Engineering Reference — EnergyPlus 8.7

<< Prev | Table of Contents | Next >>

System Design Loads and Air Flow Rates[LINK]

Overview[LINK]

The purpose of the system design calculation is to estimate design heating and cooling loads and air flow rates for each air loop in the simulation problem. The calculation sequence for system level design loads and air flow rates resembles the calculation sequence for zone loads and air flow rates. There is an update subroutine UpdateSysSizing called at the beginning, during, and end of a loop in the Sizing Manager over all the design days. The major difference is that this calculation is done at the zone time-step only. There is no idealized component calculation triggered at the system time-step as in the zone calculation. The system design calculation operates at the zone time step using the design environment weather data and the data stored in the zone sizing arrays. The results of the system design calculation are stored in the system sizing arrays described below.

System Design Data Arrays[LINK]

The system design data arrays are:

SysSizInput(i) stores the input data from the Sizing:System objects.

SysSizing(i,j) stores the results of the system design calculations for all systems and all design days. The index i is for air loops, j for design days.

CalcSysSizing(i) stores the results of the system design calculations for the peak heating and cooling cases for each air loop. The index i is for the air loops.

FinalSysSizing(i) corresponds to CalcSysSizing but includes the effect of the user specified sizing factor or user specified system design flow rate.

The data stored in SysSizing, CalcSysSizing and FinalSysSizing includes the following data items.

System Sizing Data
Name Description
All the data from SysSizInput
CoinCoolMassFlow coincident peak cooling mass flow rate [kg/s]
CoinHeatMassFlow coincident peak heating mass flow rate [kg/s]
NonCoinCoolMassFlow noncoincident peak cooling mass flow rate [kg/s]
NonCoinHeatMassFlow noncoincident peak heating mass flow rate [kg/s]
DesMainVolFlow design main supply duct volume flow [m /s]
DesHeatVolFlow design heat supply duct volume flow [m /s]
DesCoolVolFlow design cool supply duct volume flow [m /s]
SensCoolCap design sensible cooling capacity [W]
TotCoolCap design total cooling capacity [W]
HeatCap design heating capacity [W]
PreheatCap design preheat capacity [W]
MixTempAtCoolPeak mixed air temperature at the time of the cooling peak [C]
MixHumRatAtCoolPeak mixed air humidity ratio at the time of the cooling peak [kg water/kg dry air]
RetTempAtCoolPeak return air temperature at the time of the cooling peak [C]
RetHumRatAtCoolPeak return air humidity ratio at the time of the cooling peak [kg water/kg dry air]
OutTempAtCoolPeak outside air temperature at the time of the cooling peak [C]
OutHumRatAtCoolPeak outside air humidity ratio at the time of the cooling peak [kg water/kg dry air]
HeatMixTemp design mixed air temperature for heating [C]
HeatMixHumRat design mixed air humidity ratio for heating [kg water/kg dry air]
HeatRetTemp design return air temperature for heating [C]
HeatRetHumRat design return air humidity ratio for heating [kg water/kg dry air]
HeatOutTemp design outside air temperature for heating [C]
HeatOutHumRat design outside air humidity ratio for heating [kg water/kg dry air]
HeatFlowSeq(i) daily sequence of system heating air mass flow rate (zone time step) [kg/s]
CoolFlowSeq(i) daily sequence of system cooling air mass flow rate (zone time step) [kg/s]
SumZoneCoolLoadSeq(I) daily sequence of zones’ summed cooling load (zone time step) [W]
CoolZoneAvgTempSeq(I) daily sequence of zones flow weighted average temperature (zone time step) [W]
SensCoolCapSeq(I) daily sequence of system sensible cooling capacity (zone time step) [W]
TotCoolCapSeq(I) daily sequence of system total cooling capacity (zone time step) [W]
HeatCapSeq(i) daily sequence of system heating capacity (zone time step) [W]
PreHeatCapSeq(i) daily sequence of system preheat capacity (zone time step) [W]
SysCoolRetTempSeq(i) daily sequence of system cooling return temperatures (zone time step) [C]
SysCoolRetHumRatSeq(I) daily sequence of system cooling return humidity ratios (zone time step) [kg water/kg dry air]
SysHeatRetTempSeq(i) daily sequence of system heating return temperatures (zone time step) [C]
SysHeatRetHumRatSeq(I) daily sequence of system heating return humidity ratios (zone time step) [kg water/kg dry air]
SysCoolOutTempSeq daily sequence of system cooling outside temperatures (zone time step) [C]
SysCoolOutHumRatSeq daily sequence of system cooling outside humidity ratios (zone time step) [kg water/kg dry air]
SysHeatOutTempSeq daily sequence of system heating outside temperatures (zone time step) [C]
SysHeatOutHumRatSeq daily sequence of system heating outside humidity ratios (zone time step) [kg water/kg dry air]

