Thermodynamics & Heat Transfer

Formula sheet · learning zone · practice problems with answer key

Engineering thermodynamics, steam plant & heat transfer · 70 formulas · 68 practice problems · metric edition 1

The practice problems are edition-drawn: a later edition deals fresh numbers, so keep the key with the printing it came from. The Exam Room deals new numbers on every attempt.

The formula sheet

Gauge and Absolute Pressure
Pabs=Pgauge+PatmP_{abs} = P_{gauge} + P_{atm}
Sensible Heat (Q = mcΔT)
Q=mcΔTQ = m c \Delta T
Latent Heat
Q=mLQ = m L
Boyle's Law
P1V1=P2V2P_1 V_1 = P_2 V_2
Charles's Law
V1T1=V2T2\frac{V_1}{T_1} = \frac{V_2}{T_2}
Gay-Lussac's Law
P1T1=P2T2\frac{P_1}{T_1} = \frac{P_2}{T_2}
Combined Gas Law
P1V1T1=P2V2T2\frac{P_1 V_1}{T_1} = \frac{P_2 V_2}{T_2}
Ideal Gas Law
PV=nRTP V = n R T
Gas Density from Molar Mass
ρ=PMRT\rho = \frac{PM}{RT}
Compressibility Factor (Z = PV/nRT)
Z=PVnRTZ = \frac{P V}{n R T}
Van der Waals Equation of State
(P+an2V2)(Vnb)=nRT\left(P + \frac{a n^{2}}{V^{2}}\right)\left(V - n b\right) = n R T
Antoine Equation (Vapour Pressure)
log10P=ABC+T\log_{10} P = A - \frac{B}{C + T}
Clausius–Clapeyron Equation (Two-Point Form)
ln ⁣(P2P1)=ΔHvapR(1T21T1)\ln\!\left(\frac{P_2}{P_1}\right) = -\frac{\Delta H_{vap}}{R}\left(\frac{1}{T_2} - \frac{1}{T_1}\right)
Saturation Temperature and Pressure of Steam
Tsat=Ts(psat)psat=ps(Tsat)T_{sat} = T_s(p_{sat}) \qquad p_{sat} = p_s(T_{sat})
Steam Quality from Enthalpy
x=hhfhfgx = \frac{h - h_f}{h_{fg}}
Wet Steam Enthalpy from h_f and h_g
h=(1x)hf+xhgh = (1 - x) \, h_f + x \, h_g
Degrees of Superheat
ΔTsh=TTs(p)\Delta T_{sh} = T - T_s(p)
Flash Steam Percentage
%F=hf1hf2hfg2×100\%F = \frac{h_{f1} - h_{f2}}{h_{fg2}} \times 100
Flash Steam Mass Rate
m˙f=%F100m˙c\dot{m}_f = \frac{\%F}{100} \, \dot{m}_c
Steam Turbine Specific Work
w=h1h2w = h_1 - h_2
Steam Turbine Power Output
P=m˙wP = \dot{m} \, w
Turbine Isentropic Efficiency
ηisen=h1h2h1h2s\eta_{isen} = \frac{h_1 - h_2}{h_1 - h_{2s}}
Napier's Steam Leak Rate
m˙=AP70\dot{m} = \frac{A \, P}{70}
Desuperheater Water Injection Rate
m˙w=m˙1h1h2h2hw\dot{m}_w = \dot{m}_1 \, \frac{h_1 - h_2}{h_2 - h_w}
Steam Coil Condensate Load
m˙=Q˙hfg\dot m = \frac{\dot Q}{h_{fg}}
Thermal Efficiency
η=WQh\eta = \frac{W}{Q_h}
Carnot Efficiency
η=1TcTh\eta = 1 - \frac{T_c}{T_h}
Rankine Cycle Thermal Efficiency
η=(h1h2)(h4h3)h1h4\eta = \frac{\left(h_1 - h_2\right) - \left(h_4 - h_3\right)}{h_1 - h_4}
Otto Cycle Efficiency (Compression Ratio)
η=11rγ1\eta = 1 - \frac{1}{r^{\gamma - 1}}
Brayton Cycle Efficiency (Pressure Ratio)
η=11rp(γ1)/γ\eta = 1 - \frac{1}{r_p^{\left(\gamma - 1\right)/\gamma}}
Heat Conduction Rate
P=kAΔTdP = \tfrac{k A \Delta T}{d}
R-Value of an Insulation Layer (R = L/k)
R=LkR = \frac{L}{k}
Thermal Resistance of a Plane Wall
R=LkAR = \frac{L}{k A}
Thermal Resistances in Series
Rtot=R1+R2+R3R_{tot} = R_1 + R_2 + R_3
Total R-Value of an Assembly
Rtot=R1+R2+R3R_{tot} = R_1 + R_2 + R_3
U-Factor from Total R-Value (U = 1/R)
U=1RtotU = \frac{1}{R_{tot}}
Overall U from Total Resistance
U=1RtotAU = \frac{1}{R_{tot} A}
Heat Flow from Thermal Resistance
Q˙=ΔTR\dot{Q} = \frac{\Delta T}{R}
Heat Loss Through an Assembly (Q = A·ΔT/R)
Q˙=AΔTRtot\dot{Q} = \frac{A \, \Delta T}{R_{tot}}
Heat Flux Through Insulation (q = ΔT/R)
q=ΔTRq'' = \frac{\Delta T}{R}
Effective R-Value with Framing (Parallel Path)
1Reff=ffrRfr+1ffrRcav\frac{1}{R_{eff}} = \frac{f_{fr}}{R_{fr}} + \frac{1 - f_{fr}}{R_{cav}}
Conduction Through a Pipe Wall
Q˙=2πkLΔTln(r2/r1)\dot{Q} = \frac{2 \pi k L \, \Delta T}{\ln(r_2 / r_1)}
Critical Radius of Insulation
rcr=khr_{cr} = \frac{k}{h}
Newton's Law of Cooling (Q = hAΔT)
Q˙=hAΔT\dot{Q} = h A \, \Delta T
Convection Film Resistance
R=1hAR = \frac{1}{h A}
Reynolds Number
Re=ρvDμRe = \frac{\rho v D}{\mu}
Prandtl Number
Pr=μcpk\mathrm{Pr} = \frac{\mu c_p}{k}
Nusselt Number
Nu=hLk\mathrm{Nu} = \frac{h L}{k}
Grashof Number
Gr=gβΔTL3ν2\mathrm{Gr} = \frac{g \, \beta \, \Delta T \, L^{3}}{\nu^{2}}
Rayleigh Number
Ra=GrPr\mathrm{Ra} = \mathrm{Gr} \, \mathrm{Pr}
Dittus-Boelter Correlation
Nu=0.023Re0.8Prn\mathrm{Nu} = 0.023 \, \mathrm{Re}^{0.8} \, \mathrm{Pr}^{n}
Stefan-Boltzmann Law
P=εσAT4P = \varepsilon \sigma A T^4
Net Radiation Exchange Between Surfaces
Q˙=εσA(T14T24)\dot{Q} = \varepsilon \sigma A (T_1^4 - T_2^4)
Combined Convection and Radiation Coefficient
ht=hc+εσ(Ts+Tsur)(Ts2+Tsur2)h_t = h_c + \varepsilon \sigma (T_s + T_{sur})(T_s^2 + T_{sur}^2)
Overall Heat Transfer Coefficient (U)
1U=1hi+Lk+1ho\frac{1}{U} = \frac{1}{h_i} + \frac{L}{k} + \frac{1}{h_o}
Fouled Overall Coefficient
1Uf=1Uc+1hf\frac{1}{U_f} = \frac{1}{U_c} + \frac{1}{h_f}
Fouling Factor on an Overall Coefficient
1Uf=1Uc+Rf\frac{1}{U_f} = \frac{1}{U_c} + R_f
Log Mean Temperature Difference (Counterflow)
ΔTlm=ΔT1ΔT2ln(ΔT1/ΔT2)\Delta T_{lm} = \frac{\Delta T_1 - \Delta T_2}{\ln(\Delta T_1 / \Delta T_2)}
Log Mean Temperature Difference (Parallel Flow)
ΔTlm=ΔT1ΔT2ln(ΔT1/ΔT2)\Delta T_{lm} = \frac{\Delta T_1 - \Delta T_2}{\ln(\Delta T_1 / \Delta T_2)}
Heat Exchanger Duty (Q = U·A·F·LMTD)
Q˙=UAFΔTlm\dot{Q} = U A F \, \Delta T_{lm}
Stream Duty from Mass Flow (Q = ṁcΔT)
Q˙=m˙cpΔT\dot{Q} = \dot{m} \, c_p \, \Delta T
Number of Transfer Units (NTU)
NTU=UAm˙cp\mathrm{NTU} = \frac{U A}{\dot{m} \, c_p}
Capacity Rate Ratio (Cr)
Cr=m˙mincminm˙maxcmaxC_r = \frac{\dot{m}_{min} c_{min}}{\dot{m}_{max} c_{max}}
Maximum Possible Heat Transfer (Qmax)
Q˙max=m˙mincmin(Th,inTc,in)\dot{Q}_{max} = \dot{m}_{min} c_{min} (T_{h,in} - T_{c,in})
Heat Exchanger Effectiveness (ε = Q/Qmax)
ε=Q˙Q˙max\varepsilon = \frac{\dot{Q}}{\dot{Q}_{max}}
Effectiveness from NTU (Counterflow)
ε=1eNTU(1Cr)1CreNTU(1Cr)\varepsilon = \frac{1 - e^{-\mathrm{NTU}(1 - C_r)}}{1 - C_r \, e^{-\mathrm{NTU}(1 - C_r)}}
Biot Number
Bi=hLck\mathrm{Bi} = \frac{h L_c}{k}
Fourier Number
Fo=ktρcL2\mathrm{Fo} = \frac{k \, t}{\rho \, c \, L^{2}}
Lumped Capacitance Time Constant
τ=ρVchA\tau = \frac{\rho V c}{h A}
Lumped Capacitance Cooling Curve
T=T+(T0T)et/τT = T_\infty + (T_0 - T_\infty) e^{-t/\tau}

Properties of State

Gauge and Absolute Pressure

Pabs=Pgauge+PatmP_{abs} = P_{gauge} + P_{atm}
PabsPgaugePatm
Where
  • PabsP_{abs}= Absolute pressure (kPa)
  • PgaugeP_{gauge}= Gauge pressure (kPa)
  • PatmP_{atm}= Atmospheric pressure (kPa)

A tire gauge reads zero in open air, yet a barometer standing beside it says 101 kPa. Both instruments are right, because they use different zeros. A gauge has atmosphere on the back of its diaphragm, so it can only ever report the difference between what it is connected to and the air around it. Absolute pressure is measured from a sealed vacuum reference and counts everything. Pabs=Pgauge+PatmP_{abs} = P_{gauge} + P_{atm} is the bridge, and neither reading is more correct than the other — they answer different questions.

A tire inflated to 220 kPa gauge holds 321 kPa absolute. The relation matters just as much in the other direction: a "vacuum" rated at −80 kPa gauge is 21 kPa absolute, still holding a fifth of an atmosphere. And there is a floor. The most negative a gauge can ever read is −101.3 kPa at sea level, because a perfect vacuum on one side and one atmosphere on the other is the entire range available. Anyone quoting a −150 kPa vacuum is quoting something that does not exist.

The distinction dates to Torricelli's barometer in 1643, which settled two arguments at once: that the atmosphere has weight, and that a vacuum can exist. Before that there was no absolute zero of pressure to measure from. The unit conventions still carry the split — psia and psig, bara and barg, and in North American vacuum work a third convention entirely, inches of mercury below atmosphere, so that "28 inHg of vacuum" means a gauge reading of −94.8 kPa and an absolute pressure of 6.5 kPa.

Every gas law demands absolute pressure, and forgetting that is the classic error on this site. Boyle's law, Gay-Lussac's law, the combined gas law and PV=nRTPV = nRT all count molecular impacts, and a gauge has quietly subtracted an atmosphere from the count. Compressing a tire from 220 to 440 kPa gauge is not doubling the pressure — in absolute terms it goes from 321 to 541, a factor of 1.69, and a calculation done on the gauge figures is wrong by nearly 20%. The rule is simple enough to make automatic: convert to absolute before any gas law, and convert back only at the end if a gauge reading is what somebody wants.

Then the subtler trap, which is PatmP_{atm} itself. It is not 101.325 kPa where you are. That figure is the standard atmosphere at sea level, and elevation takes it away quickly — Calgary at 1045 m sits near 89 kPa, Denver near 84, and weather moves any of them by ±3 kPa on its own. Converting a gauge reading with a textbook 101.325 in a mountain city introduces a 12 kPa error, which is 12% of an atmosphere. Worse, the barometric pressure a weather service reports is usually sea-level corrected: it has been adjusted upward to what the pressure would be if the station were at sea level, specifically so that maps are comparable. It is not the local absolute pressure and must not be used as one. If you need the real figure, read a station barometer or compute it from elevation. This is not academic — it is why a pump's available NPSH is stated absolute and why a suction lift that works at sea level can cavitate at altitude, with the equation and the elevation together explaining exactly how much margin was lost.

Worked example: Tire 220 kPa gauge + 101.325 kPa atm → 321.325 kPa abs

Sensible Heat (Q = mcΔT)

Q=mcΔTQ = m c \Delta T
mcpQΔT
Where
  • QQ= Heat energy (J)
  • mm= Mass (kg)
  • cpc_p= Specific heat capacity (J/(kg·K))
  • ΔT\Delta T= Temperature change ()

Sensible heat is the energy that changes a substance's temperature without changing its phase. The specific heat capacity c is the price of each degree: how many joules one kilogram demands per kelvin of warming. Water's is famously steep at about 4186 J/(kg·K), which is why oceans moderate coastal climates and why a kettle takes its time. Heating 1.5 kg of water from 15 °C to 95 °C costs Q = 1.5 × 4186 × 80 ≈ 502 kJ.

The concept dates to Joseph Black's calorimetry experiments in 1760s Glasgow, which first pried apart the ideas of temperature and heat. Note that ΔT is a temperature difference, so a change of 80 °C equals a change of 80 K exactly — Fahrenheit differences convert by scale alone, with no offset. The formula holds as long as c stays roughly constant over the range and nothing melts or boils along the way.

Worked example: 2 kg water, c = 4186, ΔT = 30 C° → 251160 J

Latent Heat

Q=mLQ = m L
QLmm
Where
  • QQ= Heat absorbed or released (J)
  • mm= Mass changing phase (kg)
  • LL= Specific latent heat (J/kg)

Latent heat is the energy a phase change absorbs or releases while the temperature holds still. Melting 1 kg of ice at 0 °C soaks up 334 kJ — enough to heat that same water from 0 °C to 80 °C — yet the thermometer never moves until the last crystal is gone. Boiling is costlier still: vaporizing a kilogram of water takes about 2256 kJ, more than five times the energy needed to warm it from ice-cold to boiling.

Joseph Black coined the term in 1762 — latent means hidden, because the heat disappears into the phase change instead of the temperature reading. The physics runs everyday life: sweat cools you as it evaporates, steam scalds far worse than boiling water because it dumps its latent heat on condensing against skin, and every refrigerator moves heat by evaporating and condensing a working fluid in an endless loop.

One quantity, three notations, depending on whose book you are holding. Physics writes it L, as here, and splits it into Lf for fusion and Lv for vaporization. Engineering and every steam table print it hfg, where f is saturated fluid and g is saturated gas — so hfg = hg − hf is the gap between the two columns, which is exactly why it shrinks to nothing at the critical point where those columns meet. Chemistry writes it as an enthalpy of vaporization per MOLE rather than per kilogram, so its numbers look nothing like these until you divide by the molar mass. Same energy, three addresses.

Worked example: Melting 1 kg ice at 334 kJ/kg → 334000 J

Boyle's Law

P1V1=P2V2P_1 V_1 = P_2 V_2
P1V1P2V2
Where
  • P1P_1= Initial pressure (kPa)
  • V1V_1= Initial volume (L)
  • P2P_2= Final pressure (kPa)
  • V2V_2= Final volume (L)

Robert Boyle published his gas law in 1662, making it one of the oldest quantitative laws in physics: at constant temperature, the pressure and volume of a trapped gas are inversely proportional, so their product never changes. Cap a syringe and squeeze — compress 60 mL of air at 100 kPa down to 20 mL and the pressure climbs to 300 kPa. Molecularly, shrinking the space raises how often molecules hammer the walls, and pressure rises in exact proportion.

The law matters wherever gas gets squeezed. Scuba divers learn it first: air breathed at depth expands as they ascend, which is why the cardinal rule is never to hold your breath on the way up. Boyle's law assumes the temperature and the amount of gas stay fixed; change either and you need Charles's law or the combined gas law instead.

Worked example: 2 L at 1 atm → 4 atm gives 0.5 L

Charles's Law

V1T1=V2T2\frac{V_1}{T_1} = \frac{V_2}{T_2}
V1T1V2T2
Where
  • V1V_1= Initial volume (L)
  • T1T_1= Initial absolute temperature (°C)
  • V2V_2= Final volume (L)
  • T2T_2= Final absolute temperature (°C)

Hold the pressure on a gas constant and its volume follows its absolute temperature in strict proportion: warm it by 10% and it swells by 10%. The reason is worth having rather than memorising. Temperature is a measure of how fast the molecules are moving; pressure is how hard their impacts push on each square metre of wall. If the gas is to keep pushing with the same force while its molecules move faster, the only thing it can do is spread out, so that each patch of wall is struck less often. Volume rises exactly as fast as temperature to keep that balance, and V1/T1=V2/T2V_1/T_1 = V_2/T_2 is that sentence in symbols.

A 2.5 L balloon leaves a 22 °C room and goes into a −18 °C freezer. In kelvin those are 295.15 and 255.15, so V2=2.5×255.15/295.15=2.16 LV_2 = 2.5 \times 255.15/295.15 = 2.16\ \text{L} — it loses about 14% of its volume and visibly puckers. Bring it back out and it recovers. Nothing left the balloon; the same molecules simply stopped needing as much room.

Jacques Charles found this with hydrogen balloons around 1787 and never published it. Joseph Louis Gay-Lussac did the careful work and published in 1802, and generously named the result after Charles. The interesting part is what came out of extending the straight line. Plot volume against Celsius temperature and you get a line that, extrapolated backwards, hits zero volume at about −273 °C. No one in 1802 could reach anywhere near that temperature, yet the graph pointed straight at it. That extrapolation is how absolute zero was first located, and it is why William Thomson could propose an absolute scale in 1848 by simply moving the origin to where the gases were pointing. The modern value, −273.15 °C, is the zero of the kelvin scale.

Which is exactly why a ratio in Celsius is meaningless here. Going from 20 °C to 40 °C does not double the volume; in kelvin that is 293.15 to 313.15, a rise of under 7%. Worse, a Celsius ratio breaks outright when the temperature crosses zero — 0 °C in the denominator gives infinity, and a negative Celsius temperature gives a negative volume. This page converts your °C or °F entries to kelvin before it does anything, but the trap is worth recognising when you meet the equation off-screen.

Two smaller conditions do real work. The pressure must actually be constant: a gas sealed in a rigid tank obeys Gay-Lussac's law instead, where pressure rises and volume does not move at all. And the gas must stay a gas. Cool steam through 100 °C and the relation does not merely become inaccurate, it stops applying — the vapour condenses and the volume collapses by a factor of about 1600. A balloon is also only approximately constant-pressure, since the stretched rubber adds a little tension of its own, which is why the real shrinkage runs slightly under what the arithmetic predicts.

Worked example: 2 L at 300 K heated to 600 K → 4 L

Gay-Lussac's Law

P1T1=P2T2\frac{P_1}{T_1} = \frac{P_2}{T_2}
P1T1P2T2
Where
  • P1P_1= Initial pressure (kPa)
  • T1T_1= Initial absolute temperature (°C)
  • P2P_2= Final pressure (kPa)
  • T2T_2= Final absolute temperature (°C)

Seal a gas in a container that cannot change size and there is only one thing left for it to do when you heat it: push harder. Pressure then tracks absolute temperature in direct proportion, P1/T1=P2/T2P_1/T_1 = P_2/T_2. The mechanism is the same one behind Charles's law with the outcome swapped. Faster molecules strike the walls both more often and with more momentum each time, and since the walls will not move aside, all of that arrives as pressure.

