A Locking-free and Loosely Coupled Robin-Robin Scheme for Fluid-Poroelasticity Interaction

Haibiao Zheng; Jiwei Zhang; Thomas Wick; Wenlong He; Xiaohe Yue

arxiv: 2604.07033 · v2 · submitted 2026-04-08 · 🧮 math.NA · cs.NA

A Locking-free and Loosely Coupled Robin-Robin Scheme for Fluid-Poroelasticity Interaction

Wenlong He , Thomas Wick , Xiaohe Yue , Jiwei Zhang , Haibiao Zheng This is my paper

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

classification 🧮 math.NA cs.NA

keywords fluid-poroelasticity interactionRobin-Robin schemelocking-freedecoupled schemeBiot systemStokes equationsunconditional stability

0 comments

The pith

A reformulation of the Biot system into four fields enables a fully decoupled Robin-Robin scheme for fluid-poroelasticity interaction that is unconditionally stable and locking-free.

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

The paper develops a partitioned numerical method for simulating the coupled dynamics of an incompressible fluid and a poroelastic solid, as occurs in biomechanics and geomechanics. Standard coupled approaches often require sub-iterations at each time step or suffer from locking when poroelastic parameters reach extreme values such as very low permeability. The authors introduce two auxiliary variables to rewrite the Biot equations as a four-field system whose interface conditions match those of the original problem exactly. They then impose lagged Robin-Robin transmission conditions so that the fluid and poroelastic subproblems become completely independent and can be solved in parallel. The resulting scheme is proved to be unconditionally stable and to deliver optimal-order error estimates in the H1 norm that remain uniform across all parameter regimes.

Core claim

The central claim is that the four-field reformulation of the fully dynamic Biot system, paired with explicitly lagged Robin-Robin interface data, produces a fully decoupled, non-iterative time-stepping scheme for the unsteady Stokes-Biot coupling. This scheme preserves the physical interface conditions without consistency error, admits unconditional stability, and yields optimal H1 error bounds that are robust with respect to all poroelastic parameters, thereby eliminating the locking phenomena that appear in conventional formulations.

What carries the argument

The four-field reformulation of the Biot system together with lagged Robin-Robin transmission conditions, which decouples the fluid and poroelastic subproblems while preserving the original interface coupling exactly.

Load-bearing premise

The reformulation of the Biot system into four fields preserves the original interface conditions exactly, so that the lagged Robin-Robin conditions introduce no additional consistency error.

What would settle it

A sequence of computations on successively refined meshes with permeability approaching zero or Poisson ratio approaching one-half, checking whether the observed H1 error rates remain optimal and whether the discrete solution remains free of spurious oscillations.

Figures

Figures reproduced from arXiv: 2604.07033 by Haibiao Zheng, Jiwei Zhang, Thomas Wick, Wenlong He, Xiaohe Yue.

**Figure 2.** Figure 2: The comparison of numerical results obtained by different algorithms [PITH_FULL_IMAGE:figures/full_fig_p027_2.png] view at source ↗

**Figure 3.** Figure 3: Numerical results for up,y, vf,x and pf on the fluid-poroelastic interface and vf,x on the bottom fluid boundary obtained from Algorithm 1. AL 1.1, AL 1.2, and AL 1.3 refer to the parameter configurations in [PITH_FULL_IMAGE:figures/full_fig_p030_3.png] view at source ↗

**Figure 4.** Figure 4: Surface plots of pressure and fluid velocity for Algorithm 1 with Case [PITH_FULL_IMAGE:figures/full_fig_p031_4.png] view at source ↗

read the original abstract

