Pith. sign in

REVIEW 3 minor 71 references

Linearizing nonlinear poroelasticity produces a consistent model with incompressible solid phase that supports lowest-order inf-sup stable finite elements.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.3

2026-06-27 23:53 UTC pith:KV4EJTZO

load-bearing objection Linearizing nonlinear poroelasticity gives a consistent solid-incompressible model with lowest-order inf-sup stable elements and Schur solvers.

arxiv 2606.06750 v1 pith:KV4EJTZO submitted 2026-06-04 math.NA cs.NA

Linear poroelasticity with solid incompressibility: consistent formulation and scalable numerical solution

classification math.NA cs.NA
keywords poroelasticitysolid incompressibilityfinite element discretizationSchur complementlinearizationinf-sup stabilitynumerical solution
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper linearizes the equations of fully nonlinear poroelasticity to obtain a model in which only the solid phase is incompressible. This reformulation avoids inconsistency issues that arise in standard primal formulations while preserving the essential fluid-solid coupling mechanisms. It establishes a well-posed discretization using lowest-order inf-sup stable finite element spaces and develops scalable solvers via a Schur complement reduction. Numerical experiments in two and three dimensions are used to confirm the approach works as intended.

Core claim

By linearizing the equations of fully nonlinear poroelasticity, a consistent model in which only the solid phase is incompressible is proposed. This reformulation circumvents some inconsistency issues encountered in standard primal formulations of nonlinear poroelasticity while still retaining its key physical coupling mechanisms. We show a well-posed and consistent discretization strategy and also formulate scalable solvers based on a Schur complement formalism. A distinctive feature of the model is that it allows for a lowest order, inf-sup stable family of Finite Elements spaces.

What carries the argument

The linearized poroelasticity model with solid incompressibility, which permits a well-posed discretization by lowest-order inf-sup stable finite element spaces together with Schur complement solvers.

Load-bearing premise

Linearizing the fully nonlinear poroelasticity equations produces a model that remains consistent, retains the key physical couplings, and admits a well-posed discretization with lowest-order inf-sup stable finite element spaces.

What would settle it

A manufactured solution test in which the chosen lowest-order finite element spaces fail to satisfy the discrete inf-sup condition for the linearized system would falsify the well-posedness claim.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The model retains the key physical coupling mechanisms between fluid and solid phases.
  • Lowest-order inf-sup stable finite element spaces can be employed for the discretization.
  • Scalable solvers can be constructed using a Schur complement formalism.
  • Numerical tests in two and three dimensions confirm stability and consistency of the discrete solutions.

Where Pith is reading between the lines

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

  • The same linearization strategy might be applied to other nonlinear coupled systems where incompressibility leads to similar inconsistencies.
  • The retained couplings suggest the model could serve as an intermediate step toward consistent nonlinear solvers for large-deformation poroelasticity.
  • Stability of the lowest-order spaces may extend to related mixed problems involving volume constraints on one phase only.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

0 major / 3 minor

Summary. The manuscript proposes a linearized model of poroelasticity in which only the solid phase is incompressible, obtained by linearization of the fully nonlinear equations. The reformulation is claimed to avoid inconsistencies of standard primal nonlinear formulations while preserving key physical couplings. The authors present a well-posed and consistent discretization that admits a lowest-order inf-sup stable family of finite-element spaces, formulate scalable Schur-complement solvers, and report numerical validation in two and three dimensions.

Significance. If the linearization, well-posedness, and consistency results hold, the work supplies a practical route to stable low-order discretizations and scalable solvers for a class of poroelastic problems that arise in biomechanics and geomechanics. The explicit retention of solid incompressibility together with the inf-sup stable lowest-order elements constitutes a concrete computational advantage over many existing primal formulations.

