ISSN:
0271-2091
Keywords:
finite elements
;
Richards equation
;
porous flow
;
Engineering
;
Numerical Methods and Modeling
Source:
Wiley InterScience Backfile Collection 1832-2000
Topics:
Mechanical Engineering, Materials Science, Production Engineering, Mining and Metallurgy, Traffic Engineering, Precision Mechanics
Notes:
A new computational method is developed for numerical solution of the Richards equation for flow in variably saturated porous media. The new method, referred to as the mixed transform finite element method, employs the mixed formulation of the Richards equation but expressed in terms of a partitioned transform. An iterative finite element algorithm is derived using a Newton-Galerkin weak statement. Specific advantages of the new method are demonstrated with applications to a set of one- dimensional test problems. Comparisons with the modified Picard method show that the new method produces more robust solutions for a broad range of soil- moisture regimes, including flow in desiccated soils, in heterogeneous media and in layered soils with formation of perched water zones. In addition, the mixed transform finite element method is shown to converge faster than the modified Picard method in a number of cases and to accurately represent pressure head and moisture content profiles with very steep fronts. © 1997 by John Wiley & Sons, Ltd.
Additional Material:
10 Ill.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1002/fld.501
Permalink