Effects of Wall Roughness on Coupled Flow and Heat Transport in Fractured Media

Abstract: Heat transfer in fractured media is governed by the interplay between advective transport along rough-walled fractures and conductive transport, both within the fractures and in the surrounding low-permeability matrix. Flow localization induced by aperture heterogeneity, combined with matrix conduction, gives rise to anomalous thermal behavior. To capture these effects, we develop a stochastic modeling framework that couples a time-domain random walk (TDRW) representation of advective and conductive transport in the fractures with a semi-analytical model of conductive heat exchange with the matrix. Matrix trapping times follow a L\'evy-Smirnov distribution derived from first-passage theory, capturing the heavy-tailed dynamics typical of fractured systems. Heat flux at the fracture-matrix interface is computed via a nonlocal convolution integral based on Duhamel's principle, accounting for thermal memory effects. The model is validated against analytical benchmarks and finite-element simulations. Monte Carlo simulations over stochastic aperture fields quantify the influence of fracture closure, correlation length, and P\'eclet number. Results reveal a transition from superdiffusive to subdiffusive regimes, driven by the competition between advective transport along preferential paths, dispersion induced by aperture variability, and matrix-driven heat conduction. In the long-time regime, heat exchange exhibits a characteristic $t^{-1/2}$ decay. At early times, limited thermal penetration into the matrix leads to weaker interfacial fluxes, underscoring the role of matrix thermal inertia. The proposed framework enables physically consistent and computationally efficient simulations of thermal transport in complex fractured systems, with implications for geothermal energy, subsurface thermal storage, and engineered heat exchange in low-permeability environments.