pith. sign in

arxiv: 2604.21558 · v1 · submitted 2026-04-23 · 🧮 math.NA · cs.NA

A nonconforming method for a generalized Darcy-Forchheimer model

Michele Botti , Lorenzo Mascotto , Marialetizia Mosconi This is my paper

Pith reviewed 2026-05-09 21:25 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords nonconforming finite elementsDarcy-Forchheimer modelmixed methodsconvergence analysiserror estimatesbroken Sobolev spacesporous media flownonlinear permeability
0
0 comments X

The pith

A nonconforming dual mixed scheme converges for the generalized Darcy-Forchheimer model under low regularity assumptions.

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

The paper constructs and analyzes a nonconforming finite element discretization for a generalized Darcy-Forchheimer model of nonlinear porous media flow. It extends earlier quadratic-only schemes by admitting general Forchheimer nonlinearities, mixed inhomogeneous boundary conditions, and permeability tensors with reduced Lebesgue integrability. Convergence of the discrete solutions to the exact weak solution is proved without high regularity on the data or solution, using newly derived Sobolev trace inequalities that control jumps across element interfaces in broken spaces. When extra regularity is assumed, the analysis supplies error bounds of arbitrary order, and numerical tests confirm the method behaves well across different nonlinear exponents and solver choices.

Core claim

We analyze a dual mixed nonconforming discretization of a generalized Darcy-Forchheimer model. Compared to prior quadratic schemes, the method handles general nonquadratic Forchheimer nonlinearities, mixed inhomogeneous boundary conditions, and permeability tensors with lower Lebesgue regularity. General-order schemes are built, convergence to the exact solution is established under low regularity assumptions via novel Sobolev-trace inequalities for broken spaces, and general-order error estimates are derived under additional regularity of the solution and data, with numerical results assessing performance for varied nonlinearities and solvers.

What carries the argument

A dual mixed nonconforming finite element method on broken polynomial spaces that approximates the flux and pressure separately while preserving local mass conservation and accommodating the monotonicity of the generalized Forchheimer term.

If this is right

  • The discrete solutions converge to the exact solution whenever the data and permeability satisfy only the minimal integrability required for well-posedness.
  • Arbitrary-order error bounds hold once the solution and forcing terms possess sufficient additional smoothness.
  • The scheme remains stable for nonquadratic Forchheimer exponents and for mixed boundary conditions that are inhomogeneous.
  • Permeability tensors belonging to spaces with lower integrability than L-infinity are admissible without destroying convergence.

Where Pith is reading between the lines

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

  • The broken-space trace inequalities could be reused to analyze nonconforming schemes for other nonlinear elliptic problems in porous-media or filtration models.
  • The reduced regularity demands may allow the method to be applied directly to field data from reservoirs or aquifers where smooth coefficients are unavailable.
  • Coupling the scheme with adaptive mesh refinement driven by the same broken norms could produce efficient solvers for large-scale three-dimensional flows.

Load-bearing premise

The novel Sobolev-trace inequalities for broken spaces hold and supply the bounds needed to pass to the limit in the convergence argument.

What would settle it

A concrete manufactured solution with low regularity whose discrete approximations fail to converge in the broken norm as the mesh is refined, or a permeability and nonlinearity pair where the discrete problem loses uniqueness while the continuous problem remains well-posed.

Figures

Figures reproduced from arXiv: 2604.21558 by Lorenzo Mascotto, Marialetizia Mosconi, Michele Botti.

