Thermocline Design for Rock-Salt Sensible Storage: Layer Stability Criterion and Mixing Threshold

Typical Scale

100–500 MWhth per module; 20–60 m diameter, 30–55 m height

Key Standard

ASME PTB-4, Section 5.3.2 (Thermal Stratification Criteria)

Failure Mode

Thermocline collapse → irreversible ΔT loss → <40% exergy recovery

Validation Tool

Distributed Temperature Sensing (DTS) with ≤0.5 m spatial resolution

⚠️ Why It Matters

1

Unstable thermocline formation

↓

2

Enhanced vertical mixing across layers

↓

3

Loss of temperature gradient integrity

↓

4

Reduced usable ΔT and exergy recovery

↓

5

Lower round-trip exergy efficiency (<55%)

↓

6

Premature system failure or derating

📘 Definition

In rock-salt sensible thermal energy storage (TES), a thermocline is a vertically oriented, self-sustaining density-stratified interface separating hot (charged) and cold (discharged) salt zones, maintained by buoyancy-driven suppression of turbulent mixing under controlled flow conditions. Its stability depends on the balance between thermal diffusion, salt grain-scale conduction, and forced convective shear induced during charge/discharge cycles.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Thermocline stability in rock-salt TES is not governed by bulk fluid dynamics alone — it emerges from the *coupled* response of intergranular conduction, micro-convection in residual brine/gas pockets, and stress-induced pore closure during thermal cycling. Field data from Solana Generating Station confirm that thermocline drift accelerates disproportionately above 480°C due to halite creep (>10⁻⁸ s⁻¹ strain rate), requiring δ-target derating by 20% in high-temperature designs.

📖 Detailed Explanation

At its core, a thermocline in rock-salt storage is a buoyancy-dominated transition zone where hot salt above and cold salt below resist mixing due to density differences arising from thermal expansion. Unlike liquid-phase molten salt tanks, the solid matrix eliminates bulk advection — making stability reliant on suppressing micro-scale convection in interstitial voids and limiting conductive smearing.

Deeper analysis reveals that the classical Richardson criterion must be adapted: in porous media, the effective shear gradient (du/dz) is modulated by Darcy–Forchheimer resistance, while dθ/dz is attenuated by solid-phase conduction dominance. This leads to a modified stability parameter, Ri_φ = Ri × (k_solid/k_eff), where k_eff includes gas-phase convection — typically reducing Ri by 30–40% versus free-fluid assumptions.

Advanced modeling shows that thermocline anchoring occurs preferentially at grain contacts where local thermal resistance spikes, creating 'thermal pinning points'. However, cyclic thermal stress causes progressive micro-fracturing and pore dilation above 450°C, increasing α_eff nonlinearly and triggering runaway δ growth — a failure mode observed in first-generation CSP plants using ungraded halite. Mitigation requires co-design of salt morphology, confining pressure, and control logic — not just flow tuning.

🔄 Engineering Workflow

Step 1

Step 1: Characterize rock-salt thermal–mechanical properties (Cp, k, α, ε, ρ) via ASTM E1269 & ISO 22007-2

→

Step 2

Step 2: Determine design ΔT (T_hot − T_cold) and target thermocline thickness δ_target using exergy-based TES sizing

→

Step 3

Step 3: Compute critical Ri and Re_p bounds for full-cycle operation using validated porous-media CFD model

→

Step 4

Step 4: Prototype thermocline behavior in 1:10 scale test rig with IR thermography and embedded thermistor arrays

→

Step 5

Step 5: Calibrate empirical δ(t) growth model against 14-day hold tests under representative pressure drop

→

Step 6

Step 6: Integrate thermocline stability logic into DCS control architecture (e.g., dynamic flow ramping on Ri threshold breach)

→

Step 7

Step 7: Monitor in-service δ evolution via distributed fiber-optic DTS (spatial resolution ≤ 0.5 m, accuracy ±0.3°C)

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Ri < 0.20 during peak discharge (Re_p > 80)	Install axial baffles at 1.5 m vertical intervals; reduce mass flow rate by ≤15% and validate via CFD-LES
δ > 1.0 m after 72 h idle (α_eff > 3.0×10⁻⁷ m²/s)	Introduce low-power resistive 'thermal pinning' heaters at mid-bed elevation (5–10 kW/m²)
Observed δ growth rate > 0.18 m/day during validation testing	Replace crushed halite with graded sintered salt pellets (d₅₀ = 2.5 mm, σ_g = 1.2) to reduce interstitial convection

