NTU method

The number of transfer units (NTU) method is used to calculate the rate of heat transfer in heat exchangers (especially parallel flow, counter current, and cross-flow exchangers) when there is insufficient information to calculate the log mean temperature difference (LMTD). Alternatively, this method is useful for determining the expected heat exchanger effectiveness from the known geometry. In heat exchanger analysis, if the fluid inlet and outlet temperatures are specified or can be determined by simple energy balance, the LMTD method can be used; but when these temperatures are not available either the NTU or the effectiveness NTU method is used.

The effectiveness-NTU method is very useful for all the flow arrangements (besides parallel flow, cross flow, and counterflow ones) but the effectiveness of all other types must be obtained by a numerical solution of the partial differential equations and there is no analytical equation for LMTD or effectiveness.

Defining and using heat exchanger effectiveness

To define the effectiveness of a heat exchanger we need to find the maximum possible heat transfer that can be hypothetically achieved in a counter-flow heat exchanger of infinite length. Therefore one fluid will experience the maximum possible temperature difference, which is the difference of ${\displaystyle \ T_{h,i}-\ T_{c,i}}$ (the temperature difference between the inlet temperature of the hot stream and the inlet temperature of the cold stream). First, you must know the specific heat capacity of your two air streams, denoted as ${\displaystyle c_{p}}$. By definition ${\displaystyle c_{p}}$ is the derivative of enthalpy with respect to temperature:

${\displaystyle c_{p}={\frac {dh}{dT}}}$

This information can usually be found in a thermodynamics textbook, ^[1] or by using various software packages. Additionally, the mass flowrates (${\displaystyle {\dot {m}}}$) of the two streams exchanging heat must be known (here, the cold stream is denoted with subscripts 'c' and the hot stream is denoted with subscripts 'h'). The method proceeds by calculating the heat capacity rates (i.e. mass flow rate multiplied by specific heat capacity) ${\displaystyle \ C_{h}}$ and ${\displaystyle \ C_{c}}$ for the hot and cold fluids respectively. To determine the maximum possible heat transfer rate in the heat exchanger, the minimum heat capacity rate must be used, denoted as ${\displaystyle \ C_{\mathrm {min} }}$:

${\displaystyle \ C_{\mathrm {min} }=\mathrm {min} [{\dot {m}}_{c}c_{p,c},{\dot {m}}_{h}c_{p,h}]}$

Where ${\displaystyle {\dot {m}}}$ is the mass flow rate and ${\displaystyle c_{p}}$ is the fluid's specific heat capacity at constant pressure. The maximum possible heat transfer rate is then determined by the following expression:

${\displaystyle {\dot {Q}}_{\mathrm {max} }\ =C_{\mathrm {min} }(T_{h,i}-T_{c,i})}$

Here, ${\displaystyle {\dot {Q}}_{\mathrm {max} }}$ is the maximum rate of heat that could be transferred between the fluids per unit time. ${\displaystyle \ C_{\mathrm {min} }}$ must be used as it is the fluid with the lowest heat capacity rate that would, in this hypothetical infinite length exchanger, actually undergo the maximum possible temperature change. The other fluid would change temperature more slowly along the heat exchanger length. The method, at this point, is concerned only with the fluid undergoing the maximum temperature change.

The effectiveness of the heat exchanger (${\displaystyle \epsilon }$), is the ratio between the actual heat transfer rate and the maximum possible heat transfer rate:

${\displaystyle \epsilon \ ={\frac {\dot {Q}}{{\dot {Q}}_{\mathrm {max} }}}}$

where the real heat transfer rate can be determined either from the cold fluid or the hot fluid (they must provide equivalent results):

${\displaystyle {\dot {Q}}\ =C_{h}(T_{h,i}-T_{h,o})\ =C_{c}(T_{c,o}-T_{c,i})}$

Effectiveness is a dimensionless quantity between 0 and 1. If we know ${\displaystyle \epsilon }$ for a particular heat exchanger, and we know the inlet conditions of the two flow streams we can calculate the amount of heat being transferred between the fluids by:

${\displaystyle {\dot {Q}}\ =\epsilon C_{\mathrm {min} }(T_{h,i}-T_{c,i})}$

Then, having determined the actual heat transfer from the effectiveness and inlet temperatures, the outlet temperatures can be determined from the equation above.

Relating Effectiveness to the Number of Transfer Units (NTU)

For any heat exchanger it can be shown that the effectiveness of the heat exchanger is related to a non-dimensional term called the "number of transfer units" or NTU:

${\displaystyle \ \epsilon =f(NTU,{\frac {C_{\mathrm {min} }}{C_{\mathrm {max} }}})}$

For a given geometry, ${\displaystyle \epsilon }$ can be calculated using correlations in terms of the "heat capacity ratio," or ${\displaystyle C_{r}}$ and NTU:

${\displaystyle C_{r}\ ={\frac {C_{\mathrm {min} }}{C_{\mathrm {max} }}}}$

