Calculator D4

TES Tank Geometry Optimization: Aspect Ratio, Baffle Placement, and Natural Convection Suppression

Optimizing the shape and internal structure of a thermal energy storage (TES) tank helps heat move efficiently during charging and discharging — like designing a thermos that keeps coffee hot longer by controlling how the liquid swirls inside.

Typical Scale
Industrial TES tanks: 5–50 m diameter, 10–40 m height, 500–20,000 m³ volume
Key Standards
ASHRAE Guideline 36-2021, IEC 62788-6-1:2022, ASTM E2993-21
Industry Applications
CSP plants, industrial decarbonization (steel, cement), grid-scale long-duration storage

⚠️ Why It Matters

1
Excessive vertical convection in tall narrow tanks
2
Thermal short-circuiting between hot and cold zones
3
Reduced effective storage capacity (<75% usable fraction)
4
Mismatch between discharge rate and process heat demand
5
Increased parasitic pumping energy and exergy destruction
6
Premature system degradation due to thermal cycling fatigue

📘 Definition

TES tank geometry optimization is the systematic engineering design process that determines optimal aspect ratio (height-to-diameter), baffle configuration, and internal flow suppression features to minimize natural convection-driven thermal stratification loss, maximize usable energy fraction, and ensure charge/discharge rate compatibility with process heat demand profiles. It integrates fluid mechanics, transient heat transfer, and exergy-based performance validation for molten salt, PCM, and sensible TES systems operating at 200–700 °C.

🎨 Concept Diagram

HotColdH/D = ?Baffles ↓

AI-generated illustration for visual understanding

💡 Engineering Insight

Baffles are not passive 'mixing inhibitors'—they are active thermal front stabilizers. Their optimal placement aligns with the predicted plume impingement height (≈0.62·H for Ra ~ 10¹⁰), not mid-height. Over-baffling increases pumping loss more than it improves stratification—and introduces weld-seam corrosion hotspots at support brackets.

📖 Detailed Explanation

Thermal energy storage tanks rely on stable vertical temperature gradients to store energy efficiently. When hot fluid enters at the top and cold fluid exits at the bottom (or vice versa), density differences drive natural convection—causing hot and cold layers to mix, reducing usable energy. Geometry sets the stage: a tall, narrow tank amplifies vertical motion; a short, wide one promotes lateral dispersion but increases surface-area-to-volume ratio and conductive losses.

The Rayleigh number (Ra) quantifies this tendency: for molten salts like Solar Salt (60% NaNO₃ + 40% KNO₃), Ra exceeds 10⁹ even at modest ΔT (50 K) and H = 20 m—guaranteeing turbulent convection. Baffles interrupt rising thermal plumes by forcing fluid through constrained apertures, converting kinetic energy into dissipation. However, their effectiveness depends on porosity, edge geometry, and vertical spacing—not just count or height.

Advanced optimization requires coupling conjugate heat transfer with turbulence modeling (e.g., LES or k-ω SST with buoyancy corrections) and tracking exergy destruction density fields. Real-world constraints—such as weld accessibility, thermal expansion mismatch between baffle and shell, and fouling accumulation on baffle undersides—must be embedded in the CFD boundary conditions. The most robust designs use adaptive baffle heights (tapered upward) and incorporate sacrificial anode integration points at baffle-to-shell welds per NACE SP0169-2022.

🔄 Engineering Workflow

Step 1
Step 1: Define operational envelope (T_min/T_max, Δt_charge/Δt_discharge, Q_dot_process, duty cycle)
Step 2
Step 2: Compute dimensionless groups (Ra, Pr, Re_inlet, TSI_target) using worst-case ΔT and fluid properties
Step 3
Step 3: Generate parametric CFD ensemble (H/D ∈ [0.6, 3.0], baffle count ∈ [0, 5], baffle porosity ∈ [0.1, 0.6])
Step 4
Step 4: Extract usable energy fraction (UEF) and exergy efficiency (η_ex) from transient simulations across 10+ cycles
Step 5
Step 5: Validate against ASHRAE Guideline 36 Annex G thermal performance test protocol
Step 6
Step 6: Iterate geometry + insulation + nozzle design to meet UEF ≥ 0.85 and η_ex ≥ 0.72 at design point
Step 7
Step 7: Conduct scaled physical testing (1:10 water-glycerol analog) to confirm baffle-induced flow damping

📋 Decision Guide

Rock/Field Condition Recommended Design Action
High-temperature molten salt (e.g., Solar Salt, 565 °C) with ΔT ≥ 200 K and H/D > 2.2 Install stepped-height baffles at h_b/H = 0.25 and 0.75; limit max H/D to 2.0; add radial flow guides near inlet/outlet nozzles
Encapsulated PCM (paraffin, 60–80 °C) with low Prandtl number (Pr ≈ 10) and slow discharge rate (<1 kW/m²) Use H/D = 1.0–1.4; install perforated vertical baffles (25% open area) to dampen plume coalescence without impeding phase change front propagation
Sensible water/glycol TES (<120 °C) with rapid cycling (≤2 hr charge/discharge) and Ra < 10⁸ Prioritize low H/D (0.8–1.2); omit baffles; rely on controlled inlet velocity (0.05–0.15 m/s) and diffuser geometry for laminar thermal front advancement

