19–22 May 2025
US/Mountain timezone

Laminar to Turbulent Convection in Porous Media: The Role of Solid-Fluid Conductivity Ratios and Porosity Variations

22 May 2025, 09:05
15m
Oral Presentation (MS09) Pore-scale modelling MS09

Speaker

Dario Schwendener (ETH Zurich)

Description

Transition from laminar to turbulent flow within the framework of conjugate heat transfer in porous media occurs in various applications across different scales, including geothermal energy extraction, thermal energy storage systems, high-temperature gas-cooled reactors, and microchip cooling. Understanding this transition is critical for optimizing systems designed to maximize surface area for efficient heat and mass transfer.

To tackle the regime transition, we developed an advanced lattice-Boltzmann Method (LBM) solver optimized for simulating conjugate heat transfer in porous media. Built on the STLBM open-source platform, the solver integrates thermo-physical heterogeneity and supports multi-threaded, parallelized computations on both CPUs and GPUs. It accurately captures the Navier-Stokes–Fourier dynamics under the Boussinesq approximation, providing a robust framework for analyzing heat and fluid flow in porous structures, including inertial effects and thermophysical heterogeneity.

Using staggered isotropic porous media composed of cylindrical structures, we investigate the effects of porosity, solid-to-fluid conductivity ratio, and Rayleigh number on overall dynamics, including Nusselt number, Reynolds number, and boundary layer thickness. In pure hydrodynamics, increasingly confined pore space strongly influences the transition from Darcy to non-Darcy flow as porosity decreases from 45% to 30%. Extending this analysis, we explore natural convection behavior during the transition from Darcy to non-Darcy (Forchheimer) flow, focusing on Darcy numbers around 10⁻⁶. Our simulations further assess the impact of conductivity ratios (0.1 to 10) on convective dynamics across Rayleigh numbers spanning four orders of magnitude (10⁷–10¹⁰).

Our results reveal that inertial forces drive regime transitions from steady-state to oscillatory convection, as evidenced by spatial Reynolds number analysis. Nusselt–Rayleigh scaling in steady-state convection aligns with the classic Nu ~ Ra at low porosities but deviates significantly at higher porosities. For low porosities, the transition to oscillatory convection is strongly influenced by the solid-to-fluid conductivity ratio. At high Rayleigh numbers, higher kinetic energy does not necessarily enhance heat transfer, as boundary layer thickness is influenced by both velocity and the thermal conductivity of the solid and fluid phases. Furthermore, within the examined Darcy numbers, the pore-scale Prandtl number fails to reliably predict the transition to the Forchheimer regime.

References Wood, B. D., He, X., & Apte, S. V. (2020). Modeling turbulent flows in porous media. Annual Review of Fluid Mechanics, 52(1), 171–203. Published by Annual Reviews. | Latt, J., Coreixas, C., & Beny, J. (2021). Cross-platform programming model for many-core lattice Boltzmann simulations. PLOS ONE, 16(4), 1–29. DOI: 10.1371/journal.pone.0250306. Published by the Public Library of Science. | Khalifa, Z., Pocher, L., & Tilton, N. (2020). Regimes of flow through cylinder arrays subject to steady pressure gradients. International Journal of Heat and Mass Transfer, 159, 120072. Published by Elsevier. | Korba, D., & Li, L. (2022). Effects of pore scale and conjugate heat transfer on thermal convection in porous media. Journal of Fluid Mechanics, 944, A28. Published by Cambridge University Press.
Country Switzerland
Acceptance of the Terms & Conditions Click here to agree

Primary author

Dario Schwendener (ETH Zurich)

Co-authors

Dr Xiangzhao Kong Dr Jerome Noir (ETH Zurich) Dr Christophe Coreixas (University of Geneva) Dr Jonas Latt (University of Geneva)

Presentation materials

There are no materials yet.