${\displaystyle \ NTU}$ describes heat transfer across a surface

${\displaystyle NTU\ ={\frac {UA}{C_{\mathrm {min} }}}}$

Here, ${\displaystyle U}$ is the overall heat transfer coefficient, ${\displaystyle A}$ is the total heat transfer area, and ${\displaystyle C_{min}}$ is the minimum heat capacity rate. To better understand where this definition of NTU comes from, consider the follow heat transfer energy balance, which is an extension of the energy balance above:

${\displaystyle {\dot {Q}}=C_{min}(T_{o}-T_{i})_{min}=UA\Delta T_{LM}}$

From this energy balance, it is clear that NTU relates the temperature change of the flow with the minimum heat capacitance rate to the log mean temperature difference (${\displaystyle \Delta T_{LM}}$). Starting from the differential equations that describe heat transfer, several "simple" correlations between effectiveness and NTU can be made. ^[2] For brevity, below summarizes the Effectiveness-NTU correlations for some of the most common flow configurations:

For example, the effectiveness of a parallel flow heat exchanger is calculated with:

${\displaystyle \epsilon \ ={\frac {1-\exp[-NTU(1+C_{r})]}{1+C_{r}}}}$

Or the effectiveness of a counter-current flow heat exchanger is calculated with:

${\displaystyle \epsilon \ ={\frac {1-\exp[-NTU(1-C_{r})]}{1-C_{r}\exp[-NTU(1-C_{r})]}}}$

For a balanced counter-current flow heat exchanger (balanced meaning ${\displaystyle C_{r}\ =1}$, which is a scenario desirable to reduce entropy):

${\displaystyle \epsilon \ ={\frac {NTU}{1+NTU}}}$

A single-stream heat exchanger is a special case in which ${\displaystyle C_{r}\ =0}$. This occurs when ${\displaystyle C_{\mathrm {min} }=0}$ or ${\displaystyle C_{\mathrm {max} }=\infty }$ and may represent a situation in which a phase change (condensation or evaporation) is occurring in one of the heat exchanger fluids or when one of the heat exchanger fluids is being held at a fixed temperature. In this special case the heat exchanger behavior is independent of the flow arrangement and the effectiveness is given by: ^[3]

${\displaystyle \epsilon \ =1-e^{-NTU}}$

For a crossflow heat exchanger with both fluid unmixed, the effectiveness is:

${\displaystyle \epsilon \ =1-\exp(-NTU)-\exp[-(1+C_{r})NTU]\sum _{n=1}^{\infty }C_{r}^{n}P_{n}(NTU)}$

where ${\displaystyle P_{n}}$ is the polynomial function

${\displaystyle P_{n}(x)={\frac {1}{(n+1)!}}\sum _{j=1}^{n}{\frac {n+1-j}{j!}}x^{n+j}}$

If both fluids are mixed in the crossflow heat exchanger, then

${\displaystyle \epsilon \ =\left[{\frac {1}{1-\exp(-NTU)}}+{\frac {C_{r}}{1-\exp(-NTU\cdot C_{r})}}-{\frac {1}{NTU}}\right]^{-1}}$

If one of the fluids in the crossflow heat exchanger is mixed and the other is unmixed, the result depends on which one has the minimum heat capacity rate. If ${\displaystyle C_{\mathrm {min} }}$ corresponds to the mixed fluid, the result is

${\displaystyle \epsilon \ =1-\exp \left(-{\frac {1-\exp(-NTU\cdot C_{r})}{C_{r}}}\right)}$

whereas if ${\displaystyle C_{\mathrm {min} }}$ corresponds to the unmixed fluid, the solution is

${\displaystyle \epsilon \ ={\frac {1}{C_{r}}}(1-\exp\{-C_{r}[1-\exp(-NTU)]\})}$

All these formulas for crossflow heat exchangers are also valid for ${\displaystyle C_{r}=1}$.

Additional effectiveness-NTU analytical relationships have been derived for other flow arrangements, including shell-and-tube heat exchangers with multiple passes and different shell types, and plate heat exchangers. ^[4]

Effectiveness-NTU method for gaseous mass transfer

It is common in the field of mass transfer system design and modeling to draw analogies between heat transfer and mass transfer. ^[2] However, a mass transfer-analogous definition of the effectiveness-NTU method requires some additional terms. One common misconception is that gaseous mass transfer is driven by concentration gradients, however, in reality it is the partial pressure of the given gas that drive mass transfer. In the same way that the heat transfer definition includes the specific heat capacity of the fluid, which describes the change in enthalpy of the fluid with respect to change in temperature and is defined as:

${\displaystyle c_{p}={\frac {dh}{dT}}}$

then a mass transfer-analogous specific mass capacity is required. This specific mass capacity should describe the change in concentration of the transferring gas relative to the partial pressure difference driving the mass transfer. This results in a definition for specific mass capacity as follows:

${\displaystyle c_{p-x}={\frac {d\omega _{x}}{dP_{x}}}}$

