TY - JOUR

T1 - A phase-field fracture model in thermo-poro-elastic media with micromechanical strain energy degradation

AU - Liu, Yuhao

AU - Yoshioka, Keita

AU - You, Tao

AU - Li, Hanzhang

AU - Zhang, Fengshou

N1 - Publisher Copyright: © 2024 Elsevier B.V.

PY - 2024/6/24

Y1 - 2024/6/24

N2 - This work extends the hydro-mechanical phase-field fracture model to non-isothermal conditions with micromechanics based poroelasticity, which degrades Biot's coefficient not only with the phase-field variable (damage) but also with the energy decomposition scheme. Furthermore, we propose a new approach to update porosity solely determined by the strain change rather than damage evolution as in the existing models. As such, these poroelastic behaviors of Biot's coefficient and the porosity dictate Biot's modulus and the thermal expansion coefficient. For numerical implementation, we employ an isotropic diffusion method to stabilize the advection-dominated heat flux and adapt the fixed stress split method to account for the thermal stress. We verify our model against a series of analytical solutions such as Terzaghi's consolidation, thermal consolidation, and the plane strain hydraulic fracture propagation, known as the KGD fracture. Finally, numerical experiments demonstrate the effectiveness of the stabilization method and intricate thermo-hydro-mechanical interactions during hydraulic fracturing with and without a pre-existing weak interface.

AB - This work extends the hydro-mechanical phase-field fracture model to non-isothermal conditions with micromechanics based poroelasticity, which degrades Biot's coefficient not only with the phase-field variable (damage) but also with the energy decomposition scheme. Furthermore, we propose a new approach to update porosity solely determined by the strain change rather than damage evolution as in the existing models. As such, these poroelastic behaviors of Biot's coefficient and the porosity dictate Biot's modulus and the thermal expansion coefficient. For numerical implementation, we employ an isotropic diffusion method to stabilize the advection-dominated heat flux and adapt the fixed stress split method to account for the thermal stress. We verify our model against a series of analytical solutions such as Terzaghi's consolidation, thermal consolidation, and the plane strain hydraulic fracture propagation, known as the KGD fracture. Finally, numerical experiments demonstrate the effectiveness of the stabilization method and intricate thermo-hydro-mechanical interactions during hydraulic fracturing with and without a pre-existing weak interface.

KW - Fixed stress split

KW - Hydraulic fracturing

KW - Isotropic diffusion method

KW - Phase-field

KW - Thermo-hydro-mechanical coupling

KW - Thermo-poroelasticity

UR - http://www.scopus.com/inward/record.url?scp=85196637142&partnerID=8YFLogxK

U2 - 10.1016/j.cma.2024.117165

DO - 10.1016/j.cma.2024.117165

M3 - Article

VL - 429.2024

JO - Computer methods in applied mechanics and engineering

JF - Computer methods in applied mechanics and engineering

SN - 0045-7825

IS - 1 September

M1 - 117165

ER -