pith. sign in

arxiv: 2604.24731 · v1 · submitted 2026-04-27 · 🧮 math.NA · cs.NA

Error analysis for the approximation of a flow in deformable porous media with nonlinear strain-stress relation

Andrea Bonito , Vivette Girault , Diane Guignard This is my paper

Pith reviewed 2026-05-08 02:02 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords deformable porous medianonlinear strain-stressporoelasticitysemi-implicit schemefinite element methoderror analysisconvergence estimates
0
0 comments X

The pith

Semi-implicit finite element scheme for nonlinear flow in deformable porous media has unique solutions and converges when perturbations are small.

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

The paper analyzes the numerical approximation of slow fluid flow through a porous elastic solid whose stress-strain law includes nonlinear effects formulated as a small perturbation of linear elasticity. It introduces a first-order semi-implicit time integration scheme paired with standard finite element discretization for the coupled system. The authors prove existence and uniqueness of the discrete solution together with a priori error estimates that establish convergence to the continuous solution, provided the nonlinear perturbations stay sufficiently small. This matters for applications such as geomechanics where deformation alters flow patterns and reliable error control is needed to trust the computed solutions.

Core claim

For the model of slow flow in a deformable porous medium with nonlinear perturbation of linear elasticity, the first-order semi-implicit time scheme combined with finite element discretization yields a unique discrete solution and satisfies a priori convergence estimates whenever the nonlinear perturbations remain sufficiently small.

What carries the argument

First-order semi-implicit time integration combined with standard finite element spatial discretization applied to the coupled poroelastic system with nonlinear strain-stress perturbation.

If this is right

  • The discrete solution exists and is unique under the small nonlinear perturbation assumption.
  • A priori error estimates bound the difference between discrete and continuous solutions in terms of time-step and mesh sizes.
  • Numerical experiments confirm that the scheme captures the nonlinear effects present in the model.

Where Pith is reading between the lines

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

  • Larger nonlinear perturbations would likely require a fully implicit or iteratively solved scheme to retain stability and uniqueness.
  • The same error-analysis approach can be adapted to other coupled multiphysics models that contain small nonlinear perturbations of a linear core.
  • The quantitative convergence rates can be used to select mesh and time-step sizes that keep discretization error below a target tolerance in practical simulations.

Load-bearing premise

The nonlinear perturbations to the linear elasticity law must remain sufficiently small to guarantee stability of the semi-implicit scheme and to close the error estimates.

What would settle it

Increase the size of the nonlinear perturbation coefficient in the numerical experiments until the discrete solver fails to produce a unique solution or the observed errors stop decreasing at the rates predicted by the analysis.

Figures

Figures reproduced from arXiv: 2604.24731 by Andrea Bonito, Diane Guignard, Vivette Girault.

Figure 1
Figure 1. Figure 1: Solution at the final time T = 1 using ∆ t = 0.1 for the case λ1 = λ2 = 5. Left: solid displacement u (magnitude and vector field with dark blue=0.0 and dark red=1.0); middle: fluid velocity v (magnitude and vector field with dark blue=0.0 and dark red=0.22); right: fluid pressure p with isolines using 16 values uniformly distributed between maxΩ p = −0.9868 (dark blue) and maxΩ p = 0.9945 (dark red). λ1 =… view at source ↗
Figure 2
Figure 2. Figure 2: Frobenius norm |ε(u N h )| of the strain at the final time T = 1 when λ1 = λ2 = 0, 1, 2, 3, 4, 5 (from left to right and top to bottom) using 20 values uniformly distributed isolines between 0 (dark blue) and 1.5 (dark red). [3] J Arumugam, P Alagappan, J Bird, M Moreno, and KR Rajagopal. A new constitutive relation to describe the response of bones. International Journal of Non-Linear Mechanics, 161:10466… view at source ↗
read the original abstract

We study a model describing the slow flow of a fluid through a deformable, porous, elastic solid undergoing small deformations. The stress-strain relationship of the solid incorporates nonlinear effects, formulated as a perturbation of the classical linear elasticity. To approximate the coupled system, we introduce a discrete scheme based on a first order semi-implicit time integration scheme combined with a standard finite element spatial discretization. We establish the existence and uniqueness of the discrete solution and derive a priori convergence estimates under the assumption that the nonlinear perturbations remain sufficiently small. Finally, we demonstrate the efficiency of the proposed scheme through numerical experiments that also highlight the nonlinear phenomena captured by the model.

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 manuscript analyzes a model for slow fluid flow through a deformable porous elastic solid with a nonlinear perturbation to the linear stress-strain relation. It introduces a first-order semi-implicit time-stepping scheme combined with standard finite-element spatial discretization, proves existence and uniqueness of the discrete solution, and derives a priori convergence estimates under an explicit smallness assumption on the nonlinear terms. Numerical experiments illustrate the scheme and the captured nonlinear effects.

