Journal Article10.1029/WR026I002P00279
Efficient numerical methods for infiltration using Richards' equation
192
TL;DR: In this paper, two efficient finite difference methods for solving the Richards' equation in one dimension are presented, and their use in a range of soils and conditions is investigated, and the methods add points to the space grid as an infiltration or redistribution front advances.
read more
Abstract: Two efficient finite difference methods for solving Richards' equation in one dimension are presented, and their use in a range of soils and conditions is investigated. Large time steps are possible when the mass-conserving mixed form of Richards' equation is combined with an implicit iterative scheme, while a hyperbolic sine transform for the matric potential allows large spatial increments even in dry, inhomogeneous soil. Infiltration in a range of soils can be simulated in a few seconds on a personal computer with errors of only a few percent in the amount and distribution of soil water. One of the methods adds points to the space grid as an infiltration or redistribution front advances, thus gaining considerably in efficiency over the other fixed grid method for infiltration problems. In 17-s computing, this advancing front method simulated infiltration, redistribution, and drainage for 50 days in an inhomogeneous soil with nonuniform initial conditions. Only 16 space and 21 time steps were needed for the simulation, which included early ponding with the development and dissipation of a perched water table.
read more
Chat with Paper
AI Agents for this Paper
Find similar papers on Google Scholar, PubMed and Arxiv
Write a critical review of this paper
Analyze citations of this paper to find unaddressed research gaps
Citations
A dual-porosity model for simulating the preferential movement of water and solutes in structured porous media
TL;DR: In this paper, a dual-porosity model was developed for the purpose of studying variably saturated water flow and solute transport in structured soils or fractured rocks, where water in both pore systems is assumed to be mobile.
1.2K
APSIM: a novel software system for model development, model testing and simulation in agricultural systems research
TL;DR: A shell of software tools has been developed within a WINDOWS environment which includes user-installed editor, linker, compiler, testbed generator, graphics, database and version control software and the resulting product is one in which the functions are coded in the language most familiar to the developers of scientific modules but provides many of the features of object oriented programming.
951
An evaluation of lattice Boltzmann schemes for porous medium flow simulation
C. Pan,Li-Shi Luo,Cass T. Miller +2 more
TL;DR: In this article, the capability and accuracy of the lattice Boltzmann equation (LBE) for modeling flow through porous media was evaluated with the multiple-relaxation-time (MRT) and Bhatnagar-Gross-Krook (BGK) collision operators.
781
Numerical simulation of infiltration, evaporation and shallow groundwater levels with the Richards equation
J.C. van Dam,Reinder A. Feddes +1 more
TL;DR: In this paper, the authors discuss a versatile, implicit, backward finite difference scheme which is relatively easy to implement and special attention is given to the selection of a head or flux controlled top boundary condition during the iterative solution of the Richards equation.
397
A comparison of Picard and Newton iteration in the numerical solution of multidimensional variably saturated flow problems
Claudio Paniconi,Mario Putti +1 more
TL;DR: In this paper, the authors compared the performance of the Picard and Newton iterative methods in one-, two-, and three-dimensional finite element simulations involving both steady state and transient flow, and concluded that the Picard or relaxed Picard schemes are often adequate for solving Richards' equation, but in cases where these fail to converge or converge slowly, the Newton method should be used.
343
References
Numerical recipes
TL;DR: This section discusses a particular type of low-pass filter, well-adapted for data smoothing, and termed variously Savitzky-Golay, least-squares, or DISPO (Digital Smoothing Polynomial) filters.
19K
•Book
Finite Element Simulation in Surface and Subsurface Hydrology
George F. Pinder,William G. Gray +1 more
- 01 Jan 1977
742
Determining soil hydraulic properties from one-step outflow experiments by parameter estimation. I: Theory and numerical studies
TL;DR: In this article, the numerical feasibility of determining water retention and hydraulic conductivity functions simultaneously from one-step pressure outflow experiments on soil cores by a parameter estimation method is evaluated.
323
A mass-conservative procedure for time-stepping in models of unsaturated flow
TL;DR: In numerical models of moisture in the unsaturated zone, it is frequently desirable to choose the matric (or pressure) head, rather than the moisture content, as the dependent variable as discussed by the authors.
195