A hybrid numerical approach for hydraulic fracturing in a naturally fractured formation combining the XFEM and phase-field model. (August 2022)