Abstract
Fluid-driven fracture processes are central to the development of subsurface energy systems such as geothermal and hydrocarbon reservoirs. Although phase-field formulations have become a widely used tool for describing fracture initiation and growth, the diffuse representation of cracks makes it difficult to resolve flow behavior accurately inside discrete fracture networks (DFNs) and to represent hydro-mechanical coupling in a sharp-interface sense. This study develops a hybrid-dimensional iterative framework for lubrication-flow simulation in deformable fractured geomaterials. By leveraging phase-field point clouds together with non-conforming discretization schemes for both the solid matrix and fracture domains, the proposed framework enables the dynamic reconstruction of evolving fracture networks. The theoretical formulation and numerical implementation of the coupling strategy are presented in detail. Hydraulic benchmark examples verify the performance of the fluid flow solver under various physical conditions. The classical Sneddon problem and Khristianovic-Geertsma-de Klerk (KGD) model are employed to validate the solid deformation solver, confirming accurate predictions of crack opening displacement and mesh independence in fracture width calculation. Additional simulations with complex pre-existing fracture patterns further demonstrate the applicability of the framework to coupled hydro-mechanical analysis in fractured media.