System Design Flow Rate and Load Summation and Adjustment[LINK]

There is no system level subroutine corresponding to SizeZoneEquipment. Instead the system design loads and flow rates are calculated using the zone level results. The zone design flow rates for the zones served by an air loop are summed to obtain the system level design flow rates. These air flows are mixed with the system level design minimum outside air flow rate to obtain system design coil loads. These activities are all performed within the UpdateSysSizing subroutine in the SimAirServingZones module. It is called at the start of each design day (CallIndicator = BeginDay), at the zone time-step (CallIndicator = DuringDay), at the end of the design day (CallIndicator = EndDay) and at the end of the zone design calculation (CallIndicator = EndSysSizingCalc).

There is a logical flag SysSizingCalc corresponding to ZoneSizingCalc. It is used to allow the component routines to distinguish a normal simulation call from a being called during a system sizing calculation.

BeginDay[LINK]

  1. The environment (in this case, a design day) name is stored in the system sizing data structures.

  2. Loop over the zones cooled by this air loop:

  1. NonCoinCoolMassFlowsys = DesCoolMassFlowzone
  1. Loop over the zones heated by this air loop:
  1. NonCoinHeatMassFlowsys = DesHeatMassFlowzone

DuringDay[LINK]

  1. Loop over the zones cooled by this air loop:

CoolFlowSeqsys(i) = CoolFlowSeq zone(i)

SysCoolRetTemp(i) = (CoolZoneRetTempSeq(i) * CoolFlowSeqzone(i))/ CoolFlowSeqsys(i)

SysCoolRetHumRat(i) = (CoolZoneHumRatSeq(i) * CoolFlowSeqzone(i))/ CoolFlowSeqsys(i)

FracOA = rair * DesOutAirVolFlowsys/ CoolFlowSeqsys(i)

Tmix = Toutside * FracOA + SysCoolRetTemp(i)(1 - FracOA)

Wmix = Woutside * FracOA + SysCoolRetHumRat (i)(1 - FracOA)

SysCoolOutTempSeq(i) = Toutside

SysCoolOutHumRatSeq(i) = Woutside

Get the current (zone time-step) system cooling capacity:

SysSensCoolCapcur = Cp,air * CoolFlowSeqsys(i) *( Tmix-Tsup)

SensCoolCapSeq(I) = SysSensCoolCapcur

If SysSensCoolCapcur is the maximum for the day so far then save SysSensCoolCapcur as the design value:

SensCoolCap(i )sys = SysSensCoolCapcur

And save the corresponding mixed, return and outside conditions:

CoolMixTempsys = Tmix

CoolMixHumRatsys = Wmix

CoolRetTempsys = SysCoolRetTemp(i)

CoolRetHumRatsys = SysCoolRetHumRat(I)

CoolOutTempsys = Toutside

CoolOutHumRatsys = Woutside

Here rair is the density of dry air at 20C and standard elevation corrected pressure, [kg/m3];

*FracOA* is the outside air fraction;
  • Cp,air* is the specific heat of dry air at 20C, [J/kg-K];

  • Tsup* is the user specified design cooling supply temperature [C];

  • Tmix* is the current mixed air temperature [C];

  • Wmix* is the current mixed air humidity ratio [kg water / kg dry air];

  • Toutside is the current outside air temperature [C];

  • Woutside is the current outside air humidity ratio [kg water / kg dry air].

  1. Loop over the zones heated by this air loop.

HeatFlowSeqsys(i) = HeatFlowSeq zone(i)

SysHeatRetTemp(i) = (HeatZoneRetTempSeq(i) * HeatFlowSeqzone(i))/

HeatFlowSeqsys(i)

SysHeatRetHumRat(i) = (HeatZoneHumRatSeq(i) * HeatFlowSeqzone(i))/

  • HeatFlowSeqsys(i)*

