| name | stirling-engine |
| description | Stirling engine thermodynamics — Stirling cycle (isothermal expansion/compression, isochoric regeneration), ideal cycle efficiency (Carnot), Schmidt analysis (practical closed-form), regenerator effectiveness and NTU method, displacer/piston kinematics (alpha/beta/gamma configurations), dead volume penalty, power density, heat exchanger design (heater/cooler tubes), free-piston Stirling engines (FPSE), and applications in cryocoolers, solar dish Stirling, and CHP systems. |
| metadata | {"priority":7,"promptSignals":{"phrases":["Stirling engine","Stirling cycle","Stirling cooler","free piston Stirling","regenerator Stirling","Schmidt analysis"],"minScore":3}} |
Stirling Engine — Complete Skill
Stirling Cycle Thermodynamics
Ideal Stirling Cycle
Four processes (ideal):
- Isothermal expansion (T_H): heat addition; piston work done by gas at constant T_H
- Isochoric cooling (constant volume): gas transferred through regenerator; cools from T_H to T_C; heat stored in regenerator
- Isothermal compression (T_C): heat rejection; piston work done on gas at constant T_C
- Isochoric heating (constant volume): gas transferred back through regenerator; reheats from T_C to T_H; heat recovered from regenerator
Ideal Stirling efficiency:
With perfect regenerator: η_Stirling = η_Carnot = 1 − T_C/T_H [Kelvin temperatures]
This is remarkable: Stirling reaches Carnot efficiency with ideal regenerator
Without regenerator: η_no_regen = (1 − T_C/T_H) / (1 + c_v×(T_H−T_C)/(R×T_H×ln(r))) [r = compression ratio]
Specific work per cycle (ideal):
W_net = m_gas × R × (T_H − T_C) × ln(V_max/V_min) [m = mass of gas; R = specific gas constant]
Schmidt Analysis (Practical)
Sinusoidal Motion Analysis
Schmidt (1871) analysis: accounts for sinusoidal (not ideal square-wave) piston/displacer motions
Allows closed-form solution for pressure, work, and power — most widely used analytical model for Stirling engines
Alpha-configuration Stirling:
Two pistons (expansion and compression) in separate cylinders connected by regenerator + heater + cooler
Phase angle between pistons: 90° (optimal for max power; varies 60–120° in practice)
Schmidt pressure variation:
p(θ) = P_mean / (1 − B × cos(θ − θ_0))
Where: B = amplitude ratio (function of geometry); θ = crank angle; θ_0 = phase angle
Mean pressure:
P_mean = m_total × R × T_mean_harmonic / V_total [T_mean_harmonic = harmonic mean of temperatures weighted by volumes]
Power output (Schmidt):
W_cycle = π × P_mean × V_SW × sin(α) × (√((T_H/T_C)−1) / (√(T_H/T_C)+1)) × correction_factors
[approximate; V_SW = swept volume; α = phase angle; complete Schmidt analysis is algebraic]
Working gas selection:
| Gas | k [W/mK] | η [μPa·s] | cp [J/kgK] | Advantage |
|---|
| Hydrogen H₂ | 0.18 | 9.0 | 14,300 | Highest power density; safety concern |
| Helium He | 0.15 | 19.8 | 5,200 | Best compromise: high k, safe |
| Air | 0.026 | 18.5 | 1,005 | Low cost; low performance |
| Nitrogen | 0.026 | 17.8 | 1,040 | Low cost; slightly safer than air |
Hydrogen: 2.5× more power than helium at same conditions → used in Philips/Ford Stirling engines
Helium: preferred for most applications (safety, inertness; cryocoolers always He)
Regenerator Design
Regenerator Effectiveness
Regenerator function: stores heat when gas flows hot→cold; returns heat when cold→hot
Ideal: 100% heat recovery → heat input = W_net / η_Carnot only
Real: ε_regen = Q_stored / Q_max = 1 − (T_H − T_out_hot) / (T_H − T_C)
NTU (Number of Transfer Units) method:
NTU = UA / (ṁ × cp) [U = overall heat transfer coefficient; A = regenerator matrix surface area]
For counterflow matrix: ε = NTU / (NTU + 1) [simplified balance; counterflow heat exchanger analogy]
Target: ε > 0.95; NTU > 19 for ε = 0.95
Regenerator matrix materials:
Wire mesh screens (SS 316, Fe, stainless): high surface area; hydraulic diameter D_h = d_wire (if packed)
Foil: annular foil wound with spacers; very low pressure drop; lower surface area than mesh
Felt metal: sintered fiber matrix; uniform pore size; good regenerator performance
Random fiber: porosity 0.7–0.9; specific surface area 5000–20,000 m²/m³
Regenerator design tradeoff:
More matrix → better ε, higher NTU; but also higher dead volume (reduces power) and pressure drop
Dead volume penalty: V_dead reduces effective swept volume ratio → reduces work per cycle
Optimal matrix length: maximize (P × ε) / (1 + V_dead/V_swept)
Dead Volume Effect
Dead volume: regenerator volume + heat exchanger volumes (heater tubes + cooler tubes) + connecting passages
Effect: gas in dead volume does not contribute to work; effectively dilutes working gas
Power reduction factor: approx (V_swept / (V_swept + 2×V_dead_equiv))^n [empirical; n ≈ 0.5–1.0]
Design target: V_dead ≤ 0.5 × V_swept
Engine Configurations
Alpha, Beta, Gamma
Alpha configuration:
Two pistons, two cylinders; one hot, one cold; regenerator between
Higher friction (two seals); good power density; NASA FPSE type
Beta configuration:
Single cylinder; piston + displacer (coaxial); displacer moves gas between hot/cold spaces
Most compact; lower dead volume; classic RE-1000 and GPU-3 research engines
Displacer: Teflon-coated; low friction; minimal seal
Gamma configuration:
Separate cylinder for displacer and piston; connected by passage
Easy to manufacture; good power control; domestic CHP Stirling (WhisperGen, Sunpower)
Free-Piston Stirling Engine (FPSE)
Key feature: no crankshaft; piston and displacer oscillate freely at resonant frequency
No rubbing seals (gas bearings); long life (>60,000 hours demonstrated); suitable for space
Power extraction: linear alternator; piston motion directly converts to electricity
FPSE dynamics:
Spring-mass system: gas + mechanical springs set resonance frequency
f_res = 1/(2π) × √(K_eff / m_piston) [K_eff = combined gas spring + mechanical spring stiffness]
Phase between piston and displacer: determined by dynamics; controlled by spring selection
Sunpower 35W FPSE:
Power: 35 W electrical; η = 25% (Carnot = 33%); T_H = 650°C; T_C = 35°C; 20 year design life (no wear)
Used in: remote power, space (Advanced Stirling Radioisotope Generator ASRG — NASA cancelled 2013)
Applications
Solar Dish-Stirling
Dish concentrator + Stirling engine at focal point:
Concentration ratio: 1000–3000× (parabolic dish 10–50 m²)
T_H = 650–850°C (receiver cavity); T_C = ambient 20–40°C
η_system = η_dish × η_receiver × η_engine ≈ 0.88 × 0.88 × 0.35 ≈ 27% (solar-to-electricity)
Highest solar-to-electricity conversion of any technology (historically demonstrated 29.4% at Sandia, 1984)
System: WG Stirling/Kockums V160 (4 cylinders, 25 kW); control: track sun + regulate H₂ pressure
CHP Stirling (Domestic)
Micro-CHP:
Input: natural gas burner (T_H ≈ 650°C); Output: 1 kW electric + 6–8 kW heat (for space heating)
η_electrical ≈ 12–15% (low priority); η_total ≈ 90–95% (heat recovery primary benefit)
Units: Remeha, Baxi, Honda (Japan); very low noise vs. ICE-based micro-CHP
Stirling Cryocoolers
Reversed Stirling = cooler:
Input: mechanical work; output: heat pumped from cold reservoir (cryo-cooling)
Applications: IR detector cooling (FLIR cameras); cryogenic infrared sensors; LNG reliquefaction
Pulse tube refrigerator (PTR): simplified Stirling cooler; no moving part at cold end; very reliable
T_cold achievable: 20–80 K (single-stage Stirling); 2–4 K (multi-stage / JT aftercooler)
Standards and References
| Standard | Scope |
|---|
| ASME Boiler Code (Section VIII) | Pressure vessel for Stirling engine receiver |
| IEC 62282-3-100 | Fuel cell technologies (micro-CHP Stirling adjacent) |
| SAE J1297 | Alternative fuel vehicle fueling systems (Stirling context) |
| ISO 15500 | CNG vehicle components (related to H₂ Stirling) |
| Walker "Stirling Engines" (1980) | Classic reference textbook |
| Urieli & Berchowitz "Stirling Cycle Engine Analysis" | Standard analytical reference |
Output
Provide: application (power generation/cooler/CHP/space; output power [W]; T_H [°C]; T_C [°C]), ideal efficiency (η_Carnot = 1−T_C/T_H [%]; realistic η_engine ≈ 0.6–0.75×η_Carnot; realistic η [%]), working gas (He/H₂/air; charge pressure [MPa]; viscosity/conductivity trade-off), configuration (alpha/beta/gamma/FPSE; crankshaft or free-piston; phase angle [°]), swept volume (V_swept [cc]; P_mean [MPa]; W_cycle from Schmidt [J]; Power = W_cycle × f [W]; frequency f [Hz]), dead volume (V_dead [cc]; V_dead/V_swept ratio ≤ 0.5 target; power penalty factor), regenerator (material: wire mesh/foil; mesh size; D_h [mm]; length [mm]; ε ≥ 0.95; NTU; pressure drop [kPa]), heat exchangers (heater tubes: T_H [°C]; number × diameter × length; material: Inconel 625/SS310; Nusselt correlation; ΔT between gas and wall; cooler: water-cooled; Q_cooler [W]), engine mass and dimensions (L [mm] × D [mm]; mass [kg]; power density [W/kg]), and applicable standard (ASME Section VIII for receiver if pressurized; Walker/Urieli reference for verification).