📊 Key Properties & Parameters

Aspect Ratio (H/D)

0.5–3.0 (unitless)

Ratio of tank height (H) to internal diameter (D); governs vertical velocity scale and Rayleigh number magnitude.

⚡ Engineering Impact:

Ratios >2.0 amplify buoyancy-driven mixing; ratios <0.8 impede axial thermal front propagation and increase wall losses.

Baffle Height Fraction (h_b/H)

0.15–0.40 (unitless)

Vertical extent of horizontal baffles expressed as fraction of total tank height.

⚡ Engineering Impact:

Fractions <0.20 permit unimpeded plume rise; >0.35 induce excessive pressure drop and localized stagnation, increasing local corrosion risk.

Rayleigh Number (Ra)

10⁷–10¹¹ (molten salt, 300–565 °C)

Dimensionless number quantifying the relative strength of buoyancy to viscous diffusion: Ra = g·β·ΔT·H³/(ν·α).

⚡ Engineering Impact:

Ra >10⁹ triggers turbulent natural convection — requiring baffles or geometry modification to suppress convective overturning.

Thermal Stratification Index (TSI)

0.70–0.95 (dimensionless)

Normalized measure of temperature gradient stability: TSI = (T_top − T_bottom)/ΔT_design.

⚡ Engineering Impact:

TSI <0.75 indicates severe mixing; sustained values <0.65 invalidate exergy-based efficiency claims per ASHRAE Guideline 36.

📐 Key Formulas

Rayleigh Number

Ra = g·β·ΔT·H³/(ν·α)

Predicts onset and intensity of natural convection in TES fluids

Variables:
Symbol Name Unit Description
Ra Rayleigh Number dimensionless Dimensionless number predicting onset and intensity of natural convection
g Gravitational acceleration m/s² Acceleration due to gravity
β Thermal expansion coefficient 1/K Volumetric thermal expansion coefficient of the fluid
ΔT Temperature difference K Characteristic temperature difference driving convection (e.g., between hot and cold boundaries)
H Characteristic length m Typical vertical dimension (e.g., height of fluid layer)
ν Kinematic viscosity m²/s Ratio of dynamic viscosity to fluid density
α Thermal diffusivity m²/s Ratio of thermal conductivity to product of density and specific heat capacity
Typical Ranges:
Molten salt (565 °C)
1×10⁹ – 5×10¹⁰
PCM capsule array (70 °C)
1×10⁶ – 2×10⁸
Water/glycol (90 °C)
1×10⁷ – 5×10⁸
⚠️ Design target: Ra < 1×10⁹ if baffles omitted; Ra < 1×10¹¹ acceptable with optimized baffles

Usable Energy Fraction (UEF)

UEF = ∫₀^H ρ·c_p·(T(z)−T_cold) dz / [ρ·c_p·ΔT_design·H]

Fraction of theoretical storage capacity deliverable within specified temperature band

Variables:
Symbol Name Unit Description
UEF Usable Energy Fraction Fraction of theoretical storage capacity deliverable within specified temperature band
ρ Density kg/m³ Density of the storage medium
c_p Specific Heat Capacity J/(kg·K) Specific heat capacity of the storage medium
T(z) Temperature Profile K Temperature as a function of depth z
T_cold Cold Temperature K Reference cold temperature of the storage system
ΔT_design Design Temperature Difference K Temperature difference between hot and cold extremes in design
H Storage Height m Height or depth of the thermal energy storage medium
Typical Ranges:
Well-baffled molten salt
0.82 – 0.91
Unbaffled cylindrical tank
0.55 – 0.73
PCM with conductive enhancement
0.68 – 0.85
⚠️ Minimum acceptable UEF = 0.80 per DOE SunShot TES validation protocol

🏭 Engineering Example

Crescent Dunes Solar Energy Project (NV, USA)

N/A — molten salt TES system
H/D
2.3
h_b/H
0.22, 0.48, 0.74
Ra_max
3.2×10¹⁰
UEF_measured
0.79
baffle_count
3
η_ex_design
0.71

🏗️ Applications

  • Concentrated Solar Power (CSP) tower plants
  • Industrial waste heat recovery with PCM
  • District heating seasonal storage

📋 Real Project Case

Concentrated Solar Power (CSP) Integration with Cement Kiln Preheater

Heidelberg Materials plant, Morocco

Challenge: Intermittent solar input mismatched with continuous kiln heat demand (350–450°C)
CSP Integration with Cement Kiln Preheater CSP Field Hot Salt Tank Thot ≈ 565°C Cold Salt Tank Tcold ≈ 290°C Thermocline Buffer Ceramic Aggregate Kiln Preheater 350–450°C Stratification Index: 0.82 Exergy Reduction: −37% Storage Duration: 12 h CSP / Kiln Hot Salt Cold Salt Thermocline
Read full case study →

🎨 Technical Diagrams

Hot inletCold outletH/D = 2.3 → high Ra
h_b/H = 0.22h_b/H = 0.48h_b/H = 0.74

📚 References

[1]
ASHRAE Guideline 36-2021: High-Performance Sequencing Control for HVAC Systems — American Society of Heating, Refrigerating and Air-Conditioning Engineers
[3]
Thermal Energy Storage Handbook — U.S. Department of Energy (DOE)