pith. sign in

arxiv: 2604.06815 · v1 · submitted 2026-04-08 · 🧮 math.NA · cs.NA

A time-nonlocal multiphysics finite element method with Crank-Nicolson scheme for poroelasticity model with secondary consolidation

Zhihao Ge , Yanan He This is my paper

Pith reviewed 2026-05-10 18:18 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords poroelasticitysecondary consolidationfinite element methodCrank-Nicolson schemetime-nonlocalmultiphysicserror estimatesstability
0
0 comments X

The pith

A reformulated poroelasticity model with auxiliary variables yields a stable Crank-Nicolson finite element method with optimal error estimates.

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

The paper establishes that introducing auxiliary variables for fluid content and generalized pressure allows the original strongly coupled poroelasticity model with secondary consolidation to be rewritten as a generalized Stokes equation containing time-integral terms coupled to a diffusion equation. This reformulation exposes the multiphysics character and time-nonlocal features of the problem. The authors then construct a fully discrete scheme that uses Taylor-Hood mixed finite elements in space and the Crank-Nicolson method in time, while approximating the integral terms by the composite trapezoidal rule and maintaining running updates of the integrals. They prove existence and uniqueness of weak solutions, establish unconditional stability of the discrete scheme, and derive optimal-order error estimates in both space and time. Numerical experiments confirm the rates and illustrate that the Crank-Nicolson scheme exhibits superior long-time convergence behavior relative to backward Euler.

Core claim

By introducing the fluid content variable eta and the generalized pressure variable xi, the poroelasticity model is recast as a time-nonlocal multiphysics system consisting of a generalized Stokes problem with memory integrals and a diffusion equation. The fully discrete time-nonlocal multiphysics finite element method, which combines high-order Taylor-Hood elements with the Crank-Nicolson time scheme and trapezoidal quadrature for the integrals, is stable and attains optimal-order convergence rates.

What carries the argument

The reformulation that introduces auxiliary variables eta and xi to convert the original system into a generalized Stokes equation with time-integral terms plus a diffusion equation, which then admits a mixed finite-element discretization that treats the multiphysics coupling and the nonlocal memory explicitly.

If this is right

  • The scheme preserves second-order temporal accuracy while updating the integral terms in real time without repeated quadrature.
  • Long-time simulations of secondary consolidation become more accurate with Crank-Nicolson than with backward Euler, as shown by the numerical comparison.
  • The stability and error analysis rely on energy estimates and projection operators and hold under the positive-parameter assumption.
  • The method avoids storage of the full history of the integral terms and therefore remains computationally efficient for extended time intervals.

Where Pith is reading between the lines

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

  • The same auxiliary-variable reformulation could be applied to other poroelastic or viscoelastic models that contain secondary time-dependent effects.
  • The explicit time-nonlocal structure may help interpret physical mechanisms of long-term consolidation in geotechnical settings.
  • Efficiency gains from incremental integral updates could transfer to other PDEs that involve memory integrals discretized by similar quadrature rules.
  • Application to heterogeneous or nonlinear parameter regimes would test whether the stability proof extends beyond the constant-coefficient case.

Load-bearing premise

The physical parameters lambda, lambda-star, and c-zero must be finite positive constants so that the original model can be rewritten in the generalized Stokes-plus-diffusion form.

What would settle it

A numerical test on the reformulated model in which the computed solution becomes unstable for any choice of time step or in which the observed spatial or temporal convergence rate falls below the predicted optimal order would disprove the stability and error claims.

Figures

Figures reproduced from arXiv: 2604.06815 by Yanan He, Zhihao Ge.

Figure 4.1
Figure 4.1. Figure 4.1: Surface plots of displacement component u1 at the final time T. (a) Exact solution (b) Numerical solution [PITH_FULL_IMAGE:figures/full_fig_p029_4_1.png] view at source ↗
Figure 4.2
Figure 4.2. Figure 4.2: Surface plots of displacement component u2 at the final time T. (a) Exact solution (b) Numerical solution [PITH_FULL_IMAGE:figures/full_fig_p029_4_2.png] view at source ↗
Figure 4.3
Figure 4.3. Figure 4.3: Surface plots of pressure p at the final time T [PITH_FULL_IMAGE:figures/full_fig_p029_4_3.png] view at source ↗
read the original abstract

