Heat Exchanger LMTD & ε-NTU Calculator

Thermal Design and Rating of Industrial Heat Exchangers

Heat exchangers are ubiquitous thermodynamic devices engineered to transfer thermal energy between two or more fluids at different temperatures across a solid separating boundary without physical mixing of the streams. From steam surface condensers in power generation facilities and crude oil distillation preheat trains in petroleum refineries, to chilled water coils in commercial HVAC systems and plate-and-frame pasteurizers in dairy processing, thermal optimization dictates plant capital expenditure, operational efficiency, and thermodynamic irreversibility.

Thermal analysis of recuperative heat exchangers generally follows one of two primary analytical methodologies:

1. The Log Mean Temperature Difference (LMTD) Method: Primarily employed in sizing problems where process fluid flow rates and all four terminal temperatures ($T_{h,in}, T_{h,out}, T_{c,in}, T_{c,out}$) are explicitly specified by process requirements, allowing the required heat transfer surface area $A$ to be evaluated directly.
2. The Effectiveness-NTU (ε-NTU) Method: Developed by Kays and London (1984), this method is indispensable in rating and performance evaluation tasks where an existing heat exchanger of known surface area $A$ and overall heat transfer coefficient $U$ is operated under off-design inlet flow rates and temperatures, enabling explicit determination of heat duty and outlet temperatures without iterative convergence.

Heat Duty and the First Law Energy Balance

Assuming steady-state operation, negligible heat exchange with the ambient environment, negligible kinetic and potential energy changes, and constant fluid specific heat capacities ($c_p$), the conservation of energy requires that the heat released by the hot fluid must exactly equal the heat absorbed by the cold fluid:

where $C_h = dot{m}_h c_{p,h}$ and $C_c = dot{m}_c c_{p,c}$ denote the thermal heat capacity rates (in $ ext{kW/K}$ or $ ext{W/}^circ ext{C}$). The minimum and maximum heat capacity rates are identified as:

Mathematical Derivation of the Log Mean Temperature Difference (LMTD)

The local heat transfer rate across a differential area element $dA$ is governed by Newton's Law of Cooling: $$dQ = U cdot (T_h - T_c) cdot dA = U cdot Delta T cdot dA$$ Differentiating the local temperature difference $Delta T = T_h - T_c$: $$d(Delta T) = dT_h - dT_c$$ From energy balances across differential stream elements: $$dQ = -C_h dT_h = pm C_c dT_c$$ where the plus sign applies to counter-current flow (since $T_c$ increases in the opposite direction of $x$) and the minus sign applies to parallel (co-current) flow. Substituting these relations yields a separable first-order differential equation: $$rac{d(Delta T)}{Delta T} = -U left( rac{1}{C_h} pm rac{1}{C_c} ight) dA$$ Integrating over the entire exchanger area from $A = 0$ to $A_{ ext{total}}$: $$lnleft(rac{Delta T_2}{Delta T_1} ight) = -U left( rac{1}{C_h} pm rac{1}{C_c} ight) A$$ Substituting the total heat duty relations $1/C_h = (T_{h,in} - T_{h,out})/Q$ and $pm 1/C_c = (T_{c,out} - T_{c,in})/Q$ into the integrated expression yields the rigorous LMTD formula: $$Q = U cdot A cdot Delta T_{ ext{lm}}$$ $$Delta T_{ ext{lm}} = rac{Delta T_1 - Delta T_2}{ln(Delta T_1 / Delta T_2)}$$

For terminal temperature differences:

• Counter-Flow: $Delta T_1 = T_{h,in} - T_{c,out}$,   $Delta T_2 = T_{h,out} - T_{c,in}$
• Parallel-Flow: $Delta T_1 = T_{h,in} - T_{c,in}$,   $Delta T_2 = T_{h,out} - T_{c,out}$

Multipass Shell-and-Tube Geometry and the Bowman F-Factor

In practical industrial applications, pure counter-flow cannot always be maintained due to pressure drop constraints and physical equipment geometry. In shell-and-tube exchangers (TEMA designs), baffles force the shell-side fluid across the tube bundles while the tube-side fluid undergoes multiple passes (e.g., 1 shell pass, 2, 4, or 6 tube passes). This hybrid flow configuration creates localized regions of cross-flow and parallel-flow, reducing the effective mean temperature difference below that of pure counter-flow.

The actual heat transfer equation incorporates a non-dimensional correction factor $F$: $$Q = U cdot A cdot F cdot Delta T_{ ext{lm,cf}}$$ where $Delta T_{ ext{lm,cf}}$ is computed as if the streams were in pure counter-flow. The factor $F$ is calculated from thermal effectiveness $P$ and capacity ratio $R$: $$P = rac{T_{c,out} - T_{c,in}}{T_{h,in} - T_{c,in}}, quad R = rac{T_{h,in} - T_{h,out}}{T_{c,out} - T_{c,in}} = rac{C_c}{C_h}$$ For a standard 1-shell pass, even-number tube pass exchanger, the Bowman-Mueller-Nagle equation provides the exact analytical solution: $$F = rac{sqrt{R^2 + 1} cdot lnleft( rac{1 - P}{1 - P R} ight)}{(R - 1) cdot lnleft( rac{2 - P(R + 1 - sqrt{R^2 + 1})}{2 - P(R + 1 + sqrt{R^2 + 1})} ight)}$$ Design Rule: Process engineering codes (such as TEMA) require $F ge 0.80$. If $F < 0.75$, the design experiences severe thermal crossing, requiring two or more shells in series.