Here is the case everyone actually meets, worked carefully, because the careless version is the standard error. A tire is set to 220 kPa gauge on a 5 °C morning and warms to 45 °C after an hour on the highway. Convert to absolute pressure first: 220 + 101 = 321 kPa absolute. Then P2=321×318.15/278.15=367 kPaP_2 = 321 \times 318.15/278.15 = 367\ \text{kPa} absolute, which is 266 kPa on the gauge — a rise of 46 kPa, about 6.6 psi. Run the same calculation on the gauge reading alone and you get 252 kPa, understating the rise by a third. This is why tire pressures are specified cold, and why topping up a hot tire leaves it soft in the morning.

The relation is usually credited to Gay-Lussac's 1802 paper, though Guillaume Amontons had it a century earlier: around 1702 he built an air thermometer that worked on precisely this principle, and noticed that the pressure line extrapolated toward a temperature below which it could not go. Some texts call it Amontons's law for that reason. Combine it with Boyle's law and Charles's law and you have the combined gas law; add Avogadro and you have PV=nRTPV = nRT, of which this is the constant-volume slice.

Absolute pressure is the trap here, more than absolute temperature. The ratio P1/P2P_1/P_2 is only meaningful when both pressures are measured from vacuum, because the equation is counting molecular impacts and a gauge has quietly subtracted an atmosphere from the count. Temperature has the same requirement for the same reason — kelvin, not Celsius — and this page converts your entries on both fronts. But when you meet the equation on paper, ask twice whether the pressure in your hand is gauge or absolute. It usually is gauge; almost every instrument in a mechanical room reads that way.

The other honest limit is that constant volume is an idealisation. A tire is not rigid — it grows a little as it warms, which relieves some of the pressure rise, so the measured increase runs slightly below the calculation. A steel cylinder is much closer to the ideal, which is what makes this law genuinely dangerous rather than merely academic. An aerosol can left on a dashboard, a propane cylinder in a closed vehicle, or a sealed pressure vessel in a fire all follow this line with nothing to relieve them, and the pressure keeps climbing until something gives. Relief valves exist because the equation has no upper bound.

Worked example: 3 atm at 300 K heated to 400 K → 4 atm (405.3 kPa)

Combined Gas Law

P1V1T1=P2V2T2\frac{P_1 V_1}{T_1} = \frac{P_2 V_2}{T_2}
P1V1T1P2T2V2
Where
  • P1P_1= Initial pressure (kPa)
  • V1V_1= Initial volume (L)
  • T1T_1= Initial absolute temperature (°C)
  • P2P_2= Final pressure (kPa)
  • V2V_2= Final volume (L)
  • T2T_2= Final absolute temperature (°C)

The combined gas law merges Boyle's, Charles's, and Gay-Lussac's laws into a single statement: for a fixed amount of gas, PV/T is constant. A weather balloon shows all three variables moving at once. Launched with 2.0 m³ of helium at 101 kPa and 288 K, it rises to where the pressure is 30 kPa and the temperature 228 K; its new volume is V₂ = P₁V₁T₂ ÷ (P₂T₁) = 2.0 × 101 × 228 ÷ (30 × 288) ≈ 5.3 m³ — more than double.

Hold any one variable constant and the named laws drop out: fix T for Boyle's law, fix P for Charles's, fix V for Gay-Lussac's. Temperatures must be absolute — the ratio of 20 °C to 40 °C is not 1:2 but 293:313 — and Celsius or Fahrenheit inputs convert to kelvin automatically. Add Avogadro's insight about the amount of gas n and this law becomes the full ideal gas law, PV = nRT.

Worked example: 1 L at 1 atm, 273.15 K → 0.5 atm, 546.3 K gives 4 L

Ideal Gas Law

PV=nRTP V = n R T
PVTn
Where
  • PP= Pressure (kPa)
  • VV= Volume (L)
  • nn= Amount (mol)
  • TT= Temperature (°C)

PV=nRTPV = nRT says that for a gas, four quantities are not independent: fix any three and the fourth is decided. Squeeze it and the pressure rises; warm it and it pushes harder or swells; add more of it and both go up. What makes the equation remarkable is not that those things are true — anyone with a bicycle pump knows them — but that one constant serves every gas. Helium, nitrogen and steam all obey it with the same R=8.314R = 8.314 J/(mol·K), which is a strong hint that pressure has nothing to do with what the molecules are and everything to do with how many there are and how fast they are moving.

A worked case in units you would actually read off a gauge. A 20 L cylinder sits at 150 kPa absolute on a 20 °C morning. Rearranged for amount, n=PV/RT=(150000×0.020)/(8.314×293.15)=3000/24371.23 moln = PV/RT = (150\,000 \times 0.020)/(8.314 \times 293.15) = 3000/2437 \approx 1.23\ \text{mol}. Note what had to happen before the arithmetic: pascals not kilopascals, cubic metres not litres, and kelvin not Celsius. This page converts your entries for you, but the discipline is worth keeping in your head, because a scrap of paper will not.

The law arrived in pieces. Boyle established PVPV constant at fixed temperature in 1662; Charles and Gay-Lussac tied volume and pressure to temperature around 1800; Avogadro proposed in 1811 that equal volumes of gases hold equal numbers of particles. Émile Clapeyron folded them into a single expression in 1834. Kinetic theory later derived the whole thing from mechanics: treat molecules as point masses that bounce elastically and never attract one another, average over their collisions with the walls, and PV=nRTPV = nRT falls out — with RTRT revealed as a measure of the average kinetic energy per mole.

Those two assumptions are also the fine print, and here a common textbook line deserves correcting. The ideal gas law is a limit, not a fact about gases. Real molecules do occupy volume and do attract each other, so the equation is exact only as pressure approaches zero and the gas gets out of its own way. Near condensation it fails plainly: at 100 atm, or anywhere close to the boiling point, the error runs to tens of percent and you want van der Waals or a compressibility factor. Under ordinary room conditions the error is well under 1%, which is why the approximation earns its keep.

Two errors account for most wrong answers on this page, and both are unit errors rather than physics errors. The first is feeding in Celsius. Doubling a gas from 20 °C to 40 °C does not double anything — in kelvin that is 293 to 313, a rise of 7%, and a calculation that used 20 and 40 would be wrong by a factor of nearly two. The second is feeding in a gauge pressure. A tire gauge reading 220 kPa means 321 kPa absolute; PP here is absolute pressure, measured from vacuum, because the equation counts molecular impacts and vacuum is where there are none. A third, quieter trap: the familiar 22.4 L per mole belongs to 0 °C and 1 atm. IUPAC redefined standard pressure to 100 kPa in 1982, and at that pressure the molar volume is 22.71 L. Both numbers circulate, and quoting one against the other's conditions is a 1.3% error hiding inside a memorised constant.

Worked example: 1 mol at 0 C and 1 atm → 22.414 L (molar volume at STP)

Gas Density from Molar Mass

ρ=PMRT\rho = \frac{PM}{RT}
PρMT
Where
  • ρ\rho= Gas density (kg/m³)
  • PP= Pressure (kPa)
  • MM= Molar mass (g/mol)
  • TT= Absolute temperature (°C)

Start from PV=nRTPV = nRT, substitute n=m/Mn = m/M, and rearrange for mass over volume: ρ=PM/RT\rho = PM/RT. What the result says is that a gas has no density of its own. Unlike a solid or a liquid, whose density is close enough to a fixed property to tabulate, a gas takes whatever density its pressure and temperature impose, and the only thing the substance itself contributes is MM. Squeeze it and it gets denser in exact proportion; warm it and it thins in inverse proportion.

Air at 101.325 kPa and 20 °C, using M=0.028964M = 0.028964 kg/mol: ρ=(101325×0.028964)/(8.314×293.15)=1.204 kg/m3\rho = (101\,325 \times 0.028964)/(8.314 \times 293.15) = 1.204\ \text{kg/m}^3 — the figure every ventilation calculation starts from. Helium at the same conditions, with M=0.0040026M = 0.0040026, comes to 0.166 kg/m³. Subtract, and a cubic metre of helium lifts about 1.04 kg. That is the entire physics of a party balloon, and it explains why a balloon large enough to lift a person has to be the size of a house.

Rearranged for molar mass the equation becomes a measurement rather than a prediction, and a historically important one. Weigh a bulb of known volume empty, fill it with a vapour at measured temperature and pressure, weigh it again, and M=ρRT/PM = \rho RT/P hands you the molar mass of an unknown. This is the Dumas method, and through the middle of the nineteenth century it was one of the few routes to a molecular formula. It also settled arguments: measured vapour densities are what showed that many elemental gases travel as diatomic molecules rather than lone atoms.

The dominant error here is the molar mass unit, and it is a clean factor of a thousand. With R=8.314R = 8.314 J/(mol·K) the equation demands MM in kilograms per mole. Air is 0.029 kg/mol. Enter 29 and the answer comes back as 1204 kg/m³ — air denser than water — which at least announces itself. Enter 0.029 when the calculation wanted grams and you get the mirror error. The usual companions apply too: PP must be absolute, not a gauge reading, and TT must be in kelvin.

Two conceptual notes. There is no such thing as "the molar mass of air" in the strict sense — air is a mixture, and 28.96 g/mol is a mole-weighted average of nitrogen, oxygen and argon. It works precisely because an ideal gas is indifferent to what its neighbours are; only the total count matters. That same indifference produces a result most people find backwards: humid air is lighter than dry air. Water is 18 g/mol against air's 29, so at a given pressure and temperature every water molecule that joins the mixture has displaced a heavier one. Muggy days are low-density days, which is why aircraft performance charts include humidity and why a hot, humid runway is a long takeoff.

Worked example: O2 at STP (1 atm, 273.15 K) → 1.42768 kg/m3

Compressibility Factor (Z = PV/nRT)

Z=PVnRTZ = \frac{P V}{n R T}
PVn, TZ
Where
  • ZZ= Compressibility factor
  • PP= Pressure (kPa)
  • VV= Volume (L)
  • nn= Amount of gas (mol)
  • TT= Temperature (°C)

Z is the honest answer to "how wrong is the ideal gas law here?" expressed as a single multiplier. A cylinder holding 2 mol at 400 K reading 1.20 MPa in 5.00 L gives Z=6000/6651.6=0.902Z = 6000/6651.6 = 0.902, meaning the gas occupies about 10% less volume than an ideal gas would at the same pressure and temperature. Below 1, attraction is winning. Above 1, the molecules' own bulk is. At low pressure everything tends back to 1, which is why the ideal gas law survives at all.

The elegant part is the theorem of corresponding states: plot Z against reduced pressure Pr=P/PcP_r = P/P_c at fixed reduced temperature Tr=T/TcT_r = T/T_c and almost every non-polar gas falls on the same generalised chart. Nelson and Obert drew those charts in 1954 and process engineers still read them, because one chart covers nitrogen, methane, propane and argon alike. Natural gas metering leans on this heavily, where the correction is called supercompressibility and moves the invoice by percent-level amounts on a high-pressure line.

Two cautions. Z is not a property you can look up for a substance, only for a substance at a state, so quoting "the Z of methane" without a pressure and temperature says nothing. And when you use Z to correct a flow measurement, be certain which pressure and temperature the meter reports at, because applying a Z evaluated at line conditions to a volume already corrected to standard conditions double-counts the correction and is a classic custody-transfer dispute.

Worked example: An ideal mole at STP gives Z = 1.000

Van der Waals Equation of State

(P+an2V2)(Vnb)=nRT\left(P + \frac{a n^{2}}{V^{2}}\right)\left(V - n b\right) = n R T
PVa, bn, T
Where
  • PP= Pressure (kPa)
  • VV= Volume (L)
  • nn= Amount of gas (mol)
  • TT= Temperature (°C)
  • aa= Attraction constant a (Pa·m⁶/mol²)
  • bb= Excluded volume b (L/mol)

Johannes van der Waals earned the 1910 Nobel Prize for two corrections to the ideal gas law. The an2/V2a n^2/V^2 term adds back the pressure that molecular attraction steals, and the nbnb term subtracts the space the molecules themselves occupy. Put one mole of carbon dioxide in a one-litre vessel at 300 K and the ideal gas law promises 2494 kPa. Van der Waals, with a=0.3640a = 0.3640 and b=4.267×105b = 4.267\times 10^{-5}, says 2242 kPa. Attraction alone accounts for 364 kPa of that, and the real measured value sits near the van der Waals figure.

Expand the equation and it is a cubic in volume, which is the whole reason there is no clean formula for VV. Below the critical temperature the cubic has three real roots. The largest is the saturated vapour, the smallest the saturated liquid, and the middle one is physically unstable and corresponds to nothing you can put in a vessel. This page's volume brain brackets the largest root and bisects for it, so the answer it gives is the vapour branch. If you are working near or below the critical point and need the liquid root, the equation is telling you something the page cannot: use a proper equation of state.

Watch the units on aa and bb, because handbooks quote them in at least four conventions. This page wants aa in Pa·m⁶/mol² and bb in m³/mol. The common table entry of 3.640 bar·L²/mol² for CO₂ is the same number as 0.3640 Pa·m⁶/mol², and 0.04267 L/mol is 4.267×10⁻⁵ m³/mol. The factor of ten between 3.640 and 0.3640 has ruined a great many homework sets.

Worked example: 1 mol CO2 in 1.000 L at 300 K → 2241.5 kPa, not the ideal 2494.3

Antoine Equation (Vapour Pressure)

log10P=ABC+T\log_{10} P = A - \frac{B}{C + T}
TPA, B, C
Where
  • PP= Vapour pressure (kPa)
  • TT= Temperature (°C)
  • AA= Antoine constant A
  • BB= Antoine constant B (°C)
  • CC= Antoine constant C (°C)

Antoine's 1888 fit is three fitted constants standing in for a whole thermodynamic derivation, and it works startlingly well. With the water constants A=8.07131A = 8.07131, B=1730.63B = 1730.63, C=233.426C = 233.426, the equation says that at 100 °C the vapour pressure is 108.071311730.63/333.426=760.1 mmHg10^{8.07131 - 1730.63/333.426} = 760.1\ \text{mmHg}. That is one atmosphere to within a tenth of a millimetre of mercury, recovered from a curve fit rather than from steam tables.

The trap is the units, and it swallows more students than any other feature of the equation. Published constants are overwhelmingly quoted for pressure in mmHg and temperature in degrees Celsius. There is a second convention in kPa and kelvin, and one in bar, and the constants are not interchangeable: use a mmHg set with a temperature in kelvin and you will land several orders of magnitude out with no warning. This page assumes the mmHg and Celsius convention and converts your entered pressure and temperature to match, so enter the constants exactly as the handbook prints them.

The other thing worth knowing is that every constant set has a stated temperature range, and it is not a suggestion. Water has one set for 1-100 °C and a different set for 99-374 °C, and the two disagree by several percent where they overlap. Extrapolate past the range and the error grows without any hint on the page that anything is wrong, because a smooth curve keeps returning smooth numbers. Check the range before you trust the answer.

Worked example: Water at 100 degC → 760.09 mmHg (101.337 kPa)

Clausius–Clapeyron Equation (Two-Point Form)

ln ⁣(P2P1)=ΔHvapR(1T21T1)\ln\!\left(\frac{P_2}{P_1}\right) = -\frac{\Delta H_{vap}}{R}\left(\frac{1}{T_2} - \frac{1}{T_1}\right)
T1P1T2P2ΔH
Where
  • P1P_1= Vapour pressure at T₁ (kPa)
  • T1T_1= Temperature 1 (°C)
  • P2P_2= Vapour pressure at T₂ (kPa)
  • T2T_2= Temperature 2 (°C)
  • ΔHvap\Delta H_{vap}= Enthalpy of vaporisation (kJ/mol)

Integrate the Clapeyron equation with the two honest simplifications, that the vapour is ideal and that ΔHvap\Delta H_{vap} does not change over the interval, and this two-point form falls out. It is the workhorse for turning two boiling points into an enthalpy of vaporisation, or one boiling point into all the others. A liquid boiling at 373.15 K under one atmosphere and at 354.75 K under half an atmosphere has ΔHvap=Rln2/(1/354.751/373.15)=41.5 kJ/mol\Delta H_{vap} = R\ln 2 / (1/354.75 - 1/373.15) = 41.5\ \text{kJ/mol}, which is water to within a couple of percent.

Run it forward on water from 100 °C down to 80 °C with ΔHvap=40.7 kJ/mol\Delta H_{vap} = 40.7\ \text{kJ/mol} and it predicts 48.2 kPa. Steam tables say 47.4 kPa. That 1.7% gap is not an arithmetic error, it is the cost of holding ΔH\Delta H constant across 20 K, and it grows fast as you widen the interval. Near the critical point ΔHvap\Delta H_{vap} collapses towards zero and the equation fails outright.

The unglamorous mistake is temperature units. Both temperatures go in as absolute values, because they appear as 1/T1/T and the reciprocal of a Celsius reading is meaningless. This page takes any temperature unit and converts, but if you are working the equation on paper, convert to kelvin first. The unexpected use, incidentally, is in the kitchen and on mountains: the same equation, run backwards, tells you that at 3000 m water boils near 90 °C, which is why high-altitude cooking directions exist.

Worked example: Water at 80 degC → 48.20 kPa from the 100 degC point

Steam, Turbines & Power Cycles

Saturation Temperature and Pressure of Steam

Tsat=Ts(psat)psat=ps(Tsat)T_{sat} = T_s(p_{sat}) \qquad p_{sat} = p_s(T_{sat})
psatTsat
Where
  • psatp_{sat}= Saturation pressure (absolute) (kPa)
  • TsatT_{sat}= Saturation temperature (°C)

Steam has exactly one degree of freedom while it is boiling. Fix the pressure and the temperature is decided for you, and fix the temperature and the pressure is decided for you, which is why the first column of every steam table is really a single curve with two labels. At 100 kPa absolute water boils at 99.606 °C, at 1 MPa it boils at 179.886 °C, and at 100 °C the vapour pressure is 101.418 kPa. Those three numbers are the whole table in miniature, and the last one is worth staring at: standard atmosphere is 101.325 kPa, so the round 100 °C boiling point everybody learns in school is about a tenth of a kilopascal off the real one. On a steam fitter's gauge the same curve reads 15 psig for 250 °F and 100 psig for 338 °F.

The engine underneath is IAPWS-IF97 Region 4, and it is genuinely two-way. Equation (29) of the release is one implicit quadratic in β=(p/1 MPa)1/4\beta = (p/1\ \text{MPa})^{1/4} and ϑ=T+n9/(Tn10)\vartheta = T + n_9/(T - n_{10}), and Equations (30) and (31) are that same quadratic solved once for pressure and once for temperature. Nothing here iterates and nothing interpolates between printed rows, so a round trip through the pair comes back to the digit it started on. IF97 is the industrial formulation, released in 1997 and adopted by ASME, so the numbers this page returns are the numbers the ASME tables are printed from rather than an approximation of them. Two honest limits. Region 4 stops dead at both ends: 0 °C and 611.213 Pa at the bottom, and the critical point at 373.946 °C and 22.064 MPa at the top, above which liquid and vapour stop being different substances and there is no line left to sit on. And the pressure is ABSOLUTE. A 100 psig gauge is 114.7 psia, and feeding the gauge reading straight in is the single most common way to get a saturation temperature 30 degrees too low.

Worked example: R7-97 Table 36: T_sat at 0.1 MPa → 99.606 °C

Steam Quality from Enthalpy

x=hhfhfgx = \frac{h - h_f}{h_{fg}}
hxhfghfhf + hfg
Where
  • xx= Steam quality (dryness fraction)
  • hh= Enthalpy of the wet steam (J/kg)
  • hfh_f= Saturated liquid enthalpy (J/kg)
  • hfgh_{fg}= Latent heat of vaporisation (J/kg)

