arxiv: 2604.18483 · v1 · submitted 2026-04-20 · ⚛️ physics.geo-ph · physics.flu-dyn

Recognition: unknown

Steadily moving semi-infinite fracture in plane poroelasticity

Evgenii Kanin , Andreas M\"ori , Dmitry Garagash , Brice Lecampion

Authors on Pith no claims yet

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

classification ⚛️ physics.geo-ph physics.flu-dyn

keywords poroelasticityboundary integral equationssemi-infinite fractureplane strainfluid diffusionfracture propagationcoupled hydro-mechanical modelgeomechanics

0 comments

The pith

The paper develops a boundary integral formulation for steadily moving semi-infinite fractures in poroelastic media.

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

The authors construct a fully coupled boundary integral method to simulate semi-infinite fractures that propagate steadily through poroelastic rock. The technique superposes fundamental solutions for fluid sources and dislocations in time and space to obtain equations for surface tractions and pore pressure, which are then solved numerically for opening and fluid exchange when loads are prescribed. This capability is useful because it handles the interaction between rock deformation and fluid flow in a steady-state setting without needing a full transient simulation. Verification against analytical solutions for tensile and shear fractures confirms the accuracy of the implementation. The resulting framework can analyze how fluid pressure influences fracture behavior in permeable media.

Core claim

What carries the argument

The fully coupled boundary integral formulation derived from temporal and spatial superposition of poroelastic fundamental solutions for an instantaneous fluid source and for edge dislocations (normal and slip modes) in a steadily moving semi-infinite plane strain geometry.

If this is right

The numerical method computes fracture opening, slip, and cumulative fluid exchange rate from prescribed surface tractions and pore pressures.
The model captures the full coupling between mechanical stresses and pore fluid diffusion during steady propagation.
The framework provides a robust tool for analyzing coupled fracture-fluid interactions in permeable poroelastic media.
The method can be adapted to broader classes of elasto-diffusive problems by modifying the underlying physical parameters.

Where Pith is reading between the lines

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

If the steady-state assumption holds, this method could approximate segments of variable-speed fracture growth in larger models.
Applying the formulation to field data from hydraulic fracturing operations would test its practical utility for predicting fluid leak-off.
Similar superposition techniques might extend to other moving boundary problems in diffusive materials.

Load-bearing premise

The tractions and pore fluid pressures on the fracture surfaces must be known ahead of time, and the propagation must be steady and semi-infinite for the superposition of fundamental solutions to hold directly.

What would settle it

Comparing the numerically obtained fracture opening profile and fluid exchange rate to the exact analytical solution for the case of a stress-free tensile fracture with imposed exponential pore fluid pressure would directly test the correctness of the derived boundary integral equations.

Figures

Figures reproduced from arXiv: 2604.18483 by Andreas M\"ori, Brice Lecampion, Dmitry Garagash, Evgenii Kanin.

**Figure 1.** Figure 1: Fundamental loadings: (a) a fluid source and (b) an edge dislocation located at ⃗𝜒 = (𝜒1 , 𝜒2 ). The observation point is denoted by ⃗𝑥 = (𝑥1 , 𝑥2 ). 2.2.1. Instantaneous fluid source solutions A fluid source (Fig. 1a) corresponds to the injection of fluid mass at a point within the medium. The instantaneous fluid source is modeled by the term 𝛾 = 𝛿(⃗𝑥 − ⃗𝜒)𝛿(𝑡 − 𝜏), where 𝛿(⋅) denotes the Dirac delta func… view at source ↗

**Figure 2.** Figure 2: Fundamental loadings moving at constant velocity 𝑉 : (a) a fluid source and (b) an edge dislocation. Both the fixed (𝑋, 𝑌 ) and moving (𝑥, 𝑦) coordinate systems are indicated. so that the solution is steady in the moving frame (𝑥, 𝑦). Consequently, time derivatives in the fixed frame convert into spatial derivatives in the moving coordinates, and the transient problem reduces to a steady one in the co-movi… view at source ↗

