An immersed phase field fracture model for microporomechanics with Darcy–Stokes flow

This paper presents an immersed phase field model designed to predict the fracture-induced flow due to brittle fracture in vuggy porous media. Due to the multiscale nature of pores in the vuggy porous material, crack growth may connect previously isolated pores, which leads to flow conduits. This mechanism has important implications for many applications such as disposal of carbon dioxide and radioactive materials and hydraulic fracture and mining. To understand the detailed microporomechanics that causes the fracture-induced flow, we introduce a new phase field fracture framework where the phase field is not only used as an indicator function for damage of the solid skeleton but also used as an indicator of the pore space. By coupling the Stokes equation that governs the fluid transport in the voids, cavities, and cracks and Darcy’s flow in the deformable porous media, our proposed model enables us to capture the fluid–solid interaction of the pore fluid and solid constituents during crack growth. Numerical experiments are conducted to analyze how the presence of cavities affects the accuracy of predictions based on the homogenized effective medium during crack growth.