FracOA = rair * DesOutAirVolFlowsys/ HeatFlowSeqsys(i)

Tmix = Toutside * FracOA + SysHeatRetTemp(i)(1 - FracOA)

Wmix = Woutside * FracOA + SysHeatRetHumRat (i)(1 - FracOA)

SysHeatOutTempSeq(i) = Toutside

SysHeatOutHumRatSeq(i) = Woutside

Get the current (zone time-step) system heating capacity:

SysHeatCapcur = Cp,air * MinFlowRatsys *HeatFlowSeqsys(i) *( Tsup-Tmix)

HeatCapSeq(I) = SysHeatCapcur

If SysHeatCapcur is the maximum for the day so far then save SysHeatCapcur as the design value:

HeatCap(i )sys = SysHeatCapcur

And save the corresponding mixed, return and outside conditions:

HeatMixTempsys = Tmix

HeatMixHumRatsys = Wmix

HeatRetTempsys = SysHeatRetTemp(i)

HeatRetHumRatsys = SysHeatRetHumRat(I)

HeatOutTempsys = Toutside

HeatOutHumRatsys = Woutside

Here MinFlowRatsys is the user specified minimum supply flow ratio.

EndDay[LINK]

If the user has specified coincident system sizing then:

DesCoolVolFlowsys = rair *CoinCoolMassFlowsys

DesHeatVolFlowsys = rair *CoinHeatMassFlowsys

DesMainVolFlowsys = Max(DesCoolVolFlowsys, DesHeatVolFlowsys)

If the user has specified noncoincidentsystem sizing then:

DesCoolVolFlowsys = rair *NonCoinCoolMassFlowsys

DesHeatVolFlowsys = rair *NonCoinHeatMassFlowsys

DesMainVolFlowsys = Max(DesCoolVolFlowsys, DesHeatVolFlowsys)

Based on the outdoor air method selected, the DesCoolVolFlowsys and DesHeatVolFlowsys are modified based on the system ventilation effciency calculated based on the maximum outdoor air fraction.

EndSysSizingCalc[LINK]

At this point all the calculations have been done in SysSizing(i,j): we have results for each design day. Now these results need to be processed to find the heating and cooling design quantities for each system over all the design days.

For coincident sizing the task is quite easy.

  1. Loop over all of the air loops.
  1. Loop over all of the design days.

  2. If the value of DesCoolVolFlow in SysSizing for the current design day is greater than the value stored in CalcSysSizing, then move DesCoolVolFlow from SysSizing into CalcSysSizing along with CoolDesDay, CoinCoolMassFlow, SensCoolCap, CoolFlowSeq(i), SensCoolCapSeq(i), CoolMixTemp, CoolRetTemp, CoolMixHumRat, CoolRetHumRat, CoolOutTemp, CoolOutHumRat, SysCoolRetTempSeq(i), SysCoolRetHumRatSeq(i), SysCoolOutTempSeq(i) and SysCoolOutHumRatSeq(i).

  1. If the value of DesHeatVolFlow in SysSizing for the current design day is greater than the value stored in CalcSysSizing, then move DesHeatVolFlow from SysSizing into CalcSysSizing along with HeatDesDay, CoinHeatMassFlow, HeatCap, PreHeatCap, HeatFlowSeq(i), HeatCapSeq(i), PreHeatCapSeq(i), HeatMixTemp, HeatRetTemp, HeatMixHumRat, HeatRetHumRat, HeatOutTemp, HeatOutHumRat, SysHeatRetTempSeq(i), SysHeatRetHumRatSeq(i), SysHeatOutTempSeq(i) and SysHeatOutHumRatSeq(i).

At the end of each design day loop the peak cooling and the peak heating data will be stored in CalcSysSizing. At this point we set DesMainVolFlow in CalcSysSizing equal to the maximum of DesCoolVolFlow and DesHeatVolFlow.

For noncoincident sizing the task is harder since we don’t have a single time-step during which all the zone peaks occur. So there is no obvious value for outside air temperature at the peak, return air temperature at the peak and so forth. We must return to the zone sizing data and calculate average values for return and outside conditions.

  1. Loop over all of the zones cooled by this air loop

  2. In FinalZoneSizing replace the value in DesCoolCoilInTemp with the user specified CoolSupTempsys. Do the same for DesCoolCoilInHumRat and CoolSupHumRat. This ensures that zone equipment connected to an air loop will use the system design supply air conditions as coil entering conditions.

  1. NonCoinCoolMassFlowsys = DesCoolMassFlowzone