Dryness fraction, quality, x. It is the mass of vapour divided by the total mass of the mixture, so x = 0.95 means 950 g of every kilogram is actually steam and the other 50 g is still water riding along as fog. The arithmetic is a lever: saturated water at that pressure sits at hfh_f, dry saturated steam sits at hf+hfgh_f + h_{fg}, and a mixture sits proportionally between them. At 1 MPa absolute the table gives hf=762.5h_f = 762.5 and hfg=2,014.6h_{fg} = 2{,}014.6 kJ/kg, so a sample measured at 2,576 kJ/kg works out to (2,576 − 762.5)/2,014.6 = 0.900. Ninety percent dry, which in a plant is poor.

Two things matter more than the sum. The first is that both table values have to be read at the SAME pressure as the sample, because hfh_f and hfgh_{fg} each move hundreds of kilojoules across a normal working range, and mixing rows is how a perfectly good measurement turns into a nonsense answer. The second is what the answer costs. Wet steam gives up only the fraction x of the latent heat it should, so a coil fed 90% steam is short about 10% of its rating, and the water it carries is doing damage on the way through: droplets at header velocity erode valve seats and turbine blades, and slugs of it are what causes water hammer. In practice quality gets measured with a throttling calorimeter, which expands a small sample to atmospheric pressure where it becomes superheated and its enthalpy can be read from a thermometer alone. And note the honest boundary on this whole page. Quality is only defined ON the saturation line. Below hfh_f the water is subcooled and above hf+hfgh_f + h_{fg} the steam is superheated, and in both cases the answer this page hands back is a message rather than a number worth using.

Worked example: h = 2,576.0 kJ/kg at 1 MPa → x = 0.9002

Wet Steam Enthalpy from h_f and h_g

h=(1x)hf+xhgh = (1 - x) \, h_f + x \, h_g
hghfx1 − xh
Where
  • hh= Enthalpy of the wet steam (J/kg)
  • xx= Steam quality (dryness fraction)
  • hfh_f= Saturated liquid enthalpy (J/kg)
  • hgh_g= Dry saturated steam enthalpy (J/kg)

Same mixture, written from the other pair of columns. Most tables print hfh_f, hfgh_{fg} and hgh_g side by side, and since hg=hf+hfgh_g = h_f + h_{fg} you can blend the two end states directly: a fraction x of the mass arrives as dry steam at hgh_g and the remaining (1 − x) as water at hfh_f. At 1 MPa, with hf=762.5h_f = 762.5 and hg=2,777.1h_g = 2{,}777.1 kJ/kg, steam that is 95% dry carries 0.05 × 762.5 + 0.95 × 2,777.1 = 2,676.4 kJ/kg. Run it the other way and it is the lever rule from school: x = (h − hfh_f)/(hgh_ghfh_f), the distance you have travelled along the horizontal tie line divided by its full length.

Which form you use is purely about which columns are in front of you, and the answers agree to the last digit. The reason to know both is that the same weighting works on every other property in the table, not just enthalpy. Specific volume is the one that bites: at 1 MPa water occupies 0.001127 m³/kg and dry steam 0.19444, a ratio of 173, so a mixture at x = 0.95 has almost exactly 95% of the dry-steam volume and a pipe sized on the mixture is a pipe sized on the vapour. Entropy blends the same way, which is how the exhaust point of a turbine expansion gets pinned down. One thing the lever rule cannot do is average temperatures, because in the wet region there is only one temperature and both ends of the lever are already sitting at it.

Worked example: x = 0.95 between h_f = 762.5 and h_g = 2,777.1 → h = 2,676.37 kJ/kg

Degrees of Superheat

ΔTsh=TTs(p)\Delta T_{sh} = T - T_s(p)
ΔTshpT
Where
  • ΔTsh\Delta T_{sh}= Degrees of superheat ()
  • TT= Steam temperature (°C)
  • pp= Steam pressure (absolute) (kPa)

Superheat is simply how far the thermometer sits above the boiling point at the pressure the steam is actually at. Steam at 1 MPa absolute and 250 °C has 250 − 179.886 = 70.1 degrees of superheat. In trade units, 100 psig at 500 °F is 500 − 337.9 = 162 °F of superheat. The number is meaningless without the pressure beside it, which is the mistake I see most often on log sheets: somebody records 400 °F and calls it superheated without noting that the line is at 250 psig, where saturation is already 406 °F and the steam is in fact slightly wet.

Superheat buys two things and costs a third. It buys dryness, because a superheated main can lose heat through its lagging for a long way before it reaches saturation and starts making condensate, which is why long distribution runs and every steam turbine want it. It also buys turbine life, since expansion starting well above the saturation line ends with less moisture in the last stages. What it costs is heat transfer. A superheated vapour film has a poor film coefficient compared with condensing steam, so a process heater fed 150 degrees of superheat spends the first stretch of its surface doing almost nothing while the steam desuperheats, and the exchanger behaves as though it were undersized. That is exactly why plants run a desuperheating station upstream of process users while sending the same header to the turbine untouched. A negative answer here is worth taking seriously rather than shrugging at: it means the temperature is below saturation, so either the sample is wet or the pressure and temperature were taken at two different places in the system.

Worked example: 250 °C at 1 MPa → 70.114 C° of superheat

Flash Steam Percentage

%F=hf1hf2hfg2×100\%F = \frac{h_{f1} - h_{f2}}{h_{fg2}} \times 100
hf1%Fhfg2hf2
Where
  • %F\%F= Flash steam percentage (%)
  • hf1h_{f1}= Liquid enthalpy at high pressure (J/kg)
  • hf2h_{f2}= Liquid enthalpy at low pressure (J/kg)
  • hfg2h_{fg2}= Latent heat at low pressure (J/kg)

Condensate at 100 psig sits at 338 °F, and it can only do that under pressure. Open a trap to an atmospheric receiver and the water is suddenly 126 °F hotter than it is allowed to be, so it borrows its own excess sensible heat as latent heat and part of it boils instantly. The bookkeeping is pure energy balance: the surplus liquid enthalpy hf1hf2h_{f1} - h_{f2} divided by the latent heat available at the lower pressure. With hf1=309h_{f1} = 309, hf2=180.2h_{f2} = 180.2 and hfg2=970.3h_{fg2} = 970.3 BTU/lb, 128.8/970.3 = 13.3% of the condensate flashes.

Thirteen percent by mass is more than a thousand percent by volume, which is why a receiver vent that "blows steam" is usually normal flash and not a failed trap — the classic misdiagnosis on a steam survey. Tell them apart by watching the plume: flash pulses with each trap discharge, live steam blows continuously. Better still, do not vent it. A flash tank recovering that steam into a low-pressure header, or a heat exchanger putting it into the makeup, pays back in months, and the same arithmetic applied to boiler blowdown is where blowdown heat recovery gets its numbers.

Worked example: 100 psig condensate to atmosphere → 13.27% flash

Flash Steam Mass Rate

m˙f=%F100m˙c\dot{m}_f = \frac{\%F}{100} \, \dot{m}_c
cf%F
Where
  • m˙f\dot{m}_f= Flash steam rate (kg/h)
  • %F\%F= Flash steam percentage (%)
  • m˙c\dot{m}_c= Condensate rate (kg/h)

The companion arithmetic to the flash percentage, and the step that turns a percentage into something you can size a vent or a flash vessel on. Condensate at 100 psig leaving a trap into an atmospheric receiver flashes about 13.3% of its own mass, so a trap passing 2,000 lb/h of condensate releases 266 lb/h of flash steam and drops 1,734 lb/h of hot water into the receiver. In metric, 1,000 kg/h of condensate at the same conditions gives 133 kg/h of flash and 867 kg/h of water.

The mass split looks modest and the volume split does not, which is why this calculation matters more than it first appears. At atmospheric pressure that flash steam occupies 1.673 m³/kg while the condensate beside it occupies 0.00104, so the 13.3% by mass is more than 1,600 times the volume of the water it came out of. A return line sized on the condensate alone will be choked by its own flash, back-pressure the traps and stall the equipment upstream, and this is the classic reason a heat exchanger will not hold temperature on a cold morning. Two habits. Use it to argue for recovery rather than venting, since a flash vessel feeding a low-pressure header, or a heat exchanger putting that energy into makeup water, usually pays back inside a year. And use it as a diagnostic: a receiver vent plume that pulses with each trap discharge is honest flash, while one that blows steadily is a failed trap passing live steam, and running this calculation tells you which plume the numbers say you should be seeing.

Worked example: 133 kg/h of flash from 1,000 kg/h of condensate → 13.3 %

Steam Turbine Specific Work

w=h1h2w = h_1 - h_2
wh1h2
Where
  • ww= Specific work (J/kg)
  • h1h_1= Inlet enthalpy (J/kg)
  • h2h_2= Exhaust enthalpy (J/kg)

Write the steady-flow energy equation across a turbine, drop the kinetic and potential terms because they are tiny, and drop the heat loss because a lagged casing is adiabatic to about one percent, and what is left is beautifully plain: the work out per kilogram is the enthalpy in minus the enthalpy out. Steam entering at 4 MPa and 400 °C carries 3,214 kJ/kg. If it leaves at 2,600 kJ/kg, the machine has taken 614 kJ out of every kilogram that went through, and no property other than those two enthalpies is involved.

Turbine people usually flip that into steam rate, which is the mass of steam needed per unit of work: 3,600 divided by 614 kJ/kg gives 5.86 kg per kilowatt-hour. It is the number a plant engineer carries in his head because it converts a load directly into a boiler duty. Two cautions. The exhaust enthalpy h2h_2 is rarely measurable, because a wet exhaust has no unique temperature at its pressure, so in practice h2h_2 is worked back from a measured power output or from an assumed isentropic efficiency rather than read off a thermometer. And this is specific work at the steam path, not at the coupling. Bearing friction, the governor, the oil pump and windage all come off before the shaft turns anything, and on a small machine those can be a few percent.

Worked example: 3,213.6 down to 2,600.0 kJ/kg → w = 613.6 kJ/kg

Steam Turbine Power Output

P=m˙wP = \dot{m} \, w
wP
Where
  • PP= Shaft power (W)
  • m˙\dot{m}= Steam mass flow (kg/h)
  • ww= Specific work (J/kg)

Specific work is per kilogram, so multiply by the kilograms per second and you have watts. A machine passing 10 kg/s, which is 36,000 kg/h or about 79,000 lb/h, at 614 kJ/kg develops 6,140 kW at the shaft. That single multiplication is the bridge between the thermodynamics and the electrical single line, and it is also the bridge back: divide a required 5 MW by the same 614 kJ/kg and the boiler has to deliver 8.14 kg/s of steam at the throttle conditions, which is what actually sizes the plant.

Be clear about where the boundary sits. This is shaft power at the coupling, before the gearbox, before the generator and before the excitation. A generator at 96% and a gearbox at 98% turn that 6,140 kW into about 5,780 kW at the terminals, and nameplates are usually quoted at the terminals, so comparing a calculated shaft figure against a nameplate without accounting for the drivetrain will always make the machine look bad. The other thing to keep straight is which kind of turbine you are costing. A back-pressure machine exhausts into a process header, so its steam is not wasted and the work is nearly free, while a condensing machine dumps its exhaust latent heat into the cooling tower and only converts a fraction of the fuel into work. Same equation, completely different economics.

Worked example: 36,000 kg/h at 613.6 kJ/kg → 6.136 MW

Turbine Isentropic Efficiency

ηisen=h1h2h1h2s\eta_{isen} = \frac{h_1 - h_2}{h_1 - h_{2s}}
h1h2sh2ηisen
Where
  • ηisen\eta_{isen}= Isentropic efficiency
  • h1h_1= Inlet enthalpy (J/kg)
  • h2h_2= Actual exhaust enthalpy (J/kg)
  • h2sh_{2s}= Isentropic exhaust enthalpy (J/kg)

The ideal turbine expands at constant entropy. It is the best a machine can possibly do between a given inlet state and a given exhaust pressure, it is a vertical line on the Mollier chart, and it is the yardstick everything else is measured against. Work the example through. Steam at 4 MPa and 400 °C has h = 3,214 kJ/kg and s = 6.769 kJ/(kg·K). Expand to 0.1 MPa holding that entropy: at 0.1 MPa the table gives sf=1.3026s_f = 1.3026 and sfg=6.0568s_{fg} = 6.0568, so the ideal end point has quality (6.769 − 1.3026)/6.0568 = 0.9025, and its enthalpy is 417.5 + 0.9025 × 2,258.0 = 2,455 kJ/kg. The ideal drop is therefore 759 kJ/kg. If the machine really delivers 614 kJ/kg, its isentropic efficiency is 614/759 = 0.81.

Eighty-one percent is a normal answer. A modern multistage condensing turbine runs 80 to 90%, a well-matched back-pressure machine 70 to 80%, and a small single-stage unit driving a pump may only manage 50 to 65%. Friction, leakage past the blade tips, throttling at the governor valve and moisture drag all show up in that one number, and the moisture part follows the old Baumann rule of roughly one percent of efficiency lost per one percent of average exhaust wetness. If this page hands you a figure above 100%, treat it as information rather than a windfall: no expansion beats the isentropic one, because irreversibility only ever adds entropy and pushes the real end point to the RIGHT on the chart, which means to a higher enthalpy. An answer over 100% says the state points are wrong. The usual culprits are h2sh_{2s} evaluated at the wrong exhaust pressure, an entropy carried at too few figures, or an exhaust temperature measured downstream of a gland leak or a spray.

Worked example: 4 MPa/400 °C expanding to 0.1 MPa, 613.6 of 758.24 kJ/kg → η = 0.8092

Napier's Steam Leak Rate

m˙=AP70\dot{m} = \frac{A \, P}{70}
PA
Where
  • m˙\dot{m}= Steam leak rate (kg/h)
  • AA= Orifice area (mm²)
  • PP= Steam pressure (absolute) (kPa)

Read the unit convention first, because the bare 70 means nothing without it: Napier's rule is ṁ in POUNDS PER SECOND, A in SQUARE INCHES, and P in POUNDS PER SQUARE INCH ABSOLUTE. Per hour the same statement is the ASME Section I relief-capacity form ṁ = 51.5 A P, since 3,600/70 is 51.43. The solver converts for you, so you can type square millimetres and kilopascals and get kilograms an hour back, but the constant belongs to those three units and quoting the 70 against pounds per hour understates a leak by a factor of 3,600. Now the number everybody actually wants. A 1/8 inch hole is 0.0123 in². At 100 psig, which is 114.7 psia, Napier gives 0.0123 × 114.7 / 70 = 0.0201 lb/s, or 72 lb/h. Left alone for a year that is 634,000 lb of steam, about 288 tonnes, and at a typical $12 per 1,000 lb it is roughly $7,600 walking out of one hole you could cover with a thumbnail.

Napier's rule is an approximation and it is worth knowing exactly which one. It is a straight-line fit to CHOKED flow, meaning the steam is moving at its own speed of sound in the throat, which happens whenever the downstream pressure is below about 58% of the upstream absolute. Above that ratio the orifice is not choked, the flow depends on the downstream pressure too, and the rule reads high. The fit was made for saturated steam and is honest to a few percent from roughly 25 psia up to about 1,000 psia, drifting low above that. As a check on the constant itself, an ideal isentropic choked-flow calculation for that same 1/8 inch hole at 114.7 psia gives 69.8 lb/h against Napier's 72.4, which is the few percent high the rule is known for. Two field cautions. The area is the hole, not the pipe, and holes are rarely round or sharp-edged, so on a corroded gasket face or a worn valve seat you are estimating an equivalent area and the answer carries that uncertainty with it. And a plume is not a hole: plumes look the same at very different flows, so estimate the leak size from the metal, not from the cloud.

Worked example: 10 mm² hole at 1,000 kPa absolute → 52.443 kg/h

Desuperheater Water Injection Rate

m˙w=m˙1h1h2h2hw\dot{m}_w = \dot{m}_1 \, \frac{h_1 - h_2}{h_2 - h_w}
1h1whwh2
Where
  • m˙w\dot{m}_w= Spray water rate (kg/h)
  • m˙1\dot{m}_1= Superheated steam rate (kg/h)
  • h1h_1= Inlet steam enthalpy (J/kg)
  • h2h_2= Target outlet enthalpy (J/kg)
  • hwh_w= Spray water enthalpy (J/kg)

A desuperheating station is a mass balance and an energy balance solved together, and nothing more exotic. Steam comes in at m˙1\dot{m}_1 and h1h_1, spray water joins it at m˙w\dot{m}_w and hwh_w, and what leaves is the sum of the two masses carrying the sum of the two energies. Set the mixed enthalpy to the target h2h_2, rearrange, and the water you need is m˙w=m˙1(h1h2)/(h2hw)\dot{m}_w = \dot{m}_1 (h_1 - h_2)/(h_2 - h_w). The numerator is the heat you have to take out per kilogram of steam and the denominator is how much heat each kilogram of water can absorb before it too reaches the target, so the answer is one divided by the other. A worked case: 10 kg/s of steam at 3,300 kJ/kg to be brought to 2,900 kJ/kg with feedwater at 700 kJ/kg needs 10 × 400/2,200 = 1.82 kg/s of spray, which is 18% more mass leaving the station than entered it.

Enthalpies do the work here, not temperatures, and that is deliberate. Specific heat is not constant across a desuperheating range, and the outlet may end up anywhere from heavily superheated to saturated, so a mcΔTmc\Delta T balance quietly drifts while an enthalpy balance stays exact. Three practical notes. Everything at the outlet is at the DOWNSTREAM pressure, so if the station also drops pressure, read h2h_2 at the pressure the steam leaves at rather than the pressure it arrived at. Do not aim at saturation: spray control needs superheat left over to prove the water has all evaporated, and 10 to 15 degrees is the usual minimum, because a target on the saturation line means any overshoot puts liquid water into a hot pipe. And the extra mass is real, so the downstream line, the safety valve capacity and the process user all see more steam than the upstream meter reads, which is a good thing when the water is treated feedwater and a very bad thing when it is not, since every dissolved solid in that spray goes straight into the header.

Worked example: 10 kg/s from 3,213.6 to 2,800.8 kJ/kg with 632.2 kJ/kg water → 6,852.7 kg/h

Steam Coil Condensate Load

m˙=Q˙hfg\dot m = \frac{\dot Q}{h_{fg}}
Where
  • m˙\dot m= Condensate load (kg/h)
  • Q˙\dot Q= Coil duty (W)
  • hfgh_{fg}= Latent heat at coil pressure (J/kg)

A steam coil does its work by condensing, so the steam it eats is simply the duty divided by the latent heat at the pressure the coil is running at. A 150 kW preheat coil fed 15 psig steam, where hfgh_{fg} is about 2,199 kJ/kg, condenses 150/2,199 = 0.0682 kg/s, which is 246 kg/h or roughly 540 lb/h. Notice that the answer is not really a flow of steam so much as a flow of water: everything that condenses has to leave through the trap, which is why the same number is called the condensate load and why it is the figure the trap, the return line and the receiver all get sized from.

The pressure matters twice over and in opposite directions. Raising it gives a hotter coil and a bigger temperature difference, so the coil transfers more heat per square foot, but it also LOWERS the latent heat each kilogram carries, so the mass flow climbs faster than the duty does. That is the trade-off behind every steam coil selection.

Then there is the number this page does not give you, and it is the one that matters on a cold morning. The running load calculated here is a steady-state figure. Start-up is worse, sometimes far worse: the coil metal is cold, the air on it is cold, and condensate forms faster than the coil will ever produce it once warm. Traps are selected on two to three times the running load for that reason. Undersize the trap and condensate backs up into the tubes, where it subcools, and in a coil facing outdoor air in a Canadian winter it freezes and splits them — which is exactly why steam-distributing "non-freeze" coils exist, with an inner perforated tube that keeps live steam against the full length of the outer tube instead of letting a stalled slug sit in the bottom.

Worked example: A 210 kW coil on steam of 2,100 kJ/kg condenses 0.1 kg/s = 360 kg/h

Thermal Efficiency

η=WQh\eta = \frac{W}{Q_h}
QhηW
Where
  • η\eta= Efficiency
  • WW= Useful work output (J)
  • QhQ_h= Heat input (J)