The Effectiveness-NTU (ε-NTU) Relations

Thermal effectiveness $epsilon$ is defined as the ratio of actual heat duty $Q$ to the thermodynamically maximum possible heat transfer $Q_{max}$ that could be achieved in a counter-flow exchanger of infinite surface area: $$Q_{max} = C_{min} cdot (T_{h,in} - T_{c,in})$$ $$epsilon = rac{Q}{Q_{max}} = rac{C_h (T_{h,in} - T_{h,out})}{C_{min}(T_{h,in} - T_{c,in})} = rac{C_c (T_{c,out} - T_{c,in})}{C_{min}(T_{h,in} - T_{c,in})}$$ The Number of Transfer Units ($ ext{NTU}$) represents a non-dimensional measure of the heat exchanger's physical thermal size: $$ ext{NTU} = rac{U cdot A}{C_{min}}$$ Analytical $epsilon$-NTU relationships for standard flow configurations:

• Pure Counter-Flow ($C_r < 1$): $$epsilon = rac{1 - exp(- ext{NTU}(1 - C_r))}{1 - C_r cdot exp(- ext{NTU}(1 - C_r))}, quad ext{for } C_r = 1: quad epsilon = rac{ ext{NTU}}{1 + ext{NTU}}$$
• Parallel-Flow: $$epsilon = rac{1 - exp(- ext{NTU}(1 + C_r))}{1 + C_r}$$
• Cross-Flow (Both Streams Unmixed): $$epsilon approx 1 - expleft( rac{ ext{NTU}^{0.22}}{C_r} cdot left[ expleft( -C_r cdot ext{NTU}^{0.78} ight) - 1 ight] ight)$$
• Boiling or Condensation ($C_{max} o infty, C_r = 0$): $$epsilon = 1 - exp(- ext{NTU})$$

Step-by-Step Worked Numerical Examples

Example 1: Sizing a Water-to-Water Plate Heat Exchanger

Given: Hot process water flows at $dot{m}_h = 3.0 ext{ kg/s}$ entering at $T_{h,in} = 90.0^circ ext{C}$ and exiting at $T_{h,out} = 60.0^circ ext{C}$. Cold cooling tower water enters at $T_{c,in} = 25.0^circ ext{C}$ with a mass flow rate of $dot{m}_c = 3.5 ext{ kg/s}$. Specific heats are $c_{p,h} = c_{p,c} = 4.184 ext{ kJ/(kg}cdot ext{K)}$. The estimated overall heat transfer coefficient is $U = 2,200 ext{ W/(m}^2cdot ext{K)}$.

Step 1: Compute Heat Duty ($Q$):

Step 2: Determine Cold Fluid Outlet Temperature ($T_{c,out}$):

Step 3: Compute Terminal Temperature Differences for Counter-Flow:

Step 4: Compute Counter-Flow LMTD:

Step 5: Compute Required Heat Transfer Surface Area ($A$):

Example 2: Rating an Existing Lube Oil Cooler

Scenario: A marine diesel engine lube oil cooler has a fixed surface area $A = 12.0 ext{ m}^2$ and $U = 450 ext{ W/(m}^2cdot ext{K)}$. Hot oil enters at $T_{h,in} = 85.0^circ ext{C}$ with $dot{m}_h = 1.8 ext{ kg/s}$ ($c_{p,h} = 2.05 ext{ kJ/(kg}cdot ext{K)}$). Cooling seawater enters at $T_{c,in} = 18.0^circ ext{C}$ with $dot{m}_c = 2.5 ext{ kg/s}$ ($c_{p,c} = 3.98 ext{ kJ/(kg}cdot ext{K)}$). Find the oil outlet temperature.

Calculations:
1. Heat capacity rates: $C_h = 1.8 imes 2.05 = 3.69 ext{ kW/K} = 3690 ext{ W/K}$.
$C_c = 2.5 imes 3.98 = 9.95 ext{ kW/K} = 9950 ext{ W/K}$.
2. Minimum and maximum capacity: $C_{min} = 3.69 ext{ kW/K}$ (oil is the minimum fluid), $C_{max} = 9.95 ext{ kW/K}$.
Capacity ratio $C_r = 3.69 / 9.95 = 0.3709$.
3. NTU calculation: $$ ext{NTU} = rac{U cdot A}{C_{min}} = rac{450 imes 12.0}{3690} = rac{5400}{3690} = 1.463$$
4. Counter-flow effectiveness: $$epsilon = rac{1 - exp(-1.463 imes (1 - 0.3709))}{1 - 0.3709 cdot exp(-1.463 imes (1 - 0.3709))} = rac{1 - exp(-0.9204)}{1 - 0.3709 imes 0.3984} = rac{0.6016}{0.8522} = 0.7059 quad (70.59%)$$
5. Maximum possible heat transfer: $$Q_{max} = C_{min} cdot (T_{h,in} - T_{c,in}) = 3.69 imes (85.0 - 18.0) = 3.69 imes 67.0 = 247.23 ext{ kW}$$
6. Actual heat transfer: $Q = epsilon imes Q_{max} = 0.7059 imes 247.23 = 174.52 ext{ kW}$.
7. Oil outlet temperature: $$T_{h,out} = T_{h,in} - rac{Q}{C_h} = 85.0 - rac{174.52}{3.69} = 85.0 - 47.30 = 37.7^circ ext{C}$$

Comparative Analysis of Heat Exchanger Types

Typical Fouling Resistances ($R_f$) in Process Design

Engineering FAQs