Model for a steady-state or dynamic cooling tower with variable speed fan using the York calculation for the approach temperature at off-design conditions.
To compute the thermal performance, this model takes as parameters the approach temperature, the range temperature and the inlet air wet bulb temperature at the design condition. Since the design mass flow rate (of the chiller condenser loop) is also a parameter, these parameters define the rejected heat.
For off-design conditions, the model uses the actual range
temperature and a polynomial to compute the approach temperature
for free convection and for forced convection, i.e., with the fan
operating. The polynomial is valid for a York cooling tower. If the
fan input signal y is below the minimum fan revolution
yMin, then the cooling tower operates in free
convection mode, otherwise it operates in the forced convection
mode. For numerical reasons, this transition occurs in the range of
y ∈ [0.9*yMin, yMin].
The fan power consumption at the design condition can be specified as follows:
fraPFan_nominal can be used to
specify at the nominal conditions the fan power divided by the
water flow rate. The default value is 275 Watts for a water
flow rate of 0.15 kg/s.PFan_nominal can be set to the fan
power at nominal conditions. If a user does not set this parameter,
then the fan power will be PFan_nominal = fraPFan_nominal *
m_flow_nominal, where m_flow_nominal is the
nominal water flow rate.In the forced convection mode, the actual fan power is computed
as PFan=fanRelPow(y) * PFan_nominal, where the default
value for the fan relative power consumption at part load is
fanRelPow(y)=y3. In the free convection
mode, the fan power consumption is zero. For numerical reasons, the
transition of fan power from the part load mode to zero power
consumption in the free convection mode occurs in the range y
∈ [0.9*yMin, yMin].
To change the fan relative power consumption at part load in the
forced convection mode, points of fan controls signal and
associated relative power consumption can be specified. In between
these points, the values are interpolated using cubic splines.
This model is similar to the model Cooling Tower:Variable
Speed that is implemented in the EnergyPlus building energy
simulation program version 6.0. The main differences are
y.This model requires a medium that has the same computation of the enthalpy as Buildings.Media.Water, which computes
h = cp (T-T0),
where h is the enthalpy, cp = 4184
J/(kg K) is the specific heat capacity, T is the temperature
in Kelvin and T0 = 273.15 Kelvin. If this is not
the case, the simulation will stop with an error message. The
reason for this limitation is that as of January 2015, OpenModelica
failed to translate the model if Medium.temperature()
is used instead of Water.temperature().
EnergyPlus 2.0.0 Engineering Reference, April 9, 2007.
FRWat0, as this variable can
take on values that are bigger than 1.Medium.temperature() with
Water.temperature() in order for the model to work
with OpenModelica. Added an assert that stops the
simulation if an incompatible medium is used.each keyword for
fanRelPowDer. Added regularization in computation of
TAppCor. Removed intermediate states with
temperatures.Evaluate=true.allowReverseFlow=false. Avoided the use of the
conditionally enabled variables sta_a and
sta_b as this was not proper use of the Modelica
syntax.Q_flow and
mXi_flow.mWat_flow_nominal
since it is equal to m_flow_nominal, which is the
water flow rate from the chiller condenser loop.