pith. sign in

arxiv: 2602.19037 · v2 · submitted 2026-02-22 · 🧮 math.AP

Mathematical analysis for a doubly degenerate parabolic equation: Application to the Richards equation

Abderrahmane Benfanich , Yves Bourgault , Abdelaziz Beljadid This is my paper

Pith reviewed 2026-05-15 20:46 UTC · model grok-4.3

classification 🧮 math.AP
keywords doubly degenerate parabolic equationRichards equationweak solutionsmaximal monotone operatorsL-schemeweighted Sobolev spacessaturation bounds
0
0 comments X

The pith

A bounded auxiliary variable enables existence of weak solutions for the doubly degenerate Richards equation in weighted Sobolev spaces.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper analyzes a doubly degenerate parabolic equation that arises when modeling unsaturated flow through porous media with the Richards equation. The authors introduce a bounded auxiliary variable to recast the problem so that semi-implicit time discretization combined with maximal monotone operator theory can be applied in weighted Sobolev spaces. They obtain existence of weak solutions without assuming strictly positive lower bounds on the diffusivity coefficient or high regularity of the solution. The same reformulation is shown to keep the saturation variable inside its physical interval [0,1] at all times. In addition, the L-scheme linearization converges linearly and unconditionally to the semi-discrete solution.

Core claim

The Richards equation, after reformulation with a bounded auxiliary variable, admits weak solutions whose existence follows from semi-implicit time discretization and maximal monotone operator theory in weighted Sobolev spaces. The analysis requires neither a strictly positive lower bound on the diffusivity nor high regularity of the solution. The auxiliary variable guarantees that saturation remains between zero and one, and the L-scheme linearization converges linearly to the semi-discrete solution without any restriction on the time-step size.

What carries the argument

A bounded auxiliary variable that recasts the Richards equation into a form where maximal monotone operator theory applies directly in weighted Sobolev spaces.

If this is right

  • Weak solutions exist for the Richards equation under the stated monotonicity and growth conditions on the nonlinearities.
  • The saturation variable is guaranteed to remain inside the physical interval [0,1] for all time.
  • The L-scheme linearization converges linearly to the semi-discrete solution independently of the time-step size.
  • The existence result applies to a broader class of doubly degenerate parabolic equations possessing the same structural properties.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The weighted-space framework may allow adaptive numerical methods that refine meshes only where degeneracy occurs.
  • The same auxiliary-variable idea could be tested on other degenerate flow models such as those arising in two-phase filtration.
  • Unconditional convergence of the L-scheme suggests the method remains stable for very large time steps in long-term groundwater simulations.

Load-bearing premise

The hydraulic conductivity and retention curves must admit a bounded auxiliary variable that satisfies the monotonicity and growth conditions needed for the maximal monotone operator framework.

What would settle it

A concrete retention or conductivity curve that makes the auxiliary variable unbounded or destroys monotonicity would produce a case with no weak solution under this theory, or would cause the L-scheme to lose its unconditional linear convergence.

Figures

Figures reproduced from arXiv: 2602.19037 by Abdelaziz Beljadid, Abderrahmane Benfanich, Yves Bourgault.

