NTU method explained

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

Th,i-Tc,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 fluid streams, denoted as

cp

. By definition

cp

is the derivative of enthalpy with respect to temperature:
c
p=dh
dT
This information can usually be found in a thermodynamics textbook,[1] or by using various software packages. Additionally, the mass flowrates (
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)

Ch

and

Cc

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

Cmin

:

Cmin=min[

m

ccp,c,

m

hcp,h]

Where

m
is the mass flow rate and

cp

is the fluid's specific heat capacity at constant pressure. The maximum possible heat transfer rate is then determined by the following expression:
Q

max =Cmin(Th,i-Tc,i)

Here,

Q

max

is the maximum rate of heat that could be transferred between the fluids per unit time.

Cmin

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 (

\epsilon

), is the ratio between the actual heat transfer rate and the maximum possible heat transfer rate:

\epsilon =

Q
Qmax

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

Q

 =Ch(Th,i-Th,o) =Cc(Tc,o-Tc,i)

Effectiveness is a dimensionless quantity between 0 and 1. If we know

\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:
Q

 =\epsilonCmin(Th,i-Tc,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:

\epsilon=f(NTU,

Cmin
Cmax

)

For a given geometry,

\epsilon

can be calculated using correlations in terms of the "heat capacity ratio," or

Cr

and NTU:

Cr =

Cmin
Cmax

NTU

describes heat transfer across a surface

NTU =

UA
Cmin
Here,

U

is the overall heat transfer coefficient,

A

is the total heat transfer area, and

Cmin

is the minimum heat capacity rate. To better understand where this definition of NTU comes from, consider the following heat transfer energy balance, which is an extension of the energy balance above:
Q=

Cmin(To-Ti)min=UA\DeltaTLM

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 (

\DeltaTLM

). Starting from the differential equations that describe heat transfer, several "simple" correlations between effectiveness and NTU can be made. 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:

\epsilon =

1-\exp[-NTU(1+Cr)]
1+Cr

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

\epsilon =

1-\exp[-NTU(1-Cr)]
1-Cr\exp[-NTU(1-Cr)]

For a balanced counter-current flow heat exchanger (balanced meaning

Cr =1

, which is a scenario desirable to enable irreversible entropy production to be reduced given sufficient heat transfer area):

\epsilon =

NTU
1+NTU

A single-stream heat exchanger is a special case in which

Cr =0

. This occurs when

Cmin=0

or

Cmax=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:[2]

\epsilon =1-e-NTU

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

\epsilon =1-\exp(-NTU)-\exp[-(1+Cr)NTU]

infty
\sum
n=1
n
C
r

Pn(NTU)

where

Pn

is the polynomial function

Pn(x)=

1
(n+1)!
n
\sum
j=1
n+1-j
j!

xn+j

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

\epsilon =\left[

1
1-\exp(-NTU)

+

Cr
1-\exp(-NTUCr)

-

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

Cmin

corresponds to the mixed fluid, the result is

\epsilon =1-\exp\left(-

1-\exp(-NTUCr)
Cr

\right)

whereas if

Cmin

corresponds to the unmixed fluid, the solution is

\epsilon =

1
Cr

(1-\exp\{-Cr[1-\exp(-NTU)]\})

All these formulas for crossflow heat exchangers are also valid for

Cr=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.[3]

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.[4] 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:

c
p=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:

cp-x=

d\omegax
dPx
Here,

\omegax

represents the mass ratio of gas 'x' (meaning mass of gas 'x' relative to the mass of all other non-'x' gas mass) and

Px

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:

cp-x=

Mx/Mother
Pother
Here,

Mx

is the molecular weight of gas 'x' and

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:

NTUx =

UmAm
mcp-x
Here,

Um

is the overall mass transfer coefficient, which could be determined by empirical correlations,

Am

is the surface area for mass transfer (particularly relevant in membrane-based separations), and
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

Cr

becomes:

Cr =

(mcp-x)min
(mcp-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."

cp-h=

Mwv/Mair=
Pair
0.62198=
Pair
0.62198
Ptotal-Pwv,inlet
Here,

Mwv

is the molecular weight of water (vapor),

Mair

is the average molecular weight of air,

Pair

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,

Pwv,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

Notes and References

  1. Book: Çengel, Yunus A. . Thermodynamics: an engineering approach . Boles . Michael A. . 2015 . McGraw-Hill Education . 978-0-07-339817-4 . 8. edition in SI units . New York, NY.
  2. Book: J. H. Lienhard IV . J. H. Lienhard V. August 14, 2020. A Heat Transfer Textbook. Phlogiston Press. 121–127.
  3. Book: Shah, Ramesh . Fundamentals of Heat Exchanger Design . Sekulic . Dusan . 2003 . John Wiley & Sons . 0-471-32171-0 . 1 . Hoboken, NJ.
  4. 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 . 10.1016/j.applthermaleng.2023.121514 . 1359-4311.
  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 . 10.1016/j.enconman.2022.116491 . 0196-8904. free .