Thermal Engineering
Effectiveness–NTU Method: Rating a Heat Exchanger Without Knowing Outlet Temperatures
Hand an engineer a heat exchanger with a known 12 m² of tube area, a hot stream entering at 90 °C, and a cold stream entering at 20 °C — but no outlet temperatures — and the classic log-mean-temperature-difference (LMTD) approach stalls, because LMTD needs the very outlet temperatures you are trying to find. The Effectiveness–NTU method sidesteps that dead-end. It reframes performance as a single dimensionless number, effectiveness ε (0 to 1), that says what fraction of the thermodynamically maximum possible heat transfer the device actually achieves.
Developed and popularized by W. M. Kays and A. L. London in their 1955 monograph Compact Heat Exchangers, the ε–NTU method expresses ε as a closed-form function of two groups: the Number of Transfer Units, NTU = UA/C_min, and the capacity-rate ratio, C_r = C_min/C_max. Once those two numbers are known, effectiveness — and therefore heat duty and both outlet temperatures — falls out directly, with no iteration.
- TypeHeat-exchanger rating/sizing method
- Key equationε = q/q_max, NTU = UA/C_min
- Effectiveness range0 ≤ ε ≤ 1 (ε_max = 1)
- Formalized byKays & London, 1955
- Best forRating problems (outlet temps unknown)
- Typical design NTU0.5–5 (ε ≈ 0.4–0.99)
Interactive visualization
Press play, or step through manually. The visualization is yours to drive — try it before reading on.
Watch the 60-second explainer
A condensed visual walkthrough — narrated, captioned, under a minute.
What It Is and Where It's Used
The Effectiveness–NTU method is a design tool for rating and sizing heat exchangers. Its defining move is the effectiveness, ε = q / q_max, the ratio of actual heat duty to the maximum thermodynamically possible duty for the given inlet conditions. That maximum is q_max = C_min·(T_h,in − T_c,in), where C = ṁ·c_p is the capacity rate (mass flow times specific heat) and C_min is the smaller of the two streams' capacity rates.
- Rating problems: geometry and flows are fixed (UA known); you want the heat duty and outlet temperatures. ε–NTU solves these in one pass, where LMTD would require iteration.
- Compact exchangers: automotive radiators, HVAC coils, gas-turbine recuperators, and plate exchangers, where Kays & London's dimensionless framing shines.
It appears in every mechanical-engineering heat-transfer course (Incropera, Holman, Çengel) and in software from EES to aspenONE. Any device that transfers heat between two single-phase or phase-changing streams — recuperators, intercoolers, oil coolers, condensers, evaporators — can be characterized by a single ε.
How It Works: The Derivation
Start with two streams exchanging heat over area A with overall coefficient U. Write energy balances: q = C_h(T_h,in − T_h,out) = C_c(T_c,out − T_c,in). The maximum duty occurs when the stream with the smaller capacity rate undergoes the full inlet temperature difference, so q_max = C_min·ΔT_max, with ΔT_max = T_h,in − T_c,in. It must be C_min, because C_max cannot span the whole ΔT without violating the second law.
Integrating the local heat-transfer equation dq = U·dA·ΔT along the exchanger, subject to the flow arrangement, yields ε purely as a function of two groups:
- NTU = UA / C_min — a dimensionless "size" (how much conductance you have per unit of the limiting capacity).
- C_r = C_min / C_max — the capacity-rate ratio, 0 to 1.
For counterflow: ε = [1 − exp(−NTU(1 − C_r))] / [1 − C_r·exp(−NTU(1 − C_r))]. For parallel flow: ε = [1 − exp(−NTU(1 + C_r))] / (1 + C_r). When C_r → 0 (a boiling or condensing stream), both collapse to ε = 1 − exp(−NTU).
Key Quantities and a Worked Example
Consider a counterflow oil cooler. Hot oil: ṁ_h = 0.5 kg/s, c_p = 2.0 kJ/kg·K → C_h = 1000 W/K, entering at 90 °C. Cooling water: ṁ_c = 0.4 kg/s, c_p = 4.18 kJ/kg·K → C_c = 1672 W/K, entering at 20 °C. Overall UA = 1200 W/K.
- C_min = 1000 W/K (oil), C_max = 1672 W/K, so C_r = 0.598.
- NTU = UA/C_min = 1200/1000 = 1.2.
- Counterflow: exponent = NTU(1 − C_r) = 1.2·0.402 = 0.482; exp(−0.482) = 0.617. ε = (1 − 0.617)/(1 − 0.598·0.617) = 0.383/0.631 = 0.607.
Then q_max = C_min·ΔT_max = 1000·(90 − 20) = 70 kW, so q = ε·q_max = 0.607·70 = 42.5 kW. Outlet temperatures: T_h,out = 90 − q/C_h = 90 − 42.5 = 47.5 °C; T_c,out = 20 + q/C_c = 20 + 42.5/1.672 = 45.4 °C. No iteration, no LMTD guess.
Design, Selection, and Operation in Practice
Engineers use ε–NTU in two directions. For rating (known UA), read ε from the equation or the classic Kays–London charts and compute duty. For sizing (target ε known), invert the relation to get NTU, then A = NTU·C_min/U. The counterflow inverse is NTU = [1/(C_r − 1)]·ln[(ε − 1)/(ε·C_r − 1)].
- Diminishing returns: effectiveness rises steeply up to NTU ≈ 2–3, then flattens. Pushing from ε = 0.90 to 0.95 can double the area — a core capital-vs-performance trade-off.
- Arrangement matters most at high C_r: at C_r = 1, counterflow reaches ε = 0.75 at NTU = 3, but parallel flow caps near ε = 0.50 no matter how large the exchanger.
- Recuperators and regenerators (gas-turbine, cryogenic) are specified directly in ε (often 0.85–0.95) because that number sets fuel savings.
Practical UA already folds in fouling resistance and fin efficiency, so operators track ε over time as a health metric — a dropping ε flags fouling before a pressure alarm does.
Comparison to LMTD and Related Methods
The LMTD (log-mean-temperature-difference) method solves q = U·A·F·ΔT_lm, where ΔT_lm = (ΔT₁ − ΔT₂)/ln(ΔT₁/ΔT₂) and F is a correction factor for non-counterflow geometry. LMTD is elegant when all four terminal temperatures are known (a design/sizing problem) — you compute area directly. But in a rating problem the outlet temperatures are unknown, ΔT_lm can't be formed, and LMTD forces a guess-and-iterate loop.
- ε–NTU: best when outlets are unknown; algebraic, non-iterative; ε is a direct performance metric.
- LMTD: best when outlets are known and you want area; the F-factor plot handles crossflow/multipass.
- P-NTU (thermal effectiveness P): a variant used in TEMA/HEDH standards that references temperature change to one specific stream rather than C_min, convenient for multi-shell trains.
The two methods are mathematically equivalent — ε–NTU and LMTD are algebraic rearrangements of the same energy and rate balances — so they always give identical answers; the choice is about which unknowns you start with.
Limits, Failure Modes, and Significance
The ε–NTU relations rest on assumptions that, when violated, quietly corrupt results:
- Constant U and c_p: the closed-form ε expressions assume both are uniform. For viscous oils or near-critical fluids where c_p varies strongly, engineers segment the exchanger into zones, each with local ε–NTU.
- Phase change: during boiling or condensing, temperature is constant, so C_max → ∞ and C_r → 0. The device then follows ε = 1 − exp(−NTU) regardless of arrangement — a useful simplification, but you must switch zones at the saturation point.
- Fouling drift: deposits raise thermal resistance, lowering U and thus NTU and ε; a recuperator specified at ε = 0.90 may sag to 0.80, cutting recovered heat by ~10% and raising fuel burn.
- Longitudinal conduction and flow maldistribution erode effectiveness in high-NTU (>5) compact cores, an effect Kays & London themselves quantified.
Its significance is enduring: ε–NTU turns a messy coupled ODE problem into a two-number lookup, making it the default language for specifying recuperators, HVAC coils, and any exchanger where the answer, not the geometry, is the starting point.
| Configuration | NTU = 1, C_r = 1 | NTU = 3, C_r = 1 | NTU = 3, C_r = 0 |
|---|---|---|---|
| Counterflow | 50% | 75% | 95% |
| Parallel flow | 43% | 50% | 95% |
| Shell-and-tube (1 shell pass) | 46% | 58% | 95% |
| Crossflow (both unmixed) | 48% | 69% | 95% |
| Any config, C_r = 0 (boiler/condenser) | 63% | 95% | 95% |
Frequently asked questions
Why use ε–NTU instead of LMTD?
LMTD requires all four inlet and outlet temperatures to form the log-mean difference. In a rating problem the outlet temperatures are unknown, so LMTD demands iteration, while ε–NTU gives effectiveness — and thus duty and both outlets — algebraically in one pass. Use LMTD for sizing when outlets are known; use ε–NTU for rating when they are not.
What exactly is NTU and why is it dimensionless?
NTU = UA/C_min, the product of overall coefficient and area divided by the smaller capacity rate. Units are (W/K)/(W/K) = dimensionless. It represents the thermal "size" of the exchanger — how much heat-transfer conductance you have relative to the limiting stream's ability to absorb it. Higher NTU means a larger, more effective exchanger, with strongly diminishing returns above NTU ≈ 3.
Why is q_max based on C_min and not C_max?
The maximum possible heat transfer occurs when the stream with the smaller capacity rate experiences the full inlet-to-inlet temperature difference. If you tried to give C_max that full swing, C_min would have to exceed the temperature difference it can thermodynamically span, violating the second law. So q_max = C_min·(T_h,in − T_c,in) is the correct, physically attainable ceiling.
What does C_r = 0 mean physically?
C_r = C_min/C_max = 0 means one stream has an effectively infinite capacity rate — which happens during a phase change (boiling or condensing), since the stream's temperature stays constant while it absorbs or rejects latent heat. In that case every flow arrangement collapses to the same simple relation, ε = 1 − exp(−NTU). Condensers and evaporators are the classic C_r = 0 cases.
Why does counterflow beat parallel flow?
In counterflow the streams move oppositely, so the temperature difference stays more uniform along the length and the cold outlet can approach the hot inlet — allowing ε above 0.75 even at C_r = 1. In parallel flow both streams start far apart and converge toward a common temperature, capping ε near 0.50 at C_r = 1 no matter how large the area. Counterflow is thermodynamically superior for the same UA.
How do I size an exchanger if I know the target effectiveness?
Invert the ε–NTU relation. For counterflow, NTU = [1/(C_r − 1)]·ln[(ε − 1)/(ε·C_r − 1)]; for C_r = 0, NTU = −ln(1 − ε). Then the required area is A = NTU·C_min/U. Because ε flattens above NTU ≈ 3, targeting very high effectiveness (>0.95) demands disproportionately large area, so designers weigh performance against capital and pressure-drop cost.