We study a fluid-poroelasticity interaction (FPSI) problem coupling the unsteady Stokes equations with the fully dynamic Biot system. A major challenge in such problems is to design partitioned schemes that remain robust in locking-related parameter regimes while preserving the physical interface coupling structure.To address this issue, we introduce two auxiliary variables and reformulate the Biot system as a four-field problem consisting of a dynamic Stokes-like system coupled with a diffusion equation. Crucially, this reformulation preserves the original interface conditions. Based on Robin-Robin transmission conditions with explicitly lagged interface data, we construct a fully decoupled scheme in which the fluid and poroelastic subproblems can be solved independently and in parallel at each time step, without sub-iterations.We prove that the resulting method is unconditionally stable and derive optimal-order error estimates in the $H^1$-norm. The analysis further shows that the scheme is robust with respect to extreme poroelastic parameters and avoids the locking effects inherent in standard formulations. Numerical experiments confirm the theoretical convergence results and demonstrate the locking-robust performance of the proposed method.

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.

Desk Editor's Note private letter to a colleague

The paper gives a decoupled Robin-Robin scheme for fluid-poroelasticity via a four-field Biot reformulation that claims unconditional stability and locking-free behavior.

read the letter

This paper introduces a concrete way to split fluid-poroelasticity interaction into fully independent parallel solves using lagged Robin-Robin conditions. They add two auxiliary variables to recast the Biot system as a dynamic Stokes-like problem plus a diffusion equation, then apply the transmission conditions with explicit lagging so neither subproblem needs the other at the current step. The abstract states that this reformulation keeps the original interface conditions intact, which underpins the unconditional stability result and the optimal-order H1 error estimates. The analysis also targets robustness when permeability vanishes or Lamé parameters become large, and the numerics are said to confirm both the rates and the absence of locking. That combination of a new reformulation, a stability proof, and parameter-robust numerics is the actual advance here. The soft spot is exactly the one the stress test flags: whether the four-field version matches the original interface conditions exactly or only weakly. If any regularity or approximation sneaks in at the interface, the lagging step could inject consistency errors that the energy argument does not control, and those errors would be largest in the locking regimes the method is meant to handle. The paper asserts exact preservation, so the derivation details matter. Standard energy methods are used for the analysis, which is appropriate if the assumptions hold, and the citations appear to build directly on prior partitioned schemes without circularity. This is aimed at people who need practical, parallelizable schemes for FPSI in biomechanics or porous-media applications. It has enough new construction and supporting analysis to deserve a serious referee, provided the reformulation equivalence is checked carefully in review.

Referee Report

3 major / 2 minor

Summary. The manuscript proposes a partitioned time-stepping scheme for unsteady Stokes-Biot fluid-poroelasticity interaction. Two auxiliary variables are introduced to recast the Biot system as an equivalent four-field problem (dynamic Stokes-like fluid plus diffusion), after which lagged Robin-Robin transmission conditions are applied to obtain a fully decoupled, parallelizable scheme. The central claims are unconditional stability, optimal-order H¹ error estimates, and robustness to extreme poroelastic parameters (vanishing permeability, large Lamé coefficients) without locking.

Significance. If the stability and error analysis hold with parameter-independent constants, the work would provide a practical advance in the design of loosely coupled schemes for FPSI problems. The combination of exact interface preservation, unconditional stability, and locking-free behavior could reduce the need for sub-iterations or monolithic solves in applications such as biomechanics and reservoir simulation.

major comments (3)

[§3] §3 (Four-field reformulation): The assertion that the auxiliary-variable reformulation preserves the original fluid-poroelastic interface conditions exactly (paragraph following the definition of the auxiliary fields) is load-bearing for the subsequent analysis. The lagged Robin-Robin scheme then inherits these conditions without additional consistency error. The manuscript must supply the explicit verification that no extra interface regularity or weak-form modification is required, otherwise the O(Δt) lagging term can produce hidden consistency errors that grow precisely when permeability κ→0 or λ,μ→∞.
[Theorem 4.1] Theorem 4.1 (Unconditional stability): The energy estimate establishing unconditional stability must be checked for uniformity with respect to the poroelastic parameters. If any constant in the estimate depends on 1/κ or on λ, the robustness claim in the locking regimes is not supported, even if the scheme remains stable for fixed parameters.
[§5] §5 (H¹ error analysis): The optimal-order H¹ estimate is stated to hold uniformly. The proof should explicitly track the dependence of the error constant on Δt, h, κ, λ, and μ; any hidden dependence would contradict the locking-free assertion and must be removed or bounded independently of these parameters.

