Component Sizing[LINK]
In EnergyPlus each HVAC component sizes itself. Each component module contains a sizing subroutine. When a component is called for the first time in a simulation, it reads in its user specified input data and then calls the sizing subroutine. This routine checks the autosizable input fields for missing data and calculates the data when needed.
A number of highlevel variables are used in the sizing subroutines.
CurDuctType (in DataSizing) contains the information about the current duct type. The types can be main, cooling, heating or other.
CurZoneEqNum (in DataSizing) is the current zone equipment set index and indicates that the component is a piece of zone equipment and should size itself using the zone sizing data arrays.
CurSysNum (in DataSizing) is the current air loop index and indicates that the component is part of the primary air system and should size itself using the system sizing data arrays.
Fan sizing is done in subroutine SizeFan.
Max Flow Rate[LINK]
If the fan is part of the central air system then check the duct type.
For duct type = main, other or default, ˙Vfan,max=DesMainVolFlowsys
For duct type = cooling, ˙Vfan,max=DesCoolVolFlowsys
For duct type = heating, ˙Vfan,max=DesHeatVolFlowsys
If the fan is zone equipment then check whether it is part of a component that only does heating.
For heating only ˙Vfan,max=DesHeatVolFlowzone ;
Otherwise ˙Vfan,max=max(DesHeatVolFlowzone,DesCoolVolFlowzone)
If the max fan flow rate is less than SmallAirVolFlow the max flow rate is set to zero.
Coil:Cooling:Water[LINK]
The sizing is done in function SizeWaterCoil of module WaterCoils
Initial Calculations[LINK]
For central cooling coils, the first step is to determine the design air flow rate, load, and design air entering and exit conditions. The coil design air flow rate is not generally the same as the maximum system air flow rate (used to size the central fans). The cooling coil peak load (either sensible or total) can occur at a different time than the system peak flow rate. Hence the coil air entering conditions can be different than those at the peak system flow rate. Also, the method of controlling the coil’s cooling output may also affect coil design flow rate as well as the coil design exit temperature and humidity.
By choosing Type of Load to Size On = Sensible or Total in Sizing:System the user indicates to the program to save the cooling coil air flow rate and system air conditions (mixed, return, outside) at the time of either the system cooling sensible or total load peak. Note that the choice VentilationRequirement uses the time of the sensible peak.
Choosing Central Cooling Capacity Control Method =
VAV,
Bypass,
VT, or
OnOff indicates which type of cooling output control the program should assume when calculating the design air flow rate. The function
GetCoilDesFlowT in module
ReportSizingManager calculates the air flow rate and exit air temperature for each capacity control method.
Cooling coil calculations for different capacity control methods
Control Method

Calculations

VAV

Tcc,exit=Tcool,supply˙Vcc,air=˙mcc,air,peakρair

Bypass

Tcc,exit=Tcool,supply˙Vcc,air=˙Vcc,air,max⋅max(0,min(1,Tmix,at−peak−Tsup,avgTmix,at−peak−Tcc,exit))

VT

Tcc,exit=max(Tcool,supply,Tsup,avg)˙Vcc,air=˙Vcc,air,max

OnOff

Tcc,exit=Tcool,supply˙Vcc,air=˙Vsys,air,max