Thermal efficiency is a bookkeeping question with an uncomfortable answer: of all the heat you paid for, how much came back out as work? η=W/Qh\eta = W/Q_h, a ratio of two energies, so it is dimensionless and sits between 0 and 1. The uncomfortable part is that it is never close to 1, and not because engineers have failed. A heat engine cannot convert all of QhQ_h into work no matter how well it is built, because to run a cycle at all it has to dump some heat into a colder place. That requirement is the second law, and it is why the missing energy is not waste in the sense of carelessness — it is rent.

A 1000 MW thermal station delivering 380 MW to the grid runs at η=380/1000=0.38\eta = 380/1000 = 0.38. The other 620 MW is not lost in any physical sense; it leaves at low temperature through the condenser and has to go somewhere real — a river, a lake, or the plume off a cooling tower. That 620 MW is the reason large thermal plants are built beside water, and the reason a heat wave can force one to throttle back.

What sets the ceiling is on the neighbouring Carnot page: ηmax=1Tc/Th\eta_{max} = 1 - T_c/T_h, with both temperatures absolute. A steam cycle at 550 °C rejecting to a 27 °C condenser has a Carnot limit of 1300/823=63.5%1 - 300/823 = 63.5\%, so a real 38% plant is achieving about six-tenths of what thermodynamics allows. That ratio — actual over Carnot — is the honest measure of engineering quality, and it is far more flattering than the raw efficiency. It also explains why every gain in this field comes from raising ThT_h: combined-cycle plants reach 60% by putting a gas turbine at 1400 °C in front of the steam cycle and feeding the steam cycle its exhaust.

The error that matters most is comparing two efficiencies computed on different denominators. For fuel-burning equipment, QhQ_h can be the higher heating value, which counts the energy recovered when the combustion water vapour condenses, or the lower heating value, which does not. The two differ by about 10% for natural gas. This is why condensing boilers are sometimes advertised above 100% efficiency: that figure is on LHV, and it breaks no law — it is a fraction with a deliberately small denominator. A 95% AFUE furnace and a "108%" boiler may be the same machine described twice. Always ask which heating value, especially when comparing a quoted number against one you calculated yourself.

Three more distinctions worth keeping straight. WW must be net work: an engine's output minus whatever the feed pumps, compressors and auxiliaries consume, or you are counting energy that never left the building. QhQ_h is the heat that actually entered the working fluid, which for a fired system is the fuel energy less the stack loss, not the fuel energy itself. And a heat pump's coefficient of performance is not a thermal efficiency and routinely reads 3 or 4 — no contradiction, because a heat pump is not converting heat into work, it is spending work to move heat, and the heat it moves was already there. Putting a COP and an η on the same axis is comparing two different questions.

Worked example: 250 J work from 1000 J heat → η = 0.25

Carnot Efficiency

η=1TcTh\eta = 1 - \frac{T_c}{T_h}
ThηTc
Where
  • η\eta= Maximum efficiency
  • TcT_c= Cold reservoir temperature (°C)
  • ThT_h= Hot reservoir temperature (°C)

In 1824 the 28-year-old French engineer Sadi Carnot proved a startling result: no heat engine, however ingenious, can beat an efficiency ceiling set purely by the absolute temperatures of its hot source and cold sink. His short book went almost unnoticed in his lifetime, yet it founded thermodynamics and anticipated the second law. The formula demands kelvin — formula.expert converts °C and °F inputs — and it only makes sense when 0 < TcT_c < ThT_h.

A steam plant running between 800 K boiler steam and a 300 K cooling river can never exceed η = 1 − 300/800 = 62.5%; real plants manage about 40% after friction and other irreversibilities take their share. The lesson is practical: to buy efficiency, raise ThT_h or lower TcT_c. That is why power stations chase ever-hotter turbine materials — and why waste heat is not an engineering failure but a law of nature.

Worked example: 500 K hot, 300 K cold → η = 0.4

Rankine Cycle Thermal Efficiency

η=(h1h2)(h4h3)h1h4\eta = \frac{\left(h_1 - h_2\right) - \left(h_4 - h_3\right)}{h_1 - h_4}
h1h2h3h4η
Where
  • η\eta= Thermal efficiency
  • h1h_1= Turbine inlet enthalpy (J/kg)
  • h2h_2= Turbine exhaust enthalpy (J/kg)
  • h3h_3= Condensate enthalpy (J/kg)
  • h4h_4= Feedwater enthalpy (J/kg)

William Rankine set out the steam cycle in the 1850s and the arithmetic has not changed: net work over heat in, with every term read off a steam table. Take 3 MPa and 350 °C into the turbine against a 10 kPa condenser. The turbine drops the steam from 3116 to 2136 kJ/kg, the feed pump lifts the condensate by about 3 kJ/kg, and the boiler puts back 2921. That is (979.93.0)/2921=33.4%(979.9 - 3.0)/2921 = 33.4\%, and a real 1950s station would have measured something close.

The pump work is the term everybody wants to drop, and for a first pass they are right. Compressing a liquid costs vΔPv\,\Delta P, which for water lifted from 10 kPa to 3 MPa is about 3 kJ/kg against a turbine output near 980, so it is 0.3% of the gross. This is the whole reason the Rankine cycle uses a condenser at all rather than running a gas back through a compressor: the Brayton cycle spends half its turbine output driving its own compressor, and the Rankine cycle spends three parts in a thousand. Condensing the working fluid before pressurising it is the single best trick in power engineering.

What is genuinely counterintuitive is where the efficiency gains come from. Lowering the condenser pressure buys more than raising the boiler pressure, because it stretches the temperature range at the cheap end, and this is why a coastal plant with cold seawater beats an inland one on identical hardware. It is also why a station's efficiency drops measurably in August. The limit is that a lower condenser pressure means wetter steam at the turbine exhaust, and below about 88% dryness the water droplets erode the last-stage blades, which is what reheat stages exist to prevent.

Worked example: 3 MPa / 350 degC steam to a 10 kPa condenser → 33.44%

Otto Cycle Efficiency (Compression Ratio)

η=11rγ1\eta = 1 - \frac{1}{r^{\gamma - 1}}
rγη
Where
  • η\eta= Thermal efficiency
  • rr= Compression ratio
  • γ\gamma= Heat capacity ratio

The air-standard Otto cycle says something remarkable: the thermal efficiency of a spark-ignition engine depends on nothing but its compression ratio. Not fuel, not speed, not displacement, not how much heat you add. At r=8r = 8 with γ=1.4\gamma = 1.4 the answer is 180.4=56.5%1 - 8^{-0.4} = 56.5\%, and at r=10r = 10 it is 60.2%. Every engineering student meets that result and immediately asks why real petrol engines manage 25 to 35%, which is the right question.

Part of the gap is honest thermodynamics that the air-standard model throws away: real combustion is not instantaneous and not at constant volume, heat leaks into the coolant, exhaust gas leaves hot, and γ\gamma for a hot combustion mixture is nearer 1.3 than 1.4. Recompute at γ=1.3\gamma = 1.3 and the r=8r = 8 figure falls to 48%. The rest is pumping work, friction and part-load throttling, none of which appear anywhere in the equation.

The reason nobody simply raises the compression ratio is knock. Squeeze the charge harder and the end gas ahead of the flame front gets hot enough to autoignite, which produces the metallic rattle and, sustained, holes in pistons. Fuel octane rating is precisely a measure of resistance to this, and it is the ceiling that has held petrol engines near r=10r = 10 for decades. Diesels dodge it entirely by compressing air alone and injecting fuel at the top, which is why they run at 16 to 22 and why they are more efficient. Direct injection and variable valve timing have lately pushed petrol engines to 12 or 13, which is worth a few real points of efficiency.

Worked example: Compression ratio 8 at gamma 1.4 → 56.47%

Brayton Cycle Efficiency (Pressure Ratio)

η=11rp(γ1)/γ\eta = 1 - \frac{1}{r_p^{\left(\gamma - 1\right)/\gamma}}
rpγη
Where
  • η\eta= Thermal efficiency
  • rpr_p= Pressure ratio
  • γ\gamma= Heat capacity ratio

The gas turbine's air-standard efficiency depends only on pressure ratio, in the same way the Otto cycle depends only on compression ratio. At rp=8r_p = 8 with γ=1.4\gamma = 1.4 the answer is 182/7=44.8%1 - 8^{-2/7} = 44.8\%. Modern industrial machines run pressure ratios of 15 to 25 and aero engines go past 40, and the equation duly promises 55 to 65%.

The catch that the equation cannot express is that efficiency is not the only thing that matters. As pressure ratio climbs, the specific work, meaning the output per kilogram of air, rises to a peak and then falls, because the compressor is eating an ever larger share of what the turbine makes. A machine optimised purely for efficiency needs to swallow far more air for the same output, which means a physically larger and more expensive engine. Real designs sit between the pressure ratio that maximises efficiency and the one that maximises specific work, and the choice depends on whether the customer pays more for fuel or for hardware.

Turbine inlet temperature does not appear in this equation at all, which is the biggest thing the air-standard model hides. In practice it is the dominant design variable: raising it lifts both efficiency and specific work, and the entire history of gas turbine development is metallurgy, film cooling and thermal barrier coatings chasing it upward. It is also why the exhaust is still 500 to 600 °C and why combined-cycle plants bolt a Rankine steam cycle on the back, reaching over 60% overall by making the gas turbine's waste heat someone else's fuel.

Worked example: Pressure ratio 8 at gamma 1.4 → 44.80%

Conduction & the Building Envelope

Heat Conduction Rate

P=kAΔTdP = \tfrac{k A \Delta T}{d}
ΔTkAPd
Where
  • PP= Heat flow rate (W)
  • kk= Thermal conductivity (W/(m·K))
  • AA= Cross-sectional area ()
  • ΔT\Delta T= Temperature difference ()
  • dd= Thickness (mm)

Fourier's law says heat flows down a temperature gradient at a rate proportional to how steep that gradient is. Written for a flat slab, the gradient is ΔT/d and the flow is P=kAΔT/dP = kA\Delta T/d. Each term earns its place: more area gives the heat more parallel paths, a larger temperature difference drives harder, and greater thickness spreads the same difference over a longer distance and so flattens the gradient. The conductivity kk is the material's own willingness to pass heat along, and it spans four orders of magnitude — copper near 400 W/(m·K), glass 0.96, mineral wool 0.04, still air 0.026.

Run the numbers on a single window pane: 1.5 m², 3 mm thick, 20 K from inside to outside. P=0.96×1.5×20/0.003=9600 WP = 0.96 \times 1.5 \times 20 / 0.003 = 9600\ \text{W}. Nine and a half kilowatts through one window. That is obviously wrong, and being clear about why it is wrong is the most useful thing this page can teach.

The equation is fine; the model is incomplete. Heat has to reach the glass from the room air and leave it into the outdoor air, and both of those handoffs are slow. The still-air film clinging to the inside surface has a thermal resistance of about 0.12 m²·K/W and the wind-scoured outside film about 0.03, while the glass itself contributes only d/k=0.003/0.96=0.0031d/k = 0.003/0.96 = 0.0031. Resistances in series add, so the total is roughly 0.155 m²·K/W, giving U=1/R=6.5 W/(m2⋅K)U = 1/R = 6.5\ \text{W/(m}^2\text{·K)} and a real heat flow of about 195 W — fiftyfold less than the bare slab calculation. The glass was never the bottleneck. This series-resistance picture is exactly the electrical analogy: ΔT is voltage, PP is current, and d/kAd/kA is resistance, which is why building science speaks in R-values (d/kd/k per unit area) and U-values (1/Rtotal1/R_{total}) rather than in this equation directly.

So the standard error is treating one layer as the whole assembly. Add layers by summing their resistances, never by averaging their conductivities, and never forget the two air films — in a well-insulated wall they are negligible, and in a window they dominate. It is also why a double-glazed unit works: the gain comes almost entirely from the trapped gas layer and the two extra surface films, not from doubling the glass.

Three further limits. This is a steady-state relation: it tells you the flow once the temperatures have settled and says nothing whatever about how long a wall takes to get there, which is a question of thermal mass and belongs to the diffusion equation. It is one-dimensional, which means it silently assumes heat goes straight through — a steel stud bridging an insulated cavity carries far more than its share of the area, and a thermal bridge can add a third to an assembly's real loss while the calculated R-value notices nothing. And ΔT is a difference, identical in kelvin and Celsius, but imperial conductivities quoted in BTU·in/(hr·ft²·°F) are a different quantity with a different thickness convention buried in them, so convert deliberately rather than by feel.

Worked example: 2 m² glass pane, 4 mm, ΔT = 15 K → 7200 W

R-Value of an Insulation Layer (R = L/k)

R=LkR = \frac{L}{k}
kLR
Where
  • RR= R-value of the layer (RSI (m²·K/W))
  • LL= Layer thickness (mm)
  • kk= Thermal conductivity (W/(m·K))

An R-value is thermal resistance per unit area: thickness divided by conductivity, with no area term at all. That is why an insulation label can carry one number for a product sold by the roll. 140 mm of mineral wool at k = 0.040 W/(m·K) gives 0.140/0.040 = RSI 3.5, and the same layer is R-19.9 in Canada or the United States — R and RSI are one quantity in two units, and because the imperial unit (h·ft²·°F/BTU) is the smaller one, R ≈ 5.678 × RSI. An "R-20 wall" and an "RSI-3.5 wall" are the same wall. Divide by 5.678 to go the other way; the calculator will do it for you, but the factor is worth memorising because half the arguments on a job site are really about this multiplier.

Rules of thumb by product, per inch of thickness: fibreglass and mineral wool batt R-3.1 to R-4.3, blown cellulose R-3.2 to R-3.8, expanded polystyrene R-3.6 to R-4.2, polyisocyanurate R-5 to R-6, and closed-cell spray foam about R-6. Softwood framing manages only R-1.25 per inch and concrete about R-0.08 — which is why an 8 in concrete wall is R-0.64 and needs everything it can get on the outside. The traps are all in k. It rises with temperature for fibrous insulation and falls for foams as the blowing agent diffuses out, so aged polyiso is derated; it collapses when the material gets wet, because water conducts about 25 times better than the trapped air it replaced; and compressing a batt into a shallower cavity costs you resistance in proportion, so an R-19 batt stuffed into a 3.5 in bay is delivering closer to R-13.

Worked example: 140 mm mineral wool at k=0.040 → RSI 3.50

Thermal Resistance of a Plane Wall

R=LkAR = \frac{L}{k A}
AkLR
Where
  • RR= Conduction resistance (K/W)
  • LL= Wall thickness (mm)
  • kk= Thermal conductivity (W/(m·K))
  • AA= Cross-sectional area ()

Joseph Fourier spent the years after 1807 being told by Lagrange and Laplace that his heat series were not rigorous; the Théorie analytique de la chaleur finally appeared in 1822 and gave physics both the conduction law and the Fourier series. Written as a resistance, his law becomes R = L/(kA) — thickness over conductivity times area — and once heat flow wears the clothes of Ohm's law you can stack, branch and total resistances exactly as an electrician would. A 100 mm concrete wall (k = 1.4 W/(m·K)) of 12 m² has R = 0.1/(1.4 × 12) = 0.00595 K/W, so 22 K across it drives 22/0.00595 ≈ 3700 W.

Two traps. First, this is the absolute resistance in K/W, not the building-trade R-value, which is per unit area (m²·K/W, or h·ft²·°F/BTU in the US) — multiply the RSI by area and invert to compare. Second, k is not a constant: mineral wool at −20 °C is not the mineral wool on the datasheet at 24 °C, and wet insulation can lose three-quarters of its resistance because water conducts 25 times better than the trapped air it displaced. Every insulated cold line that sweats is quietly converting itself into a bare pipe.

Worked example: 100 mm concrete (k=1.4) over 12 m2 → 0.005952 K/W

Thermal Resistances in Series

Rtot=R1+R2+R3R_{tot} = R_1 + R_2 + R_3
R1R2R3Rtot
Where
  • RtotR_{tot}= Total resistance (K/W)
  • R1R_1= First layer resistance (K/W)
  • R2R_2= Second layer resistance (K/W)
  • R3R_3= Third layer resistance (K/W)

Because every layer of a composite wall passes the identical heat flow, their resistances add just like series resistors: Rtot=R1+R2+R3R_{\mathrm{tot}} = R_1 + R_2 + R_3. Set an unused layer to zero. The sum immediately shows where the money goes — an inside film of 0.02, a 6 mm steel skin of 0.006 and 50 mm of mineral wool at 0.04 K/W total 0.066 K/W, of which the steel is 9%. Adding a second steel skin changes almost nothing; adding a second inch of wool changes a great deal.

The same arithmetic exposes the most expensive mistake in insulation work, the thermal bridge. Series resistances add, but parallel paths do not: a steel stud, a through-bolt or a pipe hanger sits alongside the insulation rather than behind it, and a bridge occupying 2% of the area with 300 times the conductivity can carry a third of the heat. That is why fastener manufacturers sell thermal-break washers, and why an infrared camera pointed at a well-insulated wall on a cold morning still draws you a perfect picture of the studs.

Worked example: 0.020 + 0.006 + 0.040 → 0.066 K/W total

Total R-Value of an Assembly

Rtot=R1+R2+R3R_{tot} = R_1 + R_2 + R_3
R1R2R3Rtot
Where
  • RtotR_{tot}= Total assembly R-value (RSI (m²·K/W))
  • R1R_1= First layer R-value (RSI (m²·K/W))
  • R2R_2= Second layer R-value (RSI (m²·K/W))
  • R3R_3= Third layer R-value (RSI (m²·K/W))

Because every layer of a wall passes the same heat through the same square metre, their R-values add: Rtot=R1+R2+R3R_{\mathrm{tot}} = R_1 + R_2 + R_3. Set an unused layer to zero. A conventional 2×6 wall totals R-19 in the cavity, R-5 for an inch of exterior polyiso and about R-1.13 for the inside air film plus gypsum, giving R-25.13 — and note that the two "small" layers together are 24% of the answer, which is why the air films and finishes are in the ASHRAE tables at all rather than being waved away. In metric the same sum runs RSI 0.12 for the inside film, 3.32 of cavity and 0.08 for drywall to make RSI 3.52, the Canadian code number that everyone quotes as R-20.

Two things do not belong in this sum. Airspaces are not insulation in proportion to their thickness: a sealed 20 mm cavity is worth about RSI 0.17 (R-1) and a 100 mm one is worth barely more, because convection takes over as soon as the gap is wide enough to circulate — the exception is a low-emissivity foil facing the gap, which doubles it. And the exterior air film depends on wind, dropping from RSI 0.12 in still air to 0.03 in a 25 km/h breeze, so any R-value quoted with a generous outdoor film is quoting a calm day. The rule to hold on to is that only series layers add. Studs, joists and fasteners run alongside the insulation rather than behind it, and that is a parallel path, handled on the effective R-value page.

Worked example: RSI 3.52 wall less 0.12 film and 0.08 gypsum → cavity RSI 3.32

U-Factor from Total R-Value (U = 1/R)

U=1RtotU = \frac{1}{R_{tot}}
RtotU
Where
  • UU= U-factor (W/(m²·K))
  • RtotR_{tot}= Total assembly R-value (RSI (m²·K/W))

Insulation is sold by resistance and windows are sold by conductance, so the same envelope gets described two ways: R-20 is U-0.05 in imperial units, and RSI 3.52 is U = 0.284 W/(m²·K) in metric. They are reciprocals of each other within one unit system — and note that the conversion between systems is the same 5.678 factor running the other way, since U in W/(m²·K) equals 5.678 divided by the imperial R. Codes mix the two deliberately: opaque assemblies get a prescriptive R, fenestration gets a U-factor on the NFRC label, and the compliance path adds them up as UA products.

Here is the mistake that costs real money, and it appears in spreadsheets everywhere: U-factors add, R-values add, but never in the same direction. R-values add along the path heat takes, in series through the layers. U-factors add across the wall, area-weighted, when parallel paths share the same ΔT — which is why a wall and its windows combine as ΣUA and not as an average R. Averaging R-values across an elevation of wall and glass gives a wall that looks far better than it is: 90% at R-20 and 10% at R-3 does not average to R-18.3; it averages to U = 0.9 × 0.05 + 0.1 × 0.333 = 0.0783, which is R-12.8. The glass ate a third of the wall's performance, and only the U-side arithmetic showed it.

Worked example: R-20 wall → U = 0.2839 W/(m2.K)

Overall U from Total Resistance