minor comments (3)
  1. [§3] §3 (discretization): the statement that the chosen spaces are 'inf-sup stable' for the linearized system should be accompanied by an explicit reference to the theorem or lemma that establishes the discrete inf-sup condition, rather than only citing the continuous well-posedness result.
  2. [Table 1, Figure 4] Table 1 and Figure 4: the reported convergence rates for the 3-D test are given only for the displacement and pressure fields; inclusion of the fluid flux or Darcy velocity error would strengthen the claim that the full coupled system is consistently approximated.
  3. [§4.2] §4.2 (solver): the Schur-complement preconditioner is described at a high level; a brief statement of the spectral equivalence or eigenvalue bounds used to justify robustness with respect to the permeability parameter would clarify the scalability claim.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript, the recognition of its potential significance in biomechanics and geomechanics, and the recommendation for minor revision. No specific major comments are listed in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The central derivation proceeds by direct linearization of the fully nonlinear poroelasticity equations to obtain a model with solid-phase incompressibility. This is a standard first-principles approximation step whose output is not equivalent to its inputs by construction. Well-posedness of the discretization, inf-sup stability for lowest-order elements, and Schur-complement solvers are established as independent mathematical results. No fitted parameters are renamed as predictions, no self-citation chains carry the load-bearing claims, and no ansatz is smuggled via prior work. Numerical tests function as external validation rather than tautological confirmation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that linearization preserves consistency and key couplings, plus the existence of a well-posed lowest-order inf-sup stable discretization; no free parameters or invented entities are indicated in the abstract.

axioms (1)
  • domain assumption Linearization of fully nonlinear poroelasticity equations yields a consistent model with solid incompressibility that retains key physical coupling mechanisms
    This premise is invoked as the foundation for the proposed reformulation in the abstract.

pith-pipeline@v0.9.1-grok · 5647 in / 1271 out tokens · 28682 ms · 2026-06-27T23:53:01.453277+00:00 · methodology

0 comments
read the original abstract

In this work we propose, by linearizing the equations of fully nonlinear poroelasticity, a consistent model in which only the solid phase is incompressible. This reformulation circumvents some inconsistency issues encountered in standard primal formulations of nonlinear poroelasticity while still retaining its key physical coupling mechanisms. We show a well-posed and consistent discretization strategy and also formulate scalable solvers based on a Schur complement formalism. A distinctive feature of the model is that it allows for a lowest order, inf-sup stable family of Finite Elements (FE) spaces. Numerical tests in two and three dimensions are provided to validate the proposed method and solver framework.

Figures

Figures reproduced from arXiv: 2606.06750 by Andr\'es E. Rubiano, Nicol\'as A. Barnafi, Ricardo Ruiz-Baier.

Figure 4.1
Figure 4.1. Figure 4.1: Comparison of the pressure fields in the steady-state Mandel test. The solutions obtained from [PITH_FULL_IMAGE:figures/full_fig_p012_4_1.png] view at source ↗
Figure 4.2
Figure 4.2. Figure 4.2: Approximate solutions for the convergence test in the steady case, computed with the second [PITH_FULL_IMAGE:figures/full_fig_p013_4_2.png] view at source ↗
Figure 4.3
Figure 4.3. Figure 4.3: Approximate solutions for the convergence test in the time-dependent case computed with the [PITH_FULL_IMAGE:figures/full_fig_p014_4_3.png] view at source ↗

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

