Thermal Process Engineering & Heat Exchanger Fundamentals

Heat exchangers are thermodynamic devices engineered to transfer thermal energy between two or more fluid streams across a solid conduction boundary without intermixing. Operating in power stations, chemical refineries, district cooling plants, food pasteurization lines, and refrigeration chillers, heat exchangers are governed by the Tubular Exchanger Manufacturers Association (TEMA) and ASME Section VIII standards.

Proper sizing requires solving the combined effects of convective fluid boundary layers, conductive wall resistance, fluid fouling scale, and temperature profile convergence. Designing an exchanger with insufficient area fails to achieve process target temperatures, causing boiling instability or product ruination. Conversely, gross oversizing causes fluid velocity stagnation, accelerating fouling deposition, sedimentation, and tube corrosion.

Thermal Heat Duty Formulation

Under steady-state operation neglecting atmospheric heat leakage, the total thermal duty (\(Q\)) transferred from the hot stream to the cold stream is governed by the First Law of Thermodynamics:

$$Q = \dot{m}_h c_{p,h} (T_{h,\text{in}} - T_{h,\text{out}}) = \dot{m}_c c_{p,c} (T_{c,\text{out}} - T_{c,\text{in}}) \quad [\text{kW or BTU/hr}]$$

Where:

  • \(\dot{m}_h, \dot{m}_c\): Mass flow rates of the hot and cold streams (\(\text{kg/s}\) or \(\text{lbs/hr}\)).
  • \(c_{p,h}, c_{p,c}\): Specific heat capacities at constant pressure (\(\text{kJ/kg}\cdot\text{K}\) or \(\text{BTU/lb}\cdot^\circ\text{F}\)).
  • \(C_h = \dot{m}_h c_{p,h}\) and \(C_c = \dot{m}_c c_{p,c}\): Fluid heat capacity rates (\(\text{kW/K}\)).

The Logarithmic Mean Temperature Difference (LMTD) Method

Because fluid temperatures vary continuously along the tube bundle, arithmetic temperature differences (\(\Delta T_{\text{avg}} = (T_{h} - T_{c})\)) grossly overestimate heat transfer. The mathematically rigorous driving gradient is the Logarithmic Mean Temperature Difference (LMTD):

$$\text{LMTD} = \frac{\Delta T_1 - \Delta T_2}{\ln(\Delta T_1 / \Delta T_2)}$$

Counter-Current Flow Boundary Definitions

In a counter-flow heat exchanger, the fluids enter at opposite ends and travel in reverse directions:

$$\Delta T_1 = T_{h,\text{in}} - T_{c,\text{out}} \quad \text{and} \quad \Delta T_2 = T_{h,\text{out}} - T_{c,\text{in}}$$

Co-Current (Parallel) Flow Boundary Definitions

In parallel flow, both fluids enter at the same end and travel in parallel:

$$\Delta T_1 = T_{h,\text{in}} - T_{c,\text{in}} \quad \text{and} \quad \Delta T_2 = T_{h,\text{out}} - T_{c,\text{out}}$$

Because counter-flow maintains an almost uniform temperature difference across the entire heat exchange length, its LMTD is consistently 20% to 50% higher than parallel flow, enabling a drastically smaller required heat transfer surface area.

Shell & Tube LMTD Correction Factor (F)

When fluid flow is not purely counter-current, as occurs in multi-pass shell-and-tube exchangers (TEMA E-shell with 2, 4, or 6 tube passes) and cross-flow radiators, the effective mean temperature difference (\(\Delta T_m\)) is derated by an empirical geometry correction factor (\(F \le 1.0\)):

$$\Delta T_m = F \cdot \text{LMTD}_{\text{counter}}$$

The correction factor \(F\) is calculated as a function of the heat capacity ratio (\(R\)) and thermal effectiveness (\(P\)):

$$R = \frac{T_{h,\text{in}} - T_{h,\text{out}}}{T_{c,\text{out}} - T_{c,\text{in}}} = \frac{C_c}{C_h}, \quad P = \frac{T_{c,\text{out}} - T_{c,\text{in}}}{T_{h,\text{in}} - T_{c,\text{in}}}$$ $$F = \frac{\sqrt{R^2 + 1} \cdot \ln\left(\frac{1 - P}{1 - P \cdot R} \right)}{(R - 1) \cdot \ln\left(\frac{2 - P(R + 1 - \sqrt{R^2 + 1})}{2 - P(R + 1 + \sqrt{R^2 + 1})} \right)}$$

TEMA standards recommend that exchangers never be operated at \(F < 0.75\), as the curve drops precipitously toward zero, making the exchanger unstable and susceptible to temperature pinching.

Required Heat Transfer Surface Area & Overall U-Value

The total outside surface area (\(A\)) required to transfer heat duty \(Q\) is:

$$A = \frac{Q}{U_{\text{dirty}} \cdot \Delta T_m} \quad [\text{m}^2 \text{ or ft}^2]$$

Where \(U_{\text{dirty}}\) incorporates both internal and external convective film coefficients (\(h_i, h_o\)), tube wall metal thermal conduction (\(k_{\text{metal}}\)), and fouling scale resistance factors (\(R_{fi}, R_{fo}\)):

