Next Article in Journal
Multi-Electrode Resistivity Probe for Investigation of Local Temperature Inside Metal Shell Battery Cells via Resistivity: Experiments and Evaluation of Electrical Resistance Tomography
Previous Article in Journal
Energy, Exergy and Performance Analysis of Small-Scale Organic Rankine Cycle Systems for Electrical Power Generation Applicable in Rural Areas of Developing Countries
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Evaluation of Air-Based Heat Rejection Systems

1
Fraunhofer ISE, Fraunhofer Institute for Solar Energy Systems, Heidenhofstr. 2, 79110 Freiburg, Germany
2
ILK Dresden, Institute of Air Handling and Refrigeration, Bertolt-Brecht-Allee 20, 01309 Dresden, Germany
3
Centre of Applied Research Sustainable Energy Technology – zafh.net, Stuttgart University of Applied Sciences, Schellingstr. 24, 70174 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Energies 2015, 8(2), 714-741; https://doi.org/10.3390/en8020714
Submission received: 2 December 2014 / Accepted: 6 January 2015 / Published: 23 January 2015

Abstract

:
On the basis of the Number of Transfer Units (NTU) method a functional relation between electric power for fans/pumps and effectiveness in dry coolers and wet cooling towers is developed. Based on this relation, a graphical presentation method of monitoring and simulation data of heat rejection units is introduced. The functional relation allows evaluating the thermodynamic performance of differently sized heat rejection units and comparing performance among them. The method is used to evaluate monitoring data of dry coolers of different solar cooling field projects. The novelty of this approach is that performance rating is not limited by a design point or standardized operating conditions of the heat exchanger, but is realizable under flexible conditions.

Graphical Abstract

1. Introduction

Heat rejection is essential for a variety of processes in industrial, commercial and residential applications. In the latter two cases it is mostly used in air conditioning, cooling and refrigeration applications. The use of these technologies is ever increasing, e.g., according to [1] the market growth across the world only for air conditioning systems is expected to be above 5% for the next 3 years. One consequence of this enormous growth is an increase in electricity consumption for cooling processes. Considering the actual climate discussion and motivated by rising electricity prices around the globe, efforts are made by many stakeholders to increase the efficiency of the involved components (among which heat rejection plays a key role).
Waste heat from cold processes can be dissipated to different media such as water, soil and air. Therefore technologies used for heat rejection include night radiative cooling [2], sea-, lake- or groundwater cooling, including buried irrigation tanks [3] or swimming pools [4]. Yet the most commonly used heat sink is ambient air, as the above named options are rarely available or costly to exploit. As a result this study will concentrate on the air based technologies (e.g., dry coolers, wet cooling towers).
A first step towards the increase of efficiency is the determination of actual performance of heat rejection unit. Therefore standards for the determination of performance of heat rejection units exist. Some of the standards and performance evaluation methods will be presented below, to show the limitations of evaluation and the necessity of a more flexible method.
For dry coolers the European norm [5] defines one reference operating point at which the nominal capacity is determined (fixed temperatures for the air (25 °C) and water side (40 °C) and a fixed atmospheric pressure of 1013 hPa). Instructions on measurements and correction factors for differing conditions are given. The European Committee of Air Handling and Refrigeration Manufacturers (Eurovent) published a guideline for rating dry coolers [6]. Based on the standard defined in [5], the electric power for the fan/fans is measured. If the standard conditions are not fulfilled, instructions to correct the measurements are given. An energy ratio R , defined as nominal cooling capacity divided by the total certified power input of the fan motors at the standard rating conditions, is then calculated. Based on this value, a rating of dry coolers in an energy class (A++ for R     240   kW th / kW el to E with R   > 30   kW th / kW el ) can be performed. However the rating is limited to the standard operating point. Performance under varying conditions is left open. In [7] the German standard, according to which the energy consumption of buildings shall be determined, single average seasonal values for the electricity consumption of heat rejection equipment (including fans and pumps) are specified:
  • dry heat rejection systems: 0.045 kW el / kW th ;
  • open wet cooling towers with axial fan: 0.018 kW el / kW th ;
  • closed wet cooling towers with axial fan: 0.033 kW el / kW th .
The performance testing of series wet cooling towers must be carried out according to European Standard (EN) 13741 [8]. In this standard the procedure to verify the performance specified by the manufacturer is described. Yet the boundary conditions under which this is done are left up to the manufacturer. Therefore no comparability among products is facilitated, no labelling is possible. In [9] performance curves and tower characteristic maps for wet cooling towers are presented. The temperature range on the cold water side can be read off for different design flow rates, wet-bulb temperatures and cold water temperatures. Yet the method does not facilitate an easy comparison of different products nor does it allow the rating of different products. Also, a depiction of monitoring data in a similar way is precluded, if the available data is insufficient.
Apparently there is no description of performance or a definition of a satisfying coefficient of performance which allows the comparison of different Heat Rejection Units (HRUs). Especially as in practice HRUs rarely operate under the reference conditions defined in the standards, the question of their actual performance is open. This is even more so, as more and more frequently efforts are made to control both the air and the water volume flow according to the actual ambient conditions and load variations in order to minimize the auxiliary electricity consumption of such systems. One possibility to handle these flexible conditions is an annual and site specific simulation of the whole system including a suitable control strategy for the HRU. However, these simulations are time consuming and very specific and therefore not appropriate for evaluation of different HRUs. Consequently a method to evaluate HRUs not only in their design point or as one component in a whole system, but as an independent unit has to be elaborated.
In this paper we will develop a presentation method for the performance characterization of dry coolers by means of monitoring and simulation data at flexible operating conditions. We then extend this method to compare HRUs of different size and we illustrate the applicability of this method for monitoring and simulation data of dry coolers. The data utilized in this study originates from multiple solar cooling pilot projects which were intensively monitored. Finally we prove that an analog presentation method can be used to characterize the performance of wet cooling towers based on monitoring and simulation data at different operating conditions.

2. Comparing the Performance of Dry Coolers

