Discrete positivity and maximum principles for a finite element discretization of the Richards equation
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866910207239520256 |
|---|---|
| author | Benfanich, Abderrahmane Bourgault, Yves Beljadid, Abdelaziz |
| author_facet | Benfanich, Abderrahmane Bourgault, Yves Beljadid, Abdelaziz |
| contents | Standard finite element discretizations of the Richards equation may violate the discrete minimum principle, producing unphysical negative saturations. While existing bound-preserving methods typically rely on computationally expensive fully implicit solvers, we propose a novel semi-implicit finite element framework utilizing a bounded continuous auxiliary variable. Our approach treats the gravity-driven advective term using a linearly implicit technique, which improves the time-step restrictions required by explicit gravity methods near the degenerate limit. We provide rigorous mathematical proofs establishing sufficient geometric and algebraic constraints for discrete positivity and the discrete maximum principle, specifically a local Péclet condition and a discrete row-sum condition. When both conditions are satisfied on weakly acute meshes with mass lumping, our framework ensures that numerical solutions strictly respect physical bounds across highly degenerate conditions and initially dry soil regimes. Comprehensive numerical validation demonstrates the method across multiple flow regimes, including cases where algebraic conditions are satisfied, violated, and recovered through mesh refinement. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2605_09615 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Discrete positivity and maximum principles for a finite element discretization of the Richards equation Benfanich, Abderrahmane Bourgault, Yves Beljadid, Abdelaziz Numerical Analysis Standard finite element discretizations of the Richards equation may violate the discrete minimum principle, producing unphysical negative saturations. While existing bound-preserving methods typically rely on computationally expensive fully implicit solvers, we propose a novel semi-implicit finite element framework utilizing a bounded continuous auxiliary variable. Our approach treats the gravity-driven advective term using a linearly implicit technique, which improves the time-step restrictions required by explicit gravity methods near the degenerate limit. We provide rigorous mathematical proofs establishing sufficient geometric and algebraic constraints for discrete positivity and the discrete maximum principle, specifically a local Péclet condition and a discrete row-sum condition. When both conditions are satisfied on weakly acute meshes with mass lumping, our framework ensures that numerical solutions strictly respect physical bounds across highly degenerate conditions and initially dry soil regimes. Comprehensive numerical validation demonstrates the method across multiple flow regimes, including cases where algebraic conditions are satisfied, violated, and recovered through mesh refinement. |
| title | Discrete positivity and maximum principles for a finite element discretization of the Richards equation |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2605.09615 |