“…Yet, explicitly accounting for fractures in numerical models of poromechanical problems is a challenging task and often comes at the expenses of a very high computational cost (Lei et al., 2017; McDermott & Kolditz, 2006; Pandey et al., 2017; Salimzadeh et al., 2018; Thomas et al., 2020; Yoshioka et al., 2019). In order to facilitate the calculation for large problems, dense fractured networks are often simplified as equivalent porous media, like the ubiquitous joint model, and computational resources are used for other purposes, such as reactive transport, phase changes, and dynamics or complex plasticity models (e.g., Birdsell et al., 2018; Liu et al., 2009; Parisio, Vilarrasa, et al., 2019; Rinaldi et al., 2015; Rutqvist, Wu, et al., 2002; Rutqvist et al., 2005; Taron & Elsworth, 2009; Vallier et al., 2020; Y. Wang et al., 2019).…”