As the application and the boundary conditions for heat rejection systems differ strongly from system to system it is essential to rate a system based on a parameter which does include the operating condition (OC, e.g., ambient temperature, water/cooling fluid supply and return temperatures, water/cooling fluid mass flow rate). The two fluids in a HRU are air and a cooling fluid. The cooling fluid (e.g., water, glycol-water mixture) will be abbreviated by “cf” throughout this paper.
The electric power consumption of the circulating pump and the fans ( P el in [Wel]) of a dry cooler is the basis of this paper for a rating of different design and materials used in heat exchangers, since it is the determinant for the operating costs (at least in case of dry coolers). The electric power for fans P el , Fan includes the power input of all fan motors used for the operation on the air side of the HRU. Likewise the electric power for pumps P el , Pump includes the power input of all pump motors on the cooling fluid side of the HRU. Later this has to be specified according to the system boundary of the HRU. The sum of P el , Fan and P el , Pump yields the electric power consumption of the HRU, neglecting the power input for monitoring devices, etc. In general P el can be determined by a not yet defined function, dependent on the investigated dry cooler configuration ( DC ) and the operating conditions:
P el =   P el ( DC , OC )
The dry cooler configuration describes the geometry and materials of the heat exchanger, is therefore fixed for a given heat exchanger and can be written as a subscript. The function P el is describing in this case a not explicitly known dependency of two parameters. As the operating condition is determined by the four parameters Q ˙ ,   T cf , in ,   T cf , out , T air , in [10] and [11] the electric power consumption of a given dry cooler can be expressed as:
P el , DC =   P el , DC ( Q ˙ , T cf , in ,   T cf , out , T air , in )
To depict this function explicitly, a detailed physical model or an extensive experimental study is required. However a reduction of parameters can be achieved without major losses in accuracy, as will be shown. The reduction is based on the effectiveness-NTU method, e.g., [12] and [13], which is based on the mass and energy balance of a heat exchanger, shown schematically in Figure 1. Therein the temperature effectiveness on cooling fluid side ( ϵ cf = Q ˙ / ( C ˙ cf  Δ T max )   ) is a function of the number of transfer units on air side ( NTU air = U   A / C ˙ air ), the heat capacity rate ratio on cooling fluid side ( C ˙ cf * = C ˙ cf / C ˙ air ) and the dry cooler configuration.
Figure 1. (a) Mass and energy balance of a counterflow dry cooler. (b) Scheme of temperatures over the length of the heat exchanger. The cooling effectiveness ϵ cf is defined by ϵ cf = Q ˙ / ( C ˙ cf  Δ T max ) and can be reformulated to due to ϵ cf = Δ T cf / Δ T max the energy balance Q ˙ = C ˙ cf  Δ T cf .
Figure 1. (a) Mass and energy balance of a counterflow dry cooler. (b) Scheme of temperatures over the length of the heat exchanger. The cooling effectiveness ϵ cf is defined by ϵ cf = Q ˙ / ( C ˙ cf  Δ T max ) and can be reformulated to due to ϵ cf = Δ T cf / Δ T max the energy balance Q ˙ = C ˙ cf  Δ T cf .
Energies 08 00714 g001
The heat capacity flow rates are defined as C ˙ cf = c ¯ p , cf   m ˙ cf and C ˙ air = c ¯ p , air   m ˙ air , wherein the specific heat capacities are mean values with a dependency on fluid temperature, which has to be analyzed. The parameter Δ T max is the inlet temperature difference ( Δ T max = T cf , in T air , in in [K]). The temperature effectiveness ϵ cf will from now on be called cooling effectiveness. The cooling effectiveness ranges from zero to one. In an ideal HRU the cooling fluid would be cooled to the air inlet temperature, resulting in a value of ϵ cf equal to one with electric power as low as possible. The geometrical configuration is fixed for a given dry cooler, so ϵ cf can be written as a function of three parameters:
ϵ cf = ϵ cf ( NTU air , C ˙ cf * ,   DC )
Assuming that the overall heat transfer coefficient U is a function of C ˙ air and C ˙ cf and the medium used as cooling fluid, NTU air and C ˙ cf * are given by:
NTU air =   NTU air ( C ˙ air , C ˙ cf , DC , cf )
C ˙ cf * =   C ˙ cf * ( C ˙ air , C ˙ cf )
This assumption is reasonable as long as the operating temperatures are within a limited range. In Table 1 the relevant properties of water, ethylene glycol based water solutions (30 Vol%) and air (relative humidity of 0%) for a temperature range of 20 to 45 °C (respectively 10 to 35 °C) are depicted.
Table 1. Thermodynamic properties of ethylene glycol solution, water and dry air at different temperatures (at 1 atm) and non-dimensional ratio between properties at the minimum temperature and the maximum temperature. Data calculated with CoolProp ([14] and [15]).
Table 1. Thermodynamic properties of ethylene glycol solution, water and dry air at different temperatures (at 1 atm) and non-dimensional ratio between properties at the minimum temperature and the maximum temperature. Data calculated with CoolProp ([14] and [15]).
MediumRange and RatioThermal Conductivity (W/m K)Dynamic Viscosity (10−5 kg/m s)Density (kg/m³)Specific Heat Capacity (103 J/kg K)
ethylene glycol solutionmin (20 °C)0.465216.6451038.0463.718
max (45 °C)0.487115.2521026.1533.789
ratio1.0480.5320.9891.019
watermin (20 °C)0.598100.160998.2074.184
max (45 °C)0.63559.577990.2134.18
ratio1.0610.5950.9920.999
airmin (10 °C)0.0251.7721.2471.006
max (35 °C)0.0271.8931.1461.007
ratio1.0741.0680.9191.001
In the air-conditioning and refrigeration industry, wet cooling towers are generally used in the range of 32 to 46 °C hot water temperature. A typical standard design condition for such cooling towers is 35 °C hot water to 29.4 °C cold water and 25.6 °C wet-bulb temperature [9].
The density difference in air due to different relative humidity is small in comparison to differences due to temperature changes (cf. [16]), therefore it will not be presented here. Differences in atmospheric pressure have a high impact on air density. Other properties in Table 1 are weakly affected when the pressure changes.
The Reynolds number (Re) and Prandtl number (Pr) of air as functions of the thermodynamic properties in Table 1 are slightly dependent on the temperature for constant capacity flow rate C ˙ air and constant DC (differences of maximum to minimum Reynolds number are 6.8% with respect to maximum Reynolds number, due to temperature dependent dynamic viscosity differences). So the assumption that the heat transfer coefficient is only dependent on the capacity flow rate C ˙ air and the dry cooler configuration DC is reasonable (especially as thermal conductivity is slightly increasing with temperature as well).
For the heat transfer coefficient of the cooling fluid side the dependency on temperature is much more pronounced, due to the strong effect of the temperature on viscosity. The differences in dynamic viscosity for a low mean temperature of 20 °C and a high mean temperature of 45 °C can be depicted by the relation of the viscosity, which is 0.53 for the ethylene glycol solution and 0.60 for water (cf. Table 1). To show the differences of heat transfer coefficients at low ( α cf , low ) and high temperatures ( α cf , high ) the relation is estimated by using the Dittus–Boelter correlation for turbulent flow in a smooth duct [13]. Therein the heat transfer for cooling is described by the Nusselt number:
Nu = 0.026  Re 0.8   Pr 0.3
The relation α cf , low / α cf , high is calculated for the same geometry and equal capacity flow rates, resulting in a relation of 0.71 for the ethylene glycol solution and 0.64 for water. Although the difference seems not to be negligible, the overall heat transfer coefficient U in a gas-to-liquid heat exchanger is typically determined by the gas side due to the much lower heat transfer coefficient. Assuming a dependency of cooling fluid side of 20% on U   ( U A / ( α cf A cf   ) = 0.2 ) , a relative error of 29% (corresponding to   α cf , low / α cf , high = 0.71 ) for α cf yields an acceptable relative error in U of 5.8%. Therefore Equation (4) is only weakly dependent on the temperature of the cooling fluid and thus this dependency will be ignored. Especially as in the estimation above two extreme values are compared to each other. Comparing one extreme value with a medium value (at water temperature of 32.5 °C) the relative error in U will further decrease. The relevance of the type of cooling fluid will be preserved, as the thermodynamic properties vary strongly from one cooling fluid to another. So on the one hand ϵ cf is, by definition:
ϵ cf = Q ˙ C ˙ cf  Δ T max
On the other hand due to Equations (3)–(5) ϵ cf is given by a function   ζ :
ϵ cf = ζ (   C ˙ cf , C ˙ air , DC , cf )
Taking the difference of Equations (7) and (8) yields a function τ :
τ ( C ˙ cf , C ˙ air , DC , cf , ϵ cf ) = ζ (   C ˙ cf , C ˙ air , DC , cf ) ϵ cf
For measurement data τ is equal to zero, as the five parameters C ˙ cf ,   C ˙ air ,   DC , cf and   ϵ cf are not independent of each other. In the mathematical formulation of Equation (9) the parameters are independent of each other and τ can therefore be different from zero. Then the function τ is only equal to zero if the choice of C ˙ cf ,   C ˙ air ,   DC and cf is in such a way that is corresponds to the choice of the fifth parameter   ϵ cf in the function τ . Therefore C ˙ air might be partly given by an implicit function ϕ :
C ˙ air = ϕ (   ϵ cf ,   C ˙ cf ,   DC , cf )  
The function ϕ will not be defined for every possible choice of ( ϵ cf ,   C ˙ cf ,   DC , cf ) , but it can be shown with the implicit function theorem [17], that if there exists one value for C ˙ air and fixed values for   ( ϵ cf ,   C ˙ cf ,   DC , cf ) , such that τ ( C ˙ cf , C ˙ air , DC , cf , ϵ cf ) = 0 , the function ϕ is defined in a neighbourhood around   ( ϵ cf ,   C ˙ cf ,   DC , cf ) . A necessary condition for validity of the implicit function theorem is that τ C ˙ air | ϵ cf ,   C ˙ cf ,   DC , cf 0 at these fixed values. The proof is based on the idea, that the function τ in Equation (9) is a monotonously increasing function of   C ˙ air and only in ζ dependent on   C ˙ air .
Therefore τ C ˙ air | ϵ cf ,   C ˙ cf ,   DC , cf = ζ C ˙ air | ϵ cf ,   C ˙ cf ,   DC , cf =   ζ C ˙ air | C ˙ cf ,   DC , cf > 0 , as it is obvious, that for an increasing capacity flow rate on air side (respectively increasing air mass flow rate) and constant capacity flow rate on water side (respectively constant water mass flow rate) the function ζ for cooling effectiveness will increase in a given dry cooler. Therefore τ C ˙ a i r | ϵ cf , C ˙ cf ,   DC , cf 0 holds (cf. Equation (9)), a necessary and sufficient requirement for the implicit function theorem. For the studied monitoring data Equation (10) is always applicable, however an extrapolation of Equation (10) for other parameters might not exist.
The relationship of Equation (10) for a simulated exemplary dry cooler (with a similar configuration as the dry cooler DC1 which will be presented in the subsequent Section 3) is shown in Figure 2. The figure shows that the four parameters ϵ cf ,   C ˙ cf ,   DC and cf have strong influence on the air heat capacity flow rate C ˙ air , whereas the influence of temperature on C ˙ air is minor. The air heat capacity flow rate is an increasing function of cooling effectiveness, with a vertical asymptote less than or equal to one.
In Figure 2a two different cooling fluids are used for simulation. The cooling fluid heat capacity flow rate is constant and equal for the cooling fluids, as well as the dry cooler flow configuration and the inlet temperatures on the cooling fluid and air side. For a fixed cooling effectiveness the difference in air heat capacity flow rate is up to 10% and therefore not negligible. Figure 2b,c show a similar result, with a strong influence of the changing parameters C ˙ cf and flow configuration on the curve. The flow configuration in Figure 2c is altered just be changing the number of passes on the water side. For different inlet temperatures (Figure 2d,e) and for a fixed cooling effectiveness the difference in air heat capacity flow rate is less 4%. This approves the proceeding not to consider the dependency on temperature in the implicit function ϕ in Equation (10).
Figure 2. Dependency of cooling fluid capacity flow rate versus cooling effectiveness on different parameters. Other parameters than the changed parameters are constant. Further explanation is given in the text. (a) Changing the cooling fluid; (b) Changing the cooling fluid heat capacity flow rate; (c) Changing the dry cooler flow configuration; (d) Changing cooling fluid inlet temperature; (e) Changing air inlet temperature.
Figure 2. Dependency of cooling fluid capacity flow rate versus cooling effectiveness on different parameters. Other parameters than the changed parameters are constant. Further explanation is given in the text. (a) Changing the cooling fluid; (b) Changing the cooling fluid heat capacity flow rate; (c) Changing the dry cooler flow configuration; (d) Changing cooling fluid inlet temperature; (e) Changing air inlet temperature.
Energies 08 00714 g002

