🎓 Lesson 10
D5
Transient Resistance Networks: Modeling Charge/Discharge Dynamics
Transient resistance networks model how electrical resistance changes over time when energy is stored or released—like how a battery charges or discharges in a thermal energy storage system.
🎯 Learning Objectives
- ✓ Calculate the time constant (τ) of a transient resistance network given thermal capacitance and effective resistance
- ✓ Design a charge/discharge profile that maintains resistance drift within ±5% of nominal value over 10,000 cycles
- ✓ Analyze voltage decay curves to identify dominant resistance mechanisms (e.g., contact degradation vs. bulk heating)
- ✓ Explain how thermal runaway thresholds are influenced by transient resistance growth in phase-change material (PCM)–integrated batteries
- ✓ Apply the Thévenin–RC hybrid model to size thermal management subsystems for grid-scale thermal energy storage
📖 Why This Matters
In industrial thermal energy storage (TES) systems—such as those using molten salt, graphite-PCM composites, or resistive-heating batteries—the electrical resistance of components doesn’t stay constant. As temperature rises during charging, resistance shifts, altering power delivery, efficiency, and safety margins. Ignoring this transient behavior leads to oversized heaters, undersized cooling, or catastrophic thermal runaway. Understanding transient resistance networks ensures accurate dynamic load matching—critical for ROI-driven TES deployments in cement kilns, steel reheating furnaces, and solar-thermal peaking plants.
📘 Core Principles
Transient resistance arises from three coupled phenomena: (1) intrinsic temperature dependence of resistivity (ρ = ρ₀[1 + α(T − T₀)]), (2) interfacial resistance evolution at electrode–PCM boundaries due to microcracking or interdiffusion, and (3) geometric deformation under thermal cycling that alters current path length and cross-section. These are modeled as a time-varying resistor R(t) embedded in an RC ladder network, where C represents thermal capacitance (J/K) and the parallel RC branch captures thermal inertia. At high frequencies (>1 Hz), skin effect and eddy losses may require adding inductive elements—but for industrial TES duty cycles (minutes-to-hours), first-order RC dominance holds per IEEE 1547-2018 Annex D guidelines.
📐 Thermal Time Constant & Effective Resistance
The thermal time constant τ quantifies how quickly resistance responds to thermal excitation. It governs system stability and informs heater/coolant sizing. When resistance varies linearly with temperature, τ defines the delay between power input and steady-state resistance value.
Thermal Time Constant
τ = R_th × C_thTime required for system resistance to reach ~63.2% of its final value after a step change in power input.
Variables:
| Symbol | Name | Unit | Description |
|---|---|---|---|
| τ | Thermal time constant | s | Characteristic response time of resistance to thermal excitation |
| R_th | Effective thermal resistance | K/W | Total conductive/convective resistance between heater and storage medium |
| C_th | Thermal capacitance | J/K | Heat capacity of active thermal mass (including interface layers) |
Typical Ranges:
Molten salt TES (industrial): 500 – 5000 s
Graphite-PCM battery module: 300 – 1200 s
Resistive ceramic heater array (cement precalciner): 100 – 400 s
💡 Worked Example
Problem: A graphite-PCM thermal battery has thermal capacitance C_th = 42 kJ/K and effective thermal resistance R_th = 0.018 K/W (measured via step-power test). Calculate τ and determine if the system can respond to a 15-minute load shift without exceeding 95% of final resistance value.
1.
Step 1: Identify knowns — C_th = 42,000 J/K, R_th = 0.018 K/W
2.
Step 2: Apply τ = R_th × C_th = 0.018 × 42,000 = 756 seconds (≈12.6 minutes)
3.
Step 3: Compute response at t = 15 min (900 s): R(t)/R_∞ = 1 − e^(−t/τ) = 1 − e^(−900/756) ≈ 1 − e^(−1.19) ≈ 1 − 0.304 = 0.696 → 69.6% of final value; but since resistance *increases* asymptotically, full stabilization requires ~4τ = 50.4 min. To reach 95%, t = −τ·ln(1−0.95) = −756·ln(0.05) ≈ 2270 s ≈ 37.8 min.
Answer:
The time constant is 756 s (12.6 min); the system reaches only ~69.6% of its final resistance after 15 minutes — insufficient for tight load-matching control without predictive compensation.
🏗️ Real-World Application
At the 24 MWth Cerro Dominador solar tower plant (Chile), transient resistance modeling revealed that molten salt heater elements exhibited 12% resistance increase over 200°C–565°C range. Without correcting for this in the SCADA-based load scheduler, peak current exceeded IEC 60502-2 limits during ramp-up, causing premature tripping. Engineers retrofitted a real-time R(t) estimator using thermocouple feedback and a lookup table derived from ASTM E2847-22 calibration data — improving cycle efficiency by 4.3% and extending heater life by 3×.
🔧 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