SysCoolRetTemp = (ZoneRetTempAtCoolPeak *DesCoolMassFlowzone)

/NonCoinCoolMassFlowsys

SysCoolRetHumRat = (ZoneHumRatAtCoolPeak *

DesCoolMassFlowzone)/NonCoinCoolMassFlowsys

SysCoolOutTemp = (TOA,zonepeak *DesCoolMassFlowzone)/

  • NonCoinCoolMassFlowsys*

SysCoolOutHumRat = (WOA,zonepeak *DesCoolMassFlowzone)/

  • NonCoinCoolMassFlowsys*

At the end of the zone loop calculate mixed air conditions and the system sensible cooling capacity.

FracOA = rair * DesOutAirVolFlowsys/ NonCoinCoolMassFlowsys

Tmix = SysCoolOutTemp * FracOA + SysCoolRetTemp * (1 - FracOA)

Wmix = SysCoolOutHumRat * FracOA + SysCoolRetHumRat *

(1 - FracOA)

SysSensCoolCap = Cp,air * NonCoinCoolMassFlow *( Tmix-Tsup)

Then (for noncoincident sizing) the variables calculated in section (ii) are moved into the CalcSysSizing Array.

  1. Loop over all of the zones heated by this air loop.

  2. In FinalZoneSizing replace the value in DesHeatCoilInTemp with the user specified HeatSupTempsys. Do the same for DesHeatCoilInHumRat and HeatSupHumRat. This ensures that zone equipment connected to an air loop will use the system design supply air conditions as coil entering conditions.

  1. NonCoinHeatMassFlowsys = DesHeatMassFlowzone

SysHeatRetTemp = (ZoneRetTempAtHeatPeak *DesHeatMassFlowzone)

/NonCoinHeatMassFlowsys

SysHeatRetHumRat = (ZoneHumRatAtHeatPeak *

DesHeatMassFlowzone)/NonCoinHeatMassFlowsys

SysHeatOutTemp = (TOA,zonepeak *DesHeatMassFlowzone)/

  • NonCoinHeatMassFlowsys*

SysHeatOutHumRat = (WOA,zonepeak *DesHeatMassFlowzone)/

  • NonCoinHeatMassFlowsys*

At the end of the zone loop calculate mixed air conditions and the system sensible cooling capacity.

FracOA = rair * DesOutAirVolFlowsys/ NonCoinHeatMassFlowsys

Tmix = SysHeatOutTemp * FracOA + SysHeatRetTemp * (1 - FracOA)

Wmix = SysHeatOutHumRat * FracOA + SysHeatRetHumRat *

(1 - FracOA)

SysHeatlCap = Cp,air * NonCoinHeatlMassFlow *( Tsup-Tmix)

Then (for noncoincident sizing) the variables calculated in section (ii) are moved into the CalcSysSizing Array.

  1. We now have the calculated system sizing data. This data needs to be altered to take into account the user input system design flow rates (if any), or the fact that the user may have requested that the system flow rate be sized on the ventilation requirement. Note that user specified sizing ratios have already been applied to the zone sizing data which have been used in out preceding system sizing calculation. Thus the user specified sizing ratios do not have to be explicitly taken into account at the system level.

First we move the calculated system sizing data from CalcSysSizing array into the FinalSysSizing array. FinalSysSizing will contain the user modified system design data when we are all done.

Loop over the air loops.

  1. As in the zone case, the user specified system design flow rates are turned into sizing ratios by dividing the user input value by the calculated value. The same strategy is employed for sizing on the ventilation requirement: the design ventilation flow rate is divided by the calculated design flow rate value. For each air loop this gives us a SizRatcool and SizRatheat.

CoinCoolMassFlow = SizRatcool * CoinCoolMassFlowcalc

NonCoinCoolMassFlow = SizRatcool * NonCoinCoolMassFlowcalc

DesCoolVolFlow = SizRatcool * DesCoolVolFlowcalc

Since the flow rates have been altered the outside air fraction will change. This will alter the design mixed air conditions and lead to an altered value for the cooling capacity. This must be done for the time-step sequence and for the peak value.

  1. Loop over the zone timesteps (index = i).