Significance. Provided the smallness condition holds, the work supplies a rigorous, conditional error analysis for a practical discretization of a nonlinear poroelastic system. This extends classical linear theory in a transparent way and pairs the analysis with supporting computations, which is useful for applications in which nonlinear corrections remain modest.

minor comments (3)
  1. The dependence of the smallness threshold on physical data and discretization parameters should be stated explicitly in the statement of the main existence and error theorems so that uniformity with respect to mesh size and time step is immediately visible.
  2. In the numerical section, include a direct comparison against the corresponding linear-elasticity scheme (same mesh and time step) to quantify the additional cost and accuracy impact of the nonlinear term.
  3. A short remark on how the smallness assumption can be checked a posteriori from computed solutions would increase the practical utility of the theoretical result.

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. We appreciate the recognition that the work provides a rigorous conditional error analysis extending classical linear theory, paired with numerical experiments. Since the report raises no specific major comments, we will incorporate any minor editorial or presentational improvements in the revised version.

Circularity Check

0 steps flagged

No significant circularity; standard conditional analysis

full rationale

The central claims are existence/uniqueness of the discrete solution and a priori error estimates for a semi-implicit FEM scheme, both conditioned explicitly on a smallness assumption for the nonlinear perturbation. This is a conventional structure in nonlinear PDE analysis (fixed-point or contraction arguments for existence, then standard Galerkin estimates for convergence). No step reduces by construction to a fitted input, self-definition, or unverified self-citation chain; the smallness hypothesis is stated up-front and used transparently to close estimates without circular dependence on the derived rates. The derivation is self-contained against external benchmarks once the assumption is granted.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the smallness of the nonlinear perturbation (treated as a domain assumption) and on standard background results from linear elasticity and finite-element theory; no free parameters or new entities are introduced.

axioms (1)
  • domain assumption The nonlinear perturbations of the strain-stress relation remain sufficiently small.
    Explicitly invoked in the abstract to obtain existence, uniqueness, and convergence of the discrete scheme.

pith-pipeline@v0.9.0 · 5406 in / 1299 out tokens · 59022 ms · 2026-05-08T02:02:16.291649+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

