Journal
VADOSE ZONE JOURNAL
Volume 8, Issue 2, Pages 352-362Publisher
SOIL SCI SOC AMER
DOI: 10.2136/vzj2008.0108
Keywords
-
Categories
Funding
- GdR MoMas [CNRS-2439]
- ANDRA
- BRGM
- CEA
- EDF
Ask authors/readers for more resources
Accurate numerical simulation of infiltration in the vadose zone remains a challenge, especially when very sharp fronts are modeled. In this work, we used the mixed hybrid finite element (MHFE) method, which allows a simultaneous approximation on of both pressure head and velocity and can handle general irregular grids with highly heterogenous permeability. For many problems dealing with unsaturated water flow, however, the MHFE solutions exhibit significant unphysical oscillations. To avoid this phenomenon, we developed an efficient mass-lumping scheme with the MHFE method for solving the mixed form of the Richards equation. In this work, the standard and lumped MHFE formulations were compared and the ability of the lumped formulation to reduce unphysical oscillations was determined for one- and two-dimensional infiltration problems. Theoretical analysis based on the M-matrix property, which guarantees the discrete maximum principle, and practical test cases performed in this study underline the advantage of using an acute triangulation to completely remove the unphysical oscillations. Indeed, contrary to the standard approach, the lumped formulation satisfies the M-matrix property without any constraint on the time step size to be used.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available