CoolFlowSeqsys(i) = SizRatcool * CoolFlowSeqsys,calc(i)

FracOA = rair * DesOutAirVolFlowsys/ CoolFlowSeqsys(i)

Tmix = SysCoolOutTempSeq(i) *FracOA +

SysCoolRetTempSeq(i) *(1-FracOA)

SensCoolCapSeq(i) = Cp,air * CoolFlowSeqsys(i) *( Tmix-Tsup)

  1. Do the same calculation for peak cooling.

FracOA = rair * DesOutAirVolFlowsys/ DesCoolVolFlow

Tmix = CoolOutTempsys *FracOA + CoolRetTempsys *(1-FracOA)

Wmix = CoolOutHumRatsys *FracOA + CoolRetHumRatsys *

(1-FracOA)

SensCoolCapsys = Cp,air * DesCoolVolFlowsys *( Tmix-Tsup)

Tmix and Wmix are saved in FinalSysSizing .

  1. Do the same calculation for the heating case.

CoinHeatMassFlow = SizRatheat * CoinHeatMassFlowcalc

NonCoinHeatMassFlow = SizRatheat * NonCoinHeatMassFlowcalc

DesHeatVolFlow = SizRatheat * DesHeatVolFlowcalc

  1. Loop over the zone timesteps (index = i).

HeatFlowSeqsys(i) = SizRatHeat * HeatFlowSeqsys,calc(i)

FracOA = rair * DesOutAirVolFlowsys/ HeatFlowSeqsys(i)

Tmix = SysHeatOutTempSeq(i) * FracOA +

  • SysHeatRetTempSeq(i) ** (1-FracOA)

HeatCapSeq(i) = Cp,air * HeatFlowSeqsys(i) *(Tsup-Tmix)

  1. Do the same calculation for peak heating.

FracOA = rair * DesOutAirVolFlowsys/ DesHeatVolFlow

Tmix = HeatOutTempsys *FracOA + HeatRetTempsys * (1-FracOA)

Wmix = HeatOutHumRatsys *FracOA + HeatRetHumRatsys *

(1-FracOA)

HeatCapsys = Cp,air * DesHeatVolFlowsys *( Tsup-Tmix)

Tmix and Wmix are saved in FinalSysSizing .

(vii)DesMainVolFlowsys = MAX(DesCoolVolFlowsys,DesHeatVolFlowsys)

This concludes the system design calculation.


Scalable System HVAC Sizing

The scalable system sizing applies to system supply air flow rates and sysyem capacity in coolin and heating modes.

Scalable System Air Flow Sizing

The scalable sizing methods for supply air flow rate allowed are either, FlowPerFloorArea, FractionOfAutosizedCoolingAirflow, or FlowPerCoolingCapacity. The scalable system air flow sizing methods are defined as follows:

FlowPerFloorArea means the program calculates the cooling supply air volume flow rate from zone floor area served by the airloop and user specified Flow Per Floor Area value.

FractionOfAutosizedCoolingAirflow means the program calculates the cooling supply air volume flow rate from user specified fraction and the autosized design cooling supply air volume flow rate value determined by the simulation.

FractionOfAutosizedHeatingAirflow means the program calculates the heating supply air volume flow rate from user specified fraction and the autosized design heating supply air volume flow rate value determined by the simulation.

FlowPerCoolingCapacity means the supply air volume is calculated from user specified flow per cooling capacity and design cooling capacity determined by the simulation.

FlowPerHeatingCapacity means the supply air volume is calculated from user specified flow per heating capacity and design heating capacity determined by the simulation.

Scalable System Capacity Sizing

The scalable sizing methods for system capacity available are: CapacityPerFloorArea, FractionOfAutosizedCoolingCapacity and FractionOfAutosizedHeatingCapacity. The scalable system capacity sizing methods are defined as follows:

CapacityPerFloorArea means the program calculates the design capacity from user specified capacity per floor area and floor area of the zones served by the airloop.

FractionOfAutosizedCoolingCapacity means the program calculates the design cooling capacity from user specified fraction and the auto-sized design cooling capacity.


FractionOfAutosizedHeatingCapacity means the program calculates the design heating capacity from user specified fraction and the auto-sized design heating capacity.