Physics:NTU method

From HandWiki
Short description: Method to calculate rate of heat transfer in heat exchangers

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 [math]\displaystyle{ \ T_{h,i}- \ T_{c,i} }[/math] (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 [math]\displaystyle{ c_p }[/math]. By definition [math]\displaystyle{ c_p }[/math] is the derivative of enthalpy with respect to temperature:

[math]\displaystyle{ c_p=\frac{dh}{dT} }[/math]

This information can usually be found in a thermodynamics textbook,[1] or by using various software packages. Additionally, the mass flowrates ([math]\displaystyle{ \dot m }[/math]) 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) [math]\displaystyle{ \ C_h }[/math] and [math]\displaystyle{ \ C_c }[/math] 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 [math]\displaystyle{ \ C_\mathrm{min} }[/math]:

[math]\displaystyle{ \ C_\mathrm{min}=\mathrm{min}[ \dot m_c c_{p,c}, \dot m_h c_{p,h}] }[/math]

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

[math]\displaystyle{ \dot Q_\mathrm{max}\ = C_\mathrm{min} (T_{h,i}-T_{c,i}) }[/math]

Here, [math]\displaystyle{ \dot Q_\mathrm{max} }[/math] is the maximum rate of heat that could be transferred between the fluids per unit time. [math]\displaystyle{ \ C_\mathrm{min} }[/math] 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 ([math]\displaystyle{ \epsilon }[/math]), is the ratio between the actual heat transfer rate and the maximum possible heat transfer rate:

[math]\displaystyle{ \epsilon \ = \frac{\dot Q}{\dot Q _\mathrm{max}} }[/math]

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

[math]\displaystyle{ \dot Q \ = C_h (T_{h,i} -T_{h,o})\ = C_c (T_{c,o} - T_{c,i}) }[/math]

Effectiveness is a dimensionless quantity between 0 and 1. If we know [math]\displaystyle{ \epsilon }[/math] 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:

[math]\displaystyle{ \dot Q \ = \epsilon C_\mathrm{min} (T_{h,i} -T_{c,i}) }[/math]

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:

[math]\displaystyle{ \ \epsilon = f ( NTU,\frac{C_\mathrm{min}} {C_\mathrm{max}}) }[/math]

For a given geometry, [math]\displaystyle{ \epsilon }[/math] can be calculated using correlations in terms of the "heat capacity ratio," or [math]\displaystyle{ C_r }[/math] and NTU:

[math]\displaystyle{ C_r \ = \frac{C_\mathrm{min}}{C_\mathrm{max}} }[/math]

[math]\displaystyle{ \ NTU }[/math] describes heat transfer across a surface

[math]\displaystyle{ NTU \ = \frac{U A}{C_\mathrm{min}} }[/math]

Here, [math]\displaystyle{ U }[/math] is the overall heat transfer coefficient, [math]\displaystyle{ A }[/math] is the total heat transfer area, and [math]\displaystyle{ C_{min} }[/math] 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:

[math]\displaystyle{ \dot Q= C_{min} (T_{o} -T_{i})_{min} = UA \Delta T_{LM} }[/math]

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 ([math]\displaystyle{ \Delta T_{LM} }[/math]). 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:

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

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

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

For a balanced counter-current flow heat exchanger (balanced meaning [math]\displaystyle{ C_r \ = 1 }[/math], which is a scenario desirable to reduce entropy):

[math]\displaystyle{ \epsilon\ = \frac{NTU}{1+NTU} }[/math]

A single-stream heat exchanger is a special case in which [math]\displaystyle{ C_r \ = 0 }[/math]. This occurs when [math]\displaystyle{ C_\mathrm{min}=0 }[/math] or [math]\displaystyle{ C_\mathrm{max}=\infty }[/math] 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]

[math]\displaystyle{ \epsilon \ = 1 - e^{-NTU} }[/math]

The effectiveness-NTU relationships for crossflow heat exchangers and various types of shell and tube heat exchangers can be derived only numerically by solving a set of partial differential equations. So, there is no analytical formula for their effectiveness, but just a table of numbers or a diagram. These relationships are differentiated from one another depending (in shell and tube exchangers) on the type of the overall flow scheme (counter-current, concurrent, or cross flow, and the number of passes) and (for the crossflow type) whether any or both flow streams are mixed or unmixed perpendicular to their flow directions.

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:

[math]\displaystyle{ c_p=\frac{dh}{dT} }[/math]

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:

[math]\displaystyle{ c_{p-x}= \frac{d \omega_x}{dP_x} }[/math]

Here, [math]\displaystyle{ \omega_x }[/math] represents the mass ratio of gas 'x' (meaning mass of gas 'x' relative to the mass of all other non-'x' gas mass) and [math]\displaystyle{ P_x }[/math] 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:

[math]\displaystyle{ c_{p-x}= \frac {M_x/M_{other}} {P_{other}} }[/math]

Here, [math]\displaystyle{ M_x }[/math] is the molecular weight of gas 'x' and [math]\displaystyle{ M _{other} }[/math] 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:

[math]\displaystyle{ NTU_x \ = \frac{U_mA_m}{\dot{m}c_{p-x}} }[/math]

Here, [math]\displaystyle{ U_m }[/math] is the overall mass transfer coefficient, which could be determined by empirical correlations, [math]\displaystyle{ A_m }[/math] is the surface area for mass transfer (particularly relevant in membrane-based separations), and [math]\displaystyle{ \dot{m} }[/math] 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 [math]\displaystyle{ C_r }[/math] becomes:

[math]\displaystyle{ C_r \ = \frac{(\dot{m}c_{p-x})_{min}}{(\dot{m}c_{p-x})_{max}} }[/math]

Effectiveness-NTU method for dehumidification applications

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

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

Here, [math]\displaystyle{ M_{wv} }[/math] is the molecular weight of water (vapor), [math]\displaystyle{ M_{air} }[/math] is the average molecular weight of air, [math]\displaystyle{ P_{air} }[/math] 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, [math]\displaystyle{ P_{wv,inlet} }[/math]. 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. 2.0 2.1 2.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. https://www.sciencedirect.com/science/article/pii/S1359431123015430. 
  3. J. H. Lienhard IV; J. H. Lienhard V (August 14, 2020). A Heat Transfer Textbook. Phlogiston Press. pp. 121–127. https://ahtt.mit.edu. 
  4. 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. https://www.sciencedirect.com/science/article/pii/S0196890422012699. 
  • 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.