Core Concepts
This work extends the hydro-mechanical phase-field fracture model to non-isothermal conditions with micromechanics-based poroelasticity, where the Biot's coefficient and porosity are degraded based on the phase-field variable and the energy decomposition scheme. The model also employs an isotropic diffusion method to stabilize the advection-dominated heat flux and adapts the fixed stress split method to account for thermal stress.
Abstract
The key highlights and insights of this content are:
The authors extend the existing hydro-mechanical phase-field fracture model to thermo-poro-elastic media by introducing the thermal strain. The degradation of the thermo-poro-elastic strain energy is achieved through micromechanically derived Biot's coefficient, which depends on both the phase-field (damage) and the energy decomposition scheme.
A new approach is proposed to update the porosity that only depends on the strain change rather than the damage evolution, which is different from the existing models that rely on the phase-field profile.
An isotropic diffusion method is employed to stabilize the advective dominated heat flow in fractures, and the fixed stress split stabilization method is adapted to account for thermal stress in the numerical implementation.
The proposed model is verified against a series of analytical solutions, including Terzaghi's consolidation, thermal consolidation, and the plane strain hydraulic fracture propagation (KGD fracture). The results show that the new porosity update improves the accuracy in representing the sharp fracture behaviors compared to the existing phase-field dependent porosity models.
Numerical experiments demonstrate the effectiveness of the stabilization method and the intricate thermo-hydro-mechanical interactions during hydraulic fracturing with and without a pre-existing weak interface.
Stats
The authors provide the following key metrics and figures to support their analysis:
Dimensionless viscosity M = 3.8 × 10−7 for the toughness-dominated hydraulic fracture regime
Effective critical surface energy release rate Geff
c
adjusted based on the mesh size he and characteristic length ℓ
Comparisons of numerical results and analytical solutions for Terzaghi's consolidation, thermal consolidation, and KGD hydraulic fracture propagation