Here, ${\displaystyle \omega _{x}}$ represents the mass ratio of gas 'x' (meaning mass of gas 'x' relative to the mass of all other non-'x' gas mass) and ${\displaystyle P_{x}}$ is the partial pressure of gas 'x'. Using the ideal gas formulation for the mass ratio gives the following definition for the specific mass capacity:

${\displaystyle c_{p-x}={\frac {M_{x}/M_{other}}{P_{other}}}}$

Here, ${\displaystyle M_{x}}$ is the molecular weight of gas 'x' and ${\displaystyle M_{other}}$ is the average molecular weight of all other gas constituents. With this information, the NTU for gaseous mass transfer of gas 'x' can be defined as follows:

${\displaystyle NTU_{x}\ ={\frac {U_{m}A_{m}}{{\dot {m}}c_{p-x}}}}$

Here, ${\displaystyle U_{m}}$ is the overall mass transfer coefficient, which could be determined by empirical correlations, ${\displaystyle A_{m}}$ is the surface area for mass transfer (particularly relevant in membrane-based separations), and ${\displaystyle {\dot {m}}}$ is the mass flowrate of bulk fluid (e.g., mass flowrate of air in an application where water vapor is being separated from the air mixture). At this point, all of the same heat transfer effectiveness-NTU correlations will accurately predict the mass transfer performance, as long as the heat transfer terms in the definition of NTU have been replaced by the mass transfer terms, as shown above. Similarly, it follows that the definition of ${\displaystyle C_{r}}$ becomes:

${\displaystyle C_{r}\ ={\frac {({\dot {m}}c_{p-x})_{min}}{({\dot {m}}c_{p-x})_{max}}}}$

Effectiveness-NTU method for dehumidification applications

One particularly useful application for the above described effectiveness-NTU framework is membrane-based air dehumidification. ^[5] In this case, the definition of specific mass capacity can be defined for humid air and is termed "specific humidity capacity." ^[2]

${\displaystyle c_{p-h}={\frac {M_{wv}/M_{air}}{P_{air}}}={\frac {0.62198}{P_{air}}}={\frac {0.62198}{P_{total}-P_{wv,inlet}}}}$

Here, ${\displaystyle M_{wv}}$ is the molecular weight of water (vapor), ${\displaystyle M_{air}}$ is the average molecular weight of air, ${\displaystyle P_{air}}$ is the partial pressure of air (not including the partial pressure of water vapor in an air mixture) and can be approximated by knowing the partial pressure of water vapor at the inlet, before dehumidification occurs, ${\displaystyle P_{wv,inlet}}$. From here, all of the previously described equations can be used to determine the effectiveness of the mass exchanger.

Importance of defining the specific mass capacity

It is very common, especially in dehumidification applications, to define the mass transfer driving force as the concentration difference. When deriving effectiveness-NTU correlations for membrane-based gas separations, this is valid only if the total pressures are approximately equal on both sides of the membrane (e.g., an energy recovery ventilator for a building). This is sufficient since the partial pressure and concentration are proportional. However, if the total pressures are not approximately equal on both sides of the membrane, the low pressure side could have a higher "concentration" but a lower partial pressure of the given gas (e.g., water vapor in a dehumidification application) than the high pressure side, thus using the concentration as the driving is not physically accurate.

References

1. Çengel, Yunus A.; Boles, Michael A. (2015). Thermodynamics: an engineering approach (8. edition in SI units ed.). New York, NY: McGraw-Hill Education. ISBN   978-0-07-339817-4.
2. Fix, Andrew J.; Braun, James E.; Warsinger, David M. (2024-01-05). "A general effectiveness-NTU modeling framework for membrane dehumidification systems". Applied Thermal Engineering. 236: 121514. doi:10.1016/j.applthermaleng.2023.121514. ISSN   1359-4311.
3. J. H. Lienhard IV; J. H. Lienhard V (August 14, 2020). A Heat Transfer Textbook. Phlogiston Press. pp. 121–127.
4. Shah, Ramesh; Sekulic, Dusan (2003). Fundamentals of Heat Exchanger Design (1 ed.). Hoboken, NJ: John Wiley & Sons. ISBN   0-471-32171-0.
5. Fix, Andrew J.; Gupta, Shivam; Braun, James E.; Warsinger, David M. (2023-01-15). "Demonstrating non-isothermal vacuum membrane air dehumidification for efficient next-generation air conditioning". Energy Conversion and Management. 276: 116491. doi: 10.1016/j.enconman.2022.116491 . ISSN   0196-8904.

• Kays & London, 1955 Compact Heat Exchangers
• F. P. Incropera & D. P. DeWitt 1990 Fundamentals of Heat and Mass Transfer, 3rd edition, pp. 658–660. Wiley, New York
• Incropera, F. P.; DeWitt, D. P.; Bergman, T. L.; Lavine, A. S. (2006). Fundamentals of Heat and Mass Transfer (6th ed.). John Wiley & Sons US. pp. 686–688. ISBN   978-0471457282.