minor comments (2)

[Abstract] The abstract claims 'optimal-order error estimates in the H¹-norm' but does not state the precise rate (e.g., O(Δt + h)). Adding the rate would improve clarity.
[§4] Notation for the auxiliary variables and the lagged Robin data should be introduced once and used consistently; occasional re-definition in the analysis section creates minor confusion.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The suggestions help clarify the key aspects of the four-field reformulation, stability, and error analysis. We address each major comment below and have revised the manuscript to incorporate explicit verifications and remarks on parameter independence.

read point-by-point responses

Referee: [§3] The assertion that the auxiliary-variable reformulation preserves the original fluid-poroelastic interface conditions exactly is load-bearing... The manuscript must supply the explicit verification that no extra interface regularity or weak-form modification is required, otherwise the O(Δt) lagging term can produce hidden consistency errors that grow precisely when permeability κ→0 or λ,μ→∞.

Authors: We agree that an explicit verification strengthens the presentation. In the revised Section 3, we have added Lemma 3.1, which derives the weak form of the four-field system directly from the original Stokes-Biot interface conditions. The proof shows equivalence in the natural trace spaces without requiring additional regularity beyond the standard H^{1/2} interface traces. The auxiliary variables are constructed so that the transmission conditions (continuity of velocity and normal stress) hold identically. The O(Δt) lagging is treated as a consistency term in the error analysis of Section 5, where it is bounded independently of κ, λ, and μ using the stability of the Robin-Robin operators. revision: yes
Referee: [Theorem 4.1] The energy estimate establishing unconditional stability must be checked for uniformity with respect to the poroelastic parameters. If any constant in the estimate depends on 1/κ or on λ, the robustness claim in the locking regimes is not supported.

Authors: The proof of Theorem 4.1 already yields a stability bound whose constants are independent of κ, λ, and μ. This follows from testing the fluid and poroelastic subproblems with the velocity and displacement fields, respectively, and absorbing the interface Robin terms without introducing inverse powers of κ or λ. The four-field structure allows the poroelastic dissipation and elastic energy to be controlled uniformly. In the revised manuscript we have inserted a remark immediately after the theorem statement that explicitly lists the parameter-independent constants and sketches the key absorption steps. revision: yes
Referee: [§5] The optimal-order H¹ estimate is stated to hold uniformly. The proof should explicitly track the dependence of the error constant on Δt, h, κ, λ, and μ; any hidden dependence would contradict the locking-free assertion.

Authors: We have expanded the error analysis in Section 5 to include an explicit tracking of all constants. The proof proceeds by subtracting the continuous and discrete weak forms, applying the stability result from Theorem 4.1, and estimating the consistency terms arising from time lagging and spatial discretization. Each term is bounded using inverse inequalities and the parameter-robust coercivity of the Robin-Robin operators; no factors of 1/κ or λ appear in the final O(Δt + h) constant. A new corollary (Corollary 5.1) states the uniform bound and lists the explicit dependence (or independence) on each parameter. revision: yes

Circularity Check

0 steps flagged

No circularity: standard energy-method analysis of a reformulated partitioned scheme

full rationale

The paper's central claims rest on introducing auxiliary variables to obtain a four-field reformulation of the Biot system, then applying lagged Robin-Robin transmission conditions to produce a decoupled time-stepping scheme. Unconditional stability and optimal H1 error estimates are asserted to follow from energy estimates that exploit the exact preservation of interface conditions. No quoted equation or step reduces a derived quantity to a fitted parameter, a self-citation chain, or a definition that already encodes the target result. The reformulation is presented as an algebraic identity that leaves the original interface conditions unchanged; the subsequent stability proof is therefore independent of the result it establishes. This is the normal, non-circular structure for a numerical-analysis paper whose proofs rely on standard a-priori estimates rather than on any of the enumerated circular patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the mathematical reformulation preserving interface conditions and standard stability analysis techniques for decoupled schemes.