$$\frac{1}{U_{\text{dirty}}} = \frac{1}{U_{\text{clean}}} + R_f$$
Service & Fluid Combination Clean Coeff \(U\) (\(\text{W/m}^2\cdot\text{K}\)) Typical Fouling \(R_f\) (\(\text{m}^2\cdot\text{K/W}\)) Design Coeff \(U_{\text{dirty}}\) (\(\text{W/m}^2\cdot\text{K}\))
Water to Water (Plate & Frame) 3,000 to 7,000 0.00005 to 0.00010 2,000 to 4,500
Water to Water (Shell & Tube) 1,200 to 2,500 0.00015 to 0.00025 800 to 1,500
Steam Condenser (Water in Tubes) 2,000 to 4,000 0.00010 to 0.00020 1,500 to 2,800
Water to Lube Oil (Tubular Cooler) 250 to 500 0.00020 to 0.00035 180 to 350
Air to Water (Fin-Fan Radiator) 30 to 60 0.00020 to 0.00040 25 to 50

Step-by-Step Worked Case Study: District Heating Substation Water-Water Cooler

Design Scenario: An engineer is sizing a counter-flow shell-and-tube heat exchanger for a commercial building district heating substation:

  • Hot district water: \(T_{h,\text{in}} = 95^\circ\text{C}\), \(T_{h,\text{out}} = 65^\circ\text{C}\).
  • Hot fluid mass flow: \(\dot{m}_h = 2.0\text{ kg/s}\), \(c_{p,h} = 4.18\text{ kJ/kg}\cdot\text{K}\).
  • Cold building loop: \(T_{c,\text{in}} = 20^\circ\text{C}\), \(T_{c,\text{out}} = 45^\circ\text{C}\).
  • Overall clean heat transfer coefficient: \(U_{\text{clean}} = 850\text{ W/m}^2\cdot\text{K}\).
  • Treated water fouling resistance: \(R_f = 0.0001\text{ m}^2\cdot\text{K/W}\).

Step 1: Compute thermal heat duty (\(Q\))

$$Q = \dot{m}_h c_{p,h} (T_{h,\text{in}} - T_{h,\text{out}}) = 2.0 \cdot 4.18 \cdot (95 - 65) = 8.36 \times 30 = 250.8\text{ kW}$$

Step 2: Determine Counter-Flow LMTD

$$\Delta T_1 = T_{h,\text{in}} - T_{c,\text{out}} = 95 - 45 = 50.0^\circ\text{C}$$ $$\Delta T_2 = T_{h,\text{out}} - T_{c,\text{in}} = 65 - 20 = 45.0^\circ\text{C}$$ $$\text{LMTD} = \frac{50.0 - 45.0}{\ln(50.0 / 45.0)} = \frac{5.0}{\ln(1.1111)} = \frac{5.0}{0.10536} = 47.46^\circ\text{C}$$

Step 3: Calculate derated fouled overall coefficient (\(U_{\text{dirty}}\))

$$\frac{1}{U_{\text{dirty}}} = \frac{1}{850} + 0.0001 = 0.001176 + 0.000100 = 0.001276 \implies U_{\text{dirty}} = 783.4\text{ W/m}^2\cdot\text{K}$$

Step 4: Compute required heat transfer surface area (\(A\))

$$A = \frac{Q}{U_{\text{dirty}} \cdot \text{LMTD}} = \frac{250,800\text{ W}}{783.4\text{ W/m}^2\cdot\text{K} \cdot 47.46\text{ K}} = \frac{250800}{37180.2} = 6.75\text{ m}^2$$

Engineering Conclusion: A counter-flow exchanger with \(6.75\text{ m}^2\) of surface area provides the required 251 kW heating duty with built-in fouling reserve.

Frequently Asked Questions

What is a temperature pinch point in heat exchanger design?

A temperature pinch point occurs where the temperature difference between the hot and cold fluid streams approaches a minimal value (\(\Delta T_{\text{min}} \to 0\)). As the pinch point narrows, the LMTD approaches zero, causing the required surface area to approach infinity. In industrial process design, minimum approach temperatures are typically restricted to \(3^\circ\text{C}\) to \(5^\circ\text{C}\) for plate heat exchangers and \(5^\circ\text{C}\) to \(10^\circ\text{C}\) for shell-and-tube units.

Why are Plate and Frame heat exchangers more compact than Shell and Tube?

Gasketed plate heat exchangers feature thin, pressed metal plates with corrugation chevrons that induce intense boundary-layer turbulence even at low Reynolds numbers (\(Re < 50\)). This produces heat transfer coefficients 3 to 5 times higher than smooth cylindrical tubes (\(U \approx 4,000 - 6,000\text{ W/m}^2\cdot\text{K}\)), requiring 70% less surface area and foot-print space.

How does baffle spacing affect shell-side heat transfer and pressure drop?

Segmental baffles in shell-and-tube exchangers direct fluid across the tube bundle perpendicular to the tubes, generating high cross-flow velocities that boost heat transfer. However, closing baffle spacing increases shell-side pressure drop exponentially (\(\Delta P \propto v^2\)). TEMA standards require baffle spacing to be maintained between 20% and 100% of the inside shell diameter.

When must the Effectiveness-NTU method be used instead of LMTD?

The \(\varepsilon\)-NTU method is essential when evaluating an existing heat exchanger under new operating flow rates or inlet temperatures where the outlet temperatures are unknown. Solving for outlet temperatures using LMTD requires cumbersome, iterative numerical root-finding because LMTD depends directly on the unknown outlet temperatures. In contrast, the \(\varepsilon\)-NTU method calculates effectiveness directly from known heat capacity rates and area.