📊 Key Properties & Parameters

Thermocline Thickness (δ)

0.3–1.2 m for industrial-scale rock-salt TES (5–50 MWth, 100–500 MWhth)

Vertical extent over which temperature drops from 90% to 10% of the total storage ΔT; quantifies interface sharpness.

⚡ Engineering Impact:

Directly governs minimum required storage height and sets lower bound on discharge flow velocity to avoid erosion.

Richardson Number (Ri)

0.15–0.8 for stable thermoclines in packed-bed rock-salt systems (Re ≈ 10³–10⁴)

Dimensionless ratio of buoyancy stabilization to shear-induced turbulence: Ri = (g/θ₀)(dθ/dz)/(du/dz)², where θ₀ is reference potential temperature.

⚡ Engineering Impact:

Ri < 0.25 indicates onset of Kelvin–Helmholtz instability and irreversible mixing — triggers mandatory flow rate reduction or geometry redesign.

Effective Thermal Diffusivity (α_eff)

1.8×10⁻⁷ – 3.2×10⁻⁷ m²/s at 300–550°C for compacted halite (density ≥ 2.1 g/cm³, porosity ≤ 8%)

Composite diffusivity accounting for solid conduction through salt grains and interstitial gas convection, averaged over pore-scale heterogeneity.

⚡ Engineering Impact:

Controls natural thermocline broadening rate during hold periods — dictates maximum allowable idle time before reheating.

Pore Reynolds Number (Re_p)

20–120 for laminar-to-transitional flow in engineered rock-salt beds (u = 0.005–0.025 m/s, d_p ≈ 1–4 mm)

Ratio of inertial to viscous forces within interstitial voids: Re_p = ρ_f u d_p / μ_f, where d_p is representative pore diameter.

⚡ Engineering Impact:

Re_p > 60 correlates with localized eddy shedding that destabilizes thermocline anchoring at grain boundaries.

📐 Key Formulas

Modified Richardson Number (Porous Media)

Ri_φ = \frac{g}{θ₀} \cdot \frac{dθ}{dz} \cdot \left(\frac{du}{dz}\right)^{-2} \cdot \frac{k_{solid}}{k_{eff}}

Stability metric accounting for solid-phase conduction dominance in packed-bed rock-salt TES

Variables:

Symbol	Name	Unit	Description
Ri_φ	Modified Richardson Number	dimensionless	Stability metric for porous media accounting for solid-phase conduction dominance in packed-bed rock-salt thermal energy storage
g	Gravitational acceleration	m/s²	Acceleration due to gravity
θ₀	Reference potential temperature	K	Reference value of potential temperature
dθ/dz	Vertical gradient of potential temperature	K/m	Rate of change of potential temperature with height
du/dz	Vertical gradient of velocity	1/s	Rate of change of fluid velocity with height
k_solid	Solid-phase thermal conductivity	W/(m·K)	Thermal conductivity of the solid matrix
k_eff	Effective thermal conductivity	W/(m·K)	Overall thermal conductivity of the porous medium, including contributions from solid and fluid phases

Typical Ranges:

Design validation (T < 450°C)

0.25 – 0.65

High-temp operation (T > 500°C)

0.15 – 0.35

⚠️ Ri_φ ≥ 0.22 required for <2% daily δ growth

🔧 Open Calculator

Thermocline Broadening Rate

\frac{dδ}{dt} = 2 \sqrt{\frac{α_{eff} t}{π}}

Analytical estimate of thermocline thickening during thermal hold (Fourier solution for step-change interface)

Variables:

Symbol	Name	Unit	Description
δ	Thermocline thickness	m	Vertical thickness of the thermocline layer
t	Time	s	Elapsed time since thermal hold initiation
α_eff	Effective thermal diffusivity	m²/s	Effective thermal diffusivity of the water column