Heat Exchanger Calculator — LMTD Area Sizing, NTU-Effectiveness & Tube Sheet Layout (TEMA / Bowman)
Governing standard: TEMA· TEMA (LMTD F-correction charts) · Bowman-Mueller-Nagle 1940 (1-2 shell-and-tube closed-form F) · Incropera & DeWitt §11.4 (NTU-effectiveness) · Kern/TEMA (Bell-Delaware tube sheet)
The MechanixCalc heat exchanger calculator sizes and rates shell-and-tube, counterflow and cross-flow exchangers using the LMTD method (with the Bowman/TEMA F-correction for 1-2 shell passes and the analytical effectiveness-NTU cross-flow F) and the NTU-effectiveness method (Incropera & DeWitt §11.4). Enter the hot and cold stream temperatures, the heat duty and the overall U-value, and the tool returns the required heat-transfer area, the F correction factor, and the terminal ΔTs in a single pass.
It is built for process, mechanical and HVAC engineers who need a defensible first-principles sizing of a heat exchanger — including a fouling resistance margin analysis, a tube sheet layout to the Bell-Delaware method, and a tube-side heat-transfer enhancement estimate using the Dittus-Boelter correlation — together with a PDF engineering report that shows the full working.
What this calculator does
- LMTD method with TEMA / Bowman F-correction for shell-and-tube (1-2 pass) and analytical cross-flow (both fluids unmixed) — no constant-F shortcuts
- NTU-effectiveness method with ε-NTU charts across five heat exchanger types (counterflow, parallel, 1-2 shell-and-tube, cross-flow, condensing/evaporating C* = 0)
- Quick HX sizing by application type (water-water, air-water, oil-water, steam condensing, air-air) with application-specific U and F defaults
- Fouling resistance analysis: cleanliness factor, fouled overall U and required area margin to TEMA fouling tables
- Tube sheet layout: Bell-Delaware tube count, baffle spacing and cross-flow area for triangular and square pitch
- Tube-side heat-transfer enhancement using the Dittus-Boelter correlation (Re, Nu, h_plain vs h_enhanced) for four enhancement types
- PDF engineering report with step-by-step methodology, governing references and design assumptions
Method & formulas
LMTD area sizing with F-correction (TEMA / Bowman)
The primary sizing method is the log-mean temperature difference (LMTD) approach. For counterflow the correction factor F = 1.0; for a 1-2 shell-and-tube exchanger the Bowman-Mueller-Nagle (1940) closed form is used, which returns null when the temperature program (P, R) lies outside the physical single-shell-pass domain (P ≥ P_max(R)) — a constant-F fallback is never used. For single-pass cross-flow with both fluids unmixed, F is computed by the effectiveness-NTU inversion method (F = NTU_counterflow / NTU_crossflow), matching the published Bowman/ESDU F-charts.
The tool warns when F < 0.75 (consider adding shell passes) and reports the required area A = Q / (U · F · LMTD). The second-law check rejects any temperature program where a terminal ΔT ≤ 0, rather than clamping it to a small positive value and reporting a spuriously large but finite area.
LMTD = (ΔT₁ − ΔT₂) / ln(ΔT₁ / ΔT₂)where ΔT₁ = T_h,in − T_c,out (counterflow) or T_h,in − T_c,in (parallel); ΔT₂ = T_h,out − T_c,in (counterflow) or T_h,out − T_c,out (parallel)
A = Q / (U · F · LMTD)where Q = heat duty (W); U = overall heat-transfer coefficient (W/m²·K); F = LMTD correction factor (≤ 1.0, arrangement-dependent); LMTD = log-mean temperature difference (K)
NTU-effectiveness method (Incropera & DeWitt §11.4)
The NTU-effectiveness (ε-NTU) method rates an existing exchanger of known area A and overall U against known stream flow rates and inlet temperatures. The number of transfer units NTU = U·A / C_min, and the effectiveness ε = Q_actual / Q_max is computed from the ε-NTU correlation for the selected flow arrangement. The maximum possible duty Q_max = C_min · (T_h,in − T_c,in) is limited by the stream with the smaller heat-capacity rate C_min = min(ṁ·c_p).
For the 1-2 shell-and-tube type, the Shah & Sekulić closed form is used; for cross-flow (both fluids unmixed) the Incropera Eq. 11.32 / ESDU 86018 series is used. The capacity ratio C* = C_min / C_max; C* = 0 represents a condensing or evaporating stream.
NTU = U · A / C_minwhere 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)
ε = Q_actual / Q_max ; Q_max = C_min · (T_h,in − T_c,in)where ε = exchanger effectiveness (0–1); Q_actual = heat transferred (W); Q_max = thermodynamic upper bound (W)
Fouling resistance and tube sheet layout
The fouling analysis panel adds the hot-side and cold-side fouling resistances R_f (m²·K/W) from the TEMA fouling tables to the clean overall resistance 1/U_clean to obtain the fouled coefficient U_fouled = 1 / (1/U_clean + R_f,hot + R_f,cold). The cleanliness factor CF = U_fouled / U_clean and the extra area margin A_margin = A_fouled − A_clean quantify the design allowance needed for fouling over the service life.
The tube sheet layout follows the Kern/TEMA tube-count formula N = CTP·(π/4)·D_s² / (CL·P_t²) for triangular (CL = 0.866) and square (CL = 1.0) pitch, de-rated by the pass-partition factor CTP (0.93 / 0.90 / 0.85 for 1 / 2 / ≥4 passes). Baffle spacing is set to 40 % of shell diameter with warnings outside the 15–45 % baffle-cut range and below TEMA's 1.25 minimum pitch ratio.
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- and cold-side fouling resistances (m²·K/W, from TEMA tables)
Worked example
Size a counterflow water-to-water shell-and-tube heat exchanger to cool hot water from 90 °C to 55 °C using cold water entering at 20 °C and leaving at 45 °C, with a heat duty of 500 kW and an overall U = 800 W/m²·K.
Given
- Hot inlet T_h,in90 °C
- Hot outlet T_h,out55 °C
- Cold inlet T_c,in20 °C
- Cold outlet T_c,out45 °C
- Heat duty Q500 kW
- Overall U800 W/m²·K
- Flow arrangementCounterflow (F = 1.0)
Result
- LMTD≈ 39.8 K
- F correction factor1.000 (counterflow)
- Required heat-transfer area A≈ 15.7 m²
- Compute the counterflow terminal temperature differences: ΔT₁ = T_h,in − T_c,out = 90 − 45 = 45 K; ΔT₂ = T_h,out − T_c,in = 55 − 20 = 35 K.
- Compute the LMTD: LMTD = (ΔT₁ − ΔT₂) / ln(ΔT₁/ΔT₂) = (45 − 35) / ln(45/35) = 10 / ln(1.2857).
- ln(1.2857) ≈ 0.2513, so LMTD = 10 / 0.2513 ≈ 39.8 K.
- For pure counterflow, F = 1.0. Apply the area equation: A = Q / (U · F · LMTD) = (500 × 1000) / (800 × 1.0 × 39.8).
- A = 500 000 / 31 840 ≈ 15.7 m².
Illustrative — verify against your actual flow rates, fluid properties and fouling allowance. Add a fouling margin: A_design = A_clean / CF, where CF = U_fouled / U_clean from the TEMA fouling tables.
Frequently asked questions
Which standard does this heat exchanger calculator use?
The LMTD method uses the TEMA F-correction charts, specifically the Bowman-Mueller-Nagle (1940) closed form for 1-2 shell-and-tube exchangers and the analytical effectiveness-NTU inversion (Incropera & DeWitt §11.4; ESDU 86018) for cross-flow — no constant-F approximations. The NTU-ε method uses the same Incropera §11.4 correlations. Tube sheet layout follows Kern/TEMA. Fouling resistances come from TEMA fouling tables. The governing method and references are shown in the PDF report.
What is the LMTD F-correction factor and when does it matter?
The F factor corrects the counterflow LMTD for non-ideal flow arrangements. For a 1-2 shell-and-tube exchanger (one shell pass, two tube passes) F < 1.0 because the streams are not in pure counterflow. When F drops below 0.75 the required area becomes very sensitive to small changes in the temperature program — the calculator warns and suggests adding shell passes or switching to counterflow. For a pure counterflow arrangement F = 1.0.
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, and the tool computes the effectiveness ε and the outlet temperatures. Both methods give the same result for a given exchanger; the choice depends on which quantities are known at the design stage.
How does the fouling resistance analysis work?
The fouling panel adds hot-side and cold-side fouling resistances (m²·K/W, selected from the built-in TEMA fouling table) to the clean overall resistance 1/U_clean to obtain the fouled overall coefficient U_fouled. The cleanliness factor CF = U_fouled / U_clean (always ≤ 1) and the extra area margin A_margin = A_fouled − A_clean show how much additional surface area must be specified to maintain the duty after the exchanger has been in service.
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. Note: some sub-panels (thermal cooling-capacity and equilibrium-temperature estimates) are engineering estimates not validated against a governing standard — they carry an in-app disclaimer.
Related calculators
- 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 HX pressure drop and system resistance.
- Thermal AnalysisSteady-state conduction resistance networks and fin-array analysis for HX design.
- Pressure Vessels (ASME VIII)Verify the shell-side and tube-side pressure containment to ASME VIII.
- CompressorPair with inter- and after-cooler sizing for compressed-air and gas systems.
Run the Heat Exchanger on your own numbers
Free 30-minute preview — no sign-up. A free 14-day account trial unlocks every tool and the branded PDF report, no credit card required.
Start free