Premium
One‐dimensional infiltration with moving finite elements and improved soil water diffusivity
Author(s) -
Cox Christopher L.,
Jones Walter F.,
Quisenberry Virgil L.,
Yo Frans
Publication year - 1994
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/94wr00060
Subject(s) - vadose zone , nonlinear system , solver , richards equation , mathematics , grid , finite element method , thermal diffusivity , finite difference , partial differential equation , numerical analysis , mathematical optimization , mathematical analysis , soil water , water content , geology , geotechnical engineering , geometry , soil science , physics , quantum mechanics , thermodynamics
A problem of significant interest to environmental scientists is the flow of water and solutes through the vadose zone. The partial differential equations which govern this flow are typically time‐dependent and nonlinear. Valid solutions to these equations require (1) accurate relationships between various coefficients and variables on which they depend (e.g., coefficient of diffusivity and water content) and (2) sophisticated numerical methods which can handle complexities such as sharp moving fronts. In cases where coefficients are not known explicitly, curve‐fitting techniques are needed to smooth out scattered experimental data. Nonlinear coefficients can then be calculated. A constrained least squares spline fit is compared to empirical function fits which have appeared recently. Then, a state‐of‐the‐art numerical technique is used to accurately model transient flow through unsaturated homogeneous soils. The moving finite element method of Miller and colleagues is an adaptive approach in the sense that the grid moves so that nodes are concentrated where they are most needed. As a result, better accuracy is achieved with fewer nodes than are required for standard fixed‐grid methods. Petzold's robust Gear‐type solver DASSL is used for time‐integration. Numerical results are compared to experimental data. Mass balance errors are neglible, and accurate solutions are obtained at all time steps. Though only one‐dimensional problems are considered here, the numerical approach generalizes to heterogeneous media and problems in higher dimensions.