Continuum Porous-Medium CFD Modelling of Rock-Bed Thermal Energy Storage Systems: A Review of Pressure-Drop and Interphase Heat-Transfer Correlations

Rock-bed thermal energy storage (RTES) systems are attracting growing interest as low-cost, robust, and scalable sensible heat storage solutions for applications ranging from low-temperature building and greenhouse heating to medium- and high-temperature solar or waste-heat recovery systems. However, their thermo-hydraulic performance is strongly influenced by the complex interactions among heat-transfer-fluid flow, irregular rock morphology, porosity, pressure drop, interphase heat transfer, and transient thermal-front development. This review provides a focused evaluation of computational fluid dynamics (CFD) modelling strategies for packed beds of rocks, with particular attention to continuum porous-medium approaches and the closure correlations required for reliable simulation. First, the distinction between pore-scale and volume-averaged continuum modelling is discussed in terms of the trade-off between physical resolution and computational feasibility. The main pressure-drop and friction-factor correlations are then reviewed and compared, including classical packed-bed models and rock-bed-specific formulations. It is shown that hydraulic-resistance predictions are highly sensitive to particle shape, surface roughness, porosity, the bed-to-particle diameter ratio, and packing arrangement. Particle-fluid heat-transfer correlations are also examined and, when possible, converted into a consistent particle Nusselt-number form to enable direct comparison. Particular attention is given to generalized correlations, dispersion-corrected models, and air–rock-bed correlations applicable to thermal storage systems. Finally, a methodological framework for modelling RTES systems using local thermal equilibrium (LTE) and local thermal non-equilibrium (LTNE) formulations is proposed. Dimensionless criteria, including the interphase thermal coupling number and particle Biot number, are introduced to support the selection between LTE and LTNE formulations. The selection of pressure-drop/friction-factor and solid–fluid heat-transfer/particle Nusselt-number correlations should be based on the similarity between the original experimental conditions and the target RTES system, and system-specific validation is recommended whenever possible.