**Figure 3.** Figure 3: Characteristic shapes of the pore pressure and stresses induced by a steadily moving fluid source (a) and normal edge dislocation (b), plotted versus the dimensionless coordinate 𝜉 (Eq. (24)). The normalized undrained-drained Poisson’s ratio contrast is taken as 𝛽 = 0.125. • fluid source:  𝑠𝑚(𝜉 → −∞) ∼ 𝑒 2𝜉 √ 𝜋 2|𝜉| ,  𝑠𝑚(𝜉 → +∞) ∼ √ 𝜋 2𝜉 − 1 8 √ 𝜋 2 1 𝜉 3∕2 ;  𝑠𝑚 𝑥𝑥 (𝜉 → −∞) ∼ 1 𝜉 + 𝑒 2𝜉 2 √ 𝜋 2 1 |𝜉| … view at source ↗

**Figure 4.** Figure 4: Characteristic shapes of the stresses induced by a steadily moving normal edge dislocation, with the Cauchy singularity removed by multiplying the governing equations by the dimensionless coordinate 𝜉 (Eq. (24)). The resulting distributions are plotted as functions of 𝜉. The value of the normalized undrained-drained Poisson’s ratio contrast is 𝛽 = 0.125. Section 6, because Gauss-Chebyshev quadrature scheme… view at source ↗

**Figure 5.** Figure 5: Schematic of a semi-infinite (a) tensile and (b) shear fractures propagating steadily at a constant velocity 𝑉 along the negative 𝑋−axis. The fixed (𝑋, 𝑌 ) and moving (𝑥, 𝑦) coordinate systems are shown. The distributions of the pore fluid pressure 𝑝𝑓 (𝑥), the normal stress 𝜎𝑦𝑦 (a), and the shear stress 𝜎𝑥𝑦(𝑥) (b) on the fracture surfaces are depicted. The fracture opening 𝑤(𝑥) (a) and slip 𝑑(𝑥) (b) are in… view at source ↗

**Figure 6.** Figure 6: Dimensionless stress intensity factor for a semi-infinite tensile fracture with impermeable surfaces, loaded by an exponential normal stress, normalized by the dimensionless elastic stress intensity factor (Eq. (55)), as a function of (a) the normalized undrained-drained Poisson’s ratio contrast 𝛽 and (b) the dimensionless loading distance . The numerical solution is shown by black lines, while the analyt… view at source ↗

**Figure 7.** Figure 7: Relative difference, 𝜀 (in %), between numerical and analytical (Eq. (57)) solutions for the dimensionless stress intensity factor  of a semi-infinite tensile fracture with impermeable surfaces subjected to exponential normal loading. (a) Relative difference as a function of the parameters 𝛽 and . (b) Convergence analysis showing the dependence of the maximum (blue circles) and mean (green circles) rela… view at source ↗

**Figure 8.** Figure 8: Dimensionless (a) opening and (b) pore fluid pressure profiles along a semi-infinite tensile fracture with impermeable surfaces, loaded by an exponential normal stress, shown for  = {10−2 , 10−1 , 1, 10, 102 , 103} with 𝛽 = 0.125 and 𝜂 = 0.25. The numerical results are shown by black lines, while the analytical solutions (Eq. (58) with  given by Eq. (57) for the opening and Eq. (59) for the pressure) are… view at source ↗

**Figure 9.** Figure 9: Dimensionless pore fluid pressure around a semi-infinite tensile fracture with impermeable surfaces, loaded by an exponential normal stress. (a) Profiles along lines parallel to the fracture at 𝜁 = {0, 1, 2, 3, 4, 5} for  = 5, 𝛽 = 0.125, and 𝜂 = 0.25. The numerical results obtained from Eq. (61) are plotted as solid black lines, whereas the reference semianalytical solution from Eq. (62) is shown as dash… view at source ↗

**Figure 10.** Figure 10: Dimensionless stress intensity factor for a semi-infinite tensile fracture with permeable surfaces, subjected to prescribed zero pore fluid pressure and an exponential normal stress, normalized by the dimensionless elastic stress intensity factor (Eq. (55)), as a function of (a) the normalized undrained-drained Poisson’s ratio contrast 𝛽 and (b) the dimensionless loading distance . The poroelastic stress… view at source ↗

**Figure 11.** Figure 11: Dimensionless (a) opening and (b) fluid displacement function profiles along a semi-infinite tensile fracture with permeable surfaces, subjected to prescribed zero pore fluid pressure and an exponential normal stress, shown for  = {10−2 , 10−1 , 1, 10, 102 , 103 , 104} with 𝛽 = 0.125 and 𝜂 = 0.25. The fluid displacement function is depicted with reversed sign. The numerical results are shown by black lin… view at source ↗

