Total: 1
Planar 3D models of hydraulic fracturing provide a practical balance between models with restrictive geometric assumptions and fully 3D simulators, capturing fractures with arbitrary planar footprints at moderate computational cost. Nevertheless, applications such as treatment design optimization and mini-frac test interpretation require large ensembles of simulations, for which the cost of planar 3D models remains a significant bottleneck. This work presents acceleration strategies for the planar 3D Implicit Level Set Algorithm (ILSA) to reduce simulation runtime while preserving numerical accuracy. A unified planar 3D ILSA scheme that consolidates the nested loops of the elastohydrodynamic solver and the front tracking algorithm into a single iterative process is introduced. A matrix splitting approach is applied to the linearized elastohydrodynamic system, moving the dense part of the elasticity operator to the right-hand side and yielding a sparse system matrix that can be solved more efficiently. Anderson acceleration is incorporated into the solution of the elastohydrodynamic system to improve convergence under varying fracture geometry. Additionally, a predictor--corrector scheme is examined with the proposed methods to assess their combined effect. Each technique is evaluated individually and in combination on both the reference and unified planar 3D ILSA schemes across five benchmark cases. Numerical experiments demonstrate that the unified scheme alone delivers an average 2.5x speedup, reaching 5.7x for the sandglass geometry. The combined application of all techniques achieves an average 4x speedup and up to 11x for the sandglass case, with the relative discrepancy in fracture aperture below 5% compared with the reference scheme.