Phase-field modeling of rock fractures with roughness

TL;DR

This paper introduces a novel phase-field modeling framework that incorporates surface roughness effects into rock fracture simulations, bridging discrete and continuous approaches without explicit geometry tracking.

Contribution

It transforms a discrete rough fracture model into a continuous phase-field formulation, enabling accurate simulation of rough rock fractures within standard finite element methods.

Findings

Excellent agreement with discrete methods for shear and compression tests

Effective simulation of complex crack growth from rough discontinuities

Framework bridges discrete and continuous modeling approaches

Abstract

Phase-field modeling -- a continuous approach to discontinuities -- is gaining popularity for simulating rock fractures due to its ability to handle complex, discontinuous geometry without an explicit surface tracking algorithm. None of the existing phase-field models, however, incorporates the impact of surface roughness on the mechanical response of fractures -- such as elastic deformability and shear-induced dilation -- despite the importance of this behavior for subsurface systems. To fill this gap, here we introduce the first framework for phase-field modeling of rough rock fractures. The framework transforms a displacement-jump-based discrete constitutive model for discontinuities into a strain-based continuous model, without any additional parameter, and then casts it into a phase-field formulation for frictional interfaces. We illustrate the framework by constructing a particular phase-field form employing a rock joint model originally formulated for discrete modeling. The results obtained by the new formulation show excellent agreement with those of a well-established discrete method for a variety of problems ranging from shearing of a single discontinuity to compression of fractured rocks. It is further demonstrated that the phase-field formulation can well simulate complex crack growth from rough discontinuities. Consequently, our phase-field framework provides an unprecedented bridge between a discrete constitutive model for rough discontinuities -- common in rock mechanics -- and the continuous finite element method -- standard in computational mechanics -- without any algorithm to explicitly represent discontinuity geometry.