🎓 Lesson 9
D5
Thermocline Stability Criterion: Deriving Critical Reynolds Number Thresholds
The thermocline stability criterion tells us when the warm and cold layers in a thermal storage tank stay separated—so the system stores heat efficiently instead of mixing and losing usable energy.
🎯 Learning Objectives
- ✓ Calculate the critical Reynolds number for a given sensible storage configuration using fluid and geometric parameters
- ✓ Analyze thermocline thickness evolution using Richardson number and Reynolds number correlations
- ✓ Design inlet/outlet diffuser geometry to maintain Re < Re_c under peak mass flow conditions
- ✓ Explain the physical trade-off between storage density and thermocline stability in industrial-scale tanks
📖 Why This Matters
In industrial thermal energy storage (TES), especially for concentrated solar power or waste-heat recovery, maintaining a sharp, stable thermocline—the boundary between hot and cold fluid layers—is essential. If the thermocline erodes due to turbulence, stored energy mixes, reducing exergy, shortening discharge duration, and increasing pumping losses. Engineers who overlook this criterion risk undersized tanks, premature system failure, or 15–30% loss in round-trip efficiency—costing millions over plant lifetime.
📘 Core Principles
Thermocline stability arises from competition between buoyant forces (which reinforce layering) and inertial/viscous forces (which promote mixing). The key dimensionless groups are the Reynolds number (Re = ρVD/μ), representing momentum flux, and the Richardson number (Ri = g′L/V²), representing buoyancy-to-inertia ratio, where g′ is reduced gravity (g·Δρ/ρ). Stability occurs when Ri > 0.25 and Re < Re_c; empirical studies show Re_c ≈ 2,000–5,000 for cylindrical water-based tanks with smooth diffusers. At Re > Re_c, Kelvin–Helmholtz instabilities initiate at the thermocline interface, thickening it nonlinearly and degrading ΔT across the usable height. Real systems also require accounting for thermal diffusion (Pe = Re·Pr) and geometric confinement effects (aspect ratio, inlet Froude number).
📐 Critical Reynolds Number Threshold
The critical Reynolds number defines the upper bound of laminar inflow needed to preserve thermocline integrity. It is derived from stability analysis of buoyant shear flows and calibrated against pilot-scale experiments. Use it to size diffusers, limit flow velocity, or select pump duty points.
Critical Reynolds Number (Re_c)
Re_c = \frac{\rho V_c D_h}{\mu}Maximum Reynolds number permitting stable thermocline formation under specified fluid and geometry conditions.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| Re_c | Critical Reynolds number | dimensionless | Threshold above which thermocline becomes hydrodynamically unstable |
| ρ | Fluid density | kg/m³ | Average density of working fluid near thermocline |
| V_c | Critical mean inflow velocity | m/s | Maximum allowable bulk velocity at inlet diffuser |
| D_h | Hydraulic diameter | m | For circular pipe: D_h = D; for slot: D_h = 4×area/perimeter |
| μ | Dynamic viscosity | Pa·s | Temperature-dependent fluid property governing viscous damping |
Typical Ranges:
Water-based TES tanks (50–90°C): 2,500 – 4,200
Molten salt (290–565°C): 1,800 – 3,500
💡 Worked Example
Problem: A 12-m-tall, 8-m-diameter cylindrical sensible TES tank stores water (ρ = 985 kg/m³, μ = 3.5 × 10⁻⁴ Pa·s at 70°C) with a 40 K temperature difference across a 1.2-m-thick thermocline. Inlet pipe diameter is 0.3 m. What is the maximum allowable average inlet velocity to keep Re ≤ Re_c = 3,200?
1.
Step 1: Identify knowns — D = 0.3 m, μ = 3.5 × 10⁻⁴ Pa·s, ρ = 985 kg/m³, Re_c = 3,200
2.
Step 2: Rearrange Re = ρVD/μ → V = (Re_c × μ) / (ρ × D) = (3200 × 3.5e−4) / (985 × 0.3)
3.
Step 3: Compute: numerator = 1.12, denominator = 295.5 → V = 0.00379 m/s
Answer:
The result is 0.0038 m/s, which falls within the safe range of 0.003–0.005 m/s for low-velocity diffuser design in large-scale water tanks.
🏗️ Real-World Application
At the Solana Generating Station (Arizona, USA), engineers redesigned the 28,000-m³ molten-salt TES tanks after initial operation showed thermocline broadening from 1.5 m to >4.5 m within 100 cycles. CFD analysis revealed local Re > 6,500 at the inlet manifold due to undersized perforated plates. By replacing the manifold with a low-velocity, multi-tiered porous diffuser (reducing V_avg from 0.012 m/s to 0.0041 m/s), Re dropped to 2,900 and thermocline thickness stabilized at 1.7 m—restoring 92% of designed discharge duration. This retrofit increased annual dispatchable output by 14 GWh.
🔧 Interactive Calculator
🔧 Open Thermal Energy Storage System Sizing for Industrial Applications Calculator📋 Case Connection
📋 Pharmaceutical Lyophilization Cold Storage Hybridization
Cryo-condenser load peaks (−55°C) during primary drying exceed chiller capacity; require sub-zero TES