The paper studies a time-nonlocal multiphysics finite element method with Crank-Nicolson scheme for poroelasticity model with secondary consolidation. For the case where the physical parameters $\lambda,\lambda^*$ and $c_0$ are all finite positive constants, by introducing two auxiliary variables-the fluid content $\eta$ and the generalized pressure $\xi$ -- the original strongly coupled poroelasticity model is reformulated into a generalized Stokes equation with time integral terms and a diffusion equation. The reformulated model not only reveals the underlying multiphysics processes in the original model, but also exhibits time-nonlocal characteristics. A time-nonlocal multiphysics finite element method is designed for the reformulated model: the spatial discretization employs high order Taylor-Hood mixed finite element method, and the temporal discretization adopts the Crank-Nicolson scheme. The time integral terms are approximated using the composite trapezoidal rule, and the integral terms $J_{\xi}^n$ and $J_{\eta}^n$ are introduced for real-time updates, which not only avoids repeated calculations and improves efficiency, but also maintains second-order temporal accuracy. The existence and uniqueness of weak solutions for the reformulated model are proved via energy estimate methods, the stability of the fully discrete time-nonlocal multiphysics finite element method is established, and optimal-order error estimates are derived using projection operator techniques. Finally, numerical example verified the theoretical results and compared the long-time convergence of the Crank-Nicolson scheme and the backward Euler scheme.

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 studies a poroelasticity model with secondary consolidation. Under the assumption that the parameters λ, λ*, and c0 are finite positive constants, auxiliary variables η (fluid content) and ξ (generalized pressure) are introduced to reformulate the original strongly coupled system as a generalized Stokes equation containing time-integral terms together with a diffusion equation. A fully discrete scheme is constructed using high-order Taylor-Hood mixed finite elements in space, the Crank-Nicolson method in time, and composite trapezoidal quadrature for the integrals; auxiliary quantities J_ξ^n and J_η^n are maintained to avoid repeated integral evaluations while preserving second-order temporal accuracy. Existence and uniqueness of weak solutions for the reformulated model are proved by energy estimates, stability of the fully discrete scheme is established, and optimal-order error estimates are derived via projection operators. Numerical experiments confirm the theory and compare long-time behavior of the Crank-Nicolson and backward-Euler schemes.

Significance. If the analysis is correct, the work supplies a rigorous, second-order accurate numerical framework for time-nonlocal multiphysics poroelasticity that simultaneously reveals the underlying physical processes and improves computational efficiency. The paper credits standard but carefully applied techniques: energy-method proofs of well-posedness and stability, projection-based error analysis, and an auxiliary-variable device that retains accuracy without recomputing integrals. The inclusion of numerical verification of the theoretical rates is a further strength.

minor comments (3)
  1. [Section 3] The precise polynomial degrees employed for the Taylor-Hood elements (e.g., P_k-P_{k-1} with k=2 or higher) should be stated explicitly in the discretization section, together with the corresponding approximation properties used in the error analysis.
  2. [Section 2] A brief remark on the equivalence between the original and reformulated systems when the parameters are not all finite and positive would help delineate the scope of the analysis.
  3. [Section 5] Figure captions and axis labels in the numerical section could be expanded to indicate the specific norms or quantities plotted (displacement, pressure, fluid content) and the mesh sizes or time steps used.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive and accurate summary of our manuscript, the assessment of its significance, and the recommendation for minor revision. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No significant circularity; derivation uses standard analysis on reformulated system

full rationale