U=1RtotAU = \frac{1}{R_{tot} A}
ARtotU
Where
  • UU= Overall coefficient (W/(m²·K))
  • RtotR_{tot}= Total thermal resistance (K/W)
  • AA= Reference area ()

Two trades describe the same wall with reciprocal numbers. Heat-transfer analysis likes absolute resistance in K/W because it adds; equipment datasheets like U in W/(m²·K) because it multiplies straight into Q̇ = UAΔT. This converts between them: U=1/(RtotA)U = 1/(R_{\mathrm{tot}} A). A 12 m² assembly totalling 0.005 K/W has U = 1/(0.005 × 12) = 16.7 W/(m²·K) — a bare metal partition. Building envelopes live three orders of magnitude lower, at U = 0.15–0.35.

The trap is the area you divide by, and it has cost real money. Shell-and-tube exchangers are almost always rated on outside tube area, but the inside film resistance is physically attached to the inside area, which on a 19 mm tube with a 2 mm wall is 21% smaller. Mix the two and your U is off by a fifth, which shows up as an exchanger that will not make its approach on the hottest day of the year. Always write the area basis next to the U — "850 W/(m²·K), outside" — and the ambiguity disappears.

Worked example: 0.005 K/W over 12 m2 → U = 16.67 W/(m2.K)

Heat Flow from Thermal Resistance

Q˙=ΔTR\dot{Q} = \frac{\Delta T}{R}
ΔTQR
Where
  • Q˙\dot{Q}= Heat flow rate (W)
  • ΔT\Delta T= Temperature difference ()
  • RR= Total thermal resistance (K/W)

This is the payoff of the whole resistance network: once you have added up the layers, the heat flow is just ΔT over R, the exact analogue of I = V/R. A wall assembly totalling 0.066 K/W with 22 K across it passes 22/0.066 = 333 W, and no further physics is required. The same equation run backwards is how a thermal engineer sizes an enclosure — an electronics box dissipating 40 W that may only rise 25 K above ambient needs the path from junction to air to be under 0.625 K/W, total, including the heatsink and the interface pad.

Because resistances are additive and heat flow is not, the intermediate temperatures come free: the drop across any single layer is Q̇ × R for that layer alone. That is how you check whether a wall's dew point falls inside the insulation or safely outside it, and how a plant engineer proves the fouled side of an exchanger is the tube interior rather than the shell. Trap: ΔT here is a difference, so 40 °F of difference is 22.2 K, not 4.4 — the calculator handles the conversion, but the arithmetic in your notebook may not.

Worked example: 22 K across 0.066 K/W → 333.3 W

Heat Loss Through an Assembly (Q = A·ΔT/R)

Q˙=AΔTRtot\dot{Q} = \frac{A \, \Delta T}{R_{tot}}
ΔTAQRtot
Where
  • Q˙\dot{Q}= Heat loss rate (W)
  • AA= Assembly area ()
  • ΔT\Delta T= Inside-to-outside ΔT ()
  • RtotR_{tot}= Total assembly R-value (RSI (m²·K/W))

This is the line item that appears once per surface on every heat-loss worksheet ever drawn up: area times ΔT, divided by R. A 30 m² wall at RSI 3.5 with 22 K across it loses 30 × 22/3.5 = 189 W. In imperial the same arithmetic stays in imperial — 200 ft² of R-20 wall with 40 °F across it loses 200 × 40/20 = 400 BTU/h, which is 117 W — because BTU, hours, feet and Fahrenheit are a self-consistent set. Sum the surfaces, add the infiltration and ventilation load, and you have the design heat loss the boiler or heat pump has to meet.

Three traps, in the order they bite. Use gross area for each assembly type and take the windows and doors out of it, or you will pay for the same square metres twice. ΔT is a difference, so 40 °F of it is 22.2 K, not 4.4 — this page converts correctly, but the arithmetic in your notebook may not. And the R you divide by should be the effective R of the assembly including framing, not the number on the batt: a wall built with R-19 batts performs at about R-13 once the studs are counted, so a heat loss calculated off the label is roughly 30% optimistic. That single substitution is the most common reason a load calculation comes in under the building's measured fuel use.

Worked example: 30 m2 of RSI 3.5 wall at 22 K → 188.6 W

Heat Flux Through Insulation (q = ΔT/R)

q=ΔTRq'' = \frac{\Delta T}{R}
ΔTq″R
Where
  • qq''= Heat flux (W/m²)
  • ΔT\Delta T= Temperature difference ()
  • RR= R-value of the assembly (RSI (m²·K/W))

Strip the area out of the heat-loss equation and what is left is flux: watts per square metre, which is the fairest way to compare two assemblies because it does not care how big the building is. RSI 3.5 with 22 K across it passes 22/3.5 = 6.3 W/m². An R-19 wall with 40 °F across it passes 40/19 = 2.1 BTU/(h·ft²), the same 6.6 W/m². Codes and mechanical insulation specifications are often written directly in these terms — a personnel-protection or condensation-control spec is really a flux limit in disguise.

The reason this page matters in the field is the heat flux meter. Tape one to the inside of a wall, log the surface-to-surface or air-to-air ΔT alongside it, and R = ΔT/q″ hands you the as-built resistance of the wall in front of you — including the framing, the settled insulation, the gap the electrician left, and the moisture. ISO 9869 covers the method and asks for at least 72 hours of averaging, because the wall stores heat and a short test measures the weather rather than the wall. Two traps: mount the meter away from studs unless the stud is what you are hunting, and use air-to-air ΔT only if you intend the surface films to be part of the R you report.

Worked example: 22 K across RSI 3.5 → 6.286 W/m2

Effective R-Value with Framing (Parallel Path)

1Reff=ffrRfr+1ffrRcav\frac{1}{R_{eff}} = \frac{f_{fr}}{R_{fr}} + \frac{1 - f_{fr}}{R_{cav}}
RcavRfrffrReff
Where
  • ReffR_{eff}= Effective assembly R-value (RSI (m²·K/W))
  • RcavR_{cav}= R-value of the cavity path (RSI (m²·K/W))
  • RfrR_{fr}= R-value of the framing path (RSI (m²·K/W))
  • ffrf_{fr}= Framing factor

The R-value on a batt is measured at the centre of the cavity, in a guarded hot box, with nothing in the way. The wall you actually build has a stud every 16 inches, and the stud is a parallel path: heat picks whichever route is easier, so the two paths are added as conductances, area-weighted, not as R-values. That is this equation, the parallel-path or isothermal-planes method of ASHRAE Fundamentals. A 2×6 wall with R-19 batts, 6.875 of framing R through 5.5 inches of softwood and a 25% framing factor comes out at 1/(0.25/6.875 + 0.75/19) = R-13.2 — a wall sold as R-19 delivering R-13, and every watt of the difference is real money on a fuel bill.

Framing factors are higher than anyone guesses. Counting studs alone suggests 10–15%, but once you add plates, headers, sills, rim joists, corners and the extra studs around every opening, a conventionally framed wall lands at 23–27% at 16 in on centre, and only careful advanced framing at 24 in on centre gets you near 15%. Two consequences follow. First, you cannot fix a bridged wall by adding cavity insulation — push RcavR_{\mathrm{cav}} to infinity in the equation above and ReffR_{\mathrm{eff}} still cannot exceed Rfr/ffrR_{\mathrm{fr}}/f_{\mathrm{fr}}, which for this wall is 27.5; the studs are a hard ceiling. Second, continuous exterior insulation is added in series with the whole assembly, so it works on both paths at once, and this is exactly why modern energy codes stopped asking for thicker batts and started demanding continuous insulation instead. Steel framing makes the point brutally: steel conducts about 400 times better than wood, so a steel-stud wall filled with R-19 can measure R-7, and the correction there needs the code's tabulated factors rather than a simple two-path split.

Worked example: RSI 3.52 cavity, RSI 1.10 studs, RSI 2.75 whole wall → 12.7% framing

Conduction Through a Pipe Wall

Q˙=2πkLΔTln(r2/r1)\dot{Q} = \frac{2 \pi k L \, \Delta T}{\ln(r_2 / r_1)}
r1r2kLQΔT
Where
  • Q˙\dot{Q}= Heat flow rate (W)
  • kk= Thermal conductivity (W/(m·K))
  • LL= Pipe length (m)
  • ΔT\Delta T= Inner-to-outer ΔT ()
  • r1r_1= Inner radius (mm)
  • r2r_2= Outer radius (mm)

A pipe wall is not a flat slab: heat spreading outward keeps finding more area, so the resistance per unit thickness falls as you go, and the geometry hands you a logarithm instead of a simple L/kA. Ten metres of pipe lagged with 50 mm of k = 0.05 W/(m·K) insulation, from r₁ = 50 mm to r₂ = 100 mm with 100 K across it, loses 2π × 0.05 × 10 × 100 ÷ ln(2) = 453 W. Double the insulation thickness again, to r₂ = 150 mm, and the loss only falls to 279 W — the logarithm is a law of diminishing returns, and it is why insulation schedules stop where they do.

The trap is layering. Each layer needs its own ln(r₂/r₁) with its own k, and its own ΔT; you cannot average the conductivities across a jacketed system. In steam service the second trap is condensate: a 100 mm line at 180 °C bare loses roughly ten times what the lagged line does, and every watt of that loss is steam condensing somewhere it was not meant to, filling traps and hammering elbows. Run this page backwards on a measured surface temperature and you get the k your insulation is actually delivering, which after a decade of rain and mechanical damage is rarely the k in the catalogue.

Worked example: 10 m pipe, 50 → 100 mm lagging, k=0.05, 100 K → 453.2 W

Critical Radius of Insulation

rcr=khr_{cr} = \frac{k}{h}
rcrkh
Where
  • rcrr_{cr}= Critical radius (mm)
  • kk= Insulation conductivity (W/(m·K))
  • hh= Outside film coefficient (W/(m²·K))

Wrapping a cylinder in insulation does two opposite things: it adds conduction resistance, which cuts the loss, and it enlarges the outer surface, which raises the convective loss. Differentiate the total resistance and the two effects balance exactly at rcr=k/hr_{\mathrm{cr}} = k/h. Below that radius, the first millimetres of lagging make the loss worse. For lagging at k = 0.05 W/(m·K) in still air at h = 10 W/(m²·K), rcrr_{\mathrm{cr}} = 5 mm — so any pipe larger than a 10 mm-diameter tube is already past the peak and insulation only helps.

The number matters far more in electrical work than in piping. Wire insulation has k ≈ 0.15 W/(m·K) and sits in near-still air at h ≈ 8, giving a critical radius near 19 mm — larger than most conductors, which means the plastic jacket on a small cable genuinely helps it run cooler while doubling as insulation. Deliberately exploiting this is standard practice for fine thermocouple leads and small transistors. The trap is applying the cylindrical result to a flat wall or a sphere: a plane wall has no critical thickness at all, and for a sphere the answer is 2k/h.

Worked example: k=0.05, h=10 → critical radius 5 mm

Convection, Radiation & Heat Exchangers

Newton's Law of Cooling (Q = hAΔT)

Q˙=hAΔT\dot{Q} = h A \, \Delta T
AQΔTh
Where
  • Q˙\dot{Q}= Heat transfer rate (W)
  • hh= Convection coefficient (W/(m²·K))
  • AA= Surface area ()
  • ΔT\Delta T= Surface-to-fluid ΔT ()

Newton published this in 1701, anonymously and in Latin, as a throwaway note on how a red-hot iron bar cools: the heat leaving a surface is proportional to how far that surface is from the fluid around it. Everything hard about convection is hidden in h, the film coefficient, which is not a material property at all but a shorthand for the whole boundary layer — geometry, velocity, viscosity, whether the fluid is boiling. Still air gives h ≈ 5–25 W/(m²·K); a fan raises it to 25–250; water in a tube runs 500–10,000; and boiling or condensing water can exceed 50,000.

Worked example: a 2.5 m² transformer tank sitting 40 K above ambient with h = 25 W/(m²·K) sheds 25 × 2.5 × 40 = 2500 W. Run it backwards and it becomes the field diagnostic every service technician uses — measure the duty and the surface ΔT, and the h you compute tells you whether the airflow is what the nameplate assumed. The classic trap is using the mean fluid temperature where the correlation wanted the film temperature, or forgetting that a fouled, painted or dusty surface has quietly halved its h since commissioning day.

Worked example: 25 W/(m2.K) over 2.5 m2 at 40 K → 2500 W

Convection Film Resistance

R=1hAR = \frac{1}{h A}
hAR
Where
  • RR= Film resistance (K/W)
  • hh= Convection coefficient (W/(m²·K))
  • AA= Wetted area ()

A moving fluid does not touch a wall at the wall's temperature; it drags a thin, nearly stationary layer along with it, and all of the temperature drop happens inside that film. Treating it as a resistance, R = 1/(hA), lets you drop convection straight into the same series network as the conduction terms. A 2 m² panel in still air at h = 25 W/(m²·K) has a film resistance of 1/50 = 0.02 K/W — larger than 100 mm of concrete behind it, which is the whole reason surface films appear in every building-envelope calculation as fixed air-film allowances.

The trap is assuming the film is small enough to ignore because the coefficient is a big number. Compare resistances, not coefficients: on an air-cooled condenser the air film usually holds 80–95% of the total, so polishing the refrigerant side is wasted money while the fan speed is not. The other trap is area — in a finned or tubular geometry h and A must be quoted on the same surface, or the resistance you compute belongs to a different exchanger than the one on the drawing.

Worked example: h 25 W/(m2.K) over 2 m2 → 0.02 K/W film resistance

Reynolds Number

Re=ρvDμRe = \frac{\rho v D}{\mu}
vρμDRe
Where
  • ReRe= Reynolds number
  • ρ\rho= Fluid density (kg/m³)
  • vv= Flow velocity (m/s)
  • DD= Characteristic length (mm)
  • μ\mu= Dynamic viscosity (Pa·s)

Osborne Reynolds injected dye into pipe flow in 1883 and watched it either glide in a smooth filament or erupt into eddies — and found one dimensionless group predicted which. Below about Re = 2300 pipe flow is laminar; above roughly 4000 it is turbulent. Water at 1 m/s in a 5 cm pipe gives Re = 1000 × 1 × 0.05 / 0.001 = 50 000: solidly turbulent, like nearly all industrial water flow.

Because only the combination ρvD/μ matters, a small model in a wind tunnel can faithfully stand in for a full-size aircraft as long as the Reynolds numbers match — the principle that makes scale testing legitimate.

Worked example: Water, 1 m/s, D = 5 cm, mu = 1 mPa*s → Re = 50000

Prandtl Number

Pr=μcpk\mathrm{Pr} = \frac{\mu c_p}{k}
Where
  • Pr\mathrm{Pr}= Prandtl number
  • μ\mu= Dynamic viscosity (Pa·s)
  • cpc_p= Specific heat (J/(kg·K))
  • kk= Fluid thermal conductivity (W/(m·K))

Ludwig Prandtl's boundary-layer paper of 1904 was eight pages long and reorganised fluid mechanics; the group that carries his name asks which boundary layer is thicker, the velocity one or the thermal one. Pr = 1 means they grow together. Gases cluster tightly near 0.7 — air is 0.707 at room temperature — because momentum and heat are carried by the same wandering molecules. Water is about 7 at 20 °C and falls to 1.75 at 100 °C. Engine oil can exceed 10,000, so its thermal layer is a sliver inside a very thick velocity layer. Liquid metals sit at 0.004–0.03, heat sprinting far ahead of momentum, which is why sodium-cooled reactors need their own correlations entirely.

Worked example: water at 20 °C with μ = 1.002 × 10⁻³ Pa·s, cₚ = 4182 J/(kg·K) and k = 0.598 W/(m·K) gives Pr = 7.01. Because Pr is a pure property, the trap is temperature: water's viscosity halves between 20 °C and 55 °C, so a Prandtl number picked off the wrong row of the table poisons every correlation downstream. Evaluate properties at the film temperature, the average of wall and bulk, unless the correlation you are using explicitly says otherwise.

Worked example: Water at 20 C → Pr 7.007

Nusselt Number

Nu=hLk\mathrm{Nu} = \frac{h L}{k}
khLNu
Where
  • Nu\mathrm{Nu}= Nusselt number
  • hh= Convection coefficient (W/(m²·K))
  • LL= Characteristic length (m)
  • kk= Fluid thermal conductivity (W/(m·K))

Wilhelm Nusselt's 1915 paper on the similarity theory of heat transfer, and his 1916 analysis of laminar film condensation — still the standard result a century later — gave engineering its habit of expressing convection dimensionlessly. Nu is the convective coefficient measured against the conduction that would occur through a stagnant fluid layer of the same thickness: Nu = 1 means the fluid may as well be motionless, Nu = 100 means convection is doing a hundred times better. Nearly every convection correlation ever published has the form Nu = f(Re, Pr), and this equation is how you cash one in.

The workflow is always the same. Compute Re and Pr from the fluid and the flow, look up or apply a correlation to get Nu, then convert to a physical h using h = Nu·k/L. A Nusselt number of 20.8 in water (k = 0.6 W/(m·K)) inside a 50 mm tube means h = 20.8 × 0.6/0.05 = 250 W/(m²·K). The trap is L. It is whatever the correlation's author used — internal-flow correlations use the tube inside diameter, flat-plate correlations the distance from the leading edge, and non-circular ducts the hydraulic diameter 4A/P — and it must be the same L on both sides of the calculation. The k is the fluid's, never the wall's; using the metal's conductivity here is the most common error in the whole subject.

Worked example: h 250 in a 50 mm tube of water (k 0.6) → Nu 20.83

Grashof Number

Gr=gβΔTL3ν2\mathrm{Gr} = \frac{g \, \beta \, \Delta T \, L^{3}}{\nu^{2}}
νΔTgL
Where
  • Gr\mathrm{Gr}= Grashof number
  • gg= Gravitational acceleration (m/s²)
  • β\beta= Volumetric expansion coefficient (1/K)
  • ΔT\Delta T= Surface-to-fluid temperature difference ()
  • LL= Characteristic length (m)
  • ν\nu= Kinematic viscosity (mm²/s)

Stand near a radiator on a cold morning and you can feel air rising off it. Nobody is blowing that air; the radiator warms the layer touching it, warm air is less dense than cool air, and the difference in density is enough to lift it. That is natural convection, and it is the form of heat transfer most of us meet first — a hot mug, a window in the sun, a person in a still room — and have the least vocabulary for. The Grashof number is the vocabulary. It asks whether the buoyancy a temperature difference creates is strong enough to overcome the fluid's own stickiness, Gr=gβΔTL3/ν2\mathrm{Gr} = g\beta\,\Delta T\,L^{3}/\nu^{2}, and it plays exactly the role in natural convection that the Reynolds number plays in forced convection. Reynolds compares the inertia a pump supplies against viscosity; Grashof compares the buoyancy a temperature difference supplies against the same viscosity. Where Reynolds has a fan, Grashof has gravity and a density difference doing the work for free.

Read the pieces and the physics is legible. On top, gβΔTg\beta\Delta T is the buoyant acceleration: how much lighter the heated fluid has become, times gravity. The length appears CUBED because buoyancy acts on a volume while the viscous drag opposing it acts on a face — so a tall wall is not a little more convective than a short one, it is dramatically more. Underneath, ν2\nu^{2} is the viscosity resisting twice over, once in setting up the motion and once in damping it. Note the presence of gg, which no forced-convection group contains: turn gravity off and a hot plate in a spacecraft simply sits in a growing pocket of its own warm air, which is why cooling in orbit needs fans that a laboratory bench does not.

The characteristic length is a convention, not a measurement. This is the point at which most Grashof calculations go wrong, and the reason to be blunt about it. For a vertical plate, LL is the HEIGHT. For a horizontal plate it is the area divided by the perimeter. For a cylinder or a sphere it is the DIAMETER. These are not different measurements of the same thing; they are different agreements, adopted because the correlations that use them were fitted that way. The same wall in the same air has different Grashof numbers under different conventions, and because LL is cubed the difference is orders of magnitude — a 0.5 m wall taken on its height gives roughly a thousand times the Grashof number of the same wall taken on a 0.05 m thickness. A Grashof number without its length convention stated is not a number, it is a rumour. Before feeding one into a correlation, find the sentence in the source that says what LL is.

