17-point physics cross-check for cryogenic heat exchanger design validation.
Correlations Used
Correlation
Application
Range
Gnielinski (1976)
Primary Nusselt number
2300 < Re < 5×10⁶, 0.5 < Pr < 2000
Dittus-Boelter (1930)
Cross-check Nusselt
Re > 10000, 0.6 < Pr < 160
Petukhov (1970)
Darcy friction factor
3000 < Re < 5×10⁶
Seban-McLaughlin
Curved tube enhancement
De > 11.6, spiral HX
Žukauskas (1987)
Shell-side cross-flow
Re 1–2×10⁵, spiral HX
Fluid Properties — NIST Validated
Helium (40–300 K): Ideal gas EOS (R=2077.1 J/kg·K), Cp=5193 J/kg·K. Viscosity and conductivity via ln-polynomial fits validated against CoolProp/REFPROP within 2–5%.
Para-Hydrogen (20–31 K liquid, 31–300 K gas): CoolProp 7.2.0 ParaHydrogen HEOS. LH₂ is >99.8% para-H₂ at 20K. T_sat(P) = quadratic in ln(P), error <0.07K across 1–12 bar. h_fg(P) pressure-dependent: 446→237→143 kJ/kg from 1→10→12 bar. Gas cp includes ortho→para conversion peak (~16 kJ/kgK at 150K). 4-segment Simpson integration for enthalpy accuracy.
Validation Checks (17 points)
Reynolds numbers in turbulent range (Re > 4000) for both streams
Prandtl numbers physical (0.1 < Pr < 1000)
Nusselt numbers match Gnielinski ±5% vs Dittus-Boelter cross-check
Energy balance to machine precision (|Q_hot − Q_cold|/Q_avg < 1e−10)
Effectiveness bounded (0 < ε < 1)
Second law satisfied (no temperature crossover)
Subcooling margin ≥ 5 K (LH₂ stays liquid at inlet and outlet)