Air System Humidifiers[LINK]
Air system humidifiers are components that add moisture to the supply air stream. They fall into 2 broad categories: spray type humidifiers which act like direct evaporative coolers, cooling the supply air as well as humidifying it; and dry steam humidifiers, which humidify the supply air stream while causing almost no change to the supply air stream temperature. The EnergyPlus electric steam humidifier uses electrical energy to convert ordinary tap water to steam which it then injects into the supply air stream by means of a blower fan. The actual unit might be an electrode-type humidifier or a resistance-type humidifier.
Electric Steam Humidifier[LINK]
The electric steam humidifier model (object name: Humidifier:Steam:Electric) is based on moisture and enthalpy balance equations plus standard psychrometric relationships. The approach is similar to that taken in the ASHRAE HVAC 2 Toolkit, page 4-112 (ASHRAE 1993). EnergyPlus contains its own module of psychrometric routines; the psychrometric theory and relations are given in the 2001 edition of ASHRAE Fundamentals, Chapter 6 (ASHRAE 2001). The model contains both an ideal controller and the component. The control model assumes that there is a minimum humidity setpoint on the component air outlet node. This setpoint is established by a setpoint manager described elsewhere.
The component model is a forward model: its inputs are its inlet conditions; its outputs are its outlet conditions and its energy consumption. The inputs are the temperature, humidity ratio, and mass flow rate of the inlet air stream, which are known; and the water addition rate (kg/s) which is determined by the controller.
The controller first decides whether the humidifier is on or off. For the humidifier to be on:
- the humidifier schedule value must be nonzero;
- the inlet air mass flow must be greater than zero;
- the inlet air humidity ratio must be less than the minimum humidity ratio setpoint.
If the humidifier is off, the water addition rate is set to zero. If the humidifier is on, the water addition rate needed to meet the humidity setpoint is calculated.
˙ma⋅win+˙mw,add,needed=˙ma⋅wset
where
˙ma = the air mass flow rate [kg/s]
win = the inlet air humidity ratio [kg/kg]
˙mw,add,needed = water addition rate needed to meet the setpoint [kg/s]
wset = the humidity ratio setpoint [kg/kg]
Equation is the moisture balance equation for the component. It is solved for ˙mw,add,needed (the other variables are known) which is passed to the humidifier component model as its desired inlet water addition rate.
The inputs to the component model are the air inlet conditions and mass flow rate and the water addition rate set by the controller. The outputs are the air outlet conditions. First the desired water addition rate is checked against component capacity.
˙mw,add,needed,max=Min(˙mw,add,Capnom)
where
Capnom = the humidifier nominal capacity [kg/s], a user input.
If ˙mw,add,needed,max is zero, the outlet conditions are set to the inlet conditions and the water addition rate is set to zero. If the humidifier is scheduled on the component power consumption is set to the standby power consumption: Whum=Wstby . Otherwise Whum = 0.
If ˙mw,add,needed,max > 0, then the moisture and enthalpy balance equations
˙ma⋅win+˙mw=˙ma⋅wout
˙ma⋅hin+˙mw⋅hw=˙ma⋅hout
with ˙mw set equal to ˙mw,add,needed,max are solved for wout and hout . Here
˙ma = the air mass flow rate [kg/s]
win = the inlet air humidity ratio [kg/kg]
˙mw = the inlet water addition rate [kg/s]
wout = the outlet air humidity ratio [kg/kg]
hin = the inlet air specific enthalpy [J/kg]
hw = the steam specific enthalpy = 2676125. [J/kg] at 100 oC
hout = the outlet air specific enthalpy [J/kg]
The outlet temperature is obtained from
Tout=PsyHFnTdbW(hout,wout)
where
Tout = outlet air temperature [^o^C],
PsyHFnTdbW(hout,wout) is an EnergyPlus psychrometric function.
The humidity ratio at saturation at the outlet temperature is
wout,sat=PsyWFnTdbRhPb(Tout,1.0,Patmo)
where
Patmo = the barometric pressure [Pa],
1.0 is the relative humidity at saturation,
PsyWFnTdbRhPb is an EnergyPlus psychrometric function.
IF wout≤wout,sat then the outlet condition is below the saturation curve and the desired moisture addition rate can be met. ˙mw,add is set to ˙mw,add,needed,max and the calculation of outlet conditions is done. But if wout>wout,sat then it is assumed that this condition will be detected and the steam addition rate throttled back to bring the outlet conditions back to the saturation condition. We need to find the point where the line drawn between state 1 (inlet) and state 2 (our desired outlet) crosses the saturation curve. This will be the new outlet condition. Rather than iterate to obtain this point, we find it approximately by solving for the point where 2 lines cross: the first drawn from state 1 to state 2, the second from T1 , w1,sat to T2 , w2,sat ; where
T1 is the inlet temperature [^o^C],
w1,sat is the humidity ratio at saturation at temperature T1 [kg/kg],
T2 is the desired outlet temperature [^o^C],
w2,sat is the humidity ratio at saturation at temperature T2 [kg/kg].
The 2 lines are given by the equations:
w=w1+((w2−w1)/(T2−T1))⋅(T−T1)
w=w1,sat+((w2,sat−w1,sat)/(T2−T1))⋅(T−T1)
Solving for the point (state 3) where the lines cross:
w3=w1+((w2−w1)⋅(w1,sat−w1))/(w2−w2,sat+w1,sat−w1)
T3=T1+(w3−w1)⋅((T2−T1)/(w2−w1))
This point isn’t quite on the saturation curve since we made a linear approximation of the curve, but the temperature should be very close to the correct outlet temperature. We will use this temperature as the outlet temperature and move to the saturation curve for the outlet humidity and enthalpy. Thus we set Tout = T3 and
wout=PsyWFnTdbRhPb(Tout,1.0,Patmo)
hout=PsyHFnTdbW(Tout,wout)
where PsyHFnTdbW is an EnergyPlus psychrometric function. The water addition rate is set to
˙mw,add=˙ma⋅(wout−win)
We now have the outlet conditions and the adjusted steam addition rate for the case where the desired outlet humidity results in an outlet state above the saturation curve.
Finally, the electricity consumption is given by
Whum=(˙mw,add/Capnom)⋅Wnom+Wfan+Wstby
where
Wfan = nominal fan power [W], a user input,
Wstby = standby power [W], a user input.
and the water consumption rate is
˙Vcons=˙mw,add/ρw
where
˙Vcons = the water consumption rate [m3],
ρw = water density (998.2 kg/m3).
ASHRAE. 1993. HVAC 2 Toolkit: A Toolkit for Secondary HVAC System Energy Calculations. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.
ASHRAE. 2001. 2001 ASHRAE Handbook Fundamentals. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.
Air System Humidifiers[LINK]
Overview[LINK]
Air system humidifiers are components that add moisture to the supply air stream. They fall into 2 broad categories: spray type humidifiers which act like direct evaporative coolers, cooling the supply air as well as humidifying it; and dry steam humidifiers, which humidify the supply air stream while causing almost no change to the supply air stream temperature. The EnergyPlus electric steam humidifier uses electrical energy to convert ordinary tap water to steam which it then injects into the supply air stream by means of a blower fan. The actual unit might be an electrode-type humidifier or a resistance-type humidifier.
Electric Steam Humidifier[LINK]
The electric steam humidifier model (object name: Humidifier:Steam:Electric) is based on moisture and enthalpy balance equations plus standard psychrometric relationships. The approach is similar to that taken in the ASHRAE HVAC 2 Toolkit, page 4-112 (ASHRAE 1993). EnergyPlus contains its own module of psychrometric routines; the psychrometric theory and relations are given in the 2001 edition of ASHRAE Fundamentals, Chapter 6 (ASHRAE 2001). The model contains both an ideal controller and the component. The control model assumes that there is a minimum humidity setpoint on the component air outlet node. This setpoint is established by a setpoint manager described elsewhere.
Model[LINK]
The component model is a forward model: its inputs are its inlet conditions; its outputs are its outlet conditions and its energy consumption. The inputs are the temperature, humidity ratio, and mass flow rate of the inlet air stream, which are known; and the water addition rate (kg/s) which is determined by the controller.
Controller[LINK]
The controller first decides whether the humidifier is on or off. For the humidifier to be on:
If the humidifier is off, the water addition rate is set to zero. If the humidifier is on, the water addition rate needed to meet the humidity setpoint is calculated.
˙ma⋅win+˙mw,add,needed=˙ma⋅wset
where
˙ma = the air mass flow rate [kg/s]
win = the inlet air humidity ratio [kg/kg]
˙mw,add,needed = water addition rate needed to meet the setpoint [kg/s]
wset = the humidity ratio setpoint [kg/kg]
Equation is the moisture balance equation for the component. It is solved for ˙mw,add,needed (the other variables are known) which is passed to the humidifier component model as its desired inlet water addition rate.
Component[LINK]
The inputs to the component model are the air inlet conditions and mass flow rate and the water addition rate set by the controller. The outputs are the air outlet conditions. First the desired water addition rate is checked against component capacity.
˙mw,add,needed,max=Min(˙mw,add,Capnom)
where
Capnom = the humidifier nominal capacity [kg/s], a user input.
If ˙mw,add,needed,max is zero, the outlet conditions are set to the inlet conditions and the water addition rate is set to zero. If the humidifier is scheduled on the component power consumption is set to the standby power consumption: Whum=Wstby . Otherwise Whum = 0.
If ˙mw,add,needed,max > 0, then the moisture and enthalpy balance equations
˙ma⋅win+˙mw=˙ma⋅wout
˙ma⋅hin+˙mw⋅hw=˙ma⋅hout
with ˙mw set equal to ˙mw,add,needed,max are solved for wout and hout . Here
˙ma = the air mass flow rate [kg/s]
win = the inlet air humidity ratio [kg/kg]
˙mw = the inlet water addition rate [kg/s]
wout = the outlet air humidity ratio [kg/kg]
hin = the inlet air specific enthalpy [J/kg]
hw = the steam specific enthalpy = 2676125. [J/kg] at 100 oC
hout = the outlet air specific enthalpy [J/kg]
The outlet temperature is obtained from
Tout=PsyHFnTdbW(hout,wout)
where
Tout = outlet air temperature [^o^C],
PsyHFnTdbW(hout,wout) is an EnergyPlus psychrometric function.
The humidity ratio at saturation at the outlet temperature is
wout,sat=PsyWFnTdbRhPb(Tout,1.0,Patmo)
where
Patmo = the barometric pressure [Pa],
1.0 is the relative humidity at saturation,
PsyWFnTdbRhPb is an EnergyPlus psychrometric function.
IF wout≤wout,sat then the outlet condition is below the saturation curve and the desired moisture addition rate can be met. ˙mw,add is set to ˙mw,add,needed,max and the calculation of outlet conditions is done. But if wout>wout,sat then it is assumed that this condition will be detected and the steam addition rate throttled back to bring the outlet conditions back to the saturation condition. We need to find the point where the line drawn between state 1 (inlet) and state 2 (our desired outlet) crosses the saturation curve. This will be the new outlet condition. Rather than iterate to obtain this point, we find it approximately by solving for the point where 2 lines cross: the first drawn from state 1 to state 2, the second from T1 , w1,sat to T2 , w2,sat ; where
T1 is the inlet temperature [^o^C],
w1,sat is the humidity ratio at saturation at temperature T1 [kg/kg],
T2 is the desired outlet temperature [^o^C],
w2,sat is the humidity ratio at saturation at temperature T2 [kg/kg].
The 2 lines are given by the equations:
w=w1+((w2−w1)/(T2−T1))⋅(T−T1)
w=w1,sat+((w2,sat−w1,sat)/(T2−T1))⋅(T−T1)
Solving for the point (state 3) where the lines cross:
w3=w1+((w2−w1)⋅(w1,sat−w1))/(w2−w2,sat+w1,sat−w1)
T3=T1+(w3−w1)⋅((T2−T1)/(w2−w1))
This point isn’t quite on the saturation curve since we made a linear approximation of the curve, but the temperature should be very close to the correct outlet temperature. We will use this temperature as the outlet temperature and move to the saturation curve for the outlet humidity and enthalpy. Thus we set Tout = T3 and
wout=PsyWFnTdbRhPb(Tout,1.0,Patmo)
hout=PsyHFnTdbW(Tout,wout)
where PsyHFnTdbW is an EnergyPlus psychrometric function. The water addition rate is set to
˙mw,add=˙ma⋅(wout−win)
We now have the outlet conditions and the adjusted steam addition rate for the case where the desired outlet humidity results in an outlet state above the saturation curve.
Finally, the electricity consumption is given by
Whum=(˙mw,add/Capnom)⋅Wnom+Wfan+Wstby
where
Wfan = nominal fan power [W], a user input,
Wstby = standby power [W], a user input.
and the water consumption rate is
˙Vcons=˙mw,add/ρw
where
˙Vcons = the water consumption rate [m3],
ρw = water density (998.2 kg/m3).
References[LINK]
ASHRAE. 1993. HVAC 2 Toolkit: A Toolkit for Secondary HVAC System Energy Calculations. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.
ASHRAE. 2001. 2001 ASHRAE Handbook Fundamentals. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.
Documentation content copyright © 1996-2014 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.