Fractured rocks are the host for many engineered structures related to energy, water, waste, and transportation. For this reason, the functioning of such infrastructure, as well as optimization of its engineering efficiency, critically depends on characterizing, modeling, and monitoring fractured rock sites (National Academies of Sciences, 2015). Yet, modeling fluid flow in fractured porous media is a longstanding challenge due to the extensive numerical discretization and intractable computational cost imposed by discrete fractures and the matrix they are embedded within, where length scales span millimeter-scale fractures to the km-scale regions of the surrounding rock, all while attempting to achieve high model realism and accuracy for the correct management and reliable prediction of hydraulic and thermal behaviors of fractured rock masses. Not surprisingly, there has been a surge of interest on this research problem and it still remains interesting to many researchers due to its importance in a wide scope of applications such as geologic radioactive waste disposal (e.g.,