Convection, Radiation & Heat Exchangers
Newton's Law of Cooling (Q = hAΔT)
Where
- Q˙= Heat transfer rate (W)
- h= Convection coefficient (W/(m²·K))
- A= Surface area (m²)
- ΔT= Surface-to-fluid ΔT (C°)
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
Where
- R= Film resistance (K/W)
- h= Convection coefficient (W/(m²·K))
- A= Wetted area (m²)
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
Where
- Re= Reynolds number
- ρ= Fluid density (kg/m³)
- v= Flow velocity (m/s)
- D= Characteristic length (mm)
- μ= 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
Where
- Pr= Prandtl number
- μ= Dynamic viscosity (Pa·s)
- cp= Specific heat (J/(kg·K))
- k= 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
Where
- Nu= Nusselt number
- h= Convection coefficient (W/(m²·K))
- L= Characteristic length (m)
- k= 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
Where
- Gr= Grashof number
- g= Gravitational acceleration (m/s²)
- β= Volumetric expansion coefficient (1/K)
- ΔT= Surface-to-fluid temperature difference (C°)
- L= Characteristic length (m)
- ν= 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, 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βΔ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 is the viscosity resisting twice over, once in setting up the motion and once in damping it. Note the presence of g, 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, L 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 L 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 L 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 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, Gr⋅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
Where
- Ra= Rayleigh number
- Gr= Grashof number
- Pr= Prandtl number
The Rayleigh number is a product, Ra=GrPr, 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/(να), 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 να 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 Ra≈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 Ra≈109 for the laminar-to-turbulent transition on a vertical plate, and that one is a CONVENTION. The real transition wanders roughly between 108 and 1010 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 106, a room wall around 109 to 1010, and the Earth's mantle something near 1020 — which is why the mantle convects despite being rock. And the Nusselt numbers that come out grow slowly, roughly as Ra1/4 in the laminar range and Ra1/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
Where
- Nu= Nusselt number
- Re= Reynolds number
- Pr= Prandtl number
- n= 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 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, 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
Where
- P= Radiated power (W)
- ε= Emissivity
- A= Surface area (m²)
- T= 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 T3, and the average energy each one carries scales as T, so the product goes as T4. 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×10−8×1.8×3014≈823 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σT4, 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(T4−Tsurr4), 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)4≈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
Where
- Q˙= Net radiant heat rate (W)
- ε= Emissivity
- A= Surface area (m²)
- T1= Surface temperature (°C)
- T2= 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
Where
- ht= Combined coefficient (W/(m²·K))
- hc= Convection coefficient (W/(m²·K))
- ε= Surface emissivity
- Ts= Surface temperature (°C)
- Tsur= 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). Add it to the convective film and a single ht drives the whole surface with plain old Q̇ = htAΔ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 hr = 6.8 W/(m²·K), larger than the 5 W/(m²·K) of still-air natural convection beside it, giving ht = 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: hr 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)
Where
- U= Overall coefficient (W/(m²·K))
- hi= Inside film coefficient (W/(m²·K))
- L= Wall thickness (mm)
- k= Wall thermal conductivity (W/(m·K))
- ho= 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
Where
- Uf= Fouled overall coefficient (W/(m²·K))
- Uc= Clean overall coefficient (W/(m²·K))
- hf= 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 Rf 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/Rf, so it lives in the same W/(m²·K) picker as the rest of the shard — 0.0002 m²·K/W is hf = 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 Rf 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
Where
- Uf= Fouled (service) coefficient (W/(m²·K))
- Uc= Clean coefficient (W/(m²·K))
- Rf= 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 Rf 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 Rf 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/Uf−1/Uc 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)
Where
- ΔTlm= Log mean temperature difference (C°)
- Th,in= Hot stream inlet (°C)
- Th,out= Hot stream outlet (°C)
- Tc,in= Cold stream inlet (°C)
- Tc,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=(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 only: the four terminal temperatures sit inside a logarithm and a difference at once, so recovering an inlet temperature from a known ΔTlm 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)
Where
- ΔTlm= Log mean temperature difference (C°)
- Th,in= Hot stream inlet (°C)
- Th,out= Hot stream outlet (°C)
- Tc,in= Cold stream inlet (°C)
- Tc,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)
Where
- Q˙= Exchanger duty (kW)
- U= Overall coefficient (W/(m²·K))
- A= Heat transfer area (m²)
- F= LMTD correction factor
- ΔTlm= Log mean temperature difference (C°)
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)
Where
- Q˙= Stream duty (kW)
- m˙= Mass flow rate (kg/h)
- cp= Specific heat (J/(kg·K))
- ΔT= Temperature change (C°)
Every exchanger calculation has two halves that must agree. The transfer side says Q̇ = UAF·ΔTlm; 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)
Where
- NTU= Number of transfer units
- U= Overall coefficient (W/(m²·K))
- A= Heat transfer area (m²)
- m˙= Minimum stream mass flow (kg/h)
- cp= 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)
Where
- Cr= Capacity rate ratio
- m˙min= Minimum stream mass flow (kg/h)
- cmin= Minimum stream specific heat (J/(kg·K))
- m˙max= Maximum stream mass flow (kg/h)
- cmax= 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 ε=1−e−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 Cmin 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)
Where
- Q˙max= Maximum possible duty (kW)
- m˙min= Minimum stream mass flow (kg/h)
- cmin= Minimum stream specific heat (J/(kg·K))
- Th,in= Hot stream inlet (°C)
- Tc,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)
Where
- ε= Effectiveness
- Q˙= Actual duty (kW)
- 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)
Where
- ε= Effectiveness
- NTU= Number of transfer units
- Cr= 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=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 ε=1−e−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
Where
- Bi= Biot number
- h= Film coefficient (W/(m²·K))
- Lc= Characteristic length (m)
- k= 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) 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, Lc = 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
Where
- Fo= Fourier number
- k= Thermal conductivity (W/(m·K))
- t= Elapsed time (s)
- ρ= Density (kg/m³)
- cp= Specific heat (J/(kg·K))
- L= 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
Where
- τ= Thermal time constant (s)
- ρ= Density (kg/m³)
- V= Body volume (L)
- cp= Specific heat (J/(kg·K))
- h= Film coefficient (W/(m²·K))
- A= Surface area (m²)
ρ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
Where
- T= Temperature at time t (°C)
- T0= Initial temperature (°C)
- T∞= Fluid temperature (°C)
- t= Elapsed time (s)
- τ= 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×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 absolute — The 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 absolute — A 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 mass — A 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 mass — A 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 phase — An 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 phase — A 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 laws — An 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 laws — A 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 gas — A 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 gas — A 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 corrections — A 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 corrections — A 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 curves — A 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 curves — A 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 Final — Last 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 Final — Same 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 table — A 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 table — A 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 quality — Wet 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 quality — Steam 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 superheat — A 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 superheat — A 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 steam — A 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 steam — A 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 work — Steam 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 work — A 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 efficiency — A 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 efficiency — A 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 loads — A 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 loads — A 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 efficiencies — Over 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 efficiencies — A 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 Hall — Acceptance 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 Hall — Bonus 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 law — A 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 law — A 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 layer — A 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 layer — A 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 wall — An 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 wall — A 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 U — A 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 U — A 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 loss — A 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 loss — A 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 penalty — A 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 penalty — The 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 radius — A 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 radius — A 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 Audit — Last 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 Audit — Bonus 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 cooling — A 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 cooling — A 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–Boelter — Water 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–Boelter — Water 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 exchange — An 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 exchange — The 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 U — A 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 U — A 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 difference — An 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 difference — An 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 exchanger — On 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 exchanger — On 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–NTU — A 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–NTU — A 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 cooling — Before 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 cooling — Before 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 Final — Acceptance 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 Final — Bonus 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.