Coupling schemes for modeling hydraulic fracture propagation using the XFEM

• Published in: Computer methods in applied mechanics and engineering 2013-01, Vol.253, p.305-322
• Main Authors: Gordeliy, Elizaveta, Peirce, Anthony
• Format: Article
• Language: English
• Subjects: Computational fluid dynamics; Free boundary problems; Hydraulic fractures; XFEM
• ISSN: 0045-7825; 1879-2138
• DOI: 10.1016/j.cma.2012.08.017

► We propose two novel XFEM-based schemes to couple the equations governing hydraulic fracture propagation. ► The first scheme uses the XFEM to resolve the Neumann to Dirichlet (ND) map with O(h2) accuracy and is suitable for fluid fronts that lag behind the fracture front. ► The second scheme uses a mixed XFEM formulation to enable singular tip asymptotic behavior to be combined with the O(h2) ND map. ► The coupled XFEM algorithms yield solutions that show excellent agreement with published reference solutions. We describe coupled algorithms that use the Extended Finite Element Method (XFEM) to solve the elastic crack component of the elasto-hydrodynamic equations that govern the propagation of hydraulic fractures in an elastic medium. With appropriate enrichment, the XFEM resolves the Neumann to Dirichlet (ND) map for crack problems with O(h2) accuracy and the Dirichlet to Neumann (DN) map with O(h) accuracy. For hydraulic fracture problems with a lag separating the fluid front from the fracture front, we demonstrate that the finite pressure field makes it possible to use a scheme based on the O(h2) XFEM solution to the ND map. To treat problems in which there is a coalescence of the fluid and fracture fronts, resulting in singular tip pressures, we developed a novel mixed algorithm that combines the tip width asymptotic solution with the O(h2) XFEM solution of the ND map away from the tips. Enrichment basis functions required for these singular pressure fields correspond to width power law indices λ>12, which are different from the index λ=12 of linear elastic fracture mechanics. The solutions obtained from the new coupled XFEM schemes agree extremely well with those of published reference solutions.