The strength of RE is its ability to represent matrix drainage and capillary flows, which controlĮvapotranspiration processes, and its ability to be coupled to heat and solute transport models. Percolation processes (Beven and Germann, 2013). Practical limitations in representing the flow processes in real soils containing macropores, especially with modeling infiltration and rapid RE has been extensively reviewed (e.g., Vereecken et al., 2016 Farthing and Ogden, 2017, and references therein). Numerical solution yet RE remains challenging to solve reliably and accurately for a given set of boundary conditions and soil hydraulic properties Almost any practical application of RE requires a Richards' equation (RE) describes the movement of water in variably saturated porous media. Represent alternative properties (e.g., novel ways to parameterize the K( ψ) relationship) and processes (e.g., it is straightforward toĬouple heat or solute transport), making it ideal for testing alternative hypotheses. ThereĪre two advantages of our approach: (i) the code is concise, making it ideal for teaching purposes and (ii) the method can be easily extended to SciPy ODE solvers that are available in Python and make practical recommendations about the best way to implement them for Richards' equation. Solution can be implemented in an interpreted language, such as MATLAB or Python, making use of off-the-shelf ODE solvers.

Hydrus and better water balance performance (though both models have excellent water balance closure for all the problems we considered). Incur significant truncation errors in the simulated boundary fluxes, depending on the time steps used. We find that Celia's solution has the best water balance, but it can

We reproduce a set of benchmark solutions with all models. We compare this solution with alternatives, including the classic modified Picard iteration method and Weĭeveloped a novel method to quantify the boundary fluxes that reduce water balance errors without negative impacts on model runtimes – the

We implement such a solution in our model openRE. Provide efficient, mass-conservative solutions to the pressure-head form of Richards' equation. A simple numerical solution procedure – namely the method of lines combined with an off-the-shelf ordinary differential equation (ODE) solver – was shown in previous work to

Instal the new Hydrus Network 535