β is 1/T for an ideal gas and measured for everything else, and the T in that substitution is ABSOLUTE. For air at a film temperature of 320 K, β=1/320=0.003125\beta = 1/320 = 0.003125 per kelvin, exactly — no table needed. Use Celsius there and the answer is nonsense. For liquids there is no shortcut at all, and water is the reason it is worth belabouring: its expansion coefficient swings by more than a factor of ten between 0 and 100 °C, and it passes through ZERO at about 4 °C. Above that temperature water expands as it warms, as everything else does; below it, water CONTRACTS as it warms. So the densest water in a lake sits at 4 °C at the bottom, colder water floats above it, and ice forms at the surface rather than at the bed. That is why a lake freezes from the top down, why fish survive the winter, and — closer to this page — why a natural-convection calculation in cold water can return a Grashof number of nearly nothing and be perfectly correct.

Two practical habits. Evaluate all the properties at the FILM temperature, the average of the surface and the bulk fluid, because the boundary layer is where everything in this problem happens and it is neither at wall temperature nor at room temperature. And do not look for a Grashof threshold: unlike Reynolds, this group does not decide the flow regime on its own. The laminar-to-turbulent judgement in natural convection is made on the Rayleigh number, GrPr\mathrm{Gr}\cdot\mathrm{Pr}, which is where this calculation should go next.

Worked example: 0.5 m wall at 335 K in air at 305 K → Gr = 3.71×10⁸

Rayleigh Number

Ra=GrPr\mathrm{Ra} = \mathrm{Gr} \, \mathrm{Pr}
ΔTναL
Where
  • Ra\mathrm{Ra}= Rayleigh number
  • Gr\mathrm{Gr}= Grashof number
  • Pr\mathrm{Pr}= Prandtl number

The Rayleigh number is a product, Ra=GrPr\mathrm{Ra} = \mathrm{Gr}\,\mathrm{Pr}, and both halves of that sentence deserve saying. It is a product EXACTLY — nothing is measured to obtain it that was not already measured to obtain the Grashof and Prandtl numbers, and multiplying the two out gives gβΔTL3/(να)g\beta\Delta T L^{3}/(\nu\alpha), buoyancy divided by the product of the two diffusivities that damp it. So it is not an independent quantity. And yet it has its own name, its own thresholds and its own correlations, which would be strange if it were only shorthand.

The reason is a genuine experimental finding rather than a definition. When natural-convection data are plotted against Grashof number, fluids scatter: air, water and oil with the same Grashof number transfer heat at visibly different rates. Plot the same data against the PRODUCT and they collapse onto one curve. That collapse is a physical statement — it says buoyancy-driven transfer depends on the combination να\nu\alpha and not on the two diffusivities separately, and that a viscous fluid which also conducts well behaves like a thin fluid which conducts poorly. Once that was established there was no further reason to quote Grashof thresholds, and the literature stopped. Every natural-convection correlation you will meet — Churchill and Chu for plates and cylinders, the enclosure correlations, the horizontal-plate rules — is written in Rayleigh.

Two thresholds get quoted, and they are different kinds of thing. The first is Ra1708\mathrm{Ra} \approx 1708, below which a fluid layer heated from below simply sits there and conducts, and above which it breaks into the hexagonal Rayleigh–Bénard cells you can watch in a pan of oil on a low burner. That figure is a genuine stability result, derived exactly for an idealised infinite layer between rigid plates, and it does not transfer to any other geometry without being re-derived. The second is Ra109\mathrm{Ra} \approx 10^{9} for the laminar-to-turbulent transition on a vertical plate, and that one is a CONVENTION. The real transition wanders roughly between 10810^{8} and 101010^{10} depending on surface roughness, on how the flow was disturbed at the leading edge, and on how quiet the surrounding room is. Treat it as the centre of a range. When a calculation lands near it, the honest move is to run the correlation on both sides and see how much the answer actually moves — usually far less than the argument about which side you are on would suggest.

A working note on scale, because the numbers are unintuitive at first. A coffee cup gives a Rayleigh number around 10610^{6}, a room wall around 10910^{9} to 101010^{10}, and the Earth's mantle something near 102010^{20} — which is why the mantle convects despite being rock. And the Nusselt numbers that come out grow slowly, roughly as Ra1/4\mathrm{Ra}^{1/4} in the laminar range and Ra1/3\mathrm{Ra}^{1/3} in the turbulent one, so a hundredfold increase in Rayleigh buys only a threefold or fivefold increase in the coefficient. Natural convection is reliable and it is cheap; it is not, and cannot be made, powerful.

Worked example: Ra = 10⁹ in water at Pr = 7.0 → Gr = 1.43×10⁸

Dittus-Boelter Correlation

Nu=0.023Re0.8Prn\mathrm{Nu} = 0.023 \, \mathrm{Re}^{0.8} \, \mathrm{Pr}^{n}
Where
  • Nu\mathrm{Nu}= Nusselt number
  • Re\mathrm{Re}= Reynolds number
  • Pr\mathrm{Pr}= Prandtl number
  • nn= Prandtl exponent

Published by F. W. Dittus and L. M. K. Boelter at Berkeley in 1930 — in a university engineering bulletin, not a journal — this is the most-used correlation in heat transfer, and it is deliberately crude. It fits fully developed turbulent flow inside a smooth circular tube, with the stated validity range: Re > 10,000, 0.6 < Pr < 160, and length-to-diameter ratio above about 10 so the entrance region has stopped mattering. Properties are taken at the bulk mean temperature, and the exponent n switches between 0.4 when the fluid is being heated and 0.3 when it is being cooled, which crudely accounts for the way viscosity near the wall changes the velocity profile.

Expect ±25% scatter against experiment, and worse if the wall-to-bulk temperature difference is large or the fluid is very viscous — for those cases Sieder-Tate adds a (μ/μwall)0.14(\mu/\mu_{\mathrm{wall}})^{0.14} factor, and Gnielinski's 1976 correlation does far better across the transition region. Worked example: Re = 50,000, Pr = 4.5, heating, gives Nu=0.023×50,0000.8×4.50.4=0.023×5743×1.825=241\mathrm{Nu} = 0.023 \times 50{,}000^{0.8} \times 4.5^{0.4} = 0.023 \times 5743 \times 1.825 = 241, which in water inside a 25 mm tube is h ≈ 5800 W/(m²·K). The trap worth remembering is the exponent on Re: 0.8 means doubling the velocity buys only 74% more coefficient while the pressure drop rises about fourfold, so there is always a point beyond which pumping the tubes harder is a losing trade.

Worked example: Re 50,000, Pr 4.5, heating → Nu 241.1

Stefan-Boltzmann Law

P=εσAT4P = \varepsilon \sigma A T^4
TPAε
Where
  • PP= Radiated power (W)
  • ε\varepsilon= Emissivity
  • AA= Surface area ()
  • TT= Surface temperature (°C)

Every surface warmer than absolute zero radiates, and the power it sends out climbs with the fourth power of its absolute temperature. Double the kelvin temperature and the radiated power grows sixteenfold. The fourth power is not arbitrary: the number of photons a hot body emits per second scales roughly as T3T^3, and the average energy each one carries scales as TT, so the product goes as T4T^4. The Stefan–Boltzmann constant σ = 5.670374419 × 10⁻⁸ W/(m²·K⁴) is exact in the 2019 SI, and emissivity ε runs from 0 to 1, where 1 is a perfect black body.

Work a person. A clothed adult presents about 1.8 m² of surface at roughly 28 °C (301 K) with an emissivity near 0.98. The gross emission is 0.98×5.67×108×1.8×3014823 W0.98 \times 5.67 \times 10^{-8} \times 1.8 \times 301^4 \approx 823\ \text{W} — which is absurd, since no one eats 823 W. The resolution is in the next paragraph but one, and it is the whole practical lesson of this equation.

Josef Stefan found the fourth-power rule empirically in 1879, fitting it to John Tyndall's measurements of glowing platinum wire. Ludwig Boltzmann derived it from thermodynamics five years later, treating radiation as a gas that exerts pressure, and the joint name has stuck since. It sat as an empirical law with a thermodynamic argument behind it until 1900, when Planck's radiation formula produced it by integration over all wavelengths — one of the first things anyone checked about the new quantum picture. On the neighbouring astronomical page it does its most famous work: a star's luminosity is 4πR2σT44\pi R^2 \sigma T^4, so measuring brightness and surface temperature gives the radius of an object no telescope can resolve.

Nothing radiates into a void, and forgetting that is the standard error. That person is standing in a room whose walls are at 20 °C, and those walls are radiating back. What you feel is the net exchange, εσA(T4Tsurr4)\varepsilon\sigma A(T^4 - T_{surr}^4), which here comes to about 84 W — a believable figure and one-tenth of the gross. Use the bare form when the surroundings really are cold, as with a spacecraft radiator or a clear night sky, and use the difference form for anything in a room. It also explains why a 20 °C room feels cold beside a single-glazed window and comfortable beside an insulated wall at the same air temperature: the window's inner surface is colder, so you lose more by radiation while the thermometer reports nothing amiss.

Two further traps. Because of the fourth power, a Celsius temperature is not merely inaccurate, it is destroyed: entering 100 instead of 373.15 K is wrong by a factor of (373.15/100)4194(373.15/100)^4 \approx 194. And emissivity is a property of the surface, not the substance, and not of its colour to the eye. Polished aluminium sits near 0.05, the same aluminium oxidised near 0.2, and almost every paint — white, black, or anything between — sits around 0.9 in the thermal infrared, because visible colour says nothing about behaviour at 10 μm. This is why an infrared thermometer pointed at shiny bare pipe reads far too low: the instrument assumes an emissivity, usually 0.95, and the metal is not obliging. A strip of matte tape on the pipe fixes the reading.

Worked example: 1 m² blackbody at 1000 K radiates 56703.7 W

Net Radiation Exchange Between Surfaces

Q˙=εσA(T14T24)\dot{Q} = \varepsilon \sigma A (T_1^4 - T_2^4)
T1T2QεA
Where
  • Q˙\dot{Q}= Net radiant heat rate (W)
  • ε\varepsilon= Emissivity
  • AA= Surface area ()
  • T1T_1= Surface temperature (°C)
  • T2T_2= Surroundings temperature (°C)

Every surface both emits and absorbs radiation, and the net exchange with large surroundings goes as the difference of fourth powers, not as a simple ΔT. That non-linearity is why radiation is negligible in a chilled-water pipe and dominant in a furnace. A 1.2 m² oxidised steel panel (ε = 0.85) at 500 K facing a 300 K room radiates 0.85 × 5.670 × 10⁻⁸ × 1.2 × (6.25 × 10¹⁰ − 8.1 × 10⁹) = 3146 W — with a convective coefficient of 10 W/(m²·K) it would shed only about 2400 W by convection, so more than half the heat leaves as light you cannot see.

Emissivity is where field work goes wrong. Polished aluminium is 0.04, mill-finish steel 0.2–0.3, the same steel after a summer outdoors 0.7–0.85, and almost every paint, oxide, brick and organic surface is 0.85–0.95 regardless of colour — white paint is as good a radiator in the infrared as black. That is the trap in infrared thermography: point a camera set for ε = 0.95 at a shiny bus bar and it reports the temperature of whatever the bar is reflecting, which is usually you. Tape a square of matte tape on the target and read that instead. The formula assumes the surroundings are large enough to behave as a black enclosure; two comparable surfaces facing each other need view factors and a full radiosity network.

Worked example: eps 0.85, 1.2 m2, 500 K to 300 K → 3146 W

Combined Convection and Radiation Coefficient

ht=hc+εσ(Ts+Tsur)(Ts2+Tsur2)h_t = h_c + \varepsilon \sigma (T_s + T_{sur})(T_s^2 + T_{sur}^2)
TshcεhtTsur
Where
  • hth_t= Combined coefficient (W/(m²·K))
  • hch_c= Convection coefficient (W/(m²·K))
  • ε\varepsilon= Surface emissivity
  • TsT_s= Surface temperature (°C)
  • TsurT_{sur}= Surroundings temperature (°C)

Radiation is a fourth-power law, but T₁⁴ − T₂⁴ factors exactly into (T₁ − T₂)(T₁ + T₂)(T₁² + T₂²), so over a modest temperature range you can hide everything but the linear ΔT inside an equivalent coefficient hr=εσ(Ts+Tsur)(Ts2+Tsur2)h_r = \varepsilon\sigma(T_s + T_{\mathrm{sur}})(T_s^{2} + T_{\mathrm{sur}}^{2}). Add it to the convective film and a single hth_t drives the whole surface with plain old Q̇ = htAΔTh_t A \Delta T — which is exactly what makes building-envelope and insulation software tractable.

The size of the radiation term surprises people. A painted surface (ε = 0.9) at 350 K facing 293 K surroundings has hrh_r = 6.8 W/(m²·K), larger than the 5 W/(m²·K) of still-air natural convection beside it, giving hth_t = 11.8. That is why a bare hot pipe in a still basement loses more than half its heat by radiation, and why a low-emissivity foil wrap — ε ≈ 0.05 — kills that channel almost completely while doing nothing about convection. The two traps: hrh_r depends on both temperatures, so it is not a constant and must be recomputed if the surface moves far; and the surroundings temperature is the temperature of the walls seeing the surface, not the air temperature, which on a clear night can be 20 K colder than the air and is why cars frost over at 4 °C.

Worked example: hc 5, eps 0.9, 350 K surface in a 293 K room → ht 11.84

Overall Heat Transfer Coefficient (U)

1U=1hi+Lk+1ho\frac{1}{U} = \frac{1}{h_i} + \frac{L}{k} + \frac{1}{h_o}
khihoUL
Where
  • UU= Overall coefficient (W/(m²·K))
  • hih_i= Inside film coefficient (W/(m²·K))
  • LL= Wall thickness (mm)
  • kk= Wall thermal conductivity (W/(m·K))
  • hoh_o= Outside film coefficient (W/(m²·K))

Resistances in series add, and heat transfer borrows the electrical analogy wholesale: 1/U is the total resistance per unit area, the sum of the inside film, the metal, and the outside film. What falls out is the engineer's most useful insight — the largest resistance owns the answer. Take a steel exchanger tube, 3 mm of steel at k = 45 W/(m·K), water inside at hᵢ = 5000 and air outside at hₒ = 2000: 1/U = 0.000200 + 0.0000667 + 0.000500 = 0.000767, so U = 1304 W/(m²·K). Note that the steel — the only part you can see — contributes 9% of the resistance. Doubling the tube wall barely moves U; doubling the air-side film nearly doubles it.

That is why finned tubes exist. Air-side coefficients are typically 20–50 times worse than water-side ones, so every gas-to-liquid exchanger in the world grows fins on the gas side to buy back area rather than coefficient. The trap on real datasheets is area basis: U must be quoted against a stated area (usually the outside), and a U of 850 W/(m²·K) on inside area is not the same machine as 850 on outside area. When you see U and A on a drawing, check which surface A refers to before you believe the duty.

Worked example: hi 5000, 3 mm steel k=45, ho 2000 → U = 1304.3 W/(m2.K)

Fouled Overall Coefficient

1Uf=1Uc+1hf\frac{1}{U_f} = \frac{1}{U_c} + \frac{1}{h_f}
UchfUf
Where
  • UfU_f= Fouled overall coefficient (W/(m²·K))
  • UcU_c= Clean overall coefficient (W/(m²·K))
  • hfh_f= Fouling conductance (W/(m²·K))

Fouling is the slow accumulation of scale, biofilm, corrosion product and process gunk on a heat transfer surface, and it is the single largest source of lost capacity in operating plants. TEMA tabulates it as a fouling factor RfR_f in m²·K/W (or h·ft²·°F/BTU): 0.00018 for treated cooling-tower water, 0.00035 for river water, 0.0009 for untreated seawater, up to 0.002 for a heavy fuel oil. This page uses its reciprocal, the fouling conductance hf=1/Rfh_f = 1/R_f, so it lives in the same W/(m²·K) picker as the rest of the shard — 0.0002 m²·K/W is hfh_f = 5000 W/(m²·K), and 0.001667 h·ft²·°F/BTU is 600 BTU/(h·ft²·°F).

The arithmetic is brutal in the direction people find surprising. A clean U of 1200 W/(m²·K) with a modest RfR_f of 0.0002 falls to 1/(0.000833 + 0.0002) = 968 — a 19% loss from a deposit you could scrape off with a fingernail. On a high-U plate exchanger the same fouling factor can cost 40%, because the fouling resistance is now comparable to everything else in the stack. This is why shell-and-tube plants foul themselves into shutdowns: designers add fouling allowance as extra surface, the oversized unit runs at lower velocity than intended, low velocity deposits more solids, and the margin that was supposed to protect the exchanger is what killed it. Keep tube-side velocity above about 1 m/s and the chemistry in range, and the allowance stays an allowance.

Worked example: Clean U 1200 with R_f = 0.0002 → fouled U 967.7

Fouling Factor on an Overall Coefficient

1Uf=1Uc+Rf\frac{1}{U_f} = \frac{1}{U_c} + R_f
UcRfUf
Where
  • UfU_f= Fouled (service) coefficient (W/(m²·K))
  • UcU_c= Clean coefficient (W/(m²·K))
  • RfR_f= Fouling factor (RSI (m²·K/W))

A fouling factor is an R-value. Same quantity, same unit, same arithmetic as the insulation on a wall — TEMA simply quotes it as RfR_f in h·ft²·°F/BTU or m²·K/W and adds it to 1/U. Typical TEMA design values: treated cooling tower water 0.0002 m²·K/W (0.001 in imperial), river water 0.0004, seawater 0.0002, fuel oil 0.0005, and steam 0.00009. The insight the number hides is that fouling hurts a good exchanger far more than a bad one. Take a clean U of 2500 W/(m²·K), which is a plate unit on clean duty: 1/2500 = 0.0004, so a fouling factor of 0.0004 doubles the total resistance and halves the coefficient to 1250. Apply that same 0.0004 to an air-cooled unit at U = 50 and it costs you 2%.

That is why plate exchangers are specified with small fouling allowances and cleaned in place instead, while shell-and-tube units carry generous ones — and it is where the specification trap lives. Fouling allowance is bought as extra surface, and extra surface on a water-cooled unit means lower velocity in the tubes, which fouls faster. Over-specifying RfR_f is self-fulfilling; 20–30% excess surface is common and defensible, 100% is a fouling machine. The field version of this equation runs backwards: measure the duty and the LMTD, back out the service U, and Rf=1/Uf1/UcR_f = 1/U_f - 1/U_c is the deposit you have accumulated since commissioning. When that number reaches the design allowance, the exchanger is due for cleaning — that is what the design fouling factor was always for, and it is a maintenance trigger, not a safety factor.

Worked example: Clean U 2500 with R_f = 0.0004 → U halves to 1250

Log Mean Temperature Difference (Counterflow)

ΔTlm=ΔT1ΔT2ln(ΔT1/ΔT2)\Delta T_{lm} = \frac{\Delta T_1 - \Delta T_2}{\ln(\Delta T_1 / \Delta T_2)}
Th,inTh,outTc,outTc,inΔTlm
Where
  • ΔTlm\Delta T_{lm}= Log mean temperature difference ()
  • Th,inT_{h,in}= Hot stream inlet (°C)
  • Th,outT_{h,out}= Hot stream outlet (°C)
  • Tc,inT_{c,in}= Cold stream inlet (°C)
  • Tc,outT_{c,out}= Cold stream outlet (°C)

The driving ΔT in an exchanger is not constant along its length, so you cannot use the arithmetic mean — integrating Q̇ = UAΔT along the tube produces the logarithmic mean of the two terminal differences instead. In counterflow the streams run opposite ways, so the hot end pairs the hot inlet with the cold outlet: ΔT₁ = Th,in − Tc,out and ΔT₂ = Th,out − Tc,in. Cool 150 °C oil to 90 °C against water warming from 30 °C to 70 °C and the terminals are 80 K and 60 K, giving ΔTlm=(8060)/ln(80/60)=69.5\Delta T_{\mathrm{lm}} = (80 - 60)/\ln(80/60) = 69.5 K, not the 70 K an average would suggest. The log mean is always the smaller of the two, and the gap widens fast as the terminals diverge.

