TY - JOUR
T1 - A modified extended finite element method for fluid-driven fractures incorporating variable primary energy loss mechanisms
AU - Klimenko, Denis
AU - Dahi Taleghani, Arash
N1 - Funding Information:
Denis Klimenko would like to thank The GDL Foundation for their generous support during conducting this research. Arash Dahi Taleghani would like to thank Ivo Babuŝka for his guidance and constructive discussions on developing a stable GFEM approach for fluid flow part. Chevron's Emerging Faculty Fund has supported authors while working on this research project.
Publisher Copyright:
© 2018 Elsevier Ltd
PY - 2018/6
Y1 - 2018/6
N2 - A coupled extended finite element method (XFEM) is presented here for modeling propagation of fluid-driven fractures in different regimes including toughness- and viscosity-dominated regimes. The extended finite element method allows to model growth, and coalescence of arbitrary discontinuities (fractures) without requiring the mesh to conform to discontinuities nor significant refinement near the fractures. Fluid-driven fractures propagation is a coupled, nonlinear and non-local problem with moving boundary conditions. The proposed method is based on the extended finite element method with modifications to incorporate variable stress singularity at the crack tips for the transition between toughness-dominated and viscosity-dominated regimes. These modifications consist of enriched functions that are initially inspired by the asymptotic analytical solutions. The standard extended finite element approximation is enriched by adding near tip asymptotic solutions just for displacements, however the proposed method introduces a consistent enriched function for fluid pressure calculations close to the fracture tips to catch the singularity. Additionally, a technique is presented to remove singularity issue for required numerical integrations. Green's functions concept is proposed here to expedite calculations. To circumvent violation of partition of unity and parasitic terms in the approximation space induced by the blending elements at the edge of the enriched domain, the ramp function is utilized to improve the convergence rate. Stress intensity factors are calculated using a new contour integral method that can handle cases with different tip singularities. The proposed technique is verified with the cases that have analytical solutions. Some examples are presented to show the advantages of this technique in comparison to the regular XFEM.
AB - A coupled extended finite element method (XFEM) is presented here for modeling propagation of fluid-driven fractures in different regimes including toughness- and viscosity-dominated regimes. The extended finite element method allows to model growth, and coalescence of arbitrary discontinuities (fractures) without requiring the mesh to conform to discontinuities nor significant refinement near the fractures. Fluid-driven fractures propagation is a coupled, nonlinear and non-local problem with moving boundary conditions. The proposed method is based on the extended finite element method with modifications to incorporate variable stress singularity at the crack tips for the transition between toughness-dominated and viscosity-dominated regimes. These modifications consist of enriched functions that are initially inspired by the asymptotic analytical solutions. The standard extended finite element approximation is enriched by adding near tip asymptotic solutions just for displacements, however the proposed method introduces a consistent enriched function for fluid pressure calculations close to the fracture tips to catch the singularity. Additionally, a technique is presented to remove singularity issue for required numerical integrations. Green's functions concept is proposed here to expedite calculations. To circumvent violation of partition of unity and parasitic terms in the approximation space induced by the blending elements at the edge of the enriched domain, the ramp function is utilized to improve the convergence rate. Stress intensity factors are calculated using a new contour integral method that can handle cases with different tip singularities. The proposed technique is verified with the cases that have analytical solutions. Some examples are presented to show the advantages of this technique in comparison to the regular XFEM.
UR - http://www.scopus.com/inward/record.url?scp=85046698018&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85046698018&partnerID=8YFLogxK
U2 - 10.1016/j.ijrmms.2018.04.012
DO - 10.1016/j.ijrmms.2018.04.012
M3 - Article
AN - SCOPUS:85046698018
SN - 1365-1609
VL - 106
SP - 329
EP - 341
JO - International Journal of Rock Mechanics and Mining Sciences
JF - International Journal of Rock Mechanics and Mining Sciences
ER -