2.1. Electric Power for Fans

Assuming that the electrical consumption for the fans P el , Fan is dependent on the mass flow rate of air only (ignoring the impact of temperature dependent density and viscosity), P el , Fan is dependant on C ˙ air only, as the heat capacity variation of air in the given temperature range (cf. Table 1) is negligible. The fan is fixed as one part of the dry cooler. Due to Equation (10) it is:
P el , Fan =   P el , Fan ( ϵ cf ,   C ˙ cf ,   DC , cf )
This simplification is necessary if later electric powers have to be compared without knowledge of temperatures at operation. In general the electric power is calculated by   P el , Fan = ( Δ p air m ˙ air ) / ( ρ air η Fan ) , with fan efficiency η Fan . Therefore a more precise assumption is to include a temperature dependent density in the calculation. For two different temperatures the ratio between the electric power for fans necessary to run an equal air capacity flow rate (and assuming c p , air ,  low = c p , air ,  high ) is given by:
P el , Fan , low P el , Fan , high = f low f high ( ρ air , high ρ air , low ) 2 η Fan , high η Fan , low
where the friction factor f low (respectively   f high ) is a function of Reynolds number and therefore dependent on the dynamic viscosity of air at a given temperature. As the relationship between friction factor and Reynolds number is different for all types of heat exchangers, it is not possible to develop a generally valid expression of f low / f high in terms of Reynolds numbers. The correlations given in literature are mostly based on the formula f = c 1   Re c 2 , with c 1 and   c 2 being constants dependent on geometry. For a flow through a staggered tube bundle c 2 is in the range of 1 to 0.25 [18]. A good overview for correlations of complex geometries is given in [13] (e.g., Offset Strip Fins, Louver Fins, Individually Finned Tubes, Plain Flat Fins and Corrugated Flat Fins on a Tube Array, Cross rod geometries). For an equal air capacity flow rate at different temperatures ( C ˙ air , low = C ˙ air ,   high ) and a constant heat exchanger geometry the ratio of Reynolds numbers is equal to the inverse ratio of dynamic viscosities η , as heat capacity of air is constant for different temperatures (cf. Table 1). Using the basic formula for friction factor to Reynolds number correlations yields:
f low f high = ( Re low Re high ) c 2   = ( η high η low ) c 2
This ratio is equal to 0.95 (respectively 0.99) for c 2 = 1 (respectively   c 2 = 0.25 ). We will therefore in the following analysis neglect the ratio of friction factors in Equation (12), accepting an error of 5% in the equation. On the differences of fan efficiency η Fan at different temperatures, but equal air capacity flow rates, such a general statement is not possible. Nevertheless as ambient temperatures during operation will be in a restricted range, we will ignore fan efficiency differences due to temperature differences.
Considering a mean temperature of inlet air of 25 °C (cf. standard in [6]) Equation (12) can then be reformulated to:
P el , Fan , 25 ° C = ( ρ air ρ air , 25 ° C ) 2 P el , Fan
The measured data P el , Fan (and   ρ air ) at an inlet air temperature different to 25 °C is thus transformed to a standardized electric power   P el , Fan , 25 ° C . The negligence of density differences of air due to temperature differences in Equation (11) is therefore not necessary anymore, if the electric power of fans is calculated as follows:
P el , Fan , 25 ° C ( ϵ cf ,   C ˙ cf ,   DC , cf ) =   ( ρ air ρ air , 25 ° C ) 2 P el , Fan ( ϵ cf ,   C ˙ cf ,   DC , cf )
The density ρ air , 25   ° C is at standard ambient pressure of 1013.25 hPa for dry air. Differences in atmospheric pressure and relative humidity on density of air are included in Equation (15).

2.2. Electric Power for Pump

The electrical consumption for the pump P el , Pump can be approximated by the capacity flow rate of cooling fluid; however the temperature of cooling fluid cannot be ignored, due to strong dependency of pressure drop on viscosity. Using the Blasius equation [13] for turbulent flow in a smooth duct the dissipated power ( P Diss = m ˙ cf Δ p cf / ρ cf ) on the cooling fluid side can be calculated exemplarily. The relation of dissipated energy at low temperature P Diss , cf , low   to dissipated energy at high temperature P Diss , cf , high   of the cooling fluid (cf. Table 1) is then given by:
P Diss , cf , low   P Diss , cf , high   = ( C ˙ cf ,   low C ˙ cf ,   high ) 0.75 ( c p , cf ,   low c p , cf ,   high ) 0.75 ( η cf ,  low η cf ,  high ) 0.25 ,
as the geometric parameters are constant.
Considering equal capacity flow rates yields P Diss , cf , low   / P Diss , cf , high   = 0.84 for ethylene glycol solution and P Diss , cf , low   / P Diss , cf , high   = 0.88 for water. Following the same idea as for the electric power for fans, a mean cooling fluid temperature of 40 °C is defined (cf. standard in [6]). Neglecting the ratio of heat capacities in Equation (16) and assuming that pump efficiency due to temperature differences of the cooling fluid can be neglected, the electric power for the pump at 40 °C is:
P el , Pump , 40 ° C ( C ˙ cf , DC , cf ) =   ( η cf ,   40 ° C η cf ) 0.25 P el , Pump ( C ˙ cf , DC , cf )
The measured data P el , Pump at an inlet air temperature different to 40 °C is thus transformed to a standardized electric power P el , Pump , 40 ° C , using a correction factor based on the differences in dynamic viscosity due to temperature differences. A total standardized electric power P el , 25 ° C , 40   ° C is then calculated to be:
P el , 25 ° C , 40 ° C ( ϵ cf ,   C ˙ cf ,   DC , cf ) = P el , Fan , 25 ° C ( ϵ cf ,   C ˙ cf ,   DC , cf ) +   P el , Pump , 40 ° C ( C ˙ cf , DC , cf )
Monitored and simulated data of electric power at temperatures within a specified range (cf. Table 1) can thus be transformed with Equation (18) to a comparable value for electric power.
Summing up, a functional relationship between operating conditions (expressed by ϵ cf and   C ˙ cf ) as well as geometrical and system configuration ( DC and cf ) and the electrical power demand to operate the heat rejection unit has been developed. This method is based on two validated models describing thermodynamic relationships, first the NTU method and second the determination of electrical power demand based on fan/pump efficiency. The method above merges the two models to one method, describing the cost of heat rejection units in terms of electrical power and the benefit in terms of cooling effectiveness.

3. Monitoring Data and Simulation of Dry Coolers

