Sensible TES Sizing Using Specific Heat Integration and Temperature Gradient Constraints

Typical Scale

20–120 MWh thermal for industrial CSP integration

Key Standard

ASTM E1269 (cₚ measurement), ASME BPVC Section VIII Div 2

Industry Adoption

Used in >85% of commercial CSP plants with TES (e.g., Noor III, Ivanpah, Gemasolar)

⚠️ Why It Matters

1

Excessive ΔT across storage media

↓

2

Non-uniform thermal stress in containment vessels

↓

3

Cracking of refractory linings or salt tanks

↓

4

Unplanned downtime due to thermal fatigue failure

↓

5

Loss of process heat availability during critical production windows

↓

6

Violation of ASME BPVC Section VIII Div 2 fatigue life requirements

📘 Definition

Sensible TES sizing is the thermodynamic and process-engineering procedure that determines the minimum required mass and volume of a temperature-dependent storage medium (e.g., nitrate salts, concrete, or thermal oil), constrained by the system’s allowable temperature gradient, charge/discharge power matching, and exergy preservation during cyclic operation. It integrates specific heat capacity as a function of temperature (cₚ(T)) over the operational ΔT range and enforces thermal-hydraulic and material stability limits to avoid solidification, decomposition, or excessive pressure drop.

🎨 Concept Diagram

AI-generated illustration for visual understanding

💡 Engineering Insight

Never assume cₚ is constant—even for ‘well-characterized’ salts like Solar Salt. Between 290°C and 565°C, cₚ varies by 37%; using a mean value (1.52 kJ/kg·K) instead of integrating cₚ(T) = 0.83 + 0.0012T − 1.4×10⁻⁶T² yields a 9.3% energy deficit at 500°C discharge. Always source cₚ(T) from NIST SRM 2911a calibration data or vendor-supplied DSC reports traceable to ISO/IEC 17025 labs.

📖 Detailed Explanation

Sensible TES stores energy by raising the temperature of a material—no phase change occurs. The simplest energy balance says Q = m·cₚ·ΔT, but this assumes cₚ is constant. In reality, cₚ changes significantly with temperature, especially near melting points or decomposition thresholds. For molten salts used above 400°C, ignoring this variation leads to systematic undersizing.

Advanced sizing uses the integral form: Q = m·∫ᵀᶜᵀʰ cₚ(T) dT. This requires experimentally derived cₚ(T) functions—not textbook averages. Industry practice now mandates NIST-traceable DSC measurements per ASTM E1269, with uncertainty budgets reported to ±0.015 kJ/kg·K. The integration must also account for thermal losses, pump work, and exchanger effectiveness—each contributing multiplicative correction factors to the base mass calculation.

At the system level, temperature gradient constraints dominate mechanical design. A 50 K ΔT across a 12 m tall salt tank induces ~11 MPa thermal stress in 316L stainless steel—exceeding ASME’s fatigue-allowable stress range for 10⁵ cycles. Hence, sizing isn’t just about energy—it’s a coupled thermal-mechanical-fluid problem. Modern tools like COMSOL Multiphysics v6.2 embed cₚ(T) lookups, non-linear conduction, and creep models to resolve these interactions before fabrication begins.

🔄 Engineering Workflow

Step 1

Step 1: Extract time-resolved process heat demand profile (MW vs. hour) from DCS historian data

→

Step 2

Step 2: Define operational bounds: T_min/T_max, allowable ΔT across storage, and max charge/discharge power (MW)

→

Step 3

Step 3: Obtain cₚ(T) experimental curve (DSC/TGA) or validated NIST-certified polynomial coefficients for storage medium

→

Step 4

Step 4: Integrate ∫m·cₚ(T)·dT over ΔT to compute minimum mass; apply η_strat and safety factor (1.12–1.25) to determine final volume

→

Step 5

Step 5: Perform transient thermal-structural FEA (ANSYS Mechanical) on tank geometry using real cₚ(T) and convection BCs

→

Step 6

Step 6: Validate exergy efficiency (>78% round-trip) via second-law analysis using pinch point temperature and entropy generation rates

→

Step 7

Step 7: Commission with step-load thermal response test and IR thermography to verify ΔT distribution matches design

📋 Decision Guide

