Cooling Towers and Evaporative Fluid Coolers[LINK]
One, Two, and Variable Speed Cooling Towers and Evaporative Fluid Coolers[LINK]
The input objects CoolingTower:SingleSpeed, CoolingTower:TwoSpeed, and CoolingTower:VariableSpeed:Merkel provide models for singlespeed, twospeed, and variablespeed cooling towers that are based on Merkel’s theory (Merkel 1925), which is also the basis for the tower model included in ASHRAE’s HVAC1 Toolkit for primary HVAC system energy calculations (ASHRAE 1999, Bourdouxhe et al. 1994). Cooling tower performance is modeled using effectivenessNTU relationships for counterflow heat exchangers. The model can be used to simulate the performance of both single speed, two speed, and variable speed mechanicaldraft cooling towers. The model will also account for tower performance in the “free convection” regime, when the tower fan is off but the water pump remains on. For partload operation, the model assumes a simple linear interpolation between two steadystate regimes without accounting for any cycling losses.
For single speed cooling towers, the capacity control can be fan cycling or fluid bypass. In fluid bypass mode, portion of the water goes through the tower media and gets cooled while the remaining water flow gets bypassed, two water flows then mix together trying to meet the tower exiting water setpoint temperature. In both the free convection cooling when fan is off and normal cooling when fan is on for the entire time step, if the tower exiting water temperature is lower than the setpoint, the tower operates in fluid bypass mode. The model determines the fluid bypass fraction by iterations until the mixed water meets the tower exiting water temperature setpoint. In the fluid bypass mode, except the free convection, the tower fan runs at full speed for the entire time step. The maximum amount of tower water that can be bypassed is bounded by the freezing point of the tower water  the tower exiting water temperature cannot be lower than the freezing setpoint.
Evaporative fluid coolers are modeled very similar to cooling towers. The main difference between the two is in the “Performance input method” input field. Cooling tower has two choices for this field namely “UFactorTimesAreaAndDesignWaterFlowRate” and “Nominal capacity”. The nominal capacity is specified for the standard conditions i.e. entering water at 35C (95F), leaving water at 29.44C (85F), entering air at 25.56C (78F) wetbulb temperature and 35C (95F) drybulb temperature. On the other hand evaporative fluid cooler has three choices for “Performance input method” which are “UFactorTimesAreaAndDesignWaterFlowRate”, “StandardDesignCapacity” and “UserSpecifiedDesignCapacity”. First method is same for both tower and fluid cooler. Standard design capacity is specified for the same conditions which are used to specify nominal capacity for tower as described above. If the capacity of fluid cooler for conditions other than the standard ones is known then UserSpecifiedDesignCapacity method should be used. In this case, the conditions for which the fluid cooler capacity is known i.e. entering water temperature, entering air temperature and entering air wet bulb temperature must be specified in the input. To calculate evaporation loss for fluid cooler, spray water flow rate which is different than the process fluid flow rate must be specified for all the performance input methods. This is not required for cooling tower because cooled fluid i.e. water is in direct contact with the air so the water loss is calculated by using cooled fluid flow rate only. Unlike cooling tower, evaporative fluid cooler model does not account for free convection.
Cooling tower model is described below which holds equally good for evaporative fluid cooler. The differences are mentioned whenever required.
Model Description[LINK]
Based on Merkel’s theory, the steadystate total heat transfer between the air and water entering the tower can be defined by the following equation:
d˙Qtotal=UdAcp(hs−ha)
where
h_{s} = enthalpy of saturated air at the wettedsurface temperature, J/kg
h_{a} = enthalpy of air in the free stream, J/kg
c_{p} = specific heat of moist air, J/kg^{o}C
U = cooling tower overall heat transfer coefficient, W/m^{2}^{o}C
A = heat transfer surface area, m^{2}
Equation is based on several assumptions:
air and water vapor behave as ideal gases
the effect of water evaporation is neglected
fan heat is neglected
the interfacial air film is assumed to be saturated
the Lewis number is equal to 1
In this model, it is also assumed that the moist air enthalpy is solely a function of the wetbulb temperature and that the moist air can be treated as an equivalent ideal gas with its mean specific heat defined by the following equation:
¯cpe=ΔhΔTwb
where
Δh = enthalpy difference between the air entering and leaving the tower, J/kg
ΔT_{wb} = wetbulb temperature difference between the air entering and leaving the tower, ^{o}C
Since the liquid side conductance is much greater than the gas side conductance, the wettedsurface temperature is assumed to be equal to the water temperature. Based on this assumption and equations and , the expression for total heat transfer becomes:
d˙Qtotal=UedA(Tw−Twb)
where
Ue=U¯cpecp
T_{wb} = wetbulb temperature of the air, ^{o}C
T_{w} = temperature of the water, ^{o}C
An energy balance on the water and air sides of the air/water interface yields the following equations:
d˙Qtotal=˙mwcpwdTw
d˙Qtotal=˙ma¯cpedTwb
where
˙mw = mass flow rate of water, kg/s
˙ma = mass flow rate of air, kg/s
Assuming that the heat capacity rate ( ˙mcp ) for the cooling tower water is less than that for the air, the effectiveness of the cooling tower can be defined by analogy to the effectiveness of a simple heat exchanger:
ε=Twin−TwoutTwin−Twbin
where
ε = heat exchanger effectiveness
T_{win} = inlet water temperature, ^{o}C
T_{wout} = outlet water temperature, ^{o}C
T_{wbin} = wetbulb temperature of the inlet air, ^{o}C
Combining equations , , and and integrating over the entire heat transfer surface area, and combining the result with equation provides the following expression for cooling tower effectiveness:
ε=1−exp{−NTU[1−(˙Cw˙Ca)]}1−(˙Cw˙Ca)exp{−NTU[1−(˙Cw˙Ca)]}
where
˙Cw=˙mwcpw and ˙Ca=˙ma¯cpe
NTU=NumberofTransferUnits=UAe˙Cw
This equation is identical to the expression for effectiveness of an indirect contact (i.e., fluids separated by a solid wall) counterflow heat exchanger (Incropera and DeWitt 1981). Therefore, the cooling tower can be modeled, in the steadystate regime, by an equivalent counterflow heat exchanger as shown in the following figure.
The first fluid is water and the second fluid is an equivalent fluid entering the heat exchanger at temperature T_{wbin} and specific heat ¯cpe . The heat exchanger is characterized by a single parameter, its overall heat transfer coefficientarea product UA_{e}. The actual cooling tower heat transfer coefficientarea product is related to UA_{e} by the following expression:
UA=UAecp¯cpe
This heat transfer coefficientarea product is assumed to be a function of the air mass flow rate only and can be estimated from laboratory test results or manufacturers’ catalog data.
The model for the variable speed Merkel tower also includes Scheier’s modifications. Scheier has extended the Merkel model to also include terms that adjust UA with three factors that model how UA values change when the tower is operating away from its rated conditions. The first factor, fUA,wetbulb , adjusts UA for the current outdoor wetbulb temperature. The user enters a performance curve or lookup table that is a function of one independent variable. The independent variable is the difference between the design wetbulb temperature and the current wetbulb temperature, in degrees Celsius.
fUA,wetbulb=func((Tdes,wetbulb−Tcurrent,wetbulb))
The second factor, fUA,airflow , adjusts UA for the current air flow rate. The user enters a performance curve or lookup table that is a function of one independent variable. The independent variable is the ratio of the current air flow rate to the design air flow rate at full speed.
fUA,airflow=func(AirFlowRatio)
The third factor, fUA,waterflow , adjusts UA for the current water flow rate. The user enters a performance curve or lookup table that is a function of one independent variable. The independent variable is the ratio of the current water flow rate to the design water flow rate.
fUA,waterflow=func(WaterFlowRatio)
Then the UA value at any given time is calculated using
UAe=UAdesign∙fUA,wetbulb∙fUA,airflow∙fUA,waterflow
Method for Calculating SteadyState Exiting Water Temperature[LINK]
The objective of the cooling tower model is to predict the exiting water temperature and the fan power required to meet the exiting water setpoint temperature. Since only the inlet air and inlet water temperatures are known at any simulation time step, an iterative procedure is required to determine the exiting fluid temperatures using the equations defined in the previous section. In the case of the EnergyPlus model, the iterations are performed to determine the exiting wetbulb temperature of the air. The exiting water temperature is then calculated based on an energy balance that assumes that the energy absorbed by the air is equivalent to the energy removed from the water. The procedure for calculating the steadystate, exiting air wetbulb temperature is outlined below.
As explained previously, it is assumed that the moist air enthalpy can be defined by the wetbulb temperature alone. Therefore, the first step in the procedure is to calculate the enthalpy of moist air entering the cooling tower based on the ambient wetbulb temperature from the weather file. Since an iterative solution is required, a first guess of the outlet air wetbulb temperature is then made and the enthalpy of this estimated outlet air wetbulb temperature is calculated. Based on these inlet and outlet air conditions, the mean specific heat of the air is calculated based on equation , repeated here:
¯cpe=ΔhΔTwb
With the overall heat transfer coefficientarea product for the cooling tower entered by the user, the effective heat transfer coefficientarea product is calculated by rearranging equation :
UAe=UA¯cpecp
With ¯cpe and UAe known, the effectiveness of the heat exchanger is then calculated:
ε=1−exp{−NTU[1−(˙Cmin˙Cmax)]}1−(˙Cmin˙Cmax)exp{−NTU[1−(˙Cmin˙Cmax)]}
where
˙Cmin=Minimum(˙Cw,˙Ca) and ˙Cmax=Maximum(˙Cw,˙Ca)
˙Cw=˙mwcpw and ˙Ca=˙ma¯cpe
NTU=NumberofTransferUnits=UAe˙Cmin
The heat transfer rate is then calculated as follows:
˙Qtotal=ε˙Cmin(Twin−Twbin)
The outlet air wetbulb temperature is then recalculated:
Twbout=Twbin+˙Qtotal˙Ca
The iterative process of calculating Twbout continues until convergence is reached.
Finally, the outlet water temperature is calculated as follows:
Twout=Twin+˙Qtotal˙Cw
Calculating the Actual Exiting Water Temperature and Fan Power[LINK]
The previous section describes the methodology used for calculating the steadystate temperature of the water leaving the cooling tower. This methodology is used to calculate the exiting water temperature in the free convection regime (water pump on, tower fan off) and with the tower fan operating (including low and high fan speed for the twospeed tower). The exiting water temperature calculations use the fluid flow rates (water and air) and the UAvalues entered by the user for each regime.
The cooling tower model seeks to maintain the temperature of the water exiting the cooling tower at (or below) a setpoint. The model obtains the target temperature setpoint from the setpoints placed on either the tower outlet node or the loop’s overall setpoint node (typically set to the supply side outlet node). The model checks to see if the outlet node has a setpoint placed on it and uses that if it does. If the outlet node does not have a temperature setpoint then the model uses the looplevel outlet node specified in the input field called Loop Temperature Setpoint Node Name in the PlantLoop or CondenserLoop object. The model first checks to determine the impact of “free convection”, if specified by the user, on the tower exiting water temperature. If free convection is not specified by the user, then the exiting water temperature is initially set equal to the entering tower water temperature. If the user specifies “free convection” and the steadystate exiting water temperature based on “free convection” is at or below the setpoint, then the tower fan is not turned on.
If the exiting water temperature remains above the setpoint after “free convection” is modeled, then the tower fan is turned on to reduce the exiting water temperature to the setpoint. The model assumes that partload operation is represented by a simple linear interpolation between two steadystate regimes (e.g., tower fan on for the entire simulation time step and tower fan off for the entire simulation time step). Cyclic losses are not taken into account.
The fraction of time that the tower fan must operate is calculated based on the following equation:
ω=Tset−Twout,offTwout,on−Twout,off
where
Tset = exiting water setpoint temperature, ^{o}C
Twout,off = exiting water temperature with tower fan off, ^{o}C
Twout,on = exiting water temperature with tower fan on, ^{o}C
The average fan power for the simulation time step is calculated by multiplying ω by the steadystate fan power specified by the user.
The calculation method for the twospeed tower is similar to that for the singlespeed tower example described above. The model first checks to see if “free convection” is specified and if the resulting exiting water temperature is below the setpoint temperature. If not, then the model calculates the steadystate exiting water temperature with the tower fan at low speed. If the exiting water temperature at low fan speed is below the setpoint temperature, then the average fan power is calculated based on the result of equation and the steadystate, low speed fan power specified by the user. If lowspeed fan operation is unable to reduce the exiting water temperature below the setpoint, then the tower fan is increased to its high speed and the steadystate exiting water temperature is calculated. If this temperature is below the setpoint, then a modified version of equation is used to calculate runtime at high fan speed:
ω=Tset−Twout,lowTwout,high−Twout,low
where
Tset = exiting water setpoint temperature, ^{o}C
Twout,low = exiting water temperature with tower fan at low speed, ^{o}C
Twout,high = exiting water temperature with tower fan at high speed, ^{o}C
The average fan power for the simulation time step is calculated for the twospeed cooling tower as follows:
Pfan,avg=ω(Pfan,high)+(1−ω)(Pfan,low)
The calculation method for the variablespeed Merkel/Scheier model is different from the one and twospeed models. Whereas the other towers are setpoint driven, the variable speed Merkel/Scheier model is driven by the load dispatched from the plant operation scheme, such as PlantEquipmentOperation:CoolingLoad. The plant operation scheme provides the cooling tower model with a target load to meet and the tower is controlled to meet that load (as opposed to meeting a leaving setpoint). This allows complex configurations with more than one cooling tower to be properly dispatched using all the features available in the various plant operation schemes. The model first evaluates the load provided with the fans off using only free convection. If the free convection mode meets or exceeds the target load then the tower runs in free convection mode. The model then evaluates the load provided by running the fans at full speed. If the full fan speed load provided is less than or equal to the target load then the tower runs at full speed. The model then evaluates the load provided with the variable speed fan running at the minimum air speed ratio, which is a user input (default is 0.2). If the minimum fan speed load provided meets or exceeds the target load, then tower runs at the minimum fan speed. If the target load is between the load provided by minimum and maximum fans speeds, then the model solves for an airflow rate ratio that just meets the target load using Regula Falsi numerical method. The variable speed Merkel/Scheier model does not model bypass and may provided excess cooling at times when running in free convection or at the minimum fan speed. Fan power for the variable speed Merkel/Scheier model is then calculated from the airflow rate ratio using a design fan power and a performance curve or lookup table with one independent variable for the airflow ratio.
Pfan=Pfan,des∙func(AirFlowRatio)
Tower Basin Heater[LINK]
The tower basin heater operates in the same manner as the variable speed cooling tower basin heater. Refer to the variable speed cooling tower basin heater description in the following section.
Cooling Tower Makeup Water Usage[LINK]
The cooling tower makeup water usage is the same as the variable speed cooling tower makeup water usage. Refer to the variable speed cooling tower makeup water usage description in the following section.
Rosaler, Robert C. 1995. Standard Handbook of Plant Engineering, 2^{nd} Ed. New York, NY: McGrawHill, pp. 63637.
Variable Speed Cooling Towers Empirical Models[LINK]
The input object CoolingTower:VariableSpeed provides models for variable speed towers that are based on empirical curve fits of manufacturer’s performance data or field measurements. The user specifies tower performance at design conditions, and empirical curves are used to determine the approach temperature and fan power at offdesign conditions. The user defines tower performance by entering the inlet air wetbulb temperature, tower range, and tower approach temperature at the design conditions. The corresponding water flow rate, air flow rate, and fan power must also be specified. The model will account for tower performance in the “free convection” regime, when the tower fan is off but the water pump remains on and heat transfer still occurs (albeit at a low level). Basin heater operation and makeup water usage (due to evaporation, drift, and blowdown) are also modeled.
The cooling tower seeks to maintain the temperature of the water exiting the cooling tower at (or below) a setpoint. The setpoint temperature is defined by the setpoints placed on either the tower outlet node or the loop’s overall setpoint node (typically set to the supply side outlet node). The model checks to see if the outlet node has a setpoint placed on it and uses that if it does. If the outlet node does not have a temperature setpoint then the model uses the looplevel outlet node specified in the input field called Loop Temperature Setpoint Node Name in the PlantLoop or CondenserLoop object. The model simulates the outlet water temperature in four successive steps:
The model first determines the tower outlet water temperature with the tower fan operating at maximum speed. If the outlet water temperature is above the setpoint temperature, the fan runs at maximum speed.
If the outlet water temperature with maximum fan speed is below the setpoint temperature, then the model next determines the impact of “free convection” (water flowing through tower with fan off). If the exiting water temperature based on “free convection” is at or below the setpoint, then the tower fan is not turned on.
If the outlet water temperature remains above the setpoint after “free convection” is modeled, then the tower fan is turned on at the minimum fan speed (minimum air flow rate ratio) to reduce the leaving water temperature. If the outlet water temperature is below the setpoint at minimum fan speed, the tower fan is cycled on and off to maintain the outlet water setpoint temperature.
If the outlet water temperature remains above the setpoint after minimum fan speed is modeled, then the tower fan is turned on and the model determines the required air flow rate and corresponding fan speed to meet the desired setpoint temperature.
Model Description[LINK]
The variable speed tower model utilizes userdefined tower performance at design conditions along with empirical curves to determine tower heat rejection and fan power at offdesign conditions. Basin heater operation and makeup water usage are also modeled based on user inputs, tower entering air conditions, and tower operation. The following sections describe how each of these tower performance areas is modeled.
Tower Heat Rejection[LINK]
Heat rejection by the variable speed cooling tower is modeled based on the CoolTools correlation, YorkCalc correlation, or userdefined coefficients for either the CoolTools or YorkCalc correlations. These purelyempirical correlations model the tower approach temperature using a polynomial curve fit with a large number of terms and either three or four independent variables.
The CoolTools correlation has 35 terms with four independent variables:
Approach = Coeff(1) + Coeff(2)•FRair + Coeff(3)•(FRair)^{2} +
Coeff(4)•(FRair)^{3} + Coeff(5)•FRwater +
Coeff(6)•FRair•FRwater + Coeff(7)•(FRair)^{2}•FRwater +
Coeff(8)•(FRwater)^{2} + Coeff(9)•FRair•(FRwater)^{2} +
Coeff(10)•(FRwater)^{3} + Coeff(11)•Twb + Coeff(12)•FRair•Twb +
Coeff(13)•(FRair)^{2}•Twb + Coeff(14)•FRwater•Twb +
Coeff(15)•FRair•FRwater•Twb + Coeff(16)•(FRwater)^{2}•Twb +
Coeff(17)•(Twb)^{2} + Coeff(18)•FRair•(Twb)^{2} +
Coeff(19)•FRwater•(Twb)^{2} + Coeff(20)•(Twb)^{3} + Coeff(21)•Tr +
Coeff(22)•FRair•Tr + Coeff(23)•(FRair)^{2}•Tr +
Coeff(24)•FRwater•Tr + Coeff(25)•FRair•FRwater•Tr +
Coeff(26)•(FRwater)^{2}•Tr + Coeff(27)•Twb•Tr +
Coeff(28)•FRair•Twb•Tr + Coeff(29)•FRwater•Twb•Tr +
Coeff(30)•(Twb)^{2}•Tr + Coeff(31)•(Tr)^{2}+ Coeff(32)•FRair•(Tr)^{2}+
Coeff(33)•FRwater•(Tr)^{2}+ Coeff(34)•Twb•(Tr)^{2}+ Coeff(35)•(Tr)^{3}
where:
Approach = approach temperature (^{o}C) = outlet water temperature minus inlet air wetbulb temperature
FRair = air flow rate ratio (actual air flow rate divided by design air flow rate)
FRwater = water flow rate ratio (actual water flow rate divided by design water flow rate)
Tr = range temperature (^{o}C) = inlet water temperature minus outlet water temperature
Twb = inlet air wetbulb temperature (^{o}C)
Coeff(#) = correlation coefficients
If the user selects Tower Model Type = CoolToolsCrossFlow, then the 35 coefficients derived for the CoolTools simulation model (Benton et al. 2002) are used and these coefficients are already defined within EnergyPlus as shown in Table 63. If the user specifies Tower Model Type = CoolToolsUserDefined, then the user must enter a CoolingTowerPerformance:CoolTools object to define the 35 coefficients that will be used by the CoolTools approach temperature correlation.
Approach Temperature Correlation Coefficients
Coeff(1)

0.52049709836241

0.359741205

Coeff(2)

10.617046395344

0.055053608

Coeff(3)

10.7292974722538

0.0023850432

Coeff(4)

2.74988377158227

0.173926877

Coeff(5)

4.73629943913743

0.0248473764

Coeff(6)

8.25759700874711

0.00048430224

Coeff(7)

1.57640938114136

0.005589849456

Coeff(8)

6.51119643791324

0.0005770079712

Coeff(9)

1.50433525206692

0.00001342427256

Coeff(10)

3.2888529287801

2.84765801111111

Coeff(11)

0.0257786145353773

0.121765149

Coeff(12)

0.182464289315254

0.0014599242

Coeff(13)

0.0818947291400898

1.680428651

Coeff(14)

0.215010003996285

0.0166920786

Coeff(15)

0.0186741309635284

0.0007190532

Coeff(16)

0.0536824177590012

0.025485194448

Coeff(17)

0.00270968955115031

0.0000487491696

Coeff(18)

0.00112277498589279

0.00002719234152

Coeff(19)

0.00127758497497718

0.0653766255555556

Coeff(20)

0.0000760420796601607

0.002278167

Coeff(21)

1.43600088336017

0.0002500254

Coeff(22)

0.5198695909109

0.0910565458

Coeff(23)

0.117339576910507

0.00318176316

Coeff(24)

1.50492810819924

0.000038621772

Coeff(25)

0.135898905926974

0.0034285382352

Coeff(26)

0.152577581866506

0.00000856589904

Coeff(27)

0.0533843828114562

0.000001516821552

Coeff(28)

0.00493294869565511

N/A

Coeff(29)

0.00796260394174197

N/A

Coeff(30)

0.000222619828621544

N/A

Coeff(31)

0.0543952001568055

N/A

Coeff(32)

0.00474266879161693

N/A

Coeff(33)

0.0185854671815598

N/A

Coeff(34)

0.00115667701293848

N/A

Coeff(35)

0.000807370664460284

N/A

: Approach Temperature Correlation Coefficients
Similarly, the YorkCalc correlation has 27 terms with three independent variables:
Approach = Coeff(1) + Coeff(2)•Twb + Coeff(3)•Twb^{2} + Coeff(4)•Tr +
Coeff(5)•Twb•Tr + Coeff(6)•Twb^{2}•Tr + Coeff(7)•Tr^{2} +
Coeff(8)•Twb•Tr^{2}+ Coeff(9)•Twb^{2}•Tr^{2}+ Coeff(10)•LGRatio +
Coeff(11)•Twb•LGRatio + Coeff(12)•Twb^{2}•LGRatio +
Coeff(13)•Tr•LGRatio + Coeff(14)•Twb•Tr•LGRatio +
Coeff(15)•Twb^{2}•Tr•LGRatio + Coeff(16)•Tr^{2}•LGRatio +
Coeff(17)•Twb•Tr^{2}•LGRatio + Coeff(18)•Twb^{2}•Tr^{2}•LGRatio +
Coeff(19)•LGRatio^{2}+ Coeff(20)•Twb•LGRatio^{2} +
Coeff(21)• Twb^{2}•LGRatio^{2}+ Coeff(22)•Tr•LGRatio^{2}+
Coeff(23)•Twb•Tr•LGRatio^{2}+ Coeff(24)•Twb^{2}•Tr•LGRatio^{2}+
Coeff(25)•Tr^{2}•LGRatio^{2}+ Coeff(26)•Twb•Tr^{2}•LGRatio^{2}+
Coeff(27)•Twb^{2}•Tr^{2}•LGRatio^{2}
where:
Approach = approach temperature (^{o}C) = outlet water temperature minus inlet air wetbulb temperature
Tr = range temperature (^{o}C) = inlet water temperature minus outlet water temperature
Twb = inlet air wetbulb temperature (^{o}C)
LGratio = liquidtogas ratio = ratio of water flow rate ratio (FRwater) to air flow rate ratio (FRair)
Coeff(#) = correlation coefficients
If the user selects Tower Model Type = YorkCalc, then the 27 coefficients derived for the YorkCalc simulation model (York International Corp. 2002) are used and these coefficients are already defined within EnergyPlus as shown in Table 63. If the user specifies Tower Model Type = YorkCalcUserDefined, then the user must enter a CoolingTowerPerformance:YorkCalc object to define the 27 coefficients that will be used by the YorkCalc approach temperature correlation.
The approach temperature correlations for the CoolTools and YorkCalc simulation models are valid for a range of conditions defined in Table 64. If the user defines their own model coefficients (CoolingTowerPerformance:CoolTools or CoolingTowerPerformance:YorkCalc), then they must also define in that same object the range of conditions for which the model is valid. For all of these correlation variables, the program issues warnings if the actual values are beyond the minimum/maximum values specified for the correlation being used. For inlet air wetbulb temperature and water mass flow rate ratio, the values of these variables used in the calculation of approach temperature are limited to be within the valid minimum/maximum range. For approach, range, and liquidtogas ratio the warnings are issued if the values are beyond the specified minimum/maximum range but the actual values are still used. The warnings issued do not necessarily indicate a poor estimate of tower performance at the condition(s) which caused the warning, but are provided to identify conditions outside the defined correlation limits. Exceeding the defined limits by a small amount may not introduce significant errors, but large deviations may be problematic. It is for this reason that we recommend using a very broad range of cooling tower performance data (i.e., data covering the entire range expected during the simulation) when generating userdefined coefficients for the variable speed tower model.
Minimum and Maximum Limits for Approach Temperature Correlation Variables
Minimum Inlet Air WetBulb Temperature 
1.0°C 
34.4°C 
Maximum Inlet Air WetBulb Temperature 
26.7°C 
26.7°C 
Minimum Tower Range Temperature 
1.1°C 
1.1°C 
Maximum Tower Range Temperature 
11.1°C 
22.2°C 
Minimum Tower Approach Temperature 
1.1°C 
1.1°C 
Maximum Tower Approach Temperature 
11.1°C 
40°C 
Minimum Water Flow Rate Ratio 
0.75 
0.75 
Maximum Water Flow Rate Ratio 
1.25 
1.25 
Maximum LiquidtoGas Ratio 
N/A 
8.0 
The approach temperature correlation(s) used to simulate cooling tower heat rejection are based on water and air flow rate “ratios” and are not directly dependent on the size of the tower or the actual air and water flow rates through the tower. However, the model correlations are developed based on a reference condition. For Model Types “CoolToolsCrossFlow” and “YorkCalc”, the reference condition is a water flow rate of 0.000043 m^{3}/s per kW of heat rejected (2.4 gal/min per ton of heat rejected) with 25.6^{o}C (78^{o}F) enter air wetbulb temperature, 35^{o}C (95^{o}F) hot water inlet temperature, and 29.4^{o}C (85^{o}F) cold water outlet temperature. The reference condition may be different if the user defines tower model coefficients using CoolingTowerPerformance:CoolTools or CoolingTowerPerformance:YorkCalc.
Due to the inherent reference condition used to generate the tower performance curves, the water flow rate at the reference condition must be determined using the design performance information specified by the user and the tower model’s approach temperature correlation. This is done by using the model’s approach temperature correlation (described earlier in this section) to calculate the water flow rate ratio which yields the userdefined design approach temperature based on an air flow rate ratio of 1.0 (FR_{air} = 1.0), the design inlet air wetbulb temperature, and the design range temperature. The calculated approach temperature (using the model correlation) must satisfy the following two equations:
Twater,outlet,design=Twb,air,design+Tapproach,design{FRair=1.0} Twater,inlet,design=Twb,air,design+Tapproach,design+Trange,design{FRair=1.0}
where:
Twater,outlet,design = design outlet water temperature (^{o}C)
Twater,inlet,design = design inlet water temperature (^{o}C)
Twb,air,design = design inlet air wetbulb temperature (^{o}C)
Tapproach,design = design approach temperature (^{o}C)
Trange,design = design range temperature (^{o}C)
FRair = air flow rate ratio (actual air flow rate divided by design air flow rate)
The water flow rate ratio used in the approach temperature correlation which satisfies these two equations is the ratio of the design water flow rate (specified by the user) to the water flow rate at the reference condition. This ratio is used to calculate the reference water volumetric flow rate, which is then used throughout the simulation to determine the actual water flow rate ratio used in the approach temperature correlation for each simulation time step.
˙Vreference=˙VdesignFRwaterreference
where:
˙Vreference = water volumetric flow rate at the reference condition (m^{3}/s)
˙Vdesign = design water volumetric flow rate specified by the user (m^{3}/s)
FRwater,reference = design water flow rate divided by the reference water flow rate
The cooling tower seeks to maintain the temperature of the water exiting the cooling tower at (or below) a setpoint. The setpoint temperature is defined by the field “Condenser Loop Temperature Setpoint schedule or reference” for the CondenserLoop object. The model simulates the outlet water temperature in four successive steps:
The model first determines the tower outlet water temperature with the tower fan operating at maximum speed. If the outlet water temperature is above the setpoint temperature, the fan runs at maximum speed.
If the outlet water temperature with maximum fan speed is below the setpoint temperature, then the model next determines the impact of “free convection” (water flowing through tower with fan off). If the exiting water temperature based on “free convection” is at or below the setpoint, then the tower fan is not turned on.
If the outlet water temperature remains above the setpoint after “free convection” is modeled, then the tower fan is turned on at the minimum fan speed (minimum air flow rate ratio) to reduce the leaving water temperature. If the outlet water temperature is below the setpoint at minimum fan speed, the tower fan is cycled on and off to maintain the outlet water setpoint temperature.
If the outlet water temperature remains above the setpoint after minimum fan speed is modeled, then the tower fan is turned on and the model determines the required air flow rate and corresponding fan speed to meet the desired setpoint temperature.
For each simulation time step, the model first calculates the outlet water temperature with the tower fan operating at maximum speed (FRair = 1.0). The calculated approach temperature (using the correlations described above), inlet air wetbulb temperature (weather data), and range temperature are used to determine the tower outlet water temperature as follows:
Twater,outlet,fanMAX=Twb,air+Tapproach{FRair=1.0}=Twater,inlet+Trange
where:
Twater,outlet,fanMAX = tower outlet water temperature at maximum fan speed (^{o}C)
Twb,air = tower inlet air wetbulb temperature (^{o}C)
Tapproach = approach temperature at current operating conditions (^{o}C)
Trange = range temperature at current operating conditions (^{o}C)
Note that the approach temperature correlation as described previously is a function of range temperature, so the equations above must be solved iteratively to converge on a solution. If the resulting outlet water temperature is above the desired setpoint temperature, then the fan runs at maximum speed and does not cycle on/off (fan partload ratio = FanPLR = 1.0 and FR_{air} = 1.0).
If the outlet water temperature with maximum fan speed is below the setpoint temperature, then the model next determines the impact of “free convection” (water flowing through tower with fan off). In the free convection regime, the outlet water temperature is calculated using a fraction of the water temperature difference through the tower when the fan is at its maximum speed. This fraction is defined by the user (Fraction of Tower Capacity in Free Convection Regime).
Twater,outlet,fanOFF=Twater,inlet−[Fracfreeconv(Twater,inlet−Twater,outlet,fanMAX)]
where:
Twater,outlet,fanOFF = tower outlet water temperature in free convection regime (^{o}C)
Twater,inlet = tower inlet water temperature (^{o}C)
Fracfreeconv = fraction of tower capacity in free convection regime (user specified)
If the outlet water temperature in the free convection regime is below the setpoint temperature, the tower fan is not turned on and the fan partload ratio is set equal to 0. In addition, the air flow rate ratio through the tower is assumed to be equal to the fraction of tower capacity in the free convection regime.
FanPLR=FanPLRfreeconv=0.0
FRair=FRair,freeconv=Fracfreeconv
where:
FanPLR = fan partload ratio
FanPLRfreeconv = fan partload ratio in free convection regime
FRair,freeconv = air flow rate ratio in free convection regime
If the outlet water temperature in the free convection regime is above the setpoint temperature, then the fan is turned on at the minimum fan speed (minimum air flow rate ratio, FR_{air,min}, entered by the user) and the outlet water temperature is calculated as the inlet air wetbulb temperature plus the calculated approach temperature:
FRair=FRair,min
Twater,outlet,fanMIN=Twb,air+Tapproach{FRair=FRair,min}
where:
Twater,outlet,FanMIN = outlet water temperature at minimum fan speed (^{o}C)
FRair,min = air flow rate ratio at the minimum fan speed
If the outlet water temperature at minimum fan speed is below the setpoint temperature, the cooling tower fan cycles on and off at the minimum air flow rate ratio in order to meet the setpoint temperature.
FanPLR=(Twater,outlet,fanOFF−TsetpointTwater,outlet,fanOFF−Twater,outlet,fanMIN)
where:
Tsetpoint = outlet water setpoint temperature (^{o}C)
If the outlet water temperature at minimum fan speed is above the outlet water temperature setpoint, then the cooling tower fan speed (FR_{air}) is increased until the calculated approach temperature produces the required outlet water temperature to meet the setpoint.
Tapproach=f(FRair)
Twater,outlet=Tsetpoint=Twb,air+Tapproach
FanPLR=1.0 (i.e., fan does not cycle on/off)
When the cooling tower fan is operating, fan electric power is calculated based on the air flow rate ratio required to meet the above conditions. If the user has entered a fan power curve object (cubic curve), the output of that curve is multiplied by the design fan power. Otherwise, tower fan power is assumed to be directly proportional to the cube of the air flow rate ratio. In either case, the fan partload ratio is applied to account for times when the tower fan cycles on/off to meet the setpoint temperature. Fan energy consumption is calculated each simulation time step.
IF(FanPowerCurveObjectisDefined)THENPfan=(CurveOutput)(FanPLR)Pfan,designELSEPfan=(Pfan,design)(FanPLR)(FRair)3ENDIF
Efan=(Pfan)(TimeStepSys)3600.
where:
FanPowerCurveObject = name of fan power ratio as a function of air flow rate ratio curve
Pfan = tower fan electric power (W)
Efan = tower fan electric consumption (J)
CurveOutput = output of FanPowerCurveObject evaluated at the operating air flow rate ratio (FR_{air})
Pfan,design = design fan power at design (maximum) air flow through the tower (W)
TimeStepSys = HVAC system simulation time step (hr)
Tower Basin Heater[LINK]
Calculations are also made to estimate the electric power input to the tower basin heater. A schedule may be used to disable the basin heater during regular maintenance periods or other time periods (e.g., during summer). If a schedule is not provided, the basin heater is assumed to be available the entire simulation time period. The basin heater operates when it is scheduled on, the outdoor air drybulb temperature is below the basin heater setpoint temperature, and the cooling tower is not active (i.e., water is not flowing through the tower). The user is required to enter a basin heater capacity (watts per degree Kelvin) and a heater setpoint temperature (^{o}C) if they want to model basin heater electric power.
Pheater,basin=0.0IF(WaterNotFlowingThroughTower)THENIF(Scheduleheater,basinisDefined)THENIF(CAPheater,basin>0.AND.Scheduleheater,basin=ON)THENPheater,basin=MAX(0.0,CAPheater,basin(Tsetpoint,basin−Tdb,outdoor))ENDIFELSEIF(CAPheater,basin>0)THENPheater,basin=MAX(0.0,CAPheater,basin(Tsetpoint,basin−Tdb,outdoor))ENDIFENDIFENDIF
Eheater,basin=(Pheater,basin)(TimeStepSys)3600.
where:
Pheater,basin = tower basin heater electric power (W)
Eheater,basin = tower basin heater electric consumption (J)
T_{setpoint,basin} = basin heater setpoint temperature (^{o}C)
T_{db,outdoor} = outdoor air drybulb temperature (^{o}C)
CAP_{heater,basin} = basin heater capacity (W/K)
Schedule_{heater,basin} = basin heater schedule (schedule value > 0 means ON)
ASHRAE 1999. HVAC1 Toolkit: A Toolkit for Primary HVAC System Energy Calculations. Atlanta: American Society of Heating, Refrigerating and AirConditioning Engineers, Inc.
Benton, D.J., Bowman, C.F., Hydeman, M., Miller, P. 2002. An Improved Cooling Tower Algorithm for the CoolTools^{TM} Simulation Model. ASHRAE Transactions, Vol. 108, Part 1, pp.760768.
Bourdouxhe, J.P., M. Grodent, J. Lebrun and C. Silva. 1994. Cooling tower model developed in a toolkit for primary HVAC system energy calculation: part 1. Proceedings of the fourth international conference on system simulation in buildings, Liege (Belgium), December 57, 1994.
Incropera, F.P. and D.P. DeWitt. 1981. Fundamentals of Heat Transfer. New York: John Wiley & Sons.
Merkel, F. 1925. Verduftungskuhlung. VDI Forschungarbeiten, No 275, Berlin.
Rosaler, Robert C. 1995. Standard Handbook of Plant Engineering, 2^{nd} Ed. New York, NY: McGrawHill, pp. 63637.
Scheier, L. 2013. Personal communication.
York International Corporation, 2002. “YORKcalc^{TM} Software, ChillerPlant EnergyEstimating Program”, Form 160.00SG2 (0502).
Cooling Towers with Multiple Cells[LINK]
Many towers are constructed to be capable of being grouped together to achieve the desired capacity. Thus, many cooling towers are assemblies of two or more individual cooling towers or “cells.” The number of cells they have, e.g., an eightcell tower, often refers to such towers.
For the operation of multicell towers, the first step is to determine the number of cells n, which will be operating during the timestep using the calculation logic from DOE2.1E.
The maximum and minimum flow rates per cell are determined according to the input fractions (Minimum Water Flow Rate Fraction: kmin and Maximum Water Flow Rate Fraction: kmax) as follows:
⎧⎪⎨⎪⎩˙mw,min=˙mw,desntot×kmin˙mw,max=˙mw,desntot×kmax
where ˙mw,des is the design water flow rate through the entire cooling tower.
Then, we determine the minimum and maximum number of cells that can operate with this water flow rate:
⎧⎪
⎪⎨⎪
⎪⎩nmin=MIN(˙mw,tot˙mw,max,ntot)nmax=MIN(˙mw,tot˙mw,min,ntot)
where n_{tot} is the total number of cells of the tower, and mw,tot is the water flow rate to the tower.
The number of cells operating n is set accordingly:
If the Cell Control method is MinimalCell,
n=nmin
If the Cell Control method is MaximalCell,
n=nmax
Finally, the water mass flow rate per cell (˙mw) will be:
˙mw=˙mw,totn
Then we simulate the performance of one cell with this flow rate per cell (calling the SimSimpleTower subroutine for single and two speed cooling tower objects). As we assume that each cell is identical, the UA of one cell is calculated dividing the UA of the whole tower (obtained from the input or from the auto sizing calculations). The air flow rate per cell is also equal to the one of the whole tower divided by the number of cells operating:
UAcell=UAcellntot
˙mair,cell=˙mair,totn
At the end, the total fan power of the tower operating with a certain number of cells is given by:
FanPower=FanPowertot×nntot
If the cells operating do not meet the loads, we increase the number of cells if spare cells are available and the water flow through each cell is within the user specified minimum and maximum water flow rate fractions range. This is an iteration process.
Cooling Tower Makeup Water Usage[LINK]
Makeup water use for all types of cooling towers is made up of three components: evaporation, drift, and blowdown. The first is the amount of water evaporated to reduce the water’s temperature as it passes through the cooling tower. There are two methods that evaporation makeup water can be modeled in EnergyPlus. The first method assumes that the tower outlet air conditions are saturated (which may not always be the case for certain operating conditions). For this “Saturated Exit” mode, the enthalpy of the tower’s outlet air is calculated as the inlet air enthalpy plus the water side heat transfer divided by the air mass flow rate through the tower.
˙Qwater=˙mwaterCp(Twater,inlet−Twater,outlet)
hsatair,outlet=hair,inlet+(˙Qwater˙mair)
where:
˙Qwater = waterside heat transfer (W)
˙mwater = mass flow rate of water through the tower (kg/s)
Cp = specific heat of water (W/kgK)
hsatair,outlet = saturated outlet air enthalpy (J/kg)
hair,inlet = inlet air enthalpy (J/kg)
˙mair = mass flow rate of air through the tower (kg/s)
The saturation temperature and humidity ratio are then calculated for the tower’s outlet air.
Tsatair,outlet=PsyTsatFnHPb(hsatair,outlet,OutBaroPress)
ωsatair,outlet=PsyWFnTdbH(Tsatair,outlet,hsatair,outlet)
where:
Tsatair,outlet = saturated outlet air temperature (^{o}C)
PsyTsatFnHPb = EnergyPlus psychrometric function, returns saturation temperature given enthalpy and barometric pressure
OutBaroPress = outdoor barometric pressure (Pa)
ωsatair,outlet = saturated outlet air humidity ratio (kg_{water}/kg~dry air~)
PsyWFnTdbH = EnergyPlus psychrometric function, returns humidity ratio given drybulb temperature and enthalpy
The makeup water quantity required to replenish the water lost due to evaporation is then calculated as the product of the air mass flow rate and the difference between the entering and leaving air humidity ratio divided by the density of water.
˙Vevaporation,makeup=˙mair(ωsatair,outlet−ωair,inlet)ρwater
where:
˙Vevaporation,makeup = makeup water usage due to evaporation (m^{3}/s)
˙mair = mass flow rate of air through tower (kg/s)
ωair,inlet = humidity ratio of tower inlet air (kg_{water}/kg~dry air~)
ρwater = density of water evaluated at the tower inlet air temperature (kg/m^{3})
The second method available for calculating water makeup for evaporation is for the user to provide a value for a loss factor. The evaporation loss is then calculated as a fraction of the circulating condenser water flow and varies with the temperature change in the condenser water. The value provided by the user is in units of percentperdegree Kelvin. The evaporation rate will equal this value times each degree Kelvin of temperature drop in the condenser water. Typical values are from 0.15 to 0.27 [percent/K]. The default is 0.2. The rate of water makeup for evaporation is then calculated by multiplying this factor times the condenser water flow rate and the temperature decrease in the condenser water flow rate. For evaporative fluid coolers, a numerical value of loss factor can be entered in the same manner as for cooling towers. If this field is blank, an empirical correlation will be used to calculate the value based on current outdoor dry bulb temperature and relative humidity. The following correlation from Qureshi and Zubair (2007) is used to calculate the loss factor:
LossFactor=(113−8.417ϕ+1.6147tdb)∗10−5
where:
ϕ = relative humidity of inlet air
tdb = Drybulb temperature of inlet air
Additional makeup water usage is modeled as a percentage of design water flow rate through the tower to account for drift, and as a scheduled flow rate to model blowdown. Drift is water loss due to the entrainment of small water droplets in the air stream passing through the tower. Drift is defined by the model user as a percentage of the tower’s design water flow rate, and is assumed to vary with tower air flow rate ratio as follows:
˙Vdrift,makeup=˙mwater,design(PercentWatertoDrift100.0)FRair
where:
˙Vdrift,makeup = makeup water usage due to drift (m^{3}/s)
˙mwater,design = design (volumetric) water flow rate (m^{3}/s)
PercentWatertoDrift = percent of design water flow rate lost to drift at the tower design air flow rate
FRair = ratio of actual air flow rate to tower design air flow rate
Blowdown is water flushed from the basin on a periodic basis to purge the concentration of mineral scale or other contaminants. There are two ways that blowdown is calculated in EnergyPlus. Blowdown water rates can be scheduled so that we have:.
IF(Scheduleblowdownisdefined)THEN˙Vblowdown,makeup=ScheduleValueblowdownELSE˙Vblowdown,makeup=0.0ENDIF
where:
˙Vblowdown,makeup = makeup water usage due to blowdown (m^{3}/s)
ScheduleValue_{blowdown} = blowdown schedule value for the time step being simulated (m^{3}/s)
The second (and default) way that blowdown can be calculated is to assume that blowdown water is continually introduced at a rate that will provide a constant concentration ratio. As water evaporates it leaves behind minerals and the like causing the concentration of water impurities to be higher in the tower than in the makeup water. Acceptable concentration ratios are in the range of 3 to 5 depending on the purity of the make up water. Water lost as drift does not evaporate and decrease the water needed for blowdown. Using the “Concentration Ratio” method, the rate of blowdown can be calculated using:
˙Vblowdown,makeup=˙Vevaporation,makeupRconcentration−1−˙Vdrift,makeup
where,
Rconcentration is the concentration ratio or the ratio of solids in the blowdown water to solids in the makeup water.
The tower makeup water consumption (m^{3}) for each simulation time step is calculated as the sum of the individual components of makeup water usage multiplied by the simulation time step in hours and the conversion for hours to seconds (3600 sec/hr). Makeup water usage is only calculated when the cooling tower is active and water is flowing through the cooling tower.
Vwater,makeup=(˙Vevaporation,makeup+˙Vdrift,makeup+˙Vblowdown,makeup)TimeStepSys3600.
where:
Vwater,makeup = tower makeup water consumption (m^{3})
Robert C. Rosaler. 1995. Standard Handbook of Plant Engineering, 2^{nd} Ed. McGrawHill, New York, NY, pp 63637.
Quereshi, B.A. and S.M.Zubair. 2007. Prediction of evaporation losses in evaporative fluid coolers, Applied Thermal Engineering 27 pp. 520527
One and Two Speed Fluid Coolers[LINK]
The input objects FluidCooler:SingleSpeed and FluidCooler:TwoSpeed provide models for dry fluid coolers. Fluid cooler’s performance is modeled using effectivenessNTU relationships for cross flow heat exchanger with both streams unmixed. The model can be used to simulate the performance of both single speed and two speed mechanicaldraft fluid coolers. For partload operation, the model assumes a simple linear interpolation between two steadystate regimes without accounting for any cycling losses.
Model Description[LINK]
The expression for fluid cooler effectiveness is as follows:
ε=1−e[e(−NTUCrη)−1Crη]
Where
ε = heat exchanger effectiveness
Cr=˙Cmin˙Cmax
˙Ca=˙macpa Ca=macpaCa=macpa and ˙Cw=˙mwcpw
˙Cmax=Max(˙Ca,˙Cw) Cmax=Max (Ca,Cw)Cmax=Max(Ca,Cw); ˙Cmin=Min(˙Ca,˙Cw)
NTU=UA˙Cmin
η=NTU(−0.22) η=NTU0.22η=NTU0.22
The first fluid is water and the second fluid is air entering the heat exchanger at temperature Tdb,in TdbinTdbin and specific heat cpa cpacpa. The heat exchanger is characterized by a single parameter, its overall heat transfer coefficientarea product UA.
When the user selects the nominal capacity method, the UA is calculated as follows:
The model inputs (other than the UA) and the fluid cooler load that it must meet are specified at design conditions. Then the fluid cooler model converges to a UA value, using the regulafalsi method that will enable it to meet the design fluid cooler load given at the specified inputs.
Method for Calculating SteadyState Exiting Water Temperature[LINK]
The objective of the fluid cooler model is to predict the exiting water temperature and the fan power required to meet the exiting water setpoint temperature. The exiting water temperature is calculated based on an energy balance that assumes that the energy absorbed by the air is equivalent to the energy removed from the water. The procedure for calculating the steadystate, exiting air drybulb temperature is outlined below.
With the overall heat transfer coefficientarea product for the fluid cooler calculated by the nominal capacity information entered by the user, the effectiveness of the heat exchanger is then calculated as:
ε=1−e[e(−NTUCrη)−1Crη]
The heat transfer rate is then calculated as follows:
˙Q=ε˙Cmin(Tw,in−Tdb,in)
˙Q=ε˙Cmin(Tw,in−Tdbin)
Then the outlet air drybulb and outlet water temperature are calculated:
Tdb,out=Tdb,in+.Q.Ca
Tw,out=Tw,in−.Q.Cw
Tw,in = inlet water temperature, ^{o}C
Tw,out = outlet water temperature, ^{o}C
Tdb,in = drybulb temperature of the inlet air, ^{o}C
Tdb,out = drybulb temperature of the outlet air, ^{o}C
Calculating the Actual Exiting Water Temperature and Fan Power[LINK]
The previous section describes the methodology used for calculating the steadystate temperature of the water leaving the fluid cooler. This methodology is used to calculate the exiting water temperature with the fluid cooler fans operating (including low and high fan speed for the twospeed fluid cooler). The exiting water temperature calculations use the fluid flow rates (water and air) and the Nominal capacity information entered by the user for each regime.
The fluid cooler model seeks to maintain the temperature of the water exiting the fluid cooler at (or below) a setpoint. The setpoint schedule is defined by the field “Loop Temperature Setpoint Node or reference” for the CondenserLoop object.
The fluid cooler fans are turned on to reduce the exiting water temperature to the setpoint. The model assumes that partload operation is represented by a simple linear interpolation between two steadystate regimes (e.g., Fluid cooler fans on for the entire simulation time step and fluid cooler fans off for the entire simulation time step). Cyclic losses are not taken into account. If the outlet water temperature is less than the setpoint then the fraction of time for which the fluid cooler must operate to meet the setpoint is calculated by using the following equation:
ω=Tset−Twout,offTwout,on−Twout,off
Where
Tset = exiting water setpoint temperature, ^{o}C
Twout,off = exiting water temperature with all fluid cooler fans off, ^{o}C
Twout,on = exiting water temperature with all fluid cooler fans on, ^{o}C
The average fan power for the simulation time step is calculated by multiplying ω by the steadystate fan power specified by the user.
The calculation method for the twospeed fluid cooler is similar to that for the singlespeed fluid cooler example described above. The model first calculates the steadystate exiting water temperature with the fluid cooler fans at low speed. If the exiting water temperature at low fan speed is below the setpoint temperature, then the average fan power is calculated based on the result of previous equation and the steadystate, low speed fan power specified by the user. If lowspeed fan operation is unable to reduce the exiting water temperature below the setpoint, then the fluid cooler fans’ speed is increased to high speed and the steadystate exiting water temperature is calculated. If this temperature is below the setpoint, then a modified version of previous equation is used to calculate runtime at high fan speed:
ω=Tset−Twout,lowTwout,high−Twout,low
where
Tset = exiting water setpoint temperature, ^{o}C
Twout,low = exiting water temperature with fluid cooler fans at low speed, ^{o}C
Twout,high = exiting water temperature with fluid cooler fans at high speed, ^{o}C
The average fan power for the simulation time step is calculated for the twospeed fluid cooler as follows
Pfan,avg=ω(Pfan,high)+(1−ω)(Pfan,low)
Cooling Towers and Evaporative Fluid Coolers[LINK]
One, Two, and Variable Speed Cooling Towers and Evaporative Fluid Coolers[LINK]
Overview[LINK]
The input objects CoolingTower:SingleSpeed, CoolingTower:TwoSpeed, and CoolingTower:VariableSpeed:Merkel provide models for singlespeed, twospeed, and variablespeed cooling towers that are based on Merkel’s theory (Merkel 1925), which is also the basis for the tower model included in ASHRAE’s HVAC1 Toolkit for primary HVAC system energy calculations (ASHRAE 1999, Bourdouxhe et al. 1994). Cooling tower performance is modeled using effectivenessNTU relationships for counterflow heat exchangers. The model can be used to simulate the performance of both single speed, two speed, and variable speed mechanicaldraft cooling towers. The model will also account for tower performance in the “free convection” regime, when the tower fan is off but the water pump remains on. For partload operation, the model assumes a simple linear interpolation between two steadystate regimes without accounting for any cycling losses.
For single speed cooling towers, the capacity control can be fan cycling or fluid bypass. In fluid bypass mode, portion of the water goes through the tower media and gets cooled while the remaining water flow gets bypassed, two water flows then mix together trying to meet the tower exiting water setpoint temperature. In both the free convection cooling when fan is off and normal cooling when fan is on for the entire time step, if the tower exiting water temperature is lower than the setpoint, the tower operates in fluid bypass mode. The model determines the fluid bypass fraction by iterations until the mixed water meets the tower exiting water temperature setpoint. In the fluid bypass mode, except the free convection, the tower fan runs at full speed for the entire time step. The maximum amount of tower water that can be bypassed is bounded by the freezing point of the tower water  the tower exiting water temperature cannot be lower than the freezing setpoint.
Evaporative fluid coolers are modeled very similar to cooling towers. The main difference between the two is in the “Performance input method” input field. Cooling tower has two choices for this field namely “UFactorTimesAreaAndDesignWaterFlowRate” and “Nominal capacity”. The nominal capacity is specified for the standard conditions i.e. entering water at 35C (95F), leaving water at 29.44C (85F), entering air at 25.56C (78F) wetbulb temperature and 35C (95F) drybulb temperature. On the other hand evaporative fluid cooler has three choices for “Performance input method” which are “UFactorTimesAreaAndDesignWaterFlowRate”, “StandardDesignCapacity” and “UserSpecifiedDesignCapacity”. First method is same for both tower and fluid cooler. Standard design capacity is specified for the same conditions which are used to specify nominal capacity for tower as described above. If the capacity of fluid cooler for conditions other than the standard ones is known then UserSpecifiedDesignCapacity method should be used. In this case, the conditions for which the fluid cooler capacity is known i.e. entering water temperature, entering air temperature and entering air wet bulb temperature must be specified in the input. To calculate evaporation loss for fluid cooler, spray water flow rate which is different than the process fluid flow rate must be specified for all the performance input methods. This is not required for cooling tower because cooled fluid i.e. water is in direct contact with the air so the water loss is calculated by using cooled fluid flow rate only. Unlike cooling tower, evaporative fluid cooler model does not account for free convection.
Cooling tower model is described below which holds equally good for evaporative fluid cooler. The differences are mentioned whenever required.
Model Description[LINK]
Based on Merkel’s theory, the steadystate total heat transfer between the air and water entering the tower can be defined by the following equation:
d˙Qtotal=UdAcp(hs−ha)
where
h_{s} = enthalpy of saturated air at the wettedsurface temperature, J/kg
h_{a} = enthalpy of air in the free stream, J/kg
c_{p} = specific heat of moist air, J/kg^{o}C
U = cooling tower overall heat transfer coefficient, W/m^{2}^{o}C
A = heat transfer surface area, m^{2}
Equation is based on several assumptions:
air and water vapor behave as ideal gases
the effect of water evaporation is neglected
fan heat is neglected
the interfacial air film is assumed to be saturated
the Lewis number is equal to 1
In this model, it is also assumed that the moist air enthalpy is solely a function of the wetbulb temperature and that the moist air can be treated as an equivalent ideal gas with its mean specific heat defined by the following equation:
¯cpe=ΔhΔTwb
where
Δh = enthalpy difference between the air entering and leaving the tower, J/kg
ΔT_{wb} = wetbulb temperature difference between the air entering and leaving the tower, ^{o}C
Since the liquid side conductance is much greater than the gas side conductance, the wettedsurface temperature is assumed to be equal to the water temperature. Based on this assumption and equations and , the expression for total heat transfer becomes:
d˙Qtotal=UedA(Tw−Twb)
where
Ue=U¯cpecp
T_{wb} = wetbulb temperature of the air, ^{o}C
T_{w} = temperature of the water, ^{o}C
An energy balance on the water and air sides of the air/water interface yields the following equations:
d˙Qtotal=˙mwcpwdTw
d˙Qtotal=˙ma¯cpedTwb
where
˙mw = mass flow rate of water, kg/s
˙ma = mass flow rate of air, kg/s
Assuming that the heat capacity rate ( ˙mcp ) for the cooling tower water is less than that for the air, the effectiveness of the cooling tower can be defined by analogy to the effectiveness of a simple heat exchanger:
ε=Twin−TwoutTwin−Twbin
where
ε = heat exchanger effectiveness
T_{win} = inlet water temperature, ^{o}C
T_{wout} = outlet water temperature, ^{o}C
T_{wbin} = wetbulb temperature of the inlet air, ^{o}C
Combining equations , , and and integrating over the entire heat transfer surface area, and combining the result with equation provides the following expression for cooling tower effectiveness:
ε=1−exp{−NTU[1−(˙Cw˙Ca)]}1−(˙Cw˙Ca)exp{−NTU[1−(˙Cw˙Ca)]}
where
˙Cw=˙mwcpw and ˙Ca=˙ma¯cpe
NTU=NumberofTransferUnits=UAe˙Cw
This equation is identical to the expression for effectiveness of an indirect contact (i.e., fluids separated by a solid wall) counterflow heat exchanger (Incropera and DeWitt 1981). Therefore, the cooling tower can be modeled, in the steadystate regime, by an equivalent counterflow heat exchanger as shown in the following figure.
Cooling Tower Schematic
The first fluid is water and the second fluid is an equivalent fluid entering the heat exchanger at temperature T_{wbin} and specific heat ¯cpe . The heat exchanger is characterized by a single parameter, its overall heat transfer coefficientarea product UA_{e}. The actual cooling tower heat transfer coefficientarea product is related to UA_{e} by the following expression:
UA=UAecp¯cpe
This heat transfer coefficientarea product is assumed to be a function of the air mass flow rate only and can be estimated from laboratory test results or manufacturers’ catalog data.
The model for the variable speed Merkel tower also includes Scheier’s modifications. Scheier has extended the Merkel model to also include terms that adjust UA with three factors that model how UA values change when the tower is operating away from its rated conditions. The first factor, fUA,wetbulb , adjusts UA for the current outdoor wetbulb temperature. The user enters a performance curve or lookup table that is a function of one independent variable. The independent variable is the difference between the design wetbulb temperature and the current wetbulb temperature, in degrees Celsius.
fUA,wetbulb=func((Tdes,wetbulb−Tcurrent,wetbulb))
The second factor, fUA,airflow , adjusts UA for the current air flow rate. The user enters a performance curve or lookup table that is a function of one independent variable. The independent variable is the ratio of the current air flow rate to the design air flow rate at full speed.
fUA,airflow=func(AirFlowRatio)
The third factor, fUA,waterflow , adjusts UA for the current water flow rate. The user enters a performance curve or lookup table that is a function of one independent variable. The independent variable is the ratio of the current water flow rate to the design water flow rate.
fUA,waterflow=func(WaterFlowRatio)
Then the UA value at any given time is calculated using
UAe=UAdesign∙fUA,wetbulb∙fUA,airflow∙fUA,waterflow
Method for Calculating SteadyState Exiting Water Temperature[LINK]
The objective of the cooling tower model is to predict the exiting water temperature and the fan power required to meet the exiting water setpoint temperature. Since only the inlet air and inlet water temperatures are known at any simulation time step, an iterative procedure is required to determine the exiting fluid temperatures using the equations defined in the previous section. In the case of the EnergyPlus model, the iterations are performed to determine the exiting wetbulb temperature of the air. The exiting water temperature is then calculated based on an energy balance that assumes that the energy absorbed by the air is equivalent to the energy removed from the water. The procedure for calculating the steadystate, exiting air wetbulb temperature is outlined below.
As explained previously, it is assumed that the moist air enthalpy can be defined by the wetbulb temperature alone. Therefore, the first step in the procedure is to calculate the enthalpy of moist air entering the cooling tower based on the ambient wetbulb temperature from the weather file. Since an iterative solution is required, a first guess of the outlet air wetbulb temperature is then made and the enthalpy of this estimated outlet air wetbulb temperature is calculated. Based on these inlet and outlet air conditions, the mean specific heat of the air is calculated based on equation , repeated here:
¯cpe=ΔhΔTwb
With the overall heat transfer coefficientarea product for the cooling tower entered by the user, the effective heat transfer coefficientarea product is calculated by rearranging equation :
UAe=UA¯cpecp
With ¯cpe and UAe known, the effectiveness of the heat exchanger is then calculated:
ε=1−exp{−NTU[1−(˙Cmin˙Cmax)]}1−(˙Cmin˙Cmax)exp{−NTU[1−(˙Cmin˙Cmax)]}
where
˙Cmin=Minimum(˙Cw,˙Ca) and ˙Cmax=Maximum(˙Cw,˙Ca)
˙Cw=˙mwcpw and ˙Ca=˙ma¯cpe
NTU=NumberofTransferUnits=UAe˙Cmin
The heat transfer rate is then calculated as follows:
˙Qtotal=ε˙Cmin(Twin−Twbin)
The outlet air wetbulb temperature is then recalculated:
Twbout=Twbin+˙Qtotal˙Ca
The iterative process of calculating Twbout continues until convergence is reached.
Finally, the outlet water temperature is calculated as follows:
Twout=Twin+˙Qtotal˙Cw
Calculating the Actual Exiting Water Temperature and Fan Power[LINK]
The previous section describes the methodology used for calculating the steadystate temperature of the water leaving the cooling tower. This methodology is used to calculate the exiting water temperature in the free convection regime (water pump on, tower fan off) and with the tower fan operating (including low and high fan speed for the twospeed tower). The exiting water temperature calculations use the fluid flow rates (water and air) and the UAvalues entered by the user for each regime.
The cooling tower model seeks to maintain the temperature of the water exiting the cooling tower at (or below) a setpoint. The model obtains the target temperature setpoint from the setpoints placed on either the tower outlet node or the loop’s overall setpoint node (typically set to the supply side outlet node). The model checks to see if the outlet node has a setpoint placed on it and uses that if it does. If the outlet node does not have a temperature setpoint then the model uses the looplevel outlet node specified in the input field called Loop Temperature Setpoint Node Name in the PlantLoop or CondenserLoop object. The model first checks to determine the impact of “free convection”, if specified by the user, on the tower exiting water temperature. If free convection is not specified by the user, then the exiting water temperature is initially set equal to the entering tower water temperature. If the user specifies “free convection” and the steadystate exiting water temperature based on “free convection” is at or below the setpoint, then the tower fan is not turned on.
If the exiting water temperature remains above the setpoint after “free convection” is modeled, then the tower fan is turned on to reduce the exiting water temperature to the setpoint. The model assumes that partload operation is represented by a simple linear interpolation between two steadystate regimes (e.g., tower fan on for the entire simulation time step and tower fan off for the entire simulation time step). Cyclic losses are not taken into account.
The fraction of time that the tower fan must operate is calculated based on the following equation:
ω=Tset−Twout,offTwout,on−Twout,off
where
Tset = exiting water setpoint temperature, ^{o}C
Twout,off = exiting water temperature with tower fan off, ^{o}C
Twout,on = exiting water temperature with tower fan on, ^{o}C
The average fan power for the simulation time step is calculated by multiplying ω by the steadystate fan power specified by the user.
The calculation method for the twospeed tower is similar to that for the singlespeed tower example described above. The model first checks to see if “free convection” is specified and if the resulting exiting water temperature is below the setpoint temperature. If not, then the model calculates the steadystate exiting water temperature with the tower fan at low speed. If the exiting water temperature at low fan speed is below the setpoint temperature, then the average fan power is calculated based on the result of equation and the steadystate, low speed fan power specified by the user. If lowspeed fan operation is unable to reduce the exiting water temperature below the setpoint, then the tower fan is increased to its high speed and the steadystate exiting water temperature is calculated. If this temperature is below the setpoint, then a modified version of equation is used to calculate runtime at high fan speed:
ω=Tset−Twout,lowTwout,high−Twout,low
where
Tset = exiting water setpoint temperature, ^{o}C
Twout,low = exiting water temperature with tower fan at low speed, ^{o}C
Twout,high = exiting water temperature with tower fan at high speed, ^{o}C
The average fan power for the simulation time step is calculated for the twospeed cooling tower as follows:
Pfan,avg=ω(Pfan,high)+(1−ω)(Pfan,low)
The calculation method for the variablespeed Merkel/Scheier model is different from the one and twospeed models. Whereas the other towers are setpoint driven, the variable speed Merkel/Scheier model is driven by the load dispatched from the plant operation scheme, such as PlantEquipmentOperation:CoolingLoad. The plant operation scheme provides the cooling tower model with a target load to meet and the tower is controlled to meet that load (as opposed to meeting a leaving setpoint). This allows complex configurations with more than one cooling tower to be properly dispatched using all the features available in the various plant operation schemes. The model first evaluates the load provided with the fans off using only free convection. If the free convection mode meets or exceeds the target load then the tower runs in free convection mode. The model then evaluates the load provided by running the fans at full speed. If the full fan speed load provided is less than or equal to the target load then the tower runs at full speed. The model then evaluates the load provided with the variable speed fan running at the minimum air speed ratio, which is a user input (default is 0.2). If the minimum fan speed load provided meets or exceeds the target load, then tower runs at the minimum fan speed. If the target load is between the load provided by minimum and maximum fans speeds, then the model solves for an airflow rate ratio that just meets the target load using Regula Falsi numerical method. The variable speed Merkel/Scheier model does not model bypass and may provided excess cooling at times when running in free convection or at the minimum fan speed. Fan power for the variable speed Merkel/Scheier model is then calculated from the airflow rate ratio using a design fan power and a performance curve or lookup table with one independent variable for the airflow ratio.
Pfan=Pfan,des∙func(AirFlowRatio)
Tower Basin Heater[LINK]
The tower basin heater operates in the same manner as the variable speed cooling tower basin heater. Refer to the variable speed cooling tower basin heater description in the following section.
Cooling Tower Makeup Water Usage[LINK]
The cooling tower makeup water usage is the same as the variable speed cooling tower makeup water usage. Refer to the variable speed cooling tower makeup water usage description in the following section.
References[LINK]
Rosaler, Robert C. 1995. Standard Handbook of Plant Engineering, 2^{nd} Ed. New York, NY: McGrawHill, pp. 63637.
Variable Speed Cooling Towers Empirical Models[LINK]
Overview[LINK]
The input object CoolingTower:VariableSpeed provides models for variable speed towers that are based on empirical curve fits of manufacturer’s performance data or field measurements. The user specifies tower performance at design conditions, and empirical curves are used to determine the approach temperature and fan power at offdesign conditions. The user defines tower performance by entering the inlet air wetbulb temperature, tower range, and tower approach temperature at the design conditions. The corresponding water flow rate, air flow rate, and fan power must also be specified. The model will account for tower performance in the “free convection” regime, when the tower fan is off but the water pump remains on and heat transfer still occurs (albeit at a low level). Basin heater operation and makeup water usage (due to evaporation, drift, and blowdown) are also modeled.
The cooling tower seeks to maintain the temperature of the water exiting the cooling tower at (or below) a setpoint. The setpoint temperature is defined by the setpoints placed on either the tower outlet node or the loop’s overall setpoint node (typically set to the supply side outlet node). The model checks to see if the outlet node has a setpoint placed on it and uses that if it does. If the outlet node does not have a temperature setpoint then the model uses the looplevel outlet node specified in the input field called Loop Temperature Setpoint Node Name in the PlantLoop or CondenserLoop object. The model simulates the outlet water temperature in four successive steps:
The model first determines the tower outlet water temperature with the tower fan operating at maximum speed. If the outlet water temperature is above the setpoint temperature, the fan runs at maximum speed.
If the outlet water temperature with maximum fan speed is below the setpoint temperature, then the model next determines the impact of “free convection” (water flowing through tower with fan off). If the exiting water temperature based on “free convection” is at or below the setpoint, then the tower fan is not turned on.
If the outlet water temperature remains above the setpoint after “free convection” is modeled, then the tower fan is turned on at the minimum fan speed (minimum air flow rate ratio) to reduce the leaving water temperature. If the outlet water temperature is below the setpoint at minimum fan speed, the tower fan is cycled on and off to maintain the outlet water setpoint temperature.
If the outlet water temperature remains above the setpoint after minimum fan speed is modeled, then the tower fan is turned on and the model determines the required air flow rate and corresponding fan speed to meet the desired setpoint temperature.
Model Description[LINK]
The variable speed tower model utilizes userdefined tower performance at design conditions along with empirical curves to determine tower heat rejection and fan power at offdesign conditions. Basin heater operation and makeup water usage are also modeled based on user inputs, tower entering air conditions, and tower operation. The following sections describe how each of these tower performance areas is modeled.
Tower Heat Rejection[LINK]
Heat rejection by the variable speed cooling tower is modeled based on the CoolTools correlation, YorkCalc correlation, or userdefined coefficients for either the CoolTools or YorkCalc correlations. These purelyempirical correlations model the tower approach temperature using a polynomial curve fit with a large number of terms and either three or four independent variables.
The CoolTools correlation has 35 terms with four independent variables:
where:
Approach = approach temperature (^{o}C) = outlet water temperature minus inlet air wetbulb temperature
FRair = air flow rate ratio (actual air flow rate divided by design air flow rate)
FRwater = water flow rate ratio (actual water flow rate divided by design water flow rate)
Tr = range temperature (^{o}C) = inlet water temperature minus outlet water temperature
Twb = inlet air wetbulb temperature (^{o}C)
Coeff(#) = correlation coefficients
If the user selects Tower Model Type = CoolToolsCrossFlow, then the 35 coefficients derived for the CoolTools simulation model (Benton et al. 2002) are used and these coefficients are already defined within EnergyPlus as shown in Table 63. If the user specifies Tower Model Type = CoolToolsUserDefined, then the user must enter a CoolingTowerPerformance:CoolTools object to define the 35 coefficients that will be used by the CoolTools approach temperature correlation.
: Approach Temperature Correlation Coefficients
Similarly, the YorkCalc correlation has 27 terms with three independent variables:
where:
Approach = approach temperature (^{o}C) = outlet water temperature minus inlet air wetbulb temperature
Tr = range temperature (^{o}C) = inlet water temperature minus outlet water temperature
Twb = inlet air wetbulb temperature (^{o}C)
LGratio = liquidtogas ratio = ratio of water flow rate ratio (FRwater) to air flow rate ratio (FRair)
Coeff(#) = correlation coefficients
If the user selects Tower Model Type = YorkCalc, then the 27 coefficients derived for the YorkCalc simulation model (York International Corp. 2002) are used and these coefficients are already defined within EnergyPlus as shown in Table 63. If the user specifies Tower Model Type = YorkCalcUserDefined, then the user must enter a CoolingTowerPerformance:YorkCalc object to define the 27 coefficients that will be used by the YorkCalc approach temperature correlation.
The approach temperature correlations for the CoolTools and YorkCalc simulation models are valid for a range of conditions defined in Table 64. If the user defines their own model coefficients (CoolingTowerPerformance:CoolTools or CoolingTowerPerformance:YorkCalc), then they must also define in that same object the range of conditions for which the model is valid. For all of these correlation variables, the program issues warnings if the actual values are beyond the minimum/maximum values specified for the correlation being used. For inlet air wetbulb temperature and water mass flow rate ratio, the values of these variables used in the calculation of approach temperature are limited to be within the valid minimum/maximum range. For approach, range, and liquidtogas ratio the warnings are issued if the values are beyond the specified minimum/maximum range but the actual values are still used. The warnings issued do not necessarily indicate a poor estimate of tower performance at the condition(s) which caused the warning, but are provided to identify conditions outside the defined correlation limits. Exceeding the defined limits by a small amount may not introduce significant errors, but large deviations may be problematic. It is for this reason that we recommend using a very broad range of cooling tower performance data (i.e., data covering the entire range expected during the simulation) when generating userdefined coefficients for the variable speed tower model.
The approach temperature correlation(s) used to simulate cooling tower heat rejection are based on water and air flow rate “ratios” and are not directly dependent on the size of the tower or the actual air and water flow rates through the tower. However, the model correlations are developed based on a reference condition. For Model Types “CoolToolsCrossFlow” and “YorkCalc”, the reference condition is a water flow rate of 0.000043 m^{3}/s per kW of heat rejected (2.4 gal/min per ton of heat rejected) with 25.6^{o}C (78^{o}F) enter air wetbulb temperature, 35^{o}C (95^{o}F) hot water inlet temperature, and 29.4^{o}C (85^{o}F) cold water outlet temperature. The reference condition may be different if the user defines tower model coefficients using CoolingTowerPerformance:CoolTools or CoolingTowerPerformance:YorkCalc.
Due to the inherent reference condition used to generate the tower performance curves, the water flow rate at the reference condition must be determined using the design performance information specified by the user and the tower model’s approach temperature correlation. This is done by using the model’s approach temperature correlation (described earlier in this section) to calculate the water flow rate ratio which yields the userdefined design approach temperature based on an air flow rate ratio of 1.0 (FR_{air} = 1.0), the design inlet air wetbulb temperature, and the design range temperature. The calculated approach temperature (using the model correlation) must satisfy the following two equations:
Twater,outlet,design=Twb,air,design+Tapproach,design{FRair=1.0} Twater,inlet,design=Twb,air,design+Tapproach,design+Trange,design{FRair=1.0}
where:
Twater,outlet,design = design outlet water temperature (^{o}C)
Twater,inlet,design = design inlet water temperature (^{o}C)
Twb,air,design = design inlet air wetbulb temperature (^{o}C)
Tapproach,design = design approach temperature (^{o}C)
Trange,design = design range temperature (^{o}C)
FRair = air flow rate ratio (actual air flow rate divided by design air flow rate)
The water flow rate ratio used in the approach temperature correlation which satisfies these two equations is the ratio of the design water flow rate (specified by the user) to the water flow rate at the reference condition. This ratio is used to calculate the reference water volumetric flow rate, which is then used throughout the simulation to determine the actual water flow rate ratio used in the approach temperature correlation for each simulation time step.
˙Vreference=˙VdesignFRwaterreference
where:
˙Vreference = water volumetric flow rate at the reference condition (m^{3}/s)
˙Vdesign = design water volumetric flow rate specified by the user (m^{3}/s)
FRwater,reference = design water flow rate divided by the reference water flow rate
The cooling tower seeks to maintain the temperature of the water exiting the cooling tower at (or below) a setpoint. The setpoint temperature is defined by the field “Condenser Loop Temperature Setpoint schedule or reference” for the CondenserLoop object. The model simulates the outlet water temperature in four successive steps:
The model first determines the tower outlet water temperature with the tower fan operating at maximum speed. If the outlet water temperature is above the setpoint temperature, the fan runs at maximum speed.
If the outlet water temperature with maximum fan speed is below the setpoint temperature, then the model next determines the impact of “free convection” (water flowing through tower with fan off). If the exiting water temperature based on “free convection” is at or below the setpoint, then the tower fan is not turned on.
If the outlet water temperature remains above the setpoint after “free convection” is modeled, then the tower fan is turned on at the minimum fan speed (minimum air flow rate ratio) to reduce the leaving water temperature. If the outlet water temperature is below the setpoint at minimum fan speed, the tower fan is cycled on and off to maintain the outlet water setpoint temperature.
If the outlet water temperature remains above the setpoint after minimum fan speed is modeled, then the tower fan is turned on and the model determines the required air flow rate and corresponding fan speed to meet the desired setpoint temperature.
For each simulation time step, the model first calculates the outlet water temperature with the tower fan operating at maximum speed (FRair = 1.0). The calculated approach temperature (using the correlations described above), inlet air wetbulb temperature (weather data), and range temperature are used to determine the tower outlet water temperature as follows:
Twater,outlet,fanMAX=Twb,air+Tapproach{FRair=1.0}=Twater,inlet+Trange
where:
Twater,outlet,fanMAX = tower outlet water temperature at maximum fan speed (^{o}C)
Twb,air = tower inlet air wetbulb temperature (^{o}C)
Tapproach = approach temperature at current operating conditions (^{o}C)
Trange = range temperature at current operating conditions (^{o}C)
Note that the approach temperature correlation as described previously is a function of range temperature, so the equations above must be solved iteratively to converge on a solution. If the resulting outlet water temperature is above the desired setpoint temperature, then the fan runs at maximum speed and does not cycle on/off (fan partload ratio = FanPLR = 1.0 and FR_{air} = 1.0).
If the outlet water temperature with maximum fan speed is below the setpoint temperature, then the model next determines the impact of “free convection” (water flowing through tower with fan off). In the free convection regime, the outlet water temperature is calculated using a fraction of the water temperature difference through the tower when the fan is at its maximum speed. This fraction is defined by the user (Fraction of Tower Capacity in Free Convection Regime).
Twater,outlet,fanOFF=Twater,inlet−[Fracfreeconv(Twater,inlet−Twater,outlet,fanMAX)]
where:
Twater,outlet,fanOFF = tower outlet water temperature in free convection regime (^{o}C)
Twater,inlet = tower inlet water temperature (^{o}C)
Fracfreeconv = fraction of tower capacity in free convection regime (user specified)
If the outlet water temperature in the free convection regime is below the setpoint temperature, the tower fan is not turned on and the fan partload ratio is set equal to 0. In addition, the air flow rate ratio through the tower is assumed to be equal to the fraction of tower capacity in the free convection regime.
FanPLR=FanPLRfreeconv=0.0
FRair=FRair,freeconv=Fracfreeconv
where:
FanPLR = fan partload ratio
FanPLRfreeconv = fan partload ratio in free convection regime
FRair,freeconv = air flow rate ratio in free convection regime
If the outlet water temperature in the free convection regime is above the setpoint temperature, then the fan is turned on at the minimum fan speed (minimum air flow rate ratio, FR_{air,min}, entered by the user) and the outlet water temperature is calculated as the inlet air wetbulb temperature plus the calculated approach temperature:
FRair=FRair,min
Twater,outlet,fanMIN=Twb,air+Tapproach{FRair=FRair,min}
where:
Twater,outlet,FanMIN = outlet water temperature at minimum fan speed (^{o}C)
FRair,min = air flow rate ratio at the minimum fan speed
If the outlet water temperature at minimum fan speed is below the setpoint temperature, the cooling tower fan cycles on and off at the minimum air flow rate ratio in order to meet the setpoint temperature.
FanPLR=(Twater,outlet,fanOFF−TsetpointTwater,outlet,fanOFF−Twater,outlet,fanMIN)
where:
Tsetpoint = outlet water setpoint temperature (^{o}C)
If the outlet water temperature at minimum fan speed is above the outlet water temperature setpoint, then the cooling tower fan speed (FR_{air}) is increased until the calculated approach temperature produces the required outlet water temperature to meet the setpoint.
Tapproach=f(FRair)
Twater,outlet=Tsetpoint=Twb,air+Tapproach
FanPLR=1.0 (i.e., fan does not cycle on/off)
Fan Power[LINK]
When the cooling tower fan is operating, fan electric power is calculated based on the air flow rate ratio required to meet the above conditions. If the user has entered a fan power curve object (cubic curve), the output of that curve is multiplied by the design fan power. Otherwise, tower fan power is assumed to be directly proportional to the cube of the air flow rate ratio. In either case, the fan partload ratio is applied to account for times when the tower fan cycles on/off to meet the setpoint temperature. Fan energy consumption is calculated each simulation time step.
IF(FanPowerCurveObjectisDefined)THENPfan=(CurveOutput)(FanPLR)Pfan,designELSEPfan=(Pfan,design)(FanPLR)(FRair)3ENDIF
Efan=(Pfan)(TimeStepSys)3600.
where:
FanPowerCurveObject = name of fan power ratio as a function of air flow rate ratio curve
Pfan = tower fan electric power (W)
Efan = tower fan electric consumption (J)
CurveOutput = output of FanPowerCurveObject evaluated at the operating air flow rate ratio (FR_{air})
Pfan,design = design fan power at design (maximum) air flow through the tower (W)
TimeStepSys = HVAC system simulation time step (hr)
Tower Basin Heater[LINK]
Calculations are also made to estimate the electric power input to the tower basin heater. A schedule may be used to disable the basin heater during regular maintenance periods or other time periods (e.g., during summer). If a schedule is not provided, the basin heater is assumed to be available the entire simulation time period. The basin heater operates when it is scheduled on, the outdoor air drybulb temperature is below the basin heater setpoint temperature, and the cooling tower is not active (i.e., water is not flowing through the tower). The user is required to enter a basin heater capacity (watts per degree Kelvin) and a heater setpoint temperature (^{o}C) if they want to model basin heater electric power.
Pheater,basin=0.0IF(WaterNotFlowingThroughTower)THENIF(Scheduleheater,basinisDefined)THENIF(CAPheater,basin>0.AND.Scheduleheater,basin=ON)THENPheater,basin=MAX(0.0,CAPheater,basin(Tsetpoint,basin−Tdb,outdoor))ENDIFELSEIF(CAPheater,basin>0)THENPheater,basin=MAX(0.0,CAPheater,basin(Tsetpoint,basin−Tdb,outdoor))ENDIFENDIFENDIF
Eheater,basin=(Pheater,basin)(TimeStepSys)3600.
where:
Pheater,basin = tower basin heater electric power (W)
Eheater,basin = tower basin heater electric consumption (J)
T_{setpoint,basin} = basin heater setpoint temperature (^{o}C)
T_{db,outdoor} = outdoor air drybulb temperature (^{o}C)
CAP_{heater,basin} = basin heater capacity (W/K)
Schedule_{heater,basin} = basin heater schedule (schedule value > 0 means ON)
References[LINK]
ASHRAE 1999. HVAC1 Toolkit: A Toolkit for Primary HVAC System Energy Calculations. Atlanta: American Society of Heating, Refrigerating and AirConditioning Engineers, Inc.
Benton, D.J., Bowman, C.F., Hydeman, M., Miller, P. 2002. An Improved Cooling Tower Algorithm for the CoolTools^{TM} Simulation Model. ASHRAE Transactions, Vol. 108, Part 1, pp.760768.
Bourdouxhe, J.P., M. Grodent, J. Lebrun and C. Silva. 1994. Cooling tower model developed in a toolkit for primary HVAC system energy calculation: part 1. Proceedings of the fourth international conference on system simulation in buildings, Liege (Belgium), December 57, 1994.
Incropera, F.P. and D.P. DeWitt. 1981. Fundamentals of Heat Transfer. New York: John Wiley & Sons.
Merkel, F. 1925. Verduftungskuhlung. VDI Forschungarbeiten, No 275, Berlin.
Rosaler, Robert C. 1995. Standard Handbook of Plant Engineering, 2^{nd} Ed. New York, NY: McGrawHill, pp. 63637.
Scheier, L. 2013. Personal communication.
York International Corporation, 2002. “YORKcalc^{TM} Software, ChillerPlant EnergyEstimating Program”, Form 160.00SG2 (0502).
Cooling Towers with Multiple Cells[LINK]
Many towers are constructed to be capable of being grouped together to achieve the desired capacity. Thus, many cooling towers are assemblies of two or more individual cooling towers or “cells.” The number of cells they have, e.g., an eightcell tower, often refers to such towers.
For the operation of multicell towers, the first step is to determine the number of cells n, which will be operating during the timestep using the calculation logic from DOE2.1E.
The maximum and minimum flow rates per cell are determined according to the input fractions (Minimum Water Flow Rate Fraction: kmin and Maximum Water Flow Rate Fraction: kmax) as follows:
⎧⎪⎨⎪⎩˙mw,min=˙mw,desntot×kmin˙mw,max=˙mw,desntot×kmax
where ˙mw,des is the design water flow rate through the entire cooling tower.
Then, we determine the minimum and maximum number of cells that can operate with this water flow rate:
⎧⎪ ⎪⎨⎪ ⎪⎩nmin=MIN(˙mw,tot˙mw,max,ntot)nmax=MIN(˙mw,tot˙mw,min,ntot)
where n_{tot} is the total number of cells of the tower, and mw,tot is the water flow rate to the tower.
The number of cells operating n is set accordingly:
If the Cell Control method is MinimalCell,
n=nmin
If the Cell Control method is MaximalCell,
n=nmax
Finally, the water mass flow rate per cell (˙mw) will be:
˙mw=˙mw,totn
Then we simulate the performance of one cell with this flow rate per cell (calling the SimSimpleTower subroutine for single and two speed cooling tower objects). As we assume that each cell is identical, the UA of one cell is calculated dividing the UA of the whole tower (obtained from the input or from the auto sizing calculations). The air flow rate per cell is also equal to the one of the whole tower divided by the number of cells operating:
UAcell=UAcellntot
˙mair,cell=˙mair,totn
At the end, the total fan power of the tower operating with a certain number of cells is given by:
FanPower=FanPowertot×nntot
If the cells operating do not meet the loads, we increase the number of cells if spare cells are available and the water flow through each cell is within the user specified minimum and maximum water flow rate fractions range. This is an iteration process.
Cooling Tower Makeup Water Usage[LINK]
Makeup water use for all types of cooling towers is made up of three components: evaporation, drift, and blowdown. The first is the amount of water evaporated to reduce the water’s temperature as it passes through the cooling tower. There are two methods that evaporation makeup water can be modeled in EnergyPlus. The first method assumes that the tower outlet air conditions are saturated (which may not always be the case for certain operating conditions). For this “Saturated Exit” mode, the enthalpy of the tower’s outlet air is calculated as the inlet air enthalpy plus the water side heat transfer divided by the air mass flow rate through the tower.
˙Qwater=˙mwaterCp(Twater,inlet−Twater,outlet)
hsatair,outlet=hair,inlet+(˙Qwater˙mair)
where:
˙Qwater = waterside heat transfer (W)
˙mwater = mass flow rate of water through the tower (kg/s)
Cp = specific heat of water (W/kgK)
hsatair,outlet = saturated outlet air enthalpy (J/kg)
hair,inlet = inlet air enthalpy (J/kg)
˙mair = mass flow rate of air through the tower (kg/s)
The saturation temperature and humidity ratio are then calculated for the tower’s outlet air.
Tsatair,outlet=PsyTsatFnHPb(hsatair,outlet,OutBaroPress)
ωsatair,outlet=PsyWFnTdbH(Tsatair,outlet,hsatair,outlet)
where:
Tsatair,outlet = saturated outlet air temperature (^{o}C)
PsyTsatFnHPb = EnergyPlus psychrometric function, returns saturation temperature given enthalpy and barometric pressure
OutBaroPress = outdoor barometric pressure (Pa)
ωsatair,outlet = saturated outlet air humidity ratio (kg_{water}/kg~dry air~)
PsyWFnTdbH = EnergyPlus psychrometric function, returns humidity ratio given drybulb temperature and enthalpy
The makeup water quantity required to replenish the water lost due to evaporation is then calculated as the product of the air mass flow rate and the difference between the entering and leaving air humidity ratio divided by the density of water.
˙Vevaporation,makeup=˙mair(ωsatair,outlet−ωair,inlet)ρwater
where:
˙Vevaporation,makeup = makeup water usage due to evaporation (m^{3}/s)
˙mair = mass flow rate of air through tower (kg/s)
ωair,inlet = humidity ratio of tower inlet air (kg_{water}/kg~dry air~)
ρwater = density of water evaluated at the tower inlet air temperature (kg/m^{3})
The second method available for calculating water makeup for evaporation is for the user to provide a value for a loss factor. The evaporation loss is then calculated as a fraction of the circulating condenser water flow and varies with the temperature change in the condenser water. The value provided by the user is in units of percentperdegree Kelvin. The evaporation rate will equal this value times each degree Kelvin of temperature drop in the condenser water. Typical values are from 0.15 to 0.27 [percent/K]. The default is 0.2. The rate of water makeup for evaporation is then calculated by multiplying this factor times the condenser water flow rate and the temperature decrease in the condenser water flow rate. For evaporative fluid coolers, a numerical value of loss factor can be entered in the same manner as for cooling towers. If this field is blank, an empirical correlation will be used to calculate the value based on current outdoor dry bulb temperature and relative humidity. The following correlation from Qureshi and Zubair (2007) is used to calculate the loss factor:
LossFactor=(113−8.417ϕ+1.6147tdb)∗10−5
where:
ϕ = relative humidity of inlet air
tdb = Drybulb temperature of inlet air
Additional makeup water usage is modeled as a percentage of design water flow rate through the tower to account for drift, and as a scheduled flow rate to model blowdown. Drift is water loss due to the entrainment of small water droplets in the air stream passing through the tower. Drift is defined by the model user as a percentage of the tower’s design water flow rate, and is assumed to vary with tower air flow rate ratio as follows:
˙Vdrift,makeup=˙mwater,design(PercentWatertoDrift100.0)FRair
where:
˙Vdrift,makeup = makeup water usage due to drift (m^{3}/s)
˙mwater,design = design (volumetric) water flow rate (m^{3}/s)
PercentWatertoDrift = percent of design water flow rate lost to drift at the tower design air flow rate
FRair = ratio of actual air flow rate to tower design air flow rate
Blowdown is water flushed from the basin on a periodic basis to purge the concentration of mineral scale or other contaminants. There are two ways that blowdown is calculated in EnergyPlus. Blowdown water rates can be scheduled so that we have:.
IF(Scheduleblowdownisdefined)THEN˙Vblowdown,makeup=ScheduleValueblowdownELSE˙Vblowdown,makeup=0.0ENDIF
where:
˙Vblowdown,makeup = makeup water usage due to blowdown (m^{3}/s)
ScheduleValue_{blowdown} = blowdown schedule value for the time step being simulated (m^{3}/s)
The second (and default) way that blowdown can be calculated is to assume that blowdown water is continually introduced at a rate that will provide a constant concentration ratio. As water evaporates it leaves behind minerals and the like causing the concentration of water impurities to be higher in the tower than in the makeup water. Acceptable concentration ratios are in the range of 3 to 5 depending on the purity of the make up water. Water lost as drift does not evaporate and decrease the water needed for blowdown. Using the “Concentration Ratio” method, the rate of blowdown can be calculated using:
˙Vblowdown,makeup=˙Vevaporation,makeupRconcentration−1−˙Vdrift,makeup
where,
Rconcentration is the concentration ratio or the ratio of solids in the blowdown water to solids in the makeup water.
The tower makeup water consumption (m^{3}) for each simulation time step is calculated as the sum of the individual components of makeup water usage multiplied by the simulation time step in hours and the conversion for hours to seconds (3600 sec/hr). Makeup water usage is only calculated when the cooling tower is active and water is flowing through the cooling tower.
Vwater,makeup=(˙Vevaporation,makeup+˙Vdrift,makeup+˙Vblowdown,makeup)TimeStepSys3600.
where:
Vwater,makeup = tower makeup water consumption (m^{3})
References[LINK]
Robert C. Rosaler. 1995. Standard Handbook of Plant Engineering, 2^{nd} Ed. McGrawHill, New York, NY, pp 63637.
Quereshi, B.A. and S.M.Zubair. 2007. Prediction of evaporation losses in evaporative fluid coolers, Applied Thermal Engineering 27 pp. 520527
One and Two Speed Fluid Coolers[LINK]
Overview[LINK]
The input objects FluidCooler:SingleSpeed and FluidCooler:TwoSpeed provide models for dry fluid coolers. Fluid cooler’s performance is modeled using effectivenessNTU relationships for cross flow heat exchanger with both streams unmixed. The model can be used to simulate the performance of both single speed and two speed mechanicaldraft fluid coolers. For partload operation, the model assumes a simple linear interpolation between two steadystate regimes without accounting for any cycling losses.
Model Description[LINK]
The expression for fluid cooler effectiveness is as follows:
ε=1−e[e(−NTUCrη)−1Crη]
Where
ε = heat exchanger effectiveness
Cr=˙Cmin˙Cmax
˙Ca=˙macpa Ca=macpaCa=macpa and ˙Cw=˙mwcpw
˙Cmax=Max(˙Ca,˙Cw) Cmax=Max (Ca,Cw)Cmax=Max(Ca,Cw); ˙Cmin=Min(˙Ca,˙Cw)
NTU=UA˙Cmin
η=NTU(−0.22) η=NTU0.22η=NTU0.22
The first fluid is water and the second fluid is air entering the heat exchanger at temperature Tdb,in TdbinTdbin and specific heat cpa cpacpa. The heat exchanger is characterized by a single parameter, its overall heat transfer coefficientarea product UA.
When the user selects the nominal capacity method, the UA is calculated as follows:
The model inputs (other than the UA) and the fluid cooler load that it must meet are specified at design conditions. Then the fluid cooler model converges to a UA value, using the regulafalsi method that will enable it to meet the design fluid cooler load given at the specified inputs.
Method for Calculating SteadyState Exiting Water Temperature[LINK]
The objective of the fluid cooler model is to predict the exiting water temperature and the fan power required to meet the exiting water setpoint temperature. The exiting water temperature is calculated based on an energy balance that assumes that the energy absorbed by the air is equivalent to the energy removed from the water. The procedure for calculating the steadystate, exiting air drybulb temperature is outlined below.
With the overall heat transfer coefficientarea product for the fluid cooler calculated by the nominal capacity information entered by the user, the effectiveness of the heat exchanger is then calculated as:
ε=1−e[e(−NTUCrη)−1Crη]
The heat transfer rate is then calculated as follows:
˙Q=ε˙Cmin(Tw,in−Tdb,in)
˙Q=ε˙Cmin(Tw,in−Tdbin)
Then the outlet air drybulb and outlet water temperature are calculated:
Tdb,out=Tdb,in+.Q.Ca
Tw,out=Tw,in−.Q.Cw
Tw,in = inlet water temperature, ^{o}C
Tw,out = outlet water temperature, ^{o}C
Tdb,in = drybulb temperature of the inlet air, ^{o}C
Tdb,out = drybulb temperature of the outlet air, ^{o}C
Calculating the Actual Exiting Water Temperature and Fan Power[LINK]
The previous section describes the methodology used for calculating the steadystate temperature of the water leaving the fluid cooler. This methodology is used to calculate the exiting water temperature with the fluid cooler fans operating (including low and high fan speed for the twospeed fluid cooler). The exiting water temperature calculations use the fluid flow rates (water and air) and the Nominal capacity information entered by the user for each regime.
The fluid cooler model seeks to maintain the temperature of the water exiting the fluid cooler at (or below) a setpoint. The setpoint schedule is defined by the field “Loop Temperature Setpoint Node or reference” for the CondenserLoop object.
The fluid cooler fans are turned on to reduce the exiting water temperature to the setpoint. The model assumes that partload operation is represented by a simple linear interpolation between two steadystate regimes (e.g., Fluid cooler fans on for the entire simulation time step and fluid cooler fans off for the entire simulation time step). Cyclic losses are not taken into account. If the outlet water temperature is less than the setpoint then the fraction of time for which the fluid cooler must operate to meet the setpoint is calculated by using the following equation:
ω=Tset−Twout,offTwout,on−Twout,off
Where
Tset = exiting water setpoint temperature, ^{o}C
Twout,off = exiting water temperature with all fluid cooler fans off, ^{o}C
Twout,on = exiting water temperature with all fluid cooler fans on, ^{o}C
The average fan power for the simulation time step is calculated by multiplying ω by the steadystate fan power specified by the user.
The calculation method for the twospeed fluid cooler is similar to that for the singlespeed fluid cooler example described above. The model first calculates the steadystate exiting water temperature with the fluid cooler fans at low speed. If the exiting water temperature at low fan speed is below the setpoint temperature, then the average fan power is calculated based on the result of previous equation and the steadystate, low speed fan power specified by the user. If lowspeed fan operation is unable to reduce the exiting water temperature below the setpoint, then the fluid cooler fans’ speed is increased to high speed and the steadystate exiting water temperature is calculated. If this temperature is below the setpoint, then a modified version of previous equation is used to calculate runtime at high fan speed:
ω=Tset−Twout,lowTwout,high−Twout,low
where
Tset = exiting water setpoint temperature, ^{o}C
Twout,low = exiting water temperature with fluid cooler fans at low speed, ^{o}C
Twout,high = exiting water temperature with fluid cooler fans at high speed, ^{o}C
The average fan power for the simulation time step is calculated for the twospeed fluid cooler as follows
Pfan,avg=ω(Pfan,high)+(1−ω)(Pfan,low)
Documentation content copyright © 19962014 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.