Author(s)

Camille Poussel¹, Mehmet Ersoy¹, Yolhan Mannes¹, Aimed Ajroud², Frédéric Golay¹

¹Institut de Mathématiques de Toulon (IMATH), Université de Toulon, La Garde, France.
²École d’Ingénieurs SeaTech, Université de Toulon, La Garde, France.

In this article, we study the convergence of an IIPG (Incomplete Interior Penalty Galerkin) Discontinuous Galerkin numerical method for the Richards equation. The Richards equation is a degenerate parabolic nonlinear equation for modeling flows in porous media with variable saturation. The numerical solution of this equation is known to be difficult to calculate numerically, due to the abrupt displacement of the wetting front, mainly as a result of highly nonlinear hydraulic properties. As time scales are slow, implicit numerical methods are required, and the convergence of nonlinear solvers is very sensitive. We propose an original method to ensure convergence of the numerical solution to the exact Richards solution, using a technique of auto-calibration of the penalty parameters derived from the Galerkin Discontinuous method. The method is constructed using nonlinear 1D and 2D general elliptic problems. We show that the numerical solution converges toward the unique solution of the continuous problem under certain conditions on the penalty parameters. Then, we numerically demonstrate the efficiency and robustness of the method through test cases with analytical solutions, laboratory test cases, and large-scale simulations.

Porous Media, Richards Equation, Discontinuous Galerkin, Backward Differentiation Method, Incomplete Interior Penalty Galerkin (IIPG), Broken Sobolev Space, Picard’s Fixed Point, Minimal Regularity Solution

Poussel, C. , Ersoy, M. , Mannes, Y. , Ajroud, A. and Golay, F. (2025) Optimal Penalization and Nonlinear Solver Convergence for a DG-Based Richards’ Equation Model of Variably Saturated Flows. Journal of Applied Mathematics and Physics, 13, 4083-4127. doi: 10.4236/jamp.2025.1311227.