Counterflow is the reason so many exchangers are plumbed the way they are: it permits a temperature cross, where the cold stream leaves hotter than the hot stream leaves, which parallel flow can never do. Two traps. First, a genuine cross at either terminal makes the logarithm undefined — the calculator refuses, because the arrangement you described cannot exist. Second, this page solves for ΔTlm\Delta T_{\mathrm{lm}} only: the four terminal temperatures sit inside a logarithm and a difference at once, so recovering an inlet temperature from a known ΔTlm\Delta T_{\mathrm{lm}} is transcendental and belongs to an iterative solver, not a closed form.

Worked example: 150 → 90 C against 30 → 70 C counterflow → 69.52 K

Log Mean Temperature Difference (Parallel Flow)

ΔTlm=ΔT1ΔT2ln(ΔT1/ΔT2)\Delta T_{lm} = \frac{\Delta T_1 - \Delta T_2}{\ln(\Delta T_1 / \Delta T_2)}
Th,inTh,outTc,inTc,outΔTlm
Where
  • ΔTlm\Delta T_{lm}= Log mean temperature difference ()
  • Th,inT_{h,in}= Hot stream inlet (°C)
  • Th,outT_{h,out}= Hot stream outlet (°C)
  • Tc,inT_{c,in}= Cold stream inlet (°C)
  • Tc,outT_{c,out}= Cold stream outlet (°C)

In parallel (co-current) flow both fluids enter at the same end, so the pairing changes: ΔT₁ = Th,in − Tc,in at the inlet end and ΔT₂ = Th,out − Tc,out at the outlet end. The same 150 °C oil and 30 °C water, delivered to 90 °C and 70 °C, now give terminals of 120 K and 20 K and a log mean of only 55.8 K — a fifth less driving force than the counterflow arrangement, from identical fluids at identical temperatures. Same duty, same U, and you need 25% more surface. That is why counterflow is the default and parallel flow needs a reason.

It does have reasons. Parallel flow puts the biggest ΔT where the cold fluid is coldest, which brings a viscous fluid up to temperature fast, and it holds the hot-end wall temperature lower, which matters when a product scorches, a coating cures or a thermally sensitive fluid must never see a hot tube. The hard limit is thermodynamic: the two outlet temperatures can approach each other but can never cross, so a parallel-flow unit can never heat the cold stream above the hot stream's exit. If your process needs a cross, no amount of surface in a co-current unit will deliver it.

Worked example: Same streams in parallel flow → 55.81 K (vs 69.52 counter)

Heat Exchanger Duty (Q = U·A·F·LMTD)

Q˙=UAFΔTlm\dot{Q} = U A F \, \Delta T_{lm}
QUAFΔTlm
Where
  • Q˙\dot{Q}= Exchanger duty (kW)
  • UU= Overall coefficient (W/(m²·K))
  • AA= Heat transfer area ()
  • FF= LMTD correction factor
  • ΔTlm\Delta T_{lm}= Log mean temperature difference ()

This is the equation on which the world's exchangers are bought and sold. Compute the counterflow log mean, multiply by U, A and the correction factor F, and you have the duty. F answers a single question: how much worse than pure counterflow is this geometry? A 1-2 shell-and-tube unit — one shell pass, two tube passes — has half its tubes running the wrong way, so F falls below 1; crossflow coils with one or both fluids unmixed sit somewhere between. F comes from charts plotted against the parameters P = (Tc,out − Tc,in)/(Th,in − Tc,in) and R = (Th,in − Th,out)/(Tc,out − Tc,in), and it is a factor, never a bonus: F ≤ 1 always, and true counterflow is F = 1.

The design rule handed down since Bowman, Mueller and Nagle published the F charts in 1940 is: never design below F = 0.80. Not because the physics fails, but because the chart goes vertical there — a one-degree measurement error in a terminal temperature swings F by a tenth, and your exchanger's duty becomes a guess. Cross that line and the answer is more shells in series, not more tubes. Worked example: U = 850 W/(m²·K) on 24 m² with F = 0.95 and a 30 K log mean gives 850 × 24 × 0.95 × 30 = 581 kW. Run it backwards from the observed duty and the F you compute is a fouling alarm — F does not degrade with time, so if the equation only balances at F = 0.6, the real culprit is U.

Worked example: U 850, A 24 m2, F 0.95, LMTD 30 K → 581.4 kW

Stream Duty from Mass Flow (Q = ṁcΔT)

Q˙=m˙cpΔT\dot{Q} = \dot{m} \, c_p \, \Delta T
cpQΔT
Where
  • Q˙\dot{Q}= Stream duty (kW)
  • m˙\dot{m}= Mass flow rate (kg/h)
  • cpc_p= Specific heat (J/(kg·K))
  • ΔT\Delta T= Temperature change ()

Every exchanger calculation has two halves that must agree. The transfer side says Q̇ = UAF·ΔTlm\Delta T_{\mathrm{lm}}; the process side says Q̇ = ṁcₚΔT for each stream. Write both, set them equal, and the whole problem closes. Because the heat leaving the hot stream must arrive in the cold one, ṁcₚΔT for the hot side equals ṁcₚΔT for the cold side — so the stream with the smaller ṁcₚ, the smaller heat capacity rate, always shows the larger temperature swing. That single observation lets you sanity-check a datasheet from across the room: if both streams change by the same amount, their capacity rates are equal.

Worked example: 2.5 kg/s of water (cₚ = 4186 J/(kg·K)) heated 12 K takes 2.5 × 4186 × 12 = 125.6 kW. In North American units the same physics reads 500,000 BTU/h into 20,000 lb/h of oil at cₚ = 0.5 across 50 °F. Traps: cₚ is not constant — water is flat enough to ignore, but oils and glycols vary 10–20% over a working range, so use the value at the mean temperature. And this equation is sensible heat only; the moment anything boils or condenses, the temperature stops moving and you need ṁ times the latent heat instead.

Worked example: 2.5 kg/s water, cp 4186, 12 K → 125.58 kW

Number of Transfer Units (NTU)

NTU=UAm˙cp\mathrm{NTU} = \frac{U A}{\dot{m} \, c_p}
UAcpNTU
Where
  • NTU\mathrm{NTU}= Number of transfer units
  • UU= Overall coefficient (W/(m²·K))
  • AA= Heat transfer area ()
  • m˙\dot{m}= Minimum stream mass flow (kg/h)
  • cpc_p= Minimum stream specific heat (J/(kg·K))

NTU is an exchanger's size measured in the only currency that matters — conductance compared with the thermal inertia of the fluid flowing through it. UA is how well the machine can transfer; ṁcₚ is how much heat the limiting stream can carry per degree. Their ratio is dimensionless and tells you immediately what class of equipment you are holding: NTU below 0.5 is a trim heater that barely touches the fluid, 1–3 is normal process duty, 3–5 is a close-approach unit, and above 5 you are in regenerator and cryogenic territory where surface becomes very expensive per degree gained.

The effectiveness-NTU method was developed by W. M. Kays and A. L. London for the compact heat exchangers of gas-turbine regenerators and published in Compact Heat Exchangers (1955) — a book still on working desks seventy years later. Their motivation was practical: LMTD design demands all four terminal temperatures, but a regenerator problem usually gives you the two inlets and asks what comes out, which drives LMTD into iteration and NTU straight to an answer. Worked example: U = 500 W/(m²·K) on 8 m² against 1.5 kg/s of water gives NTU = 4000/6279 = 0.64, a small unit. Trap: ṁcₚ must be the minimum stream's, not whichever stream you measured first.

Worked example: UA 4000 W/K against 1.5 kg/s of water → NTU 0.637

Capacity Rate Ratio (Cr)

Cr=m˙mincminm˙maxcmaxC_r = \frac{\dot{m}_{min} c_{min}}{\dot{m}_{max} c_{max}}
mincminmaxcmaxCr
Where
  • CrC_r= Capacity rate ratio
  • m˙min\dot{m}_{min}= Minimum stream mass flow (kg/h)
  • cminc_{min}= Minimum stream specific heat (J/(kg·K))
  • m˙max\dot{m}_{max}= Maximum stream mass flow (kg/h)
  • cmaxc_{max}= Maximum stream specific heat (J/(kg·K))

Heat capacity rate, C = ṁcₚ in watts per kelvin, is how much heat a stream absorbs for each degree it warms. Divide the smaller by the larger and you get Cr, which by construction runs from 0 to 1 and controls how effectiveness responds to size. Example: 1.2 kg/s of air (cₚ ≈ 1005) is 1206 W/K; 0.8 kg/s of water (cₚ ≈ 4186) is 3349 W/K; Cr = 0.36, and the air — despite the higher flow — is the limiting stream, because water carries four times the heat per kilogram per degree.

The two ends of the range are the interesting ones. Cr = 0 means one stream's capacity rate is effectively infinite, which is exactly what happens when a fluid boils or condenses: it absorbs heat at constant temperature, and every exchanger arrangement — counterflow, parallel, crossflow — collapses to the same ε=1eNTU\varepsilon = 1 - e^{-\mathrm{NTU}}. Cr = 1 is the balanced exchanger, hardest to make effective, and the case where counterflow's advantage over parallel flow is largest. The trap is bookkeeping: identify CminC_{\mathrm{min}} from ṁcₚ, not from flow rate alone. Steam-to-water and refrigerant-to-air units are Cr = 0 problems no matter what the flow meters read.

Worked example: Air 1206 W/K over water 3348.8 W/K → Cr 0.3601

Maximum Possible Heat Transfer (Qmax)

Q˙max=m˙mincmin(Th,inTc,in)\dot{Q}_{max} = \dot{m}_{min} c_{min} (T_{h,in} - T_{c,in})
cminTh,inTc,inQmax
Where
  • Q˙max\dot{Q}_{max}= Maximum possible duty (kW)
  • m˙min\dot{m}_{min}= Minimum stream mass flow (kg/h)
  • cminc_{min}= Minimum stream specific heat (J/(kg·K))
  • Th,inT_{h,in}= Hot stream inlet (°C)
  • Tc,inT_{c,in}= Cold stream inlet (°C)

Before asking how well an exchanger performs, you have to know what perfection would look like. Q̇max is that reference: an infinitely long counterflow unit in which the limiting stream is brought all the way to the other stream's inlet temperature. Only the minimum capacity rate may be used — if you used the larger one, the smaller stream would have to overshoot past the other stream's inlet, which is the second law being violated in plain sight.

Worked example: 0.9 kg/s of air (cₚ ≈ 1005 J/(kg·K)) entering at 25 °C against exhaust gas at 200 °C. Q̇max = 0.9 × 1005 × 175 = 158.3 kW. Whatever the real recuperator recovers, it is a fraction of that number, and that fraction is the effectiveness. This is also the quickest audit tool on a plant walkdown: measure the two inlet temperatures and the limiting flow, compute Q̇max, compare it with the duty the process is actually getting, and you have a percentage that tells you whether a cleaning, a re-pass or a new exchanger is the honest recommendation. Trap: the inlet-to-inlet difference, never the inlet-to-outlet difference of one stream.

Worked example: 0.9 kg/s air, 200 C gas against 25 C air → Qmax 158.3 kW

Heat Exchanger Effectiveness (ε = Q/Qmax)

ε=Q˙Q˙max\varepsilon = \frac{\dot{Q}}{\dot{Q}_{max}}
Where
  • ε\varepsilon= Effectiveness
  • Q˙\dot{Q}= Actual duty (kW)
  • Q˙max\dot{Q}_{max}= Maximum possible duty (kW)

Effectiveness is the one exchanger number a non-specialist can read without a chart: 0.82 means the unit captured 82% of everything thermodynamics allowed. Unlike efficiency in a boiler sense it has no fuel in it, and unlike U it needs no area basis to be meaningful. Typical values: a plate heat exchanger with a close approach runs 0.85–0.95, a shell-and-tube process cooler 0.6–0.8, an air-to-air plate recovery core in an HRV 0.6–0.75, a rotary wheel 0.75–0.85, and a fouled unit somewhere well below where it started.

Its real power is diagnostic. Because ε for a given geometry depends only on NTU and Cr, a drop in measured effectiveness at unchanged flows can only mean UA has fallen — which on a water side means scale or biofilm, and on an air side usually means a plugged coil face. Worked example: an economiser recovering 420 kW where the inlet-to-inlet ceiling is 600 kW runs at ε = 0.70; six months later the same flows return 480 kW against a 750 kW ceiling, ε = 0.64, and the cleaning is overdue. The trap is comparing effectiveness values measured at different flow rates — dropping the flow raises NTU and flatters ε, so an exchanger can look better simply because the pump is throttled.

Worked example: 420 kW recovered of a 600 kW ceiling → eps 0.70

Effectiveness from NTU (Counterflow)

ε=1eNTU(1Cr)1CreNTU(1Cr)\varepsilon = \frac{1 - e^{-\mathrm{NTU}(1 - C_r)}}{1 - C_r \, e^{-\mathrm{NTU}(1 - C_r)}}
Where
  • ε\varepsilon= Effectiveness
  • NTU\mathrm{NTU}= Number of transfer units
  • CrC_r= Capacity rate ratio

Effectiveness is the fraction of the thermodynamically possible heat an exchanger actually moves, and for counterflow it depends on nothing but NTU and Cr. The shape of the curve is the practical lesson: at NTU = 1 with Cr = 0.5 you get ε ≈ 0.56; doubling the surface to NTU = 2 buys ε ≈ 0.75; doubling again to NTU = 4 buys only 0.88. Surface is bought at a rising price per point of effectiveness, which is why the argument in every exchanger review is about where on that curve the money stops making sense. Worked example: NTU = 1.5, Cr = 0.5 gives x=e0.75=0.4724x = e^{-0.75} = 0.4724 and ε = 0.5276/0.7638 = 0.691.

Two limits deserve memorising. When Cr → 0 — a boiling or condensing stream — the expression collapses to ε=1eNTU\varepsilon = 1 - e^{-\mathrm{NTU}}, the best any arrangement can do. When Cr = 1, the balanced exchanger, the algebra goes 0/0 and the true limit is ε = NTU/(1 + NTU); this page detects that case and uses the limit rather than returning a NaN. The trap is inverting in your head: given ε and Cr the NTU comes back cleanly, but given ε and NTU there is no closed form for Cr, because it sits in the exponent and the denominator at once — that one needs an iterative solver, so this page does not offer it.

Worked example: NTU 1.5, Cr 0.5 counterflow → effectiveness 0.691

Biot Number

Bi=hLck\mathrm{Bi} = \frac{h L_c}{k}
kLchBi
Where
  • Bi\mathrm{Bi}= Biot number
  • hh= Film coefficient (W/(m²·K))
  • LcL_c= Characteristic length (m)
  • kk= Solid thermal conductivity (W/(m·K))

Named for Jean-Baptiste Biot, whose 1804 experiments on heated bars preceded Fourier's theory by two decades, this is the number that decides whether you may treat a cooling object as a single lump. It compares the resistance to getting heat out of the body (Lc/k)(L_c/k) with the resistance to getting it off the surface (1/h). Below Bi = 0.1 the interior is within a few percent of uniform and the lumped-capacitance method is legitimate; above it, the centre lags the skin badly and you need a chart, a series solution or a finite-element model.

The characteristic length is volume divided by surface area, not the diameter — for a sphere that is r/3 and for a long cylinder r/2, and forgetting the factor is the classic error that pushes a valid problem out of the lumped regime on paper. Worked example: a 20 mm steel billet, LcL_c = 0.01 m, k = 45 W/(m·K), quenched in air at h = 100 gives Bi = 0.022, comfortably lumped. Quench the same billet in agitated water at h = 5000 and Bi = 1.1, so the surface transforms while the core is still glowing — which is precisely the metallurgy that hardening exploits, and precisely why quench cracks happen.

Worked example: Steel billet, h 100, Lc 10 mm, k 45 → Bi 0.0222

Fourier Number

Fo=ktρcL2\mathrm{Fo} = \frac{k \, t}{\rho \, c \, L^{2}}
kρcptLFo
Where
  • Fo\mathrm{Fo}= Fourier number
  • kk= Thermal conductivity (W/(m·K))
  • tt= Elapsed time (s)
  • ρ\rho= Density (kg/m³)
  • cpc_p= Specific heat (J/(kg·K))
  • LL= Characteristic length (m)

The Fourier number is dimensionless time: how far a thermal disturbance has diffused compared with the size of the object. It is normally written Fo = αt/L² with the thermal diffusivity α = k/(ρc), but diffusivity has no entry in this calculator's unit picker, so the group is spelled out as kt/(ρcL²) — identical physics, and it has the pleasant side effect of showing where the diffusivity comes from. High-k, low-ρc materials diffuse heat fast: copper's α is 1.1 × 10⁻⁴ m²/s, steel's 1.2 × 10⁻⁵, brick's 5 × 10⁻⁷, and that thousand-fold spread is why a copper pan responds instantly and a masonry wall takes half a day.

Fo ≈ 1 is the rough marker for "the disturbance has crossed the body". Below Fo = 0.2 the one-term approximations in the textbook charts are not valid and you need the full series; above about 1 the transient is essentially over. Worked example: a 100 mm steel plate (L = 0.05 m half-thickness, k = 45, ρ = 7850, c = 480) after 10 minutes has Fo = 45 × 600/(7850 × 480 × 0.0025) = 2.87, thoroughly soaked through. The same plate in firebrick would need most of a day. This is the number behind cooking times, heat-treat soak schedules and the thermal-mass lag that lets a stone building coast through an afternoon.

Worked example: Steel, 50 mm, 600 s → Fo 2.866 (alpha = 1.19e-5 m2/s)

Lumped Capacitance Time Constant

τ=ρVchA\tau = \frac{\rho V c}{h A}
ρVcphAτ
Where
  • τ\tau= Thermal time constant (s)
  • ρ\rho= Density (kg/m³)
  • VV= Body volume (L)
  • cpc_p= Specific heat (J/(kg·K))
  • hh= Film coefficient (W/(m²·K))
  • AA= Surface area ()

ρVc is the heat a body stores per kelvin — its thermal capacitance — and hA is the conductance draining it. Their ratio is a time constant in exactly the sense an electrical engineer means, and the analogy is complete: the body is a capacitor, the surface film is a resistor, and the temperature decays exponentially. In one τ the gap to ambient closes by 63%, in three τ by 95%, in five τ by over 99%.

Worked example: an aluminium block, ρ = 2700 kg/m³, V = 100 cm³, c = 900 J/(kg·K), cooling in air at h = 30 W/(m²·K) over A = 0.02 m². τ = 2700 × 0.0001 × 900/(30 × 0.02) = 243/0.6 = 405 s, so it is essentially at room temperature in about twenty minutes. This is the number that sizes thermocouple response — a fine bead responds in milliseconds, a 6 mm thermowell in a stagnant pocket can lag a minute, and a control loop tuned without knowing which one you have will hunt forever. Trap: τ assumes the lumped regime, so check Bi < 0.1 before trusting it, and remember h is not constant during a violent transient — free convection off a hot block starts strong and weakens as the block cools.

Worked example: 100 cm3 aluminium block, h 30 over 0.02 m2 → tau 405 s

Lumped Capacitance Cooling Curve

T=T+(T0T)et/τT = T_\infty + (T_0 - T_\infty) e^{-t/\tau}
TTT0τt
Where
  • TT= Temperature at time t (°C)
  • T0T_0= Initial temperature (°C)
  • TT_\infty= Fluid temperature (°C)
  • tt= Elapsed time (s)
  • τ\tau= Thermal time constant (s)

When Bi < 0.1 the interior of a body stays essentially uniform, an energy balance gives ρVc dT/dt = −hA(T − T∞), and the solution is a pure exponential approach to ambient. Every point on the curve is the same fraction of the remaining gap, which is why the answer never depends on how you got there. Worked example: an aluminium block at 200 °C with τ = 405 s, cooling in 25 °C air, after 300 s sits at 25+175×e0.741=108.425 + 175 \times e^{-0.741} = 108.4 °C.

