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

Alessandro Lenci; Daniel M. Tartakovsky; Maria Klepikova; Vittorio Di Federico; Yves M\'eheust

arxiv: 2510.07294 · v2 · submitted 2025-10-08 · ⚛️ physics.geo-ph · physics.flu-dyn

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

Alessandro Lenci , Yves M\'eheust , Maria Klepikova , Vittorio Di Federico , Daniel M. Tartakovsky This is my paper

Pith reviewed 2026-05-18 09:21 UTC · model grok-4.3

classification ⚛️ physics.geo-ph physics.flu-dyn

keywords fractured mediaheat transportwall roughnessanomalous diffusiontime-domain random walkmatrix conductiongeothermal energystochastic modeling

0 comments

The pith

A stochastic framework reveals that fracture wall roughness drives a shift from superdiffusive to subdiffusive thermal transport.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a modeling framework that links advective flow in rough fractures with conductive heat exchange to the rock matrix using random walks and a nonlocal integral. This captures how aperture variations create preferential paths and how matrix conduction affects overall heat movement. Readers interested in geothermal systems or subsurface storage would care because the model shows efficient ways to simulate these complex interactions without full numerical grids. It demonstrates a regime transition in diffusion behavior and characteristic decay rates for heat flux. The approach is validated against benchmarks and used to explore effects of different fracture properties.

Core claim

The central claim is that coupling a time-domain random walk model for transport within rough-walled fractures to a semi-analytical convolution model for heat exchange with the matrix, based on Levy-Smirnov trapping times, enables simulation of anomalous thermal transport. This reveals a transition from superdiffusive to subdiffusive regimes due to the competition between advective flow in preferential paths, dispersion from aperture variability, and matrix-driven conduction. In the long-time limit, interfacial heat exchange decays as t to the power of negative one-half.

What carries the argument

Time-domain random walk (TDRW) for fracture transport coupled with a nonlocal convolution integral for matrix heat exchange derived from first-passage theory and Duhamel's principle.

Load-bearing premise

The load-bearing premise is that matrix trapping times follow a Levy-Smirnov distribution from first-passage theory, which allows the semi-analytical nonlocal model to describe conductive heat exchange.

What would settle it

A mismatch between the model's predicted regime transition and direct numerical simulations using finite elements on the same stochastic aperture field would falsify the claim that the framework accurately captures the superdiffusive to subdiffusive shift.

read the original 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.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper couples TDRW fracture transport with a Lévy-Smirnov matrix convolution to explore roughness effects on heat diffusion regimes, with decent validation but a modeling assumption worth checking.

read the letter

The main point is a stochastic framework that links time-domain random walks for advective-conductive transport inside rough fractures to a nonlocal Duhamel convolution for heat exchange with the matrix, where trapping times follow a Lévy-Smirnov distribution from first-passage theory. Monte Carlo runs over random aperture fields then map out how closure, correlation length, and Péclet number push the system from superdiffusive to subdiffusive behavior, with the expected long-time t to the minus one-half flux decay showing up from matrix conduction. They back the setup with comparisons to analytical solutions and finite-element simulations, which keeps the claims grounded and shows the method is faster than full-grid numerics for these problems. That combination of efficiency and parameter exploration is the useful part for applied work on geothermal or thermal storage systems. The soft spot is the fixed Lévy-Smirnov kernel for matrix interaction. It comes from a flat-interface, homogeneous-matrix assumption, and when aperture varies spatially the local contact area and diffusion lengths change, which could alter early-time fluxes or shift the exact balance that produces the reported regime transition. Their overall validation probably keeps the discrepancy small, but it would help to see a direct test of that sensitivity rather than leaving it implicit. This is aimed at researchers who model subsurface heat flow and already work with random-walk or first-passage methods. A reader looking for a practical way to run many realizations of rough-fracture thermal transport without heavy computation will find value here. The methods are coherent enough on their own terms that the paper deserves a serious referee to examine the coupling details and Monte Carlo statistics. I would send it for peer review.

Referee Report

1 major / 2 minor

Summary. The manuscript develops a stochastic modeling framework coupling a time-domain random walk representation of advective and conductive transport in rough-walled fractures with a semi-analytical nonlocal convolution model for matrix heat exchange. Matrix trapping times are drawn from a Lévy-Smirnov distribution obtained via first-passage theory; heat flux at the interface is computed with Duhamel's principle. Monte Carlo simulations over stochastic aperture fields quantify the roles of fracture closure, correlation length, and Péclet number, revealing a transition from superdiffusive to subdiffusive thermal regimes driven by preferential advection, aperture-induced dispersion, and matrix conduction. The model is validated against analytical benchmarks and finite-element simulations, with long-time interfacial flux exhibiting a t^{-1/2} decay.

Significance. If the modeling assumptions remain valid under roughness, the framework supplies a computationally efficient route to simulate anomalous heat transport in fractured media. The Monte Carlo quantification of parameter effects and the identification of competing mechanisms that produce regime transitions would be of direct interest for geothermal energy, thermal storage, and low-permeability heat-exchange applications.

major comments (1)