71 extracted references · 3 canonical work pages

  1. [1]

    2004 , publisher=

    Coussy, O , title=. 2004 , publisher=

  2. [2]

    Bansal, A and Barnafi, N A and Pandey, D N and Ruiz-Baier, R , journal=. A. 2026 , publisher=

  3. [3]

    Calcolo , volume=

    Error estimates of a new mixed finite element method for arbitrary element pair for a nonlinear thermo-poroelasticity model , author=. Calcolo , volume=. 2026 , publisher=

  4. [4]

    SIAM Journal on Applied Mathematics , volume =

    Barnafi Wittwer, NA and Di Gregorio, S and Dede', L and Zunino, P and Vergara, C and Quarteroni, AM , title =. SIAM Journal on Applied Mathematics , volume =. 2022 , doi =

  5. [5]

    Journal of Scientific Computing , volume=

    Finite Element Methods for Large-Strain Poroelasticity/Chemotaxis Models Simulating the Formation of Myocardial Oedema , author=. Journal of Scientific Computing , volume=. 2022 , publisher=

  6. [6]

    Ruiz-Baier, R and Taffetani, M and Westermeyer, HD and Yotov, I , journal=. The. 2022 , publisher=

  7. [7]

    Journal of the Mechanics and Physics of Solids , volume=

    A transversely isotropic thermo-hyperelastic constitutive model of myocardial tissue with a three-state cell death dynamics for cardiac radiofrequency ablation , author=. Journal of the Mechanics and Physics of Solids , volume=. 2022 , publisher=

  8. [8]

    International Journal for Numerical Methods in Biomedical Engineering , volume=

    A computational model of open-irrigated radiofrequency catheter ablation accounting for mechanical properties of the cardiac tissue , author=. International Journal for Numerical Methods in Biomedical Engineering , volume=. 2019 , publisher=

  9. [9]

    World Journal of Surgery , volume=

    Pathogenesis of edema formation in burn injuries , author=. World Journal of Surgery , volume=. 1992 , publisher=

  10. [10]

    European Journal of Mechanics, B/Fluids , year=

    Chapelle, D and Moireau, P , title=. European Journal of Mechanics, B/Fluids , year=

  11. [11]

    Communications in Computational Physics , year =

    Ruiz-Baier, R and Gizzi, A and Loppini, A Loppini and Cherubini, C and Filippi, S , title =. Communications in Computational Physics , year =. doi:https://doi.org/10.4208/cicp.OA-2018-0253 , url =

  12. [12]

    Journal of Elasticity , keywords =

    Effect of laminar orthotropic myofiber architecture on regional stress and strain in the canine left ventricle , volume =. Journal of Elasticity , keywords =. 2000 , author=. doi:10.1023/A:1010883920374 , issn =

  13. [13]

    Computing and Visualization in Science , volume=

    Computational model of three-dimensional cardiac electromechanics , author=. Computing and Visualization in Science , volume=. 2002 , publisher=

  14. [14]

    Frontiers in physiology , volume=

    Multiscale and Multiphysics Modeling of Anisotropic Cardiac RFCA: Experimental-Based Model Calibration via Multi-Point Temperature Measurements , author=. Frontiers in physiology , volume=. 2022 , publisher=

  15. [15]

    Journal of Interventional Cardiac Electrophysiology , pages=

    Spatial temperature reconstructions in myocardial tissues undergoing radiofrequency ablations by performing high-resolved temperature measurements , author=. Journal of Interventional Cardiac Electrophysiology , pages=. 2022 , publisher=

  16. [16]

    ACM Transactions on Mathematical Software (TOMS) , volume=

    Firedrake: automating the finite element method by composing abstractions , author=. ACM Transactions on Mathematical Software (TOMS) , volume=. 2016 , publisher=

  17. [17]

    2015 , publisher=

    The paraview guide: a parallel visualization application , author=. 2015 , publisher=

  18. [18]

    2007 , publisher=

    The porous medium equation: mathematical theory , author=. 2007 , publisher=

  19. [19]

    Springer Science , volume=

    Numerical optimization , author=. Springer Science , volume=

  20. [20]

    2000 , organization=

    Amestoy, PR and Duff, IS and L’Excellent, J-Y and Koster, J , booktitle=. 2000 , organization=

  21. [21]

    SIAM Journal on Numerical Analysis , volume=

    Anderson acceleration for fixed-point iterations , author=. SIAM Journal on Numerical Analysis , volume=. 2011 , publisher=

  22. [22]

    Journal of the ACM (JACM) , volume=

    Iterative procedures for nonlinear integral equations , author=. Journal of the ACM (JACM) , volume=. 1965 , publisher=

  23. [23]

    Computer Methods in Applied Mechanics and Engineering , volume=

    Computational methods for inverse finite elastostatics , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 1996 , publisher=

  24. [24]

    Basic research in cardiology , volume=

    Model of the coronary circulation based on pressure dependence of coronary resistance and compliance , author=. Basic research in cardiology , volume=. 1988 , publisher=

  25. [25]

    Journal of biomechanics , volume=

    A novel porous mechanical framework for modelling the interaction between coronary perfusion and myocardial mechanics , author=. Journal of biomechanics , volume=. 2012 , publisher=

  26. [26]

    Biomechanics and Modeling in Mechanobiology , volume=

    In silico coronary wave intensity analysis: application of an integrated one-dimensional and poromechanical model of cardiac perfusion , author=. Biomechanics and Modeling in Mechanobiology , volume=. 2016 , publisher=

  27. [27]

    1996 , publisher=

    Solving ordinary differential equations II , author=. 1996 , publisher=

  28. [28]

    Nonstiff problems , author=

    Solving ordinary differential equations I. Nonstiff problems , author=. 1993 , publisher=

  29. [29]

    Mathematics , volume=

    Coupling chemotaxis and growth poromechanics for the modelling of feather primordia patterning , author=. Mathematics , volume=. 2022 , publisher=

  30. [30]

    Physiology masterclass: Extremes of age: newborn and infancy , volume =

    LoMauro, A and Aliverti, A , doi =. Physiology masterclass: Extremes of age: newborn and infancy , volume =. Breathe , number =

  31. [31]

    Biomechanics and Modeling in Mechanobiology , volume=

    A quasi-static poromechanical model of the lungs , author=. Biomechanics and Modeling in Mechanobiology , volume=. 2022 , publisher=

  32. [32]

    Journal of Biomechanical Engineering , volume =

    Guccione, JM and McCulloch, AD and Waldman, LK , title = ". Journal of Biomechanical Engineering , volume =. 1991 , issn =

  33. [33]

    Biomechanics and modeling in mechanobiology , volume=

    The importance of the pericardium for cardiac biomechanics: from physiology to computational modeling , author=. Biomechanics and modeling in mechanobiology , volume=. 2019 , publisher=

  34. [34]

    arXiv preprint arXiv:2305.01199 , year=

    Modeling of cardiac fibers as oriented liquid crystals , author=. arXiv preprint arXiv:2305.01199 , year=

  35. [35]

    American Journal of Physiology-Heart and Circulatory Physiology , volume=

    Porous medium finite element model of the beating left ventricle , author=. American Journal of Physiology-Heart and Circulatory Physiology , volume=

  36. [36]

    Computational Mechanics , volume=

    A poroelastic model valid in large strains with applications to perfusion in cardiac modeling , author=. Computational Mechanics , volume=. 2010 , publisher=

  37. [37]

    Journal of biomechanics , volume=

    Prestressing in finite deformation abdominal aortic aneurysm simulation , author=. Journal of biomechanics , volume=. 2009 , publisher=

  38. [38]

    Journal of biomechanics , volume=

    An augmented iterative method for identifying a stress-free reference configuration in image-based biomechanical modeling , author=. Journal of biomechanics , volume=. 2017 , publisher=

  39. [39]

    Theory and Applications

    A simple introduction to the mixed finite element method , author=. Theory and Applications. Springer Briefs in Mathematics. Springer, London , year=

  40. [40]

    Tissue poromechanical deformation effects on steam pop likelihood in

    Wongchadakul, P and Datta, AK and Rattanadecho, P , journal=. Tissue poromechanical deformation effects on steam pop likelihood in. 2023 , publisher=

  41. [41]

    1987 , publisher=

    Theory and applications of transport in porous media , author=. 1987 , publisher=

  42. [42]

    American Journal of Physiology-Heart and Circulatory Physiology , volume=

    Mechanisms of myocardium-coronary vessel interaction , author=. American Journal of Physiology-Heart and Circulatory Physiology , volume=. 2010 , publisher=

  43. [43]

    Journal of the mechanical behavior of biomedical materials , volume=

    A modular inverse elastostatics approach to resolve the pressure-induced stress state for in vivo imaging based cardiovascular modeling , author=. Journal of the mechanical behavior of biomedical materials , volume=. 2018 , publisher=

  44. [44]

    Computers in Biology and Medicine , volume=

    Mechanoelectric effects in healthy cardiac function and under Left Bundle Branch Block pathology , author=. Computers in Biology and Medicine , volume=. 2023 , publisher=

  45. [45]

    Computer Methods in Applied Mechanics and Engineering , volume=

    Fully nonlinear inverse poroelasticity: Stress-free configuration recovery , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2024 , publisher=

  46. [46]

    Computer Methods in Applied Mechanics and Engineering , volume=

    Reconstructing relaxed configurations in elastic bodies: Mathematical formulations and numerical methods for cardiac modeling , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2024 , publisher=

  47. [47]

    Journal of biomechanics , volume=

    Heterogeneous growth-induced prestrain in the heart , author=. Journal of biomechanics , volume=. 2015 , publisher=

  48. [48]

    Journal of Biomechanical Engineering , volume=

    Estimation of regional pulmonary compliance in idiopathic pulmonary fibrosis based on personalized lung poromechanical modeling , author=. Journal of Biomechanical Engineering , volume=. 2022 , publisher=

  49. [49]

    Biomechanics and Modeling in Mechanobiology , pages=

    Comparison of optimization parametrizations for regional lung compliance estimation using personalized pulmonary poromechanical modeling , author=. Biomechanics and Modeling in Mechanobiology , pages=. 2023 , publisher=

  50. [50]

    SIAM journal on control and optimization , volume=

    Monotone operators and the proximal point algorithm , author=. SIAM journal on control and optimization , volume=. 1976 , publisher=

  51. [51]

    2003 , publisher=

    Numerical methods for image registration , author=. 2003 , publisher=

  52. [52]

    SIAM Journal on Imaging Sciences , volume=

    Primal and mixed finite element methods for deformable image registration problems , author=. SIAM Journal on Imaging Sciences , volume=. 2018 , publisher=

  53. [53]

    SIAM Journal on Imaging Sciences , volume=

    Adaptive mesh refinement in deformable image registration: a posteriori error estimates for primal and mixed formulations , author=. SIAM Journal on Imaging Sciences , volume=. 2021 , publisher=

  54. [54]

    Mathematical Models and Methods in Applied Sciences , volume=

    New primal and dual-mixed finite element methods for stable image registration with singular regularization , author=. Mathematical Models and Methods in Applied Sciences , volume=. 2021 , publisher=

  55. [55]

    Parameter-robust methods for the

    Boon, WM and Hornkj. Parameter-robust methods for the. Journal of Computational Physics , volume=. 2022 , publisher=

  56. [56]

    Hybridizable discontinuous

    Cesmelioglu, A and Lee, JJ and Rhebergen, S , journal=. Hybridizable discontinuous. 2023 , publisher=

  57. [57]

    A numerical solution of the

    Taylor, C and Hood, P , journal=. A numerical solution of the. 1973 , publisher=

  58. [58]

    Spe Journal , volume=

    A coupled reservoir and geomechanical simulation system , author=. Spe Journal , volume=. 1998 , publisher=

  59. [59]

    Robust fixed stress splitting for

    Both, JW and Borregales, M and Nordbotten, JM and Kumar, K and Radu, FA , journal=. Robust fixed stress splitting for. 2017 , publisher=

  60. [60]

    Computer Methods in Applied Mechanics and Engineering , volume=

    Iterative splitting schemes for a soft material poromechanics model , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2022 , publisher=

  61. [61]

    Computer Methods in Applied Mechanics and Engineering , volume=

    Block-partitioned solvers for coupled poromechanics: A unified framework , author=. Computer Methods in Applied Mechanics and Engineering , volume=. 2016 , publisher=

  62. [62]

    Journal of Computational Physics , volume=

    Scalable algorithms for three-field mixed finite element coupled poromechanics , author=. Journal of Computational Physics , volume=. 2016 , publisher=

  63. [63]

    Numerical solution of partial differential equations on parallel computers , pages=

    The design and implementation of hypre, a library of parallel high performance preconditioners , author=. Numerical solution of partial differential equations on parallel computers , pages=. 2006 , organization=

  64. [64]

    Balay, S and Abhyankar, S and Adams, M and Brown, J and Brune, P and Buschelman, K and Dalcin, L and Dener, A and Eijkhout, V and Gropp, W , year=

  65. [65]

    Graham, RL and Woodall, TS and Squyres, JM , booktitle=. Open. 2006 , organization=

  66. [66]

    Numerical solution of the

    Kumar, S and Mora, D and Ruiz-Baier, R and Verma, N , journal=. Numerical solution of the. 2024 , publisher=

  67. [67]

    The fixed-stress splitting scheme for

    Storvik, E and Both, JW and Nordbotten, JM and Radu, FA , booktitle=. The fixed-stress splitting scheme for. 2020 , organization=

  68. [68]

    PAMM , volume=

    Benchmark computations of dynamic poroelasticity , author=. PAMM , volume=. 2023 , publisher=

  69. [69]

    SIAM Journal on Matrix Analysis and Applications , volume=

    Reducing complexity in parallel algebraic multigrid preconditioners , author=. SIAM Journal on Matrix Analysis and Applications , volume=. 2006 , publisher=

  70. [70]

    Numerical Linear Algebra with Applications , volume=

    Distance-two interpolation for parallel algebraic multigrid , author=. Numerical Linear Algebra with Applications , volume=. 2008 , publisher=

  71. [71]

    2013 , doi =

    Boffi, D and Brezzi, F and Fortin, M , title =. 2013 , doi =