**Figure 12.** Figure 12: Dimensionless pore fluid pressure around a semi-infinite tensile fracture with permeable surfaces, subjected to prescribed zero pore fluid pressure and an exponential normal stress. (a) Profiles along lines parallel to the fracture at 𝜁 = {0, 1, 2, 3, 4, 5} for  = 5, 𝛽 = 0.125, and 𝜂 = 0.25. The numerical results obtained from Eq. (67) are plotted as solid black lines, whereas the reference semi-analytic… view at source ↗

**Figure 13.** Figure 13: Dimensionless tensile stress profiles ahead (𝜉 < 0) of a semi-infinite tensile fracture with permeable surfaces, subjected to prescribed zero pore fluid pressure and an exponential normal stress. Results are shown for  = {10−2 , 10−1 , 1, 10, 102 , 103 , 104} with 𝛽 = 0.125 and 𝜂 = 0.25. Both the tensile stress and the 𝜉−coordinate are plotted with reversed sign. The numerical results are plotted by blac… view at source ↗

**Figure 14.** Figure 14: Dimensionless stress intensity factor for a semi-infinite tensile fracture with permeable surfaces, stress-free and subjected to an imposed exponential pore pressure, as a function of (a) the normalized undrained-drained Poisson’s ratio contrast 𝛽, (b) the dimensionless loading distance , and (c) the poroelastic stress coefficient 𝜂. The dimensionless stress intensity factor is shown with reversed sign. … view at source ↗

**Figure 15.** Figure 15: Dimensionless (a) opening and (b) fluid displacement function profiles along a semi-infinite tensile fracture with permeable surfaces, stress-free and subjected to an imposed exponential pore pressure, shown for  = {10−2 , 10−1 , 1, 10, 102 , 103 , 104} with 𝛽 = 0.125 and 𝜂 = 0.25. The opening profile is shown with reversed sign. The numerical results are shown by black lines, while the analytical solut… view at source ↗

**Figure 16.** Figure 16: Dimensionless stress intensity factor for a semi-infinite shear fracture with permeable surfaces, subjected to prescribed zero pore fluid pressure and a uniform shear stress applied over a finite region, normalized by the dimensionless elastic stress intensity factor (Eq. (76)), as a function of (a) the normalized undrained-drained Poisson’s ratio contrast 𝛽 and (b) the dimensionless loading distance . T… view at source ↗

**Figure 17.** Figure 17: Dimensionless slip profiles along a semi-infinite shear fracture with permeable surfaces, subjected to prescribed zero pore fluid pressure and a uniform shear stress applied over a finite region, shown for  = {10−2 , 10−1 , 1, 10, 102 , 103 , 104} with 𝛽 = 0.125. The numerical results are shown by black lines, while the analytical solution (Eq. (79)), are indicated by dashed red lines. with  given by E… view at source ↗

read the original abstract

We present a fully coupled boundary integral formulation for modeling steadily propagating semi-infinite plane strain fractures in poroelastic media. By combining fundamental solutions of plain strain poroelasticity for instantaneous fluid source and edge dislocations (normal and slip modes) with temporal and spatial superposition principles, we derive boundary integral equations governing the tractions (normal and shear stresses) and pore fluid pressure on the fracture surfaces. Assuming prescribed tractions and pore fluid pressure profiles, we develop a numerical methodology to solve the governing equations for fracture opening, slip, and cumulative fluid exchange rate. The formulation is systematically verified on several relevant problems, including the case of a tensile fracture with exponential normal loading, a stress-free tensile fracture with an imposed exponential pore fluid pressure, and a shear fracture under uniform shear loading over a finite region, demonstrating excellent agreement with analytical solutions. The framework provides a robust tool for analyzing coupled fracture-fluid interactions in permeable poroelastic media and can be adapted to broader classes of elasto-diffusive problems by modifying the underlying physical parameters.

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

This paper gives a clean boundary-integral method for steady semi-infinite poroelastic fractures, verified against three analytical cases.

read the letter

