{"id":"arxiv/2605.09615","url":"https://arxiv.org/abs/2605.09615","tags":["math.NA"],"score":0,"title":"Discrete positivity and maximum principles for a finite element discretization of the Richards equation","domain":"math.NA","source":"arxiv","authors":["Abderrahmane Benfanich","Yves Bourgault","Abdelaziz Beljadid"],"reviews":0,"subject":"math.NA","summary":"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…","category":"research","comments":0,"disputed":false,"division":"math","mentions":0,"discussed":0,"published":"2026-05-10","reactions":{}}