[Model description (semi-analytical matrix exchange section)] Model description (semi-analytical matrix exchange section): The nonlocal convolution kernel is fixed by the Lévy-Smirnov first-passage distribution derived for a homogeneous semi-infinite matrix bounded by a locally flat interface. With spatially varying aperture (finite correlation length and roughness), local contact area, effective diffusion length, and instantaneous heat-flux boundary conditions become position-dependent; the global kernel therefore omits a heterogeneous memory effect that can quantitatively alter early-time flux decay and the long-time t^{-1/2} tail, thereby changing the balance that produces the reported super-to-subdiffusive transition.

minor comments (2)

[Abstract] Abstract: The statement that results reveal a transition 'driven by the competition between advective transport along preferential paths, dispersion induced by aperture variability, and matrix-driven heat conduction' should be cross-referenced to the specific Monte Carlo diagnostics (e.g., mean-squared displacement scaling or effective diffusivity) that isolate each contribution.
[Validation section] Validation section: While analytical benchmarks and finite-element comparisons are mentioned, explicit error norms, grid-convergence data, or the precise definition of the Péclet number used in the coupled problem would strengthen reproducibility.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive review and positive overall assessment of our manuscript. We address the major comment on the semi-analytical matrix exchange model below, offering clarifications supported by our validation results and indicating the revisions we will make.

read point-by-point responses

Referee: The nonlocal convolution kernel is fixed by the Lévy-Smirnov first-passage distribution derived for a homogeneous semi-infinite matrix bounded by a locally flat interface. With spatially varying aperture (finite correlation length and roughness), local contact area, effective diffusion length, and instantaneous heat-flux boundary conditions become position-dependent; the global kernel therefore omits a heterogeneous memory effect that can quantitatively alter early-time flux decay and the long-time t^{-1/2} tail, thereby changing the balance that produces the reported super-to-subdiffusive transition.

Authors: We appreciate the referee highlighting this modeling assumption. Our semi-analytical matrix exchange indeed employs the Lévy-Smirnov first-passage distribution and Duhamel convolution under the approximation of a homogeneous semi-infinite matrix with a locally flat interface. Spatial variations arising from aperture heterogeneity and wall roughness are incorporated via the TDRW discretization of the fracture, where local aperture values determine advective velocities, particle paths, and interface locations on a realization-by-realization basis. The nonlocal kernel is then applied to compute ensemble-averaged interfacial fluxes. This effective treatment is justified by direct validation against finite-element simulations that explicitly resolve the full rough geometry, position-dependent conduction lengths, and heterogeneous boundary conditions. The quantitative agreement between the stochastic model and these resolved simulations—including reproduction of the characteristic t^{-1/2} long-time decay and the superdiffusive-to-subdiffusive regime transition—indicates that the leading-order effects of matrix memory and competing transport mechanisms are captured. We acknowledge that a fully position-dependent heterogeneous kernel could refine early-time flux predictions. In the revised manuscript we will expand the model description section to explicitly discuss the homogeneous-matrix approximation, its range of validity, and the supporting evidence from the FEM benchmarks. revision: partial

Circularity Check

0 steps flagged

No circularity: derivations apply external first-passage theory and Duhamel's principle to rough-fracture problem

full rationale

The paper couples a TDRW representation of advective-conductive transport in fractures to a semi-analytical matrix model whose trapping-time kernel is taken from Lévy-Smirnov first-passage theory and Duhamel's principle. These are standard external results applied to the heterogeneous-aperture setting; Monte Carlo results over stochastic fields then quantify the super-to-subdiffusive transition. No equation or claim in the abstract or described framework reduces a reported prediction or regime transition to a quantity fitted inside the same study or to a self-citation chain that itself lacks independent verification. The central outputs therefore remain independent of the paper's own inputs.

Axiom & Free-Parameter Ledger

3 free parameters · 1 axioms · 0 invented entities

The central modeling approach rests on standard stochastic transport theory and first-passage results rather than new fitted constants or invented physical entities; the Lévy-Smirnov distribution is treated as a derived input from prior theory.

free parameters (3)

fracture closure
Varied parametrically in Monte Carlo runs to quantify influence; not fitted to match a specific target data set within the study.
correlation length
Varied parametrically in Monte Carlo runs to quantify influence; not fitted to match a specific target data set within the study.
Péclet number
Varied parametrically in Monte Carlo runs to quantify influence; not fitted to match a specific target data set within the study.

axioms (1)

domain assumption Matrix trapping times follow a Lévy-Smirnov distribution derived from first-passage theory.
Invoked to capture heavy-tailed dynamics typical of fractured systems and to enable the semi-analytical heat-exchange model.

pith-pipeline@v0.9.0 · 5842 in / 1542 out tokens · 43416 ms · 2026-05-18T09:21:32.869052+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Matrix trapping times follow a Lévy-Smirnov distribution derived from first-passage theory... Heat flux... via a nonlocal convolution integral based on Duhamel’s principle
IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

Results reveal a transition from superdiffusive to subdiffusive regimes, driven by... matrix-driven heat conduction... t^{-1/2} decay

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.