The main contribution is a set of boundary integral equations for a steadily propagating semi-infinite fracture in plane poroelasticity. The authors combine the usual fundamental solutions for instantaneous fluid sources and normal/slip edge dislocations, then apply spatial and temporal superposition in the moving frame to relate the prescribed surface tractions and pore pressure to the resulting opening, slip, and fluid exchange rate. They also outline a numerical scheme to solve those equations. That specific combination for the steady semi-infinite case appears new. The verification is the part that stands out. They compare the numerical output to closed-form solutions on three problems: exponential normal loading, imposed exponential pressure on a stress-free fracture, and uniform shear over a finite patch. The agreement is reported as excellent across the board, which gives reasonable that the kernels and discretization are handled correctly. The derivation rests on standard linear superposition and established poroelastic kernels, so there is little circularity. The assumptions are stated plainly: steady motion, semi-infinite geometry, and known surface tractions plus pressure. Those limits mean the method cannot be dropped straight into problems with self-determined loading, finite length, or strong transients without further work. The paper is aimed at researchers who already use boundary integrals for fracture problems in geomechanics or energy applications. It supplies a practical tool for the steady semi-infinite regime that was not previously available in this coupled form. I would send it for peer review. The core formulation is internally consistent and the checks are direct enough to justify referee time.

Referee Report

2 major / 2 minor

Summary. The paper presents a fully coupled boundary integral equation (BIE) formulation for steadily propagating semi-infinite plane-strain fractures in poroelastic media. It combines fundamental solutions for instantaneous fluid sources and normal/slip edge dislocations with spatial-temporal superposition in the moving frame to obtain integral equations relating surface tractions, pore pressure, opening, slip, and fluid exchange. Assuming prescribed traction and pressure profiles on the fracture faces, a numerical scheme is developed and verified against closed-form analytical solutions on three test problems (exponential normal loading on a tensile fracture, imposed exponential pore pressure on a stress-free tensile fracture, and uniform shear over a finite patch on a shear fracture), with reported pointwise agreement.

Significance. If the central formulation and verifications hold, the work supplies a parameter-free, reproducible numerical tool for coupled hydro-mechanical analysis of propagating fractures in permeable poroelastic solids. This is directly relevant to geophysical applications such as hydraulic fracturing and fault slip. The explicit use of established poroelastic fundamental solutions plus linear superposition, together with systematic comparison to analytical benchmarks on three distinct loading configurations, constitutes a clear strength that supports extensibility to other elasto-diffusive problems.

major comments (2)

[Numerical methodology and verification] The numerical methodology section does not specify the discretization (collocation vs. Galerkin, element type, quadrature rule for singular kernels) or the convergence criteria employed when solving the resulting system for opening and slip; without these details the reported excellent agreement with the three analytical solutions cannot be independently reproduced or assessed for robustness under mesh refinement.
[Formulation and governing equations] The derivation assumes that the fracture propagates at constant speed with all fields steady in the moving frame; the manuscript should explicitly state the range of validity (e.g., when transient diffusion length scales remain comparable to the process zone) because violation of this assumption would invalidate the superposition used to obtain the BIE kernels.

minor comments (2)

A short table summarizing the three verification cases (loading type, analytical solution reference, reported error metric) would improve readability and allow quick comparison of accuracy across normal, pressure-driven, and shear modes.
Notation for the poroelastic parameters (e.g., Biot coefficient, diffusivity) should be collected in a single nomenclature table or defined at first appearance to avoid ambiguity when the formulation is adapted to other materials.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive comments and positive recommendation. We address each major point below and will revise the manuscript to improve reproducibility and clarity.

read point-by-point responses

Referee: [Numerical methodology and verification] The numerical methodology section does not specify the discretization (collocation vs. Galerkin, element type, quadrature rule for singular kernels) or the convergence criteria employed when solving the resulting system for opening and slip; without these details the reported excellent agreement with the three analytical solutions cannot be independently reproduced or assessed for robustness under mesh refinement.

Authors: We agree that these implementation details are required for independent reproduction and assessment of robustness. The manuscript currently lacks an explicit description of the discretization approach, quadrature rules, and convergence criteria. In the revised version we will add a dedicated paragraph in the numerical methodology section specifying the collocation scheme, element type, treatment of singular kernels, and the residual-based convergence tolerance used for the reported solutions. revision: yes
Referee: [Formulation and governing equations] The derivation assumes that the fracture propagates at constant speed with all fields steady in the moving frame; the manuscript should explicitly state the range of validity (e.g., when transient diffusion length scales remain comparable to the process zone) because violation of this assumption would invalidate the superposition used to obtain the BIE kernels.

Authors: The steady-state assumption in the moving frame is fundamental to the superposition that yields the time-independent kernels. We will revise the formulation section to include an explicit statement of the validity range, noting that the approach holds when the fracture velocity is sufficiently high that diffusive transients remain negligible relative to the process-zone length scale (i.e., when the diffusion length is comparable to or smaller than the characteristic length over which the prescribed tractions or pressures vary). revision: yes