Run it the other way and it becomes the field measurement everyone actually uses: log a cooling curve, read the time to fall from 300 °F to 150 °F in 70 °F air, and τ = 634 ÷ ln(230/80) drops out — no need to know ρ, V, c, h or A at all. Two traps. First, temperature differences are what matter, so the ratio inside the logarithm works in any consistent scale, but the temperatures you type must be absolute or Celsius, not differences. Second, the model dies quietly when h is not constant: a body that starts by boiling its quench fluid and finishes in ordinary convection has two different time constants, and forcing one exponential through that data gives an h that describes neither regime.

Worked example: 200 C block, tau 405 s, 25 C air, after 300 s → 108.4 C

Practice problems

Answer key at the back. Work in the units each problem states.

Properties of State

1. Gauge versus absoluteThe drum gauge on a package boiler reads 550 kPa(g). The plant barometer stands at 99.5 kPa. Calculate the absolute pressure in the drum.

2. Gauge versus absoluteA surface condenser runs under vacuum. Its shell gauge reads -80 kPa(g) with the barometer at 101.3 kPa. Determine the absolute pressure inside the condenser shell.

3. Heating a massA plant heats 300 kg of a 50 % glycol loop charge in the freeze-protected outdoor loop through a rise of 20 K. Its specific heat capacity is 3.4 kJ/(kg·K). Calculate the heat the charge absorbs.

4. Heating a massA plant heats 200 kg of a mild-steel casting batch in the stress-relief oven through a rise of 50 K. Its specific heat capacity is 0.5 kJ/(kg·K). Calculate the heat the charge absorbs.

5. Changing phaseAn ice-storage tank holds 400 kg of ice at 0 °C. The plant wants all of it melted to water at 0 °C. The latent heat of fusion for water is 334 kJ/kg. Calculate the heat the melt demands.

6. Changing phaseA flash tank receives 900 kg of saturated water at 100 °C and evaporates all of it at atmospheric pressure. The latent heat of vaporisation there is 2257 kJ/kg. Calculate the heat the evaporation demands.

7. The gas lawsAn accumulator holds 3 m³ of gas at 400 kPa absolute. It is compressed isothermally to 240 kPa absolute, with the temperature held constant throughout. Calculate the volume the gas then occupies.

8. The gas lawsA gasholder bell keeps its gas at constant pressure. It contains 2 m³ at 300 K, and the gas is then brought to 450 K. Calculate the volume at the new temperature.

9. The ideal gasA 0.02 m³ cylinder of nitrogen stands at 15000 kPa absolute and 280 K. Take R as 8.314 J/(mol·K). Determine the amount of gas in the cylinder.

10. The ideal gasA purge system must store 200 mol of nitrogen at 2000 kPa absolute and 350 K. Take R as 8.314 J/(mol·K). Determine the receiver volume that holds it.

11. Real gas correctionsA test cell holds 2 mol of gas in 0.01 m³ at 350 K, and the transducer reads 500 kPa absolute. An ideal gas at those conditions would read something else. Calculate the compressibility factor of the gas.

12. Real gas correctionsA 0.05 m³ bottle of gas stands at 15000 kPa absolute and 300 K. At those conditions the generalised chart gives a compressibility factor of 0.85. Determine the amount of gas actually in the bottle.

13. Vapour pressure curvesA closed feedwater tank holds water at 80 °C. The Antoine constants for water on the 1–100 °C fit are A = 8.07131, B = 1730.63 and C = 233.426, quoted for pressure in millimetres of mercury and temperature in degrees Celsius. Calculate the vapour pressure above the water.

14. Vapour pressure curvesA vacuum deaerator is held at 50 kPa absolute. Water in it boils at whatever temperature matches that pressure. The Antoine constants for water are A = 8.07131, B = 1730.63 and C = 233.426, for pressure in mmHg and temperature in °C. Determine the temperature at which the water boils in the vessel.

15. The State FinalLast job of the shift. A rigid instrument-air receiver reads 900 kPa(g) at 250 K, and the compressor room heats it to 300 K with the outlet valve shut. The barometer reads a round 100 kPa today, and the relief valve on the receiver is set at 1050 kPa(g). Work each line — every answer feeds the next. Determine whether the relief valve lifts, one line at a time.

16. The State FinalSame shift, the deaerator. 2000 kg of condensate arrives at 40 °C and must reach 100 °C, after which 200 kg of it flashes to steam. No calculator: call water's specific heat a round 4 kJ/(kg·K) and its latent heat of vaporisation a round 2250 kJ/kg — the numbers an operator keeps in their head. Determine the heat each stage of that duty demands.

Steam, Turbines & Power Cycles

17. Reading the steam tableA steam header carries a Bourdon gauge reading 298.7 kPa, and local atmospheric pressure is 101.3 kPa. The plant's saturated-steam table lists: 200 kPa → 120.2 °C · 400 kPa → 143.6 °C · 700 kPa → 165.0 °C. Determine the saturation temperature in the header.

18. Reading the steam tableA data logger on a saturated header reports the steam temperature as 416.75 K. The plant's saturated-steam table lists: 200 kPa → 120.2 °C · 400 kPa → 143.6 °C · 700 kPa → 165.0 °C. Determine the absolute pressure in the header.

19. Wet steam and qualityWet steam of quality 0.8 is measured at 2,349.7 kJ/kg. The saturated liquid enthalpy at the same pressure is h_f = 697.1 kJ/kg. Determine the latent heat of vaporisation at that pressure.

20. Wet steam and qualitySteam leaving a separator at 2,000 kPa is 95 % dry. At that pressure the table gives h_f = 908.6 kJ/kg and h_g = 2,798.4 kJ/kg. Calculate the specific enthalpy of the steam leaving the separator.

21. Degrees of superheatA superheater outlet header runs at 700 kPa absolute, and the thermowell in it reads 225.0 °C. The plant's saturated-steam table lists: 400 kPa → 143.6 °C · 700 kPa → 165.0 °C · 1,000 kPa → 179.9 °C. Determine the degrees of superheat carried by that steam.

22. Degrees of superheatA superheater outlet header runs at 1,000 kPa absolute, and the thermowell in it reads 259.9 °C. The plant's saturated-steam table lists: 700 kPa → 165.0 °C · 1,000 kPa → 179.9 °C · 2,000 kPa → 212.4 °C. Determine the degrees of superheat carried by that steam.

23. Flash steamA trap set passes 500 kg/h of condensate from a 700 kPa main into a 100 kPa flash vessel. The table gives h_f = 697.1 kJ/kg at the higher pressure, and h_f = 417.5 kJ/kg with h_fg = 2,258 kJ/kg at the lower. Determine the percentage that flashes and the flash steam rate it produces.

24. Flash steamA trap set passes 1,500 kg/h of condensate from a 1,000 kPa main into a 100 kPa flash vessel. The table gives h_f = 762.5 kJ/kg at the higher pressure, and h_f = 417.5 kJ/kg with h_fg = 2,258 kJ/kg at the lower. Determine the percentage that flashes and the flash steam rate it produces.

25. Turbine workSteam enters a turbine stop valve at 3,116.1 kJ/kg and leaves the last stage at 2,416.1 kJ/kg. Calculate the specific work the machine takes from the steam.

26. Turbine workA turbine passes 20 kg/s of steam and takes 650 kJ/kg out of it between the stop valve and the exhaust. Determine the shaft power the machine develops.

27. Isentropic efficiencyA turbine takes steam at 3,213.6 kJ/kg and exhausts it at 2,533.6 kJ/kg. An expansion at constant entropy to the same exhaust pressure would have ended at 2,413.6 kJ/kg. Determine the isentropic efficiency of the machine.

28. Isentropic efficiencyA supply contract guarantees an isentropic efficiency of 82 %. On the acceptance test the machine takes steam at 3,213.6 kJ/kg and exhausts it at 2,533.6 kJ/kg; the isentropic end point at that exhaust pressure is 2,413.6 kJ/kg. Determine whether the machine meets its guarantee.

29. Steam losses and loadsA steam-trap survey finds a blowing element on a header whose gauge reads 1,898.7 kPa; local atmospheric pressure is 101.3 kPa. The leak is treated as a sharp orifice of 12 mm². Calculate the steam lost through the orifice, in kilograms per hour.

30. Steam losses and loadsA spray desuperheater takes 30,000 kg/h of steam at 3,240 kJ/kg and must bring it down to 2,790 kJ/kg. The feedwater at the nozzle carries 540 kJ/kg. Determine the spray water rate the station requires.

31. Cycle efficienciesOver one hour of steady running, a plant burns fuel worth 100 MJ of heat into the working fluid and delivers 35 MJ of work at the coupling. Calculate the plant's thermal efficiency.

32. Cycle efficienciesA steam plant takes heat in at 280 °C and rejects it to cooling water at 38 °C. Determine the highest thermal efficiency any engine could reach between those two temperatures.

33. The Turbine HallAcceptance run, turbine hall. Steam reaches the stop valve at 3,250 kJ/kg and 20 kg/s. It leaves the last stage into the condenser at a quality of 0.9, where the rounded table gives h_f = 240 kJ/kg and h_fg = 2,400 kJ/kg. The condensate leaves at 240 kJ/kg and the feed pump lifts it to 250 kJ/kg before the boiler. Work each line — every answer feeds the next. Determine the shaft power and the cycle efficiency, one line at a time.

34. The Turbine HallBonus mark, no calculator. A vendor's brochure claims 32 % thermal efficiency for a cycle taking heat in at 600 K and rejecting it at 360 K. (Both temperatures are already absolute — the boss does not do conversions.) Determine whether that claim can be true.

Conduction & the Building Envelope

35. Fourier's lawA heat-flux survey measures 480 W crossing 50 mm of mineral wool on a flat duct wall. The slab is 10 m² in area and its conductivity is 0.04 W/(m·K). Determine the temperature difference across the slab.

36. Fourier's lawA heat-flux survey measures 9600 W crossing 150 mm of refractory in the boiler setting. The slab is 3 m² in area and its conductivity is 1.2 W/(m·K). Determine the temperature difference across the slab.

37. R-value of a layerA submittal sheet lists 50 mm of polyisocyanurate board with a declared thermal conductivity of 0.025 W/(m·K). Calculate the R-value of the layer.

38. R-value of a layerA roof specification calls for an insulation layer of RSI 3. The board on the truck is extruded polystyrene, conductivity 0.03 W/(m·K). Determine the thickness of board the specification demands.

39. Stacking the wallAn energy model lists a plant-office wall as three resistances in series: the inside and outside air films together at RSI 0.2, the cavity insulation at RSI 4, and the sheathing, gypsum and cladding together at RSI 0.4. Calculate the total R-value of the assembly.

40. Stacking the wallA wall assembly is measured at RSI 3 overall. Its air films are worth RSI 0.15 and its sheathing, gypsum and cladding together RSI 0.35. Determine the R-value of the cavity insulation inside it.

41. From R to UA wall assembly is calculated at RSI 1.25 overall, films included. The energy code that governs the job is written in U-factors, not R-values. Determine the U-factor of the assembly.

42. From R to UA window schedule quotes an opaque spandrel panel at U-0.25 W/(m²·K). The insulation submittal that has to match it is written in RSI. Determine the panel's total R-value.

43. Whole-wall heat lossA design-day calculation covers 15 m² of exterior wall at RSI 5, films included. Indoors is held at 20 °C and the design outdoor temperature is -10 °C. Calculate the heat lost through that wall at design conditions.

44. Whole-wall heat lossA thermographer measures RSI 4 across an assembly with 30 K from the warm face to the cold face. No area is recorded — the report is written per square metre. Calculate the heat flux through the assembly.

45. The framing penaltyA wall is insulated to RSI 5 between the studs. The path straight through a stud is worth RSI 1.2, and the framing takes 25% of the wall's area. Calculate the effective R-value of the whole wall.

46. The framing penaltyThe same build-up again: RSI 3.5 in the cavity, RSI 1 through a stud, framing factor 0.25. The design team wants to know how much of the batt's rating survives the framing. Determine the effective R-value, then the percentage of the cavity rating the wall actually delivers.

47. Pipes and the critical radiusA 8 m run of steam main is lagged from an inner radius of 40 mm to an outer radius of 120 mm. The lagging's conductivity is 0.05 W/(m·K), and there is 110 K from the inner face of the lagging to its outer face. Calculate the heat conducted through the lagging.

48. Pipes and the critical radiusA lagging with a conductivity of 0.04 W/(m·K) is being considered for a small 32 mm-radius line, and the outside film coefficient in that space is 5 W/(m²·K). Determine the critical radius for this lagging, then decide what a first thin wrap does to the line's heat loss.

49. The Envelope AuditLast audit of the day. One exterior wall: air films worth RSI 0.2, cavity insulation RSI 4.5, and sheathing, gypsum and cladding together RSI 0.3. The wall is 45 m² and the design condition puts 40 K across it. The code the job is built to caps walls at U-0.28 W/(m²·K). No calculator — the numbers are chosen to fit in your head, and every answer feeds the next. Determine whether this wall meets the code's U-factor cap, one line at a time.

50. The Envelope AuditBonus mark, on the way out: the estimator wants the same kind of wall — RSI 5 under 40 K — quoted per square metre rather than for the whole elevation. Determine the heat flux through the assembly.

Convection, Radiation & Heat Exchangers

51. Newton's law of coolingA lagged condensate line loses 4500 W through 2 m² of jacket, and the jacket's outside film coefficient is 45 W/(m²·K). Determine how far the jacket surface sits above the surrounding air.

52. Newton's law of coolingA bare section of steam main in a plant room presents 4 m² of surface. The still-air film coefficient on it is 25 W/(m²·K), and the surface stands 50 K above the room air. Calculate the convective heat rate off the bare section.

53. Dittus–BoelterWater flows turbulently through the 25 mm tubes of a shell-and-tube exchanger at Re = 50 000 and Pr = 6, and hot water is being COOLED inside the tubes. Take the water's thermal conductivity as 0.6 W/(m·K). Determine the inside film coefficient.

54. Dittus–BoelterWater flows turbulently through the 20 mm tubes of a shell-and-tube exchanger at Re = 100 000 and Pr = 6, and hot water is being COOLED inside the tubes. Take the water's thermal conductivity as 0.6 W/(m·K). Determine the inside film coefficient.

55. Radiation exchangeAn uninsulated boiler-front panel presents 1.5 m² at 127 °C. Its oxidised paint has an emissivity of 0.8. Calculate the total power the panel radiates.

56. Radiation exchangeThe same 1 m² panel, emissivity 0.95, runs at 527 °C while the boiler-house walls around it sit at 27 °C. Calculate the net radiant heat the panel loses to the room.

57. Building the UA rating sheet for lube oil against cooling water lists an inside film coefficient of 1000 W/(m²·K), an outside film coefficient of 1250 W/(m²·K), and a 3 mm wall of conductivity 15 W/(m·K). Calculate the clean overall heat transfer coefficient.

58. Building the UA rating sheet for boiler feedwater behind a heavy-wall tube lists an inside film coefficient of 5000 W/(m²·K), an outside film coefficient of 2500 W/(m²·K), and a 6 mm wall of conductivity 15 W/(m·K). Calculate the clean overall heat transfer coefficient.

59. Log mean temperature differenceAn exchanger on turbine lube oil against cooling water cools the hot stream from 150 °C to 90 °C while the cold stream rises from 30 °C to 70 °C. The streams run COUNTERFLOW — they enter at opposite ends of the shell. Calculate the log mean temperature difference for this arrangement.

60. Log mean temperature differenceAn exchanger on hot condensate against feedwater cools the hot stream from 160 °C to 100 °C while the cold stream rises from 40 °C to 60 °C. The streams run COUNTERFLOW — they enter at opposite ends of the shell. Calculate the log mean temperature difference for this arrangement.

61. Rating the exchangerOn a one-shell two-tube-pass cooler, water at 4 kg/s is to be heated through 15 K; take cₚ = 4.2 kJ/(kg·K). The design overall coefficient is 800 W/(m²·K), the log mean temperature difference is 25 K, and the arrangement gives a correction factor F = 0.9. Determine the heat transfer area the unit requires.

62. Rating the exchangerOn a counterflow plate pack, water at 5 kg/s is to be heated through 20 K; take cₚ = 4.2 kJ/(kg·K). The design overall coefficient is 1200 W/(m²·K), the log mean temperature difference is 25 K, and the arrangement gives true counterflow, so F = 1.0. Determine the heat transfer area the unit requires.

63. Effectiveness–NTUA counterflow plate exchanger carries hot water at 1 kg/s on one side and cold water at 2 kg/s on the other; take cₚ = 4.2 kJ/(kg·K) for both. Its overall coefficient is 900 W/(m²·K) over 7 m² of plate. Determine the exchanger's number of transfer units.

64. Effectiveness–NTUA counterflow plate exchanger carries hot water at 2 kg/s on one side and cold water at 2 kg/s on the other; take cₚ = 4.2 kJ/(kg·K) for both. Its overall coefficient is 1200 W/(m²·K) over 7 m² of plate. Determine the exchanger's number of transfer units.

65. Transient coolingBefore any transient calculation, a fireclay refractory block is checked for a uniform internal temperature. Its characteristic length is 50 mm, the quench film coefficient on it is 50 W/(m²·K), and the solid's own conductivity is 1 W/(m·K). Determine whether the lumped-capacitance method may be used on this body.

66. Transient coolingBefore any transient calculation, a stainless baffle plate is checked for a uniform internal temperature. Its characteristic length is 15 mm, the quench film coefficient on it is 300 W/(m²·K), and the solid's own conductivity is 15 W/(m·K). Determine whether the lumped-capacitance method may be used on this body.

67. The Exchanger FinalAcceptance day on a counterflow water-to-water heat exchanger. Both streams run at 2.5 kg/s, and today cₚ = 4 kJ/(kg·K) for both — equal flows and equal specific heats, so the two capacity rates are equal. The hot side enters at 110 °C and leaves at 70 °C; the cold side enters at 60 °C. The overall coefficient is 800 W/(m²·K). Work each line — every answer feeds the next. Determine the heat transfer surface this unit must carry, one line at a time.

68. The Exchanger FinalBonus mark. The same balanced pair of streams — equal flows, equal specific heats — is offered in a cheaper PARALLEL-FLOW shell instead: both fluids enter at the same end. The hot side still enters at 150 °C and the cold side at 50 °C, and the duty specification still calls for cold water leaving at 95 °C. Determine whether the parallel-flow unit can meet that outlet specification.

Answer key

  1. 649.5 kPa
  2. 21.3 kPa
  3. 20.4 MJ
  4. 5 MJ
  5. 133.6 MJ
  6. 2031.3 MJ
  7. 5 m³
  8. 3 m³
  9. 128.9 mol
  10. 0.291 m³
  11. 0.859 (no unit)
  12. 353.7 mol
  13. 47.27 kPa
  14. 81.39 °C
  15. 1000 kPa
  16. 480 MJ
  17. 143.6 °C
  18. 400 kPa absolute
  19. 2065.75 kJ/kg
  20. 2703.91 kJ/kg
  21. 165 °C
  22. 179.9 °C
  23. 12.3826 %
  24. 15.279 %
  25. 700 kJ/kg
  26. 13000 kW
  27. 680 kJ/kg
  28. 0.85 (no unit)
  29. 125.862 kg/h
  30. 6000 kg/h
  31. 0.35 (no unit)
  32. 0.437494 (no unit)
  33. 2400 kJ/kg
  34. 0.4 (no unit)
  35. 60 K
  36. 400 K
  37. 2 m²·K/W
  38. 90 mm
  39. 4.6 m²·K/W
  40. 2.5 m²·K/W
  41. 0.8 W/(m²·K)
  42. 4 m²·K/W
  43. 90 W
  44. 7.5 W/m²
  45. 2.791 m²·K/W
  46. 2.154 m²·K/W
  47. 251.6 W
  48. 8 mm
  49. 5 m²·K/W
  50. 8 W/m²
  51. 50 K
  52. 5000 W
  53. 226.125 —
  54. 393.707 —
  55. 1741.94 W
  56. 21628.2 W
  57. 500 W/(m²·K)
  58. 1000 W/(m²·K)
  59. 69.5212 K
  60. 78.3046 K
  61. 252 kW
  62. 420 kW
  63. 4.2 kW/K
  64. 8.4 kW/K
  65. 2.5 —
  66. 0.3 —
  67. 400 kW
  68. 100 °C