The paper begins with a parameter assumption (λ, λ*, c0 finite positive constants) that permits an exact algebraic reformulation of the original coupled system into a generalized Stokes problem with time integrals plus a diffusion equation, via the explicit introduction of auxiliary fields η (fluid content) and ξ (generalized pressure). Existence/uniqueness follows from standard energy estimates on this equivalent weak form; stability of the fully discrete scheme (Taylor-Hood mixed FEM + Crank-Nicolson + trapezoidal quadrature) is obtained by the same energy method applied to the discrete equations; optimal error bounds are derived via standard projection operators and consistency estimates. The auxiliary accumulators J_ξ^n and J_η^n are introduced solely for efficient real-time quadrature updates and do not alter the underlying continuous or discrete operators. No step reduces a claimed result to a fitted parameter, a self-referential definition, or a load-bearing self-citation; all load-bearing arguments are self-contained textbook techniques applied to an explicitly equivalent reformulation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 2 invented entities

The approach depends on assuming constant positive parameters to justify the reformulation and on introducing two auxiliary variables whose physical interpretation is tied to the original model.

axioms (1)
  • domain assumption The parameters λ, λ* and c0 are finite positive constants
    Explicitly stated as the case considered for the reformulation and analysis.
invented entities (2)
  • fluid content η no independent evidence
    purpose: Auxiliary variable to reformulate the strongly coupled system
    Introduced to obtain the generalized Stokes equation with time integrals.
  • generalized pressure ξ no independent evidence
    purpose: Auxiliary variable to reveal multiphysics and time-nonlocal features
    Used alongside η to decouple the original poroelasticity model.

pith-pipeline@v0.9.0 · 5587 in / 1411 out tokens · 62247 ms · 2026-05-10T18:18:34.708196+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