Figure 1
Figure 1. Figure 1: Extension of the saturation function θ(η) for parameters c = 5 3 and b = 3 5 . For a homogeneous medium, the constant ϕ = θs − θr. To apply these results to the full physical model over the domain η ∈ [0, u∗ ], the conductivity is defined as: K(η) = C ϕ KsKr(θ(η))θ(η) −a , (127) where Ks is the hydraulic conductivity, C is a constant related to the hy￾draulic properties, Kr is the relative permeability, an… view at source ↗
Figure 2
Figure 2. Figure 2: Extension of the diffusivity function K(η) with parameters KsC = 1, m = 0.6, and a = 5 3 . Finally, we identify the convective term coefficient K¯ and its factorization required by Hypothesis (H3). We define the auxiliary scaling function K¯ 1 as: K¯ 1(η) = 1 C θ(η) aez, for η ∈ [0, u∗ ]. (129) We extend K¯ 1 to R similarly to K: 1. For η > u∗ , we set K¯ 1(η) = 1 C ez. 2. For η < 0, we set K¯ 1(η) = K¯ 1(… view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of the extended functions K¯ z and K¯ 1z. We suppose that the source term S = 0. Verification of Hypotheses We now verify that the extended saturation function θ and the hydraulic conductivity K constructed above satisfy the structural hypotheses (H1)– (H3) required for the existence theory. Proposition 7.1 (Verification of H1). Let b ∈ [0, 1) and c ≥ 1. The extended saturation function θ : R → … view at source ↗
read the original abstract

This paper presents a mathematical analysis of a doubly degenerate parabolic equation and its application to the Richards equation using a bounded auxiliary variable. We establish the existence of weak solutions using semi-implicit time discretization combined with maximal monotone operator theory. The analysis is conducted within weighted Sobolev spaces, allowing for a rigorous treatment of the equation's strict degeneracy and strong nonlinearities. A key feature of this study is the derivation of convergence results without imposing strictly positive lower bounds on the diffusivity or requiring high regularity of the solution. Furthermore, we prove that the Richards equation using the introduced auxiliary variable preserves the physical bounds of the saturation and demonstrate the unconditional linear convergence of the L-scheme linearization to the semi-discrete solution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript establishes existence of weak solutions to a doubly degenerate parabolic equation by means of semi-implicit time discretization combined with maximal monotone operator theory in weighted Sobolev spaces. It applies the framework to the Richards equation via a bounded auxiliary variable, proving that saturation remains in the physical interval [0,1] and that the L-scheme converges unconditionally and linearly to the semi-discrete solution without imposing strictly positive lower bounds on the diffusivity.

Significance. If the auxiliary-variable construction and the weighted-space estimates hold rigorously, the work would supply a useful rigorous justification for numerical schemes applied to degenerate parabolic models arising in porous-media flow. The avoidance of artificial positive lower bounds on diffusivity and the unconditional convergence result are potentially valuable for practical computations.

major comments (2)
  1. [§3.2] §3.2 (semi-implicit discretization and operator reformulation): the boundedness of the auxiliary variable is asserted to guarantee coercivity and maximal monotonicity in the weighted spaces, yet no explicit construction or uniform bound independent of the time step is supplied; this assumption is load-bearing for both the existence proof and the passage to the limit.
  2. [§5] §5 (Richards-equation application): the claims that saturation bounds are preserved and that the L-scheme converges linearly without discretization-dependent constants rest on the auxiliary variable satisfying the required monotonicity and growth conditions for the chosen retention and conductivity curves; the manuscript provides no verification that these conditions hold when diffusivity vanishes at the degeneracy points.
minor comments (2)
  1. [Abstract] The abstract states 'unconditional linear convergence' but does not specify the norm; this should be clarified in the statement of the main theorem.
  2. [§2] Notation for the weight function in the Sobolev spaces is introduced without an explicit formula; a displayed definition would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback on our manuscript. The points raised highlight areas where additional clarification and explicit verification will strengthen the presentation. We address each major comment below and will incorporate the necessary revisions.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (semi-implicit discretization and operator reformulation): the boundedness of the auxiliary variable is asserted to guarantee coercivity and maximal monotonicity in the weighted spaces, yet no explicit construction or uniform bound independent of the time step is supplied; this assumption is load-bearing for both the existence proof and the passage to the limit.

    Authors: We agree that the construction and bound should be made fully explicit. In the revised manuscript we will add a dedicated paragraph in §3.2 that constructs the auxiliary variable directly from the saturation via the inverse of the retention curve and then proves, via a discrete maximum principle applied to the semi-implicit scheme, that the auxiliary variable remains uniformly bounded by constants depending only on the initial data and the physical interval [0,1], independently of the time-step size. This bound is then used to verify the required coercivity and maximal monotonicity in the weighted Sobolev spaces, thereby closing the existence argument and the passage to the limit. revision: yes

  2. Referee: [§5] §5 (Richards-equation application): the claims that saturation bounds are preserved and that the L-scheme converges linearly without discretization-dependent constants rest on the auxiliary variable satisfying the required monotonicity and growth conditions for the chosen retention and conductivity curves; the manuscript provides no verification that these conditions hold when diffusivity vanishes at the degeneracy points.

    Authors: We acknowledge that explicit verification for the degeneracy case was omitted. In the revised §5 we will insert a short subsection that checks the monotonicity and growth conditions for standard retention curves (van Genuchten) and conductivity functions (Mualem) at the points where diffusivity vanishes. The verification uses the already-established uniform bound on the auxiliary variable together with the explicit algebraic form of the nonlinearities; it confirms both the preservation of saturation in [0,1] and that the Lipschitz constant of the L-scheme linearization remains independent of the time step, yielding the unconditional linear convergence result. revision: yes

Circularity Check

0 steps flagged

No significant circularity; relies on external maximal monotone theory

full rationale

The paper's core existence result for the doubly degenerate parabolic equation is obtained via semi-implicit time discretization recast as a maximal monotone operator problem in weighted Sobolev spaces. This framework is drawn from standard external theory rather than self-referential definitions or fitted parameters. The auxiliary variable is introduced with an independent argument for boundedness that preserves saturation bounds [0,1] for the Richards equation; no step equates a prediction to its own input by construction, and no load-bearing uniqueness theorem or ansatz is smuggled via self-citation. The unconditional L-scheme convergence likewise follows from the monotonicity properties without reducing to the discretization parameters themselves. This yields a low circularity score consistent with a self-contained analysis against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on standard assumptions from the theory of maximal monotone operators applied to the specific nonlinear structure of the Richards equation and on the introduction of a bounded auxiliary variable to manage strict degeneracy.

axioms (2)
  • domain assumption The nonlinear functions satisfy monotonicity and growth conditions that allow application of maximal monotone operator theory.
    Invoked when establishing existence via the semi-implicit scheme in weighted Sobolev spaces.
  • ad hoc to paper The auxiliary variable remains bounded for the given hydraulic functions.
    Central to handling strict degeneracy and proving preservation of saturation bounds.
invented entities (1)
  • bounded auxiliary variable no independent evidence
    purpose: To reformulate the doubly degenerate equation so that maximal monotone theory applies and physical saturation bounds are preserved.
    New device introduced in the paper to bypass the need for strictly positive lower bounds on diffusivity.

pith-pipeline@v0.9.0 · 5419 in / 1583 out tokens · 46581 ms · 2026-05-15T20:46:13.909614+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    L.A.Richards, Capillaryconductionofliquidsthroughporousmediums, Physics 1 (5) (1931) 318–333

    work page 1931

  2. [2]

    H. W. Alt, S. Luckhaus, Quasilinear elliptic-parabolic differential equa- tions, Math. Z. 183 (3) (1983) 311–341

    work page 1983

  3. [3]

    DiBenedetto, Degenerate Parabolic Equations, Universitext, Springer New York, New York, NY, 1993

    E. DiBenedetto, Degenerate Parabolic Equations, Universitext, Springer New York, New York, NY, 1993. 45

    work page 1993

  4. [4]

    Otto, L1-contraction and uniqueness for quasilinear elliptic–parabolic equations, Journal of Differential Equations 131 (1) (1996) 20–38

    F. Otto, L1-contraction and uniqueness for quasilinear elliptic–parabolic equations, Journal of Differential Equations 131 (1) (1996) 20–38

    work page 1996

  5. [5]

    J.Carrillo, Entropysolutionsfornonlineardegenerateproblems, Archive for Rational Mechanics and Analysis 147 (4) (1999) 269–361

    work page 1999

  6. [6]

    F. List, F. A. Radu, A study on iterative methods for solving Richards’ equation, Computational Geosciences 20 (2016) 341–353

    work page 2016

  7. [7]

    J. E. Jones, C. S. Woodward, Newton-krylov-multigrid solvers for large- scale, highly heterogeneous, variably saturated flow problems, Advances in Water Resources 24 (7) (2001) 763–774

    work page 2001

  8. [8]

    H. J. G. Diersch, P. Perrochet, On the primary variable switching tech- nique for simulating unsaturated-saturated flows, Adv. Water Resour. 23 (3) (1999) 271–301

    work page 1999

  9. [9]

    Forsyth, Y

    P. Forsyth, Y. Wu, K. Pruess, Robust numerical methods for saturated- unsaturated flow with dry initial conditions in heterogeneous media, Adv. Water Resour. 18 (1) (1995) 25–38

    work page 1995

  10. [10]

    Y. S. Wu, P. A. Forsyth, On the selection of primary variables in nu- merical formulation for modeling multiphase flow in porous media, J. Contam. Hydrol. 48 (3–4) (2001) 277–304

    work page 2001

  11. [11]

    Brunner, C

    P. Brunner, C. T. Simmons, HydroGeoSphere: a fully integrated, phys- ically based hydrological model, Ground Water 50 (2) (2012) 170–176

    work page 2012

  12. [12]

    Krabbenhøft, An alternative to primary variable switching in saturated-unsaturated flow computations, Adv

    K. Krabbenhøft, An alternative to primary variable switching in saturated-unsaturated flow computations, Adv. Water Resour. 30 (3) (2007) 483–492

    work page 2007

  13. [13]

    Y. Zha, J. Yang, L. Yin, Y. Zhang, W. Zeng, L. Shi, A modified Pi- card iteration scheme for overcoming numerical difficulties of simulating infiltration into dry soil, J. Hydrol. 551 (2017) 56–69

    work page 2017

  14. [14]

    F. H. Maina, P. Ackerer, Ross scheme, Newton–Raphson iterative meth- ods and time-stepping strategies for solving the mixed form of Richards’ equation, Hydrol. Earth Syst. Sci. 21 (6) (2017) 2667–2683

    work page 2017

  15. [15]

    C. E. Kees, C. T. Miller, Higher order time integration methods for two-phase flow, Adv. Water Resour. 25 (2) (2002) 159–177. 46

    work page 2002

  16. [16]

    Y. Zha, J. Yang, J. Zeng, C.-H. M. Tso, W. Zeng, L. Shi, Review of nu- merical solution of Richardson–Richards equation for variably saturated flow in soils, WIREs Water 6 (2019) e1364

    work page 2019

  17. [17]

    Schweizer, Regularization of outflow problems in unsaturated porous media with dryregions, J

    B. Schweizer, Regularization of outflow problems in unsaturated porous media with dryregions, J. Differential Equations237 (2) (2007) 278–306

    work page 2007

  18. [18]

    I. S. Pop, B. Schweizer, Regularization schemes for degenerate Richards equations and outflow conditions, Math. Models Methods Appl. Sci. 21 (8) (2011) 1685–1712

    work page 2011

  19. [19]

    Févotte, A

    F. Févotte, A. Rappaport, M. Vohralík, Adaptive regularization for the Richards equation, Comput. Geosci. 28 (2024) 1371–1388

    work page 2024

  20. [20]

    I. Pop, F. Radu, P. Knabner, Mixed finite elements for the Richards’ equation: linearization procedure, Journal of Computational and Ap- plied Mathematics 168 (1) (2004) 365–373, selected Papers from the Second International Conference on Advanced Computational Methods in Engineering (ACOMEN 2002)

    work page 2004

  21. [21]

    M. Slodička, A robust and efficient linearization scheme for doubly non- linearanddegenerateparabolicproblemsarisinginflowinporousmedia, SIAM Journal on Scientific Computing 23 (5) (2002) 1593–1614

    work page 2002

  22. [22]

    Mitra, I

    K. Mitra, I. S. Pop, A modified L-scheme to solve nonlinear diffusion problems, Comput. Math. Appl. 77 (6) (2019) 1722–1738

    work page 2019

  23. [23]

    Kamil, A

    H. Kamil, A. Beljadid, A. Soulaïmani, Y. Bourgault, Semi-implicit schemes for modeling water flow and solute transport in unsaturated soils, Advances in Water Resources 193 (2024) 104835

    work page 2024

  24. [24]

    Paniconi, M

    C. Paniconi, M. Putti, G. F. Pinder, Numerical evaluation of iterative and noniterative methods for the solution of the nonlinear Richards equation, Water Resources Research 27 (6) (1991) 1147–1163

    work page 1991

  25. [25]

    Keita, A

    S. Keita, A. Beljadid, Y. Bourgault, Implicit and semi-implicit second- order time stepping methods for the richards equation, Advances in Water Resources 148 (2021) 103841

    work page 2021

  26. [26]

    M. A. Celia, E. T. Bouloutas, R. L. Zarba, A general mass-conservative numerical solution for the unsaturated flow equation, Water Resour. Res. 26 (7) (1990) 1483–1496. 47

    work page 1990

  27. [27]

    Teubner Verlagsgesellschaft, 1982, pp

    K.Rektorys, Themethodofdiscretizationintimeandpartialdifferential equations, in: Equadiff 5, BSB B.G. Teubner Verlagsgesellschaft, 1982, pp. 293–296

    work page 1982

  28. [28]

    Kačur, Method of rothe in evolution equations, in: J

    J. Kačur, Method of rothe in evolution equations, in: J. Vosmanský, M. Zlámal (Eds.), Equadiff 6, Springer Berlin Heidelberg, Berlin, Hei- delberg, 1986, pp. 23–34

    work page 1986

  29. [29]

    Benfanich, Y

    A. Benfanich, Y. Bourgault, A. Beljadid, A finite element method using a bounded auxiliary variable for solving the richards equation (2025). arXiv:2510.13012. URLhttps://arxiv.org/abs/2510.13012

    work page arXiv 2025

  30. [30]

    A. C. Cavalheiro, Weighted sobolev spaces and degenerate elliptic equa- tions, Boletim da Sociedade Paranaense de Matemática (3s.) 26 (1-2) (2008) 117–132

    work page 2008

  31. [31]

    H.Brezis, Ope¦rateursmaximauxmonotonesetsemi-groupesdecontrac- tions dans les espaces de Hilbert, North-Holland Mathematics Studies, North Holland, 1973

    work page 1973

  32. [32]

    P.D.Lax, A.N.Milgram, IX.ParabolicEquations, PrincetonUniversity Press, Princeton, 1955, pp. 167–190

    work page 1955

  33. [33]

    Agarwal, M

    P. Agarwal, M. Jleli, B. Samet, Fixed Point Theory in Metric Spaces: RecentAdvancesandApplications, 1stEdition, SpringerSingapore, Sin- gapore, 2018

    work page 2018

  34. [34]

    Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer New York, New York, NY, 2011

    H. Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer New York, New York, NY, 2011

    work page 2011

  35. [35]

    J. L. Lions, E. Magenes, Non-Homogeneous Boundary Value Problems and Applications, Vol. I of Die Grundlehren der mathematischen Wis- senschaften, Springer-Verlag Berlin Heidelberg New York, 1972

    work page 1972

  36. [36]

    Aubin, Un théorème de compacité, Comptes Rendus de l’Académie des Sciences de Paris 256 (1963) 5042–5044

    J.-P. Aubin, Un théorème de compacité, Comptes Rendus de l’Académie des Sciences de Paris 256 (1963) 5042–5044

    work page 1963

  37. [37]

    H. L. Royden, Real Analysis, Prentice Hall, 1988

    work page 1988

  38. [38]

    Ern, J.-L

    A. Ern, J.-L. Guermond, Finite Elements III: First-Order and Time- Dependent PDEs, Springer International Publishing, Cham, 2021. 48

    work page 2021

  39. [39]

    Moussa, Some variants of the classical Aubin–Lions Lemma, Journal of Evolution Equations 16 (1) (2016) 65–93

    A. Moussa, Some variants of the classical Aubin–Lions Lemma, Journal of Evolution Equations 16 (1) (2016) 65–93

    work page 2016

  40. [40]

    M. T. van Genuchten, A closed-form equation for predicting the hy- draulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J. 44 (5) (1980) 892–898

    work page 1980

  41. [41]

    Mualem, A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resour

    Y. Mualem, A new model for predicting the hydraulic conductivity of unsaturated porous media, Water Resour. Res. 12 (3) (1976) 513–522. 49

    work page 1976