23 extracted references · 23 canonical work pages

  1. [1]

    Elsevier, 2003

    Robert A Adams and John JF Fournier.Sobolev spaces, volume 140. Elsevier, 2003

    work page 2003

  2. [2]

    The deal.II finite element library: Design, features, and insights.Computers & Mathematics with Appli- cations, 81:407–422, 2021

    Daniel Arndt, Wolfgang Bangerth, Denis Davydov, Timo Heister, Luca Heltai, Martin Kro- nbichler, Matthias Maier, Jean-Paul Pelteret, Bruno Turcksin, and David Wells. The deal.II finite element library: Design, features, and insights.Computers & Mathematics with Appli- cations, 81:407–422, 2021. ERROR ANALYSIS FOR THE APPROXIMATION OF A FLOW IN A POROUS ME...

    work page 2021

  3. [3]

    A new constitutive relation to describe the response of bones.International Journal of Non-Linear Mechanics, 161:104664, 2024

    J Arumugam, P Alagappan, J Bird, M Moreno, and KR Rajagopal. A new constitutive relation to describe the response of bones.International Journal of Non-Linear Mechanics, 161:104664, 2024. 30 ANDREA BONITO, VIVETTE GIRAULT, AND DIANE GUIGNARD

    work page 2024

  4. [4]

    Continuum theories of mixtures: applications.IMA Journal of Applied Mathematics, 17(2):153–207, 1976

    RJ Atkin and RE Craine. Continuum theories of mixtures: applications.IMA Journal of Applied Mathematics, 17(2):153–207, 1976

    work page 1976

  5. [5]

    Kitware, Inc., USA, 2015

    Utkarsh Ayachit.The ParaView Guide: A Parallel Visualization Application. Kitware, Inc., USA, 2015

    work page 2015

  6. [6]

    Theory of mixtures.Continuum physics, 3:2–127, 1976

    Ray M Bowen. Theory of mixtures.Continuum physics, 3:2–127, 1976

    work page 1976

  7. [7]

    Springer Science & Business Me- dia, 2012

    Franck Boyer and Pierre Fabrie.Mathematical Tools for the Study of the Incompressible Navier-Stokes Equations and Related Models, volume 183. Springer Science & Business Me- dia, 2012

    work page 2012

  8. [8]

    A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles.Flow, Turbulence and Combustion, 1(1):27–34, 1949

    Hendrik C Brinkman. A calculation of the viscous force exerted by a flowing fluid on a dense swarm of particles.Flow, Turbulence and Combustion, 1(1):27–34, 1949

    work page 1949

  9. [9]

    On the permeability of media consisting of closely packed porous particles.Flow, Turbulence and Combustion, 1(1):81–86, 1949

    Hendrik C Brinkman. On the permeability of media consisting of closely packed porous particles.Flow, Turbulence and Combustion, 1(1):81–86, 1949

    work page 1949

  10. [10]

    Ueber diffusion.Annalen der physik, 170(1):59–86, 1855

    Adolf Fick. Ueber diffusion.Annalen der physik, 170(1):59–86, 1855

    work page

  11. [11]

    Error analysis of a BDF2-FEM discretization for flow in deformable porous media

    Vivette Girault, Mar´ ıa Gonz´ ales Taboada, Fr´ ed´ eric Hecht, and Kumbakonam R Rajagopal. Error analysis of a BDF2-FEM discretization for flow in deformable porous media

    work page

  12. [12]

    A model for flow in deformable porous media.Communications in optimization Theory, 20:1–44, 2023

    Vivette Girault, Mar´ ıa Gonz´ ales Taboada, Fr´ ed´ eric Hecht, and Kumbakonam R Rajagopal. A model for flow in deformable porous media.Communications in optimization Theory, 20:1–44, 2023

    work page 2023

  13. [13]

    SIAM, 2011

    Pierre Grisvard.Elliptic problems in nonsmooth domains. SIAM, 2011

    work page 2011

  14. [14]

    American Mathematical Soc., 2017

    Giovanni Leoni.A first course in Sobolev spaces. American Mathematical Soc., 2017

    work page 2017

  15. [15]

    On implicit constitutive theories.Applications of Mathematics, 48(4):279–319, 2003

    Kumbakonam R Rajagopal. On implicit constitutive theories.Applications of Mathematics, 48(4):279–319, 2003

    work page 2003

  16. [16]

    The elasticity of elasticity.Zeitschrift f¨ ur angewandte Mathe- matik und Physik, 58(2):309–317, 2007

    Kumbakonam R Rajagopal. The elasticity of elasticity.Zeitschrift f¨ ur angewandte Mathe- matik und Physik, 58(2):309–317, 2007

    work page 2007

  17. [17]

    On a hierarchy of approximate models for flows of incompress- ible fluids through porous solids.Mathematical Models and Methods in Applied Sciences, 17(02):215–252, 2007

    Kumbakonam R Rajagopal. On a hierarchy of approximate models for flows of incompress- ible fluids through porous solids.Mathematical Models and Methods in Applied Sciences, 17(02):215–252, 2007

    work page 2007

  18. [18]

    Conspectus of concepts of elasticity.Mathematics and Mechanics of Solids, 16(5):536–562, 2011

    Kumbakonam R Rajagopal. Conspectus of concepts of elasticity.Mathematics and Mechanics of Solids, 16(5):536–562, 2011

    work page 2011

  19. [19]

    Kumbakonam R Rajagopal. An implicit constitutive relation for describing the small strain response of porous elastic solids whose material moduli are dependent on the density.Math- ematics and Mechanics of Solids, 26(8):1138–1146, 2021

    work page 2021

  20. [20]

    World scien- tific, 1995

    Kumbakonam R Rajagopal and Luoyi Tao.Mechanics of mixtures, volume 35. World scien- tific, 1995

    work page 1995

  21. [21]

    Thermodynamics of irreversible processes in fluid mixtures (approached by rational thermodynamics).(No Title), 1987

    Ivan Samoh` yl. Thermodynamics of irreversible processes in fluid mixtures (approached by rational thermodynamics).(No Title), 1987

    work page 1987

  22. [22]

    Navier-Stokes equations: Theory and numerical analysis.Amsterdam, North- Holland Publishing Co.(Studies in Mathematics and Its Applications, 2:510, 1977

    Roger Temam. Navier-Stokes equations: Theory and numerical analysis.Amsterdam, North- Holland Publishing Co.(Studies in Mathematics and Its Applications, 2:510, 1977

    work page 1977

  23. [23]

    Sulle basi della termomeccanica.Rend

    Clifford Truesdell. Sulle basi della termomeccanica.Rend. Lincei, 22(8):33–38, 1957. (Andrea Bonito) Texas A&M University, Department of Mathematics, College Sta- tion, TX 77843, USA Email address:bonito@tamu.edu (Vivette Girault) Sorbonne-Universit´e, CNRS, Universit´e de Paris, Laboratoire Jacques- Louis Lions (LJLL), F-75005 Paris, France Email address...

    work page 1957