Rock/Field Condition	Recommended Design Action
Process demand requires >4 h discharge at 420°C with ±3°C tolerance	Use dual-tank stratified molten salt with cₚ(T)-integrated sizing; enforce ΔT ≤ 22 K across bed; specify ASTM A312 TP347H piping
Site has limited footprint (<1,200 m²) and must serve intermittent 1.8 MW steam load	Adopt high-density concrete (ρ = 2,350 kg/m³, cₚ = 0.87 kJ/kg·K) with vertical finned tube heat exchangers; limit ΔT to 18 K; validate via ANSYS Fluent transient CFD
Retrofit into existing 30-year-old refinery with 12 bar saturated steam header	Size single-tank oil-based TES (Therminol VP-1) using cₚ(T) spline fit from DSC data; cap max T_hot at 390°C; install inline particle filter to prevent coking-induced fouling

📊 Key Properties & Parameters

cₚ(T) integration error

±2.1% to ±8.7% for 220–565°C NaNO₃–KNO₃ (60:40)

Relative deviation in stored energy calculation when cₚ is assumed constant versus integrated over measured cₚ(T) curve

⚡ Engineering Impact:

Directly causes 3–12% under- or over-sizing of tank volume, risking thermal runaway or insufficient discharge duration

Maximum allowable ΔT across storage bed

15–45 K for stainless steel 316L tanks; 8–22 K for ceramic-packed beds

Largest permissible temperature difference between inlet and outlet fluid during full-rate discharge, governed by thermal stress limits in structural materials

⚡ Engineering Impact:

Dictates minimum mass flow rate and maximum bed height — violating this triggers ASME Code Case N-900 fatigue assessments

Thermal stratification efficiency (η_strat)

0.68–0.89 for well-designed packed-bed systems; 0.52–0.75 for single-tank forced-convection designs

Ratio of usable sensible energy in a stratified tank to total theoretical energy, quantifying mixing losses during charge/discharge

⚡ Engineering Impact:

Reduces effective storage capacity by up to 38%, forcing oversized tanks or auxiliary heating to meet process setpoints

Material decomposition onset temperature (T_decomp)

540–580°C for Solar Salt (60NaNO₃–40KNO₃); 420–450°C for Hitec XL

Temperature at which the storage medium begins irreversible chemical breakdown, measured via TGA under inert atmosphere

⚡ Engineering Impact:

Sets absolute upper bound on hot-side operating temperature — exceeding it degrades heat transfer fluid compatibility and corrodes piping

📐 Key Formulas

Energy Storage Capacity (Integrated)

Q = m \int_{T_c}^{T_h} c_p(T)\, dT

Total sensible energy stored, accounting for temperature-dependent specific heat

Variables:

Symbol	Name	Unit	Description
Q	Energy Storage Capacity	J	Total sensible energy stored
m	Mass	kg	Mass of the storage material
c_p(T)	Specific Heat Capacity	J/(kg·K)	Temperature-dependent specific heat capacity
T_c	Cold Temperature	K	Lower temperature limit of integration
T_h	Hot Temperature	K	Upper temperature limit of integration

Typical Ranges:

Molten salt (290–565°C)

0.92–1.68 kJ/kg·K

Concrete (100–400°C)

0.78–0.94 kJ/kg·K

⚠️ cₚ(T) uncertainty < ±0.015 kJ/kg·K per ASTM E1269

🔧 Open Calculator

Minimum Mass with Stratification Efficiency

m_{min} = \frac{Q_{demand}}{\eta_{strat} \int_{T_c}^{T_h} c_p(T)\, dT}

Corrected mass accounting for thermal mixing losses in stratified tanks

Variables:

Symbol	Name	Unit	Description
m_{min}	Minimum Mass	kg	Corrected mass accounting for thermal mixing losses in stratified tanks
Q_{demand}	Thermal Demand	J	Required thermal energy
\eta_{strat}	Stratification Efficiency	dimensionless	Efficiency factor accounting for thermal mixing losses in stratified thermal storage tanks
T_c	Cold Temperature	K	Lower temperature limit of the integration range
T_h	Hot Temperature	K	Upper temperature limit of the integration range
c_p(T)	Specific Heat Capacity	J/(kg·K)	Temperature-dependent specific heat capacity of the storage medium