Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

Deep beneath the Earth's surface, in the porous rocks that hold oil, gas, and geothermal heat, a silent battle is constantly being waged between solid stone and the fluids trapped within its tiny pores. When engineers try to extract these resources, they often pump fluids into the rock to force it apart, a process known as hydraulic fracturing. For decades, scientists have modeled this process by treating the rock as a simple, elastic sponge that stretches and snaps back. However, real rocks are more complex; they can also bend permanently, crumble, and slide along internal cracks, especially when squeezed by high pressures. Understanding exactly how these rocks break when both the solid skeleton and the fluid inside are pushing and pulling is crucial for safely drilling wells, managing waste storage, and even predicting natural geological events like icebergs cracking or magma pushing through the crust. The challenge lies in creating a mathematical description that captures the messy reality of rock breaking without getting lost in impossible calculations.

A team of researchers has revisited the fundamental equations used to describe this breaking process, finding that a common shortcut used in previous models leads to a significant error. In the world of rock mechanics, scientists use a "phase-field" approach, which treats a crack not as a sharp, invisible line, but as a blurry zone where the material is gradually weakening. This allows computers to simulate how a crack grows and changes shape. To do this accurately, the model must account for two distinct behaviors: when the rock is being pulled apart (tension) and when it is being squeezed together (compression). In tension, tiny cracks open up; in compression, they close and the rough surfaces inside them slide against each other, creating friction. The researchers discovered that many existing models, which are widely used in industry and academia, mix up the variables used to describe the energy of the system. Specifically, they found that treating the fluid pressure as a direct input in the energy equation, rather than the amount of fluid in the pores, creates a mathematical "jump" or discontinuity. This jump means the model predicts the rock's strength changes abruptly and unrealistically as it switches from being pulled to being squeezed, leading to inaccurate predictions of how and where a fracture will form.

To fix this, the researchers developed a new, consistent formulation that carefully tracks the relationship between the rock's deformation and the fluid content. They built their model on a microscopic view of the rock, imagining it as a collection of tiny, penny-shaped cracks filled with fluid. By analyzing how these micro-cracks open and close under different stresses, they derived a set of rules that ensure the rock's strength surface—the mathematical boundary that defines when it will break—remains smooth and continuous, no matter how the pressure changes. They proved that by using a specific type of flow rule that links the rock's permanent deformation directly to the forces acting on it, and by correctly including a missing coupling term in the energy equation, the model can accurately describe the transition from tension to compression. This correction is not just a minor tweak; it fundamentally changes how the model calculates the driving force that pushes a crack forward.

The team tested their new model against two very different scenarios to see if it held up. First, they simulated a classic hydraulic fracturing benchmark where fluid is injected into a rock to create a crack. They compared their results with a known analytical solution, which is a precise mathematical answer derived for ideal conditions. The new model matched this solution much more closely than the older, mixed formulations, correctly predicting the pressure needed to start the fracture and the width of the crack as it grew. The older models had tended to overestimate the fluid pressure required, a flaw that could lead to inefficient or unsafe drilling operations. In the second test, they simulated a rock sample being squeezed from two sides, a condition that mimics the deep underground environment where rocks are under immense pressure. Here, the rock did not just snap; it deformed plastically, forming a shearing band where the material slid past itself before breaking. The new model successfully reproduced this complex behavior, capturing both the shear fractures caused by mechanical squeezing and the tension fractures driven by fluid injection.

The results of these simulations reveal that the plastic deformation of the rock acts as a significant barrier to fracture growth. When the rock is allowed to deform permanently, it absorbs more energy, which means the fluid pressure inside the crack must be higher to keep the fracture open and growing. This extra resistance changes the entire profile of the fracture, making it narrower and requiring more force to propagate than a purely elastic model would suggest. Furthermore, in the compression tests, the model showed how the pore pressure within the rock changes as the fracture forms. As the rock deforms and creates space for the fluid to move, the pressure drops locally, drawing more fluid into the damaged zone. This interaction creates a feedback loop where the fluid flow and the rock's mechanical failure are tightly coupled, a dynamic that the new model captures with high fidelity.

By correcting the mathematical foundation of these models, the researchers have provided a more reliable tool for engineers and scientists working with porous media. The new formulation ensures that the transition between different types of rock failure is smooth and physically realistic, eliminating the artificial jumps that plagued previous attempts. This accuracy is vital for applications ranging from optimizing oil and gas production to designing safe underground storage for carbon capture. The study demonstrates that even in complex systems where solid and fluid interact, a careful, consistent approach to the underlying physics can reveal the true nature of how materials break. The model now stands ready to help predict the behavior of the Earth's crust with greater confidence, bridging the gap between theoretical mechanics and the messy, plastic reality of the ground beneath our feet.