Efficient numerical methods for infiltration using Richards' equation

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.

Impacts of elevated CO 2 , climate change and their interactions on water budgets in four different catchments in Australia

“…The WAVES model was chosen because (1) it has a coupled 160 water-energy-carbon modelling structure with a good balance of complexity and accuracy in 161 both hydrological and physiological processes and thus can capture the dynamic coupling of 162 water and carbon fluxes within soil-plant-atmosphere continuum and suitable for this study 163 In WAVES, soil water movement in both the unsaturated and saturated zones is 167 simulated by numerically solving the Richards equation using a finite difference method 168 (Ross, 1990; Dawes and Short, 1993). For each soil type, an analytical soil model proposed 169…”

MT3D-USGS version 1: A U.S. Geological Survey release of MT3DMS updated with new and expanded transport capabilities for use with MODFLOW

“…This approach neglects the diffusive term found in Richards' equation (Richards, 1931), resulting in sharp wetting fronts and precludes the ability to simulate capillary pressures that draw water upward. However, the approach is numerically efficient and stable; the method has provided reasonable results where Richards' equation has historically struggled (Harman and others, 2011;McGrath and others, 2008a;McGrath and others, 2008b;Ross, 1990;Struthers and others, 2006;Van Dam and Feddes, 2000) and therefore is a viable alternative.…”

Case Study: Finite Element Method and Artificial Neural Network Models for Flow through Jeziorsko Earthfill Dam in Poland

“…The value of C͑h͒ determined from Eq. ͑3͒ leads to large errors in the numerical model ͑Celia et Ross 1990;Paniconi et al 1991;Li 1993;Rathfelder and Abriola 1994;and Tocci et al 1997͒. Therefore C͑h͒ was evaluated effectively by following Cooley ͑1983͒ and Abriola and Rathfelder ͑1993͒ as…”