For the comparison of various types of dry coolers and the evaluation of measurement data this approach enables us to describe the electric power needed to run pumps and fans for a given dry cooler and a cooling fluid in use by solely knowing the value of C ˙ cf and ϵ cf as described in Equation (18). Further knowledge on temperatures to determine air density and cooling fluid viscosity will increase the accuracy. As the systems for solar cooling differ, so does the method to determine monitoring data. This results in a non-uniform definition of electric power for the HRU. In Figure 3 different system boundaries for electric powers are defined.
The first boundary (B1) includes only the electric power P el , Fan necessary to run the fan. If simulation data is used, the electric power for the fan is calculated by the quotient of dissipated energy ( P Diss , air = V ˙ air Δ p air ) and an assumed average fan efficiency of η Fan = 0.5 (this value is equal to 150% of the minimum fan efficiency in the energy efficiency requirements for axial fans with electrical power input of 0.8 kW (cf. [19]). The second boundary (B2) includes the electric powers from B1 plus the fraction of electric power for the cooling fluid pump needed to overcome the pressure drop on the fluid side of the dry cooler (not the tubes connecting chiller and dry cooler). For simulation data this fraction is given by   V ˙ cf Δ p cf / η Pump , with   η Pump = 0.3 . In B3 the whole electric power for the cooling fluid pump is taken into account, additional to the electric power for the fan (B1). Some dry coolers operate with a second heat exchanger, separating the cooling fluid (e.g., ethylene glycol solution) from the water in the chiller.
Figure 3. System boundaries for electric power in monitoring projects. Boundary 1 includes the electric power input of the fans on the air side only. All other boundaries include in addition a fraction of or the total electric power input of the pumps on the cooling fluid and water side.
Figure 3. System boundaries for electric power in monitoring projects. Boundary 1 includes the electric power input of the fans on the air side only. All other boundaries include in addition a fraction of or the total electric power input of the pumps on the cooling fluid and water side.
Energies 08 00714 g003
The electric power needed to run this second pump is added in B4 to the power needed in B3. An increased cooling water temperature after the heat exchanger is measured in this case. However determination of cooling effectiveness in this paper will always be based on temperatures in the cooling fluid circle, not in the water circle to the chiller. If this data is not available an estimation of temperature reduction in the second heat exchanger takes place. The monitoring data which will be analyzed does not contain information on all of these boundaries. A comparison of different dry coolers should take place only within the same boundary. In the following, data from seven sources were used to depict the performance graphically, according to Equations (2), (15) and (18):
  • DC1: Data calculated with the software CoilDesigner [20] for a dry cooler ( A HTS = 152   m ² ) ; air and cooling fluid (water) mass flow rate as well as both inlet temperatures are varied; data for B1 and B2 is given;
  • DC2: set of monitoring data from the German project SolCoolSys [21] ( A HTS unknown); air mass flow rate and both inlet temperatures are varied; The dry cooler is operated with an ethylene glycol solution (approximately 30 Vol%), the heat capacity is measured in the water circuit after a separating heat exchanger; temperatures for the cooling fluid side are calculated by assuming an effective temperature difference of 1.5 K in the separating heat exchanger; data for B1, B3 (temperatures calculated) and B4 is given;
  • DC3: set of monitoring data from the German project SolCoolSys [21] with a more recent dry cooler model than DC2 ( A HTS = 221   m ² ) ; same system as DC2; data for B1, B3 (temperatures calculated) and B4 is given;
  • DC4: Lab data measured by German manufacturer Thermofin [22] for a dry cooler ( A HTS = 46   m ² ) ; air and cooling fluid (water) mass flow rate as well as both inlet temperatures are varied; data for B1 and B2 is given;
  • DC5: set of monitoring data from the German project SolCoolSys [21] ( A HTS = 221   m ² ) ; air mass flow rate and both inlet temperatures are varied; the dry cooler is operated with an ethylene glycol solution (30 Vol%) as cooling fluid; temperatures are measured in the cooling fluid circuit and in the water circuit after a separating heat exchanger; data for B1, B3 and B4 is given;
  • DC6: Data calculated with the software Güntner Product Calculator [23] for a dry cooler with axial fans and multiple-row line-up (type GFW, A HTS = 1072   m ² ) ; air and cooling fluid (water) mass flow rate as well as both inlet temperatures are varied; data for B1 is given;
  • DC7: set of monitoring data by the Bavarian Center for Applied Energy Research (ZAE Bayern) [24]; The dry cooler with axial fans is assembled horizontally and produced by Güntner (type GFH, A HTS = 197   m ² ) ; air and cooling fluid (water) mass flow rate as well as both inlet temperatures are varied; data for B1 is given.
The simulated and monitored HRU were largely operated at a constant cooling fluid mass flow rate and therefore at a constant cooling fluid capacity flow rate. The ratio between C ˙ cf and the heat transfer surface (HTS) of the heat exchangers describes the type of operation. A low value of cooling fluid capacity flow rate at a high value of HTS may be an indication for high investment cost, compared to operation cost, as a larger heat exchanger has to be installed if a certain heat flow rate shall be achieved. However the cost for electric power for fans is reduced, as the HTS is large. In Figure 4 the cooling fluid capacity flow rate and the heat transfer surface of DC1 to DC7 are shown.
Figure 4. Differences of size and operating conditions of simulated and monitored dry coolers, expressed by heat transfer surface (HTS) and cooling fluid capacity flow rate ( C ˙ cf ) . Areas of same grey coloration correspond to same ratio of C ˙ cf and HTS.
Figure 4. Differences of size and operating conditions of simulated and monitored dry coolers, expressed by heat transfer surface (HTS) and cooling fluid capacity flow rate ( C ˙ cf ) . Areas of same grey coloration correspond to same ratio of C ˙ cf and HTS.
Energies 08 00714 g004
It is to be expected that DC4 has a high electricity demand, as the capacity flow rate is high in relation to the available HTS. Further DC5 und DC7 should show a better thermodynamic performance. In Figure 5 a general depiction of performance is shown. The electric power for fans is plotted versus heat flow rate. According to Equation (2) the electric power is further dependent on temperatures. This dependency is not included in this plot. The good performance of DC1 and DC6 (high Q ˙ at low P el , Fan ) is at least partly due to a larger driving temperature difference between the two fluids. An evaluation of performance based on this plot is not possible.
In Figure 6 the electric power for fans at 25 °C is normalized with the cooling fluid capacity flow rate. This specific electric power demand is plotted versus cooling effectiveness, according to Equation (15). The reason for normalizing the electric power demand is to compare different sized HRU with each other. A performance evaluation based on Figure 6 yields for every cooling ratio a best HRU. For example, at ϵ cf = 0.3 DC7 has the lowest specific electric power demand, DC4 the highest. At ϵ cf = 0.2 DC1 is better than DC7. Ignoring the scattering of monitoring data, each dry cooler can be represented by a curve. The courses of the curves are similar, ranging from a straight line in the log-plot to a slightly bended curve, as an increase in cooling effectiveness can only be achieved when increasing the air volume flow and therefore the electric power for fans.
Figure 5. Electric power for fans ( P el , Fan ) versus rejected heat ( Q ˙ ) for simulated and monitored dry coolers (in boundary B1). This depiction is not suitable to compare performance of HRU.
Figure 5. Electric power for fans ( P el , Fan ) versus rejected heat ( Q ˙ ) for simulated and monitored dry coolers (in boundary B1). This depiction is not suitable to compare performance of HRU.
Energies 08 00714 g005
Figure 6. Electric power for fans at 25 °C ( P el , Fan , 25 ° C ) normalized with cooling fluid capacity flow rate versus cooling effectiveness ( ϵ cf ) for simulated and monitored dry coolers (in boundary B1).
Figure 6. Electric power for fans at 25 °C ( P el , Fan , 25 ° C ) normalized with cooling fluid capacity flow rate versus cooling effectiveness ( ϵ cf ) for simulated and monitored dry coolers (in boundary B1).
Energies 08 00714 g006
The information given in the plot is diverse; however two findings shall be touched upon. First DC2 has been replaced by DC3 during monitoring. The improvement is evident: a three-fold decrease of the specific electric power for fans at a cooling effectiveness of ϵ cf = 0.4 (from 0.1 W el / ( W th / K ) for DC2 to 0.03 W el / ( W th / K ) for DC3) and further an increase in cooling effectiveness during operation from a range of ϵ cf = (0.1,0.5) to ϵ cf = (0.4,0.7). Second the dry coolers DC3 and DC5 are the same type of cooler, but operated at different locations and different cooling fluid mass flow rates. DC5 shows a better performance within the total range of cooling effectiveness, e.g., the specific electric power for fans at ϵ cf = 0.4 is only 66% of DC3.
To determine the electric power for fans for a specific rate of heat flow Q ˙ and a specific temperature driving force Δ T max , the cooling effectiveness has to be calculated first with Equation (7), based on the cooling fluid capacity flow rate used in Figure 6. Then the electric power demand can be determined from the curves in the plot.
The following sample calculation is intended to clarify the application of Figure 6. Assuming   Q ˙ = 20  kW , ΔT max = 10   K and a cooling fluid capacity flow rate of   C ˙ cf = 4000   W / K Equation (7) yields a cooling effectiveness of   ϵ cf = 0.5 . The required cooling fluid capacity flow rate can, e.g., be realized with DC2 or with two devices of DC7 (based on the data). The specific electric power demand for DC2 is approximately 0.09   W el / ( W th / K ) , for DC7 it is 0.03   W el / ( W th / K ) . As C ˙ cf is equal for both applications the electric power for fans is three times higher for DC2 (360 W) than for two devices of DC7 (120W). However the heat transfer surface is doubled for the two devices of DC7, resulting probably in higher investment cost for application of DC7 at this operating point.
The difference between plotting electric power for fans at 25 °C (Equation (15)) and at the monitored ambient temperatures (Equation (11)) is marginal with ± 2.6 % . The scattering of monitoring data is due to a transient heat transfer process and a non-constant value for cooling fluid capacity flow rate (limited to a certain range) to increase the usable data points.
In Figure 7 the specific electric power for the pump versus the cooling effectiveness (cf. Equation (17)) is depicted. As the cooling fluid capacity flow rate is constant the electric power is almost constant (cf. Equation (16)). The values for specific electric power for the pump range from 0.2   W el / ( kW th / K ) to   10   W el / ( kW th / K ) . The difference in demand of the dry coolers is due to the specific fraction of electric power that is used in this calculation. For DC1 and DC4 the fraction of B2 is taken (without B1). For DC2, DC3 and DC5 boundary B3 (without B1) is taken into account. No information on electric power for pumps was available for DC6 and DC7. Further the differences in cooling fluid capacity flow rate yield a higher value of specific electric power for the pump for DC1 (compared to DC4).
Adding the values in Figure 6 and Figure 7 yields a specific electric power for the HTU (cf. Equation (18)). This is depicted in Figure 8. The fraction of electric power for the pump of the entire electrical demand is between 50% and 80% for DC2, DC3 and DC5 (dependent on cooling effectiveness). Therefore at lower cooling effectiveness the differences between the three DCs is less dominant than in Figure 6. A Comparison to DC1 and DC4 is possible only to a limited extent as the electric power for pumps are within different system boundaries for electric power.
The validity of the method is shown in particular by the unique functional relationship for the simulated data in Figure 6 and Figure 8, which reveals no scattering at all.
Figure 7. Electric power for pump at 40 °C ( P el , Pump , 40 ° C ) normalized with cooling fluid capacity flow rate versus cooling effectiveness ( ϵ cf ) for simulated and monitored dry coolers (DC1 and DC4 in boundary B2 and DC2, DC3 and DC5 in B3; Excluding B1).
Figure 7. Electric power for pump at 40 °C ( P el , Pump , 40 ° C ) normalized with cooling fluid capacity flow rate versus cooling effectiveness ( ϵ cf ) for simulated and monitored dry coolers (DC1 and DC4 in boundary B2 and DC2, DC3 and DC5 in B3; Excluding B1).
Energies 08 00714 g007
Figure 8. Electric power for fans and pump at 25 °C and 40 °C ( P el , 25 ° C , 40 ° C ) normalized with cooling fluid capacity flow rate versus cooling effectiveness ( ϵ cf ) for simulated and monitored dry coolers (DC1 and DC4 in boundary B2 and DC2, DC3 and DC5 in B3).
Figure 8. Electric power for fans and pump at 25 °C and 40 °C ( P el , 25 ° C , 40 ° C ) normalized with cooling fluid capacity flow rate versus cooling effectiveness ( ϵ cf ) for simulated and monitored dry coolers (DC1 and DC4 in boundary B2 and DC2, DC3 and DC5 in B3).
Energies 08 00714 g008
In this study we will not present a method to combine investment cost (related strongly to size of heat transfer surface, cf. Figure 4) and operation cost (related strongly to electricity demand) to one parameter characterizing total cost, as additional information on operating conditions over the entire operational lifetime have to be given. These conditions vary strongly from system to system and total cost has to be calculated in the individual case. The advantage of the depiction in Figure 6 to Figure 8 is that all necessary information for changing operating conditions are included and operation cost for pumps and fans can easily be calculated if the operating conditions over the entire operational lifetime are given and converted to cooling efficiency. Different designs and configuration of heat rejection units can only be judged on the level of combined investment and operating costs, otherwise if ignoring investment cost, a heat exchanger with very large heat transfer surface would be preferable. Nevertheless the depiction method described above shows the difference in operation cost due to electricity demand for different designs and offers a possibility to compare the investment cost due to size of heat transfer surface (cf. Figure 4).

4. Wet Cooling Towers

The effectiveness-NTU method is generally applied for heat exchanger based on sensible heat transfer. Yet in the following chapter it will be shown, that a similar method as the one used for dry coolers above can be derived for wet cooling towers. Due to the lack of appropriate monitoring data, the method is only developed theoretically and not yet applied to measured data.
To start with, the following assumptions are adopted from [25]:
  • Heat and mass transfer is in a direction normal to the flows only;
  • Negligible heat and mass transfer through the tower walls to the environment;
  • Negligible heat transfer from the tower fans to the air or water streams;
  • Constant water and dry air specific heats;
  • Constant heat and mass transfer coefficients throughout the tower;
  • Value of Lewis number equal to unity;
  • Water loss due to evaporation is negligible (constant m ˙ w );
  • Uniform temperature throughout the water stream at each cross section;
  • Uniform cross sectional area of the tower.
In Figure 9 a counterflow wet cooling tower is shown schematically, including important states and dimensions. The steady state energy and mass balances on an incremental volume ( d V ) of a wet cooling tower yield the equations (cf. Figure 9, [25,26]):
d q = m ˙ w   d h w = m ˙ air   d h air =   h D A V   ( h air h s , w ) d V
Therein the heat, transferred to air ( d q ) is given by the product of water mass flow rate ( m ˙ w ) and incremental change in enthalpy of water ( d h w ). Assumption 7 is used here to consider a constant water flow rate at each point of the tower, which is equal to the inlet flow. The air mass flow rate ( m ˙ air ) and the incremental change in specific enthalpy of moist air per mass of dry air ( d h air ) yield the same   d q . Assumption 6 of Le = 1 is very strong and often investigated (and criticized) in literature [27]. Due to this assumption the enthalpy difference h air h s , w between the moist air and the saturated air at temperature T w represents the driving force for the cooling process. The convective mass transfer coefficient ( h D ) and the surface area of water droplets ( A V ) define the factor in front of the driving force.
Equation (19) is used as a basis for the NTU method for cooling towers. To develop the method it is necessary to define effectiveness, NTU and capacity flow rates appropriately.
In [26] and [28] this is done, by defining the averaged saturation specific heat:
c p , sat = Δ h s , w Δ T w = h s , w , in  h s , w , out T w , in   T w , out
with units of specific heat. This parameter allows the estimate of the derivative d h w d h s , w   (as a function of temperature) for temperatures in the range of T w , in to T w , out by c p , w c p , sat as   c p , w c p , sat = Δ h w Δ h s , w d h w d h s , w , with Δ h w = h w , in h w , out .
Figure 9. Mass and energy balance of a wet counter flow cooling tower (cf. [25]). The term m ˙ air   d X has been neglected on the water side.
Figure 9. Mass and energy balance of a wet counter flow cooling tower (cf. [25]). The term m ˙ air   d X has been neglected on the water side.
Energies 08 00714 g009
The effectiveness for sensible heat exchangers is defined as the quotient of actual heat transfer to maximum possible heat transfer. For heat exchangers with evaporation it is reasonable to define the effectiveness based on the maximum possible air-side heat transfer [26]. This air-side heat transfer effectiveness ϵ air is defined as:
ϵ air = Q m ˙ air  Δ h max ˙
with Δ h max =   h s , w , in h air , in and the rate of heat flow in the cooling tower   Q ˙ .
Using Assumption 5 of constant mass transfer coefficients throughout the tower and Assumption 8 of a uniform temperature throughout the water stream at each uniform cross section, Equations (19) and (20) yield:
d q = m ˙ w c p , w c p , sat   d h s , w =   m ˙ air   d h air   =   h D A   ( h air h s , w ) d z / L
Following the analysis for differential equations in [13] one can obtain the relationship between ϵ air ,   m * = m ˙ air m ˙ w c p , w c p , sat   and NTU = h D A   m ˙ air as:
ϵ air = 1 exp [ NTU   ( 1 m * ) ] 1 m * exp [ NTU   ( 1 m * ) ]
for m * < 1 . If m * > 1 holds true, a similar analysis as presented can be performed. For flow geometries different to counter-flow an explicit function comparable to Equation (23) is not always possible to develop.
The following analysis uses the reduction of parameters in Equation (23) to find the necessary parameters to describe the electric power of fans and pumps of the cooling tower. Therefore an exact description of the dependencies among the parameters has to be presented. Due to Equation (22) the water mass flow rate is given after integrating by:
m ˙ w = Q ˙ Δ h s , w   c p , sat c p , w
The water-side heat transfer effectiveness for a cooling tower is defined as:
ϵ w , wet = Δ h s , w Δ h max
Assuming that the mass transfer coefficient h D is only dependent on m ˙ air ,   m ˙ w ,   c p , sat and the cooling tower flow arrangement CT (this includes the design, materials, ...) (cf. [26]) NTU is:
NTU = NTU ( m ˙ air , m ˙ w , c p , sat , CT )
Hence ϵ air is on the one hand due to Equations (23) and (26) and m * = m * ( m ˙ air ,   m ˙ w ,   c p , sat ) a not explicitly known function ζ :
ϵ air   = ζ ( m ˙ air , m ˙ w ,   c p , sat ,   CT )
On the other hand due to its definition in Equation (21) ϵ air is a known function φ :
ϵ air   = φ ( m ˙ air , Q ˙ Δ h max )
Using:
Q ˙ Δ h max = m ˙ w c p , w c p , sat ϵ w , w e t
and taking the difference of Equations (27) and (28) yields a function   τ , that solely equals zero if the parameters   m ˙ air , m ˙ w ,   c p , sat , ϵ w , wet ,   CT are choosen appropriate:
τ ( m ˙ air , m ˙ w ,   c p , sat , ϵ w , wet ,   CT )   = ζ ( m ˙ air , m ˙ w ,   c p , sat ,   CT ) φ ( m ˙ air , m ˙ w ,   c p , sat ,   ϵ w , wet )
Therefore m ˙ air is given implicitly as a function ϕ :
m ˙ air = ϕ ( m ˙ w ,   c p , sat ,   ϵ w , wet , CT )  
for some values of m ˙ w ,   c p , sat ,   ϵ w , wet , CT . If there exist one value for m ˙ air and fixed values for ( m ˙ w ,   c p , sat ,   ϵ w , w e t , CT ) fixed , such that τ ( m ˙ air , m ˙ w ,   c p , sat , ϵ w , w e t ,   CT ) = 0 is fulfilled, the function ϕ is defined in a neighborhood around   ( m ˙ w ,   c p , sat ,   ϵ w , wet , CT ) fixed . Mathematically this can be proven again with the implicit function theorem [17], showing that τ m ˙ air | m ˙ w ,   c p , sat , ϵ w , wet ,   CT 0 :
ζ m ˙ air | m ˙ w ,   c p , sat ,   CT = Q ˙ m ˙ air  Δ h max m ˙ air | m ˙ w ,   c p , sat ,   CT = Q ˙ Δ h max m ˙ air | m ˙ w ,   c p , sat ,   CT 1 m ˙ air Q ˙ Δ h max   m ˙ air 2
And by using Equation (29):
φ m ˙ air | m ˙ w ,   c p , sat ,   ϵ w , w e t = m ˙ w m ˙ air c p , w c p , sat ϵ w , wet m ˙ air | m ˙ w ,   c p , sat ,   ϵ w , wet = m ˙ w c p , w ϵ w , wet c p , sat   m ˙ air 2 =   Q ˙ Δ h max   m ˙ air 2
with Equation (24) it follows that:
ζ m ˙ air | m ˙ w ,   c p , sat ,   CT φ m ˙ air | m ˙ w ,   c p , sat ,   ϵ w , wet   =   Q ˙ Δ h max m ˙ air | m ˙ w ,   c p , sat ,   CT 1 m ˙ air = Δ h s Δ h max m ˙ air | m ˙ w , c p , sat , CT m ˙ w m ˙ air c p , w c p , sat > 0
as Δ h s Δ h max m ˙ air = ϵ w , wet m ˙ air > 0 for constant m ˙ w , c p , sat , CT ; which is a necessary and sufficient condition for the implicit function theorem. Therefore Equation (31) is valid.
Following the approach for dry coolers, we assume that the power for fans and pumps P el is only dependent on the mass flow rates m ˙ air , m ˙ w and on the density of air, viscosity of water and possibly on c p , sat , as well as on the cooling tower flow arrangement   CT .
If friction factor differences on the air side due to different inlet temperatures (same air mass flow rate) can be neglected, due to Equation (31) the dry cooler Equation (15) can be adapted:
P el , Fan , 25 ° C ( m ˙ w ,   ϵ w , wet ,   c p , sat ,   CT ) =   ( ρ air ρ air , 25 ° C ) 2 P el , Fan ( m ˙ w ,   ϵ w , wet ,   c p , sat ,   CT )
The density ρ air , 25 ° C is again calculated at standard ambient pressure of 1013.25 hPa for dry air. The density ρ air shall represent air density at the CT inlet. If the friction factor Reynoldsnumber correlation for the water side is still given by the Blasius Equation [13], Equation (17) can be adapted to:
P el , Pump , 40 ° C ( m ˙ w ,   CT ) =   ( η cf , 40 ° C η cf ) 0.25 P el , Pump ( m ˙ w , CT )
This is a strong simplification of pressure drop origination. Valves or other devices might increase the pressure drop significantly and the exponent of 0.25 in the Equation (36) has to be adapted.
Combining Equation (35) and a potentially adapted Equation (36) results in:
P el , 25 ° C , 40 ° C ( m ˙ w ,   ϵ w , wet ,   c p , sat ,   CT ) = P el , Fan , 25 ° C ( m ˙ w ,   ϵ w , wet ,   c p , sat ,   CT ) +   P el , Pump , 40 ° C ( m ˙ w , CT )
Comparison of Equations (18) and (37) shows that the only difference in the rating of dry coolers and wet cooling towers is that in the latter case one has to add the saturation specific heat c p , sat in the examination and the definition of effectiveness is based additionally to temperature on humidity. The specification for the cooling fluid in a wet cooling tower with water is not necessary.

5. Simulation of Wet Cooling Towers

In order to validate and show the benefit of the mathematical unification of NTU method and electrical power demand in Section 4, we examined the cooling tower performance of the Type 51 cooling tower model in the transient systems simulation program TRNSYS 17. This validated model is based on the work of Braun [26,29] and the Equipment Guide of the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) [28]. A detailed description of TRNSYS and the model can be found in [30,31]. Coefficients of mass transfer correlation and overall performance data have to be given as an input to the model. In a next step the operating conditions (temperatures and mass/volumetric flow rates) have to be specified. The objective of the following application of the model is to show the functional relationship of Equation (35) for a variety of data. The range of input parameters ( V ˙ air ,   T air , in , T air , in , wb , m ˙ w and T w , in ) is given in Table 2. The water outlet temperature T w , out and the electric power for the fan P el , Fan are outputs of the model. Further overall performance data is oriented on the Axima EWK 036/06 cooling tower with a rated cooling capacity of 45 W for the cooling down of water from 32 to 26 °C at a wetbulb temperature of T air , in , wb = 20 °C.
Table 2. Range of input parameters to Type 51 cooling tower model in TRNSYS 17.
Table 2. Range of input parameters to Type 51 cooling tower model in TRNSYS 17.
VariableDescriptionLower BoundaryUpper BoundaryUnits
V ˙ air air volumetric flow rate1,3004,500m³/hr
T air , in air dry bulb temperature2030°C
T air , in , wb air wet bulb temperature929°C
m ˙ w water mass flow rate3,0006,000kg/hr
T w , in water inlet temperature2238°C
The results of the simulation are shown in Figure 10 by plotting heat flow rate Q ˙ versus the electric power for fans P el , Fan . A functional relationship between these two quantities alone is not identifiable. With this depiction no information on performance can be worked out. However Equation (35) is suited for an evaluation of monitoring data for cooling towers. As the water mass flow rate for most cooling tower is constant a normalization of electric power for fans and pumps (share of cooling tower) with water mass flow rate is possible and necessary to compare different sized cooling towers analog to the normalization method for dry coolers. We will follow this procedure even though the size of the simulated cooling tower (CT) is constant and we will filter the data for constant water mass flow rates.
Figure 10. Electric power for fans ( P el , Fan ) versus rejected heat ( Q ˙ ) for simulated wet cooling tower Type 51 in TRNSYS 17. This depiction is not suitable to compare performance of heat rejection units.
Figure 10. Electric power for fans ( P el , Fan ) versus rejected heat ( Q ˙ ) for simulated wet cooling tower Type 51 in TRNSYS 17. This depiction is not suitable to compare performance of heat rejection units.
Energies 08 00714 g010
As the performance of CT additionally depends on c p , sat , which approximately represents the absolute temperature at which the process takes place and thus the capacity of the ambient air to take up humidity, data is filtered for c p , sat as well. For the unfiltered data c p , sat is in the range of 3250 J/(kg K) to 7000 J/(kg K) due to the defined input parameter to the Type 51 model. The applied filters are described in Table 3.
Table 3. Description of filtered data for simulated wet cooling tower Type 51 in TRNSYS 17.
Table 3. Description of filtered data for simulated wet cooling tower Type 51 in TRNSYS 17.
Filter m ˙ w in kg/s c p , sat in J/(kg K)
A1.4 ± 0.34,000 ± 100
B1.4 ± 0.34,900 ± 100
C0.9 ± 0.34,000 ± 100
D0.9 ± 0.34,900 ± 100
In Figure 11 the quotient of electric power for the fan and water capacity flow rate (specific electric power) is plotted versus cooling effectiveness for cooling tower ϵ w , wet . The unfiltered data (light blue) shows a large scattering of specific electric power with a factor of approximately 10 between the lowest and highest values. However a trend of increasing specific electric power with increasing cooling effectiveness is identifiable. Using the filter A to D from Table 3 yields a clear functional relationship (with only small scattering) which has been proven before in Equation (35). A further decrease of the filtered interval (e.g., from ±0.3 kg/s to ±0.1 kg/s for   m ˙ w ) would on the one hand reduce the small scattering of data, but on the other hand reduce the applicable data and is therefore not performed. The specific electric power of the cooling tower is lower for the low water mass flow rates of 0.9 kg/s (at equal cooling effectiveness), with the disadvantage of a lower heat flow rate. A higher c p , sat yields in the four studied cases a lower specific electric power. This is in general not true and can easily be checked by combining Equation (23), Equations (28) and (29) and solving for ϵ w , wet as a function of c p , sat . The course of the curves in Figure 11 is close to a straight line and due to the logarithmic scaling of the y-axis an exponential relationship between specific electric power and cooling effectiveness is proposed in this range of cooling effectiveness. This direct correlation shows the validity of Equation (31) with respect to Equation (35) exemplarily.
Figure 11. Electric power for fans at 25 °C ( P el , Fan , 25 ° C ) normalized with water capacity flow rate ( C ˙ w ) versus cooling effectiveness for cooling towers ( ϵ w , wet ) for simulated wet cooling tower Type 51 in TRNSYS 17. Unfiltered data (light blue) is based on input parameters in Table 2. Data A to D is filtered for water mass flow rate and averaged saturation specific heat. The filter is described in Table 3.
Figure 11. Electric power for fans at 25 °C ( P el , Fan , 25 ° C ) normalized with water capacity flow rate ( C ˙ w ) versus cooling effectiveness for cooling towers ( ϵ w , wet ) for simulated wet cooling tower Type 51 in TRNSYS 17. Unfiltered data (light blue) is based on input parameters in Table 2. Data A to D is filtered for water mass flow rate and averaged saturation specific heat. The filter is described in Table 3.
Energies 08 00714 g011
In the future, this approach shall be applied to monitoring data of wet cooling towers collected within the IEA SHC Task 48. The respective report will be published in [32].