Figure 1
Figure 1. Figure 1: Exact solution in (97); h-convergence of the flux ( [PITH_FULL_IMAGE:figures/full_fig_p026_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Exact solution in (97); h-convergence of the flux ( [PITH_FULL_IMAGE:figures/full_fig_p026_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Exact solution in (97); h-convergence of the flux ( [PITH_FULL_IMAGE:figures/full_fig_p027_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Exact solution in (97); h-convergence of the flux ( [PITH_FULL_IMAGE:figures/full_fig_p028_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Exact solution in (97); h-convergence of the flux ( [PITH_FULL_IMAGE:figures/full_fig_p028_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Exact solution u2 in (98); h-convergence of the flux (left-panel) and the potential (right-panel) for k = 1, 2, α = 3, and β = 10. Standard fixed point iterative scheme (94). different directions: • (α-2)-type Forchheimer nonlinearities are allowed, moving beyond the standard quadratic model; • the proposed method accommodates mixed, inhomogeneous boundary conditions, and per￾mits permeability tensors with… view at source ↗
read the original abstract

We analyze a dual mixed nonconforming discretization of a generalized Darcy-Forchheimer model. Compared to the analogous scheme proposed by Girault and Wheeler, we consider general, i.e., nonquadratic, Forchheimer nonlinearities; we admit mixed, inhomogeneous boundary conditions; we allow for more general, i.e., with lower Lebesgue regularity, permeability tensors; we construct general-order schemes; we prove convergence to the exact solution under low regularity assumptions, based on novel Sobolev-trace inequalities for broken spaces; we derive error estimates of general-order assuming extra regularity of the exact solution and data; we present numerical results assessing the performance of the proposed schemes for different types of nonlinearity and nonlinear solvers.

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

0 major / 3 minor

Summary. The paper analyzes a dual-mixed nonconforming finite element discretization for a generalized Darcy-Forchheimer model. It extends the scheme of Girault and Wheeler to general (non-quadratic) Forchheimer nonlinearities, inhomogeneous mixed boundary conditions, permeability tensors with lower Lebesgue regularity, and arbitrary-order approximations. Convergence to the exact solution is proved under low regularity assumptions using novel Sobolev-trace inequalities on broken spaces; general-order a priori error estimates are derived under additional regularity of the data and solution; numerical experiments assess performance across nonlinearity types and solvers.

Significance. If the novel trace inequalities and monotonicity arguments hold, the work meaningfully broadens the applicability of nonconforming methods to nonlinear porous-media flows with realistic boundary data and reduced regularity. The general-order error analysis and numerical validation provide a solid foundation for further development in the field.

minor comments (3)
  1. §3.2: the statement of the discrete inf-sup condition for the nonconforming space would benefit from an explicit reference to the precise mesh assumptions (e.g., shape-regularity constant) used in the proof.
  2. Table 1: the reported L2 velocity errors for the cubic nonlinearity appear to converge at a rate slightly below the predicted order; a brief remark on whether this is due to the solver tolerance or the nonlinearity would clarify the results.
  3. The notation for the broken Sobolev spaces (e.g., the precise definition of the jump operator across edges) is introduced in §2 but used without re-statement in the error analysis of §5; a short reminder would improve readability.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and the recommendation for minor revision. The report accurately summarizes the contributions of our work on the nonconforming dual-mixed scheme for the generalized Darcy-Forchheimer model, including the novel trace inequalities, convergence under low regularity, and general-order error estimates. We will incorporate minor revisions to improve the presentation and clarity of the manuscript.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The manuscript establishes convergence of the dual-mixed nonconforming scheme to an external exact solution via novel Sobolev-trace inequalities on broken spaces that are proved inside the paper, together with monotonicity/coercivity arguments on the generalized Forchheimer term. Error estimates of arbitrary order are derived under stated extra regularity assumptions on the data and solution. No derivation step reduces by construction to a fitted parameter, a self-referential definition, or a load-bearing self-citation whose content is merely renamed or assumed without independent verification. The argument chain remains self-contained against external solution benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only the abstract is available; the ledger is therefore limited to the structural assumptions explicitly required by the model statement.

axioms (2)
  • domain assumption The Forchheimer nonlinearity satisfies standard monotonicity and growth conditions that guarantee well-posedness of the continuous problem.
    Invoked to ensure existence and uniqueness before discretization.
  • domain assumption The permeability tensor belongs to L^∞ or weaker Lebesgue spaces as stated.
    Explicitly relaxed compared with prior work.

pith-pipeline@v0.9.0 · 5413 in / 1342 out tokens · 21573 ms · 2026-05-09T21:25:49.603964+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

40 extracted references · 40 canonical work pages

  1. [1]

    Ainsworth and R

    M. Ainsworth and R. Rankin. Fully computable bounds for t he error in nonconforming finite element approx- imations of arbitrary order on triangular elements. SIAM J. Numer. Anal. , 46(6):3207–3232, 2008

    work page 2008

  2. [2]

    J. A. Almonacid, H. S. D ´ ıaz, G. N. Gatica, and A. M´ arquez . A fully mixed finite element method for the coupling of the Stokes and Darcy-Forchheimer problems. IMA J. Numer. Anal. , 40(2):1454–1502, 2020. 30

    work page 2020

  3. [3]

    Amigo, F

    D. Amigo, F. Lepe, E. Ot´ arola, and G. Rivera. A virtual el ement method for a convective Brinkman- Forchheimer problem coupled with a heat equation. Comput. Math. Appl. , 191:1–23, 2025

    work page 2025

  4. [4]

    Badia, C

    S. Badia, C. Carstensen, A. F. Martin, R. Ruiz-Baier, and S. Villa-Fuentes. A velocity-vorticity-pressure formulation for the steady Navier-Stokes-Brinkman-Forch heimer problem. Comput. Methods Appl. Mech. Engrg., 447:Paper No. 118343, 27, 2025

    work page 2025

  5. [5]

    A. Kh. Balci, Ch. Ortner, and J. Storn. Crouzeix-Raviart finite element method for non-autonomous variational problems with Lavrentiev gap. Numer. Math. , 151(4):779–805, 2022

    work page 2022

  6. [6]

    Baran and G

    ´A. Baran and G. Stoyan. Gauss-Legendre elements: a stable, h igher order non-conforming finite element family. Computing, 79(1):1–21, 2007

    work page 2007

  7. [7]

    J. W. Barrett and W. B. Liu. Finite element approximation of the p-Laplacian. Math. Comp. , 61(204):523–537, 1993

    work page 1993

  8. [8]

    V. Berinde. Iterative Approximation of Fixed Points , volume 1912 of Lecture Notes in Mathematics . Springer, Berlin, second edition, 2007

    work page 1912

  9. [9]

    Bonetti, M

    S. Bonetti, M. Botti, and P. F. Antonietti. Conforming an d discontinuous discretizations of non-isothermal Darcy-Forchheimer flows. http://arxiv.org/abs/2508.21630, 2025

    work page arXiv 2025

  10. [10]

    Botti and L

    M. Botti and L. Mascotto. Trace inequalities for piecew ise W 1,p functions over general polytopic meshes. http://arxiv.org/abs/2512.09752, 2025

    work page arXiv 2025

  11. [11]

    Botti and L

    M. Botti and L. Mascotto. Sobolev-Poincar´ e inequalit ies for piecewise W 1,p functions over general polytopic meshes. Accepted on SIAM J. Numer. Anal. , 2026

    work page 2026

  12. [12]

    S. C. Brenner and L. R. Scott. The Mathematical Theory of Finite Element Methods , volume 15 of Texts in Applied Mathematics . Springer, New York, third edition, 2008

    work page 2008

  13. [13]

    Bressan, L

    A. Bressan, L. Mascotto, and M. Mosconi. New Crouzeix-R aviart elements of even degree: theoretical aspects, numerical performance, and applications to the St okes’ equations. IMA J. Numer. Anal , 2026. https://doi.org/10.1093/imanum/draf091

    work page doi:10.1093/imanum/draf091 2026

  14. [14]

    Burman and P

    E. Burman and P. Hansbo. Stabilized Crouzeix-Raviart e lement for the Darcy-Stokes problem. Numer. Methods Partial Differential Equations , 21(5):986–997, 2005

    work page 2005

  15. [15]

    Carstensen and S

    C. Carstensen and S. A. Sauter. Critical functions and i nf-sup stability of Crouzeix-Raviart elements. Comput. Math. Appl. , 108:12–23, 2022

    work page 2022

  16. [16]

    Caucao, G

    S. Caucao, G. N. Gatica, and L. F. Gatica. A Banach spaces -based mixed finite element method for the stationary convective Brinkman-Forchheimer problem. Calcolo, 60(4):Paper No. 51, 32, 2023

    work page 2023

  17. [17]

    Caucao, G

    S. Caucao, G. N. Gatica, and F. Sandoval. A fully-mixed fi nite element method for the coupling of the Navier- Stokes and Darcy-Forchheimer equations. Numer. Methods Partial Differential Equations , 37(3):2550–2587, 2021

    work page 2021

  18. [18]

    Crouzeix and P.-A

    M. Crouzeix and P.-A. Raviart. Conforming and nonconfo rming finite element methods for solving the station- ary Stokes equations. I. Rev. Fran. ¸ Automat. Informat. Rech. Op´ er. S´ er. Rouge, 7:33–75, 1973

    work page 1973

  19. [19]

    Douglas, Jr., P

    J. Douglas, Jr., P. J. Paes-Leme, and T. Giorgi. General ized Forchheimer flow in porous media. In Boundary value problems for partial differential equations and appli cations, volume 29 of RMA Res. Notes Appl. Math. , pages 99–111. Masson, Paris, 1993

    work page 1993

  20. [20]

    Ern and J.-L

    A. Ern and J.-L. Guermond. Finite Elements I: Approximation and Interpolation , volume 72. Springer Nature, 2021

    work page 2021

  21. [21]

    P. H. Forchheimer. W asserbewegung durch Boden. Z. Ver. Deutsch. Ing. , 45:1782–1788, 1901

    work page 1901

  22. [22]

    Fortin and M

    M. Fortin and M. Soulie. A nonconforming piecewise quad ratic finite element on triangles. Internat. J. Numer. Methods Engrg., 19(4):505–520, 1983

    work page 1983

  23. [23]

    Girault and M

    V. Girault and M. F. Wheeler. Numerical discretization of a Darcy-Forchheimer model. Numer. Math. , 110(2):161–198, 2008

    work page 2008

  24. [24]

    Kim and E.-J

    M.-Y. Kim and E.-J. Park. Fully discrete mixed finite ele ment approximations for non-Darcy flows in porous media. Comput. Math. Appl. , 38(11-12):113–129, 1999

    work page 1999

  25. [25]

    L´ opez, B

    H. L´ opez, B. Molina, and J. J. Salas. Comparison betwee n different numerical discretizations for a Darcy- Forchheimer model. Electron. Trans. Numer. Anal. , 34:187–203, 2008/09

    work page 2008

  26. [26]

    E.-J. Park. Mixed finite element methods for generalize d Forchheimer flow in porous media. Numer. Methods Partial Differential Equations , 21(2):213–228, 2005

    work page 2005

  27. [27]

    Raviart and J

    P.-A. Raviart and J. M. Thomas. Primal hybrid finite elem ent methods for 2nd order elliptic equations. Math. Comp., 31(138):391–413, 1977

    work page 1977

  28. [28]

    Rui and H

    H. Rui and H. Pan. A block-centered finite difference meth od for the Darcy-Forchheimer model. SIAM J. Numer. Anal. , 50(5):2612–2631, 2012

    work page 2012

  29. [29]

    Rui and H

    H. Rui and H. Pan. A block-centered finite difference meth od for slightly compressible Darcy-Forchheimer flow in porous media. J. Sci. Comput. , 73(1):70–92, 2017

    work page 2017

  30. [30]

    J. J. Salas, H. L´ opez, and B. Molina. An analysis of a mix ed finite element method for a Darcy-Forchheimer model. Math. Comput. Modelling , 57(9-10):2325–2338, 2013. 31

    work page 2013

  31. [31]

    Sayah, G

    T. Sayah, G. Semaan, and F. Triki. Finite element method s for the Darcy-Forchheimer problem coupled with the convection-diffusion-reaction problem. ESAIM Math. Model. Numer. Anal. , 55(6):2643–2678, 2021

    work page 2021

  32. [32]

    R. E. Showalter. Monotone Operators in Banach Space and Nonlinear Partial Di fferential Equations , volume 49 of Mathematical Surveys and Monographs . American Mathematical Society, Providence, RI, 1997

    work page 1997

  33. [33]

    Stoyan and ´A

    G. Stoyan and ´A. Baran. Crouzeix-Velte decompositions for higher-order finite elements. Comput. Math. Appl. , 51(6-7):967–986, 2006

    work page 2006

  34. [34]

    Triki, T

    F. Triki, T. Sayah, and G. Semaan. A posteriori error est imates for Darcy-Forchheimer’s problem coupled with the convection-diffusion-reaction equation. Int. J. Numer. Anal. Model. , 21(1):65–103, 2024

    work page 2024

  35. [35]

    W ang and H

    Y. W ang and H. Rui. Stabilized Crouzeix-Raviart elemen t for Darcy-Forchheimer model. Numer. Methods Partial Differential Equations , 31(5):1568–1588, 2015

    work page 2015

  36. [36]

    Zhang and H

    J. Zhang and H. Rui. A stabilized Crouzeix-Raviart elem ent method for coupling Stokes and Darcy-Forchheimer flows. Numer. Methods Partial Differential Equations , 33(4):1070–1094, 2017

    work page 2017

  37. [37]

    L. Zhao, E. T. Chung, E.-J. Park, and G. Zhou. Staggered D G method for coupling of the Stokes and Darcy- Forchheimer problems. SIAM J. Numer. Anal. , 59(1):1–31, 2021. A Properties of the nonlinear operator This section is devoted with showing three properties of A(·):

    work page 2021

  38. [38]

    monotonicity, in Section A.1

    work page

  39. [39]

    coercivity, in Section A.2

    work page

  40. [40]

    A.1 Monotonicity Here, we show that A(·) is monotone

    hemi-continuity, in Section A.3. A.1 Monotonicity Here, we show that A(·) is monotone. Proposition A.1. Let y be in V. Then, recalling that λ min is the smallest eigenvalue of K − 1, the following estimate holds true: µ ρλ min∥w − v∥L2(Ω) ≤ ∫ Ω (A(w + y) − A(v + y)) ·(w − v) dx ∀ w, v ∈ V. Proof. We define J : Lα (Ω) → R as J (v) := 1 2 µ ρ ∫ Ω (K− 1)v ·v ...

    work page