39 extracted references · 39 canonical work pages

  1. [1]

    Sekhdaria

    A. Sekhdaria. Irreversibility of soil skeleton deformations: the pedagogical limitations of terzaghi’s celebrated model for soil consolidation. Computers and Geotechnics , 137:104–137, 2021

    work page 2021

  2. [2]

    Józefiak, A

    K. Józefiak, A. Zbiciak, K. Brzozziński, and M. Maślakowski. A novel approach to the analysis of the soil consolidation problem by using non-classical rheological schemes. Applied Rheology, 26(1):137–156, 2016

    work page 2016

  3. [3]

    Gaspar, J.L

    F.J. Gaspar, J.L. Gracia, F.J. Lisbona, and P.N. Vabishchevich. A stabilized method for a secondary consolidation biot’s model. Numerical Methods for Partial Differential Equations , 24(1):60–78, 2010

    work page 2010

  4. [4]

    A new multiphysics finite element method for a biot model with secondary consolidation

    Zhihao Ge and Wenlong He. A new multiphysics finite element method for a biot model with secondary consolidation. CSIAM Transactions on Applied Mathematics , 5(3):515–550, 2024

    work page 2024

  5. [5]

    V. Mow, M. Kwan, W. Lai, and M. Holmes. A finite deformation theory for nonlinearly permeable soft hydrated biological tissues. In Biomechanics: Current Interdisciplinary Research , pages 153–179. Springer-Verlag, 1986. OPTIMAL ERROR ESTIMATES OF A NEW MFEM FOR NONLINEAR POROELASTICITY MODEL 33

    work page 1986

  6. [6]

    Y. C. Fung. Biomechanics: Mechanical Properties of Living Tissues . Springer-Verlag, New York, second edition edition, 1993

    work page 1993

  7. [7]

    Yang and L

    M. Yang and L. Taber. The possible role of poroelasticity in the apparent viscoelastic behavior of passive cardiac muscle. Journal of Biomechanics , 24(7):587–597, 1991

    work page 1991

  8. [8]

    Bociu, R

    L. Bociu, R. Guidoboni, R. Sacco, and J.T. Webster. Analysis of nonlinear poro-elastic and poro-visco-elastic models. Archive for Rational Mechanics and Analysis , 222(3):1445–1519, 2016

    work page 2016

  9. [9]

    W. Pao, R. Lewis, and I. Masters. A fully coupled hydro-thermo-poro-mechanical model for black oil reservoir simulation. International Journal for Numerical and Analytical Methods in Geomechanics , 25(13):1229–1256, 2001

    work page 2001

  10. [10]

    Poromechanics

    Olivier Coussy. Poromechanics. John Wiley & Sons, Chichester, 2004

    work page 2004

  11. [11]

    Rutqvist

    J. Rutqvist. The geomechanics of co2 storage in deep sedimentary formations. Geotechnical and Geological Engineering, 30(3):525–551, 2012

    work page 2012

  12. [12]

    Nelenne, P

    D. Nelenne, P. Ranga, M. Veble, P. Andersson, C. Tsang, and L. Jing. Coupled t-h-m issues related to radioactive waste repository design and performance. International Journal of Rock Mechanics and Mining Sciences , 143:104751, 2021

    work page 2021

  13. [13]

    Jacob Bear and Alexander H.-D. Cheng. Modeling Groundwater Flow and Contaminant Transport , volume 23 of Theory and Applications of Transport in Porous Media . Springer, Dordrecht, 2010

    work page 2010

  14. [14]

    A new mixed finite element method for a swelling clay model with secondary consolidation

    Wenlong Ge and Zhihao He. A new mixed finite element method for a swelling clay model with secondary consolidation. Applied Mathematical Modelling , 112:391–414, 2022

    work page 2022

  15. [15]

    Murad and John H

    Márcio A. Murad and John H. Cushman. Multiscale flow and deformation in hydrophilic swelling porous media. International Journal of Engineering Science , 34(3):313–338, 1996

    work page 1996

  16. [16]

    Showalter

    R.E. Showalter. Diffusion in poro-elastic media. Journal of Mathematical Analysis and Applications , 251(1):310–340, 2000

    work page 2000

  17. [17]

    Carleman estimate for biot consolidation system in poro-elasticity and appli- cation to inverse problems

    Mourad Bellassoued and Bochra Riahi. Carleman estimate for biot consolidation system in poro-elasticity and appli- cation to inverse problems. Mathematical Methods in the Applied Sciences , 39(18):5281–5301, 2016

    work page 2016

  18. [18]

    P. A. Vermeer and A. Verruijt. An accuracy condition for consolidation by finite elements. International Journal for Numerical and Analytical Methods in Geomechanics , 5(1):1–14, 1981

    work page 1981

  19. [19]

    O. C. Zienkiewicz and T. Shiomi. Static and dynamic behaviour of saturated porous media: The biot theory. Engi- neering Computations , 1(1):3–12, 1984

    work page 1984

  20. [20]

    M. A. Murad and A. F. D. Loula. On stability and convergence of finite element approximations of biot’s consolidation problem. International Journal for Numerical Methods in Engineering , 37(4):645–667, 1994

    work page 1994

  21. [21]

    F. J. Gaspar, F. J. Lisbona, and P. N. Vabishchevich. A finite difference analysis of biot’s consolidation model. Applied Numerical Mathematics , 44(4):487–506, 2003

    work page 2003

  22. [22]

    O. C. Zienkiewicz, A. H. C. Chan, M. Pastor, D. K. Paul, and T. Shiomi. Computational Geomechanics with Special Reference to Earthquake Engineering . John Wiley & Sons, Chichester, UK, 1999

    work page 1999

  23. [23]

    Lewis and Bernhard A

    Roland W. Lewis and Bernhard A. Schrefler. The Finite Element Method in the Static and Dynamic Deformation and Consolidation of Porous Media . John Wiley and Sons, Chichester, 1998

    work page 1998

  24. [24]

    Bause, F

    M. Bause, F. A. Radu, and U. K”ocher. Space–time finite element approximation of the biot poroelasticity system with iterative coupling. Computer Methods in Applied Mechanics and Engineering , 320:745–768, 2017

    work page 2017

  25. [25]

    Kunwar, H

    H. Kunwar, H. Lee, and K. Seelman. Second-order time discretization for a coupled quasi-newtonian fluid-poroelastic system. International Journal for Numerical Methods in Fluids , 92(6):687–702, 2020

    work page 2020

  26. [26]

    An arbitrary high-order discontinuous galerkin method for elastic waves on unstructured meshes - iv

    Josep de la Puente, Martin K”aser, Michael Dumbser, and Heiner Igel. An arbitrary high-order discontinuous galerkin method for elastic waves on unstructured meshes - iv. anisotropy. Geophysical Journal International, 169(3):1210– 1228, June 2007

    work page 2007

  27. [27]

    D. Yang, J. Peng, M. Lu, et al. A nearly analytical discrete method for wave-field simulations in 2d porous media. Communications in Computational Physics , 1(3):528–547, 2006

    work page 2006

  28. [28]

    Egger and M

    H. Egger and M. Sabouri. On the structure preserving high-order approximation of quasistatic poroelasticity. Math- ematics and Computers in Simulation , 189:237–252, 2021

    work page 2021

  29. [29]

    Stabilized multiphysics finite element method with Crank–Nicolson scheme for a poroelasticity model

    Zhihao Ge, Yanan He, and Tingting Li. Stabilized multiphysics finite element method with Crank–Nicolson scheme for a poroelasticity model. Numerical Methods for Partial Differential Equations , 35(4):1412–1428, 2019

    work page 2019

  30. [30]

    Crank and P

    J. Crank and P. Nicolson. A practical method for numerical evaluation of solutions of partial differential equations of the heat–conduction type. Mathematical Proceedings of the Cambridge Philosophical Society , 43(1):50–67, 1947

    work page 1947

  31. [31]

    Lawrence C. Evans. Partial Differential Equations . American Mathematical Society, Providence, RI, second edition edition, 2016

    work page 2016

  32. [32]

    Analysis of a multiphysics finite element method for a poroelasticity model

    Xiaobing Feng, Zhihao Ge, and Yukun Li. Analysis of a multiphysics finite element method for a poroelasticity model. IMA Journal of Numerical Analysis , 38(1):330–359, 2018

    work page 2018

  33. [33]

    Phillip Joseph Phillips and Mary F. Wheeler. A coupling of mixed and continuous galerkin finite element methods for poroelasticity i: the continuous in time case. Computational Geosciences, 11(2):131–144, 2007

    work page 2007

  34. [34]

    Phillip Joseph Phillips and Mary F. Wheeler. A coupling of mixed and continuous galerkin finite element methods for poroelasticity ii: the discrete in time case. Computational Geosciences, 11(2):145–158, 2007

    work page 2007

  35. [35]

    A lowest equal-order stabilized mixed finite element method based on multiphysics approach for a poroelasticity model

    Zhihao Ge, Yanan He, and Yinnian He. A lowest equal-order stabilized mixed finite element method based on multiphysics approach for a poroelasticity model. Applied Numerical Mathematics , 153:1–14, 2020

    work page 2020

  36. [36]

    Brenner and L.R

    S.C. Brenner and L.R. Scott. The Mathematical Theory of Finite Element Methods . Springer-Verlag New York Inc, third edition edition, 2008

    work page 2008

  37. [37]

    Navier–Stokes Equations: Theory and Numerical Analysis

    Roger Temam. Navier–Stokes Equations: Theory and Numerical Analysis . North-Holland, 1984

    work page 1984

  38. [38]

    Finite Element Method for Navier-Stokes Equations: theory and algo- 34 YANAN HE, ZHIHAO GE rithms

    Vivette Girault and Pierre-Arnaud Raviart. Finite Element Method for Navier-Stokes Equations: theory and algo- 34 YANAN HE, ZHIHAO GE rithms. Springer-Verlag, Berlin, 1986

    work page 1986

  39. [39]

    Mixed and Hybrid Finite Element Methods

    Franco Brezzi and Michel Fortin. Mixed and Hybrid Finite Element Methods . Springer-Verlag, New York, 1991

    work page 1991