6. Conclusions

We have developed a method to depict monitoring data of heat rejection units for dry cooler and wet cooling towers with varying operating conditions. The method allows an evaluation of heat rejection unit performance and a rating of units, with reduced amount of input parameters, representing the relevant operation data. For evaluation the normalized electric power is plotted versus the cooling effectiveness for fixed cooling fluid capacity flow rate. Independent of size the heat rejection unit performance can be depicted in one figure. Based on this plot heat rejection units and their operation can be directly compared, resulting in clear illustration of differences. In future the method could be used to compare devices in operation with design points and it could be extended to hybrid coolers and coolers with evaporative pre-cooling. A graphical comparison of dry and wet coolers is in progress.

Acknowledgments

The authors acknowledge the financial support from the German Federal Ministry for the Environment, Nature Conservation and Nuclear Safety (BMU) for the SolaRück Project [33].

Nomenclature

A V
surface area of water droplets per unit volume of tower (m2/m3)
A
Heat transfer surface for DC and surface area of water droplets for CT (m2)
c p , sat
Averaged saturation specific heat cp,sat=hs,w,in hs,w,outTw,in Tw,out (J/(kg K))
C ˙ *
Heat capacity rate ratio (-)
C ˙
Heat capacity flow rate C ˙ = m ˙   c p (W/K)
c p
Specific heat capacity (J/(kg K))
f
Friction factor (-)
h air
Enthalpy of moist air per mass of dry air ( J / kg dry air )
h w
Enthalpy of water evaluated at T w ( J / kg dry air )
h s , w
Enthalpy of saturated moist air evaluated at T w ( J / kg dry air )
Δ h
Enthalpy difference ( J / kg dry air )
h D
Convective mass transfer coefficient (kg/(m2 s))
L
Length of heat exchanger (m)
Le
Lewis number (-)
m ˙
Mass flow rate (kg/s)
Nu
Nusselt number (-)
NTU
Number of transfer units NTU = U A / C ˙ (-)
Δ p
Pressure difference (Pa)
Pr
Prandtl number (-)
P el
Electric power (W)
P el , 25 ° C , 40 ° C
Electric power for fans and pump at inlet temperatures of 25 °C for dry air and 40 °C for cooling fluid at 1 atm (W)
P el , Fan
Electric power input of all fan motors used for the operation on the air side of the heat rejection unit (W)
P el , Pump
Electric power input of all pump motors used for the operation on the cooling fluid side of the heat rejection unit within the system boundary (W)
P Diss
Dissipated power P Diss = m   ˙ Δ p / ρ (W)
Q ˙
Heat flow rate (W)
R
Energy ratio R = Q ˙ / P el , Fan (-)
Re
Reynolds number (-)
T
Temperature (K)
Δ T max
Temperature difference between both inlet fluids (K)
U
overall heat transfer coefficient (W/(m²K))
V ˙
Volume flow rate (m3/s)
X
Humidity ratio of moist air ( kg w / kg dry air )
α
heat transfer coefficient (W/(m2K))
ϵ cf
Cooling effectiveness for dry cooler ϵ cf = Q ˙ / ( C ˙ cf  Δ T max ) , see Figure 1 (-)
ϵ w , wet
Cooling effectiveness for cooling tower ϵ w , wet = Δ h s , w / Δ h max (-)
ϵ air
Air-side heat transfer effectiveness (-)
η
Dynamic viscosity (kg/(m s))
η Fan
Fan efficiency (-)
η Pump
Pump efficiency (-)
λ
Thermal conductivity (W/(m K))
ρ
Density (kg/m3)
τ ,  ζ ,   ϕ , φ
Not explicitly known functions