Where:
Tsup,avg=Tzones,avg−∑zones˙Qsens,at−peakρaircp,air˙Vcool,air,max
and:
Cp,air : the specific heat of air (J/kgC)
˙mcc,air,peak : the air mass flow rate through the cooling coil at the sensible or total system peak cooling load (m3/s)
∑zones˙Qsens,at−peak : sum of the zone sensible cooling loads at the time of the peak system cooling load
ρair : the density of air (kg/m3)
Tcc,exit : the design cooling coil exit temperature (c)
Tcool,supply : the supply air temperature for cooling specified in Sizing:System (C)
Tmix,at−peak : the mixed air temperature at the time of the system peak cooling load (C)
Tzones,avg : the average zone temperature at the time of the system peak cooling load (C)
˙Vcc,air : the design volumetric air flow rate through the cooling coil (m3/s). This is the flow rate at either the sensible or total cooling load peak from the design period calculations.
˙Vcool,air,max : the maximum cooling volumetric air flow rate from the design calculations (m3/s). This flow rate occurs at the maximum zone cooling demand.
˙Vsys,air,max : the maximum volumetric air flow rate from the design calculations (m3/s). This flow rate occurs at either the maximum zone cooling or heating demand.
Design Coil Load  System Coils[LINK]
Design coil load (cooling capacity) is not an input for Coil:Cooling:Water. It is used for calculating the design water flow rate.
The design load is calculated as:
˙Qcoil,des=˙ma,coil,des(ha,coil,des,in−ha,coil,des,out)
Where:
ha,coil,des,in : is the coil design inlet air enthalpy (J/kg)
ha,coil,des,out : is the coil design outlet air enthalpy (J/kg), and
˙ma,coil,des : is the coil design air mass flow rate (kg/s)
The design air mass flow rate depends on the location of the coil. If the coil is in the outside air stream, the flow rate is set to
ρair˙Va,coil,oa,des where ˙Va,coil.oa,des is the design outside air volumetric flow rate for the system. Otherwise, it is set to:
ρair˙Vcc,air
where ˙Vcc,air is calculated above in the Initial Calculations section.
To obtain the inlet and outlet enthalpies, we need the inlet and outlet temperatures and humidity ratios. The inlet and outlet conditions depend on whether the coil is in the outside air stream and if it is not, whether or not there is outside air preconditioning.
Coil in outside air stream:
[LINK]
Tair,in,des=Tout,cool,at−peak (the outside air temperature at the design cooling peak)
Tair,out,des=Tsys,precool (the specified Precool Design Temperature from the System:Sizing object)
Wair,in,des=Wout,cool,at−peak (the outside humidity ratio at the design cooling peak)
Wair,out,des=Wsys,precool (the specified Precool Design Humidity Ratio from the System:Sizing object)
Coil in main air stream, no preconditioning of outside air
[LINK]
Tair,in,des=Tmix,cool,at−peak (the mixed air temperature at the design cooling peak)
Wair,in,des=Wmix,cool,at−peak (the mixed humidity ratio at the design cooling peak)
Tair,out,des=Tcc,exit (calculated above in the Initial Calculation section)
Wair,out,des=Wsup,cool (the specified Central Cooling Design Supply Air Humidity Ratio from the Sizing:System object)
Coil in main air stream, outside air preconditioned
[LINK]
The oustide air fraction is calculated as (where V_{cc,air} is calculated as above)
foa=˙Vair,out,des˙Vcc,air
Tair,in,des=foaTprecool+(1−foa)Tret,cool,at−peak (Precool temperature is the specified Precool Design Temperature from System:Sizing Manager; T_ret_cool_atpeak is the return temperature at the system cooling peak load)
Wair,in,des=foaWprecool+(1−foa)Wret,cool,at−peak (Precool humidity ratio is the specified Precool Design Humidity Ratio from System:Sizing Manager; W_ret_cool_atpeak is the return humidity ratio at the system cooling peak load)
Tair,out,des=Tcc,exit (calculated above in the Initial Calculation section)
Wair,out,des=Wsup,cool (the specified Central Cooling Design Supply Air Humidity Ratio from the Sizing:System object)
With the inlet and outlet conditions established, we can obtain the inlet and outlet enthalpies:
hair,coil,des,in=PsyHFnTdbW(Tair,in,des,Wair,in,des)hair,coil,des,out=PsyHFnTdbW(Tair,out,des,Wair,out,des)
Where PsyHFnTdbW is the EnergyPlus function for calculation air specific enthalpy given the air temperature and humidity ratio. We now have all we need to calculate the design coil capacity, ˙Qcoil,des .
Design Coil Load  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit or an ZoneHVAC:FourPipeFanCoil, the cooling load (cooling capacity) is passed down from the terminal unit or fan coil sizing calculations. Otherwise the load is defined as:
˙Qcoil,des=˙ma,coil,des(ha,coil,des,in−ha,coil,des,out)
Where:
ha,coil,des,in : is the coil design inlet air enthalpy (J/kg)
ha,coil,des,out : is the coil design outlet air enthalpy (J/kg), and
˙ma,coil,des : is the coil design air mass flow rate (kg/s)
The enthalpies are given by:
hair,coil,des,in=PsyHFnTdbW(Tair,in,des,Wair,in,des)hair,coil,des,out=PsyHFnTdbW(Tair,out,des,Wair,out,des)
Where the inputs to those functions are the coil inlet design conditions. For coils in terminal units these are set at the system level to the system design supply air temperature. For zonal units they are set to design return air, mixed air, or outside air as appropriate to the unit. T_{air,out,des} is set to the zone cooling design supply air temperature as specified in the Zone:Sizing inputs. W_{air,out,des} is set to the zone cooling design supply air humidity ratio as specified in the Zone:Sizing inputs.
Design Water Flow Rate (m^{3}/s)  System Coils[LINK]
The design water volumetric flow rate is calculated using:
˙Vw,coil,des=˙Qcoil,desρwcp,wΔTw,des
Where ΔTw,des is just the Loop Design Temperature Difference user input from Sizing:Plant (if the coil is in the outside air stream, 1/2 the Loop Design Temperature Difference is used). The design coil load Load_{coil,des} is calculated from:
Design Water Flow Rate (m^{3}/s)  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit or an ZoneHVAC:FourPipeFanCoil, the chilled water flow rate is passed down from the terminal unit or fan coil sizing calculations. Otherwise the flow is set to:
˙Vw,coil,des=˙Qcoil,desρwcp,wΔTw,des
Where ΔTw,des is just the Loop Design Temperature Difference user input from Sizing:Plant.
Design Air Flow Rate  System Coils[LINK]
The design air volumetric flow rate for the system cooling coil is set to:
the design outside air flow rate if the coil is in the outside air stream;
the coil design flow rate from function GetCoilDesFlowT described in section “Initial Calculations”;
the design flow rate set by the parent component (such as a unitary system) containing the cooling coil.
Design Air Flow Rate  Zone Coils[LINK]
Zone chilled water coils are always part of a zone HVAC component. In almost all cases the design flow rate is passed down from the design flow rate of the parent component. Otherwise if the parent component does cooling only the flow rate for the coil is set to the zone design cooling flow rate. And if the parent component does both cooling and heating, the coil flow rate is set to the maximum of the zone design cooling and heating flow rates.
Design Air Inlet Temperature  System Coils[LINK]
The inlet air temperature depends on whether the coil is in the outside air stream and if it is not, whether or not there is outside air preconditioning.
Coil in outside air stream: Tair,in,des=Tout,cool,at−peak (the outside air temperature set at the design cooling peak).
Coil in main air stream, no preconditioning of outside air: Tair,in,des=Tmix,cool,at−peak (the mixed air temperature at the cooling design peak).
Coil in main air stream, outside air preconditioned. The outside air fraction is calculated as foa=˙Vair,out,des/˙Vcc,air , where ˙Vcc,air is calculated above. Then Tair,in,des=foaTprecool+(1−foa)Tret,cool,at−peak , where Tprecool is the specified Precool Design Temperature from System:Sizing, and Tret,cool,at−peak is the return temperature at the system cooling peak load.
Design Air Inlet Temperature  Zone Coils[LINK]
The design inlet temperature depends on whether the coil is in a terminal unit or a zonal unit, and where the coil is positioned within the unit.
For the AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction terminal unit the design inlet temperature is set to the zone temperature at the time of the zone cooling peak, since the coil is located in the induced air stream.
For fan coil units the design inlet temperature is set to the mixed air temperature: Tair,in,des=foaToa,coolpeak+(1−foa)Tz,coolpeak , where foa=ρa˙Vz,oa,des/˙mz,cool,des
In all other cases the design inlet temperature is set to the zone design cooling coil inlet temperature which is calculated in the zone sizing simulation and is basically the same calculation as the fan coil unit.
Design Air Outlet Temperature  System Coils[LINK]
The outlet air temperature depends on whether the coil is in the outside air stream.
Coil in outside air stream: T_{air,out,des} = T_{sys,des,precool} (the specified Precool Design Temperature from the Sizing:System object).
Coil in main air stream: the design outlet air temperature is set to the temperature calculated in the Initial Calculation section above.
Design Air Outlet Temperature  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit, then:
˙Qcoil,des=cp,airρair˙Vw,coil,desΔTw,desT1=Tair,in,des−˙Qcoil,des/(ρaircp,air˙Vair,coil,des)T2=Tw,out,des+2Tair,out,des=max(T1,T2)
For all other cases T_{air,out,des} is set to T_{z,sup,des} (the zone design supply air temperature as specified in Sizing:Zone).
Design Inlet Air Humidity Ratio  System Coils[LINK]
The design inlet humidity ratio depends on whether the coil is in the outside air stream and if it is not, whether or not there is outside air preconditioning.
Coil in outside air stream: Wair,in,des=Wout,cool,at−peak (the outside air humidity ratio at the design cooling peak).
Coil in main air stream, no preconditioning of outside air: Wair,in,des=Wmix,cool,at−peak (the mixed air humidity ratio at the cooling design peak).
Coil in main air stream, outside air preconditioned. The outside air fraction is calculated as foa=˙Vair,out,des/˙Vcc,air , where ˙Vcc,air is calculated above. Then Wair,in,des=foaWprecool+(1−foa)Wret,cool,at−peak , where Wprecool is the specified Precool Design Humidity Ratio from System:Sizing, and Wret,cool,at−peak is the return humidity ratio at the system cooling peak load.
Design Air Inlet Humidity Ratio  Zone Coils[LINK]
The design inlet humidity ratio depends on whether the coil is in a terminal unit or a zonal unit, and where the coil is positioned within the unit.
For the AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction terminal unit the design inlet humidity ratio is set to the zone humidity ratio at the time of the zone cooling peak, since the coil is located in the induced air stream.
For fan coil units the design inlet humidity ratio is set to the mixed air humidity ratio: Wair,in,des=foaWoa,coolpeak+(1−foa)Wz,coolpeak , where foa=ρa˙Vz,oa,des/˙mz,cool,des
In all other cases the design inlet humidity ratio is set to the zone design cooling coil inlet hunidity ratio which is calculated in the zone sizing simulation and is basically the same calculation as the fan coil unit.
Design Outlet Air Humidity Ratio  System Coils[LINK]
The outlet air humidity ratio depends on whether the coil is in the outside air stream.
Coil in outside air stream: W_{air,out,des} = W_{sys,des,precool} (the specified Precool Design Humidity Ratio from the Sizing:System object)
Coil in main air stream: W_{air,out,des} = PsyWFnTdbRhPb(T_{air,out,des},0.9,P_{air,std}), where PsyWFnTdbRhPb is the EnergyPlus psychrometric function to calculate humidity ratio from drybulb temperature, relative humidity, and atmospheric pressure. The design outlet humidity ratio is being set to the humidity ratio at 90% relative humidity and design outlet temperature.
Design Outlet Air Humidity Ratio  Zone Coils[LINK]
Design Inlet Water Temperature  System Coils[LINK]
The Design Inlet Water Temperature is set to the Design Loop Exit Temperature specified in the Sizing:Plant object for the water loop serving this coil.
Design Inlet Water Temperature  Zone Coils[LINK]
The Design Inlet Water Temperature is set to the Design Loop Exit Temperature specified in the Sizing:Plant object for the water loop serving this coil.
Coil:Cooling:Water:DetailedGeometry Sizing[LINK]
The sizing is done in subroutine SizeWaterCoil
Max Water Flow Rate of Coil[LINK]
The calculation is identical to that done for Coil:Cooling:Water.
Number of Tubes per Row[LINK]
Ntube/row=Int(13750⋅˙Vcoil,water,max)
N_{tube/row} = Max(N_{tube/row},3)
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙mair,des=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙mair,des=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙mair,des=ρair⋅DesHeatVolFlowsys
Dfin=0.335⋅˙mair,des
Minimum Air Flow Area[LINK]
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙mair,des=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙mair,des=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙mair,des=ρair⋅DesHeatVolFlowsys
AMinAirFlow=0.44⋅˙mair,des
Fin Surface Area[LINK]
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙mair,des=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙mair,des=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙mair,des=ρair⋅DesHeatVolFlowsys
AFinSurf=78.5⋅˙mair,des
Total Tube Inside Area[LINK]
A_{tube,total\ inside} = 4.4 *D_{tube,inside} *N_{tube\ rows} *N_{tubes/row}
Where D_{tube,inside} is the tube inside diameter.
Tube Outside Surf Area[LINK]
A_{tube,outside} = 4.1 *D_{tube,outside} *N_{tube\ rows} *N_{tubes/row}
Where D_{tube,outside} is the tube outside diameter.
Depth_{coil} = Depth_{tube\ spacing} * N_{tube\ rows}
Coil:Cooling:WaterToAirHeatPump:EquationFit Sizing[LINK]
The sizing is done in subroutine SizeHVACWaterToAir
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water, which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. If there is a companion heating coil, the heating coil design load is used so that both modes will have the same rated water flow rate. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Cooling Capacity[LINK]
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:Water. The following calculations are then performed to determine the rated total cooling capacity.
TWB,ratio={\left( {{T_{WB,air,in,des}} + 273.15\,C} \right)}/(TWB,air,in,des+273.15C)283.15C{283.15\,C}
TS,ratio={\left( {29.44\,C + 273.15\,C} \right)}/(29.44C+273.15C)283.15C{283.15\,C}
where:
${T_{WB,ratio}} = $ ratio of loadside inlet air wetbulb temperature in Kelvin to a reference temperature
${T_{S,ratio}} = $ ratio of sourceside inlet water temperature in Kelvin to a reference temperature
TCC1 = user input for Total Cooling Capacity Coefficient 1
TCC2 = user input for Total Cooling Capacity Coefficient 2
TCC3 = user input for Total Cooling Capacity Coefficient 3
TCC4 = user input for Total Cooling Capacity Coefficient 4
TCC5 = user input for Total Cooling Capacity Coefficient 5
TotCapTempModFac=TCC1+TCC2(TWB,ratio)+TCC3(TS,ratio)+TCC4+TCC5
The 4^{th} and 5^{th} coefficient (TCC4 and TCC5) used in the above equation are multipliers for the loadside and sourceside flow ratios, respectively. For sizing, these ratios are assumed to be 1.
The enthalpy of the entering air is then compared with the enthalpy of the exiting air. The calculations for air enthalpy are identical to that done for Coil:Cooling:Water. If the entering air enthalpy is less than the exiting air enthalpy, a reference value of 48,000 J/kg is used as the entering air enthalpy. If the TotCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the design air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil total cooling capacity is set equal to 0.
IF(Hin>Hout)THEN⋅Qcoil,des,total={{m_{air,des}}\left( {{H_{in}}  {H_{out}}} \right)}/mair,des(Hin−Hout)TotCapTempModFac{TotCapTempModFac}ELSE⋅Qcoil,des,total={{m_{air,des}}\left( {48000  {H_{out}}} \right)}/mair,des(48000−Hout)TotCapTempModFac{TotCapTempModFac}ENDIF
Rated Sensible Cooling Capacity[LINK]
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:Water. The following calculations are then performed to determine the rated sensible cooling capacity.
TDB,ratio={\left( {{T_{DB,air,in,des}} + 273.15\,C} \right)}/(TDB,air,in,des+273.15C)283.15C{283.15\,C}
TS,ratio={\left( {29.44\,C + 273.15\,C} \right)}/(29.44C+273.15C)283.15C{283.15\,C}
where:
${T_{DB,ratio}} = $ ratio of loadside inlet air drybulb temperature in Kelvin to a reference temperature
SCC1 = user input for Sensible Cooling Capacity Coefficient 1
SCC2 = user input for Sensible Cooling Capacity Coefficient 2
SCC3 = user input for Sensible Cooling Capacity Coefficient 3
SCC4 = user input for Sensible Cooling Capacity Coefficient 4
SCC5 = user input for Sensible Cooling Capacity Coefficient 5
SCC6 = user input for Sensible Cooling Capacity Coefficient 6
SensCapTempModFac=SCC1+SCC2(TDB,ratio)+SCC3(TWB,ratio)+SCC4(TS,ratio)+SCC5+SCC6
The 5^{th} and 6^{th} coefficient (SCC5 and SCC6) used in the above equation are multipliers for the loadside and sourceside flow ratios, respectively. For sizing, these ratios are assumed to be 1.
The drybulb temperature of the entering air is then compared with the drybulb temperature of the exiting air. The calculations for air drybulb temperature are identical to that done for Coil:Cooling:Water. If the entering air drybulb temperature is less than the exiting air drybulb temperature, a reference value of 24 C is used as the entering air drybulb temperature. If the SensCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the design air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil sensible cooling capacity is set equal to 0.
IF(TDB,in>TDB,out)THEN⋅Qcoil,des,sensible={{m_{air,des}}C{p_{air,des}}\left( {{T_{DB,in}}  {T_{DB,out}}} \right)}/mair,desCpair,des(TDB,in−TDB,out)SensCapTempModFac{SensCapTempModFac}ELSE⋅Qcoil,des,sensible={{m_{air,des}}C{p_{air,des}}\left( {24  {T_{DB,out}}} \right)}/mair,desCpair,des(24−TDB,out)SensCapTempModFac{SensCapTempModFac}ENDIF
Coil:Cooling:WaterToAirHeatPump:VariableSpeedEquationFit Sizing[LINK]
For the cooling coil of VS WSHP, we specify a nominal speed level. During the sizing calculation, the Rated Air Volume Flow Rate, the Rated Water Volume Flow Rate and the Rated Total Cooling Capacity at the Selected Nominal Speed Level are determined in the same way as the Coil:Cooling:WaterToAirHeatPump:EquationFit object. The sensible heat transfer rate is not allowed for autosizing, instead, it is a function of the rated air and water flow rates, rated total cooling capacity and the Reference Unit SHR at the nominal speed level. The default nominal speed level is the highest speed. However, the model allows the user to select a nominal speed level rather than the highest.
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit , which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. If there is a companion heating coil, the heating coil design load is used so that both modes will have the same rated water flow rate. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Cooling Capacity[LINK]
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit. The calculations for air enthalpy are similar to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit. The difference is in calculating the total cooling capacity temperature modifier function at the selected nominal speed level, as below:
TotCapTempModFracNominalSpeed=a+b∗WBi+c∗WB2i+d∗EWT+e∗EWT2+f∗WBi∗EWT
where
WB_{i} = wetbulb temperature of the air entering the heating coil, degC
EWT = entering water temperature, degC
af = regression curvefit coefficients.
If the entering air enthalpy is less than the exiting air enthalpy, a reference value of 48,000 J/kg is used as the entering air enthalpy. If the TotCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the rated air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil total cooling capacity is set equal to 0.
If H_{in} > H_{out} Then
˙Qcoil,rated,total=mair,rated(Hin−Hout)/TotCapTempModFracNominalSpeed
Else
˙Qcoil,rated,total=mair,rated(48000−Hout)/TotCapTempModFracNominalSpeed
End If
Coil:Heating:WaterToAirHeatPump:EquationFit Sizing[LINK]
The sizing is done in subroutine SizeHVACWaterToAir.
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water , which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Heating Capacity[LINK]
The rated total heating capacity is set equal to the rated total cooling capacity.
Coil:Heating:WaterToAirHeatPump:VariableSpeedEquationFit Sizing[LINK]
For the heating coil of VS WSHP, we specify a nominal speed level. During the sizing calculation, the Rated Air Volume Flow Rate and the Rated Water Volume Flow Rate are determined in the same way as the Coil:Heating:WaterToAirHeatPump:EquationFit object. On the other hand, the Rated Heating Capacity at the Selected Nominal Speed Level should be the same as the total cooling capacity of its corresponding cooling coil, which has to be sized first. The default nominal speed level will be the highest speed. However, the model allows the user to select a nominal speed level rather than the highest.
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit, which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Heating Capacity[LINK]
The rated total heating capacity is set equal to the rated total cooling capacity.
Coil:Heating:Water Sizing[LINK]
The sizing is done in subroutine SizeWaterCoil.
Max Water Flow Rate of Coil[LINK]
With the coil load from the system design data array and the user specified (in a Sizing:Plant object) design hot water temperature fall, calculate the max water flow rate:
˙Vcoil,water,max=HeatCapsys/(Cp,water⋅ρwater⋅ΔTplt,hw,des)
Using the zone design coil inlet and supply air conditions calculate the design coil load.
If the coil is not part of an induction unit then obtain the coil inlet temperature from the zone design data array;
T_{in,air} = DesHeatCoilInTemp_{zone}
If the coil is part of an induction unit take into account the induced air:
Frac_{minflow} = MinFlowFrac_{zone}
T_{in,air} = DesHeatCoilInTemp_{zone} * Frac_{minflow} +
ZoneTempAtHeatPeak_{ zone} *(1 Frac_{minflow})
T_{out,air} = HeatDesTemp_{zone}
W_{out,air} = HeatDesHumRat_{zone}
If the coil is part of a terminal unit the mass flow rate is determined by the volumetric flow rate of the terminal unit:
˙mair,des=ρair⋅˙mair,des,tu
Otherwise the design flow is obtained from the zone design data array:
˙mair,des=DesHeatMassFlowzone
Qcoil,des=cp,air˙mair,des⋅(Tout,air−Tin,air)
Here c_{p,air} is calculated at the outlet humidity and the average of the inlet and outlet temperatures.
With the coil load and the user specified (in a Sizing:Plant object) design hot water temperature decrease, calculate the max water flow rate:
˙Vcoil,water,max=Qcoil,des/(Cp,water⋅ρwater⋅ΔTplt,hw,des)
UA of the Coil[LINK]
To obtain the UA of the coil, we specify the model inputs (other than the UA) at design conditions and the design coil load that the coil must meet. Then we numerically invert the coil model to solve for the UA that will enable the coil to meet the design coil load given the specified inputs.
The design coil load is the system design sensible cooling capacity;
Q_{coil,des} = HeatCap_{sys}
The required inputs for the simple coil model are:
T_{in,air} = HeatMixTemp_{sys}
W_{in,air} = HeatMixHumRat_{sys}
T_{in,water} = ExitTemp_{plt,hw,des}
˙min,water=ρwater⋅˙Vcoil,water,max
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙min,air=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙min,air=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙min,air=ρair⋅DesHeatVolFlowsys
We now have all the data needed to obtain UA. The numerical inversion is carried out by calling subroutine SolveRegulaFalsi. This is a general utility routine for finding the zero of a function. In this case it finds the UA that will zero the residual function  the difference between the design coil load and the coil output divided by the design coil load. The residual is calculated in the function SimpleHeatingCoilUAResidual.
If the coil is not part of an induction unit then obtain the coil inlet temperature from the zone design data array;
T_{in,air} = DesHeatCoilInTemp_{zone}
If the coil is part of an induction unit take into account the induced air:
Frac_{minflow} = MinFlowFrac_{zone}
T_{in,air} = DesHeatCoilInTemp_{zone} * Frac_{minflow} +
ZoneTempAtHeatPeak_{ zone} *(1 Frac_{minflow})
W_{in,air} = DesHeatCoilInHumRat_{zone}
T_{in,water} = ExitTemp_{plt,hw,des}
˙min,water=ρwater⋅˙Vcoil,water,max
T_{out,air} = HeatDesTemp_{zone}
W_{out,air} = HeatDesHumRat_{zone}
If the coil is part of a terminal unit the mass flow rate is determined by the volumetric flow rate of the terminal unit:
˙mair,des=ρair⋅˙mair,des,tu
Otherwise the design flow is obtained from the zone design data array:
˙mair,des=DesHeatMassFlowzone
˙Qcoil,des=cp,air⋅˙mair,des⋅(Tout,air−Tin,air)
Here c_{p,air} is calculated at the outlet humidity and the average of the inlet and outlet temperatures.
We now have all the data needed to obtain UA. The numerical inversion is carried out by calling subroutine SolveRegulaFalsi. This is a general utility routine for finding the zero of a function. In this case it finds the UA that will zero the residual function  the difference between the design coil load and the coil output divided by the design coil load. The residual is calculated in the function SimpleHeatingCoilUAResidual.
Coil:Heating:Steam Sizing[LINK]
The sizing is done in subroutine SizeSteamCoil.
Maximum Steam Flow Rate[LINK]
The maximum steam volumetric flow rate is calculated using:
˙Vcoil,steam,max=Loadcoil,desρsteam(hfg+cp,w⋅ΔTsc)
The steam density (ρsteam ) is for saturated steam at 100 degC (101325.0 Pa) and h_{fg} is the latent heat of vaporization of water at 100 degC (101325.0 Pa). C_{p,w} is the heat capacity of saturated water (condensate) at 100 degC (101325.0 Pa) and ΔTsc is the Degree of Subcooling defined in the Coil:Heating:Steam object input. The design coil load Load_{coil,des} is calculated from:
Loadcoil,des=˙mair,des(cp,air)(Tair,coil,des,out−Tair,coil,des,in)
The design air mass flow rate depends on the location of the coil (duct type). For duct type = main, the flow rate is set to r_{air} *DesMainVolFlow_{sys} *MinSysAirFlowRatio. If the coil is in a cooling duct the flow rate is set to r_{air} *DesCoolVolFlow_{sys} *MinSysAirFlowRatio. If the coil is in a heating duct the flow rate is set to r_{air} *DesHeatVolFlow_{sys}. If the coil is in any other kind of duct, the flow rate is set to r_{air} *DesMainVolFlow_{sys}.
For sizing, the design outlet air temperature (T_{air,coil,des,out}) is the Central Heating Design Supply Air Temperature specified in the Sizing:System object.
The design inlet air temperature depends on whether the coil is being sized for 100% outdoor air or minimum outdoor air flow (per 100% Outdoor Air in Heating input field in the Sizing:System object).
 Sizing based on 100% Outdoor Air in Heating
