Quantitative phase field modeling of hydraulic fracture branching in heterogeneous formation under anisotropic in-situ stress