Circularity Check

0 steps flagged

Derivation relies on standard fundamental solutions and superposition; no circularity

full rationale

The paper constructs its boundary integral formulation by superposing established fundamental solutions (instantaneous fluid source and normal/slip edge dislocations) from plane poroelasticity together with temporal-spatial superposition in the moving frame. These kernels are external, pre-existing results rather than derived or fitted inside the paper. The governing equations are solved numerically under explicitly prescribed surface tractions and pore pressure, then verified pointwise against independent closed-form analytical solutions on three separate test problems. No step equates a derived quantity to its own input by definition, renames a fitted parameter as a prediction, or rests on a self-citation chain whose validity is presupposed by the present work. The derivation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard linear superposition and pre-existing fundamental solutions of poroelasticity; no new free parameters, ad-hoc axioms, or invented physical entities are introduced.

axioms (2)

standard math Linear superposition in space and time holds for the governing poroelastic equations
Invoked to combine instantaneous fluid-source and edge-dislocation solutions into the moving-fracture problem.
domain assumption Fundamental solutions exist for an instantaneous fluid source and for normal/slip edge dislocations in plane-strain poroelasticity
These known solutions are taken as given from prior poroelasticity literature and superposed.

pith-pipeline@v0.9.0 · 5490 in / 1503 out tokens · 57869 ms · 2026-05-10T02:42:47.611602+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

6 extracted references

[1]

Pergamon Press

Fundamentals of poroelasticity, in: Comprehensive Rock Engineering. Pergamon Press. Detournay,E.,Garagash,D.,2003. Thenear-tipregionofafluid-drivenfracturepropagatinginapermeableelasticsolid. JournalofFluidMechanics 494, 1–32. Dontsov, E.,

2003
[2]

Computer Methods in Applied Mechanics and Engineering 313, 53–84

A multiscale implicit level set algorithm (ilsa) to model hydraulic fracture propagation incorporating combined viscous, toughness, and leak-off asymptotics. Computer Methods in Applied Mechanics and Engineering 313, 53–84. Dontsov,E.V.,2016. Anapproximatesolutionforapenny-shapedhydraulicfracturethataccountsforfracturetoughness,fluidviscosityandleak-off. ...

2016
[3]

International Journal of Rock Mechanics and Mining Sciences 108, 198–208

Influence of pore pressure on the development of a hydraulic fracture in poroelastic medium. International Journal of Rock Mechanics and Mining Sciences 108, 198–208. Hills,D.A.,Kelly,P.,Dai,D.,Korsunsky,A.,2013. Solutionofcrackproblems: thedistributeddislocationtechnique.volume44. SpringerScience & Business Media. Irwin, G.R.,

2013
[4]

Computer Methods in Preprint submitted to Journal of the Mechanics and Physics of SolidsPage 35 of 36 Applied Mechanics and Engineering 283, 881–908

Modeling multi-scale processes in hydraulic fracture propagation using the implicit level set algorithm. Computer Methods in Preprint submitted to Journal of the Mechanics and Physics of SolidsPage 35 of 36 Applied Mechanics and Engineering 283, 881–908. Peirce,A.,Detournay,E.,2008. Animplicitlevelsetmethodformodelinghydraulicallydrivenfractures. Computer...

2008
[5]

International Journal for Numerical and Analytical Methods in Geomechanics 13, 215–224

A two-dimensional poroelastic displacement discontinuity method for hydraulic fracture simu- lation. International Journal for Numerical and Analytical Methods in Geomechanics 13, 215–224. Viesca,R.C.,Garagash,D.I.,2018. Numericalmethodsforcoupledfractureproblems. JournaloftheMechanicsandPhysicsofSolids113,13–34. Weertman, J.,

2018
[6]

Journal of the Mechanics and Physics of Solids 118, 15–39

Steady-state crack growth in polymer gels: a linear poroelastic analysis. Journal of the Mechanics and Physics of Solids 118, 15–39. Yu,Y.,Landis,C.M.,Huang,R.,2020b. Poroelasticeffectsonsteadystatecrackgrowthinpolymergelsunderplanestress. MechanicsofMaterials 143, 103320. Zhou,S.,Zhuang,X.,Rabczuk,T.,2018. Aphase-fieldmodelingapproachoffracturepropagatio...

2018