T_{air,coil,des,in} = HeatOutTemp_{sys} (the outdoor air temperature at the design heating peak)
 Sizing based on minimum outdoor air flow. The outdoor air fraction is calculated as Frac_{oa} = DesOutAirVolFlow_{sys} / DesVolFlow. DesVolFlow is ∙mair,des/∙mair,desρairρair .
T_{air,coil,des,in} = Frac_{oa}* * HeatOutTemp_{sys}* + (1. Frac_{oa}) *HeatRetTemp_{sys} (see Table 41. System Sizing Data)
If the coil is part of an AirTerminal:SingleDuct:* unit (e.g., AirTerminal:SingleDuct:ConstantVolume:Reheat, AirTerminal:SingleDuct:VAV:Reheat, AirTerminal:SingleDuct:SeriesPIU:Reheat, etc.), the maximum steam flow rate is set equal to the terminal unit’s maximum steam flow rate. Otherwise (e.g., the zonelevel coil is part of ZoneHVAC:PackagedTerminalAirConditioner, ZoneHVAC:UnitVentilator, ZoneHVAC:UnitHeater or ZoneHVAC:VentilatedSlab) the calculation is similar to that at the system level. A design load is calculated:
Loadcoil,des=˙mair,des(cp,air)(Tair,coil,des,out−Tair,coil,des,in)
where:
˙mair,des = DesHeatMassFlow_{zone} (see Table 40. Zone Sizing Data)
T_{air,coil,des,in} = DesHeatCoilInTemp_{zone} (see Table 40)
T_{air,coil,des,out} = HeatDesTemp_{zone} (user input from Sizing:Zone object)
cp,air = Specific heat of air (evaluated at the average of inlet and outlet air temperatures, and at the zone heating design supply air humidity ratio HeatDesHumRat_{zone} [user input from Sizing:Zone object])
˙Vcoil,steam,max=Loadcoil,desρsteam(hfg+cp,w⋅ΔTsc)
The terms in the denominator of this equation (ρ_{steam}, h_{fg}, etc.) are evaluated in the same way as described above for steam System Coils.
Sizing of Gas and Electric Heating Coils[LINK]
The sizing calculation is done in subroutine SizeHeatingCoil in module HeatingCoils.
Nominal Capacity of the Coil[LINK]
The value is obtained from the system design array.
Cap_{nom} = HeatCap_{sys}
The capacity is calculated from the design coil inlet and outlet conditions.
If the coil is not part of an induction unit then obtain the coil inlet temperature from the zone design data array;
T_{in,air} = DesHeatCoilInTemp_{zone}
If the coil is part of an induction unit take into account the induced air:
Frac_{minflow} = MinFlowFrac_{zone}
T_{in,air} = DesHeatCoilInTemp_{zone} * Frac_{minflow} +
ZoneTempAtHeatPeak_{ zone} *(1 Frac_{minflow})
T_{out,air} = HeatDesTemp_{zone}
W_{out,air} = HeatDesHumRat_{zone}
Q_{coil,des} = C_{p,air} * DesHeatMassFlow_{zone} *(T_{out,air}T_{in,air})
Here c_{p,air} is calculated at the outlet humidity and the average of the inlet and outlet temperatures.
DX Coil Sizing[LINK]
The sizing calculations are done in subroutine SizeDXCoil in module DXCoils. This section covers the sizing of the objects
Coil:Cooling:DX:SingleSpeed
Coil:Heating:DX:SingleSpeed
Coil:Cooling:DX:TwoSpeed
Rated Air Volume Flow Rate[LINK]
The rated air flow rate is obtained from the system design array.
˙Vair,rated=DesMainVolFlowsys
The rated air flow rate is the maximum of the heating and cooling design flow rates from the zone design array.
˙Vair,rated=Max(DesCoolVolFlowzone,DesHeatVolFlowzone)
Rated Total Cooling Capacity[LINK]
The rated cooling capacity is obtained by dividing the peak cooling capacity by the Cooling Capacity Modifier Curve evaluated at peak mixed wetbulb and outdoor drybulb temperatures.
T_{mix} = CoolMixTemp_{sys}
W_{mix} = CoolMixHumRat_{sys}
T_{sup} = CoolSupTemp_{sys}
W_{sup} = CoolSupHumRat_{sys}
T_{outside} = CoolOutTemp_{sys}
r_{air} = PsyRhoAirFnPbTdbW(p_{air,std}, T_{mix},W_{mix})
h_{mix} = PsyHFnTdbW(T_{mix},W_{mix})
h_{sup} = PsyHFnTdbW(T_{sup},W_{sup})
T_{mix,wb} = PsyTwbFnTdbWPb(T_{mix},W_{mix}, p_{air,std})
CapModFac = CurveValue(CCapFTemp,T_{mix,wb},T_{outside})
CCappeak=ρair⋅˙Vair,rated⋅(hmix−hsup)
CCap_{rated} = CCap_{peak}/ CapModFac
We check that the design volume flow per total capacity is within the prescribed range:
FlowCapRatio=˙Vair,rated/CCaprated
If FlowCapRatio < FlowCapRatio_{min} then
CCaprated=˙Vair,rated/FlowCapRatiomin
If FlowCapRatio > FlowCapRatio_{max} then
CCaprated=˙Vair,rated/FlowCapRatiomax
where
FlowCapRatio_{min} = 0.00004027 m^{3}/s per watt (300 cfm/ton)
And
FlowCapRatio_{max} = 0.00006041 m^{3}/s per watt (450 cfm/ton)
The sizing calculation for DX cooling coils for 100% dedicated outdor air system (DOAS) are identical to regular DX cooling coils. However, they operate operate at different flow to capacity ratio ranges and are within the prescribed range below:
FlowCapRatio_{min} = 0.00001677 m^{3}/s per Watt (125 cfm/ton)
And
FlowCapRatio_{max} = 0.00003355 m^{3}/s per Watt (250 cfm/ton)
The rated cooling capacity for zone coils is calculated in the same manner as for system coils.
T_{mix} = DesCoolCoilInTemp_{zone}
W_{mix} = DesCoolCoilInHumRat_{zone}
T_{sup} = CoolDesTemp_{zone}
W_{sup} = CoolDesHumRat_{zone}
T_{outside} = T_{outside},_{desday,peak}
r_{air} = PsyRhoAirFnPbTdbW(p_{air,std}, T_{mix},W_{mix})
h_{mix} = PsyHFnTdbW(T_{mix},W_{mix})
h_{sup} = PsyHFnTdbW(T_{sup},W_{sup})
T_{mix,wb} = PsyTwbFnTdbWPb(T_{mix},W_{mix}, p_{air,std})
CapModFac = CurveValue(CCapFTemp,T_{mix,wb},T_{outside})
CCappeak=ρair⋅˙Vair,rated⋅(hmix−hsup)
CCap_{rated} = CCap_{peak}/ CapModFac
We check that the design volume flow per total capacity is within the prescribed range:
FlowCapRatio=˙Vair,rated/CCaprated
If FlowCapRatio < FlowCapRatio_{min} then
CCaprated=˙Vair,rated/FlowCapRatiomin
If FlowCapRatio > FlowCapRatio_{max} then
CCaprated=˙Vair,rated/FlowCapRatiomax
where
FlowCapRatio_{min} = 0.00004027 m^{3}/s per watt (300 cfm/ton)
And
FlowCapRatio_{max} = 0.00006041 m^{3}/s per watt (450 cfm/ton)
We check the design flow to the total cooling capacity rato for dedicated zone outdoor unit DX cooling coils to be within the limits prescribed below:
FlowCapRatio_{min} = 0.00001677 m^{3}/s per Watt (125 cfm/ton)
And
FlowCapRatio_{max} = 0.00003355 m^{3}/s per Watt (250 cfm/ton)
Rated Total Heating Capacity[LINK]
For Coil:Heating:DX:SingleSpeed the rated heating capacity is set equal to the cooling capacity.
The rated sensible heat ratio is calculated to be the sensible cooling (from rated inlet conditions to user specified supply conditions) divided by the total cooling (from rated inlet to specified supply).
T_{in,rated} = 26.6667 ^{o}C (80 ^{o}F)
W_{in,rated} = 0.01125 (corresponds to 80 ^{o}F drybulb, 67 ^{o}F wetbulb)
C_{p,air} = PsyCpAirFnWTdb(W_{in,rated}, T_{in,rated})
For system coils
T_{sup} = CoolSupTemp_{sys}
W_{sup} = CoolSupHumRat_{sys}
For zone coils
T_{sup} = CoolDesTemp_{zone}
W_{sup} = CoolDesHumRat_{zone}
Then
h_{rated} = PsyHFnTdbW(T_{in,rated}, W_{in,rated})
h_{sup} = PsyHFnTdbW(T_{sup}, W_{sup})
Dh_{rated,sup} = h_{rated}h_{sup}
DQs_{rated,sup} = C_{p,air} *(T_{in,rated}T_{sup})
SHR_{rated} = DQs_{rated,sup}/Dh_{rated,sup}
Evaporative Condenser Air Volume Flow Rate[LINK]
The evaporative condenser air volume flow rate (m^{3}/s) is set to 0.000114 m^{3}/s per watt (850 cfm/ton) times the total rated cooling capacity.
Evaporative Condenser Air Volume Flow Rate, Low Speed[LINK]
The evaporative condenser air volume flow rate, low speed (m^{3}/s) is set to 1/3 times 0.000114 m^{3}/s per watt (850 cfm/ton) times the total rated cooling capacity.
Evaporative Condenser Pump Rated Power Consumption[LINK]
The evaporative condenser pump rated power consumption is set equal to the total cooling capacity times 0.004266 watts pump power per watt capacity (15 W/ton).
Evaporative Condenser Pump Rated Power Consumption, Low Speed[LINK]
The evaporative condenser pump rated power consumption, low speed, is set equal to 1/3 times the total cooling capacity times 0.004266 watts pump power per watt capacity (15 W/ton).
Rated Air Volume Flow Rate, low speed[LINK]
The rated air volume flow rate, low speed, is set equal to 1/3 times the full rated air volume flow rate.
Rated Total Cooling Capacity, Low Speed[LINK]
The rated total cooling capacity, low speed, is set equal to 1/3 times the full rated total cooling capacity.
Rated SHR, low speed[LINK]
The rated sensible heat ratio, low speed, is set equal to the full speed SHR.
Resistive Defrost Heater Capacity[LINK]
For the heat pump the resistive defrost heat capacity is set equal to the cooling capacity.
DX MultiSpeed Coil Sizing[LINK]
The sizing calculations are done in subroutine SizeDXCoil in module DXCoils. This section covers the sizing of the objects
The rated air volume flow rate, rated total cooling capacity, rated heating capacity, rated SHR, evaporative condenser air volume flow rate, evaporative condenser pump rated power consumption at the highest speed are sized in the same ways as DX Coil Sizing.
After the sizes are determined at the highest speed, the sizes in the rest of speeds are assumed to
Valuen=nNumberOfSpeed∗ValueNumberOfSpeed
where
Value_{n} = Any autosizable variable at Speed n, except SHR
SHR_{n} = SHR_{NumberOfSpeed}
n = Speed Index number from 1 to NumberOfSpeed1
NumberOfSpeed = The highest speed number
Coil:Cooling:DX:VariableSpeed Sizing[LINK]
For the variablespeed DX cooling coil, we specify a nominal speed level. During the sizing calculation, the Rated Total Cooling Capacity at the Selected Nominal Speed Level is determined in the same way as the Coil:Cooling:DX:SingleSpeed object. If the user chooses to autosize the Rated Air Volume Flow Rate, the flow rate, as compared to the Rated Total Cooling Capacity, is sized to have the same ratio as the air volume flow rate to the total cooling capacity at the nominal speed, of the Reference Unit. The sensible heat transfer rate is not allowed for autosizing, instead, it is a function of the rated air flow, rated total cooling capacity and the Reference Unit SHR at the nominal speed level. The default nominal speed level is the highest speed. However, the model allows the user to select a nominal speed level rather than the highest.
Rated Total Cooling Capacity
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:DX:SingleSpeed. The calculations for air enthalpy are similar to that done for Coil:Cooling:DX:SingleSpeed. The difference is in calculating the total cooling capacity temperature modifier function at the selected nominal speed level, as below:
TotCapTempModFracNominalSpeed=a+b∗WBi+c∗WB2i+d∗DBo+e∗DBoT2+f∗WBi∗DBo
where
WB_{i} = wetbulb temperature of the air entering thecooling coil, degC
DB_{o} = condenser entering air temperature, degC
af = regression curvefit coefficients.
If the entering air enthalpy is less than the exiting air enthalpy, a reference value of 48,000 J/kg is used as the entering air enthalpy. If the TotCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the rated air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil total cooling capacity is set equal to 0.
If H_{in} > H_{out} Then
˙Qcoil,rated,total=mair,rated(Hin−Hout)/TotCapTempModFracNominalSpeed
Else
˙Qcoil,rated,total=mair,rated(48000−Hout)/TotCapTempModFracNominalSpeed
End If
The other sizing procedures, e.g. evaporative condenser pump, etc., are the same as Coil:Cooling:DX:SingleSpeed.
Coil:Heating:DX:VariableSpeed Sizing[LINK]
For the variablespeed DX heating coil, we specify a nominal speed level. During the sizing calculation, the Rated Heating Capacity at the Selected Nominal Speed Level should be the same as the total cooling capacity of its corresponding cooling coil, which has to be sized first. The default nominal speed level will be the highest speed. However, the model allows the user to select a nominal speed level rather than the highest. If the user chooses to autosize the Rated Air Volume Flow Rate, the flow rate, as compared to the Rated Heating Capacity, is sized to have the same ratio as the air volume flow rate to the heating capacity at the nominal speed, of the Reference Unit. The other sizing procedures are the same as Coil:Heating:DX:SingleSpeed.
The loop pumps’ autosizable inputs are nominal volumetric flow rate and nominal power consumption. We have
Eff_{tot} = Eff_{mot} *Eff_{impeller}
The motor efficiency is an input. Since we need the total efficiency to calculate the nominal power consumption we assume an impeller efficiency of 0,78 for purposes of sizing.
Rated Volumetric Flow Rate[LINK]
This is just set equal to the design loop demand obtained from summing the needs of the components on the demand side of the loop.
Rated Power Consumption[LINK]
˙Qnom=Hnom⋅˙Vnom/Efftot
H_{nom}, the nominal head, is an input.
Electric Chiller Sizing[LINK]
Generally chillers will need nominal cooling capacity, evaporator flow rate and condenser flow rate. All 3 quantities can be straightforwardly obtained using the user specified loop sizing data and the loop design flow rates.
All chillers on a loop are sized to meet the full loop load. If there are multiple chillers on a loop that call for autosizing, they will all be assigned the same cooling capacity and evaporator flow rate.
Nominal Cooling Capacity[LINK]
˙Qchiller,nom=Cp,w⋅ρw⋅ΔTloop,des⋅˙Vlo
Component Sizing[LINK]
Introduction[LINK]
In EnergyPlus each HVAC component sizes itself. Each component module contains a sizing subroutine. When a component is called for the first time in a simulation, it reads in its user specified input data and then calls the sizing subroutine. This routine checks the autosizable input fields for missing data and calculates the data when needed.
A number of highlevel variables are used in the sizing subroutines.
CurDuctType (in DataSizing) contains the information about the current duct type. The types can be main, cooling, heating or other.
CurZoneEqNum (in DataSizing) is the current zone equipment set index and indicates that the component is a piece of zone equipment and should size itself using the zone sizing data arrays.
CurSysNum (in DataSizing) is the current air loop index and indicates that the component is part of the primary air system and should size itself using the system sizing data arrays.
Fan Sizing[LINK]
Fan sizing is done in subroutine SizeFan.
Max Flow Rate[LINK]
If the fan is part of the central air system then check the duct type.
For duct type = main, other or default, ˙Vfan,max=DesMainVolFlowsys
For duct type = cooling, ˙Vfan,max=DesCoolVolFlowsys
For duct type = heating, ˙Vfan,max=DesHeatVolFlowsys
If the fan is zone equipment then check whether it is part of a component that only does heating.
For heating only ˙Vfan,max=DesHeatVolFlowzone ;
Otherwise ˙Vfan,max=max(DesHeatVolFlowzone,DesCoolVolFlowzone)
If the max fan flow rate is less than SmallAirVolFlow the max flow rate is set to zero.
Coil:Cooling:Water[LINK]
The sizing is done in function SizeWaterCoil of module WaterCoils
Initial Calculations[LINK]
For central cooling coils, the first step is to determine the design air flow rate, load, and design air entering and exit conditions. The coil design air flow rate is not generally the same as the maximum system air flow rate (used to size the central fans). The cooling coil peak load (either sensible or total) can occur at a different time than the system peak flow rate. Hence the coil air entering conditions can be different than those at the peak system flow rate. Also, the method of controlling the coil’s cooling output may also affect coil design flow rate as well as the coil design exit temperature and humidity.
By choosing Type of Load to Size On = Sensible or Total in Sizing:System the user indicates to the program to save the cooling coil air flow rate and system air conditions (mixed, return, outside) at the time of either the system cooling sensible or total load peak. Note that the choice VentilationRequirement uses the time of the sensible peak.
Choosing Central Cooling Capacity Control Method = VAV, Bypass, VT, or OnOff indicates which type of cooling output control the program should assume when calculating the design air flow rate. The function GetCoilDesFlowT in module ReportSizingManager calculates the air flow rate and exit air temperature for each capacity control method.Where:
Tsup,avg=Tzones,avg−∑zones˙Qsens,at−peakρaircp,air˙Vcool,air,max
and:
Cp,air : the specific heat of air (J/kgC)
˙mcc,air,peak : the air mass flow rate through the cooling coil at the sensible or total system peak cooling load (m3/s)
∑zones˙Qsens,at−peak : sum of the zone sensible cooling loads at the time of the peak system cooling load
ρair : the density of air (kg/m3)
Tcc,exit : the design cooling coil exit temperature (c)
Tcool,supply : the supply air temperature for cooling specified in Sizing:System (C)
Tmix,at−peak : the mixed air temperature at the time of the system peak cooling load (C)
Tzones,avg : the average zone temperature at the time of the system peak cooling load (C)
˙Vcc,air : the design volumetric air flow rate through the cooling coil (m3/s). This is the flow rate at either the sensible or total cooling load peak from the design period calculations.
˙Vcool,air,max : the maximum cooling volumetric air flow rate from the design calculations (m3/s). This flow rate occurs at the maximum zone cooling demand.
˙Vsys,air,max : the maximum volumetric air flow rate from the design calculations (m3/s). This flow rate occurs at either the maximum zone cooling or heating demand.
Design Coil Load  System Coils[LINK]
Design coil load (cooling capacity) is not an input for Coil:Cooling:Water. It is used for calculating the design water flow rate.
The design load is calculated as:
˙Qcoil,des=˙ma,coil,des(ha,coil,des,in−ha,coil,des,out)
Where:
ha,coil,des,in : is the coil design inlet air enthalpy (J/kg)
ha,coil,des,out : is the coil design outlet air enthalpy (J/kg), and
˙ma,coil,des : is the coil design air mass flow rate (kg/s)
The design air mass flow rate depends on the location of the coil. If the coil is in the outside air stream, the flow rate is set to
ρair˙Va,coil,oa,des where ˙Va,coil.oa,des is the design outside air volumetric flow rate for the system. Otherwise, it is set to:
ρair˙Vcc,air
where ˙Vcc,air is calculated above in the Initial Calculations section.
To obtain the inlet and outlet enthalpies, we need the inlet and outlet temperatures and humidity ratios. The inlet and outlet conditions depend on whether the coil is in the outside air stream and if it is not, whether or not there is outside air preconditioning.
Coil in outside air stream: [LINK]
Tair,in,des=Tout,cool,at−peak (the outside air temperature at the design cooling peak)
Tair,out,des=Tsys,precool (the specified Precool Design Temperature from the System:Sizing object)
Wair,in,des=Wout,cool,at−peak (the outside humidity ratio at the design cooling peak)
Wair,out,des=Wsys,precool (the specified Precool Design Humidity Ratio from the System:Sizing object)
Coil in main air stream, no preconditioning of outside air [LINK]
Tair,in,des=Tmix,cool,at−peak (the mixed air temperature at the design cooling peak)
Wair,in,des=Wmix,cool,at−peak (the mixed humidity ratio at the design cooling peak)
Tair,out,des=Tcc,exit (calculated above in the Initial Calculation section)
Wair,out,des=Wsup,cool (the specified Central Cooling Design Supply Air Humidity Ratio from the Sizing:System object)
Coil in main air stream, outside air preconditioned [LINK]
The oustide air fraction is calculated as (where V_{cc,air} is calculated as above)
foa=˙Vair,out,des˙Vcc,air
Tair,in,des=foaTprecool+(1−foa)Tret,cool,at−peak (Precool temperature is the specified Precool Design Temperature from System:Sizing Manager; T_ret_cool_atpeak is the return temperature at the system cooling peak load)
Wair,in,des=foaWprecool+(1−foa)Wret,cool,at−peak (Precool humidity ratio is the specified Precool Design Humidity Ratio from System:Sizing Manager; W_ret_cool_atpeak is the return humidity ratio at the system cooling peak load)
Tair,out,des=Tcc,exit (calculated above in the Initial Calculation section)
Wair,out,des=Wsup,cool (the specified Central Cooling Design Supply Air Humidity Ratio from the Sizing:System object)
With the inlet and outlet conditions established, we can obtain the inlet and outlet enthalpies:
hair,coil,des,in=PsyHFnTdbW(Tair,in,des,Wair,in,des)hair,coil,des,out=PsyHFnTdbW(Tair,out,des,Wair,out,des)
Where PsyHFnTdbW is the EnergyPlus function for calculation air specific enthalpy given the air temperature and humidity ratio. We now have all we need to calculate the design coil capacity, ˙Qcoil,des .
Design Coil Load  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit or an ZoneHVAC:FourPipeFanCoil, the cooling load (cooling capacity) is passed down from the terminal unit or fan coil sizing calculations. Otherwise the load is defined as:
˙Qcoil,des=˙ma,coil,des(ha,coil,des,in−ha,coil,des,out)
Where:
ha,coil,des,in : is the coil design inlet air enthalpy (J/kg)
ha,coil,des,out : is the coil design outlet air enthalpy (J/kg), and
˙ma,coil,des : is the coil design air mass flow rate (kg/s)
The enthalpies are given by:
hair,coil,des,in=PsyHFnTdbW(Tair,in,des,Wair,in,des)hair,coil,des,out=PsyHFnTdbW(Tair,out,des,Wair,out,des)
Where the inputs to those functions are the coil inlet design conditions. For coils in terminal units these are set at the system level to the system design supply air temperature. For zonal units they are set to design return air, mixed air, or outside air as appropriate to the unit. T_{air,out,des} is set to the zone cooling design supply air temperature as specified in the Zone:Sizing inputs. W_{air,out,des} is set to the zone cooling design supply air humidity ratio as specified in the Zone:Sizing inputs.
Design Water Flow Rate (m^{3}/s)  System Coils[LINK]
The design water volumetric flow rate is calculated using:
˙Vw,coil,des=˙Qcoil,desρwcp,wΔTw,des
Where ΔTw,des is just the Loop Design Temperature Difference user input from Sizing:Plant (if the coil is in the outside air stream, 1/2 the Loop Design Temperature Difference is used). The design coil load Load_{coil,des} is calculated from:
Design Water Flow Rate (m^{3}/s)  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit or an ZoneHVAC:FourPipeFanCoil, the chilled water flow rate is passed down from the terminal unit or fan coil sizing calculations. Otherwise the flow is set to:
˙Vw,coil,des=˙Qcoil,desρwcp,wΔTw,des
Where ΔTw,des is just the Loop Design Temperature Difference user input from Sizing:Plant.
Design Air Flow Rate  System Coils[LINK]
The design air volumetric flow rate for the system cooling coil is set to:
the design outside air flow rate if the coil is in the outside air stream;
the coil design flow rate from function GetCoilDesFlowT described in section “Initial Calculations”;
the design flow rate set by the parent component (such as a unitary system) containing the cooling coil.
Design Air Flow Rate  Zone Coils[LINK]
Zone chilled water coils are always part of a zone HVAC component. In almost all cases the design flow rate is passed down from the design flow rate of the parent component. Otherwise if the parent component does cooling only the flow rate for the coil is set to the zone design cooling flow rate. And if the parent component does both cooling and heating, the coil flow rate is set to the maximum of the zone design cooling and heating flow rates.
Design Air Inlet Temperature  System Coils[LINK]
The inlet air temperature depends on whether the coil is in the outside air stream and if it is not, whether or not there is outside air preconditioning.
Coil in outside air stream: Tair,in,des=Tout,cool,at−peak (the outside air temperature set at the design cooling peak).
Coil in main air stream, no preconditioning of outside air: Tair,in,des=Tmix,cool,at−peak (the mixed air temperature at the cooling design peak).
Coil in main air stream, outside air preconditioned. The outside air fraction is calculated as foa=˙Vair,out,des/˙Vcc,air , where ˙Vcc,air is calculated above. Then Tair,in,des=foaTprecool+(1−foa)Tret,cool,at−peak , where Tprecool is the specified Precool Design Temperature from System:Sizing, and Tret,cool,at−peak is the return temperature at the system cooling peak load.
Design Air Inlet Temperature  Zone Coils[LINK]
The design inlet temperature depends on whether the coil is in a terminal unit or a zonal unit, and where the coil is positioned within the unit.
For the AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction terminal unit the design inlet temperature is set to the zone temperature at the time of the zone cooling peak, since the coil is located in the induced air stream.
For fan coil units the design inlet temperature is set to the mixed air temperature: Tair,in,des=foaToa,coolpeak+(1−foa)Tz,coolpeak , where foa=ρa˙Vz,oa,des/˙mz,cool,des
In all other cases the design inlet temperature is set to the zone design cooling coil inlet temperature which is calculated in the zone sizing simulation and is basically the same calculation as the fan coil unit.
Design Air Outlet Temperature  System Coils[LINK]
The outlet air temperature depends on whether the coil is in the outside air stream.
Coil in outside air stream: T_{air,out,des} = T_{sys,des,precool} (the specified Precool Design Temperature from the Sizing:System object).
Coil in main air stream: the design outlet air temperature is set to the temperature calculated in the Initial Calculation section above.
Design Air Outlet Temperature  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit, then:
˙Qcoil,des=cp,airρair˙Vw,coil,desΔTw,desT1=Tair,in,des−˙Qcoil,des/(ρaircp,air˙Vair,coil,des)T2=Tw,out,des+2Tair,out,des=max(T1,T2)
For all other cases T_{air,out,des} is set to T_{z,sup,des} (the zone design supply air temperature as specified in Sizing:Zone).
Design Inlet Air Humidity Ratio  System Coils[LINK]
The design inlet humidity ratio depends on whether the coil is in the outside air stream and if it is not, whether or not there is outside air preconditioning.
Coil in outside air stream: Wair,in,des=Wout,cool,at−peak (the outside air humidity ratio at the design cooling peak).
Coil in main air stream, no preconditioning of outside air: Wair,in,des=Wmix,cool,at−peak (the mixed air humidity ratio at the cooling design peak).
Coil in main air stream, outside air preconditioned. The outside air fraction is calculated as foa=˙Vair,out,des/˙Vcc,air , where ˙Vcc,air is calculated above. Then Wair,in,des=foaWprecool+(1−foa)Wret,cool,at−peak , where Wprecool is the specified Precool Design Humidity Ratio from System:Sizing, and Wret,cool,at−peak is the return humidity ratio at the system cooling peak load.
Design Air Inlet Humidity Ratio  Zone Coils[LINK]
The design inlet humidity ratio depends on whether the coil is in a terminal unit or a zonal unit, and where the coil is positioned within the unit.
For the AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction terminal unit the design inlet humidity ratio is set to the zone humidity ratio at the time of the zone cooling peak, since the coil is located in the induced air stream.
For fan coil units the design inlet humidity ratio is set to the mixed air humidity ratio: Wair,in,des=foaWoa,coolpeak+(1−foa)Wz,coolpeak , where foa=ρa˙Vz,oa,des/˙mz,cool,des
In all other cases the design inlet humidity ratio is set to the zone design cooling coil inlet hunidity ratio which is calculated in the zone sizing simulation and is basically the same calculation as the fan coil unit.
Design Outlet Air Humidity Ratio  System Coils[LINK]
The outlet air humidity ratio depends on whether the coil is in the outside air stream.
Coil in outside air stream: W_{air,out,des} = W_{sys,des,precool} (the specified Precool Design Humidity Ratio from the Sizing:System object)
Coil in main air stream: W_{air,out,des} = PsyWFnTdbRhPb(T_{air,out,des},0.9,P_{air,std}), where PsyWFnTdbRhPb is the EnergyPlus psychrometric function to calculate humidity ratio from drybulb temperature, relative humidity, and atmospheric pressure. The design outlet humidity ratio is being set to the humidity ratio at 90% relative humidity and design outlet temperature.
Design Outlet Air Humidity Ratio  Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:ConstantVolume:FourPipeInduction unit, then:
Get the dewpoint temperature at W_{air,in,des}: Tdp,in=PsyTdpFnWPb(Wair,in,des,Pair,std)
If T_{dp,in} < = T_{w,in,des} set W_{air,out,des} = W_{air,in,des}. Otherwise set W_{air,out,des} = min(PsyWFnTdbRhPb(T_{air,out,des},0.9,P_{air,std}),W_{air,in,des})
Design Inlet Water Temperature  System Coils[LINK]
The Design Inlet Water Temperature is set to the Design Loop Exit Temperature specified in the Sizing:Plant object for the water loop serving this coil.
Design Inlet Water Temperature  Zone Coils[LINK]
The Design Inlet Water Temperature is set to the Design Loop Exit Temperature specified in the Sizing:Plant object for the water loop serving this coil.
Coil:Cooling:Water:DetailedGeometry Sizing[LINK]
The sizing is done in subroutine SizeWaterCoil
Max Water Flow Rate of Coil[LINK]
The calculation is identical to that done for Coil:Cooling:Water.
Number of Tubes per Row[LINK]
Ntube/row=Int(13750⋅˙Vcoil,water,max)
N_{tube/row} = Max(N_{tube/row},3)
Fin Diameter[LINK]
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙mair,des=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙mair,des=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙mair,des=ρair⋅DesHeatVolFlowsys
Dfin=0.335⋅˙mair,des
Minimum Air Flow Area[LINK]
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙mair,des=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙mair,des=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙mair,des=ρair⋅DesHeatVolFlowsys
AMinAirFlow=0.44⋅˙mair,des
Fin Surface Area[LINK]
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙mair,des=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙mair,des=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙mair,des=ρair⋅DesHeatVolFlowsys
AFinSurf=78.5⋅˙mair,des
Total Tube Inside Area[LINK]
A_{tube,total\ inside} = 4.4 *D_{tube,inside} *N_{tube\ rows} *N_{tubes/row}
Where D_{tube,inside} is the tube inside diameter.
Tube Outside Surf Area[LINK]
A_{tube,outside} = 4.1 *D_{tube,outside} *N_{tube\ rows} *N_{tubes/row}
Where D_{tube,outside} is the tube outside diameter.
Coil Depth[LINK]
Depth_{coil} = Depth_{tube\ spacing} * N_{tube\ rows}
Coil:Cooling:WaterToAirHeatPump:EquationFit Sizing[LINK]
The sizing is done in subroutine SizeHVACWaterToAir
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water, which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. If there is a companion heating coil, the heating coil design load is used so that both modes will have the same rated water flow rate. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Cooling Capacity[LINK]
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:Water. The following calculations are then performed to determine the rated total cooling capacity.
TWB,ratio={\left( {{T_{WB,air,in,des}} + 273.15\,C} \right)}/(TWB,air,in,des+273.15C)283.15C{283.15\,C}
TS,ratio={\left( {29.44\,C + 273.15\,C} \right)}/(29.44C+273.15C)283.15C{283.15\,C}
where:
${T_{WB,ratio}} = $ ratio of loadside inlet air wetbulb temperature in Kelvin to a reference temperature
${T_{S,ratio}} = $ ratio of sourceside inlet water temperature in Kelvin to a reference temperature
TCC1 = user input for Total Cooling Capacity Coefficient 1
TCC2 = user input for Total Cooling Capacity Coefficient 2
TCC3 = user input for Total Cooling Capacity Coefficient 3
TCC4 = user input for Total Cooling Capacity Coefficient 4
TCC5 = user input for Total Cooling Capacity Coefficient 5
TotCapTempModFac=TCC1+TCC2(TWB,ratio)+TCC3(TS,ratio)+TCC4+TCC5
The 4^{th} and 5^{th} coefficient (TCC4 and TCC5) used in the above equation are multipliers for the loadside and sourceside flow ratios, respectively. For sizing, these ratios are assumed to be 1.
The enthalpy of the entering air is then compared with the enthalpy of the exiting air. The calculations for air enthalpy are identical to that done for Coil:Cooling:Water. If the entering air enthalpy is less than the exiting air enthalpy, a reference value of 48,000 J/kg is used as the entering air enthalpy. If the TotCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the design air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil total cooling capacity is set equal to 0.
IF(Hin>Hout)THEN⋅Qcoil,des,total={{m_{air,des}}\left( {{H_{in}}  {H_{out}}} \right)}/mair,des(Hin−Hout)TotCapTempModFac{TotCapTempModFac}ELSE⋅Qcoil,des,total={{m_{air,des}}\left( {48000  {H_{out}}} \right)}/mair,des(48000−Hout)TotCapTempModFac{TotCapTempModFac}ENDIF
Rated Sensible Cooling Capacity[LINK]
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:Water. The following calculations are then performed to determine the rated sensible cooling capacity.
TDB,ratio={\left( {{T_{DB,air,in,des}} + 273.15\,C} \right)}/(TDB,air,in,des+273.15C)283.15C{283.15\,C}
TS,ratio={\left( {29.44\,C + 273.15\,C} \right)}/(29.44C+273.15C)283.15C{283.15\,C}
where:
${T_{DB,ratio}} = $ ratio of loadside inlet air drybulb temperature in Kelvin to a reference temperature
SCC1 = user input for Sensible Cooling Capacity Coefficient 1
SCC2 = user input for Sensible Cooling Capacity Coefficient 2
SCC3 = user input for Sensible Cooling Capacity Coefficient 3
SCC4 = user input for Sensible Cooling Capacity Coefficient 4
SCC5 = user input for Sensible Cooling Capacity Coefficient 5
SCC6 = user input for Sensible Cooling Capacity Coefficient 6
SensCapTempModFac=SCC1+SCC2(TDB,ratio)+SCC3(TWB,ratio)+SCC4(TS,ratio)+SCC5+SCC6
The 5^{th} and 6^{th} coefficient (SCC5 and SCC6) used in the above equation are multipliers for the loadside and sourceside flow ratios, respectively. For sizing, these ratios are assumed to be 1.
The drybulb temperature of the entering air is then compared with the drybulb temperature of the exiting air. The calculations for air drybulb temperature are identical to that done for Coil:Cooling:Water. If the entering air drybulb temperature is less than the exiting air drybulb temperature, a reference value of 24 C is used as the entering air drybulb temperature. If the SensCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the design air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil sensible cooling capacity is set equal to 0.
IF(TDB,in>TDB,out)THEN⋅Qcoil,des,sensible={{m_{air,des}}C{p_{air,des}}\left( {{T_{DB,in}}  {T_{DB,out}}} \right)}/mair,desCpair,des(TDB,in−TDB,out)SensCapTempModFac{SensCapTempModFac}ELSE⋅Qcoil,des,sensible={{m_{air,des}}C{p_{air,des}}\left( {24  {T_{DB,out}}} \right)}/mair,desCpair,des(24−TDB,out)SensCapTempModFac{SensCapTempModFac}ENDIF
Coil:Cooling:WaterToAirHeatPump:VariableSpeedEquationFit Sizing[LINK]
For the cooling coil of VS WSHP, we specify a nominal speed level. During the sizing calculation, the Rated Air Volume Flow Rate, the Rated Water Volume Flow Rate and the Rated Total Cooling Capacity at the Selected Nominal Speed Level are determined in the same way as the Coil:Cooling:WaterToAirHeatPump:EquationFit object. The sensible heat transfer rate is not allowed for autosizing, instead, it is a function of the rated air and water flow rates, rated total cooling capacity and the Reference Unit SHR at the nominal speed level. The default nominal speed level is the highest speed. However, the model allows the user to select a nominal speed level rather than the highest.
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit , which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. If there is a companion heating coil, the heating coil design load is used so that both modes will have the same rated water flow rate. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Cooling Capacity[LINK]
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit. The calculations for air enthalpy are similar to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit. The difference is in calculating the total cooling capacity temperature modifier function at the selected nominal speed level, as below:
TotCapTempModFracNominalSpeed=a+b∗WBi+c∗WB2i+d∗EWT+e∗EWT2+f∗WBi∗EWT
where
WB_{i} = wetbulb temperature of the air entering the heating coil, degC
EWT = entering water temperature, degC
af = regression curvefit coefficients.
If the entering air enthalpy is less than the exiting air enthalpy, a reference value of 48,000 J/kg is used as the entering air enthalpy. If the TotCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the rated air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil total cooling capacity is set equal to 0.
If H_{in} > H_{out} Then
˙Qcoil,rated,total=mair,rated(Hin−Hout)/TotCapTempModFracNominalSpeed
Else
˙Qcoil,rated,total=mair,rated(48000−Hout)/TotCapTempModFracNominalSpeed
End If
Coil:Heating:WaterToAirHeatPump:EquationFit Sizing[LINK]
The sizing is done in subroutine SizeHVACWaterToAir.
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:Water , which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Heating Capacity[LINK]
The rated total heating capacity is set equal to the rated total cooling capacity.
Coil:Heating:WaterToAirHeatPump:VariableSpeedEquationFit Sizing[LINK]
For the heating coil of VS WSHP, we specify a nominal speed level. During the sizing calculation, the Rated Air Volume Flow Rate and the Rated Water Volume Flow Rate are determined in the same way as the Coil:Heating:WaterToAirHeatPump:EquationFit object. On the other hand, the Rated Heating Capacity at the Selected Nominal Speed Level should be the same as the total cooling capacity of its corresponding cooling coil, which has to be sized first. The default nominal speed level will be the highest speed. However, the model allows the user to select a nominal speed level rather than the highest.
Rated Air Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit.
Rated Water Flow Rate[LINK]
The calculation is identical to that done for Coil:Cooling:WaterToAirHeatPump:EquationFit, which is the coil design load divided by the Loop Design Temperature Difference user input from Sizing:Plant. For sizing the plant loop serving this coil, only one half of this flow rate is used since both the cooling and heating coil will save a flow rate but only one of these coils will operate at a time.
Rated Total Heating Capacity[LINK]
The rated total heating capacity is set equal to the rated total cooling capacity.
Coil:Heating:Water Sizing[LINK]
The sizing is done in subroutine SizeWaterCoil.
Max Water Flow Rate of Coil[LINK]
System Coils[LINK]
With the coil load from the system design data array and the user specified (in a Sizing:Plant object) design hot water temperature fall, calculate the max water flow rate:
˙Vcoil,water,max=HeatCapsys/(Cp,water⋅ρwater⋅ΔTplt,hw,des)
Zone Coils[LINK]
Using the zone design coil inlet and supply air conditions calculate the design coil load.
If the coil is not part of an induction unit then obtain the coil inlet temperature from the zone design data array;
T_{in,air} = DesHeatCoilInTemp_{zone}
If the coil is part of an induction unit take into account the induced air:
Frac_{minflow} = MinFlowFrac_{zone}
T_{in,air} = DesHeatCoilInTemp_{zone} * Frac_{minflow} +
ZoneTempAtHeatPeak_{ zone} *(1 Frac_{minflow})
T_{out,air} = HeatDesTemp_{zone}
W_{out,air} = HeatDesHumRat_{zone}
If the coil is part of a terminal unit the mass flow rate is determined by the volumetric flow rate of the terminal unit:
˙mair,des=ρair⋅˙mair,des,tu
Otherwise the design flow is obtained from the zone design data array:
˙mair,des=DesHeatMassFlowzone
Qcoil,des=cp,air˙mair,des⋅(Tout,air−Tin,air)
Here c_{p,air} is calculated at the outlet humidity and the average of the inlet and outlet temperatures.
With the coil load and the user specified (in a Sizing:Plant object) design hot water temperature decrease, calculate the max water flow rate:
˙Vcoil,water,max=Qcoil,des/(Cp,water⋅ρwater⋅ΔTplt,hw,des)
UA of the Coil[LINK]
To obtain the UA of the coil, we specify the model inputs (other than the UA) at design conditions and the design coil load that the coil must meet. Then we numerically invert the coil model to solve for the UA that will enable the coil to meet the design coil load given the specified inputs.
System Coils[LINK]
The design coil load is the system design sensible cooling capacity;
Q_{coil,des} = HeatCap_{sys}
The required inputs for the simple coil model are:
T_{in,air} = HeatMixTemp_{sys}
W_{in,air} = HeatMixHumRat_{sys}
T_{in,water} = ExitTemp_{plt,hw,des}
˙min,water=ρwater⋅˙Vcoil,water,max
Depending on the duct type, get the coil design air flow rate.
For duct type = main, other or default
˙min,air=ρair⋅DesMainVolFlowsys
for duct type = cooling
˙min,air=ρair⋅DesCoolVolFlowsys
for duct type = heating
˙min,air=ρair⋅DesHeatVolFlowsys
We now have all the data needed to obtain UA. The numerical inversion is carried out by calling subroutine SolveRegulaFalsi. This is a general utility routine for finding the zero of a function. In this case it finds the UA that will zero the residual function  the difference between the design coil load and the coil output divided by the design coil load. The residual is calculated in the function SimpleHeatingCoilUAResidual.
Zone Coils[LINK]
If the coil is not part of an induction unit then obtain the coil inlet temperature from the zone design data array;
T_{in,air} = DesHeatCoilInTemp_{zone}
If the coil is part of an induction unit take into account the induced air:
Frac_{minflow} = MinFlowFrac_{zone}
T_{in,air} = DesHeatCoilInTemp_{zone} * Frac_{minflow} +
ZoneTempAtHeatPeak_{ zone} *(1 Frac_{minflow})
W_{in,air} = DesHeatCoilInHumRat_{zone}
T_{in,water} = ExitTemp_{plt,hw,des}
˙min,water=ρwater⋅˙Vcoil,water,max
T_{out,air} = HeatDesTemp_{zone}
W_{out,air} = HeatDesHumRat_{zone}
If the coil is part of a terminal unit the mass flow rate is determined by the volumetric flow rate of the terminal unit:
˙mair,des=ρair⋅˙mair,des,tu
Otherwise the design flow is obtained from the zone design data array:
˙mair,des=DesHeatMassFlowzone
˙Qcoil,des=cp,air⋅˙mair,des⋅(Tout,air−Tin,air)
Here c_{p,air} is calculated at the outlet humidity and the average of the inlet and outlet temperatures.
We now have all the data needed to obtain UA. The numerical inversion is carried out by calling subroutine SolveRegulaFalsi. This is a general utility routine for finding the zero of a function. In this case it finds the UA that will zero the residual function  the difference between the design coil load and the coil output divided by the design coil load. The residual is calculated in the function SimpleHeatingCoilUAResidual.
Coil:Heating:Steam Sizing[LINK]
The sizing is done in subroutine SizeSteamCoil.
Maximum Steam Flow Rate[LINK]
System Coils[LINK]
The maximum steam volumetric flow rate is calculated using:
˙Vcoil,steam,max=Loadcoil,desρsteam(hfg+cp,w⋅ΔTsc)
The steam density (ρsteam ) is for saturated steam at 100 degC (101325.0 Pa) and h_{fg} is the latent heat of vaporization of water at 100 degC (101325.0 Pa). C_{p,w} is the heat capacity of saturated water (condensate) at 100 degC (101325.0 Pa) and ΔTsc is the Degree of Subcooling defined in the Coil:Heating:Steam object input. The design coil load Load_{coil,des} is calculated from:
Loadcoil,des=˙mair,des(cp,air)(Tair,coil,des,out−Tair,coil,des,in)
The design air mass flow rate depends on the location of the coil (duct type). For duct type = main, the flow rate is set to r_{air} *DesMainVolFlow_{sys} *MinSysAirFlowRatio. If the coil is in a cooling duct the flow rate is set to r_{air} *DesCoolVolFlow_{sys} *MinSysAirFlowRatio. If the coil is in a heating duct the flow rate is set to r_{air} *DesHeatVolFlow_{sys}. If the coil is in any other kind of duct, the flow rate is set to r_{air} *DesMainVolFlow_{sys}.
For sizing, the design outlet air temperature (T_{air,coil,des,out}) is the Central Heating Design Supply Air Temperature specified in the Sizing:System object.
The design inlet air temperature depends on whether the coil is being sized for 100% outdoor air or minimum outdoor air flow (per 100% Outdoor Air in Heating input field in the Sizing:System object).
T_{air,coil,des,in} = HeatOutTemp_{sys} (the outdoor air temperature at the design heating peak)
T_{air,coil,des,in} = Frac_{oa}* * HeatOutTemp_{sys}* + (1. Frac_{oa}) *HeatRetTemp_{sys} (see Table 41. System Sizing Data)
Zone Coils[LINK]
If the coil is part of an AirTerminal:SingleDuct:* unit (e.g., AirTerminal:SingleDuct:ConstantVolume:Reheat, AirTerminal:SingleDuct:VAV:Reheat, AirTerminal:SingleDuct:SeriesPIU:Reheat, etc.), the maximum steam flow rate is set equal to the terminal unit’s maximum steam flow rate. Otherwise (e.g., the zonelevel coil is part of ZoneHVAC:PackagedTerminalAirConditioner, ZoneHVAC:UnitVentilator, ZoneHVAC:UnitHeater or ZoneHVAC:VentilatedSlab) the calculation is similar to that at the system level. A design load is calculated:
Loadcoil,des=˙mair,des(cp,air)(Tair,coil,des,out−Tair,coil,des,in)
where:
˙mair,des = DesHeatMassFlow_{zone} (see Table 40. Zone Sizing Data)
T_{air,coil,des,in} = DesHeatCoilInTemp_{zone} (see Table 40)
T_{air,coil,des,out} = HeatDesTemp_{zone} (user input from Sizing:Zone object)
cp,air = Specific heat of air (evaluated at the average of inlet and outlet air temperatures, and at the zone heating design supply air humidity ratio HeatDesHumRat_{zone} [user input from Sizing:Zone object])
˙Vcoil,steam,max=Loadcoil,desρsteam(hfg+cp,w⋅ΔTsc)
The terms in the denominator of this equation (ρ_{steam}, h_{fg}, etc.) are evaluated in the same way as described above for steam System Coils.
Sizing of Gas and Electric Heating Coils[LINK]
The sizing calculation is done in subroutine SizeHeatingCoil in module HeatingCoils.
Nominal Capacity of the Coil[LINK]
System Coils[LINK]
The value is obtained from the system design array.
Cap_{nom} = HeatCap_{sys}
Zone Coils[LINK]
The capacity is calculated from the design coil inlet and outlet conditions.
If the coil is not part of an induction unit then obtain the coil inlet temperature from the zone design data array;
T_{in,air} = DesHeatCoilInTemp_{zone}
If the coil is part of an induction unit take into account the induced air:
Frac_{minflow} = MinFlowFrac_{zone}
T_{in,air} = DesHeatCoilInTemp_{zone} * Frac_{minflow} +
ZoneTempAtHeatPeak_{ zone} *(1 Frac_{minflow})
T_{out,air} = HeatDesTemp_{zone}
W_{out,air} = HeatDesHumRat_{zone}
Q_{coil,des} = C_{p,air} * DesHeatMassFlow_{zone} *(T_{out,air}T_{in,air})
Here c_{p,air} is calculated at the outlet humidity and the average of the inlet and outlet temperatures.
DX Coil Sizing[LINK]
The sizing calculations are done in subroutine SizeDXCoil in module DXCoils. This section covers the sizing of the objects
Coil:Cooling:DX:SingleSpeed
Coil:Heating:DX:SingleSpeed
Coil:Cooling:DX:TwoSpeed
Rated Air Volume Flow Rate[LINK]
System Coils[LINK]
The rated air flow rate is obtained from the system design array.
˙Vair,rated=DesMainVolFlowsys
Zone Coils[LINK]
The rated air flow rate is the maximum of the heating and cooling design flow rates from the zone design array.
˙Vair,rated=Max(DesCoolVolFlowzone,DesHeatVolFlowzone)
Rated Total Cooling Capacity[LINK]
System Coils[LINK]
The rated cooling capacity is obtained by dividing the peak cooling capacity by the Cooling Capacity Modifier Curve evaluated at peak mixed wetbulb and outdoor drybulb temperatures.
T_{mix} = CoolMixTemp_{sys}
W_{mix} = CoolMixHumRat_{sys}
T_{sup} = CoolSupTemp_{sys}
W_{sup} = CoolSupHumRat_{sys}
T_{outside} = CoolOutTemp_{sys}
r_{air} = PsyRhoAirFnPbTdbW(p_{air,std}, T_{mix},W_{mix})
h_{mix} = PsyHFnTdbW(T_{mix},W_{mix})
h_{sup} = PsyHFnTdbW(T_{sup},W_{sup})
T_{mix,wb} = PsyTwbFnTdbWPb(T_{mix},W_{mix}, p_{air,std})
CapModFac = CurveValue(CCapFTemp,T_{mix,wb},T_{outside})
CCappeak=ρair⋅˙Vair,rated⋅(hmix−hsup)
CCap_{rated} = CCap_{peak}/ CapModFac
We check that the design volume flow per total capacity is within the prescribed range:
FlowCapRatio=˙Vair,rated/CCaprated
If FlowCapRatio < FlowCapRatio_{min} then
CCaprated=˙Vair,rated/FlowCapRatiomin
If FlowCapRatio > FlowCapRatio_{max} then
CCaprated=˙Vair,rated/FlowCapRatiomax
where
FlowCapRatio_{min} = 0.00004027 m^{3}/s per watt (300 cfm/ton)
And
FlowCapRatio_{max} = 0.00006041 m^{3}/s per watt (450 cfm/ton)
The sizing calculation for DX cooling coils for 100% dedicated outdor air system (DOAS) are identical to regular DX cooling coils. However, they operate operate at different flow to capacity ratio ranges and are within the prescribed range below:
FlowCapRatio_{min} = 0.00001677 m^{3}/s per Watt (125 cfm/ton)
And
FlowCapRatio_{max} = 0.00003355 m^{3}/s per Watt (250 cfm/ton)
Zone Coils[LINK]
The rated cooling capacity for zone coils is calculated in the same manner as for system coils.
T_{mix} = DesCoolCoilInTemp_{zone}
W_{mix} = DesCoolCoilInHumRat_{zone}
T_{sup} = CoolDesTemp_{zone}
W_{sup} = CoolDesHumRat_{zone}
T_{outside} = T_{outside},_{desday,peak}
r_{air} = PsyRhoAirFnPbTdbW(p_{air,std}, T_{mix},W_{mix})
h_{mix} = PsyHFnTdbW(T_{mix},W_{mix})
h_{sup} = PsyHFnTdbW(T_{sup},W_{sup})
T_{mix,wb} = PsyTwbFnTdbWPb(T_{mix},W_{mix}, p_{air,std})
CapModFac = CurveValue(CCapFTemp,T_{mix,wb},T_{outside})
CCappeak=ρair⋅˙Vair,rated⋅(hmix−hsup)
CCap_{rated} = CCap_{peak}/ CapModFac
We check that the design volume flow per total capacity is within the prescribed range:
FlowCapRatio=˙Vair,rated/CCaprated
If FlowCapRatio < FlowCapRatio_{min} then
CCaprated=˙Vair,rated/FlowCapRatiomin
If FlowCapRatio > FlowCapRatio_{max} then
CCaprated=˙Vair,rated/FlowCapRatiomax
where
FlowCapRatio_{min} = 0.00004027 m^{3}/s per watt (300 cfm/ton)
And
FlowCapRatio_{max} = 0.00006041 m^{3}/s per watt (450 cfm/ton)
We check the design flow to the total cooling capacity rato for dedicated zone outdoor unit DX cooling coils to be within the limits prescribed below:
FlowCapRatio_{min} = 0.00001677 m^{3}/s per Watt (125 cfm/ton)
And
FlowCapRatio_{max} = 0.00003355 m^{3}/s per Watt (250 cfm/ton)
Rated Total Heating Capacity[LINK]
For Coil:Heating:DX:SingleSpeed the rated heating capacity is set equal to the cooling capacity.
Rated SHR[LINK]
The rated sensible heat ratio is calculated to be the sensible cooling (from rated inlet conditions to user specified supply conditions) divided by the total cooling (from rated inlet to specified supply).
T_{in,rated} = 26.6667 ^{o}C (80 ^{o}F)
W_{in,rated} = 0.01125 (corresponds to 80 ^{o}F drybulb, 67 ^{o}F wetbulb)
C_{p,air} = PsyCpAirFnWTdb(W_{in,rated}, T_{in,rated})
For system coils
T_{sup} = CoolSupTemp_{sys}
W_{sup} = CoolSupHumRat_{sys}
For zone coils
T_{sup} = CoolDesTemp_{zone}
W_{sup} = CoolDesHumRat_{zone}
Then
h_{rated} = PsyHFnTdbW(T_{in,rated}, W_{in,rated})
h_{sup} = PsyHFnTdbW(T_{sup}, W_{sup})
Dh_{rated,sup} = h_{rated}h_{sup}
DQs_{rated,sup} = C_{p,air} *(T_{in,rated}T_{sup})
SHR_{rated} = DQs_{rated,sup}/Dh_{rated,sup}
Evaporative Condenser Air Volume Flow Rate[LINK]
The evaporative condenser air volume flow rate (m^{3}/s) is set to 0.000114 m^{3}/s per watt (850 cfm/ton) times the total rated cooling capacity.
Evaporative Condenser Air Volume Flow Rate, Low Speed[LINK]
The evaporative condenser air volume flow rate, low speed (m^{3}/s) is set to 1/3 times 0.000114 m^{3}/s per watt (850 cfm/ton) times the total rated cooling capacity.
Evaporative Condenser Pump Rated Power Consumption[LINK]
The evaporative condenser pump rated power consumption is set equal to the total cooling capacity times 0.004266 watts pump power per watt capacity (15 W/ton).
Evaporative Condenser Pump Rated Power Consumption, Low Speed[LINK]
The evaporative condenser pump rated power consumption, low speed, is set equal to 1/3 times the total cooling capacity times 0.004266 watts pump power per watt capacity (15 W/ton).
Rated Air Volume Flow Rate, low speed[LINK]
The rated air volume flow rate, low speed, is set equal to 1/3 times the full rated air volume flow rate.
Rated Total Cooling Capacity, Low Speed[LINK]
The rated total cooling capacity, low speed, is set equal to 1/3 times the full rated total cooling capacity.
Rated SHR, low speed[LINK]
The rated sensible heat ratio, low speed, is set equal to the full speed SHR.
Resistive Defrost Heater Capacity[LINK]
For the heat pump the resistive defrost heat capacity is set equal to the cooling capacity.
DX MultiSpeed Coil Sizing[LINK]
The sizing calculations are done in subroutine SizeDXCoil in module DXCoils. This section covers the sizing of the objects
Coil:Heating:DX:MultiSpeed
Coil:Cooling:DX: MultiSpeed
The rated air volume flow rate, rated total cooling capacity, rated heating capacity, rated SHR, evaporative condenser air volume flow rate, evaporative condenser pump rated power consumption at the highest speed are sized in the same ways as DX Coil Sizing.
After the sizes are determined at the highest speed, the sizes in the rest of speeds are assumed to
Valuen=nNumberOfSpeed∗ValueNumberOfSpeed
where
Value_{n} = Any autosizable variable at Speed n, except SHR
SHR_{n} = SHR_{NumberOfSpeed}
n = Speed Index number from 1 to NumberOfSpeed1
NumberOfSpeed = The highest speed number
Coil:Cooling:DX:VariableSpeed Sizing[LINK]
For the variablespeed DX cooling coil, we specify a nominal speed level. During the sizing calculation, the Rated Total Cooling Capacity at the Selected Nominal Speed Level is determined in the same way as the Coil:Cooling:DX:SingleSpeed object. If the user chooses to autosize the Rated Air Volume Flow Rate, the flow rate, as compared to the Rated Total Cooling Capacity, is sized to have the same ratio as the air volume flow rate to the total cooling capacity at the nominal speed, of the Reference Unit. The sensible heat transfer rate is not allowed for autosizing, instead, it is a function of the rated air flow, rated total cooling capacity and the Reference Unit SHR at the nominal speed level. The default nominal speed level is the highest speed. However, the model allows the user to select a nominal speed level rather than the highest.
Rated Total Cooling Capacity
The calculation for coil operating temperatures (inlet and outlet) are identical to that done for Coil:Cooling:DX:SingleSpeed. The calculations for air enthalpy are similar to that done for Coil:Cooling:DX:SingleSpeed. The difference is in calculating the total cooling capacity temperature modifier function at the selected nominal speed level, as below:
TotCapTempModFracNominalSpeed=a+b∗WBi+c∗WB2i+d∗DBo+e∗DBoT2+f∗WBi∗DBo
where
WB_{i} = wetbulb temperature of the air entering thecooling coil, degC
DB_{o} = condenser entering air temperature, degC
af = regression curvefit coefficients.
If the entering air enthalpy is less than the exiting air enthalpy, a reference value of 48,000 J/kg is used as the entering air enthalpy. If the TotCapTempModFac calculation above yields 0 as the result, a value of 1 is used in the following calculation. If the rated air mass flow rate is determined to be less than a very small flow value (0.001 kg/s) or the capacity calculated here is less than 0, the coil total cooling capacity is set equal to 0.
If H_{in} > H_{out} Then
˙Qcoil,rated,total=mair,rated(Hin−Hout)/TotCapTempModFracNominalSpeed
Else
˙Qcoil,rated,total=mair,rated(48000−Hout)/TotCapTempModFracNominalSpeed
End If
The other sizing procedures, e.g. evaporative condenser pump, etc., are the same as Coil:Cooling:DX:SingleSpeed.
Coil:Heating:DX:VariableSpeed Sizing[LINK]
For the variablespeed DX heating coil, we specify a nominal speed level. During the sizing calculation, the Rated Heating Capacity at the Selected Nominal Speed Level should be the same as the total cooling capacity of its corresponding cooling coil, which has to be sized first. The default nominal speed level will be the highest speed. However, the model allows the user to select a nominal speed level rather than the highest. If the user chooses to autosize the Rated Air Volume Flow Rate, the flow rate, as compared to the Rated Heating Capacity, is sized to have the same ratio as the air volume flow rate to the heating capacity at the nominal speed, of the Reference Unit. The other sizing procedures are the same as Coil:Heating:DX:SingleSpeed.
Pump Sizing[LINK]
The loop pumps’ autosizable inputs are nominal volumetric flow rate and nominal power consumption. We have
Eff_{tot} = Eff_{mot} *Eff_{impeller}
The motor efficiency is an input. Since we need the total efficiency to calculate the nominal power consumption we assume an impeller efficiency of 0,78 for purposes of sizing.
Rated Volumetric Flow Rate[LINK]
This is just set equal to the design loop demand obtained from summing the needs of the components on the demand side of the loop.
Rated Power Consumption[LINK]
˙Qnom=Hnom⋅˙Vnom/Efftot
H_{nom}, the nominal head, is an input.
Electric Chiller Sizing[LINK]
Generally chillers will need nominal cooling capacity, evaporator flow rate and condenser flow rate. All 3 quantities can be straightforwardly obtained using the user specified loop sizing data and the loop design flow rates.
All chillers on a loop are sized to meet the full loop load. If there are multiple chillers on a loop that call for autosizing, they will all be assigned the same cooling capacity and evaporator flow rate.
Nominal Cooling Capacity[LINK]
˙Qchiller,nom=Cp,w⋅ρw⋅ΔTloop,des⋅˙Vlo