The finite element method is an attractive numerical method for modeling water flow in variably saturated porous media due to its flexibility in dealing with complicated geometries. It is well known that the conventional mass-distributed finite element method suffers from numerical oscillations at the wetting front, especially for very dry initial conditions. Routinely, mass-lumped procedures are used to eliminate them. This paper proposes a physical interpretation of the finite element method applied to the water flow problem. With the finite element method, mass conservation is applied at the element level. The water storage and the flux within each element are split into several components in the function space, each of which corresponds to one component of the boundary flux of the element. However, even though physical laws are correctly applied at the element level, it is shown that the traditional mass-distributed scheme can still generate an incorrect neighboring node response due to the highly nonlinear properties of water flow in unsaturated soil and cause numerical oscillation. We propose two new mass-distributed schemes which are free of the numerical oscillation, and reduce the smearing near the wetting front, at slightly increased CPU time.
[1]
Temam.
Survey of the status of finite element methods for partial differential equations. Final report
,
1986
.
[2]
M. Celia,et al.
A General Mass-Conservative Numerical Solution for the Unsaturated Flow Equation
,
1990
.
[3]
P. Milly,et al.
A mass-conservative procedure for time-stepping in models of unsaturated flow
,
1985
.
[4]
M. Th. van Genuchten,et al.
A comparison of numerical solutions of the one-dimensional unsaturated—saturated flow and mass transport equations
,
1982
.
[5]
Van Genuchten,et al.
A closed-form equation for predicting the hydraulic conductivity of unsaturated soils
,
1980
.
[6]
Linda M. Abriola,et al.
Mass conservative numerical solutions of the head‐based Richards equation
,
1994
.