Subscripts and shortcuts

DC
Dry Cooler
Diss
Dissipation
HTS
Heat Ttransfer Surface
HRU
Heat Rejection Unit
OC
Operating Conditions
cf
Cooling fluid
air
Ambient air
high
High temperature (45 °C for cooling fluid, 35 °C for air)
in
Inlet
low
Low temperature (20 °C for cooling fluid, 10 °C for air)
out
Outlet
w
Water
wb
Wet bulb

Author Contributions

This study is part of the SolaRück Project [33], wherein all authors developed within several meetings the theoretical background of the method of presentation. The measurement data has been collected by Björn Nienborg, Hannes Fugmann and Gregor Trommler from different internal and external projects. The simulations have been conducted by Hannes Fugmann for dry coolers and by Björn Nienborg for wet cooling towers. These two authors are primarily responsible for data analysis and preparation of data for graphical representation, however all authors had input in this process. The mathematical proofs have been carried out by Hannes Fugmann, supported with physical background information on wet cooling towers by Antoine Dalibard. Hannes Fugmann accomplished the writing of the paper in close and steady cooperation with all authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alfonso, O. BSRIA in the US and AC Trends. Presented at AHR Expo, New York, NY; 2014. Available online: http://de.slideshare.net/BSRIA/air-conditioning-usa-2014-ahr-expo-edited-market-research (accessed on 13 January 2015).
  2. Eicker, U.; Dalibard, A. Photovoltaic–thermal collectors for night radiative cooling of buildings. Sol. Energy 2011, 85, 1322–1335. [Google Scholar] [CrossRef]
  3. Palacín, F.; Monné, C.; Alonso, S. Improvement of an existing solar powered absorption cooling system by means of dynamic simulation and experimental diagnosis. Energy 2011, 36, 4109–4118. [Google Scholar] [CrossRef]
  4. Eicker, U.; Pietruschka, D.; Pesch, R. Heat rejection and primary energy efficiency of solar driven absorption cooling systems. Int. J. Refrig. 2012, 35, 729–738. [Google Scholar] [CrossRef]
  5. Heat Exchangers—Air Cooled Liquid Coolers (“dry coolers”)—Test Procedures for Establishing the Performance; German version EN 1048:2014. Available online: http://www.beuth.de/de/norm/din-en-1048/198087229?SearchID=833923486 (accessed on 13 January 2015).
  6. Eurovent Rating Standard for DX Air Coolers, Air Cooled Condensers, Dry Coolers; Eurovent Certification Company: Paris, France, 2013.
  7. Energy Efficiency of Buildings—Calculation of The Net, Final and Primary Energy Demand for Heating, Cooling, Ventilation, Domestic Hot Water and Lighting—Part 1: General Balancing Procedures, Terms and Definitions, Zoning and Evaluation of Energy Sources. Available online: http://www.beuth.de/de/vornorm/din-v-18599-1/142650957?SearchID=833917916 (accessed on 13 January 2015).
  8. Thermal Performance Acceptance Testing of Mechanical Draught Series Wet Cooling Towers; German Version EN 13741:2003. Available online: http://www.beuth.de/de/norm/din-en-13741/64656756?SearchID=833921369 (accessed on 13 January 2015).
  9. American Society of Heating, Refrigerating and Air-Conditioning Engineers ASHRAE Handbook 2008: Heating, Ventilating, and Air-Conditioning Systems and Equipment; ASHRAE: Atlanta, GA, USA, 2008. Available online: https://www.ashrae.org/resources--publications/handbook (accessed on 13 January 2015).
  10. Webb, R.L.; Kim, N.-H. Principl of Enhanced Heat Transfer; Taylor Francis: New York, NY, USA, 1994. [Google Scholar]
  11. VDI-Guideline 2049: Thermal Acceptance and Performance Tests on Dry Cooling Towers. VDI-Verlag GmbH: Düsseldorf, Germany, 1981. Available online: https://www.vdi.de/uploads/tx_vdirili/pdf/1361735.pdf (accessed on 13 January 2015).
  12. Kröger, D.G. Air-cooled Heat Exchangers and Cooling Towers; PennWell Books: Tulsa, OK, USA, 2004. [Google Scholar]
  13. Shah, R.K.; Sekulić, D.P. Fundamentals of Heat Exchanger Design; Wiley: Hoboken, NJ, USA, 2003. [Google Scholar]
  14. Bell, I.H.; Wronski, J.; Quoilin, S.; Lemort, V. Pure and pseudo-pure fluid thermophysical property evaluation and the open-source thermophysical property library coolprop. Ind. Eng. Chem. Res. 2014, 53, 2498–2508. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. CoolProp: Fluid Properties for the Masses—CoolProp 4.2.5 Documentation. Available online: http://www.coolprop.org/ (accessed on 18 August 2014).
  16. Giacomo, P. Equation for the determination of the density of moist air. Metrologia 1982, 18. [Google Scholar] [CrossRef]
  17. Amann, H.; Escher, J. Analysis 2. Deutsche Ausgabe; Birkhäuser: Basel, Switzerland, 1999. [Google Scholar]
  18. Bommes, L. Ventilatoren; Vulkan-Verlag GmbH: Essen, Germany, 2003. [Google Scholar]
  19. European Parliament and Council. European Commision Regulation (EU) No 327/2011. Series L 90; Official Journal of the European Union, 2011; pp. 8–21. Available online: http://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32011R0327&from=en (accessed on 13 January 2015).
  20. Jiang, H.; Aute, V.; Radermacher, R. CoilDesigner: A general-purpose simulation and design tool for air-to-refrigerant heat exchangers. Int. J. Refrig. 2006, 29, 601–610. [Google Scholar] [CrossRef]
  21. SolCoolSys. Available online: https://www.solcoolsys.de/ (accessed on 29 January 2014).
  22. Thermofin. Available online: http://www.thermofin.de/ (accessed on 4 June 2014).
  23. Güntner GmbH & Co. KG: Product Calculator (GPC). Available online: http://www.guentner.de/know-how/product-calculator-gpc/ (accessed on 4 June 2014).
  24. Publicly Funded Project: Latent Heat Storage for Solar Heating and Cooling with Absorption Chiller—System Design, ZAE Bayern, Funding number FKZ 0329605O. 2013. Available online: http://www.zae-bayern.de/uploads/media/tb_2013.pdf (accessed on 13 January 2015).
  25. Khan, J.-U.-R.; Yaqub, M.; Zubair, S.M. Performance characteristics of counter flow wet cooling towers. Energy Convers. Manag. 2003, 44, 2073–2091. [Google Scholar] [CrossRef]
  26. Braun, J.E.; Klein, S.A.; Mitchell, P. Effectiveness models for cooling towers and cooling coils. ASHRAE Trans. 1989, 95, 154–174. [Google Scholar]
  27. Kloppers, J.C.; Kröger, D.G. Cooling tower performance: A critical evaluation of the merkel assumptions. R & D J. 2004, 20, 24–29. [Google Scholar]
  28. American Society of Heating, Refrigerating and Air-Conditioning Engineers ASHRAE Handbook 2005: Fundamentals. ASHRAE: Atlanta, GA, USA, 2005. Available online: https://www.ashrae.org/resources--publications/handbook (accessed on 13 January 2015).
  29. Braun, J.E. Methodologies for the Design and Control of Central Cooling Plants. Ph.D. Thesis; University of Wisconsin-Madison: Madison, WI, USA, 1988. Available online: http://minds.wisconsin.edu/handle/1793/46694 (accessed on 13 January 2015).
  30. TRNSYS—A Transient Systems Simulation Program, Official Website. 2013. Available online: http://sel.me.wisc.edu/trnsys/ (accessed on 13 January 2015).
  31. TRNSYS 17. Mathematical Reference. University of Wisconsin-Madison: Madison, WI, USA, 2009, 4th ed. Available online: web.mit.edu/parmstr/Public/TRNSYS/04-MathematicalReference.pdf (accessed on 13 January 2015).
  32. IEA SHC || Task 48. Available online: http://task48.iea-shc.org/ (accessed on 25 September 2014).
  33. SolaRück. Available online: https://www.solarueck.de/ (accessed on 25 September 2014).

Share and Cite

MDPI and ACS Style

Fugmann, H.; Nienborg, B.; Trommler, G.; Dalibard, A.; Schnabel, L. Performance Evaluation of Air-Based Heat Rejection Systems. Energies 2015, 8, 714-741. https://doi.org/10.3390/en8020714

AMA Style

Fugmann H, Nienborg B, Trommler G, Dalibard A, Schnabel L. Performance Evaluation of Air-Based Heat Rejection Systems. Energies. 2015; 8(2):714-741. https://doi.org/10.3390/en8020714

Chicago/Turabian Style

Fugmann, Hannes, Björn Nienborg, Gregor Trommler, Antoine Dalibard, and Lena Schnabel. 2015. "Performance Evaluation of Air-Based Heat Rejection Systems" Energies 8, no. 2: 714-741. https://doi.org/10.3390/en8020714

APA Style

Fugmann, H., Nienborg, B., Trommler, G., Dalibard, A., & Schnabel, L. (2015). Performance Evaluation of Air-Based Heat Rejection Systems. Energies, 8(2), 714-741. https://doi.org/10.3390/en8020714

Article Metrics

Back to TopTop