“…In addition, the assumptions that we considered in the current work for the sake of simplicity (such as linear, elastic, isotropic, and continuum media) or may result in non‐accurate predictions for the real hydraulic fracture conditions. For such cases, some better numerical methods such as finite element with continuum‐discontinuum elements, 83 cohesive zone modeling, 84 extended finite element (XFEM), 85 displacement discontinuity/boundary collocation method, 86 conjugated bond‐pair‐based peridynamics 87–89 and the nonlocal general particle dynamics 90–93 may probably result in better estimations and predictions for the onset of hydraulic fracture, or the trajectory and coalescence of natural fractures around the hydraulic fracture. For example, Zhou et al, 88 presented hydromechanical bond based peridynamic model for investigating the cracking phenomenon in porous rocks containing single or multiple fissures.…”