axioms (1)

domain assumption Standard assumptions on solution regularity for the unsteady Stokes and fully dynamic Biot equations
Required to obtain the stated H1 error estimates and stability bounds in numerical PDE analysis.

pith-pipeline@v0.9.0 · 5501 in / 1166 out tokens · 28903 ms · 2026-05-10T17:52:35.095227+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We introduce two auxiliary variables and reformulate the Biot system as a four-field problem consisting of a dynamic Stokes-like system coupled with a diffusion equation. Crucially, this reformulation preserves the original interface conditions.
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We prove that the resulting method is unconditionally stable and derive optimal-order error estimates in the H1-norm. The analysis further shows that the scheme is robust with respect to extreme poroelastic parameters and avoids the locking effects inherent in standard formulations.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

[1]

Alnæs, J

[1]M. Alnæs, J. Blechta, J. Hake, A. Johansson, B. Kehlet, A. Logg, C. Richardson, J. Ring, M. E. Rognes, and G. N. Wells, The fenics project version 1.5, Archive of Numerical Software, 3 (2015). 34 [2]M. S. Alnæs, A. Logg, K. B. Ølgaard, M. E. Rognes, and G. N. Wells,Unified form language: A domain-specific language for weak for- mulations of partial dif...

work page 2015
[2]

Buka ˇc, S

[10]M. Buka ˇc, S. ˇCani´c, B. Muha, and Y. Wang,A computational algo- rithm for optimal design of a bioartificial organ scaffold architecture, Plos Computational Biology, 20 (2024), p. e1012079. [11]M. Buka ˇc, I. Yotov, R. Zakerzadeh, and P. Zunino,Partitioning strategies for the interaction of a fluid with a poroelastic material based on a nitsche’s co...

work page 2024
[3]

Cesmelioglu, H

[14]A. Cesmelioglu, H. Lee, A. Quaini, K. Wang, and S.-Y. Yi, Optimization-based decoupling algorithms for a fluid-poroelastic system, in Topics in numerical partial differential equations and scientific computing, Springer, 2016, pp. 79–98. 35 [15]W. Chen, M. Gunzburger, F. Hua, and X. Wang,A parallel robin- robin domain decomposition method for the stok...

work page 2016
[4]

Cuisiat, M

[17]F. Cuisiat, M. Gutierrez, R. Lewis, and I. Masters,Petroleum reservoir simulation coupling flow and deformation, in SPE Europec Fea- tured at EAGE Conference and Exhibition?, SPE, 1998, pp. SPE–50636. [18]A. Dalal, R. Durst, A. Quaini, and I. Yotov,A robin-robin split- ting method for the stokes-biot fluid-poroelastic structure interaction model, Math...

work page 1998
[5]

[22]P. H. Feenstra and C. A. Taylor,Drug transport in artery walls: a se- quential porohyperelastic-transport approach, Computer Methods in Biome- chanics and Biomedical Engineering, 12 (2009), pp. 263–276. [23]S. Guo, Y. Sun, Y. Wang, X. Yue, and H. Zheng,A fully parallelizable loosely coupled scheme for fluid-poroelastic structure interaction problems, ...

work page 2009
[6]

Oyekole and M

[29]O. Oyekole and M. Buka ˇc,Second-order, loosely coupled methods for fluid-poroelastic material interaction, Numerical Methods for Partial Dif- ferential Equations, 36 (2020), pp. 800–822. [30]C. Parrow and M. Buka ˇc,Stability and error analysis of a partitioned robin-robin method for fluid poroelastic structure interaction, Journal of Scientific Comp...

work page 2020
[7]

[31]P. J. Phillips and M. F. Wheeler,Overcoming the problem of locking in linear elasticity and poroelasticity: an heuristic approach, Computational Geosciences, 13 (2009), pp. 5–12. [32]T. Richter,Fluid-structure interactions: models, analysis and finite ele- ments, vol. 118, Springer,

work page 2009
[8]

Richter and T

