Search

In this work, we present a domain-based algorithm to simulate the propagation of a plane-strain hydraulic fracture in a zero-toughness permeable elastic medium. The algorithm utilizes a domain-based method to solve the elasticity equation and integrates a multi-scale tip asymptote, which is particular to hydraulic fractures, into this framework. This integration is key to accurately model the energy dissipation and the fluid leak-off in the fracture tip region. The algorithm combines a 2D finite volume method (FVM) for solving the elasticity equation with a 1D FVM for solving the nonlinear lubrication equation. Incorporating the far-field asymptotics and using a moving-mesh scheme reduces the computational burden while improving the accuracy of the scheme. The paper concludes with an analysis of the numerical results. This study demonstrates the potential of this domain-based approach for modeling hydraulic fractures in poroelastic media. [ABSTRACT FROM AUTHOR]

The contraction or expansion of a cylindrical cavity in an elastoplastic medium is usually analyzed from a continuum based approach with a plasticity constitutive model. However, localized deformations, which are rooted in the post‐failure softening response of geomaterials, are observed in the form of spiral‐shaped fractures in laboratory tests. An alternative approach based on dislocation theory is introduced in this paper for modeling cavity contraction or expansion. In this model, several equally spaced spiral‐shaped shear fractures initiate and propagate away from the cavity within the linearly elastic medium. The Mohr‐Coulomb criterion and a dilatancy rule are imposed on the shear fractures to constrain the stresses and the displacement jumps. The direction of fracture propagation is determined by minimizing plastic dissipation. The displacement discontinuity method is used to discretize the shear and normal displacement jumps along the fracture and solve the problem numerically. The calculated crack path follows a logarithmic‐like spiral, similar to the slip lines predicted by plasticity theory. The relationship between the pressure and radial displacement at the cavity boundary converge towards the classical elastoplastic solution as the number of fracture branches increases. [ABSTRACT FROM AUTHOR]

Fluid injection‐induced deformation around a cylindrical cavity is of particular interest in the area of subsurface energy extraction. In this study, a model is proposed to analyze the time‐dependent expansion of a cavity caused by fluid injection in an elastoplastic dry porous medium. This problem is characterized by the existence of two moving boundaries, a permeation front and an elastoplastic interface, which leads to distinct time‐dependent zones governed by different sets of equations. The interplay between these two boundaries leads to three phases of solution. The a priori unknown partitioning of the injection rate into the rate of change of cavity volume and infiltration in the porous medium necessitates the introduction of a time‐dependent permeation coefficient as one of the primary variables. The method of solution takes advantage of the quasi‐static and quasi‐stationary nature of the problem, which makes it possible to treat time as a parameter rather than a variable. It follows that the problem can be solved in two steps. In a first step, closed‐form expressions for the pore pressure, displacement, and stress fields in each zone are derived, with parameters in these expressions depending explicitly on four time‐dependent variables, namely the positions of the two interfaces, the cavity radius, and the permeation coefficient. In a second step, the rate equations governing the evolution of these variables during different phases are derived and solved numerically. The paper concludes with a parametric analysis of the influence of the stiffness and strength of the material on the solution. [ABSTRACT FROM AUTHOR]

This paper describes a numerical algorithm for solving the classic problem of a plane strain (KGD) fracture propagating in an impermeable elastic medium with zero toughness. The method, which takes advantage of the self-similar nature of the solution, combines a domain-based scheme to solve the elasticity equations and a finite volume method to solve the nonlinear lubrication equation. This work represents a first step towards developing a model able to account for pore pressure diffusion in the medium and corresponding poroelastic effects, noting that these processes are more efficiently solved using a domain-based rather than a boundary integral method. To enhance the efficiency and accuracy of the numerical scheme, the far-field crack asymptotics is embedded in the discretized elastic relationship between the fluid pressure and the crack opening, while the coupled fluid-solid tip asymptote is enforced in a weak form when solving the nonlinear lubrication equation. The proposed technique yields results that closely match the analytical solution, even with a coarse mesh. This approach offers potential for addressing more complex hydraulic fracturing problems in the future. [ABSTRACT FROM AUTHOR]