Heat Exchanger Sizing: LMTD and Effectiveness-NTU Methods
Two methods dominate heat exchanger engineering: the log-mean temperature difference (LMTD) method, which sizes a new exchanger when all four terminal temperatures are known, and the effectiveness-NTU (ε-NTU) method, which rates an existing exchanger from its area, flow rates and inlet conditions. Both methods are thermodynamically identical — they give the same result for a given exchanger — but they differ in which quantities you know at the design stage, and choosing the wrong method for the task leads to avoidable iteration.
This guide covers the LMTD formula, the F-correction factor for non-counterflow arrangements, the ε-NTU correlations for the most common flow types, and a fully worked counterflow sizing example. Every formula below matches the engine used in the MechanixCalc heat exchanger calculator, so the numbers you work by hand are directly reproducible in the tool.
The LMTD method — sizing from four terminal temperatures
The log-mean temperature difference is the driving-force average of the two terminal temperature differences across the exchanger. For a pure counterflow arrangement the LMTD is the highest achievable for a given temperature program, which is why counterflow is the thermodynamically preferred layout: less area is needed for the same duty. The required heat-transfer area follows directly from the overall heat-transfer coefficient U and the heat duty Q.
For any arrangement other than pure counterflow (parallel flow, 1-2 shell-and-tube, or cross-flow), the counterflow LMTD is multiplied by a correction factor F ≤ 1.0 that accounts for the departure from ideal counterflow. When F drops below 0.75 the required area becomes very sensitive to small changes in the temperature program, and adding shell passes or switching to counterflow should be considered. The LMTD engine in the calculator rejects temperature programs where a terminal ΔT ≤ 0 — a second-law violation — rather than clamping it to a small value and reporting a spuriously large area.
LMTD = (ΔT₁ − ΔT₂) / ln(ΔT₁ / ΔT₂)where ΔT₁ = T_h,in − T_c,out (hot-end temperature difference); ΔT₂ = T_h,out − T_c,in (cold-end temperature difference); ln = natural logarithm. When ΔT₁ = ΔT₂ (isothermal streams crossing), LMTD = ΔT₁.
A = Q / (U · F · LMTD)where Q = heat duty (W); U = overall heat-transfer coefficient (W/m²·K); F = LMTD correction factor (dimensionless, 0 < F ≤ 1.0); LMTD = log-mean temperature difference (K)
The F-correction factor for non-counterflow arrangements
For a 1-2 shell-and-tube exchanger (one shell pass, two tube passes) the correction factor F is given by the Bowman-Mueller-Nagle (1940) closed form, which is a function of the dimensionless temperature parameters P (cold-stream thermal effectiveness) and R (heat-capacity rate ratio). P and R are defined from the four terminal temperatures alone, without knowing the flow rates. The Bowman formula returns no valid result when the temperature program (P, R) lies outside the physical single-shell-pass domain — the calculator reports this as a domain violation and suggests adding shell passes rather than silently returning a floored constant.
For a single-pass cross-flow exchanger with both fluids unmixed, there is no simple closed-form F expression analogous to the Bowman formula. The calculator uses the effectiveness-NTU inversion method (Incropera & DeWitt §11.4): it computes the exchanger effectiveness ε and capacity ratio C* from (P, R), finds NTU_counterflow and NTU_crossflow for the same ε and C*, and returns F = NTU_counterflow / NTU_crossflow. This matches the published F-correction charts (TEMA / ESDU 86018). For pure counterflow F = 1.0 exactly — no calculation needed.
P = (T_c,out − T_c,in) / (T_h,in − T_c,in) ; R = (T_h,in − T_h,out) / (T_c,out − T_c,in)where P = cold-stream thermal effectiveness (0 < P < 1); R = heat-capacity rate ratio of hot to cold stream (also equal to C_c / C_h); T_h,in, T_h,out = hot-stream inlet and outlet temperatures; T_c,in, T_c,out = cold-stream inlet and outlet temperatures
The effectiveness-NTU method — rating an existing exchanger
The effectiveness-NTU method is used for rating: you know the exchanger (area A, overall U) and the stream conditions (inlet temperatures and flow rates), and want to find the heat duty and outlet temperatures without guessing outlet temperatures and iterating. The number of transfer units NTU = U·A / C_min is a dimensionless measure of exchanger size relative to the weaker stream. The capacity ratio C* = C_min / C_max (0 ≤ C* ≤ 1) describes the stream balance; C* = 0 represents a condensing or evaporating stream (infinite capacity rate on one side).
The effectiveness ε is the ratio of actual to maximum possible duty. Q_max is limited by the stream with the smaller heat-capacity rate C_min = min(ṁ·c_p), because that stream undergoes the largest temperature change. Once ε is found from the ε-NTU correlation for the flow arrangement, the actual duty Q_actual = ε · Q_max, and the outlet temperatures follow from Q_actual / C_h and Q_actual / C_c. The counterflow and 1-2 shell-and-tube ε-NTU expressions have compact closed forms; the both-unmixed cross-flow expression (Incropera Eq. 11.32) is a series approximation used in the calculator.
NTU = U · A / C_min ; C* = C_min / C_maxwhere U = overall heat-transfer coefficient (W/m²·K); A = heat-transfer area (m²); C_min = min(ṁ_h·c_p,h , ṁ_c·c_p,c) (W/K); C_max = max(ṁ_h·c_p,h , ṁ_c·c_p,c) (W/K); C* = capacity ratio (0 ≤ C* ≤ 1)
ε = (1 − exp[−NTU·(1 − C*)]) / (1 − C*·exp[−NTU·(1 − C*)]) ; (C* < 1)where ε = exchanger effectiveness (0 ≤ ε ≤ 1); NTU = number of transfer units; C* = C_min / C_max. For C* = 1 (balanced counterflow): ε = NTU / (1 + NTU).
Q_max = C_min · (T_h,in − T_c,in) ; Q_actual = ε · Q_maxwhere Q_max = thermodynamic upper bound on heat transfer (W); Q_actual = heat actually transferred (W); T_h,in − T_c,in = maximum available temperature driving force (K)
Fouling resistance and the design area margin
A clean exchanger seldom stays clean in service. Fouling deposits on the heat-transfer surfaces add thermal resistance and reduce the effective overall coefficient below U_clean. The fouled coefficient U_fouled is found by adding the hot-side and cold-side fouling resistances from the TEMA fouling tables to the clean overall resistance 1/U_clean. Because A = Q / (U·F·LMTD), a lower U_fouled means a larger required area — the extra area A_margin = A_fouled − A_clean is the design allowance for fouling over the service life.
The cleanliness factor CF = U_fouled / U_clean (always ≤ 1) is a compact way to express the severity of the fouling allowance. A CF of 0.80 means the exchanger must be 25% larger than the clean-area calculation to maintain the design duty after fouling. TEMA fouling resistances for common service fluids are built into the calculator's fouling panel.
1 / U_fouled = 1 / U_clean + R_f,hot + R_f,coldwhere U_clean = clean overall heat-transfer coefficient (W/m²·K); R_f,hot, R_f,cold = hot-side and cold-side fouling resistances (m²·K/W, from TEMA fouling tables); U_fouled = fouled overall coefficient (W/m²·K)
CF = U_fouled / U_clean ; A_margin = A_fouled − A_cleanwhere CF = cleanliness factor (0 < CF ≤ 1); A_clean = area required with clean surfaces (m²); A_fouled = area required after fouling (m²); A_margin = additional area to specify for fouling allowance (m²)
Worked example
Size a counterflow shell-and-tube heat exchanger to cool hot water from 80 °C to 50 °C using cold water entering at 25 °C and leaving at 45 °C, with a heat duty of 250 kW and an overall U = 600 W/m²·K.
Given
- Hot inlet T_h,in80 °C
- Hot outlet T_h,out50 °C
- Cold inlet T_c,in25 °C
- Cold outlet T_c,out45 °C
- Heat duty Q250 kW
- Overall U600 W/m²·K
- Flow arrangementCounterflow (F = 1.0)
Result
- LMTD≈ 29.7 K
- F correction factor1.000 (counterflow)
- Required heat-transfer area A≈ 14.0 m²
- Compute the counterflow terminal temperature differences: ΔT₁ = T_h,in − T_c,out = 80 − 45 = 35 K; ΔT₂ = T_h,out − T_c,in = 50 − 25 = 25 K.
- Compute the LMTD: LMTD = (ΔT₁ − ΔT₂) / ln(ΔT₁/ΔT₂) = (35 − 25) / ln(35/25) = 10 / ln(1.4).
- ln(1.4) ≈ 0.3365, so LMTD = 10 / 0.3365 ≈ 29.7 K.
- For pure counterflow, F = 1.0 exactly. Apply the area equation: A = Q / (U · F · LMTD) = (250 × 1 000) / (600 × 1.0 × 29.7).
- A = 250 000 / 17 820 ≈ 14.0 m².
Illustrative — verify against your actual flow rates, fluid properties and fouling allowance. For a 1-2 shell-and-tube arrangement, F < 1.0 (Bowman-Mueller-Nagle formula) and the required area will be larger. Add a fouling margin: A_design = A_clean / CF, where CF = U_fouled / U_clean from the TEMA fouling tables.
Do it on your own numbers
Run the full LMTD area sizing with Bowman F-correction, NTU-effectiveness rating, fouling resistance analysis and tube sheet layout for your exchanger. Free 30-minute preview, no sign-up.
Open the Heat ExchangerFrequently asked questions
What is the difference between the LMTD method and the NTU-effectiveness method?
The LMTD method is used for sizing — you specify all four terminal temperatures and the heat duty to find the required area. The NTU-effectiveness method is used for rating — you specify an existing exchanger (area A, overall U) and the stream flow rates and inlet temperatures to find the outlet temperatures and actual duty. Both are thermodynamically equivalent; the choice depends on which quantities you know at the design stage. Knowing all four outlet temperatures? Use LMTD. Only know inlet conditions? Use NTU.
What is the F-correction factor and when is it less than 1?
The F factor corrects the counterflow LMTD for a non-counterflow arrangement. Pure counterflow gives F = 1.0. A 1-2 shell-and-tube exchanger (one shell pass, two tube passes) gives F < 1.0 via the Bowman-Mueller-Nagle (1940) formula because the streams are not in pure counterflow. Single-pass cross-flow (both fluids unmixed) also gives F < 1.0, found by the effectiveness-NTU inversion method. When F drops below 0.75, the required area becomes very sensitive to the temperature program — adding shell passes or using counterflow is advisable.
How do I calculate the overall heat-transfer coefficient U?
U combines the tube-side film coefficient h_i, the wall conduction resistance, the shell-side film coefficient h_o, and fouling resistances in series: 1/U = 1/h_i + (t_wall/k_wall) + 1/h_o + R_f,hot + R_f,cold. Film coefficients are typically estimated from the Dittus-Boelter correlation (Nu = 0.023·Re⁰·⁸·Pr^n, n = 0.4 for heating / 0.3 for cooling) for turbulent tube-side flow, and from Bell-Delaware or Kern correlations on the shell side. In practice, U is often taken from published ranges for the service type (water-water, air-water, steam-water) and then refined from operating data.
Why does a temperature cross prevent sizing a 1-2 shell-and-tube exchanger?
A temperature cross occurs when the cold stream outlet temperature exceeds the hot stream outlet temperature (T_c,out > T_h,out), meaning the streams have passed each other thermally. In a single shell pass of a 1-2 exchanger, the two tube passes run in both counterflow and parallel-flow directions simultaneously, so the stream with the lower outlet temperature can never recover — the exchanger cannot meet the program, and the F-correction formula returns no valid result. The solution is to use multiple shell passes in series (each shell approaches pure counterflow) or to switch to a true counterflow arrangement.
Is the heat exchanger calculator free?
You can use it during a free 30-minute preview with no sign-up, and a free 14-day account trial unlocks every calculator with no credit card. The branded PDF engineering report and saved calculations are part of a paid plan.
Related
- Pipe Flow (Darcy-Weisbach)Size the connecting pipework and check pressure drop on both HX streams.
- Pump & Fan SelectionSelect the circulation pump to overcome the heat exchanger pressure drop and system resistance.
- Pressure Vessels (ASME VIII)Verify the shell-side and tube-side pressure containment to ASME VIII.
- NPSH and pump cavitation (ISO 9906)
- Pipe pressure drop — Darcy-Weisbach