[33]T. Richter and T. Wick,Finite elements for fluid-structure interac- tion in ale and fully eulerian coordinates, Computer Methods in Applied Mechanics and Engineering, 199 (2010), pp. 2633–2642. [34]A. Seboldt, O. Oyekole, J. Tamba ˇca, and M. Buka ˇc,Numerical modeling of the fluid-porohyperelastic structure interaction, SIAM journal on Scientific Com...

work page 2010

[1] [1]

Alnæs, J

[1]M. Alnæs, J. Blechta, J. Hake, A. Johansson, B. Kehlet, A. Logg, C. Richardson, J. Ring, M. E. Rognes, and G. N. Wells, The fenics project version 1.5, Archive of Numerical Software, 3 (2015). 34 [2]M. S. Alnæs, A. Logg, K. B. Ølgaard, M. E. Rognes, and G. N. Wells,Unified form language: A domain-specific language for weak for- mulations of partial dif...

work page 2015

[2] [2]

Buka ˇc, S

[10]M. Buka ˇc, S. ˇCani´c, B. Muha, and Y. Wang,A computational algo- rithm for optimal design of a bioartificial organ scaffold architecture, Plos Computational Biology, 20 (2024), p. e1012079. [11]M. Buka ˇc, I. Yotov, R. Zakerzadeh, and P. Zunino,Partitioning strategies for the interaction of a fluid with a poroelastic material based on a nitsche’s co...

work page 2024

[3] [3]

Cesmelioglu, H

[14]A. Cesmelioglu, H. Lee, A. Quaini, K. Wang, and S.-Y. Yi, Optimization-based decoupling algorithms for a fluid-poroelastic system, in Topics in numerical partial differential equations and scientific computing, Springer, 2016, pp. 79–98. 35 [15]W. Chen, M. Gunzburger, F. Hua, and X. Wang,A parallel robin- robin domain decomposition method for the stok...

work page 2016

[4] [4]

Cuisiat, M

[17]F. Cuisiat, M. Gutierrez, R. Lewis, and I. Masters,Petroleum reservoir simulation coupling flow and deformation, in SPE Europec Fea- tured at EAGE Conference and Exhibition?, SPE, 1998, pp. SPE–50636. [18]A. Dalal, R. Durst, A. Quaini, and I. Yotov,A robin-robin split- ting method for the stokes-biot fluid-poroelastic structure interaction model, Math...

work page 1998

[5] [5]

[22]P. H. Feenstra and C. A. Taylor,Drug transport in artery walls: a se- quential porohyperelastic-transport approach, Computer Methods in Biome- chanics and Biomedical Engineering, 12 (2009), pp. 263–276. [23]S. Guo, Y. Sun, Y. Wang, X. Yue, and H. Zheng,A fully parallelizable loosely coupled scheme for fluid-poroelastic structure interaction problems, ...

work page 2009

[6] [6]

Oyekole and M

[29]O. Oyekole and M. Buka ˇc,Second-order, loosely coupled methods for fluid-poroelastic material interaction, Numerical Methods for Partial Dif- ferential Equations, 36 (2020), pp. 800–822. [30]C. Parrow and M. Buka ˇc,Stability and error analysis of a partitioned robin-robin method for fluid poroelastic structure interaction, Journal of Scientific Comp...

work page 2020

[7] [7]

[31]P. J. Phillips and M. F. Wheeler,Overcoming the problem of locking in linear elasticity and poroelasticity: an heuristic approach, Computational Geosciences, 13 (2009), pp. 5–12. [32]T. Richter,Fluid-structure interactions: models, analysis and finite ele- ments, vol. 118, Springer,

work page 2009

[8] [8]

Richter and T

[33]T. Richter and T. Wick,Finite elements for fluid-structure interac- tion in ale and fully eulerian coordinates, Computer Methods in Applied Mechanics and Engineering, 199 (2010), pp. 2633–2642. [34]A. Seboldt, O. Oyekole, J. Tamba ˇca, and M. Buka ˇc,Numerical modeling of the fluid-porohyperelastic structure interaction, SIAM journal on Scientific Com...

work page 2010