Laboratory heat transport experiments reveal grain size and flow velocity dependent local thermal non-equilibrium effects

Abstract. Heat transport in porous media is crucial for gaining earth science process understanding and engineering applications such as geothermal system design. While heat transport models are commonly simplified by assuming local thermal equilibrium (LTE, solid and fluid phases are averaged), local thermal non-equilibrium (LTNE, solid and fluid phases are considered separately) heat transport has long been hypothesized and reports have emerged. However, experiments with realistic grain sizes and flow conditions are still lacking in the literature. To detect LTNE effects, we conducted comprehensive laboratory heat transport experiments at Darcy velocities ranging from 3 to 23 m d⁻¹ and measured the temperature of fluid and solid phases separately for glass spheres with diameters of 5, 10, 15, 20, 25 and 30 mm. Four replicas of each size were embedded at discrete distances along the flow path in small glass beads to stabilize the flow field. Our sensors were meticulously calibrated and measurements were post-processed to reveal LTNE, expressed as the difference between solid and fluid temperature during the passing of a thermal step input. To gain insight into the heat transport properties and processes, we simulated our experimental results in 1D using commonly accepted analytical solutions for LTE and a numerical solution of LTNE equations. Our results demonstrate significant LTNE effects with increasing grain size and water flow velocity. Surprisingly, some temperature differences were negative indicating that the heat front propagates non-uniformly likely caused by spatial variations of the flow field. The fluid temperature modeled by the LTE analytical solution exhibited relatively good agreement with experimental fluid temperature only for grain sizes from 5 mm to 15 mm. However, for larger grain sizes (between 20 mm and 30 mm), the temperature difference between fluid and solid phases became too significant to be represented by an LTE model. Additionally, for larger grain sizes (≥ 20 mm), the LTNE model failed to predict the magnitude of LTNE for all tested flow velocities due to experimental conditions being inadequately represented by the 1D model with ideal step input. Future studies should employ more sophisticated numerical models to examine the heat transport processes and accurately analyze LTNE effects, considering non-uniform flow effects and multi-dimensional solution. This is essential to determine the validity limits of LTE conditions for heat transport in natural systems such as gravel aquifers with grain sizes larger than 20 mm.