Coils[LINK]

Chilled-Water-Based Air Cooling Coil[LINK]

The input object Coil:Cooling:Water is simpler than the detailed geometry model. The simple model provides a good prediction of the air and water outlet conditions without requiring the detailed geometric input required for the detailed model. A greatly simplified schematic of enthalpy and temperature conditions in a counter flow cooling/dehumidifying coil is shown in the schematic Figure. The input required to model the coil includes only a set of thermodynamic design inputs, which require no specific manufacturer’s data. The coil simulation model is essentially a modification of one presented by Elmahdy and Mitalas (1977), TRNSYS, 1990 and Threlkeld, J.L. 1970. The model calculates the UA values required for a Dry, Wet and Part Wet & Part Dry Coil and iterates between the Dry and Wet Coil to output the fraction wet. There are two modes of flow operation for this model: Cross Flow, which is widely applicable in HVAC systems and the second being Counter flow mode. The default value in program is set up for Counter Flow. In addition the coil has two modes of analysis: Simple Analysis and Detailed Analysis. The Simple analysis mode operates the coil as either wet or dry while the detailed mode simulates the coil as part wet part-dry. While the detailed mode provides more accurate results, it is significantly slower than the simple model. The simple mode gives good results for an annual simulation but will not be adequate for a time step performance analysis.

Simplified Schematic of Cooling/Dehumidifying Coil

Heat Transfer and Energy Balance[LINK]

The cooling coil may be completely dry, completely wet with condensation, or it may have wet and dry sections. The actual condition of the coil surface depends on the humidity and temperature of the air passing over the coil and the coil surface temperature. The part-dry part-wet case represents the most general scenario for the coil surface conditions. There are subroutines present in the model for both the dry and wet regions of the coil, and a subroutine that iterates between the dry and wet subroutines to calculate the fraction of the coil surface that is wet. For each region the heat transfer rate from air to water may be defined by the rate of enthalpy change in the air and in the water. The rates must balance between each medium for energy to be conserved.

Model Description[LINK]

The Model has two blocks: 1 $^{s t}$ = Design Block with the Design Inputs. This block calculates the Design U-Factor Times Area Value (UA) values required by the model. Using these UA values the model simulates the operating conditions. The operating block is the one containing the operating conditions, the conditions at which the coil operates. Following is the list of Design and Operating inputs and subsequently the Design and Operating variables used in the model.

Design Inputs (User Inputs)
Input Field	Description
DesWaterVolFlowRate:	Maximum Water Volume Flow Rate
DesAirVolFlowRate:	Maximum Air Volume Flow Rate
DesInletWaterTemp:	Inlet Water Temperature at Design Condition
DesInletAirTemp:	Inlet Air Temperature at Design Condition
DesOutletAirTemp:	Outlet Air Temperature at Design Condition
DesInletAirHumRat:	Inlet Air Humidity Ratio at Design Conditions
DesOutletAirHumRat:	Outlet Air Humidity Ratio at Design Conditions.

Operating Conditions (From Nodes – not user inputs)
Condition Variable	Description
InletWaterMassFlowRate:	Entering Water Mass Flow Rate at operating condition
InletWaterTemp:	Inlet Water Temperature at operating condition
InletAirMassFlowRate:	Entering Air Mass Flow Rate at operating condition
InletAirTemp:	Inlet Air Temperature at operating condition
InletAirHumRat:	Entering air humidity ratio at operating conditions

Intermediate calculated U-Factor Times Area Values: The Crux of the Model[LINK]

The various U-Factor Times Area values (UA) required by this model are calculated from the above inputs, which are explained later in the document. The various UA are:

UA Descriptions of Model
UA Variable Name	Description
CoilUATotal:	Overall heat transfer coefficient (W/C)
CoilUAInternal:	Overall internal UA (W/C)
CoilUAExternal:	Overall external UA (W/C)
CoilUInternal:	Internal overall heat transfer coefficient (W/m ∙C)
CoilUWetExternal:	Wet part external overall heat transfer coefficient (W/m ∙C)
CoilUDryExternal:	Dry part external overall heat transfer coefficient (W/m ∙C)

The UA values are calculated assuming a wet coil at the design conditions. Following are a few important calculations to understand the working of the model. The model is basically divided into two blocks: the Design Block and the Operating Block.

The Design Block is a one time calculation. The aim of the Design Block is to calculate the Coil UA for use in the operating Block.

Design Block Calculations:[LINK]

The design block has the code for calculating the six Coil UA values required by the operating block. Reasonable assumptions have been made in the calculations to maintain the simplicity of the model.

Heat transfer in a wet coil model is based on enthalpy rather than temperature to take into account latent effects. While heat transfer rates are commonly expressed as the product of an overall heat transfer coefficient, UA, and a temperature difference, the use of enthalpy-based heat transfer calculations requires an enthalpy-based heat transfer coefficient which we denote as DesUACoilTotalEnth and hence the equation.

Q = DesUACoilTotalEnth * (H $_{a i r, m e a n}$ - H $_{w a t e r, m e a n}$ ). The value of Q is calculated using product of air mass flow rate and difference in inlet and outlet air enthalpies at design conditions.

The relation between the enthalpy-based UA and the temperature-based UA is DesUACoilTotalEnth = CoilUA / C $_{p}$ . CoilUA is the conventional heat transfer coefficient and C $_{p}$ = specific heat of the air.

We need the following quantities for our design calculations. The Psy functions are the EnergyPlus built-in psychrometric functions.

${˙ m}_{a i r} = ρ_{a i r} {˙ V}_{a i r}$

$h_{a i r, i n} = PsyHFnTdbW (T_{a i r, i n}, w_{a i r, i n})$

$h_{a i r, o u t} = PsyHFnTdbW (T_{a i r, o u t}, w o u t)$

$h_{w, s a t, i n} = P s y H F n T d b W (T_{w, i n}, P s y W F n T d p P b (T_{w, i n}, P_{a t m}))$

${˙ Q}_{c o i l} = {˙ m}_{a i r} (h_{a i r, i n} - h_{a i r, o u t})$

$T_{w, o u t} = T_{w, i n} + {˙ Q}_{c o i l} / ({˙ m}_{w, max} C_{p, w})$

$h_{w, s a t, o u t} = P s y H F n T d b W (T_{w, o u t}, P s y W F n T d p P b (T_{w, o u t}, P_{a t m}))$

We now calculate the design coil bypass factor. The bypass factor is not used in subsequent calculations. It is calculated solely to use as check on the reasonableness of the user-input design inlet and outlet conditions. First we make an initial estimate of the apparatus dew point temperature:

$T_{a i r, d p, a p p} = P s y T d p F n W P b (w_{a i r, o u t}, P_{a t m})$

we also need the “slope” of temperature versus humidity ratio on the psych chart between the inlet and outlet air conditions:

$S_{T, w} = (T_{a i r, i n} - T_{a i r, o u t}) / (w_{a i r, i n} - w_{a i r, o u t})$

We now obtain the actual design apparatus dewpoint temperature by iterating over the following two equations:

$w_{a i r, d p, a p p} = P s y W F n T d p P b (T_{a i r, d p, a p p}, P_{a t m})$

$T_{a i r, d p, a p p} = T_{a i r, i n} - S_{T, w} (w_{a i r, i n} - w_{a i r, d p, a p p})$

The apparatus dewpoint enthalpy is then:

$h_{a i r, d p, a p p} = P s y H F n T d b W (T_{a i r, d p, a p p}, w_{a i r, d p, a p p})$

and the coil bypass factor is:

$F_{c o i l b y p a s s} = (h_{a i r, o u t} - h_{a i r, d p, a p p}) / (h_{a i r, i n} - h_{a i r, d p, a p p})$

If the iterative procedure doesn’t converge, or the coil bypass factor is too large (greater than 0.5), or the apparatus dewpoint enthalpy is less than the saturated air enthalpy at the water inlet temperature, the design outlet air conditions are reset to 90% relative humidity at the same outlet enthalpy. The above design calculations are then repeated.

We are now ready to calculate the design coil UA. This will be accomplished by inverting the simple coil calculation routine CoolingCoil using the Regula Falsi method. First we make an initial estimate of the coil UA.

$Δ h_{l m d} = ((h_{a i r, i n} - h_{w, s a t, o u t}) - (h_{a i r, o u t} - h_{w, s a t, i n})) / l o g ((h_{a i r, i n} - h_{w, s a t, o u t}) - (h_{a i r, o u t} - h_{w, s a t, i n}))$

$U A_{c o i l, e n t h a l p y - b a s e d} = {˙ Q}_{c o i l} / Δ h_{l m d}$

$U A_{c o i l, e x t} = C_{p, a i r} U A_{c o i l, e n t h a l p y - b a s e d}$

We set the internal UA to 3.3 times the external UA (as a typical value for a coil). Then the total UA is:

$U A_{c o i l, t o t} = \frac{1}{(1 / U A_{c o i l, i n t} + 1 / U A_{c o i l, e x t})}$

The next step is to estimate the coil external heat transfer surface area. This is done in the function EstimateHEXSurfaceArea:

$A r e a_{c o i l, e x t} = E s t i m a t e H E X S u r f a c e A r e a$

using the following assumptions:

Tube inside diameter = 0.0122 (m)
Tube side water velocity = 2.0 (m/s)
Inside to outside coil surface area ratio (Ai/Ao) = 0.07 (-)
Fins overall efficiency = 0.92 (-)
Aluminum fins, 12 fins per inch with fins to total outside surface area ratio of 90%.
Airside combined heat and mass transfer coefficient = 140 (W/m2∙°C)

Interior and exterior U values (really UA’s per unit exterior surface area) are calculated by dividing the above UA’s by the area. The resulting U $_{c o i l, e x t}$ is assumed to be U $_{c o i l, e x t, w e t}$ ; U $_{c o i l, e x t, d r y}$ is set equal to U $_{c o i l, e x t, w e t}$ . We now have all the starting values needed for inverting the simple coil model using the chosen Regula Falsi iterative method. Once the iteration is completed, we have coil UA’s and U’s that yield the design outlet air and water enthalpies given the inlet design conditions and flow rates. Note that the simple coil model can not exactly match the specified design outlet air temperature and humidity ratio. It can only match the design air outlet enthalpy. Generally the simple coil model will yield outlet conditions near the saturation curve if any dehumidification is occuring. Typical outlet relative humidities are around 95%.

Variable UA[LINK]

The above calculations yield coil UA’s for the design inlet conditions and air and water flow rates. As the flow rates vary during the time step calculations, the UA’s need to be adjusted, since coil UA’s are a rather strong function of air and water side flow rates. Each time step the coil UA’s are modified using the same formulas as are used in the hot water coil model. Refer to that model for the flow dependences.

Operating Block Calculations:[LINK]

There are two modes of coil analysis in the operating block. They are the Simple analysis mode and the detailed analysis mode. The simple analysis mode assumes the coil to be either all wet or either all dry and execute the model , on the other hand the detailed mode checks for part wet part dry mode of operation and reports surface area wet fraction of coil, however the program execution time in detailed mode is noticeably higher.

The operating block for Detailed Mode Analysis of this coil model is divided into three modes of coil performance. The modes being

Coil is completely dry: There is no moisture condensation on the coil surface and the coil is a dry coil. This is an extreme condition when the entering air has very low humidity ratio or is dry air.
Coil is completely wet: The entire coil is wet due to complete condensation on the surface of the coil.
Part Wet Part Dry Mode: This is the usual/frequent mode of operation of coil, as shown in Figure, where part of the coil at entry of air is dry and as air cools condensation occurs and part of the coil becomes wet.

The Part Wet Part Dry Mode of operation is essentially a function the Coil Completely Dry and Coil Completely Wet mode. This subroutine iterates between the Dry Coil and the Wet Coil to give outputs, a detailed explanation is given later in the document. The operating block requires 5 inputs, which are mentioned earlier in the document. These inputs are automatically generated from the node connections in Energy Plus. The user does not have to input any information to run this coil model.

The option to identify which mode of operation the coil should perform ie, for a given set of inputs would the coil be Dry, Wet or Part Wet Part Dry, is decided by set of conditions described below.

IF (Temperature Dewpoint Air < Water Inlet Temperature) THEN the coil is Dry and we call the Subroutine Coil Completely Dry. In this case outlet temperature of air would be higher than the air dewpoint and hence there would be no condensation.
IF (Temperature Dewpoint Air > Water Inlet Temperature) THEN the coil is completely wet, call subroutine Coil Completely Wet, it is assumed that moisture condensation occurs over completely surface of the coil. However we go ahead and check for the coil being partially wet with the following condition.
IF (AirDewPointTemp < AirInletCoilSurfTemp) THEN, the coil is Partially Wet because there is possibility that air temperature will go below its dewpoint and moisture will condense on latter part of the cooling coil.

The Operating Block for Simple Mode Analysis is divided into two modes of coil performance, the two modes being

Coil is completely dry: There is no moisture condensation on the coil surface and the coil is a dry coil.
Coil is completely wet: The entire coil is wet due to complete condensation on the surface of the coil.

The option to identify which mode of operation the Simple mode analysis should perform ie, for a given set of inputs would the coil be Dry or Wet is decided by set of conditions described below.

IF (Temperature Dewpoint Air < Water Inlet Temperature) THEN the coil is Dry and we call the Subroutine Coil Completely Dry. In this case outlet temperature of air would be higher than the air dewpoint and hence there would be no condensation.
IF (Temperature Dewpoint Air > Water Inlet Temperature) THEN the coil is completely wet, call subroutine Coil Completely Wet, it is assumed that moisture condensation occurs over completely surface of the coil. However we go ahead and check for the coil being partially wet with the following condition.

The above is a simple mode of analysis and the results are very slightly different from the detailed mode of analysis. The algorithms used in Simple mode and the Detailed mode are identically similar. The surface area wet fraction in the coil is reported as 1.0 or 0.0 for wet or dry coil respectively. The program defaults to simple mode of analysis for enabling higher execution speed.

Effectiveness Equations:[LINK]

There are two modes of flow for the coil, Counter Flow mode or the Cross Flow mode, default set up is as cross flow since most air condition applications have cross flow heat exchangers. According to the mode of flow the following NTU - Effectiveness relationships are used to calculate coil effectiveness, which is used later by all the three modes (Dry, Wet, Part Wet) for calculating air outlet conditions and heat transfer.

Following are the relations used for calculating effectiveness equation for the Heat exchangers.

Counter Flow Heat Exchanger: Effectiveness Equation:

$η_{C o u n t e r F l o w} = \frac{(1 - E x p (- N T U \times (1 - R a t i o_{S t r e a m C a p a c i t y})))}{1 - R a t i o_{S t r e a m C a p a c i t y} \times E x p (- N T U \times (1 - R a t i o_{S t r e a m C a p a c i t y}))}$

In Equation the variable Ratio_StreamCapacity is defined as below

$R a t i o_{S t r e a m C a p a c i t y} = \frac{M i n C a p a c i t y_{S t r e a m}}{M a x C a p a c i t y_{S t r e a m}}$

In equation capacity of stream is defined as below in equation

$(M i n, M a x) C a p a c i t y_{S t r e a m} = (M a s s F l o w R a t e \times C p)_{a i r, w a t e r}$

NTU in equation , is defined as the Number of Transfer Units, it is a function of Coil UA and the Minimum Capacity of Stream. The Coil UA is a variable in this equation and depends on which mode of the coil operation (Dry, Wet, Part Wet) is calling upon equation , i.e., if it is Coil Completely Dry calling upon the effectiveness equation with the value of Dry UA total, which in our case is defined as CoilUA_total. Equation gives definition for NTU.

$N T U = \frac{C o i l U A}{M i n S t r e a m_{C a p a c i t y}}$

Cross Flow Heat Exchanger: Effectiveness Equation:

$η_{C r o s s F l o w} = 1 - E X P {\frac{E x p (- N T U \times R a t i o_{S t r e a m C a p a c i t y} \times N T U^{- 0.22}) - 1}{R a t i o_{S t r e a m C a p a c i t y} \times N T U^{- 0.22}}}$

The variables in the above equation have already been defined earlier. Depending on the mode of operation of the coil model the cross or the counter flow equations are used to calculate the effectiveness.

Coil Outlet Conditions:[LINK]

Calculating the Outlet Stream Conditions using the effectiveness value from equation or depending on the mode of flow. The energy difference between the outlet and inlet stream conditions gives the amount of heat transfer that has actually take place. Temperature of air and water at outlet to the coil is given as in following equations

$T e m p A i r_{O u t} = T e m p A i r_{i n l e t} - η_{c r o s s, c o u n t e r} \times \frac{M a x H e a t T r a n s f e r}{S t r e a m C a p a c i t y_{A i r}}$

$T e m p W a t e r_{O u t} = T e m p W a t e r_{I n l e t} + η_{C r o s s, c o u n t e r} \times \frac{M a x H e a t T r a n s f e r}{S t r e a m C a p a c i t y_{W a t e r}}$

In the above equations and the maximum heat transfer is calculated as shown in the following equation

$M a x H e a t T r a n s f e r = M i n S t r e a m_{C a p a c i t y} \times (T e m p A i r_{I n l e t} - T e m p W a t e r_{I n l e t})$

Coil Completely Dry Calculations: (operating block)[LINK]

Since the coil is dry, the sensible load is equal to total load and the same with the humidity ratios at inlet and outlet, as in equations and .

$Q S e n s i b l e_{D r y C o i l} = Q T o t a l_{D r y C o i l}$

$H u m R a t i o_{I n l e t} = H u m R a t i o_{O u t l e t}$

Total Heat Transfer in dry coil is as follows:

$Q_{T o t a l D r y C o i l} = C a p a c i t y A i r \times (A i r T e m p_{I n} - A i r T e m p_{O u t l e t})$

The variables in the above equation are calculated earlier in equations and to give the total cooling load on the coil.

Coil Completely Wet Calculations: (operating block)[LINK]

In wet coil we need to account for latent heat transfer, hence calculations are done using enthalpy of air and water instead of stream temperatures Hence we need to define coil UA for the wet coil based on enthalpy of the operating streams and not design streams.

Similar to equations and we calculate the air outlet enthalpy and water outlet enthalpy ie by replacing temperature with enthalpy of the respective streams. The input variable for Coil UA in equation for calculating NTU, in this case it would be enthalpy based and is given as shown in equation

$C o i l U A_{E n t h a l p y B a s e d} = \frac{1}{(\frac{C p S a t_{I n t e r m e d i a t e}}{C o i l U A_{I n t e r n a l}} + \frac{C p_{A i r}}{C o i l U A_{E x t e r n a l}})}$

Total Coil Load in case of Wet Coil is the product of mass flow rate of air and enthalpy difference between the inlet and outlet streams as given in the following equation

$Q_{T o t a l} = * M a i r \times (E n t h A i r_{I n l e t} - E n t h A i r_{O u t l e t})$

Once the enthalpy is known the outlet temperatures and outlet humidity ratios of the wet coil are calculated as in equations below.

IF (TempCondensation < PsyTdpFnWPb(InletAirHumRat ,Patm)) THEN

$A i r T e m p_{O u t} = A i r T e m p_{i n l e t} - (A i r T e m p_{i n l e t} - C o n d e n s a t i o n_{T e m p}) \times η$

and

OutletAirHumdityRatio = PsyWFnTdbH(OutletAirTemp,EnthAirOutlet)

ELSE

There is no condensation and hence the inlet and outlet Hum Ratios are equal , and outlet temperature is a function of outlet air enthalpy as below

OutletAirTemp = PsyTdbFnHW (EnthalpyAirOutlet, OutletAirHumRat)

and

OutletAirHumRat = InletAirHumRat

ENDIF

Effectiveness η used in equation is defined in equation and Condensation Temperature is calculated using psychrometric function as in equation .

$η = 1 - E x p {- \frac{C o i l U A_{E x t e r n a l}}{C a p a c i t a n c e_{A i r}}}$

$C o n d e n s a t i o n_{T e m p} = P s y T s a t F n H P b (E n t h_{A i r C o n d e n s a t e T e m p}, P a t m)$

$E n t h_{A i r C o n d e n s a t e T e m p} = E n t h_{A i r I n l e t} - \frac{(E n t h_{A i r I n l e t} - E n t h_{A i r O u t l e t})}{η}$

Once the air outlet temperature are known, then sensible load is calculated as a product of capacitance of air and temperature difference at inlet and outlet, as in equation

$Q_{S e n s i b l e} = C a p a c i t a n c e_{A i r} \times (A i r T e m p_{I n l e t} - A i r T e m p_{O u t l e t})$

Coil Part Wet Part Dry Calculations: (operating block)[LINK]

The Coil would perform under part wet part dry conditions when Air Dewpoint Temperature is less than Coil surface temperature at inlet to air. In this case part of the coil used value of Dry UA for heat transfer and part the coil used Wet UA value for heat transfer.

This problem is solved utilizing the fact that the Exit conditions from the Dry Part of the Coil would become the inlet conditions to the wet part of the coil (see Figure) and the coil model determines by iteration what fraction of the coil is wet and based on that it calculates the areas and subsequently the UA values of that dry and wet part, based on the area of the dry and wet part respectively. Explained below are the steps followed to the estimating the wet dry behavior of the coil.

Iterate between the Dry Coil and the Wet Coil. First calculate Coil Completely Dry performance by estimating the wet dry interface water temperature using equation and inputting this variable as the water inlet temperature to dry Coil.

$W e t D r y I n t e r f a c e_{W a t e r T e m p} = W a t e r T e m p_{I n l e t} + A r e a_{W e t F r a c t i o n} * (W a t e r T e m p_{O u t l e t} - W a t e r T e m p_{I n l e t})$

The value of Surface Area Wet fraction is estimated initially as follows

$A r e a_{W e t F r a c t i o n E s t i m a t e} = \frac{A i r D e w P t_{T e m p} - I n l e t W a t e r_{T e m p}}{O u t l e t W a t e r_{T e m p} - I n l e t W a t e r_{T e m p}}$

For the above mentioned iteration the value of Coil UA for Wet and Dry part need to be varied according to the new respective area of the wet and dry parts. This estimate of Wet and Dry area is a product of the estimated Surface Area Fraction and total coil external area, which keeps varying as will be explained further in the document.

UA value for Dry part of the Coil is estimated as below.

$C o i l U A_{D r y E x t e r n a l} = \frac{S u r f A r e a_{D r y}}{\frac{1}{C o i l U_{D r y E x t e r n a l}} + \frac{1}{C o i l U_{I n t e r n a l}}}$

Where Surface Area Dry = (Total Coil Area – Wet Part Area), where the Wet part area is the product of Surface fraction Wet and Total Coil Area.

UA value for the Wet part of the Coil requires Wet UA external and Wet UA Internal, which are calculated as below.

$W e t P a r t U A_{E x t e r n a l} = C o i l U_{W e t E x t e r n a l} \times S u r f a c e A r e a W e t$

$W e t P a r t U A_{I n t e r n a l} = C o i l U_{I n t e r n a l} \times S u r f a c e A r e a W e t$

It is essential to remember that the mode of calculation for the coils remains the same as in completely wet and completely dry mode, only the UA values and water, air outlet and inlet values change.

Now Iterate between the Dry Coil and wet Coil with the above respective UA, and usual operating inputs except the variable water inlet temperature for dry Coil is replaced with Wet Dry Interface Water temperature, and in the Wet Coil the Outlet Air Temperature from dry Coil is the inlet air temperature to Wet Coil. The iteration proceeds till the Outlet Water Temperature from Wet Coil equals the Wet Dry Interface Water Temp, which is the input to Dry Coil.

Dry Part Inputs: (changed operating inputs) :Iteration Case 1: Explained In Programming Fashion:

CALL CoilCompletelyDry (WetDryInterfcWaterTemp, InletAirTemp, DryCoilUA, &

OutletWaterTemp, WetDryInterfcAirTemp, WetDryInterfcHumRat, &

DryCoilHeatTranfer).

Input the calculated values calculated by Dry Coil above into Wet Coil below. The variables have been highlighted in color red and blue.

CALLCoilCompletelyWet (InletWaterTemp, WetDryInterfcAirTemp, WetDryInterfcHumRat

WetPartUAInternal,WetPartUAExternal, &

EstimateWetDryInterfcWaterTemp, OutletAirTemp, OutletAirHumRat, &

WetCoilTotalHeatTransfer, WetCoilSensibleHeatTransfer, &

EstimateSurfAreaWetFraction, WetDryInterfcSurfTemp)

Iterate Between the above two Wet and Dry Coil calls until the two variables in blue ie WetDryInterfcWaterTemp = EstimateWetDryInterfcWaterTemp.The key is to have the difference between the variables (WetDryInterfcWaterTemp – OutletWaterTemp) in Dry Coil equal to (InletWaterTemp-EstimatedWetDryInterfcWaterTemp) in Wet Coil. This equality quantized the relative part of coil that is dry and part that is wet on the basis of heat transfer that has occurred.

After the above convergence check for the coil being dry otherwise iterate to calculate surface fraction area wet.

${(A r e a F r a c t i o n_{W e t} \leq 0.0) a n d (W e t D r y I n t e r f a c e_{S u r f T e m p} > A i r D e w P t)}$

THEN CoilCompletelyDry

If equation is satisfied then Coil is Dry and simply output the value for Dry Coil calculated else the coil is partially wet and then iterate to find the surface fraction area wet. Start with the initially guess value of surface area fraction (equation wet and iterate on the entire loop starting from until the Wet Dry Interface Temperature equals the Air Dewpoint Temperature. The value of Surface Area fraction wet at which the interface air temperature equals is dewpoint is the transition point from wet to dry and gives the % of coil that is dry and % that is wet.

Graphs Showing the Performance of the coil model at optimum operating conditions are shown below. All values of variable used have been normalized.

Air Outlet Temperature Vs Air Mass Flow Rate

Sensible Load variations Vs Air mass Flow Rate

Total and Sensible Load variations Vs Air Mass Flow Rate

Surface Area Fraction Wet Vs Air Mass Flow Rate

References[LINK]

IBPSA BuildSim-2004. 2004. Colarado Boulder: An Improvement of Ashrae Secondary HVAC toolkit Simple Cooling Coil Model for Building Simulation, Rahul J Chillar, Richard J Liesen M&IE ,UIUC.

Stoecker, W.F. <dates unspecified> Design of Thermal Systems,: ME 423 Class Notes , M& IE Dept UIUC.

Brandemeuhl, M. J. 1993. HVAC2 Toolkit: Algorithms and Subroutines for Secondary HVAC Systems Energy Calculations, ASHRAE.

Elmahdy, A.H. and Mitalas, G.P. 1977. “A Simple Model for Cooling and Dehumidifying Coils for Use In Calculating Energy Requirements for Buildings ASHRAE Transactions, Vol.83 Part 2, pp. 103-117.

Threlkeld, J.L. 1970. Thermal Environmental Engineering, 2nd Edition, Englewood Cliffs: Prentice-Hall,Inc. pp. 254-270.

ASHRAE Secondary HVAC Toolkit TRNSYS. 1990. A Transient System Simulation Program: Reference Manual. Solar Energy Laboratory, Univ. Wisconsin-Madison, pp. 4.6.8-1 - 4.6.8-12.

Kays, W.M. and A.L. London. 1964. Compact Heat Exchangers, 2nd Edition, New York: McGraw-Hill.

Clark, D.R.. 1985. HVACSIM+ Building Systems and Equipment Simulation Program Reference Manual, Pub. No. NBSIR 84-2996, National Bureau of Standards, U.S. Department of Commerce, January, 1985

Elmahdy, A.H. 1975. Analytical and Experimental Multi-Row Finned-Tube Heat Exchanger Performance During Cooling and Dehumidifying Processes, Ph.D. Thesis, Carleton University, Ottawa, Canada, December, 1975.

Elmahdy, A.H., and Mitalas, G.P. 1977. “A Simple Model for Cooling and Dehumidifying Coils for Use in Calculating Energy Requirements for Buildings,” ASHRAE Transactions, Vol. 83, Part 2, pp. 103-117.

Chilled-Water-Based Detailed Geometry Air Cooling Coil[LINK]

The input object Coil:Cooling:Water:DetailedGeometry provides a coil model that predicts changes in air and water flow variables across the coil based on the coil geometry. A greatly simplified schematic of enthalpy and temperature conditions in a counterflow cooling/dehumidifying coil is shown in the following schematic figure. In addition, the variables required to model a cooling/dehumidifying coils and their definitions are extensively listed in “Table. Coil Geometry and Flow Variables for Coils”. The input required to model the coil includes a complete geometric description that, in most cases, should be derivable from specific manufacturer’s data. The coil simulation model is essentially the one presented by Elmahdy and Mitalas (1977) and implemented in HVACSIM+ (Clark 1985), a modular program also designed for energy analysis of building systems. The model solves the equations for the dry and wet sections of the coil using log mean temperature and log mean enthalpy differences between the liquid and the air streams. Elmahdy and Mitalas state that crossflow counterflow coils with at four rows or more are approximated well by this model. This does not constitute a major limitation since cooling and dehumidifying coils typically have more than four rows.

Simplified Schematic of Cooling/Dehumidifying Coil

Heat Transfer and Energy Balance[LINK]

The cooling coil may be completely dry, completely wet with condensation, or it may have wet and dry sections. The actual condition of the coil surface depends on the humidity and temperature of the air passing over the coil and the coil surface temperature. The partly wet-partly dry case represents the most general scenario for the coil surface conditions. The all dry and all wet cases can be considered as limiting solutions of the wet or dry areas respectively going to zero. In the general case, equations are written for both the dry and wet regions of the coil. For each region the heat transfer rate from air to water may be defined by the rate of enthalpy change in the air and in the water. The rates must balance between each medium for energy to be conserved. Equations through express the energy balance between the water and the air for the case of dry and wet coils respectively. Equations and represent the heat transfer rate between water and air based on the actual performance of the coil. The UA parameter can be calculated from the parameters in the following table.

Coil Geometry and Flow Variables for Coils
A	area	LMHD	log mean enthalpy difference
A	air, air side	LMTD	log mean temperature difference
aa, bb	coeff. in enthalpy approximation	$˙ m$	mass flow rate
C1, C2	coeff. in air side film coeff.	mf	metal and fouling
Cp	specific heat	$μ$	viscosity
D	diameter, effective diameter	o	outside (air side)
Dhdr	hydraulic diameter on air side	Pr	Prandtl number
D	dry region	$˙ Q$	heat transfer rate
$δ$	thickness	R	overall thermal resistance
$Δ$	spacing	Re	Reynolds number
F	heat transfer film coefficient	$ρ$	ratio of diameters
Fai	variable in fin eff. calculation	s	surface, outside of metal
fin, fins	air side fin geometry	St	Stanton number
H	enthalpy	T	temperature
$η$	efficiency	tube	water tube
I0()	mod Bessel fn, 1st kind, ord 0	UAdry	dry heat xfer coeff. * dry area
I1()	mod Bessel fn, 1st kind, ord 1	UcAw	wet heat xfer coeff. * wet area
K0()	mod Bessel fn, 2nd kind, ord 0	ub, ue	variables in fin eff. calculation
K1()	mod Bessel fn, 2nd kind, ord 1	V	average velocity
I	inside (water side)	w	water, water side, or wet region
K1	variable in sol’n form of eq.	wa	humidity ratio
K	thermal conductivity	Z	variables in sol’n form of eq.
L	length	1, 2, 3	positions (see diagram)

Equations through represent two sets of three equations with 7 unknowns: ${˙ Q}_{d}$ , Ta,1, Ta,2, Tw,2, Tw,3, ${˙ m}_{a}$ , ${˙ m}_{w}$ . However, normally at least four of these variables are specified, for example: inlet water temperature, outlet air temperature, water flow rate, air flow rate, so that the system of equations is effectively closed.

${˙ Q}_{d} = m_{a} C p_{a} (T_{a, 1} - T_{a, 2})$

${˙ Q}_{d} = m_{w} C p_{w} (T_{w, 3} - T_{w, 2})$

${˙ Q}_{d} = (U A_{d r y}) (L M T D)$

${˙ Q}_{w} = m_{a} (H_{a, 2} - H_{a, 3})$

${˙ Q}_{w} = m_{w} C p_{w} (T_{w, 2} - T_{w, 1})$

${˙ Q}_{w} = (U_{c} A_{w}) (L M H D)$

In order to manipulate these equations, the log mean temperature and enthalpy differences are expanded as shown in Equations and . Finally, a linear approximation of the enthalpy of saturated air over the range of surface temperature is made using Equation . Note that in Equation Hw refers to the enthalpy of saturated air at the water temperature.

$L M T D = \frac{(T_{a, 1} - T_{w, 3}) - (T_{a, 2} - T_{w, 2})}{ln \frac{T_{a, 1} - T_{w, 3}}{T_{a, 2} - T_{w, 2}}}$

$L M H D = \frac{(H_{a, 2} - H_{w, 2}) - (H_{a, 3} - H_{w, 1})}{ln \frac{H_{a, 2} - H_{w, 2}}{H_{a, 3} - H_{w, 1}}}$

$H_{w} = a a + b b T_{w}$

Equation is derived from the above equations and is used to solve for the coil conditions when all of the inlet conditions are given as input. Operating in this manner, the coil does not have a controlled outlet air temperature.

$T_{w, 2} = \frac{(1 - Z) (H_{a, 1} - a a - K 1 C p_{a} T_{a, 1}) + Z T_{w, 1} (b b - \frac{m_{w} C p_{w}}{m_{a}})}{b b - Z \frac{m_{w} C p_{w}}{m_{a}} - (1 - Z) K 1 C p_{a}}$

An alternative solution method is to define the coil leaving air temperature as an input with a variable water flow rate. In this case Equations and are more convenient. Equations through define terms that are used to simplify Equations , and .

$T_{w, 2} = \frac{(1 - Z) (H_{a, 3} - a a) + T_{w, 1} (\frac{m_{w} C p_{w}}{m_{a}} - b b Z)}{\frac{m_{w} C p_{w}}{m_{a}} - b b}$

$T_{w, 2} = \frac{(Z_{d} - 1) T_{a 1} C p_{a} + T_{w, 3} (C p_{a} - Z_{d} \frac{m_{w} C p_{w}}{m_{a}})}{Z_{d} (C p_{a} - \frac{m_{w} C p_{w}}{m_{a}})}$

$Z = e x p (U_{c} A_{w} (\frac{1}{m_{a}} - \frac{b b}{m_{w} C p_{w}}))$

$K 1 = \frac{Z_{d} - 1}{Z_{d} - \frac{m_{a} C p_{a}}{m_{w} C p_{w}}}$

$Z_{d} = e x p (U_{c} A_{d r y} (\frac{1}{m_{a} C p_{a}} - \frac{1}{m_{w} C p_{w}}))$

Underlying Correlations, Properties, and Assumptions[LINK]

Overall heat transfer coefficients are calculated from the specified coil geometry and by using empirical correlations from fluid mechanics and heat transfer. For the water side, Equation gives the film heat transfer coefficient in SI units:

$f_{i} = 1.429 (1 + 0.0146 T_{w}) V_{w}^{0.8} D_{i}^{- 0.2}$

This is valid for Reynolds numbers greater than 3100 based on water flow velocity and pipe inside diameter and is given in Elmahdy and Mitalas(1977) as recommended in the standard issued by the Air-Conditioning and Refrigeration Institute (1972) for air-cooling coils. The definition of overall inside thermal resistance follows directly as shown in Equation.

$R_{i} = \frac{1}{f_{i} A_{i}}$

Equation gives the film coefficient for the air side. Another form of the same equation is Equation , which is familiar from the data presented in Kays and London (1984). For coil sections that have a wet surface due to condensation, the air side film coefficient is modified according to Equation . The correction term, a function of air Reynolds number, is valid for Reynolds numbers between 400 and 1500. The coefficients in Equation and are calculated by Equations and that are functions of the coil geometry. Elmahdy (1977)explains the modifier for the wet surface and coefficients for the film coefficient. Equations through show definitions and values of common parameters and properties.

$f_{o} = C_{1} R e_{a}^{C_{2}} \frac{m_{a}}{A_{a_{m} i n_{f} l o w}} C_{p, a} P r_{a}^{2 / 3}$

$C_{1} R e_{a}^{C_{2}} = S t_{a} P r_{a}^{2 / 3}$

$f_{o, w} = f_{o} (1.425 - 5.1 \times 10^{- 4} R e_{a} + 2.63 \times 10^{- 7} R e_{a}^{2})$

$C_{1} = 0.159 {(\frac{δ_{f i n}}{D_{h d r}})}^{- 0.065} {(\frac{δ_{f i n}}{L_{f i n}})}^{0.141}$

$C_{2} = - 0.323 {(\frac{Δ_{f i n s}}{L_{f i n}})}^{0.049} {(\frac{D_{f i n}}{Δ_{t u b e_{r} o w s}})}^{0.549} {(\frac{δ_{f i n}}{Δ_{f i n s}})}^{- 0.028}$

$D_{h d r} = \frac{4 A_{a_{m} i n_{f} l o w} δ_{c o i l}}{A_{s_{t} o t a l}}$

$R e_{a} = \frac{4 δ_{c o i l} (1 + w_{a}) m_{a}}{A_{s_{t} o t a l} μ_{a}}$

$P r_{a} = 0.733$

$μ_{a} = 1.846 \times 10^{- 5}$

The film coefficients above act on the extended surface of the air side, that is the area of the fins and the tubes. Therefore, the fin efficiency must also be considered in calculating the overall thermal resistance on the outside. Gardner(1945) gives the derivation of Equation , used as a curve fit to find the fin efficiency as a function of film coefficient. This equation is based on circular fins of constant thickness. To model a coil with flat fins, an effective diameter – that of circular fins with the same fin area – is used. Equations through define variables used in Equation . The overall efficiency of the surface is shown by Equation . Note that the efficiency is found by the same equations for the wet surface using the wet surface film coefficient.

$η_{f i n} = \frac{- 2 ρ}{f a i (1 + ρ)} [\frac{I_{1} (u_{b}) K_{1} (u_{e}) - K_{1} (u_{b}) I_{1} (u_{e})}{I_{0} (u_{b}) K_{1} (u_{e}) + K_{0} (u_{b}) I_{1} (u_{e})}]$

$f a i = \frac{(D_{f i n} - D_{t u b e})}{2} \sqrt{\frac{2 f_{o}}{k_{f i n} δ_{f i n}}}$

$ρ = \frac{D_{t u b e}}{D_{f i n}}$

$u_{e} = \frac{f a i}{1 - ρ}$

$u_{b} = u_{e} ρ$

$η_{o} = 1 - (1 - η_{f i n}) \frac{A_{f i n s}}{A_{s_{t} o t a l}}$

The definition of overall outside thermal resistance is given in Equation as a function of fin efficiency and film coefficient. For a wet coil surface the resistance must be defined differently because the heat transfer equations are based on enthalpy rather than temperature differences, as shown in Equation .

$R_{o} = \frac{1}{f_{o} η_{o} A_{s, t o t a l}}$

$R_{o, w} = \frac{C_{p, a} / b b}{f_{o, w} η_{o, w} A_{s, t o t a l}}$

Equation gives the last two overall components of thermal resistance. They represent the metal tube wall and internal fouling. The fouling factor, due to deposits of dirt and corrosion of the tube inside surfaces, is assumed to be 5x10-5 m2·K/W. All components of thermal resistance are added in series to produce the overall heat transfer coefficients shown in Equations and .

$R_{m f} = \frac{δ_{t u b e}}{k_{t u b e} A_{i}} + \frac{F l}{A_{i}}$

$U A_{d r y} = \frac{A_{d r y}}{A_{s, t o t a l}} [\frac{1}{R_{i} + R_{m f} + R_{o}}]$

$U_{c} A_{w} = \frac{A_{w}}{A_{s, t o t a l}} [\frac{1 / b b}{R_{i} + R_{m f} + R_{o, w}}]$

Solution Method of Model[LINK]

The complicated equations derived above were implemented in a successive substitution solution procedure to calculate the coil performance based on the input parameters. The MODSIM implementation of a cooling coil, the TYPE12 subroutine, was the motivation for this approach; the method used there has been retained with modifications for the uncontrolled coil model. Clark (1985) contains notes about the MODSIM routine.

In the general case, the cooling coil is only partially wet. For an uncontrolled coil, Equation is used to find the water temperature at the boundary. Several simple equations in the loop adjust the boundary point until the dry surface temperature at the boundary is equal to the dewpoint of the inlet air. For the controlled coil, Equations and give two calculations of the boundary temperature, and the water flow rate and boundary position are adjusted until the two equations agree.

Special cases occur when the coil is all wet or all dry. The coil is solved as if it were all wet before the general case is attempted. If the wet surface temperatures at the coil inlet and outlet are both below the dewpoint, no further solution is required. However, to ensure a continuous solution as flow variables are changed, when the surface is all dry or when it is wet with only the dry surface equations yielding a surface temperature below the dewpoint at the water outlet, the general solution is used to calculate the unknowns. In the solution of the controlled coil the outlet air enthalpy, given some resulting dehumidification, must correspond to the enthalpy at the specified outlet air temperature.

Application of Cooling Coil Model to Heating Coils[LINK]

The implementation of detailed heating coil models in IBLAST was another important aspect of the system/plant integration. The same kind of loops exist to provide hot water to the heating coils from the boilers as exist to supply the cooling coils with chilled water from the chillers. Some simplifications can be made, however, since the enthalpy change of the air flowing over a heating coil is entirely sensible. There is no condensation in a heating coil. In order to allow heating and cooling coils to be specified using the same geometric parameters, a heating coil simulation was developed from the cooling coil model described above by eliminating the wet surface analysis.

In addition, it was concluded that, since much simpler and less computationally expensive heating coil simulations are possible, an option was provided in IBLAST for a heating coil design using only the UA value of the coil, the product of heat transfer coefficient and coil area. This model was largely based on the TYPE10 subroutine implemented in MODSIM. The equations used to model the performance of the TYPE10 heating coil are as follows:

$\begin{matrix} T_{a, o u t} & = T_{a, i n} + (T_{w, i n} - T_{a i n}) ε (\frac{min (C_{p, a} {˙ m}_{a}, C_{p, w} {˙ m}_{w})}{C_{p, a} {˙ m}_{a}}) T_{w, o u t} & = T_{w, i n} - (T_{a, o u t} - T_{a i n}) (\frac{C_{p, a} {˙ m}_{a}}{C_{p, w} {˙ m}_{w}}) \end{matrix}$

where the coil effectiveness is given by:

$ε = 1 - exp ⎛ ⎜ ⎜ ⎜ ⎜ ⎝ \frac{{exp [- (\frac{min {C_{p, a} {˙ m}_{a}, C_{p, w} {˙ m}_{w}}}{max {C_{p, a} {˙ m}_{a}, C_{p, w} {˙ m}_{w}}}) {N T U}^{0.78}] - 1}}{(\frac{min {C_{p, a} {˙ m}_{a}, C_{p, w} {˙ m}_{w}}}{max {C_{p, a} {˙ m}_{a}, C_{p, w} {˙ m}_{w}}}) {N T U}^{- .22}} ⎞ ⎟ ⎟ ⎟ ⎟ ⎠$

The parameter NTU is the number of transfer units and is defined as a function of the UA value of the coil as follows:

$N T U = \frac{U A}{m i n (C_{p, a} {˙ m}_{a}, C_{p, w} {˙ m}_{w})}$

Hot-Water-Based Air Heating Coil[LINK]

Overview[LINK]

The input object Coil:Heating:Water provides a model that uses an NTU–effectiveness model of a static heat exchanger. The model is an inlet – outlet model: given the inlet conditions and flow rates and the UA, the effectiveness is calculated using the formula for the effectiveness of a cross-flow heat exchanger with both fluid streams unmixed. The effectiveness then allows the calculation of the outlet conditions from the inlet conditions.

The inputs to the model are: (1) the current inlet temperatures and flow rates of the air and water fluid streams and (2) the UA of the coil. Note that the UA is fixed in this model and is not a function of the flow rates.

There are 2 alternative user inputs for the component: the user may input the design water volumetric flow rate and the UA directly; or the user may choose to input the more familiar design heating capacity plus design inlet & outlet temperatures and let the program calculate the design UA. These alternative user inputs are fully described in the EnergyPlus Input Output Reference document.

Model Description[LINK]

The air and water capacitance flows are defined as:

${˙ C}_{a i r} = c_{p, a i r} \cdot {˙ m}_{a i r}$

${˙ C}_{w a t e r} = c_{p, w a t e r} \cdot {˙ m}_{w a t e r}$

The minimum and maximum capacity flows are then:

${˙ C}_{min} = min ({˙ C}_{a i r}, {˙ C}_{w a t e r})$

${˙ C}_{max} = max ({˙ C}_{a i r}, {˙ C}_{w a t e r})$

The capacitance flow ratio is defined as:

$Z = {˙ C}_{min} / {˙ C}_{max}$

The number of transfer units (NTU) is:

$N T U = U A / {˙ C}_{min}$

The effectiveness is:

$ε = 1 - exp (\frac{e^{- N T U \cdot Z \cdot η} - 1}{Z \cdot η})$

Where $η = N T U^{- 0.22}$ .

The outlet conditions are then:

$T_{a i r, o u t} = T_{a i r, i n} + ε \cdot {˙ C}_{min} \cdot (T_{w a t e r, i n} - T_{a i r, i n}) / {˙ C}_{a i r}$

$T_{w a t e r, o u t} = T_{w a t e r, i n} - {˙ C}_{a i r} \cdot (T_{a i r, o u t} - T_{a i r, i n}) / {˙ C}_{w a t e r}$

The output of the coil in watts is:

${˙ Q}_{c o i l} = {˙ C}_{w a t e r} \cdot (T_{w a t e r, i n} - T_{w a t e r, o u t})$

The UA value is recalculated for each timestep. A nominal UA, UA $_{0}$ , at the rating point is calculated by the program using the input for rated conditions and a search routine called regula falsi.

User input for the ratio of convective heat transfers at the nominal or rated operating point, “r,” is used in the model. This ratio is defined as

$r = \frac{η_{f} {(h A)}_{a i r}}{{(h A)}_{w a t e r}}$

where,

$η_{f}$ is the fin efficiency, (dimensionless)

h is the surface convection heat transfer coefficient

A is the surface area

The value calculated for UA $_{0}$ is used with the input for r to characterize the convective heat transfer on the water sides at the nominal rating operation point using

${({(h A)}_{w})}_{0} = U A_{0} (\frac{r + 1}{r})$

and on the air side at the nominal rating point using

${(η_{f} {(h A)}_{a})}_{0} = r {(h A)}_{w, 0}$

Then the following equations are used to calculate a new UA as a function of the flow rates and inlet temperatures at each timestep.

$x_{a} = 1 + 4.769 \cdot 10^{- 3} (T_{a i r, i n} - T_{a i r, i n, 0})$

$η_{f} {(h A)}_{a} = x_{a} {(\frac{{˙ m}_{a}}{{˙ m}_{a, 0}})}^{0.8} {(η_{f} {(h A)}_{a})}_{0}$

$x_{w} = 1 + (\frac{0.014}{1 + 0.014 T_{w a t e r, i n, 0}}) (T_{w a t e r, i n} - T_{w a t e r, i n, 0})$

${(h A)}_{w} = x_{w} {(\frac{{˙ m}_{w}}{{˙ m}_{w, 0}})}^{0.85} {(h A)}_{w, 0}$

$U A = {(\frac{1}{{(h A)}_{w}} + \frac{1}{η_{f} {(h A)}_{a}})}^{- 1}$

The above formulas are from the following reference, along with further references. The equation for x $_{w}$ was modified from that published in Wetter (1999) to correct a small error.

References[LINK]

Wetter, M. 1999. Simulation Model: Finned Water-to-Air Coil Without Condensation. LBNL-42355. This document can be downloaded from .

Single-Speed Electric DX Air Cooling Coil[LINK]

Overview[LINK]

This model (object names Coil:Cooling:DX:SingleSpeed and
Coil:Cooling:DX:TwoStageWithHumidityControlMode, with CoilPerformance:DX:Cooling) simulates the performance of an air-cooled or evaporative-cooled direct expansion (DX) air conditioner. The model uses performance information at rated conditions along with curve fits for variations in total capacity, energy input ratio and part-load fraction to determine the performance of the unit at part-load conditions (Henderson et al. 1992, ASHRAE 1993). Sensible/latent capacity splits are determined by the rated sensible heat ratio (SHR) and the apparatus dewpoint (ADP)/bypass factor (BF) approach. This approach is analogous to the NTU-effectiveness calculations used for sensible-only heat exchanger calculations, extended to a cooling and dehumidifying coil.

This model simulates the thermal performance of the DX cooling coil and the power consumption of the outdoor condensing unit (compressor, fan, crankcase heater and evap condenser water pump). The total amount of heat rejected by the condenser is also calculated and stored for use by other waste heat recovery models (e.g., Coil:Heating:Desuperheater). The performance of the indoor supply air fan varies widely from system to system depending on control strategy (e.g., constant fan vs. AUTO fan, constant air volume vs. variable air volume, etc.), fan type, fan motor efficiency and pressure losses through the air distribution system. Therefore, this DX system model does not account for the thermal effects or electric power consumption of the indoor supply air fan. EnergyPlus contains separate models for simulating the performance of various indoor fan configurations, and these models can be easily linked with the DX system model described here to simulate the entire DX air conditioner being considered (e.g., see AirLoopHVAC:Unitary:Furnace:HeatCool, AirLoopHVAC:UnitaryHeatCool, ZoneHVAC:WindowAirConditioner or AirLoopHVAC:UnitaryHeatPump:AirToAir).

Model Description[LINK]

The user must input the total cooling capacity, sensible heat ratio (SHR), coefficient of performance (COP) and the volumetric air flow rate across the cooling coil at rated conditions. The capacity, SHR and COP inputs should be “gross” values, excluding any thermal or energy impacts due to the indoor supply air fan. The rated conditions are considered to be air entering the cooling coil at 26.7°C drybulb/19.4°C wetbulb and air entering the outdoor condenser coil at 35°C drybulb/23.9°C wetbulb. The rated volumetric air flow should be between 0.00004027 m $^{3}$ /s and 0.00006041 m $^{3}$ /s per watt of rated total cooling capacity (300 – 450 cfm/ton). The rated volumetric air flow to total cooling capacity ratio for 100% dedicated outdoor air (DOAS) application DX cooling coils should be between 0.00001677 (m3/s)/W (125 cfm/ton) and 0.00003355 (m3/s)/W (250 cfm/ton).

The user must also input five performance curves or performance tables that describe the change in total cooling capacity and efficiency at part-load conditions:

1) Total cooling capacity modifier curve or table (function of temperature)

2) Total cooling capacity modifier curve or table (function of flow fraction)

3) Energy input ratio (EIR) modifier curve or table (function of temperature)

4) Energy input ratio (EIR) modifier curve or table (function of flow fraction)

5) Part load fraction correlation curve or table (function of part load ratio)

The total cooling capacity modifier curve (function of temperature) is a curve with two independent variables: wet-bulb temperature of the air entering the cooling coil, and dry-bulb temperature of the air entering the air-cooled condenser coil (wet-bulb temperature if modeling an evaporative-cooled condenser). The output of this curve is multiplied by the rated total cooling capacity to give the total cooling capacity at the specific entering air temperatures at which the DX coil unit is operating (i.e., at temperatures different from the rating point temperatures). This curve is typically a biquadratic but any curve or table with two independent variables can be used.

Note: The data used to develop the total cooling capacity modifier curve (function of temperature) should represent performance when the cooling coil is ‘wet’ (i.e., coil providing sensible cooling and at least some dehumidification). Performance data when the cooling coil is ‘dry’ (i.e., not providing any dehumidification) should not be included when developing this modifier curve. This model automatically detects and adjusts for ‘dry coil’ conditions (see section “Dry Coil Conditions” below).

$TotCapTempModFac = Func (T_{w b, i}, T_{c}, i)$

where

$T_{w b, i}$ = x values = wet-bulb temperature of the air entering the cooling coil, °C

$T_{c, i}$ = y values = dry-bulb temperature of the air entering an air-cooled condenser or wet-bulb temperature of the air entering an evaporative-cooled condenser, °C

* The total cooling capacity modifier curve (function of flow fraction) is a curve with one independent variable being the ratio of the actual air flow rate across the cooling coil to the rated air flow rate (i.e., fraction of full load flow). The output of this curve is multiplied by the rated total cooling capacity and the total cooling capacity modifier curve (function of temperature) to give the total cooling capacity at the specific temperature and air flow conditions at which the DX unit is operating. This curve is typically a quadratic but any curve or table with one independent variable can be used.

$TotCapFlowModFac = Func (f f)$

where

$ff=flow fraction=(Actual air mass flow rateRated air mass flow rate)=x value$

Note: The actual volumetric air flow rate through the cooling coil for any simulation time step where the DX unit is operating must be between 0.00002684 m $^{3}$ /s and .00006713 m $^{3}$ /s per watt of rated total cooling capacity (200 - 500 cfm/ton). The simulation will issue a warning message if this air flow range is exceeded.

The energy input ratio (EIR) modifier curve (function of temperature) is a curve with two independent variables: wet-bulb temperature of the air entering the cooling coil, and dry-bulb temperature of the air entering the air-cooled condenser coil (wet-bulb temperature if modeling an evaporative-cooled condenser). The output of this curve is multiplied by the rated EIR (inverse of the rated COP) to give the EIR at the specific entering air temperatures at which the DX coil unit is operating (i.e., at temperatures different from the rating point temperatures). This curve is typically a biquadratic but any curve or table with two independent variables can be used.

Note: The data used to develop the energy input ratio (EIR) modifier curve (function of temperature) should represent performance when the cooling coil is ‘wet’ (i.e., coil providing sensible cooling and at least some dehumidification). Performance data when the cooling coil is ‘dry’ (i.e., not providing any dehumidification) should not be included when developing this modifier curve. This model automatically detects and adjusts for ‘dry coil’ conditions (see section “Dry Coil Conditions” below).

$EIRTempModFac = Func (T_{w b, i}, T_{c}, i)$

where

$T_{w b, i}$ = x values = wet-bulb temperature of the air entering the cooling coil, °C

$T_{c, i}$ = y values = dry-bulb temperature of the air entering an air-cooled condenser or wet-bulb temperature of the air entering an evaporative-cooled condenser, °C

The energy input ratio (EIR) modifier curve (function of flow fraction) is a curve with one independent variable being the ratio of the actual air flow rate across the cooling coil to the rated air flow rate (i.e., fraction of full load flow). The output of this curve is multiplied by the rated EIR (inverse of the rated COP) and the EIR modifier curve (function of temperature) to give the EIR at the specific temperature and air flow conditions at which the DX unit is operating. This curve is typically a quadratic but any curve or table with one independent variable can be used.

$EIRFlowModFrac = Func (f f)$

where

$ff=flow fraction=(Actual air mass flow rateRated air mass flow rate)=x value$

The part load fraction correlation (function of part load ratio) is a curve with one independent variable being part load ratio (sensible cooling load / steady-state sensible cooling capacity). The output of this curve is used in combination with the rated EIR and EIR modifier curves to give the “effective” EIR for a given simulation time step. The part load fraction (PLF) correlation accounts for efficiency losses due to compressor cycling. This curve is typically a linear, quadratic, or cubic but any curve or table with one independent variable can be used.

$PartLoadFrac = PLF = Func (PLR)$

where

$PLR = part load ratio = (\frac{actual sensible cooling load}{steady-state sensible cooling load}) = x values$

The part-load fraction correlation should be normalized to a value of 1.0 when the part load ratio equals 1.0 (i.e., no efficiency losses when the compressor(s) run continuously for the simulation time step). For PLR values between 0 and 1 (0 < = PLR < 1), the following rules apply:

$\begin{matrix} PLF & \geq 0.7 PLF & \geq PLR \end{matrix}$

If PLF < 0.7 a warning message is issued, the program resets the PLF value to 0.7, and the simulation proceeds. The runtime fraction of the coil is defined as PLR/PLF. If PLF < PLR, then a warning message is issued and the runtime fraction of the coil is limited to 1.0.

A typical part load fraction correlation for a conventional, single-speed DX cooling coil (e.g., residential or small commercial unit) would be:

PLF = 0.85 + 0.15(PLR)

All five part-load curves are accessed through EnergyPlus’ built-in performance curve equation manager (curve: quadratic, curve:cubic and curve:biquadratic). It is not imperative that the user utilize all coefficients shown in equations (449) through (453) if their performance equation has fewer terms (e.g., if the user’s PartLoadFrac performance curve is linear instead of quadratic, simply enter the values for a and b, and set coefficient c equal to zero).

For any simulation time step, the total (gross) cooling capacity of the DX unit is calculated as follows:

$∙ Q t o t a l = ∙ Q t o t a l, r a t e d (T o t C a p T e m p M o d F a c) (T o t C a p F l o w M o d F a c)$

In a similar fashion, the electrical power consumed by the DX unit (compressors plus outdoor condenser fans) for any simulation time step is calculated using the following equation:

$P o w e r = ({∙ Q}_{t o t a l}) (E I R) (R T F)$

where

$∙ Q t o t a l$ = Total cooling capacity, W – ref. equation (455)

$E I R = E n e r g y i n p u t r a t i o = (\frac{1}{C O P_{r a t e d}}) (E I R T e m p M o d F a c) (E I R F l o w M o d F a c)$

$C O P_{r a t e d}$ = Coefficient of performance at rated conditions (user input)

$RTF=({PLR}/PLRPartLoadFrac{PartLoadFrac})=runtimefractionofthecoolingcoil$

The total amount of heat rejected by the condenser is then calculated and stored for use by other waste heat recovery models (e.g., Coil:Heating:Desuperheater).

$∙ Q c o n d = ∙ Q t o t a l (1 + E I R)$

where

$∙ Q c o n d$ = total amount of heat rejected by the condenser (W)

The crankcase heater is assumed to operate when the cooling coil’s compressor is OFF and the outdoor dry-bulb temperature is below the maximum outdoor temperature for crankcase heater operation. The average crankcase heater power for the simulation time step is calculated as follows:

$P_{c r a n k c a s e} = ∙ Q c a p, c r a n k c a s e (1 - R T F)$

where

$P_{c r a n k c a s e}$ = DX cooling coil crankcase heater power, W

$∙ Q c a p, c r a n k c a s e$ = crankcase heater capacity, W

If this cooling coil is used as part of an air-to-air heat pump (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir), the crankcase heater defined for this DX cooling coil is disregarded and the associated output variable is omitted. Instead, the crankcase heater defined for the DX heating coil (Coil:Heating:DX:SingleSpeed) is enabled during the time that the compressor is not running for either heating or cooling. In this instance, RTF in the above equations would be the runtime fraction of the heat pump’s heating coil or cooling coil, whichever is greater.

In addition to calculating the total cooling capacity provided by the DX air conditioner, it is important to properly determine the break down of total cooling capacity into its sensible (temperature) and latent (dehumidification) components. The model computes the sensible/ latent split using the rated SHR and the ADP/BF approach (Carrier et al. 1959). When the DX coil model is initially called during an EnergyPlus simulation, the rated total capacity and rated SHR are used to calculate the coil bypass factor (BF) at rated conditions. The rated total capacity and rated SHR are first used to determine the ratio of change in air humidity ratio to air dry-bulb temperature:

$S l o p e R a t e d = (\frac{ω_{i n} - ω_{o u t}}{T_{d b, i n} - T_{d b, o u t}}) r a t e d$

where

$ω_{i n}$ = humidity ratio of the air entering the cooling coil at rated conditions, kg/kg

$ω_{o u t}$ = humidity ratio of the air leaving the cooling coil at rated conditions, kg/kg

$T_{d b, i n}$ = dry-bulb temperature of the air entering the cooling coil at rated conditions, °C

$T_{d b, o u t}$ = dry-bulb temperature of the air leaving the cooling coil at rated conditions, °C

Along with the rated entering air conditions, the algorithm then searches along the saturation curve of the psychrometric chart until the slope of the line between the point on the saturation curve and the inlet air conditions matches SlopeRated. Once this point, the apparatus dewpoint, is found on the saturation curve the coil bypass factor at rated conditions is calculated as follows:

$B F_{r a t e d} = \frac{h_{o u t, r a t e d} - h_{A D P}}{h_{i n, r a t e d} - h_{A D P}}$

where

$h_{o u t, r a t e d}$ = enthalpy of the air leaving the cooling coil at rated conditions, J/kg

$h_{i n, r a t e d}$ = enthalpy of the air entering the cooling coil at rated conditions, J/kg

$h_{A D P}$ = enthalpy of saturated air at the coil apparatus dewpoint, J/kg

The coil bypass factor is analogous to the “ineffectiveness” (1-ε) of a heat exchanger, and can be described in terms of the number of transfer of unit (NTU).

$B F = e^{- N T U} = e^{- (\frac{U A}{c p}) / ˙ m} = e^{- A o / ˙ m}$

For a given coil geometry, the bypass factor is only a function of air mass flow rate. The model calculates the parameter A $_{o}$ in equation (460) based on BF $_{r a t e d}$ and the rated air mass flow rate. With A $_{o}$ known, the coil BF can be determined for non-rated air flow rates.

For each simulation time step when the DX air conditioner operates to meet a cooling load, the total cooling capacity at the actual operating conditions is calculated using equation (455) and the coil bypass factor is calculated based on equation (460). The coil bypass factor is used to calculate the operating sensible heat ratio (SHR) of the cooling coil using equations (461) and (462).

$h_{A D P} = h_{i n} - \frac{({∙ Q}_{t o t a l} / ˙ m)}{1 - B F}$

$S H R = M i n i m u m ((\frac{h_{T i n, w A D P} - h_{A D P}}{h_{i n} - h_{A D P}}), 1)$

where

$h_{i n}$ = enthalpy of the air entering the cooling coil, J/kg

$h_{A D P}$ = enthalpy of air at the apparatus dewpoint condition, J/kg

$h_{T i n, w A D P}$ = enthalpy of air at the entering coil dry-bulb temperature and humidity ratio at ADP, J/kg

$˙ m$ = air mass flow rate, kg/s

With the SHR for the coil at the current operating conditions, the properties of the air leaving the cooling coil are calculated using the following equations:

$h_{o u t} = h_{i n} - \frac{{˙ Q}_{t o t a l} \cdot P L R}{˙ m}$

$h_{T i n, ω o u t} = h_{i n} - (1 - S H R) (h_{i n} - h_{o u t})$

$ω_{o u t} = P s y W F n T d b H (T_{i n}, h_{T i n, ω o u t})$

$T_{d b, o u t} = P s y T d b F n H W (h_{o u t}, ω_{o u t})$

where

$h_{o u t}$ = enthalpy of the air leaving the cooling coil, J/kg

$h_{T i n, ω o u t}$ = enthalpy of air at the entering coil dry-bulb temperature and leaving air humidity ratio, J/kg

$ω_{o u t}$ = leaving air humidity ratio, kg/kg

$T_{d b, o u t}$ = leaving air dry-bulb temperature, °C

PsyWFnTdbH = EnergyPlus psychrometric function, returns humidity ratio given dry-bulb temperature and enthalpy

PsyTdbFnHW = EnergyPlus psychrometric function, returns dry-bulb temperature given enthalpy and humidity ratio

Dry Coil Conditions[LINK]

If the model determines that the cooling coil is dry ( $ω_{i n}$ < $ω_{A D P}$ ), then equations (455) and (456) are invalid since they are functions of entering wet-bulb temperature. Under dry-coil conditions, coil performance is a function of dry-bulb temperature rather than wet-bulb temperature. In this case, the model recalculates the performance of the DX cooling unit using the calculation procedure described above but with $ω_{i n =}$ $ω_{d r y,}$ where $ω_{d r y}$ is the inlet air humidity ratio at the coil dry-out point (SHR = 1.0).

SHR Calculation Using User Specified SHR Modifier Curves[LINK]

This alternative SHR calculation method is based on user specified SHR modifying curves for temperature and flow fractions. The modifying curves correct the rated SHR value for a given DX cooling coil entering air temperatures and air mass flow fraction. These SHR modifying curves are optional additional curve inputs to the DX cooling coil objects. These two curves are a biquadratic SHR modifier curve for temperature (SHRFT), and a quadratic SHR correction curve for flow fraction (SHRFFF).

· Biquadratic SHR modifier normalized curve for DX cooling coil entering air (outdoor) wet-bulb and dry-bulb temperatures. The coil entering conditions can be outdoor air or pretreated outdoor air.

· Quadratic SHR modifier curve for flow fraction.

The SHR is given by:

$S H R = S H R_{r a t e d} \cdot S H R F T (T_{w b, i}, T_{d b, i}) \cdot S H R F F F (F F)$

The cooing coil outlet air enthalpy is given by:

$h_{o u t} = h_{i n} - \frac{{˙ Q}_{t o t a l}}{˙ m}$

The cooling coil outlet air enthalpy at the coil enlet air temperature and coil outlet humidity ratio is given by:

$h_{T_{i n} ω_{o u t}} = h_{i n} - (1.0 - S H R) \frac{{˙ Q}_{t o t a l}}{˙ m}$

The DX cooling coil outlet air humidity ratio is calculated from the psychometric function as follows:

$ω_{o u t} = P s y W F n T d b H (T_{d b, i}, h_{T_{i n} ω_{o u t}})$

The DX cooling coils leaving humidity ratio cannot exceed the DX coil entering air humidity ratio. And for dry air condition entering the cooling coil, the above equation may yield unrealistic (in some cases negative values) coil outlet humidity ratio. In this case the coil outlet air humidity is set a small value as follows:

$ω_{o u t} = M I N (0.000010, ω_{i n})$

The DX cooling coil outlet air dry-bulb temperature is determined as follows:

$ω_{o u t, C o n t F a n C y c C o m p} = ω_{o u t} (C y c l i n g R a t i o) + ω_{i n} (1 - C y c l i n g R a t i o)$

where

$S H R F T$ = sensible heat ratio modifier normalized biquadratic curve as a function of coil entering air wet-bulb and dry-bulb temperatures, (-). The DX cooling coil entering air temperature can be the outdoor air condition or pretreated outdoor air when the DX cooling coil is placed after an exhaust heat recovery heat exchangers.

$S H R F F F$ = sensible heat ratio modifier normalized quadratic curve as a function of air mass flow fraction. Flow fraction is the ratio of actual to rated mass flow rate of air through the DX cooling coil, (-).

$S H R_{r a t e d}$ = sensible heat ratio at rated condition, (-).

Condenser Options: AirCooled vs. EvaporativelyCooled[LINK]

As described previously, this model can simulate the performance of air-cooled or evaporative-cooled DX air conditioners. The following paragraphs describe three modeling options.

If the user wants to model an air-cooled condenser, they should simply specify AirCooled in the field Condenser Type. In this case, the Total Cooling Capacity Modifier Curve (function of temperature) and the Energy Input Ratio Modifier Curve (function of temperature) (equations (449) and (451) above) will utilize the outdoor dry-bulb temperature.

If the user wishes to model an evaporatively-cooled condenser AND they have performance curves that are a function of the wet-bulb temperature of air entering the condenser coil, then the user should specify Condenser Type = EvaporativlyCooled and the evaporative condenser effectiveness value should be entered as 1.0. In this case, the Total Cooling Capacity Modifier Curve (function of temperature) and the Energy Input Ratio Modifier Curve (function of temperature) (equations (449) and (451) above) will utilize the outdoor wet-bulb temperature.

If the user wishes to model an air-cooled condenser that has evaporative media placed in front of it to cool the air entering the condenser coil, then the user should specify Condenser Type = EvaporativelyCooled. The user must also enter the appropriate evaporative effectiveness for the media. In this case, the Total Cooling Capacity Modifier Curve (function of temperature) and the Energy Input Ratio Modifier Curve (function of temperature) will utilize the condenser inlet air temperature as calculated below:

$T_{c, i} = (T_{w b, o}) + (1 - E v a p C o n d E f f e c t i v e n e s s) (T_{d b, o} - T_{w b, o})$

where

$T_{c, i}$ = the temperature of the air entering the condenser coil, °C

$T_{w b, o}$ = the wet-bulb temperature of the outdoor air, °C

$T_{d b, o}$ = the dry-bulb temperature of the outdoor air, °C

In this case, the Total Cooling Capacity Modifier Curve (function of temperature) and the Energy Input Ratio Modifier Curve (function of temperature) input fields for this object should reference performance curves that are a function of outdoor dry-bulb temperature. Be aware that the evaporative media will significantly reduce the dry-bulb temperature of the air entering the condenser coil, so the Total Cooling Capacity and EIR Modifier Curves must be valid for the expected range of dry-bulb temperatures that will be entering the condenser coil.

If an evaporatively-cooled condenser is modeled, the power requirements for the water pump are calculated as follows:

$P_{e v a p c o n d p u m p} = ∙ Q c a p, e v a p c o n d p u m p (R T F)$

where

$P_{e v a p c o n d p u m p}$ = DX cooling coil evap condenser pump electric power, W

$∙ Q c a p, e v a p c o n d p u m p$ = evaporative condenser pump rated power consumption, W

Water consumption for the evaporative-cooled condenser is calculated using the difference in air humidity level across the evaporative media and the condenser air mass flow rate:

$V_{w a t e r} = \frac{{∙ m}_{_{a i r}} (ω_{e v a p c o n d, o u t} - ω_{e v a p c o n d, i n})}{ρ_{w a t e r}} (R T F) (T i m e S t e p S y s) (3600.)$

where

$V_{w a t e r}$ = DX cooling coil evap condenser water consumption, m $^{3}$

$∙ m_{a i r}$ = evaporative condenser air mass flow rate, kg/s

$ω_{e v a p c o n d, i n}$ = humidity ratio of outdoor air entering the evap condenser, kg/kg

$ω_{e v a p c o n d, o u t}$ = humidity ratio of air leaving the evap condenser, kg/kg

$ρ_{w a t e r}$ = density of water at the outdoor dry-bulb temperature, kg/m $^{3}$

TimeStepSys = HVAC system simulation time step, hr

Condenser Inlet Air Temperature[LINK]

The air temperature entering the outdoor condenser is based on the weather data used during a simulation. This temperature is either taken directly from the weather data, or can be adjusted based on the height of the outdoor condenser. The input for Condenser Air Inlet Node Name can be used to control this optional feature. If this input is left blank, the air temperature entering the condenser is based solely on the weather data. If this input is not blank, then the node name specified must also be listed in an OutdoorAir:Node or OutdoorAir:NodeList object. When the node name is listed in an OutdoorAir:NodeList object, the air temperature entering the condenser is based solely on the weather data. When the node name is listed in an OutdoorAir:Node object, the height of the node determines the air temperature entering the outdoor condenser (see description of Local Outdoor Air Temperature Calculation in the Atmospheric Variation section of this document for further details).

Supply Air Fan Control: Cycling vs. Continuous[LINK]

The DX cooling coil model simulates two specific supply air fan operation modes: cycling fan, cycling compressor or continuous fan, cycling compressor. The first operation mode is frequently referred to as “AUTO fan”, where the compressor(s) and supply air fan operate in unison to meet the zone cooling load, and cycle off together when the cooling load has been met. The second operation mode is often referred to as “fan ON”, where the compressor(s) cycle on and off to meet the zone cooling load but the supply air fan operates continuously regardless of compressor operation.

Since this DX heating coil can only be used in conjunction with a DX cooling coil (i.e. heat pumps), and these coils are used in AC equipment that specifies a fan operation mode schedule (e.g AirLoopHVAC:UnitaryHeatPump:AirToAir), the fan operation mode schedule value determines the fan operation mode for each time step throughout the simulation. A fan operation mode schedule value of 0 specifies AUTO fan mode operation while values other than 0 (a 1 is usually used) specify fan ON operation. The use of a schedule allows the fan operation mode to change based on time-of-day or with changes in season.

The EnergyPlus methodology for determining the impact that HVAC equipment has on an air stream is to calculate the mass flow rate and air properties (e.g., enthalpy, dry-bulb temperature, humidity ratio) exiting the equipment. These exiting conditions are passed along as inlet conditions to the next component model in the air stream. Eventually the flow rate and properties of the air being supplied to the conditioned zone are used in the zone energy balance to determine the resulting zone air temperature and humidity ratio. With this methodology, the determination of the air mass flow rate and air properties for the two different supply air fan operation modes is slightly different.

Cycling Fan Mode:

For the case of cycling fan/cycling compressor when humidity control is not specified, the conditions of the air leaving the cooling coil are the steady-state values calculated using equations (463), (465) and (466) above. However the air mass flow rate passed along to the next component (and eventually to the conditioned zone) is the average air mass flow rate for the system simulation time step (determined by the cooling system; see ZoneHVAC:WindowAirConditioner, AirLoopHVAC:Unitary:Furnace:HeatCool, AirLoopHVAC:UnitaryHeatCool, or AirLoopHVAC:UnitaryHeatPump:AirToAir).

For the case of cycling fan/cyling compressor when humidity control is specified, the conditions of the air leaving the cooling coil are calculated as the average conditions during the fan operating period. When the compressor operates in tandem with the fan (i.e., compressor part-load ratio [PLR] is equal to the fan PLR), the outlet conditions are identical to the calculations described above. When the compressor operates for a shorter duration than the fan (i.e., the compressor PLR is less than the heating/fan PLR), the air properties leaving the cooling coil are calculated as the average conditions during the fan operating period. In this case the calculation of exiting air conditions is analogous to the calculations for continuous fan mode described below except that PLR in the equations represents the ratio of the compressor to the fan operating period. For cycling fan systems, the fan will only operate longer than the compressor, and therefore latent degradation may be modeled (user input), when humidity control is specified, a moisture load exists (i.e., the zone air humidistat senses a moisture load), and a heating load exists where the heating PLR is greater than the cooling PLR.

Continuous Fan Mode:

For the case of continuous fan/cycling compressor, the air mass flow rate is constant. However, the air properties leaving the cooling coil are calculated as the average conditions during the system simulation time step. The model assumes that the exiting air conditions are the steady-state values calculated using equations (463), (465) and (466) above when the compressor(s) operate. For the remainder of the system simulation time step, it is assumed that the air exiting the DX coil has the same properties as the air entering the coil. For this supply air fan operating strategy, the leaving air properties are calculated as follows:

$h_{o u t, C o n t F a n C y c C o m p} = h_{o u t} (P L R) + h_{i n} (1 - P L R)$

$ω_{o u t, C o n t F a n C y c C o m p} = ω_{o u t} (P L R) + ω_{i n} (1 - P L R)$

$T_{d b, o u t, C o n t F a n C y c C o m p} = P s y T d b F n H W (h_{o u t, C o n t F a n C y c C o m p}, ω_{o u t, C o n t F a n C y c C o m p})$

Latent Capacity Degradation[LINK]

The latent (dehumidification) capacity of a direct-expansion (DX) cooling coil is strongly affected by part-load, or cyclic, operation. This is especially true in applications where the supply air fan operates continuously while the cooling coil cycles on and off to meet the cooling load. During constant fan operation, moisture condenses on the cooling coil when the compressor operates, but part or all of the moisture that is held by the coil evaporates back into the airstream when the cooling coil is deactivated (Figure). The net effect is that the amount of moisture removed from the air is degraded at part-load conditions as compared to steady-state conditions when the compressor operates continuously (Figure).

EnergyPlus is able to model latent capacity degradation based on algorithms developed by Henderson and Rengarajan (1996). The model is applicable to single-stage cooling units, like residential and small commercial air conditioners or heat pumps with less than 19 kW of nominal cooling capacity. The model inputs are described in the EnergyPlus Input/Output Reference for the object Coil:Cooling:DX:SingleSpeed. The model is enabled only if the four numerical inputs are defined (values greater than zero, see IO Reference).

The following discussion applies to both cycling fan and continuous fan operation when the fan operates for a longer period of time than the compressor and air continues to flow over the moisture laden cooling coil after compressor operation has terminated.

Transient Sensible and Latent Capacity of a Cooling Coil Over an Operating Cycle

Field Data Showing the Net Impact of Part-Load Operation on Sensible Heat Ratio

Concepts of Moisture Buildup and Evaporation

Figure graphically depicts the latent degradation concepts and defines several key model parameters. After the cooling coil starts to operate, the coil temperature is eventually reduced below the dewpoint temperature of the entering air. Moisture from the air then builds on the surface of the coil until time t $_{o}$ has elapsed and the total moisture mass on the coil is M $_{o}$ . After this time (t $_{o}$ ), moisture begins to fall from the coil and all of the latent capacity provided by the coil is “useful” since this condensate is collected and removed from the unit. When the coil cycles off and the supply air fan continues to operate, the initial moisture mass buildup on the coil (M $_{o}$ ) evaporates back into the supply air stream. If the cooling coil cycles back on before all of the moisture has evaporated, then the time until the first condensate removal (t $_{o}$ ) is shorter for this cooling cycle since the coil is already partially wetted.

Figure also shows several parameters that are used in the latent degradation model. The ratio of the coil’s moisture holding capacity (M $_{o}$ ) and the steady-state latent capacity ( $∙ Q L$ ) is defined as $t_{w e t}$ : the nominal time for moisture to fall from the coil (ignoring transient effects at startup and starting with a dry coil). The ratio of the initial moisture evaporation rate ( $∙ Q e$ ) and the steady-state latent capacity ( $∙ Q L$ ) is defined as $γ$ . Both $t_{w e t}$ and $γ$ at the rated air volume flow rate and temperature conditions are required model inputs. Two other model inputs are the Maximum ON/OFF Cycling Rate (cycles per hour, $N_{m a x}$ ) and the time constant ( $τ$ , in seconds) for the cooling coil’s latent capacity to reach steady state after startup. The development of the latent degradation model is fully described by Henderson and Rengarajan (1996). The model implemented in EnergyPlus is for their “linear decay” evaporation model.

During the simulation, all of the steady-state calculations described previously in equations (449) through (462) are completed. The latent degradation model then modifies the steady-state sensible heat ratio for the coil as shown below. The value of $t_{w e t}$ at the current air volume flow rate and entering air conditions is first calculated based on the rated value of $t_{w e t}$ entered by the user:

$t_{w e t} = M i n i m u m ⎛ ⎜ ⎝ t_{w e t, r a t e d} ⎛ ⎜ ⎝ \frac{{∙ Q}_{l a t e n t, r a t e d}}{{∙ Q}_{l a t e n t}} ⎞ ⎟ ⎠, t_{w e t, max} ⎞ ⎟ ⎠$

where

$t_{w e t}$ = nominal time for condensate removal to begin at the current airflow and entering air conditions, starting with a dry coil (sec)

$t_{w e t, r a t e d}$ = nominal time for condensate removal to begin at the coil’s rated airflow and entering air conditions, starting with a dry coil (sec)

$∙ Q l a t e n t, r a t e d$ = cooling coil latent capacity at the rated airflow and temperature conditions, W

$∙ Q l a t e n t$ = cooling coil latent capacity at the current airflow and temperature conditions, W

$t_{w e t, m a x}$ = maximum allowed value for $t_{w e t}$ (9999.0 sec)

Likewise, the value of $γ$ at the current air volume flow rate and entering air conditions is calculated based on the rated value of $γ$ entered by the user:

$γ = γ_{r a t e d} ⎛ ⎜ ⎝ \frac{{∙ Q}_{l a t e n t, r a t e d}}{{∙ Q}_{l a t e n t}} ⎞ ⎟ ⎠ (\frac{T_{d b, i} - T_{w b, i}}{T_{d b, r a t e d} - T_{w b, r a t e d}})$

where:

$γ$ = ratio of the initial moisture evaporation rate from the cooling coil (when the compressor first turns off, in Watts) and the coil’s steady-state latent capacity (Watts) at the current air volume flow rate and entering air conditions

$γ_{r a t e d}$ = $γ$ at rated air flow and entering air conditions

$T_{d b, i}$ = dry-bulb temperature of the air entering the cooling coil, °C

$T_{w b, i}$ = wet-bulb temperature of the air entering the cooling coil, °C

$T_{d b, r a t e d}$ = dry-bulb temperature of air entering the cooling coil at rated conditions (26.7°C)

$T_{w b, r a t e d}$ = wet-bulb temperature of air entering the cooing coil at rated conditions (19.4°C)

The cooling coil on and off times are then calculated based on the maximum number of cycles per hour and the calculated run-time fraction for the coil.

$t_{o n} = \frac{3600}{4 N_{max} (1 - X)}$

$t_{o f f} = \frac{3600}{4 N_{m a x} X}$

where

$t_{o n}$ = duration of cooling coil on-cycle (sec)

$N_{m a x}$ = maximum on/off cycles per hour (cph)

X = cooling coil runtime fraction (-)

$t_{o f f}$ = duration of cooling coil off-cycle (sec)

For cycling fan systems, the duration of the heating coil on and off cycles are also calculated. When the heating coil operates for a longer duration of time than the cooling coil, latent degradation can also occur. For this case, the off-cycle time (i.e., the amount of time the fan operates longer than the cooling coil) is recalculated and based on the difference between the heating coil on-cycle time and the cooling coil on-cycle time. Ton and Toff for the heating coil are calculated in the same manner as shown above except that the heating run-time fraction is used for the calculations. This model assumes that the cycling rate of the cooling and heating coils are the same. In addition, since the heating coil cycling rate may be different than the cooling coil (based on run-time fractions), the heating coil may also turn on again before the next cooling coil on-cycle. The following equations are used when the heating coil operates for a longer time period than does the cooling coil.

$t o n, h e a t i n g = \frac{3600}{4 N m a x (1 - X h e a t i n g)}$

$t o f f, h e a t i n g = \frac{3600}{4 N m a x (X h e a t i n g)}$

$t o n, h e a t i n g = t o n, h e a t i n g + M A X (0, M I N (t o n, h e a t i n g, ((t o n + t o f f) - (t o n, h e a t i n g + t o f f, h e a t i n g))))$

$t o f f = M I N (t o f f, t o n, h e a t i n g - t o n)$

where

$t_{o n, h e a t i n g}$ = duration of heating coil on-cycle (sec)

$X_{, h e a t i n g}$ = heating coil runtime fraction (-)

$t_{o f f, h e a t i n g}$ = duration of heating coil off-cycle (sec)

The equation for calculating the time $t_{o}$ when moisture first begins to fall from the cooling coil is shown below, and is solved iteratively by EnergyPlus:

$t_{o}^{j + 1} = γ t_{o f f} - (\frac{γ^{2}}{4 t_{w e t}}) t_{o f f}^{2} - τ (e^{\frac{t_{o}^{j}}{τ}} - 1), t_{o f f} \leq (\frac{2 t_{w e t}}{γ})$

where

$t_{o}$ = time where condensate removal begins (sec)

$τ$ = latent capacity time constant at start-up (sec)

j = iteration number

The part-load latent heat ratio of the cooling coil is then calculated with $t_{o}$ , t $_{o n}$ and $τ$ , which is in turn used to calculate the “effective” sensible heat ratio of the cooling including part-load latent degradation effects.

$\frac{L H R}{L H R_{s s}} = M a x i m u m ⎛ ⎜ ⎜ ⎜ ⎜ ⎝ ⎛ ⎜ ⎜ ⎜ ⎜ ⎝ \frac{t_{o n} - t_{o}}{t_{o n} + τ (e^{(\frac{- t_{o n}}{τ})} - 1)} ⎞ ⎟ ⎟ ⎟ ⎟ ⎠, 0.0 ⎞ ⎟ ⎟ ⎟ ⎟ ⎠$

$S H R_{e f f} = 1 - (1 - S H R_{s s}) (\frac{L H R}{L H R_{s s}})$

where

$L H R$ = part-load latent heat ratio

$L H R_{s s}$ = latent heat ratio at steady-state conditions ( $1 - S H R_{s s}$ with $S H R_{s s}$ from eqn. (462))

$S H R_{e f f}$ = part-load sensible heat ratio ( $S H R_{s s}$ ≤ $S H R_{e f f}$ ≤ 1.0)

$S H R_{s s}$ = steady-state sensible heat ratio (from eqn. (462))

With the “effective” SHR for the coil at the current operating conditions, including the impacts of latent degradation, equations (463) through (466) are then used to calculate the properties of the air leaving the cooling coil when it operates. Finally, equations (467) through (469) are used to calculate the average leaving air conditions (average when the coil is on and off) for the simulation time step.

Basin Heater For Single-Speed DX Coil[LINK]

Calculations are also made to estimate the electric power input to the DX coil basin heater. A schedule may be used to disable the basin heater during regular maintenance periods or other time periods (e.g., during summer). If a schedule is not provided, the basin heater is assumed to be available the entire simulation time period. The basin heater operates when it is scheduled on, the outdoor air dry-bulb temperature is below the basin heater setpoint temperature, and the DX coil is not active. The user is required to enter a basin heater capacity (watts per degree Kelvin) and a heater setpoint temperature ( $^{o}$ C) if they want to model basin heater electric power.

$P_{b a s i n h e a t e r} = max (0.0, C A P_{b a s i n h e a t e r} (T_{s e t p o i n t, b a s i n} - T_{d b, o u t d o o r}) (1 - R T F))$

$E_{b a s i n h e a t e r} = {(P}_{b a s i n h e a t e r}) (T i m e S t e p S y s) 3600$

where:

$P_{b a s i n h e a t e r}$ = DX coil basin heater electric power (W)

$E_{b a s i n h e a t e r}$ = DX coil basin heater electric consumption (J)

$T_{s e t p o i n t, b a s i n}$ = Basin heater setpoint temperature ( $^{o}$ C)

$T_{d b, o u t d o o r}$ = Outdoor air dry-bulb temperature ( $^{o}$ C)

${C A P}_{b a s i n h e a t e r}$ = Basin heater capacity (W/K)

$R T F$ = Run time fraction of DX cooling coil

Special Calculations for Coil:Cooling:DX:TwoStageWithHumidityControlMode with CoilPerformance:DX:Cooling[LINK]

The multimode DX coil uses 1 to 4 set of performance specifications which all feed into the main DX coil model as described above. The multimode DX coil may have 1 or 2 capacity stages and may have 1 enhanced dehumidification mode in addition to its “normal” or base mode performance. Any mode may specify that a portion of the total airflow is bypassed.

Bypass Operation[LINK]

When a mode has coil bypass, the non-bypassed air flow fraction is used to model the DX coil performance as described above. Then the bypassed air stream is mixed with the conditioned air stream and these conditions are place on the coil outlet node.

Multi-Stage Operation[LINK]

Multi-stage control is modeled with the DX coil. If 2-stage operation has been specified, stage 1 is assumed to operate first and continue to operate if stage 2 is required. If latent degradation is active, only one of the stages will be degraded. If stage 1 operates for less than the full time step, then latent degradation is applied. If stage 2 operates, then stage 1 is running for the full time step and latent degradation is applied only to the portion of the latent load which is attributed to stage 2 operation.

Enhanced Dehumidification Mode[LINK]

If enhanced dehumidification mode is available, this is controlled by the parent object of the DX coil, such as DXSystem:AirLoop. A dehumidification mode flag is passed to the coil model which is used to select a different set of performance specifications from a CoilPerformance:DX:Cooling input object. Then the simulation of the DX coil proceeds as described above including any bypass or multi-stage operation.

Standard Rating of Single-Speed DX Cooling Coils[LINK]

For small single-speed direct expansion (DX) cooling coils, the industry standard ratings of Standard Rating Cooling Capacity and Seasonal Energy Efficiency Ratio (SEER) are calculated according to ANSI/AHRI Standard 210/240 (AHRI 2008). These ratings apply to unitary air conditioners and air-source unitary heat pumps with air-cooled condensers with standard rating cooling capacities under 19 kW (<65,000 Btu/hr). For larger DX cooling coils, the industry standard ratings of Standard Rating Cooling Capacity, Energy Efficiency Ratio (EER), and Integrated Energy Efficiency Ratio (IEER) are calculated according to ANSI/AHRI Standard 340/360 (AHRI 2007). These ratings apply to unitary air conditioners and air-source unitary heat pumps with standard rating cooling capacities from 19 kW to below 73.2 kW (65,000 Btu/hr to <250,000 Btu/hr).

For the Coil:Cooling:DX:SingleSpeed object in EnergyPlus, these standard ratings are not direct inputs to the model. However, these standard ratings can be calculated using user-entered information for the Coil:Cooling:DX:SingleSpeed object. Since users sometimes lump the performance of several smaller DX cooling units into a single larger cooling coil object for simulation purposes, EnergyPlus outputs the Standard Rating Cooling Capacity, SEER, EER, and IEER regardless of the magnitude of the standard rating cooling capacity of the coil. It is up to the user to determine which standard ratings are applicable to the cooling coil(s) they are modeling. These standard rating values are provided in the eplusout.eio output file (Ref. OutputDetailsAndExamples.pdf) and also in the predefined tabular output reports (Output:Table:SummaryReports object, Equipment Summary). Currently, the standard ratings are only calculated and output for single-speed DX cooling coils with air-cooled condensers. If the single-speed DX coling coil is specified with an evaporatively-cooled condenser, then no standard ratings are output from EnergyPlus at this time.

Note: The standard ratings described in this section require that the DX cooling coil model be evaluated at specific operating conditions (i.e., specific wet-bulb temperatures for air entering the cooling coil and dry-bulb temperatures for air entering the air-cooled [outdoor] condenser). If the cooling coil performance curves can not be evaluated at the required test conditions, then a standard rating value calculated at the curves limit will be output and a warning message will written to eplusout.err. For example, if the curve object (Curve:Biquadratic) for Total Cooling Capacity Function of Temperature Curve has a minimum value of 21C for dry-bulb temperature entering the air-cooled condenser coil, the IEER calculation requires that EER $_{D}$ be calculated at 18.3C – so, this would result in IEER calculatd at user specified curve limit as an output and a warning message in the eplusout.err file.

The standard rating cooling capacity (AHRI 2007, AHRI 2008) is calculated as follows:

$\begin{matrix} {˙ Q}_{S t a n d a r d R a t i n g} & = {˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{R a t e d}) (T o t C a p F l o w M o d F a c_{R a t e d}) - (F a n P o w e r P e r V o l F l o w R a t e \cdot {˙ V}_{R a t e d}) \end{matrix}$

where,

${˙ Q}_{S t a n d a r d R a t i n g}$ = Standard Rating (Net) Cooling Capacity (W)

${˙ Q}_{T o t a l, R a t e d}$ = Rated Total (Gross) Cooling Capacity, user input (W)

$T o t C a p T e m p M o d F a c_{R a t e d}$ = Total Cooling Capacity Function of Temperature Curve evaluated with 19.44°C wet-bulb temperature air entering the cooling coil and 35.0°C dry-bulb temperature air entering the air-cooled (outdoor) condenser (dimensionless)

$T o t C a p F l o w M o d F a c_{R a t e d}$ = Total Cooling Capacity Function of Flow Fraction Curve evaluated at a flow fraction of 1.0 (dimensionless)

$F a n P o w e r P e r V o l F l o w R a t e$ = Rated Evaporator Fan Power Per Volume Flow Rate, user input ( W/(m $^{3}$ /s) )

${˙ V}_{R a t e d}$ = Rated Air Volume Flow Rate, user input (m $^{3}$ /s)

The Rated Evaporator Fan Power Per Volume Flow rate is a user-entered value, with a default of 773.3 W/(m $^{3}$ /s)) if the user leaves this input field blank. The default value is taken from ANSI/AHRI Standards 210/240 and 340/360 where it is defined for systems which do not have a cooling coil fan furnished as part of the system (e.g., a DX cooling coil mounted in the ductwork downstream of a gas furnace where the furnace contains the fan used for air distribution across the gas heating coil and the downstream DX cooling coil). The test conditions in ANSI/AHRI Standards 210/240 and 340/360 vary the external static pressure (i.e., pressure drop associated with ductwork and other devices external to the indoor fan/coil section) seen by the supply air fan based on the standard rating cooling capacity. Note, however, that external static pressure in actual installations is typically much higher. Further details regarding indoor fan power per volume flow rate can be found in Walker and Lutz (2005) and Walker (2007), including differences between Permanent Split Capacitor (PSC) and Brushless Permanent Magnet (BPM) fan motors. Especially at the low external static pressures defined in the ANSI/AHRI Standards, BPM motors (e.g., Electronically Commutated Motors (ECMs)) can draw significantly less power (e.g., 50-75% less) than PSC motors.

The seasonal energy efficiency ratio (SEER) is calculated as follows:

$S E E R = P L F_{0.5} \cdot E E R_{T e s t B}$

$E E R_{T e s t B} = \frac{{˙ Q}_{T o t a l, N e t, T e s t B}}{P o w e r_{T o t a l, T e s t B}}$

$\begin{matrix} {˙ Q}_{T o t a l, N e t, T e s t B} & = {˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{T e s t B}) (T o t C a p F l o w M o d F a c_{R a t e d}) - (F a n P o w e r P e r V o l F l o w R a t e \cdot {˙ V}_{R a t e d}) \end{matrix}$

$\begin{matrix} P o w e r_{T o t a l, T e s t B} & = {˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{T e s t B}) (T o t C a p F l o w M o d F a c_{R a t e d}) \cdot E I R_{T e s t B} + (F a n P o w e r P e r V o l F l o w R a t e \cdot {˙ V}_{R a t e d}) \end{matrix}$

$E I R_{T e s t B} = \frac{1}{C O P_{r a t e d}} (E I R T e m p M o d F a c_{T e s t B}) (E I R F l o w M o d F a c_{R a t e d})$

where,

$P L F_{0.5}$ = Part Load Fraction Correlation Curve evaluated at a part load ratio (PLR) of 0.5 (dimensionless)

$E E R_{T e s t B}$ = Energy efficiency ratio with 19.44°C wet-bulb temperature air entering the cooling coil, 27.78°C dry-bulb temperature air entering the air-cooled (outdoor) condenser, and rated air volume flow through the cooling coil (W/W)

${˙ Q}_{T o t a l, N e t, T e s t B}$ = Net total cooling capacity with 19.44°C wet-bulb temperature air entering the cooling coil, 27.78°C dry-bulb temperature air entering the air-cooled (outdoor) condenser, and rated air volume flow through the cooling coil (W)

$T o t C a p T e m p M o d F a c_{T e s t B}$ = Total Cooling Capacity Function of Temperature Curve evaluated with 19.44°C wet-bulb temperature air entering the cooling coil and 27.78°C dry-bulb temperature air entering the air-cooled (outdoor) condenser (dimensionless)

$P o w e r_{T o t a l, T e s t B}$ = Total electric power (compressors, condenser fans and evaporator fan) with 19.44°C wet-bulb temperature air entering the cooling coil, 27.78°C dry-bulb temperature air entering the air-cooled (outdoor) condenser, and rated air volume flow through the cooling coil (W)

$C O P_{r a t e d}$ = Coefficient of Performance at Rated Conditions, user input (W/W)

$E I R T e m p M o d F a c_{T e s t B}$ = Energy Input Ratio Function of Temperature Curve evaluated with 19.44°C wet-bulb temperature air entering the cooling coil and 27.78°C dry-bulb temperature air entering the air-cooled (outdoor) condenser (dimensionless)

$E I R F l o w M o d F a c_{R a t e d}$ = Energy Input Ratio Function of Flow Fraction Curve evaluated at a flow fraction of 1.0 (dimensionless).

Energy Efficiency Ratio (EER) is another standard rating (AHRI 2007), and it is defined as the ratio of the total cooling capacity to the total power input at any given set of rating conditions, expressed in W/W (or Btu/W-h). For this class of air-cooled DX cooling coils, EER is calculated at rated test conditions as follows:

$E E R = \frac{{˙ Q}_{S t a n d a r d R a t i n g}}{P o w e r_{T o t a l, R a t e d}}$

$\begin{matrix} P o w e r_{T o t a l, R a t e d} & = {˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{T e s t B}) (T o t C a p F l o w M o d F a c_{R a t e d}) \cdot E I R_{R a t e d} + (F a n P o w e r P e r V o l F l o w R a t e \cdot {˙ V}_{R a t e d}) \end{matrix}$

$E I R_{R a t e d} = (\frac{1}{C O P_{R a t e d}}) (E I R T e m p M o d F a c_{R a t e d}) (E I R F l o w M o d F a c_{R a t e d})$

where,

EER = Energy Efficiency Ratio (W/W)

$P o w e r_{T o t a l, R a t e d}$ = Total electric power (compressors, condenser fans and evaporator fan) with 19.44°C wet-bulb temperature air entering the cooling coil, 35.0°C dry-bulb temperature air entering the air-cooled (outdoor) condenser, and air flow rate across the evaporator at the Rated Air Volume Flow Rate (W).

$E I R T e m p M o d F a c_{R a t e d}$ = Energy Input Ratio Function of Temperature Curve evaluated with 19.44°C wet-bulb temperature air entering the cooling coil and 35.0°C dry-bulb temperature air entering the air-cooled (outdoor) condenser (dimensionless)

The Integrated Energy Efficiency Ratio (IEER) is intended to be a measure of merit for the cooling coil’s part-load efficiency. IEER replaced Integrated Part-Load Value (IPLV) as the part-load performance metric in Std. 340/360 as of January 1, 2010. Full details regarding the IEER calculation are available in ANSI/AHRI Std. 340/360 (AHRI 2007). A summary of the IEER calculations made by EnergyPlus for single-speed air-cooled DX cooling coils is provided below:

$I E E R = (0.02 \cdot E E R_{A}) + (0.617 \cdot E E R_{B}) + (0.238 \cdot E E R_{C}) + (0.125 \cdot E E R_{D})$

where,

$E E R_{A}$ = EER at 100% net capacity at AHRI standard rating conditions (same as EER calculation shown above)

$E E R_{B}$ = EER at 75% net capacity and reduced outdoor air temperature

$E E R_{C}$ = EER at 50% net capacity and reduced outdoor air temperature

$E E R_{D}$ = EER at 25% net capacity and reduced outdoor air temperature

$E E R_{B - D} = \frac{L F \cdot {˙ Q}_{T o t a l, N e t, P a r t L o a d}}{P o w e r_{T o t a l, P a r t L o a d}}$

$P o w e r_{T o t a l, P a r t L o a d} = L F (C_{D} \cdot P o w e r_{C o m p C o n d F a n, P a r t L o a d}) + (F a n P o w e r P e r V o l F l o w R a t e \cdot {˙ V}_{R a t e d})$

$\begin{matrix} P o w e r_{C o m p C o n d F a n, P a r t L o a d} = {˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{P a r t L o a d}) \cdot (T o t C a p F l o w M o d F a c_{R a t e d}) \cdot E I R_{P a r t L o a d} \end{matrix}$

$E I R_{P a r t L o a d} = (\frac{1}{C O P_{R a t e d}}) (E I R T e m p M o d F a c_{P a r t L o a d}) (E I R F l o w M o d F a c_{R a t e d})$

where,

${˙ Q}_{T o t a l, N e t, P a r t L o a d}$ = Net total cooling capacity with 19.44°C wet-bulb temperature air entering the cooling coil rated air volume flow through the cooling coil (W). The dry-bulb temperature of air entering the air-cooled condenser varies (B = 27.5°C, C = 20.0°C, D = 18.3°C).

$P o w e r_{T o t a l, P a r t L o a d}$ = Total electric power (compressors, condenser fans and evaporator fan) with 19.44°C wet-bulb temperature air entering the cooling coil and air flow rate across the evaporator at the Rated Air Volume Flow Rate (W). The dry-bulb temperature of air entering the air-cooled condenser varies (B = 27.5°C, C = 20.0°C, D = 18.3°C).

$P o w e r_{C o m p C o n d F a n, P a r t L o a d}$ = Electric power of the compressor and condenser fan at the various part-load ratios, with 19.44°C wet-bulb temperature air entering the cooling coil and rated supply air volume flow rate (W). The dry-bulb temperature of air entering the air-cooled condenser varies per the part-load ratio (B = 27.5°C, C = 20.0°C, D = 18.3°C).

$E I R T e m p M o d F a c_{P a r t L o a d}$ = Energy Input Ratio Function of Temperature Curve evaluated with 19.44°C wet-bulb temperature air entering the cooling coil and dry-bulb temperature of air entering the air-cooled condenser corresponding to the reduced part-load ratio (B = 27.5°C, C = 20.0°C, D = 18.3°C) (dimensionless).

$C_{D}$ = degradation coefficient to account for cycling of the compressor = 1.13 – 0.13LF.

$L F$ = fractional “on” time at the desired load point

The load factor (LF) is the fractional “on” time for the desired reduced load points (75%, 50%, or 25%) calculated from the following equation:

$LF=({\% Load}/%Load100{100})⋅˙QStandardRating˙QTotal,Net,PartLoad$

where,

$% L o a d$ = Part-load operating points, i.e., 75% (B), 50% (C), 25% (D)

The calculations for $Q_{T o t a l, N e t, P a r t L o a d}$ and $P o w e r_{T o t a l, P a r t L o a d}$ are calculated in nearly the same way as $Q_{T o t a l, N e t, T e s t B}$ and $P o w e r_{T o t a l, T e s t B}$ are calculated for SEER (defined above). The only difference is that these cooling capacity and power values, used for calculating $E E R_{B} / E E R_{C} / E E R_{D}$ for IEER, are calculated for a series of dry-bulb temperatures of air entering the air-cooled condenser (B = 27.5°C, C = 20.0°C, D = 18.3°C) and part-load performance degradiation correction is also applied to the condensing unit electric power calculation.

ANSI/ASHRAE 127 - Standard Ratings of Single-Speed DX Cooling Coils[LINK]

For computer and data processing room unitary air conditioners single-speed direct expansion (DX) cooling coils, the standard ratings net total cooling capacity and total cooling electric power inputs are calculated according to ANSI/AHRI Standard 127 (ASHRAE 2012). These ratings apply to unitary air conditioners with air-cooled. If the single-speed DX cooling coil is specified with an evaporatively-cooled condenser, then no standard ratings are output from EnergyPlus at this time. These standard ratings are not direct inputs to the model. However, these standard ratings can be calculated using user-entered information for the Coil:Cooling:DX:SingleSpeed object. These standard rating values are provided in the eplusout.eio output file and also in the predefined tabular output reports (Output:Table:SummaryReports object, Equipment Summary).

Note: The standard ratings described in this section require that the DX cooling coil model be evaluated at sixteen different test conditions (i.e., specific wet-bulb temperatures for air entering the cooling coil and dry-bulb temperatures for air entering the air-cooled [outdoor] condenser) for each of the four standard tests and four application classes (ASHRAE – 2012). The four test conditions: A, B, C and D are provided in the ANSI/ASHRAE Standard 127. And the test conditions are different for each application classes described in the standard. In total sixteen performance data of net cooling capacity and total electric power inputs are reported. The total cooling electric power includes the supply fan power.


	Application Class	A	B	C	D
	Class 1	23.9°C(75.0°F)	23.9°C(75.0°F)	23.9°C(75.0°F)	23.9°C(75.0°F)
Tdb, Indoor	Class 2	29.4°C(85.0°F)	29.4°C(85.0°F)	29.4°C(85.0°F)	29.4°C(85.0°F)
	Class 3	35.0°C(95.0°F)	35.0°C(95.0°F)	35.0°C(95.0°F)	35.0°C(95.0°F)
	Class 4	40.5°C(105°F)	40.5°C(105°F)	40.5°C(105°F)	40.5°C(105°F)
Tdp, Indoor		35.0°C(95.0°F)	26.7°C(80.0°F)	18.3°C(65.0°F)	4.4°C(40.0°F)
Tdb, Outdoor		35.0°C(95.0°F)	26.7°C(80.0°F)	18.3°C(65.0°F)	4.4°C(40.0°F)

The standard rating net total cooling capacity is calculated as follows:

${˙ Q}_{S t a n d a r d R a t i n g} = {˙ Q}_{T o t a l R a t e d} ({TotCapTempModFac}_{R a t e d}) ({TotCapFlowModFac}_{R a t e d}) - (FanPowerPerVolFlowRate \cdot {˙ V}_{R a t e d})$

where:

${˙ Q}_{S t a n d a r d R a t i n g}$ = Standard Rating (Net) Cooling Capacity (W)
${˙ Q}_{T o t a l, R a t e d}$ = Rated Total (Gross) Cooling Capacity, user input (W)
${TotCapTempModFac}_{R a t e d}$ = Total Cooling Capacity Function of Temperature Curve evaluated at wet-bulb temperature of air entering the cooling coil and dry-bulb temperature of air entering the air-cooled (outdoor) condenser (dimensionless)
${TotCapFlowModFac}_{R a t e d}$ = Total Cooling Capacity Function of Flow Fraction Curve evaluated at a flow fraction of 1.0 (dimensionless)
$FanPowerPerVolFlowRate$ = Rated Evaporator Fan Power Per Volume Flow Rate, user input ( W/(m3/s) )
${˙ V}_{R a t e d}$ = Rated Air Volume Flow Rate, user input (m3/s)

The standard rating net total cooling electric power input is calculated as follows:

${P o w e r}_{T o t a l} = [{˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c) ({T o t C a p F l o w M o d F a c}_{R a t e d})] \cdot E I R + (F a n P o w e r P e r V o l F l o w R a t e \cdot {˙ V}_{R a t e d})$

$\frac{E I R T e m p M o d F a c E I R F l o w M o d F a c_{R a t e d}}{C O P_{R a t e d}}$

where:

$E E R$ = Energy efficiency ratio at wet-bulb temperature of air entering the cooling coil, dry-bulb temperature of air entering the air-cooled (outdoor) condenser, and rated air volume flow through the cooling coil (W/W)
$T o t C a p T e m p M o d F a c$ = Total Cooling Capacity Function of Temperature Curve evaluated at the test condition of wet-bulb temperature of air entering the cooling coil and dry-bulb temperature of air entering the air-cooled (outdoor) condenser (dimensionless)
${P o w e r}_{T o t a l}$ = Total electric power (compressors, condenser fans and evaporator fan) at the test conditions of wet-bulb temperature of air entering the cooling coil, and dry-bulb temperature of air entering the air-cooled (outdoor) condenser, and rated air volume flow through the cooling coil (W)
$C O P_{R a t e d}$ = Coefficient of Performance at Rated Conditions, user input (W/W)
$E I R T e m p M o d F a c$ = Energy Input Ratio Function of Temperature Curve evaluated at the test condition of wet-bulb temperature of air entering the cooling coil and dry-bulb temperature of air entering the air-cooled (outdoor) condenser (dimensionless)
${E I R F l o w M o d F a c}_{R a t e d}$ = Energy Input Ratio Function of Flow Fraction Curve evaluated at a flow fraction of 1.0 (dimensionless).

Reference:

ASHRAE 2012. ANSI/ASHRAE Standard 127-2012 Method of Testing for Rating Computer and Data Processing Room Unitary Air Conditioners.

Basin Heater For Two-Stage DX Coil[LINK]

Basin heater for the object Coil:Cooling:DX:TwoStageWithHumidityControlMode operates in the same manner as for Coil:Cooling:DX:SingleSpeed. Refer to the “Basin Heater For Single-Speed DX Coil” description above.

If the number of capacity stages is equal to 1 and the CondenserType for that stage is EvapCooled, then the basin heater power is calculated for (1 - RunTimeFractionStage1) of DX coil.

If the number of capacity stages is greater than 1, then

1) If the CondenserType for stage 1 is EvapCooled, then the basin heater power is calculated for (1 - RunTimeFractionStage1) of DX coil

2) Elseif the CondenserType for stage 2 is EvapCooled, then the basin heater power is calculated for (1 - RunTimeFractionStage2) of DX coil

References[LINK]

AHRI 2008. ANSI/AHRI Standard 210/240: 2008 Standard for Performance Rating of Unitary Air-Conditioning & Air-Source Heat Pump Equipment. Arlington, VA: Air-Conditioning, Heating, and Refrigeration Institute.

AHRI 2007. ANSI/AHRI Standard 340/360: 2007 Standard for Performance Rating of Commercial and Industrial Unitary Air-Conditioning and Heat Pump Equipment. Arlington, VA: Air-Conditioning, Heating, and Refrigeration Institute.

ASHRAE. 1993. HVAC2 Toolkit: A Toolkit for Secondary HVAC System Energy Calculation. Atlanta: American Society of Heating, Refrigerating and Air-Conditioning Engineers, Inc.

Carrier, W.H., R.E. Cherne, W.A. Grant, and W.H. Roberts. 1959. Modern air conditioning, heating and ventilating, 3d ed. New York: Pitman Publishing Corporation.

Henderson, H.I. Jr., K. Rengarajan, and D.B. Shirey III. 1992. The impact of comfort control on air conditioner energy use in humid climates. ASHRAE Transactions 98(2): 104-113.

Henderson, H.I. Jr. and K. Rengarajan. 1996. A Model to Predict the Latent Capacity of Air Conditioners and Heat Pumps at Part-Load Conditions with Constant Fan Operation. ASHRAE Transactions 102(2): 266-274.

Henderson, H. 1998. The impact of part load air conditioner operation on dehumidification performance: Validating a latent capacity degradation model. Presented at ASHRAE’s IAQ & Energy ‘98 conference. New Orleans, LA. October.

Henderson, H. and D. Shirey. 2003. Understanding the Dehumidification Performance of Air-Conditioning Equipment at Part-Load Conditions. Paper presented at Joint CIBSE/ASHRAE Conference, Building Sustainability, Value and Profit. September. Edinburgh, Scotland.

Walker, I.S. and Lutz, J.D. 2005. Laboratory Evaluation of Residential Furnace Blower Performance. Berekely, CA: Lawrence Berkeley National Laboratory, Report LBNL 58752.

Walker, I.S. 2007. Comparing Residential Furnace Blowers for Rating and Installed Performance. Berkeley, CA: Lawrence Berkeley National Laboratory, Environmental Energy Technologies Division Report LBNL 62344.

Multi-Speed Electric DX Air Cooling Coil[LINK]

Overview[LINK]

This model (object name Coil:Cooling:DX:MultiSpeed) simulates the performance of an air-to-air direct expansion (DX) cooling system. The main difference compared to the other cooling coil models, such as Coil:Cooling:DX:SingleSpeed, is that this cooling coil allows modeling of two to four discrete compressor speeds. Each speed has a set of corresponding performance information at rated conditions along with curve fits for variations in total capacity, SHR, energy input ratio and part-load fraction to determine the performance of the unit at part-load conditions (DOE 1982). The full load supply airflow rate is dependent on the speed number and provided by its parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed). The part-load impact on coil energy use is automatically applied to the lowest speed. A choice is provided to determine whether part-load impacts on coil energy use are applied when the coil is operating at speeds greater than speed 1.

This model simulates the thermal performance of the indoor DX cooling coil, and the power consumption of the outdoor unit (multispeed compressor, fans, and crankcase heaters). The performance of the indoor supply air fan varies widely from system to system depending on control strategy (e.g., constant fan vs. AUTO fan), fan type, fan motor efficiency and pressure losses through the air distribution system. Therefore, this DX system model does not account for the thermal effects or electric power consumption of the indoor supply air fan. EnergyPlus contains separate models for simulating the performance of various indoor fan configurations, and these models can be easily linked with the DX system model described here to simulate the entire DX system being considered. For the time being, this coil model can only be called by the parent object AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed.

When the model determines performance at Speed 1 (the lowest speed) or cycling between OFF and Speed 1, its performance is almost the same as the performance for the Coil:Cooling:DX:SingleSpeed model. However, the outlet conditions are calculated slightly differently. Therefore, the Coil:Cooling:DX:SingleSpeed model may be considered as a subset of the model described here. When the multispeed coil model determines performance at higher speeds (above 1), the model linearly interpolates the performance at two consecutive speeds (n-1 and n) as needed to meet the cooling load, with the fraction of time at each speed established by the speed ratio.

When single mode operation is specified at higher speeds (above 1), defined in the UnitarySystemPerformance:Multispeed object, a child object of AirLoopHVAC:UnitarySystem, its performance is almost the same as the performance for the Coil:Cooling:DX:SingleSpeed model at different flow rate and capacity with given speed number. No linear interpolation is performed between two adjacent speeds. The highest speed which does not exceed the sensible load is used for each HVAC timestep. This operation does not apply to speed 1 operation and the case when the sensible load is above the highest sensible capacity.

Model Inputs[LINK]

The model inputs are also very similar to the inputs of the Coil:Cooling:DX:SingleSpeed object. The main difference is that this multispeed model requires a set of fields at each speed, such as rated capacity, rated SHR, rated COP, two capacity modifiers, two energy input ratio modifiers, part-load correction, and latent degradation inputs. The inputs also include waste heat fraction at the rated conditions and modifier as a function of temperature to calculate recoverable waste heat for heat recovery, which are not available in the similar Coil:Cooling:DX:SingleSpeed object

Speed 1 Operation[LINK]

The calculation procedures in this model, including defrost and crankcase heater, are indentical to the Coil:Heating:DX:SingleSpeed object (Ref: Coil:Heating:DX:SingleSpeed) with one exception: outlet node condition calculation when the supply air fan operates continuously (i.e., supply air fan operating mode schedule value is not equal to 0; Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed).

The following procedure provides the detailed description of the exception.

Total delivered cooling capacity

The total delivered cooling capacity for speed 1 operating at the cycling ratio needed to meet the requested cooling load is:

$Q_{c o i l, c y c l i n g} = \cdot m S p e e d 1 * C y c R a t i o * (h_{i n l e t} - h_{o u t l e t, f u l l})$

where,

Q $_{c o i l, c y c l i n g l}$ = delivered total cooling capacity for Speed 1 operating at a specific cycling ratio [W]

$\cdot m S p e e d 1$ = air mass flow rate through cooling coil at Speed 1 as set by the parent object [kg/s]

$h_{o u t l e t, f u l l}$ = specific enthalpy of the coil outlet air during full-load operation at Speed 1 (no cycling) [J/kg]

$h_{i n l e t}$ = specific enthalpy of the coil inlet air [J/kg]

CycRatio = cycling ratio at Speed 1, ratio of requested heating load to the full-load capacity of the coil at Speed 1 [dimensionless]

It is assumed that the coil provides no cooling capacity when the coil is OFF, even if the supply air fan continues to operate.

Outlet air specific enthalpy

The average specific enthalpy of the coil outlet air is then calculated based on the total delivered cooling capacity and the average air mass flow rate entering the coil:

$h_{o u t l e t, a v e r a g e} = h_{i n l e t} - \frac{Q_{c o i l, c y c l i n g}}{{\cdot m}_{i n l e t}}$

where

h $_{o u t l e t, a v e r a g e}$ = averaged specific enthalpy at the coil outlet [J/kg]

h $_{i n l e t}$ = specific enthalpy at the coil inlet [J/kg]

Q $_{c o i l, c y c l i n g}$ = total capacity at full load [W]

$\cdot m i n l e t$ = mass flow rate at the inlet to the coil as established by the parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed, Mass Flow Rate Calculation). This flow rate is the average value determined by the parent object, accounting for the specified flow rate when the cooling coil is ON and the specified flow rate when the cooling coil is OFF for the time step being simulated.

Sensible capacity

The minimum humidity ratio (HR $_{m i n}$ ) is based on humidity ratios between inlet and full load outlet as:

HR $_{m i n}$ = Minimum(HR $_{i n l e t}$ , HR $_{f u l l}$ )

where

HR $_{i n l e t}$ = Humidity ratio at the inlet [kg/kg]

HR $_{f u l l}$ = Full load humidity ratio at the outlet [kg/kg]

The coil sensible capacity may be calculated as:

$Q_{c o i l, s e n s} = \cdot m S p e e d 1 * C y c R a t i o * [h_{i n l e t} (T_{i n l e t}, H R_{min}) - h_{o u t l e t, f u l l} (T_{o u t l e t, f u l l}, H R_{min})]$

where

Q $_{c o i l, s e n s}$ = delivered sensible cooling capacity [W]

h $_{o u t l e t, f u l l}$ = full load specific enthalpy at the coil outlet as a function of outlet dry-bulb temperature at the full load, and the minimum humidity ratio [J/kg]

h $_{i n l e t}$ = specific enthalpy at the coil inlet [J/kg]

Latent capacity

The latent capacity is the difference between total and sensible capacities

$Q_{c o i l, l a t e n t} = Q_{c o i l, c y c l i n g} - Q_{c o i l, s e n s}$

where

Q $_{c o i l, l a t e n t}$ = delivered latent cooling capacity [W]

Average outlet air humidity ratio

The averaged outlet HR can be calculated as:

$H R_{o u t l e t, a v e r a g e} = H R_{i n l e t} - \frac{Q_{c o i l, l a t e n t}}{\cdot λ m i n l e t}$

where

$λ$ = heat of vaporization as a function of HR $_{m i n}$ and CycRatio * T $_{o u t l e t, f u l l}$ +(1-CycRatio) * T $_{i n l e t}$ [J/kg]

Average outlet air temperature

Using the above averaged outlet humidity ratio and specific enthalpy, the averaged outlet temperature can be calculated using the psych function of PsyTdbFnHW.

The main reason using the above approach is that outlet conditions are calculated in the same way in low and high speed operation.

The crankcase heater defined for this DX cooling coil is enabled during the time that the compressor is not running for either heating or cooling. The crankcase heater power use from either heating or cooling is reported in the heating coil (Coil:Heating:DX:MultiSpeed).

Higher Speed Operation[LINK]

This section describes how higher speed operation is simulated. When the required sensible load is less than the full load sensible capacity at Speed n (Speed Number > 1), the following calculations are performed:

Bypass factor at Speed n-1 and Speed n

$B y p a s s F a c t o r_{n} = f (R a t e d B y p a s s F a c t o r_{n}, R a t e d F l o w R a t e_{n}, A c t u a l F o w R a t e_{n})$

$B y p a s s F a c t o r_{n - 1} = f (R a t e d B y p a s s F a c t o r_{n - 1}, R a t e d F l o w R a t e_{n - 1}, A c t u a l F o w R a t e_{n - 1})$

where

BypassFactor $_{i}$ = bypass factor at actual flow rate conditions at Speed i [dimensionless]

RatedBypassFactor $_{i}$ = bypass factor at the rated conditions at Speed i [dimensionless]

RatedFowRate $_{i}$ = air mass flow rate at the rated conditions at Speed i [kg/s]

ActualFowRate $_{i}$ = actual air mass flow rate at Speed i [kg/s]

i = Speed n or Speed n-1

The bypass factor at Speed n is a function of the bypass factor at the rated conditions, rated airflow rate, and actual flow rate at Speed n. The calculation is performed by a function, called AdjustCBF in the DXCoil module.

Total capacity at Speed n-1 and Speed n

$T o t C a p_{n - 1} = f (R a t e d C a p_{n - 1}, T o t C a p T e m p M o d F a c_{n - 1}, T o t C a p F l o w M o d F a c_{n - 1}, B y p a s s F a c t o r_{n - 1})$

$T o t C a p_{n} = f (R a t e d C a p_{n}, T o t C a p T e m p M o d F a c_{n}, T o t C a p F l o w M o d F a c_{n}, B y p a s s F a c t o r_{n})$

where

TotCap $_{i}$ = total cooling capacity at given temperatures and flow rates at Speed i [w]

RatedCap $_{i}$ = cooling capacity at the rated conditions at Speed i [W]

TotCapTempModFac $_{i}$ = total cooling capacity modifier as a function of indoor web-bulb temperature and outdoor air dry-bulb temperature at Speed i

TotCapFlowModFac $_{i}$ = total cooling capacity modifier as a function of the ratio of the actual flow rate across the cooling coil to the rated airflow rate at Speed i

i = Speed n or Speed n-1

The calculation is performed by a subroutine, called CalcTotCapSHR in the DXCoil module.

EIR at Speed n-1 and Speed n

$E I R_{n - 1} = R a t e d E I R_{n - 1} * E I R T e m p M o d F a c_{n - 1} * E I R F l o w M o d F a c_{n - 1}$

$E I R_{n} = R a t e E I R_{n} * E I R T e m p M o d F a c_{n} * E I R F l o w M o d F a c_{n}$

where

EIR $_{i}$ = Energy input ratio at given temperatures and flow rates at Speed i [w]

RatedEIR $_{i}$ = Energy input ratio at the rated conditions at Speed i [W]

EIRTempModFac $_{i}$ = Energy input ratio modifier as a function of indoor and outdoor air dry-bulb temperature at Speed i

EIRFlowModFac $_{i}$ = Energy input ratio modifier as a function of ratio of the actual flow rate across the heating coil to the rated airflow rate at Speed i

i = n or n-1

Full load outlet conditions at Speed n-1 and Speed n

The calculation procedure of full load outlet conditions at Speed n-1 and Speed n is the same as the calculation procedure used in the Coil:Cooling:DX:SingleSpeed model (Ref. Coil:Cooling:DX:SingleSpeed). The difference is that the outlet conditions at Speed n-1 are calculated based on the total cooling capacity and mass flow rate at Speed n-1, while the outlet conditions at Speed n are calculated based on the total cooling capacity and mass flow rate at Speed n.

Effective total cooling capacity

$\begin{matrix} Q_{c o i l, S p e e d R a t i o} = (S p e e d R a t i o) \cdot m S p e e d n (h_{i n l e t} - h_{o u t l e t, f u l l_{S} p e e d n}) + (1 - S p e e d R a t i o) \cdot m S p e e d n - 1 (h_{i n l e t} - h_{o u t l e t, f u l l_{S} p e e d n - 1}) \end{matrix}$

where

$Q_{c o i l, S p e e d R a t i o}$ = delivered sensible cooling capacity at a given speed ratio between two consecutive speeds [W]

$\cdot m S p e e d n$ = air mass flow rate through cooling coil at Speed n as set by the parent object [kg/s]

$\cdot m S p e e d n - 1$ = air mass flow rate through cooling coil at Speed 1 as set by the parent object [kg/s]

h $_{i n l e t}$ = specific enthalpy at the coil inlet [J/kg]

h $_{o u t l e t, f u l l_S p e e d n}$ = full load specific enthalpy at the coil outlet at Speed n [J/kg]

h $_{o u t l e t, f u l l_S p e e d n - 1}$ = full load specific enthalpy at the coil outlet at Speed n-1 [J/kg]

Average outlet air specific enthalpy

$h_{o u t l e t, a v e r a g e} = h_{i n l e t} - \frac{Q_{c o i l, S p e e d R a t i o}}{{\cdot m}_{i n l e t}}$

where

h $_{o u t l e t, a v e r a g e}$ = averaged specific enthalpy at the coil outlet [J/kg]

h $_{i n l e t}$ = specific enthalpy at the coil inlet [J/kg]

$\cdot m i n l e t$ = mass flow rate at the inlet to the coil as established by the parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed, Mass Flow Rate Calculation). This flow rate is the average value determined by the parent object, accounting for the specified flow rate when the heating coil is at Speed n and the specified flow rate when the heating coil is at Speed n-1 for the time step being simulated.

Effective sensible cooling capacity

The minimum humidity ratio (HR $_{m i n}$ ) is calculated as

HR $_{m i n}$ = Minimum[HR $_{i n l e t}$ , (SpeedRatio)HR $_{f u l l, n}$ +(1.0-SpeedRatio)HR $_{f u l l, n - 1}$ )

The effective sensible cooling capacity is expressed as:

$\begin{matrix} Q_{c o i l, s e n s} = \cdot m S p e e d n (S p e e d R a t i o) [h_{i n l e t} (T_{i n l e t}, H R_{min}) - h_{o u t l e t, f u l l_{S} p e e d n} (T_{o u t l e t, n}, H R_{min})] + \cdot m S p e e d n - 1 (1 - S p e e d R a t i o) [h_{i n l e t} (T_{i n l e t}, H R_{min}) - h_{o u t l e t, f u l l_{S} p e e d n - 1} (T_{o u t l e t, n - 1}, H R_{min})] \end{matrix}$

where

Q $_{c o i l, s e n s}$ = effective sensible cooling capacity [W]

h $_{o u t l e t, f u l l_S p e e d n}$ = full load specific enthalpy at the coil outlet at Speed n as a function of outlet dry-bulb temperature at the full load, and the minimum humidity ratio [J/kg]

h $_{o u t l e t, f u l l_S p e e d n - 1}$ = full load specific enthalpy at the coil outlet at Speed n-1 as a function of outlet dry-bulb temperature at the full load, and the minimum humidity ratio [J/kg]

h $_{i n l e t}$ = specific enthalpy at the coil inlet [J/kg]

Average outlet air humidity ratio and temperature

The effective latent cooling capacity is the difference between the total and sensible capacity:

$Q_{c o i l, l a t e n t} = Q_{c o i l, S p e e d R a t i o} - Q_{c o i l, s e n s}$

Q $_{c o i l, l a t e n t}$ = effective latent cooling capacity [W]

The average outlet air HR can be calculated as:

$H R_{o u t l e t, a v e r a g e} = H R_{i n l e t} - \frac{Q_{c o i l, l a t e n t}}{{\cdot λ m}_{i n l e t}}$

where

$λ$ = heat of vaporization as a function of HR $_{m i n}$ and SpeedRatio*T $_{o u t l e t, n}$ +(1-SpeedRatio)*T $_{o u t l e t, n - 1}$ [J/kg]

At the given averaged outlet humidity ratio and specific enthalpy, the averaged outlet temperature can be calculated using the psych function of PsyTdbFnHW.

Calculate combined energy input

When the input for the field ‘Apply Part Load Fraction to Speeds Greater than 1’ is No in the object (equivalent to a single compressor), the combined energy output is calculated as follows:

$C o o l i n g P o w e r = (T o t C a p_{n}) (E I R_{n}) (S p e e d R a t i o) + (T o t C a p_{n - 1}) (E I R_{n - 1}) (1.0 - S p e e d R a t i o)$ When the input for the field ’Apply Part Load Fraction to Speeds Greater than 1’ is Yes in the object (equivalent to multiple compressors), the combined energy output is calculated as follows:

$C o o l i n g P o w e r = (T o t C a p_{n}) (E I R_{n}) (R T F) + (T o t C a p_{n - 1}) (E I R_{n - 1}) (1.0 - R T F)$

where

CoolingPower = Power used in Watt

RTF = Run time fraction at Speed n

Latent degradation

When the supply fan operation mode is ContinuousFanWithCyclingCompressorand the input of the Apply Latent Degradation to Speeds Greater than 1 is Yes, the latent degradation is included at Speed n. The calculation procedure is the same as one in the Coil:Cooling:DX:SingleSpeed object. The difference is that the rated values and run time fraction at Speed n are used. The adjusted SHR is used to calculate full load outlet conditions at Speed n.

It is expected to have less latent degradation at Speed n than Speed 1. Therefore, smaller values of the latent degradation inputs at Speed n than those at Speed 1 are recommended.

Crankcase heater

There is no power need at higher speed operation.

Higher Speed Operation with Single Mode Operation[LINK]

Waste heat calculation[LINK]

The waste heat generated by this coil object is calculated as:

$Q_{W a s t e H e a t} = (F r a c t i o n) (T e m p M o d i f i e r) (C o o l i n g P o w e)$

where

Fraction = rated waste heat fraction of the energy input

TempModifer = waste heat modifier as a function of indoor and outdoor air dry-bulb temperature

Basin Heater For Multi-Speed DX Coil[LINK]

Basin heater for the object Coil:Cooling:DX:MultiSpeed operates in the same manner as for Coil:Cooling:DX:SingleSpeed. Refer to the Basin Heater For Single-Speed DX Coil description above.

Standard Rating of Multi-Speed DX Cooling Coils[LINK]

For multi-speed direct expansion cooling coils, the industry standard ratings of Standard Rating Seasonal Energy Efficiency Ratio (SEER) are calculated according to ANSI/AHRI Standard 210/240 (AHRI 2008). These standard ratings can be calculated using the user-entered information for the Coil:Cooling:DX:MultiSpeed object. According to Standard 210/240, the Standard Rating SEER applies to air-to-air unitary heat pumps and air conditions with rated cooling capacities under 19,000 Watts (<65,000 Btu/h). The Cooling Mode Test Conditions for Units Having two-capacity Compressor standard tests A $_{2}$ , B $_{2}$ , B $_{1}$ , and F $_{1}$ are also provided in Table. Cooling Mode Test Conditions for Units Having a Two-Capacity Compressor. The equations required to calculate the net cooling capacity and SEER value are described next.

Standard Rating Cooling Capacity[LINK]

The standard net cooling capacity of multi-speed DX cooling coil is reported at the maximum speed and full load A2 test condition only. The Standard Rating (Net) Cooling Capacity calculation is the same as single speed DX cooling coil procedure when calculated at maximum speed.The procedure for calculating the Standard Rating Cooling Capacity is given by:

${˙ Q}_{S t a n d a r d R a t i n g} = {˙ Q}_{T, C R a t e d}^{k = 2} \cdot C A P F T_{T e s t A 2}^{k = 2} \cdot C A P F F F_{T e s t A 2}^{k = 2}$

${˙ Q}_{S t a n d a r d R a t i n g N e t} = {˙ Q}_{S t a n d a r d R a t i n g} - F a n_{P o w e r, C R a t e d}$

$F a n_{P o w e r, C R a t e d} = F a n P o w e r P e r V o l F l o w R a t e^{k = 2} \cdot {˙ V}_{C R a t e d}^{k = 2}$

where,

${˙ Q}_{S t a n d a r d R a t i n g N e t}$ = total standard (net) cooling capacity (W) of the air-conditioning or heat pump equipment in cooling mode determined from ANSI/AHRI Standard 210/240 and A2 test conditions shown in Table. The standard cooling test conditions for air-cooled condenser are: indoor coil entering air dry-bulb and wet-bulb temperatures 26.7°C/19.4°C and outdoor coil entering air dry-bulb and wet-bulb temperatures 35.0°C/23.9°C.

${˙ Q}_{S t a n d a r d R a t i n g}$ = total standard cooling capacity (W) of the air-conditioning or heat pump equipment determined from ANSI/AHRI Standard 210/240 and A2 test conditions shown in Table. The standard cooling test conditions for air-cooled are: indoor coil entering air dry-bulb and wet-bulb temperatures 26.7°C/19.4°C and outdoor coil entering air dry-bulb and wet-bulb temperatures 35.0°C/23.9°C.

${˙ Q}_{T, C R a t e d}^{k = 2}$ = Rated total cooling capacity at maximum speed, user input (W)

$C A P F T_{T e s t A 2}^{k = 2}$ = User-specified bi-quadratic curve evaluated at the indoor coil entering air wet-bulb temperature (19.4°C) and outdoor coil entering air dry-bulb temperature (35.0°C) for air-cooled condenser per A2 test condition as specified in Table, (dimensionless).

$C A P F F F_{T e s t A 2}^{k = 2}$ = User-specified quadratic or cubic curve modifying the total cooling capacity as function of flow fraction, (dimensionless). This curve is evaluated at a flow fraction of 1.0.

$F a n_{P o w e r, R a t e d}$ = Supply air fan power at rated conditions at high (maximum) compressor speed, (W).

${˙ V}_{_{C R a t e d}}^{k = 2}$ = Rated Air Volume Flow Rate at high (maximum) compressor speed, user input (m $^{3}$ /s)

$F a n P o w e r P e r V o l F l o w R a t e^{k = 2}$ = The Rated Indoor Coil Fan Power Per Volume Flow rate is a user-entered value, with a default of 773.3 W/(m $^{3}$ /s)) if the user leaves this input field blank. The default value is taken from ANSI/ASHRAE Standard 210/240 -2008 where it is defined for systems which do not have an indoor coil fan furnished as part of the system. See the description given below how this value is calculated. User input ( W/(m $^{3}$ /s)).

The Rated Evaporator (Indoor Coil) Fan Power Per Volume Flow rate is a user-entered value, with a default of 773.3 W/(m $^{3}$ /s)) if the user leaves this input field blank. The default value is taken from ANSI/ASHRAE Standard 210/240-2008 where it is defined for systems which do not have an Indoor Coil (Evaporator) fan furnished as part of the system. The test conditions in ANSI/AHRI Standard 210/240 vary the external static pressure (i.e., pressure drop associated with ductwork and other devices external to the indoor fan/coil section) seen by the supply air fan based on the standard rating cooling capacity.

Seasonal Energy Efficiency Ratio (SEER) for Multi-Speed DX Coil[LINK]

The SEER value for multi-speed compressor air conditioner or air-to-air heat pumps per AHRI/ANSI Std. 210/240 – 2008 is calculated as follows:

$S E E R = \frac{8 \sum j = 1 \frac{q_{c} (T_{j})}{N}}{8 \sum j = 1 \frac{e_{c} (T_{j})}{N}}$

For multi-speed compressor the SEER value is weighted average performance at different outdoor air temperature bins. The eight outdoor air temperature bins and the corresponding weight are provided in Table. Distribution of Fractional Hours with in Cooling Season Temperature Bins.

Where,

$q_{c}$ (T $_{j}$ )/N = = the ratio of space cooling capacity provided by the unit during periods of the space cooling season when the outdoor temperature fell within the range represented by bin temperature $T_{j}$ to the total number of hours in the cooling season (N), (W)

$e_{c}$ (T $_{j}$ )/N = the ratio of the electrical energy consumed by the unit during periods of the space cooling season when the outdoor temperature fell within the range represented by bin temperature Tj to the total number of hours in the cooling season (N), W.

j = the bin number, dimensionless. For cooling season j ranges from 1 to 8.

$T_{j}$ = outdoor air bin temperature, °C. Outdoor temperatures are “binned” such that calculations are only performed based one temperature within the bin. Bins of 2.8 °C with 8 cooling season bin temperatures being 19.44°C, 22.22°C, 25.0°C, 27.78°C, 30.56°C, 33.33°C, 36.11°C, 38.89°C.

The steady-state cooling capacity delivered and the electric power inputs when the DX coil is operating at minimum compressor speed (k = 1), and outdoor air temperature $T_{j}$ , are determined by linear interpolation using the values of B1 and F1 tests as follows:

${˙ Q}_{c}^{k = 1} (T_{j}) = {˙ Q}_{c}^{k = 1} (19.44) + \frac{{˙ Q}_{c}^{k = 1} (27.78) - {˙ Q}_{c}^{k = 1} (19.44)}{27.78 - 19.44} (T_{j} - 19.44)$

${˙ E}_{c}^{k = 1} (T_{j}) = {˙ E}_{c}^{k = 1} (19.44) + \frac{{˙ E}_{c}^{k = 1} (27.78) - {˙ E}_{c}^{k = 1} (19.44)}{27.78 - 19.44} (T_{j} - 19.44)$

The steady-state cooling capacities and electric power inputs at the minimum speed and test condition $B_{1}$ and $F_{1}$ are calculated from the minimum (low) speed performance curves as follows:

${˙ Q}_{c}^{k = 1} (27.78) = {˙ Q}_{T, C R a t e d}^{k = 1} \cdot C A P F T_{T e s t B 1}^{k = 1} \cdot C A P F F F_{T e s t B 1}^{k = 1} - F a n P o w e r_{C R a t e d}^{k = 1}$

${˙ E}_{c}^{k = 1} (27.78) = ⎡ ⎣ \frac{{˙ Q}_{c}^{k = 1} (27.78)}{C O P_{C R a t e d}^{k = 1}} ⎤ ⎦ \cdot E I R F T_{T e s t B 1}^{k = 1} \cdot E I R F F F_{T e s t B 1}^{k = 1} + F a n P o w e r_{C R a t e d}^{k = 1}$

${˙ Q}_{c}^{k = 1} (19.44) = {˙ Q}_{T, C R a t e d}^{k = 1} \cdot C A P F T_{T e s t F 1}^{k = 1} \cdot C A P F F F_{T e s t F 1}^{k = 1} - F a n P o w e r_{C R a t e d}^{k = 1}$

${˙ E}_{c}^{k = 1} (19.44) = ⎡ ⎣ \frac{{˙ Q}_{c}^{k = 1} (19.44)}{C O P_{C R a t e d}^{k = 1}} ⎤ ⎦ \cdot E I R F T_{T e s t F 1}^{k = 1} \cdot E I R F F F_{T e s t F 1}^{k = 1} + F a n P o w e r_{C R a t e d}^{k = 1}$

$F a n P o w e r_{C R a t e d}^{k = 1} = F a n P o w e r P e r V o l F l o w R a t e^{k = 1} \cdot {˙ V}_{C R a t e d}^{k = 1}$

Where,

${˙ Q}_{T, C R a t e d}^{k = 1}$ = rated total cooling capacity at minimum compressor speed specified by users, (W)

$C O P_{C R a t e d}^{k = 1}$ = rated gross COP at minimum compressor speed specified by users, (-)

$C A P F T_{T e s t B 1}^{k = 1}$ = cooling capacities modifier curve for temperature at minimum compressor speed and B1 test condition, (-)

$C A P F F F_{T e s t B 1}^{k = 1}$ = cooling capacities modifier curve for flow fraction at minimum compressor speed and B1 test condition, (-)

$E I R F T_{T e s t B 1}^{k = 1}$ = EIR modifier curve for temperature at minimum compressor speed and B1 test condition, (-)

$E I R F F F_{T e s t B 1}^{k = 1}$ = EIR modifier curve for flow fraction at minimum compressor speed and B1 test condition, (-)

$C A P F T_{T e s t F 1}^{k = 1}$ = cooling capacities modifier curve for temperature at minimum compressor speed and F1 test condition, (-)

$C A P F F F_{T e s t F 1}^{k = 1}$ = cooling capacities modifier curve for flow fraction at minimum compressor speed and F1 test condition, (-)

$E I R F T_{T e s t F 1}^{k = 1}$ = EIR modifier curve for temperature at minimum compressor speed and F1 test condition, (-)

$E I R F F F_{T e s t F 1}^{k = 1}$ = EIR modifier curve for flow fraction at minimum compressor speed and F1 test condition, (-)

$F a n P o w e r_{C R a t e d}^{k = 1}$ = the rated supply air fan power when the unit is operating at minimum compressor speed, (W)

$V_{C R a t e d}^{k = 1}$ = rated cooling supply air volume flow rate at minimum compressor speed specified by users, (-)

$F a n P o w e r P e r V o l F l o w R a t e^{k = 1}$ = the rated Indoor Coil fan power per volume flow rate at low (minimum) compressor speed specified value by the user, (W/(m3/s))

The steady-state cooling capacity delivered and the electric power inputs when the DX cooling coil is operating at maximum (high) compressor speed (k = 2), and outdoor air temperature $T_{j}$ , are determined as follows:

${˙ Q}_{c}^{k = 2} (T_{j}) = {˙ Q}_{c}^{k = 2} (27.78) + \frac{{˙ Q}_{c}^{k = 2} (35.0) - {˙ Q}_{c}^{k = 2} (27.78)}{35.0 - 27.78} (T_{j} - 27.78)$

${˙ E}_{c}^{k = 2} (T_{j}) = {˙ E}_{c}^{k = 2} (27.78) + \frac{{˙ E}_{c}^{k = 2} (35.0) - {˙ E}_{c}^{k = 2} (27.7)}{35.0 - 27.78} (T_{j} - 27.78)$

The steady-state cooling capacities and electric power input at the maximum speed and test condition $A_{2}$ and $B_{2}$ are determined from the maximum (high) speed performance curves as follows:

${˙ Q}_{c}^{k = 2} (35.0) = {˙ Q}_{T, C R a t e d}^{k = 2} \cdot C A P F T_{T e s t A 2}^{k = 2} \cdot C A P F F F_{T e s t A 2}^{k = 2} - F a n P o w e r_{C R a t e d}^{k = 2}$

${˙ E}_{c}^{k = 2} (35.0) = ⎡ ⎣ \frac{{˙ Q}_{c}^{k = 2} (35.0)}{C O P_{C R a t e d}^{k = 2}} ⎤ ⎦ \cdot E I R F T_{T e s t A 2}^{k = 2} \cdot E I R F F F_{T e s t A 2}^{k = 2} + F a n P o w e r_{C R a t e d}^{k = 2}$

${˙ Q}_{c}^{k = 2} (27.78) = {˙ Q}_{T, C R a t e d}^{k = 2} \cdot C A P F T_{T e s t B 2}^{k = 2} \cdot C A P F F F_{T e s t B 2}^{k = 2} - F a n P o w e r_{C R a t e d}^{k = 2}$

${˙ E}_{c}^{k = 2} (27.78) = ⎡ ⎣ \frac{{˙ Q}_{c}^{k = 2} (19.44)}{C O P_{C R a t e d}^{k = 2}} ⎤ ⎦ \cdot E I R F T_{T e s t B 2}^{k = 2} \cdot E I R F F F_{T e s t B 2}^{k = 2} + F a n P o w e r_{C R a t e d}^{k = 2}$

$F a n P o w e r_{C R a t e d}^{k = 2} = F a n P o w e r P e r V o l F l o w R a t e^{k = 2} \cdot {˙ V}_{C R a t e d}^{k = 2}$

Where,

${˙ Q}_{T, C R a t e d}^{k = 2}$ = rated total cooling capacity at maximum (high) compressor speed specified by users, (W)

$C O P_{C R a t e d}^{k = 2}$ = rated gross COP at maximum (high) compressor speed specified by users, (-)

$C A P F T_{T e s t A 2}^{k = 2}$ = cooling capacity modifier curve for temperature at maximum (high) compressor speed and A $_{2}$ test condition, (-)

$C A P F F F_{T e s t A 2}^{k = 2}$ = cooling capacity modifier curve for flow fraction at maximum (high) compressor speed and A $_{2}$ test condition, (-)

$E I R F T_{T e s t A 2}^{k = 2}$ = EIR modifier curve for temperature at maximum (high) compressor speed and A $_{2}$ test condition, (-)

$E I R F F F_{T e s t A 2}^{k = 2}$ = EIR modifier curve for flow fraction at maximum (high) compressor speed and A $_{2}$ test condition, (-)

$C A P F T_{T e s t B 2}^{k = 2}$ = cooling capacity modifier curve for temperature at maximum (high) compressor speed and B $_{2}$ test condition, (-)

$C A P F F F_{T e s t B 2}^{k = 2}$ = cooling capacity modifier curve for flow fraction at maximum (high) compressor speed and B $_{2}$ test condition, (-)

$E I R F T_{T e s t B 2}^{k = 2}$ = EIR modifier curve for temperature at maximum (high) compressor speed and B $_{2}$ test condition, (-)

$E I R F F F_{T e s t B 2}^{k = 2}$ = EIR modifier curve for flow fraction at maximum (high) compressor speed and B $_{2}$ test condition, (-)

$F a n P o w e r_{C R a t e d}^{k = 2}$ = the rated supply air fan power when the unit is operating at maximum (high) compressor speed, (W)

$V_{C R a t e d}^{k = 2}$ = rated supply air volume flow rate at maximum (high) compressor speed specified by users, (-)

$F a n P o w e r P e r V o l F l o w R a t e k^{k = 2}$ = the Rated Evaporator (Indoor Coil) Fan Power Per Volume Flow rate at maximum (high) compressor speed specified value by the user, (W/(m3/s))

The above steps show how the cooling capacity and electric power inputs are determined when the DX cooling coil is operating at minimum (low) and maximum (high) compressor speeds. But the unit may operate at minimum (low) speed capacity, cycle on–off, cycle between successiave lower and higher compressor speed capacity, or operate at maximum (high) speed capacity depending on the building cooling load. The operating range of the DX cooling coil is determined based on the building cooling load for each binned outside air temperature. The building cooling load at an outdoor air temperature $T_{j}$ , is calculated as follows:

$B L (T_{j}) = (\frac{T_{j} - 18.3}{35.0 - 18.3}) ⎛ ⎝ \frac{{˙ Q}_{c}^{k = 2} (35.0)}{1.1} ⎞ ⎠$

Where,

BL(T $_{j}$ ) = the building space cooling load corresponding to outdoor temperature of $T_{j}$ , (W).

${˙ Q}_{c}^{k = 2} (35.0)$ = the cooling capacity determined from the standard A2 test, (W)

The temperatures 35.0 °C and 18.3 °C in the building load calculation equation represent the outdoor design air temperature, and zero-load base temperature, respectively (ANSI/ASHRAE, 2008). 1.1 is a sizing factor.

The cooling capacity delivered and the electric power inputs calculations when the DX cooling coil is cycling on-off, operating at minimum (low) compressor speed, cycling between successive minimum (low) and maximum (high) compressor speed, or operating at maximum (high) compressor speed are described next.

Case 1: The steady state cooling capacity when the unit is operating at or below the minimum (low) speed compressor capacity, i.e., when the building cooling load is less or equal to the minimum (low) compressor speed capacity, is calculated as follows:

${˙ Q}_{c}^{k = 1} (T_{j}) \geq B L (T_{j})$

$\frac{q_{c} (T_{j})}{N} = [X^{k = 1} (T_{j}) \cdot {˙ Q}_{c}^{k = 1} (T_{j})] \cdot \frac{n_{j}}{N}$

$\frac{e_{c} (T_{j})}{N} = [\frac{X^{k = 1} (T_{j})}{P L F_{j}} {˙ E}_{c}^{k = 1} (T_{j})] \cdot \frac{n_{j}}{N}$

$X^{k = 1} (T_{j}) = \frac{B L (T_{j})}{{˙ Q}_{c}^{k = 1} (T_{j})}$

$P L F_{j} = 1 - C_{D}^{c} [1 - X^{k = 1} (T_{j})]$

Where,

X(T $_{j}$ ) = the cooling mode load factor or part-load ratio for temperature bin j, (-)

n $_{i}$ /N = fractional bin hours for the cooling season; the ratio of the number of hours during the cooling season when the outdoor temperature fell within the range represented by bin temperature T $_{j}$ to the total number of hours in the cooling season, (-). (see Table)

$C_{D}^{c}$ = cooling coefficient of degradation, default value is 0.25.

Case 2: The unit cycles between successive the minimum (low) and maximum (high) compressor speed capacity to meet the building cooling load at outdoor air temperature T $_{j}$ . That is, the cooling building load is between the units successive minimum (low) and maximum (high) compressor speed capacities:

${˙ Q}_{c}^{k = 1} (T_{j}) < B L (T_{j}) < {˙ Q}_{c}^{k = 2} (T_{j})$

$\frac{q_{c} (T_{j})}{N} = [X^{k = 1} (T_{j}) \cdot {˙ Q}_{c}^{k = 1} (T_{j}) + (1 - X^{k = 1} (T_{j})) \cdot {˙ Q}_{c}^{k = 2} (T_{j})] \cdot \frac{n_{j}}{N}$

$\frac{e_{c} (T_{j})}{N} = [X^{k = 1} (T_{j}) \cdot {˙ E}_{c}^{k = 1} (T_{j}) + (1 - X^{k = 1} (T_{j})) \cdot {˙ E}_{c}^{k = 2} (T_{j})] \cdot \frac{n_{j}}{N}$

$X^{k = 1} (T_{j}) = \frac{{˙ Q}_{c}^{k = 2} (T_{j}) - B L (T_{j})}{{˙ Q}_{c}^{k = 2} (T_{j}) - {˙ Q}_{c}^{k = 1} (T_{j})}$

Case 3: The steady-state cooling capacity when the unit is operating continuously at maximum (high) compressor speed capacity at outdoor air temperature $T_{j}$ . That is the building cooling load is greater or equal to the available capacity at maximum (high) compressor speed:

$B L (T_{j}) \geq {˙ Q}_{c}^{k = 2} (T_{j})$

For units when operating continuously at maximum (high) compressor speed (k = 2) at temperature $T_{j}$ , the delivered cooling capacity and electric power inputs are calculated as follows:

$\frac{q_{c} (T_{j})}{N} = [{˙ Q}_{c}^{k = 2} (T_{j})] \cdot \frac{n_{j}}{N}$

$\frac{e_{c} (T_{j})}{N} = [{˙ E}_{c}^{k = 2} (T_{j})] \cdot \frac{n_{j}}{N}$

Cooling Mode Test Conditions for Units Having a Two-Capacity Compressor
Test
(r)2-3 (l)4-5 Description	Dry Bulb	Wet Bulb	Dry Bulb	Wet Bulb	Compressor Capacity	Cooling Air Volume Rate
A2 Test-required (steady, wet coil)	80.0/26.7	67.0/19.4	95.0/35.0	75.0/23.9	High	Cooling Full-load
B2 Test-required (steady, wet coil)	80.0/26.7	67.0/19.4	82.0/27.8	65.0/18.3	High	Cooling Full-load
B1 Test-required (steady, wet coil)	80.0/26.7	67.0/19.4	82.0/27.8	65.0/18.3	Low	Cooling Minimum
C2 Test-optional (steady, dry coil)	80.0/26.7		82.0/27.8		High	Cooling Full-load
D2 Test-required (cyclic, dry coil)	80.0/26.7		82.0/27.8		High
C1 Test-optional (steady, dry coil)	80.0/26.7		82.0/27.8		Low	Cooling Minimum
D1 Test-optional (cyclic, dry coil)	80.0/26.7		82.0/27.8		Low
F1 Test-optional (steady, dry coil)	80.0/26.7	67.0/19.4	67.0/19.4	53.5/11.9	High	Cooling Minimum

Distribution of Fractional Hours with in Cooling Season Temperature Bins
Source: AHRI Standard 210-240, 2008; Table 16, Page 94
Bin Number j	Bin Temperature Range (C)	Representative Temperature for bin (C)	Fraction of Total Temperature Bin Hours Nj/N
1	Dry Bulb	Wet Bulb	Dry Bulb
2	18.33 - 20.56	19.44	0.214
3	21.11 - 23.33	22.22	0.231
4	23.89 - 26.11	25	0.216
5	26.67 - 28.89	27.78	0.161
6	29.44 - 31.67	30.56	0.104
7	32.22 - 34.44	33.33	0.052
8	35.00 - 37.22	36.11	0.018

References[LINK]

See the references for the single speed DX cooling coil earlier in this document.

Two-Speed Electric DX Air Cooling Coil[LINK]

Overview[LINK]

The input object Coil:Cooling:DX:TwoSpeed provides a model that is meant to represent either a 2 speed (fan and compressor) DX cooling unit (compressor, evaporator, condenser) or a variable speed DX cooling unit in which the variation of compressor speed and air flow rate is nearly linear as a function of cooling part load. In EnergyPlus Coil:Cooling:DX:TwoSpeed is used as the cooling coil in packaged variable volume systems (PVAV). The model is based upon the single speed DX unit Coil:Cooling:DX:SingleSpeed. Basically the model for Coil:Cooling:DX:TwoSpeed employs a separate single speed model at high speed (full load) and low speed (minimum load) and interpolates between these 2 states to obtain the needed cooling output. Below minimum load the unit cycles on/off, just like the single speed unit.

Inputs and Data[LINK]

The input is similar to that for Coil:Cooling:DX:SingleSpeed except there needs to be two complete sets of performance data. The user inputs a rated total cooling capacity, rated SHR, rated COP, and rated air volumetric flow rate for the high speed and low speed states. Performance curves – cooling capacity as a function of entering air wet-bulb temperature and outside dry-bulb temperature (wet-bulb if there is an evaporatively cooled condenser), EIR as a function of the same two temperatures – must be defined and referenced for both high and low speed states. The performance characteristics of the evaporative condenser, if present, also need to be given at high and low speed states: effectiveness, air volumetric flow rate, and pump power consumption. The full list of inputs is given in the Input/Output Reference document.

The data for Coil:Cooling:DX:TwoSpeed is stored in the array structure DXCoil in the program module DXCoils

Calculation[LINK]

Calculation is performed in routine CalcMultiSpeedDXCoil in module DXCoils. The inputs to the calculation are the coil index (denoting a specific coil with its particular performance parameters), the speed ratio, and the cycling ratio. The speed ratio is an artificial parameter between 0 and 1. If the speed ratio (SR) is greater than zero, the performance of the unit will be:

$S R \cdot H i g h S p e e d P e r f o r m a n c e + (1 - S R) \cdot L o w S p e e d P e r f o r m a n c e$

Here HighSpeedPerformance means the electricity consumption, cooling output, and outlet conditions of the unit if the unit were operating at high speed (full load). Similarly LowSpeed Performance means the electricity consumption, cooling output, and outlet conditions if the unit were operating at low speed (minimum non-cycling load). The calculations for each state are the same as for the single speed, cycling DX unit Coil:Cooling:DX:SingleSpeed.

If the speed ratio is zero and the cycling ratio is > 0, the unit will be in cycling mode. The unit will be on for cycling ratio fraction of the time step, off for the remainder of the time step. While on, the unit will perform according to the low speed performance parameters.

Simulation and Control[LINK]

Coil:Cooling:DX:TwoSpeed is not used by itself; it is used as part of an encompassing component or system that provides control for the unit. In setting up a PVAV system, for instance, Coil:Cooling:DX:TwoSpeed is part of the component CoilSystem:Cooling:DX, which controls Coil:Cooling:DX:TwoSpeed. When used in the CoilSystem:Cooling:DX component, the simulation allows for either temperature or temperature and dehumidification control based on a dehumidification control selection in the CoilSystem:Cooling:DX component. If None is selected the simulation runs the DX unit to satisfy a supply air temperature setpoint. It typically varies the speed ratio in an iterative solution process to establish the operating point that will give the desired unit outlet air temperature. If CoolReheat is selected, the simulation first runs the DX unit to satisfy a supply air temperature setpoint. However, in this case if the DX unit outlet air humidity ratio is above the humidity ratio setpoint, the simulation will increase the speed and/or cycling ratio to meet the desired unit outlet air humidity ratio. The humidity ratio setpoint is specified through the use of a humidistat (ref. ZoneControl:Humidistat) and a setpoint manager (ref. SetpointManager:SingleZone:Humidity:Maximum, SetpointManager:MultiZone:MaximumHumidity:Average, SetpointManager:MultiZone:Humidity:Maximum or SetpointManager:OudoorAirPretreat). When used with the CoilSystem:Cooling:DX component, the dehumidification control type Multimode is not allowed.

Basin Heater For Two-Speed DX Coil[LINK]

Basin heater for the object Coil:Cooling:DX:TwoSpeed operates in the same manner as for Coil:Cooling:DX:SingleSpeed. Refer to the “Basin Heater For Single-Speed DX Coil” description above.

Standard Ratings[LINK]

The Coil:Cooling:DX:TwoSpeed object includes industry standard ratings of Standard Rating Cooling Capacity, Energy Efficiency Ratio (EER), and Integrated Energy Efficiency Ratio (IEER) which are calculated according to ANSI/AHRI Standard 340/360-2007 with Addenda 1 and 2 (AHRI 2007). These ratings apply to unitary air conditioners with air-cooled condensers with standard rated cooling capacities from 19 kW (65,000 Btu/h) to below 220 kW (760,000 Btu/h).

For the Coil:Cooling:DX:TwoSpeed object in EnergyPlus, these standard ratings are not direct inputs to the model. However, these standard ratings can be calculated using user-entered information for the Coil:Cooling:DX:TwoSpeed object. These standard rating values are provided in the eplusout.eio output file (Ref. OutputDetailsAndExamples.pdf) and also in the predefined tabular output reports (Output:Table:SummaryReports object, Equipment Summary report, DX Cooling Coil table and VAV DX Cooling Standard Rating Details table). Currently, the standard ratings are only calculated and output for Packaged VAV cooling coils with air-cooled condensers and variable-air-volume fans. If the two-speed DX coling coil is specified with an evaporatively-cooled condenser or along with a constant volume fan, then no standard ratings are output from EnergyPlus at this time.

The standard rating cooling capacity is calculated as follows:

$\begin{matrix} {˙ Q}_{S t a n d a r d R a t i n g} = {˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{R a t e d}) (T o t C a p F l o w M o d F a c_{R a t e d}) - (F a n H e a t C o r r e c t i o n) \end{matrix}$

where,

${˙ Q}_{S t a n d a r d R a t i n g}$ = Standard Rating (Net) Cooling Capacity (W)

${˙ Q}_{T o t a l, R a t e d}$ = Rated Total (Gross) Cooling Capacity, user input (W)

$T o t C a p T e m p M o d F a c_{R a t e d}$ = Total Cooling Capacity Function of Temperature Curve evaluated with 19.4°C wet-bulb temperature air entering the cooling coil and 35.0°C dry-bulb temperature air entering the air-cooled (outdoor) condenser (dimensionless)

$T o t C a p F l o w M o d F a c_{R a t e d}$ = Total Cooling Capacity Function of Flow Fraction Curve evaluated at a flow fraction of 1.0 (dimensionless)

$F a n H e a t C o r r e c t i o n$ = The rate that air is heated by the supply fan ( W )

$F a n H e a t C o r r e c t i o n$ is calculated in one of two ways depending on the user’s input for the field called Unit Internal Static Air Pressure. If this field was left blank then the fan heat is calculated using

$F a n H e a t C o r r e c t i o n = {˙ V}_{R a t e d H i g h S p e e d} * 773.3$

where

${˙ V}_{R a t e d H i g h S p e e d}$ = Rated Air Volume Flow Rate at high speed, user input (m $^{3}$ /s)

The value of 773.3 W/(m $^{3}$ /s) (365 W/1000cfm) is specified by ANSI/AHRI 340/360 for indoor-coil-only units that are not provided with a supply fan.

If an internal static pressure is provided, then the fan heat is modeled by evaluating the full model for Fan:VariableVolume. This is different that the ratings calculated for single-speed DX equipment which does not use data from an associated fan object. The program detects the VAV fan associated with the two-speed coil and uses the input data for that fan but with a different total air pressure rise. The total pressure is the sum of the internal pressure and the external pressure (i.e., pressure drop associated with ductwork and other devices external to the indoor fan/coil section). The test conditions in ANSI/AHRI Standard 340/360 vary the (minimum) external static pressure seen by the supply air fan based on the standard rating cooling capacity (see Table # 5 in ANSI/AHRI Standard 340/360). For the part load test points the supply air flow is reduced and the external pressure varies with the square of the flow ratio (see Note 1 in Table # 6 in ANSI/AHRI Standard 340/360-2007 with Addenda 1 and 2). The entire VAV fan model in EnergyPlus is evaluated at the specified air flow rate and the total static pressure for rating and the fan heat is calculated using

$F a n H e a t C o r r e c t i o n = E n t h a l p y_{l e a v i n g F a n} - E n t h a l p y_{e n t e r i n g F a n}$

$E E R = \frac{{˙ Q}_{S t a n d a r d R a t i n g}}{P o w e r_{T o t a l, R a t e d}}$

$\begin{matrix} P o w e r_{T o t a l, R a t e d} = ({˙ Q}_{T o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{R a t e d}) (T o t C a p F l o w M o d F a c_{R a t e d})) \cdot E I R_{R a t e d} + (F a n P o w e r C o r r e c t i o n) \end{matrix}$

$E I R_{R a t e d} = (\frac{1}{C O P_{R a t e d}}) (E I R T e m p M o d F a c_{R a t e d}) (E I R F l o w M o d F a c_{R a t e d})$

where,

EER = Energy Efficiency Ratio (W/W)

$P o w e r_{T o t a l, R a t e d}$ = Total electric power (compressors, condenser fans and evaporator fan) with 19.4°C wet-bulb temperature air entering the cooling coil, 35.0°C dry-bulb temperature air entering the air-cooled (outdoor) condenser, and air flow rate across the evaporator at the rated high speed air flow rate. Similar to the fan heat correction, the fan power correction is calculated in one of two ways depending on the input for internal static pressure. If no internal static is provided, fan power is calculated using:

$F a n P o w e r C o r r e c t i o n = {˙ V}_{R a t e d H i g h S p e e d} * 773.3$

If the internal static pressure is provided, then the entire VAV fan model in EnergyPlus is evaluated at the rated high speed air flow rate and the total static pressure for rating and the then the fan power calculated by the fan model is used for the fan power correction.

The Integrated Energy Efficiency Ratio (IEER) is intended to be a measure of merit for the cooling coil’s part-load efficiency. IEER replaced Integrated Part-Load Value (IPLV) as the part-load performance metric in Std. 340/360 as of January 1, 2010. Full details regarding the IEER calculation are available in ANSI/AHRI Std. 340/360 (AHRI 2007). A summary of the IEER calculations made by EnergyPlus for two-speed air-cooled DX cooling coils is provided below:

$I E E R = (0.02 \cdot E E R_{A}) + (0.617 \cdot E E R_{B}) + (0.238 \cdot E E R_{C}) + (0.125 \cdot E E R_{D})$

where,

$E E R_{A}$ = EER at 100% net capacity at AHRI standard rating conditions (same as EER calculation shown above but in Btu/h)

$E E R_{B}$ = EER at 75% net capacity and reduced outdoor air temperature Btu/h

$E E R_{C}$ = EER at 50% net capacity and reduced outdoor air temperature Btu/h

$E E R_{D}$ = EER at 25% net capacity and reduced outdoor air temperature Btu/h

The following table summarizes the test points.

Test points for Two-Speed DX Coil Ratings
Point	Net Capacity	Coil inlet wetbulb	Coil inlet drybulb	Condenser inlet drybulb	Supply air flow rate
A	100%	19.4C	26.7C	35.0C	Rated High speed flow rate
B	75%	19.4C	26.7C	27.5C	Regula falsi result where coil leaving drybulb matches Point A
C	50%	19.4C	26.7C	20.0C
D	25%	19.4C	26.7C	18.3C

The units are assumed to be VAV and Standard 340/360 is be applied as for VAV indoor supply fan. Because the standard stipulates “the airflow rate at part load should be adjusted to maintain the full load measured leaving dry-bulb temperature,” the part load rating test points “B”, “C”, and “D” are evaluated by using the entire DX coil model calculations to obtain the supply air conditions leaving the coil and iterating on supply air flow rate to find the supply air flow rate. The numerical method called Regula Falsi is used to find the supply flow rate at each of the part load test points. Once the supply air flow rate is known, the two-speed DX coil model results are used to determine EER at the part load test points.

When evaluating the two-speed DX coil model, the speed ratio and cycling ratio are specified based on the desired or target net capacity. When the low speed net capacity is lower than the target part load net capacity, the speed ratio is calculated using

$S p e e d R a t i o = \frac{D e s i r e d N e t C a p a c i t y - L o w S p e e d N e t C a p a c i t y}{H i g h S p e e d N e t C a p a c i t y - L o w S p e e d N e t C a p a c i t y}$

When the low speed net capacity is higher than a target part load net capacity then the unit must cycle to meet the lower load. The speed ratio is then set to 0.0 (minimum unloading) and the cycling ratio is calculated using

$C y c l i n g R a t i o = \frac{D e s i r e d N e t C a p a c i t y}{L o w S p e e d N e t C a p a c i t y}$

The results for EER, COP, and supply air mass flow rate at the part load test points are reported to the EIO file and a predefined table report called VAV DX Cooling Standard Rating Details.

References[LINK]

See the references for the single speed DX cooling coil earlier in this document.

Variable Speed DX Cooling Coil[LINK]

Overview[LINK]

The Coil:Cooling:DX:VariableSpeed object will simulate the performance of a DX cooling coil used in combination with a variable-speed unitary air conditioner and air-to-air heat pump. It fits into the parent objects of AirLoopHVAC:UnitaryHeatCool, ZoneHVAC:PackagedTerminalAirConditioner, AirLoopHVAC:UnitaryHeatPump:AirToAir and ZoneHVAC:PackagedTerminalHeatPump, etc.

The rated conditions for obtaining the capacities, COPs and SHR are indoor dry-bulb temperature at 26.67 ˚C (80 ˚F), wet bulb temperature at 19.44 ˚C (67 ˚F), and the source side entering air temperature at 35 ˚C (95 ˚F).

Variable-speed cooling coils lead to varied dehumidification behaviors, that the Bypass Factor (BF) is not only dependent on the indoor air flow rate, but also on the refrigerant mass flow rate, i.e. the compressor speed. The methods of calculating Bypass factor and Sensible Heat Transfer Ratio at each speed are the same as the water source variable-speed cooling coil (object name Coil:Cooling:WaterToAirHeatPump:VariableSpeedEquationFit).

Model Description[LINK]

The Coil:Cooling:DX:VariableSpeed object is modeled in a manner similar to Coil:Cooling:WaterToAirHeatPump:VariableSpeedEquationFit. Of course, rather than referencing a water-cooled condenser, the new coil object references to outdoor air-cooled condenser, and has the air entering temperature in place of the condenser water entering temperature.

It shall be noted for the total cooling capacity and flow rate inputs, two fields are autosizable, which are Rated Total Cooling Capacity at Selected Nominal Speed Level and Rated Volumetric Air Flow Rate at Selected Nominal Speed Level. They are used to scale the performances of an actual unit and correlate with the actual loop flow. Except the two fields, all other total cooling capacity and flow rate inputs at individual speed levels should be directly obtained from Reference Unit catalog.

The Rated Total Cooling Capacity at Selected Nominal Speed Level contains the rated total cooling capacity to match the building sensible or latent cooling load. The rated cooling capacity is used to determine a capacity scaling factor, as compared to the Reference Unit catalog capacity at the nominal speed level.

$C a p a c i t y S c a l e F a c t o r = \frac{R a t e d T o t a l C o o l i n g C a p a c i t y}{R e f e r e n c e U n i t T o t a l C o o l i n g C a p a c i t y @ N o m i n a l S p e e d L e v e l}$

And then, this scaling factor is used to determine capacities at rated condition for other speed levels, as below,

$R a t e d T o t a l C o o l i n g C a p a c i t y @ S p e e d L e v e l (x) = C a p a c i t y S c a l e F a c t o r \times R e f e r e n c e U n i t T o t a l C o o l i n g C a p a c i t y @ S p e e d L e v e l (x)$

The Rated Volumetric Air Flow Rate is used to determine an internal scaling factor, and calculate the air flow rates in the parent objects, as follows:

$\begin{matrix} D e s i g n A i r M a s s F l o w R a t e @ S p e e d L e v e l (x) = R e f e r e n c e U n i t A i r M a s s F l o w R a t e @ S p e e d L e v e l (x) \times C a p a c i t y S c a l e F a c t o r \end{matrix}$

And the loop volumetric air flow rates in the parent objects are calculated as below,

$L o o p V o l u m e t r i c A i r F l o w R a t e @ S p e e d L e v e l (x) = A i r F l o w S c a l e F a c t o r \times R e f e r e n c e U n i t V o l A i r F l o w R a t e @ S p e e d L e v e l (x) \times C a p a c i t y S c a l e F a c t o r$

If the volumetric air flow rate at one speed level is higher than the flow rate allowed by the fan in the parent object, the flow rate will be set back to the fan flow rate.

If $A i r F l o w S c a l e F a c t o r$ equals to unity, the loop air flow rate become the design flow rates of the original unit (after scaled by the rated total cooling capacity). The Rated Volumetric Air Flow Rate is introduced here to correlate with the actual flow rate in the air loop, in case that these differ from the design specifications. Certainly, it is recommended that the Rated Volumetric Air Flow Rate is selected in the way that $A i r F l o w S c a l e F a c t o r$ is unity, so as to get more accurate results from the performance curves.

If the condenser is evaporatively cooled, it is necessary to know the condenser air flow rate, so as to calculate the water evaporation rate. We will have fields for specifying the Reference Unit condenser volumetric air flow rates at individual speed levels, and these inputs are optional. If the condenser air flow rates are not inputted by the user, default values of condenser volumetric air flow rate as a function of the rated total cooling capacity will be used. Condenser air flow rates of the simulated unit are calculated as below,

$\begin{matrix} C o n d e n s e r V o l u m e t r i c A i r F l o w R a t e @ S p e e d L e v e l (x) = R e f e r e n c e U n i t C o n d e n s e r V o l A i r F l o w R a t e @ S p e e d L e v e l (x) \times C a p a c i t y S c a l e F a c t o r \end{matrix}$

Performance Curves:[LINK]

This object includes 4 curve objects at each individual speed level.

1) Total cooling capacity modifier curve (function of temperature)

2) Total cooling capacity modifier curve (function of air flow fraction)

3) Energy input ratio (EIR) modifier curve (function of temperature)

4) Energy input ratio (EIR) modifier curve (function of air flow fraction)

The flow fraction modifier curves are used as a placeholder, to account for off-design flow rates if needed. If the manufacturer doesn’t provide the off-design performances, we can simply use a default modification multiplier of 1.0.

At the lowest speed, there will be one additional performance curve to correlate the part-load condition, i.e.

Part load fraction correlation (function of part load ratio)

1) Total cooling capacity modifier curve (function of temperature)

The total cooling capacity modifier as a function of temperature curve (CAP-FT) is a biquadratic curve with two independent variables: wet-bulb temperature of the air entering the cooling coil and the condenser entering air temperature. The output of this curve is multiplied by the rated total cooling capacity at the speed, to give the total cooling capacity at the specific entering air WB and condenser inlet air DB at which the DX unit is operating (i.e., at temperatures different from the rating point temperatures).

Note: The data used to develop the total cooling capacity modifier curve (function of temperature) should represent performance when the cooling coil is ‘wet’ (i.e., coil providing sensible cooling and at least some dehumidification). Performance data when the cooling coil is ‘dry’ (i.e., not providing any dehumidification) should not be included when developing this modifier curve. This model automatically detects and adjusts for ‘dry coil’ conditions.

$T o t C a p T e m p M o d F a c = a + b * W B_{i} + {c * W B}_{i}^{2} + d * D B_{o} + e * D {B_{o}}^{2} + f * W B_{i} * D B_{o}$

where

WB $_{i}$ = wet-bulb temperature of the air entering the indoor cooling coil, °C

DB $_{o}$ = dry-bulb temperature of the air entering the condenser coil, °C

a-f = regression curve-fit coefficients.

2) Total cooling capacity modifier curve (function of air flow fraction)

The total cooling capacity modifier curve (function of air flow fraction) is a cubic curve with the independent variable being the ratio of the actual air flow rate across the cooling coil to the design air flow rate (i.e., fraction of design flow at the speed).

$T o t C a p A i r F l o w M o d F a c = a + b * f f_{a} + {c * f f}_{a}^{2} + {d * f f}_{a}^{3}$

where

ff $_{a}$ = actual air mass flow rate/design air mass flow rate, at one speed level.

a-d = regression curve fit coefficients, if no data for correction, the user can simply put a = 1.0, and the other coefficients as 0.0.

Energy input ratio (EIR) modifier curve (function of temperature)

The energy input ratio modifier curve as a function of temperature (EIR-FT) is a biquadratic curve with two independent variables: wet-bulb temperature of the air entering the cooling coil and the condenser entering air temperature. The output of this curve is multiplied by the rated EIR (inverse of the rated COP) at the speed level, to give the EIR at the specific entering air temperatures at which the DX unit is operating (i.e., at temperatures different from the rating point temperatures).

$E I R T e m p M o d F a c = a + b * W B_{i} + {c * W B}_{i}^{2} + d * D B_{o} + e * D {B_{o}}^{2} + f * W B_{i} * D B_{o}$

where

a-f = regression curve fit coefficients

Energy input ratio (EIR) modifier curve (function of air flow fraction)

$E I R A i r F l o w M o d F a c = a + b * f f_{a} + {c * f f}_{a}^{2} + {d * f f}_{a}^{3}$

where

a-d = regression curve fit coefficients, if no data available for correction, the user can simply put a = 1.0, and the other coefficients as 0.0.

Part load fraction correlation (function of part load ratio)

This field defines the name of a quadratic or cubic performance curve (Ref: Performance Curves) that parameterizes the variation of electrical power input to the unit as a function of the part load ratio (PLR, Sensible or Latent Load/Steady-State Sensible or Latent Capacity for Speed 1). The description of the part load fraction correlation for the variable-speed DX cooling coil is similar to the variable-speed DX heating coil.

Lowest Speed Operation:[LINK]

The lowest speed operation of the variable-speed DX cooling coil is similar to the single speed DX cooling coil. The total (gross) cooling capacity of the variable-speed DX coil is calculated as follows:

$\begin{matrix} Q_{t o t a l, 1} = R a t e d T o t a l C o o l i n g C a p a c i t y @ S p e e d L e v e l (1) * T o t C a p T e m p M o d F a c_{1} * T o t C a p A i r F l o w M o d F a c_{1} \end{matrix}$

And the EIR is calculated as:

$E I R_{1} = \frac{1.0}{R e f e r e n c e U n i t C O P @ S p e e d (1)} E I R T e m p M o d F a c_{1} * E I R A i r F l o w M o d F a c_{1}$

And the power consumption excluding the indoor fan is,

$P o w e r_{1} = Q_{t o t a l, 1} * E I R_{1} * R T F$

At the lowest speed, the dehumidification calculation is exactly the same as the single speed DX coil. That is to use the rated SHR and the design air flow rate at the lowest speed to calculate rated bypass factor of BF $_{r a t e d, 1}$ , and the corresponding effective surface area of A $_{o, 1}$ . With A $_{o, 1}$ known, the coil BF can be adjusted for non-rated air flow rates.

And the part load ratio for sensible cooling is,

$P L R = S e n s i b l e C o o l i n g L o a d / (Q_{t o t a l, 1} * S H R_{1})$

For latent cooling,

$P L R = L a t e n t C o o l i n g L o a d / [Q_{t o t a l, 1} * (1.0 - S H R_{1})]$

Higher Speed Operation:[LINK]

At the speed level between the lowest and the highest, there is no part-load loss. A parameter of speed ratio (SpeedRatio) is used to define the capacity partition between Speed x-1 and Speed x.

The design air flow rate at the speed ratio are given as following:

$\begin{matrix} D e s i g n A i r F l o w R a t e S p e e d R a t i o = R e f e r e n c e U n i t A i r M a s s F l o w R a t e @ S p e e d L e v e l (x - 1) * C a p a c i t y S c a l e F a c t o r * (1 - S p e e d R a t i o) + R e f e r e n c e U n i t A i r M a s s F l o w R a t e @ S p e e d L e v e l (x) * C a p a c i t y S c a l e F a c t o r * S p e e d R a t i o \end{matrix}$

And the fraction of air flow is given:

$f f_{a, x - 1}$ = $f f_{a, x}$ = actual air mass flow rate/DesignAirFlowRateSpeedRatio

The total cooling capacities and EIRs at Speed x-1 and Speed x are calculated:

$\begin{matrix} Q_{t o t a l, x - 1} = R a t e d T o t a l C o o l i n g C a p a c i t y @ S p e e d L e v e l (x - 1) * T o t C a p T e m p M o d F a c_{x - 1} * T o t C a p A i r F l o w M o d F a c_{x - 1} \end{matrix}$

$\begin{matrix} Q_{t o t a l, x} = R a t e d T o t a l C o o l i n g C a p a c i t y @ S p e e d L e v e l (x) * T o t C a p T e m p M o d F a c_{x} * T o t C a p A i r F l o w M o d F a c_{x} \end{matrix}$

$Q_{t o t a l, S p e e d R a t i o} = H e a t i n g L o a d$

$E I R_{x} = \frac{1.0}{R e f e r e n c e U n i t C O P @ S p e e d (x)} E I R T e m p M o d F a c_{x} * E I R A i r F l o w M o d F a c_{x}$

The total cooling capacity at the corresponding speed ratio is:

$Q_{t o t a l, S p e e d R a t i o} = (1.0 - S p e e d R a t i o) * Q_{t o t a l, x - 1} + S p e e d R a t i o * Q_{t o t a l, x}$

And the power consumption, excluding the indoor fan, is

$P o w e r_{S p e e d R a t i o} = (1.0 - S p e e d R a t i o) * Q_{t o t a l, x - 1} * E I R_{x - 1} + S p e e d R a t i o * Q_{t o t a l, x} * E I R_{x}$

And the net heat discharged from the condenser:

$Q_{c o n d, S p e e d R a t i o} = P o w e r_{S p e e d R a t i o} + Q_{t o t a l, S p e e d R a t i o}$

And the effective surface area in the correlations of BF factor is calculated as below:

$A_{o, S p e e d R a t i o} = (1 - S p e e d R a t i o) * A_{o, x - 1} + S p e e d R a t i o * A_{o, x}$

Using A $_{o, S p e e d R a t i o}$ in the same BF and SHR calculation procedure as the single speed DX cooling coil, we can get BF $_{S p e e d R a t i o}$ , and SHR $_{S p e e d R a t i o}$ . And the sensible cooling capacity is calculated:

$Q_{s e n s i b l e, S p e e d R a t i o} = Q_{t o t a l, S p e e d R a t i o} * S H R_{S p e e d R a t i o}$

If the variable-speed DX cooling coil used to match the sensible cooling load,

$Q_{s e n s i b l e, S p e e d R a t i o} = S e n s i b l e L o a d$

If it intends to match the latent cooling load,

$Q_{t o t a l, S p e e d R a t i o} * (1.0 - S H R_{S p e e d R a t i o}) = L a t e n t L o a d$

If the speed reaches the highest level, the speed ratio becomes 1.0, and Speed n represents the highest speed.

Evaporatively Cooled Condenser:[LINK]

If using evaporative pre-cooling pad upstream of the condenser, the user must also enter the appropriate evaporative effectiveness for the media at each compressor speed, so as to correlate variation of the effectiveness as a function of compressor speed, i.e. the variation of condenser air flow rate. The evaporative effectiveness is interpolated between speed levels, i.e.

$E f f e c t i v e n e s s_{S p e e d R a t i o} = (1.0 - S p e e d R a t i o) * E f f e c t i v e n e s s_{x - 1} + S p e e d R a t i o * E f f e c t i v e n e s s_{x}$

The interpolated effectiveness at the speed ratio is used to calculate the air temperature reduction across the evaporative cooling pad, and the decreased air temperature entering the condenser coil is used in place of the ambient temperature, i.e.

$T_{c, i} = T_{w b, o} + (1.0 - E f f e c t i v e n e s s_{S p e e d R a t i o}) * (T_{d b, o} - T_{w b, o})$

Where

$T_{c, i}$ is the condenser entering air temperature, °C.

$T_{w b, o}$ is the outdoor air wet bulb temperature, °C.

$T_{d b, o}$ is the outdoor air dry bulb temperature, °C.

And the water consumption rate is calculated based on the interpolated condenser air flow rate at the speed ratio, i.e. the condenser flow rate at the speed ratio,

${˙ m}_{c o n d a i r, S p e e d R a t i o} = (1.0 - S p e e d R a t i o) * {˙ m}_{c o n d a i r, x - 1} + S p e e d R a t i o * {˙ m}_{c o n d a i r, x}$

Electric Air Heating Coil[LINK]

The electric air heating coil (object name: Coil:Heating:Electric) is a simple capacity model with a user-input efficiency. In many cases, this efficiency for the electric coil will be 100%. This coil only has air nodes to connect it in the system. The coil can be used in the air loop simulation or in the zone equipment as a reheat coil. Depending on where it is used determines if this coil is temperature or capacity controlled. If used in the air loop simulation it will be controlled to a specified temperature on the setpoint node by the SetpointManager. If it is used in zone equipment, it will be controlled from the zone thermostat by meeting the zone demand. See Input Output Reference for additional information.

For a coil that is controlled to meet the zone demand and will meet the capacity necessary for that zone unless it exceeds the capacity of the coil specified by the user.

! Control output to meet load QCoilReq

IF((AirMassFlow .GT. 0.0 .AND. HeatingCoil(CoilNum)%NominalCapacity > 0.0) .and. &

(GetCurrentScheduleValue(HeatingCoil(CoilNum)%SchedPtr) .gt. 0.0) .and. &

(QCoilReq .gt. 0.0) .and. (TempSetPoint = = 0.0)) THEN

!check to see if the Required heating capacity is greater than the user! specified capacity.

IF(QCoilReq > HeatingCoil(CoilNum)%NominalCapacity) Then

QCoilCap = HeatingCoil(CoilNum)%NominalCapacity

Else

QCoilCap = QCoilReq

End IF

TempAirOut = TempAirIn + QCoilCap/CapacitanceAir

HeatingCoilLoad = QCoilCap

!The HeatingCoilLoad is the change in the enthalpy of the Heating

HeatingCoil(CoilNum)%ElecUseLoad = HeatingCoilLoad/Effic

For a temperature setpoint coil the delta temperature from the coil inlet temperature to the setpoint is determined and the capacity of the coil is calculated and is met if less than the user specified capacity.

! Control coil output to meet a setpoint temperature.

Else IF((AirMassFlow .GT. 0.0 .AND. HeatingCoil(CoilNum)%NominalCapacity > 0.0) .and. &

(GetCurrentScheduleValue(HeatingCoil(CoilNum)%SchedPtr) .gt. 0.0) .and. &

(QCoilReq = = 0.0) .and. &

(ABS(TempSetPoint-TempAirIn) .gt. TempControlTol) ) THEN

QCoilCap = CapacitanceAir*(TempSetPoint - TempAirIn)

! check to see if setpoint above enetering temperature. If not, set

! output to zero.

IF(QCoilCap .LE. 0.0) THEN

QCoilCap = 0.0

TempAirOut = TempAirIn

!check to see if the Required heating capacity is greater than the user

! specified capacity.

Else IF(QCoilCap > HeatingCoil(CoilNum)%NominalCapacity) Then

QCoilCap = HeatingCoil(CoilNum)%NominalCapacity

TempAirOut = TempAirIn + QCoilCap/CapacitanceAir

Else

TempAirOut = TempSetPoint

End IF

HeatingCoilLoad = QCoilCap

!The HeatingCoilLoad is the change in the enthalpy of the Heating

HeatingCoil(CoilNum)%ElecUseLoad = HeatingCoilLoad/Effic

Gas Air Heating Coil[LINK]

The gas-powered heating coil (object name: Coil:Heating:Gas) is a simple capacity model with user inputted gas burner efficiency. The default for the gas burner efficiency is 80%. This coil only has air nodes to connect it in the system. The coil can be used in the air loop simulation or in the zone equipment as a reheat coil. Depending on where it is used determines if this coil is temperature or capacity controlled. If used in the air loop simulation it will be controlled to a specified temperature scheduled from the SetpointManager. If it is used in zone equipment, it will be controlled from the zone thermostat by meeting the zone demand. The gas coil has additional features that can add a part load correction and a parasitic gas or electric load.

The parasitic electric load associated with the gas coil operation, such as an inducer fan, etc. This will be modified by the PLR (or coil runtime fraction if a part-load fraction correlation is provided in the next input field) to reflect the time of operation in a simulation time step.

The parasitic gas load associated with the gas coil’s operation (Watts), such as a standing pilot light. The model assumes that this parasitic load is consumed only for the portion of the simulation time step where the gas heating coil is not operating.

Field: Part Load Fraction Correlation (function of part load ratio)[LINK]

The part load correction defines the name of a quadratic or cubic performance curve (Ref: Performance Curves) that parameterizes the variation of gas consumption rate by the heating coil as a function of the part load ratio (PLR, sensible heating load/nominal capacity of the heating coil). For any simulation time step, the nominal gas consumption rate (heating load/burner efficiency) is divided by the part-load fraction (PLF) if a part-load curve has been defined. The part-load curve accounts for efficiency losses due to transient coil operation.

The part-load fraction correlation should be normalized to a value of 1.0 when the part load ratio equals 1.0 (i.e., no efficiency losses when the heating coil runs continuously for the simulation time step). For PLR values between 0 and 1 ( 0 < = PLR < 1), the following rules apply:

PLF > = 0.7 and PLF > = PLR

If PLF < 0.7 a warning message is issued, the program resets the PLF value to 0.7, and the simulation proceeds. The runtime fraction of the heating coil is defined a PLR/PLF. If PLF < PLR, then a warning message is issues and the runtime fraction of the coil is limited to 1.0.

A typical part load fraction correlation for a conventional gas heating coil (e.g., residential furnace) would be:

$P L F = 0.8 + 0.2 \cdot P L R$

For a better understanding of how the coil meets the temperature setpoint in the air loop or the zone demand as zone equipment, see Coil:Electric:Heating for additional information. Also see Input Output Reference for additional input information.

Multi-Stage Electric and Gas Air Heating Coil[LINK]

Overview[LINK]

These models (object names Coil:Heating:Electric:MultiStage and Coil:Heating:Gas:MultiStage) simulate the performance of multi stage electric and gas heating systems with two to four discrete air flow stages. Each stage has a set of user inputted gas burner efficiency, Nominal capacity, and Parasitic electric load (for gas coil only). The full load supply airflow rate is dependent on the stage number and is set by its parent object (Ref: AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed). The part-load impact on coil energy use is automatically applied to the lowest stage.

When the model determines performance at Stage 1 (the lowest stage) or cycling between OFF and Stage 1, its performance is almost the same as the performance for the Coil:Heating:Electric and Coil:Heating:Gas models. However, the outlet conditions are calculated slightly differently. Therefore, the Coil:Heating:Electric and Coil:Heating:Gas model may be considered as subsets of the multistage model described here. When the multistage coil model determines performance at higher stages (above 1), the model linearly interpolates the performance at two consecutive stages (n-1 and n) as needed to meet the heating load, with the fraction of time at each speed established by the speed ratio. For the time being, this coil model can only be called by the parent object AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed.

Model Inputs[LINK]

The model inputs are also very similar to the inputs of the Coil:Heating:Electric and Coil:Heating:Gas objects. The main difference is that this multistage model requires a set of fields at each speed, such as nominal capacity, and efficiency.

Stage 1 Operation[LINK]

The following procedure provides the detailed description of the exception.

Total delivered heating capacity

The total delivered heating capacity for Stage 1 operating at the cycling ratio needed to meet the requested heating load is:

$Q_{c o i l, c y c l i n g} = \cdot m S p e e d 1 (C y c R a t i o) (h_{i n l e t} - h_{o u t l e t, f u l l})$

where,

Q $_{c o i l, c y c l i n g}$ = delivered sensible heating capacity for Stage 1 operating at a specific cycling ratio [W]

$\cdot m S p e e d 1$ = air mass flow rate through heating coil at Speed 1 as set by the parent object [kg/s]

h $_{o u t l e t, f u l l}$ = specific enthalpy of the coil outlet air during full-load operation at Stage 1 (no cycling) [J/kg]

h $_{i n l e t}$ = specific enthalpy of the coil inlet air [J/kg]

CycRatio = cycling ratio at Stage 1, ratio of requested heating load to the full-load capacity of the coil at Stage 1 [dimensionless]

It is assumed that the coil provides no heating capacity when the coil is OFF, even if the supply air fan continues to operate.

Outlet air specific enthalpy

The average specific enthalpy of the coil outlet air is then calculated based on the delivered sensible heating capacity and the average air mass flow rate entering the coil:

$h_{o u t l e t, a v e r a g e} = h_{i n l e t} - \frac{Q_{c o i l, c y c l i n g}}{{\cdot m}_{i n l e t}}$

where,

h $_{o u t l e t, a v e r a g e}$ = average specific enthalpy at the coil outlet [J/kg]

$\cdot m i n l e t$ = mass flow rate at the inlet to the coil as established by the parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiStage, Mass Flow Rate Calculation). This flow rate is the average value determined by the parent object, accounting for the specified flow rate when the heating coil is ON and the specified flow rate when the heating coil is OFF for the time step being simulated.

Outlet air temperature

The heating coil’s outlet air humidity ratio equals the inlet air humidity ratio since the coil does not change the moisture content of the air. So the average outlet air temperature is calculated based on the inlet air humidity ratio and the average outlet air enthalpy using the psychrometric function PsyTdbFnHW.

The main reason for using the above approach is that outlet air conditions are calculated in the same way for all operating Stages.

Higher Stage Operation[LINK]

This section describes how higher Stage operation is simulated. When the required sensible load is less than the full load sensible capacity at Stage n (Stage Number > 1), the following calculations are performed:

Total delivered heating capacity at Stage n-1 and Stage n

$T o t C a p_{n - 1} = M S N o m i n a l C a p_{n - 1}$

$T o t C a p_{n} = M S N o m i n a l C a p_{n}$

where,

TotCap $_{i}$ = total delivered heating capacity at given temperatures and flow rates at Stage i [W]

MSNominalCap $_{i}$ = heating capacity at the rated conditions at Stage i [W]

i = Stage n or Stage n-1

Full load outlet air specific enthalpy at Stage n-1 and Stage n

$h_{o u t l e t, f u l l_{S} p e e d_{n}} = h_{i n l e t} - (M S N o m i n a l C a p_{n} / . m i n l e t)$

$h_{o u t l e t, f u l l_{S} p e e d e_{n} - 1} = h_{i n l e t} - (M S N o m i n a l C a p_{n - 1} / . m i n l e t)$

where,

h $_{o u t l e t, f u l l_S t a g e_n}$ = specific enthalpy of the coil outlet air during full-load operation at Stage n (no cycling) [J/kg]

h $_{o u t l e t, f u l l_S t a g e_n - 1}$ = specific enthalpy of the coil outlet air during full-load operation at Stage n-1 (no cycling) [J/kg]

Effective total heating capacity

where,

$Q_{c o i l, S p e e d R a t i o}$ = delivered sensible heating capacity at a given Stage ratio between two consecutive Stages [W]

$\cdot m S p e e d n$ = air mass flow rate through heating coil at Stage n as set by the parent object [kg/s]

$\cdot m S p e e d n - 1$ = air mass flow rate through heating coil at Stage 1 as set by the parent object [kg/s]

Average outlet air enthalpy

$h_{o u t l e t, a v e r a g e} = h_{i n l e t} - \frac{Q_{c o i l, S p e e d R a t i o}}{{\cdot m}_{i n l e t}}$

where,

h $_{o u t l e t, a v e r a g e}$ = average specific enthalpy at the coil outlet [J/kg]

h $_{i n l e t}$ = specific enthalpy of the coil inlet air [J/kg]

${˙ m}_{i n l e t}$ = Mass flow rate at the inlet to the coil as established by the parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiStage, Mass Flow Rate Calculation). This flow rate is the average value determined by the parent object, accounting for the specified flow rate when the heating coil is at Stage n and the specified flow rate when the heating coil is at Stage n-1 for the time step being simulated.

Average outlet air temperature

Full load energy inputs at Stage n-1 and Stage n

$H e a t i n g P o w e r_{n} = M S N o m i n a l C a p_{n} / M S E f f i c i e n c y_{n}$

$H e a t i n g P o w e r_{n - 1} = M S N o m i n a l C a p_{n - 1} / M S E f f i c i e n c y_{n - 1}$

where,

MSEfficiency $_{n}$ = Efficiency at stage n

MSEfficiency $_{n - 1}$ = Efficiency at stage n-1

Calculate combined energy input

$H e a t i n g P o w e r = H e a t i n g P o w e r_{n} (S p e e d R a t i o) + H e a t i n g P o w e r_{n - 1} (1.0 - S p e e d R a t i o)$

where,

HeatingPower = Power used in Watt

Single-Speed Electric Heat Pump DX Air Heating Coil[LINK]

Overview[LINK]

This model (object name Coil:Heating:DX:SingleSpeed) simulates the performance of an air-to-air direct expansion (DX) heating system. The model uses performance information at rated conditions along with curve fits for variations in total capacity, energy input ratio and part-load fraction to determine the performance of the unit at part-load conditions (DOE 1982). Adjustment factors are applied to total capacity and input power to account for frost formation on the outdoor coil.

This model simulates the thermal performance of the indoor DX heating coil, and the power consumption of the outdoor unit (compressors, fans, crankcase heaters and defrost heaters). The performance of the indoor supply air fan varies widely from system to system depending on control strategy (e.g., constant fan vs. AUTO fan, constant air volume vs. variable air volume, etc.), fan type, fan motor efficiency and pressure losses through the air distribution system. Therefore, this DX system model does not account for the thermal effects or electric power consumption of the indoor supply air fan. EnergyPlus contains separate models for simulating the performance of various indoor fan configurations, and these models can be easily linked with the DX system model described here to simulate the entire DX system being considered (e.g., see AirLoopHVAC:UnitaryHeatPump:AirToAir).

Model Inputs[LINK]

The user must input the total heating capacity, coefficient of performance (COP) and the volumetric airflow rate across the heating coil at rated conditions. The capacity and COP inputs should be “gross” values, excluding any thermal or energy impacts due to the indoor supply air fan. The rating condition is considered to be outdoor air at 8.33C dry-bulb and 6.11C wet-bulb temperatures (i.e., air entering the outdoor coil), with air entering the indoor DX heating coil at 21.11C dry-bulb and 15.55C wet-bulb temperatures. The rated volumetric air flow across the DX heating coil should be between 0.00004027 m $^{3}$ /s and 0.00006041 m $^{3}$ /s per watt of rated total heating capacity (300 – 450 cfm/ton).

Depending on the defrost strategy that is selected, the user must also input up to six performance curves that describe the change in total heating capacity and efficiency at part-load conditions, and efficiency during reverse-cycle defrosting:

The total heating capacity modifier curve (function of temperature) can be a function of both the outdoor and indoor air dry-bulb temperature or only the outdoor air dry-bulb temperature. User has the choice of a bi-quadratic curve with two independent variables or a quadratic curve as well as a cubic curve with a single independent variable. The bi-quadratic curve is recommended if sufficient manufacturer data is available as it provides sensitivity to the indoor air dry-bulb temperature and a more realistic output. The output of this curve is multiplied by the rated total heating capacity to give the total heating capacity at specific temperature operating conditions (i.e., at an outdoor or indoor air temperature different from the rating point temperature).

$T o t C a p T e m p M o d F a c = a + b (T_{d b, o}) + c {(T_{d b, o})}^{2}$

or

$T o t C a p T e m p M o d F a c = a + b (T_{d b, o}) + c {(T_{d b, o})}^{2} + d {(T_{d b, o})}^{3}$

$T o t C a p T e m p M o d F a c = a + b (T_{d b, i}) + c {(T_{d b, i})}^{2} + d (T_{d b, o}) + e {(T_{d b, o})}^{2} + f (T_{d b, i}) (T_{d b, o})$

where

$T_{d b, i}$ = dry-bulb temperature of the air entering the indoor coil, °C

$T_{d b, o}$ = dry-bulb temperature of the air entering the outdoor coil, °C
The total heating capacity modifier curve (function of flow fraction) is a quadratic or cubic curve with the independent variable being the ratio of the actual air flow rate across the heating coil to the rated air flow rate (i.e., fraction of full load flow). The output of this curve is multiplied by the rated total heating capacity and the total heating capacity modifier curve (function of temperature) to give the total heating capacity at the specific temperature and air flow conditions at which the coil is operating.

$T o t C a p F l o w M o d F a c = a + b (f f) + c {(f f)}^{2}$

or

$T o t C a p F l o w M o d F a c = a + b (f f) + c {(f f)}^{2} + d {(f f)}^{3}$

where

$ff=flowfraction=(ActualairmassflowrateRatedairmassflowrate)$

Note: The actual volumetric airflow rate through the heating coil for any simulation time step where the DX unit is operating should be between 0.00002684 m $^{3}$ /s and .00008056 m $^{3}$ /s per watt of rated total heating capacity (200 - 600 cfm/ton). The simulation will issue a warning message if this airflow range is exceeded.
The energy input ratio (EIR) modifier curve (function of temperature) can be a function of both the outdoor and indoor air dry-bulb temperature or only the outdoor air dry-bulb temperature. User has the choice of a bi-quadratic curve with two independent variables or a quadratic curve as well as a cubic curve with a single independent variable. The bi-quadratic curve is recommended if sufficient manufacturer data is available as it provides sensitivity to the indoor air dry-bulb temperature and a more realistic output. The output of this curve is multiplied by the rated EIR (inverse of the rated COP) to give the EIR at specific temperature operating conditions (i.e., at an outdoor or indoor air temperature different from the rating point temperature).

$E I R T e m p M o d F a c = a + b (T_{d b, o}) + c {(T_{d b, o})}^{2}$

or

$E I R T e m p M o d F a c = a + b (T_{d b, o}) + c {(T_{d b, o})}^{2} + d {(T_{d b, o})}^{3}$

or

$E I R T e m p M o d F a c = a + b (T_{d b, i}) + c {(T_{d b, i})}^{2} + d (T_{d b, o}) + e {(T_{d b, o})}^{2} + f (T_{d b, i}) (T_{d b, o})$
The energy input ratio (EIR) modifier curve (function of flow fraction) is a quadratic or cubic curve with the independent variable being the ratio of the actual air flow rate across the heating coil to the rated air flow rate (i.e., fraction of full load flow). The output of this curve is multiplied by the rated EIR (inverse of the rated COP) and the EIR modifier curve (function of temperature) to give the EIR at the specific temperature and air flow conditions at which the coil is operating.

$E I R F l o w M o d F a c = a + b (f f) + c {(f f)}^{2}$

or

$E I R F l o w M o d F a c = a + b (f f) + c {(f f)}^{2} + d {(f f)}^{3}$
The part-load fraction correlation (function of part-load ratio) is a quadratic or cubic curve with the independent variable being part-load ratio (sensible heating load / steady-state heating capacity). The output of this curve is used in combination with the rated EIR and EIR modifier curves to give the “effective” EIR for a given simulation time step. The part-load fraction (PLF) correlation accounts for efficiency losses due to compressor cycling.

$P a r t L o a d F r a c = P L F = a + b (P L R) + c {(P L R)}^{2}$

or

$P a r t L o a d F r a c = a + b (P L R) + c {(P L R)}^{2} + d {(P L R)}^{3}$

where

$P L R = p a r t - l o a d r a t i o = (\frac{s e n s i b l e h e a t i n g l o a d}{s t e a d y - s t a t e s e n s i b l e h e a t i n g c a p a c i t y})$

The part-load fraction correlation should be normalized to a value of 1.0 when the part load ratio equals 1.0 (i.e., no efficiency losses when the compressor(s) run continuously for the simulation time step). For PLR values between 0 and 1 ( 0 < = PLR < 1), the following rules apply:

PLF > = 0.7 and PLF > = PLR

If PLF < 0.7 a warning message is issued, the program resets the PLF value to 0.7, and the simulation proceeds. The runtime fraction of the coil is defined a PLR/PLF. If PLF < PLR, then a warning message is issued and the runtime fraction of the coil is limited to 1.0.

A typical part load fraction correlation for a conventional, single-speed DX heating coil (e.g., residential heat pump) would be:

PLF = 0.85 + 0.15(PLR)
The defrost energy input ratio (EIR) modifier curve (function of temperature) is a bi-quadratic curve with two independent variables: outdoor air dry-bulb temperature and the heating coil entering air wet-bulb temperature. The output of this curve is multiplied by the heating coil capacity, the fractional defrost time period and the runtime fraction of the heating coil to give the defrost power at the specific temperatures at which the coil is operating. This curve is only required when a reverse-cycle defrost strategy is specified.

$D e f r o s t E I R T e m p M o d F a c = a + b (T_{w b, i}) + c {(T_{w b, i})}^{2} + d (T_{d b, o}) + e {(T_{d b, o})}^{2} + f (T_{w b, i}) (T_{d b, o})$

where

$T_{w b, i}$ = wet-bulb temperature of the air entering the indoor heating coil, °C

$T_{d b, o}$ = dry-bulb temperature of the air entering the outdoor coil, °C

All six curves are accessed through EnergyPlus’ built-in performance curve equation manager (curve:quadratic, curve:cubic and curve:biquadratic). It is not imperative that the user utilize all coefficients shown in the preceding equations { through } in items (1) through (6) if their performance equation has fewer terms (e.g., if the user’s PartLoadFrac performance curve is linear instead of quadratic or cubic, simply enter the appropriate values for the coefficients a and b, and set the remaining coefficients to zero).

The next input item for the Heating DX single speed coil is the supply air fan operation mode. Either the supply air fan runs continuously while the DX coil cycles on/off, or the fan and coil cycle on/off together. The next two inputs define the minimum outdoor dry-bulb temperature that the heat pump compressor will operate and the maximum outdoor dry-bulb temperature for defrost operation. Crankcase heater capacity and crankcase heater cutout temperature are entered in the following two inputs. The final four inputs cover the type of defrost strategy (reverse-cycle or resistive), defrost control (timed or on-demand), the fractional defrost time period (timed defrost control only), and the resistive defrost heater capacity if a resistive defrost strategy is selected.

Model Description[LINK]

The general flow of the model is as follows:

1) If the outdoor air dry-bulb temperature is below the specified minimum temperature for compressor operation or the DX heating coil is not scheduled to operate, simply pass through the heating coil inlet air conditions as the coil outlet conditions, set power and heating rates equal to zero, and set crankcase heater power equal to the crankcase heater capacity value specified by the input file.

2) If the outdoor air dry-bulb temperature is above the specified minimum temperature for compressor operation and the DX heating coil is scheduled to operate, then:

a. If the outdoor dry-bulb temperature is below the specified maximum outdoor dry-bulb temperature for defrost operation, calculate a heating capacity multiplier, input power multiplier and fractional defrost time period depending on the defrost strategy and defrost control type specified for the heating coil.

b. Using the rated heating capacity and COP, the part-load curves specified for the DX heating coil, the defrost multipliers calculated above (if applicable), and the part-load ratio that is being requested of the heating coil, determine the following: heating coil exiting air conditions (dry-bulb temperature, humidity ratio and enthalpy), total DX coil heating rate, electric power during heating (compressors and outdoor fans), electric power during defrost, and crankcase heater power.

The following paragraphs give a detailed description of the model calculations that are performed when the DX heating coil is operating (i.e., scenario # 2 above).

Frost Adjustment Factors[LINK]

Frost formation on the outdoor coil, and the need to periodically defrost this coil, has a significant impact on heating capacity and energy use by the DX heating system. If the outdoor air dry-bulb temperature is below the specified maximum temperature for defrost operation, then the model calculates adjustment factors for heating capacity and input power due to frost formation, and the fractional defrost time period, depending on the defrost strategy and defrost control type specified for the heating coil. This method of accounting for the impacts of frosting/defrost was taken from the model used in DOE-2.1E (ESTSC 2001, Miller and Jaster 1985).

The model first estimates the outdoor coil temperature according to a linear empirical relationship with outdoor air dry-bulb temperature as the independent variable.

$T_{c o i l, o u t} = 0.82 T_{d b, o} - 8.589$

The difference between the outdoor air humidity ratio (from the weather file) and the saturated air humidity ratio at the estimated outdoor coil temperature is then calculated, and this value is used as an indication of frost formation on the outdoor coil.

$Δ ω_{c o i l, o u t} = M A X [1.0 E - 6, ω_{o u t d o o r} - ω_{s a t} (T_{c o i l, o u t}, O u t B a r o P r e s s)]$

Frost formation on the outdoor coil must be periodically removed. The fraction of compressor runtime when the coil is being defrosted is either entered by the user (for timed defrost) or is calculated by the model (for on-demand defrost) using an empirical equation and $Δ ω_{c o i l, o u t}$ . Adjustment factors to total heating coil capacity and input power due to frost formation on the outdoor coil are also calculated by empirical models with $Δ ω_{c o i l, o u t}$ or fractional defrost time period as the independent variable. The defrost time period fraction and adjustment factors due to frost formation on the outdoor coil vary depending on the defrost control type as shown below.

Timed Defrost:

$FractionalDefrostTime=timeperiodspecifiedbyuser=tfrac,defrost$

$H e a t i n g C a p a c i t y M u l t i p l i e r = 0.909 - 107.33 (Δ ω_{c o i l, o u t})$

$I n p u t P o w e r M u l t i p l i e r = 0.9 - 36.45 (Δ ω_{c o i l, o u t})$

On-Demand Defrost:

$FractionalDefrostTime=11+(0.01446Δωcoil,out)=tfrac,defrost$

$HeatingCapacityMultiplier=0.875(1−tfrac,defrost)$

$InputPowerMultiplier=0.954(1−tfrac,defrost)$

If the outdoor air dry-bulb temperature is above the specified maximum temperature for defrost operation, the fractional defrost time period is set to zero and the heating capacity/input power multipliers are set to unity.

Defrost Operation[LINK]

If the fractional defrost time period is greater than zero for the simulation time step, then the model calculates the electrical power used during defrost. The method for calculating defrost power varies based on the defrost strategy specified (i.e., reverse-cycle or resistive). In the case of reverse-cycle defrost, the additional heating load due to defrost (indoor cooling during defrost) is also calculated so that it may be added to the existing heating load when calculating input power for the compressor(s) and outdoor coil fan(s).

Reverse-Cycle:

$Qdefrost=0.01(tfrac,defrost)(7.222−Tdb,o)(Qtotal,rated1.01667)$

$Pdefrost=DefrostEIRTempModFac(Qtotal,rated1.01667)(tfrac,defrost)(RTF)$

Resistive:

$Q_{d e f r o s t} = 0.0$

$Pdefrost=(Qcap,defrost)(tfrac,defrost)(RTF)$

where:

$Q_{d e f r o s t}$ = additional indoor heating load due to reverse-cycle defrost (W)

$Q_{t o t a l, r a t e d}$ = total full-load heating capacity of the coil at rated conditions (W)

$P_{d e f r o s t}$ = average defrost power for the simulation time step (W)

$Q_{c a p, d e f r o s t}$ = capacity of the resistive defrost heating element (W)

DefrostEIRTempModFac = energy input ratio modifier curve applicable during defrost

$RTF=({{\rm{PLR}}}/PLRPartLoadFrac{{\rm{PartLoadFrac}}})=runtimefractionoftheheatingcoil$

Heating Operation[LINK]

For any simulation time step, the total heating capacity of the DX unit is calculated as follows:

$Q t o t a l = Q t o t a l, r a t e d (T o t C a p T e m p M o d F a c) (T o t C a p F l o w M o d F a c)$

If the outdoor air dry-bulb temperature is below the maximum temperature for defrost operation, then the total heating capacity is further adjusted due to outdoor coil frost formation based on the results of Equation and Equation or .

$Q t o t a l = Q t o t a l (H e a t i n g C a p a c i t y M u l t i p l i e r)$

In a similar fashion, the electrical power draw by the DX unit (compressors plus outdoor coil fans) for any simulation time step is calculated. For a reverse-cycle defrost strategy, the additional heating load (Q $_{d e f r o s t}$ ) generated during defrost operation is added to the heating load being requested by adjusting the part-load ratio. If a resistive defrost strategy is selected, Q $_{d e f r o s t}$ = 0. The part-load fraction correlation for the heating coil (user input, Equation or ) is used in the calculation of electrical power draw to account for efficiency losses due to compressor cycling.

$P L R = M I N (1.0, P L R + (\frac{Q_{d e f r o s t}}{Q_{t o t a l}}))$

$P a r t L o a d F r a c = a + b (P L R) + c {(P L R)}^{2} + d {(P L R)}^{3}$

$P_{h e a t i n g} = \frac{(Q_{t o t a l}) (E I R) (P L R)}{P a r t L o a d F r a c} \times I n p u t P o w e r M u l t i p l i e r$

where

$P_{h e a t i n g}$ = average compressor and outdoor fan power for the simulation time step(W)

$Q_{t o t a l}$ = total heating capacity W, Eqn.

$E I R = E n e r g y i n p u t r a t i o = (\frac{1}{C O P_{r a t e d}}) (E I R T e m p M o d F a c) (E I R F l o w M o d F a c)$

$C O P_{r a t e d}$ = coefficient of performance at rated conditions (user input)

InputPowerMultiplier = power adjustment due to frost if applicable -Eqn. or

The crankcase heater is assumed to operate when the heating coil’s compressor is OFF, and the average crankcase heater power for the simulation time step is calculated as follows:

$P_{c r a n k c a s e} = Q_{c a p, c r a n k c a s e} (1 - R T F)$

$RTF=({{\rm{PLR}}}/PLRPartLoadFrac{{\rm{PartLoadFrac}}})=runtimefractionoftheheatingcoil$

where

$P_{c r a n k c a s e}$ = average crankcase heater power for the simulation time step (W)

$Q_{c a p, c r a n k c a s e}$ = crankcase heater capacity (W)

If this heating coil is used as part of an air-to-air heat pump (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir), the crankcase heater defined for this DX heating coil is enabled during the time that the compressor is not running for either heating or cooling (and the crankcase heater power defined in the DX cooling coil object is disregarded in this case). In this instance, RTF in the above equations would be the runtime fraction of the heat pump’s heating coil or cooling coil, whichever is greater.

The properties of the air leaving the heating coil at full-load operation are calculated using the following equations:

$h_{o u t l e t} = h_{i n l e t} + \frac{Q_{t o t a l}}{∙ m}$

$ω_{o u t l e t} = ω_{i n l e t}$

$T_{d b, o u t l e t} = P s y T d b F n H W (h_{o u t l e t}, ω_{o u t l e t})$

where

$h_{o u t l e t}$ = enthalpy of the air leaving the heating coil (J/kg)

$ω_{o u t l e t}$ = leaving air humidity ratio (kg/kg)

$T_{d b, o u t l e t}$ = leaving air dry-bulb temperature (°C)

PsyTdbFnHW = EnergyPlus psychrometric function, returns dry-bulb temp given enthalpy and humidity ratio

Condenser Inlet Air Temperature[LINK]

The air temperature entering the outdoor condenser is based on the weather data used during a simulation. This temperature is either taken directly from the weather data, or can be adjusted based on the height of the outdoor condenser. Since this DX heating coil can only be used in conjunction with a DX cooling coil, the input for Condenser Air Inlet Node Name in the DX cooling coil object can be used to control this optional feature. If this input is left blank, the air temperature entering the condenser is based solely on the weather data. If this input is not blank, then the node name specified must also be listed in an OutdoorAir:Node or OutdoorAir:NodeList object. When the node name is listed in an OutdoorAir:NodeList object, the air temperature entering the condenser is based solely on the weather data. When the node name is listed in an OutdoorAir:Node object, the height of the node determines the air temperature entering the outdoor condenser (see description of Local Outdoor Air Temperature Calculation in the Atmospheric Variation section of this document for further details).

Supply Air Fan Control: Cycling vs. Continuous[LINK]

One of the inputs to the DX coil model is the supply air fan operation mode: cycling fan, cycling compressor (CyclingFanAndCompressor) or continuous fan, cycling compressor (ContinuousFanWithCyclingCompressor). The first operation mode is frequently referred to as “AUTO fan”, where the compressor(s) and supply air fan operate in unison to meet the zone heating load, and cycle off together when the heating load has been met. The second operation mode is often referred to as “fan ON”, where the compressor(s) cycle on and off to meet the zone heating load but the supply air fan operates continuously regardless of compressor operation.

Since this DX heating coil can only be used in conjunction with a DX cooling coil (i.e. heat pumps), and these coils are used in AC equipment that specifies a fan operation mode schedule (e.g AirLoopHVAC:UnitaryHeatPump:AirToAir), the selection made in the DX heating coil is currently ignored and the fan operation mode schedule value determines the fan operation mode for each time step throughout the simulation. A fan operation mode schedule value of 0 specifies AUTO fan mode operation while values other than 0 specify fan ON operation. The use of a schedule allows the fan operation mode to change based on time-of-day or with changes in season.

With this methodology, the determination of the air mass flow rate and air properties for the two different supply air fan operation modes is slightly different. For the case of cycling fan/cycling compressor, the conditions of the air leaving the heating coil are the steady-state values calculated using equations , and above. However the air mass flow rate passed along to the next component (and eventually to the conditioned zone) is the average air mass flow rate for the system simulation time step (determined by the heating system; see AirLoopHVAC:UnitaryHeatPump:AirToAir). For this fan control type, the heating coil part-load fraction (Equation or ) is also passed to Fan:OnOff (if used) to properly calculate the supply air fan power and associated fan heat.

For the case of continuous fan/cycling compressor, the air mass flow rate is constant. However, the air properties leaving the heating coil are calculated as the average conditions during the system simulation time step. The model assumes that the exiting air conditions are the steady-state values calculated using equations , and above when the compressor(s) operate. For the remainder of the system simulation time step, it is assumed that the air exiting the DX coil has the same properties as the air entering the coil. For this supply air fan operating strategy, the leaving air properties are calculated as follows:

$h_{o u t l e t, C o n t F a n C y c C o m p} = h_{o u t l e t} (P L R) + h_{i n l e t} (1 - P L R)$

$ω_{o u t l e t, C o n t F a n C y c C o m p} = ω_{o u t l e t} (P L R) + ω_{i n l e t} (1 - P L R)$

$T_{d b, o u t l e t, C o n t F a n C y c C o m p} = P s y T d b F n H W (h_{o u t l e t, C o n t F a n C y c C o m p}, ω_{o u t l e t, C o n t F a n C y c C o m p})$

References[LINK]

DOE. 1982. DOE-2 engineers manual, version 2.1A. LBL-11353. Berkeley, CA: Lawrence Berkeley National Laboratory.

ESTSC. 2001. DOE-2.1E Version 110 (source code).Oak Ridge, TN: Energy Science and Technology Software Center.

Miller, R.L. and Jaster, H. 1985. Performance of Air-Source Heat Pumps. EM-4226. Palo Alto, CA: Electric Power Research Institute.

Single-Speed DX Heating Coil Standard Ratings[LINK]

For single-speed direct expansion (DX) heating coils, the industry standard ratings of High Temperature Heating Standard (Net) Rating Capacity, Low Temperature Heating Standard (Net) Rating Capacity and Heating Seasonal Performance Factor (HSPF) are calculated according to ANSI/AHRI Standard 210/240 (AHRI 2008). The rated Energy Efficiency Ratio (EER) is not calculated for any DX heatings coils at this time.

For the Coil:Heating:DX:SingleSpeed object in EnergyPlus, these standard ratings are not direct inputs to the model. However, these standard ratings can be calculated using user-entered information for the Coil:Heating:DX:SingleSpeed object. These standard rating values are provided in the eplusout.eio output file (Ref. OutputDetailsAndExamples.pdf) and also in the predefined tabular output reports (Output:Table:SummaryReports object, Equipment Summary). Currently, the standard ratings are only calculated and output for single-speed DX heating coils..

Note: The standard ratings described in this section require that the DX heating coil model be evaluated at specific operating conditions (i.e., specific temperatures for air entering the heating coil and for air entering the air-cooled [outdoor] condenser). If the heating coil performance curves can not be evaluated at the required test conditions, then a standard rating value will be calculated at user specified curve limit as an output and a warning message will written to eplusout.err. For example, if the curve object (Curve:Biquadratic) for Total Heating Capacity Function of Temperature Curve has a minimum value of -5C for dry-bulb temperature entering the condenser coil, the HSPF calculation requires that heating capacity and EIR be calculated at -8.33C, so this would result in HSPF value calculated at -5C as an output and a warning message in the eplusout.err file.

High Temperature Heating Standard (Net) Rating Capacity[LINK]

The procedure for calculating the Standard Rating Heating Capacity is given by:

${˙ Q}_{H i g h S t a n d a r d R a t i n g} = {˙ Q}_{t o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{R a t e d})$

${˙ Q}_{H i g h S t a n d a r d R a t i n g N e t} = {˙ Q}_{H i g h S t a n d a r d R a t i n g} + F a n_{P o w e r, R a t e d}$

where,

${˙ Q}_{H i g h S t a n d a r d R a t i n g N e t}$ = Total standard (net) heating capacity (W) of the air-source heat pump equipment in heating mode determined from ANSI/AHRI Standard 210/240 and test conditions shown in Table. The standard heating test conditions for air-cooled condenser are: indoor coil entering air dry-bulb 21.1°C and outdoor coil entering air dry-bulb 8.33°C.

${˙ Q}_{H i g h S t a n d a r d R a t i n g}$ = Total standard heating capacity (W) of the air-source heat pump equipment in heating mode determined from ANSI/AHRI Standard 210/240 and test conditions shown in Table. The standard heating test conditions for air-cooled condenser are: indoor coil entering air dry-bulb 21.1°C and outdoor coil entering air dry-bulb 8.33°C.

${˙ Q}_{t o t a l, R a t e d}$ = Rated total heating capacity, user input (W)

$T o t C a p T e m p M o d F a c_{R a t e d}$ = User-specified bi-quadratic curve evaluated at the indoor coil entering air dry-bulb 21.1°C and outdoor coil entering air dry-bulb 8.33°C for air-cooled condenser as specified in Table, (dimensionless).

$F a n_{P o w e r, R a t e d}$ = Supply air fan power at rated conditions, (W). The Rated Indoor Coil Fan Power Per Volume Flow rate is a user-entered value, with a default of 773.3 W/(m3/s)) if the user leaves this input field blank. The default value is taken from ANSI/ASHRAE Standard 210/240-2008 where it is defined for systems which do not have an Indoor Coil fan furnished as part of the system. See the description given at the end of this section for how this value is calculated.

Low Temperature Heating Standard (Net) Rating Capacity[LINK]

The procedure for calculating the Low Temperature Standard Rating Heating Capacity is given by:

${˙ Q}_{L o w S t a n d a r d R a t i n g} = {˙ Q}_{t o t a l, R a t e d} (T o t C a p T e m p M o d F a c_{H 3 T e s t})$

${˙ Q}_{L o w S t a n d a r d R a t i n g N e t} = {˙ Q}_{L o w S t a n d a r d R a t i n g} + F a n_{P o w e r, R a t e d}$

where,

${˙ Q}_{L o w S t a n d a r d R a t i n g N e t}$ = Total standard (net) heating capacity (W) of the air-source heat pump equipment in heating mode determined from ANSI/AHRI Standard 210/240 and test conditions shown in Table. The standard heating test conditions for air-cooled condenser are: indoor coil entering air dry-bulb 21.1°C and outdoor coil entering air dry-bulb (-8.33) °C.

${˙ Q}_{L o w S t a n d a r d R a t i n g}$ = Total standard heating capacity (W) of the air-source heat pump equipment in heating mode determined from ANSI/AHRI Standard 210/240 and test conditions shown in Table. The standard heating test conditions for air-cooled condenser are: indoor coil entering air dry-bulb 21.1°C and outdoor coil entering air dry-bulb (-8.33) °C.

${˙ Q}_{t o t a l, R a t e d}$ = Rated total heating capacity, user input (W)

$T o t C a p T e m p M o d F a c_{H 3 T e s t}$ = User-specified bi-quadratic curve evaluated at the indoor coil entering air dry-bulb 21.1°C and outdoor coil entering air dry-bulb (-8.33) °C for air-cooled condenser as specified in Table, (dimensionless).

Heating Seasonal Performance Factor (HSPF)[LINK]

Heating Seasonal Performance Factor (HSPF) is defined as the total space heating required during the space heating season, divided by the total electrical energy consumed by the heat pump system during the same season.

Calculations of HSPF of a heat pump having a single-speed compressor that was tested with a fixed speed indoor fan installed, a constant-air-volume-rate indoor fan installed, or with no indoor fan installed is given below.

$H S P F = \frac{J \sum j [\frac{n_{j}}{N} \cdot B L (T_{j})]}{J \sum j [\frac{e_{h} (T_{j})}{N}] + J \sum j [\frac{R H (T_{j})}{N}]} \cdot F_{d e f}$

Where,

$\frac{e_{h} (T_{j})}{N}$ = The ratio of the electrical energy consumed by the heat pump during periods of the space heating season when the outdoor temperature fell within the range represented by bin temperature $T_{j}$ to the total number of hours in the heating season (N), W.

$\frac{R H (T_{j})}{N}$ = The ratio of the electrical energy used for resistive space heating during periods when the outdoor temperature fell within the range represented by bin temperature $T_{j}$ to the total number of hours in the heating season (N),W.

$T_{j}$ = The outdoor bin temperature, °C. Outdoor temperatures are “binned” such that calculations are only performed based one temperature within the bin. Bins of 2.78 °C are used.

$\frac{n_{j}}{N}$ = Fractional bin hours for the heating season; the ratio of the number of hours during the heating season when the outdoor temperature fell within the range represented by bin temperature $T_{j}$ to the total number of hours in the heating season, dimensionless. Obtain $\frac{n_{j}}{N}$ values from Table # 17 of AHRI Std. 210/240-2008 (Table below).

$j$ = The bin number, dimensionless.

$J$ = For each generalized climatic region, the total number of temperature bins, dimensionless. Referring to Table # 17, $J$ is the highest bin number ( $j$ ) having a nonzero entry for the fractional bin hours for the generalized climatic region of interest.

$F_{d e f}$ = The demand defrost credit described in section 3.9.2 of AHRI Std. 210/240-2008, dimensionless. For simplification, assigned a value of 1 for timed defrost control and a value of 1.03 for demand defrost control.

$B L (T_{j})$ = The building space conditioning load corresponding to an outdoor temperature of $T_{j}$ ; the heating season building load also depends on the generalized climatic region’s outdoor design temperature and the design heating requirement, W.

Building heating load can be evaluated by using the following equation

$B L (T_{j}) = \frac{(18.33 - T_{j})}{(18.33 - T_{O D})} \cdot C \cdot D H R$

Where,

$T_{O D}$ = The outdoor design temperature, °C. An outdoor design temperature is specified for each generalized climatic region in Table # 17 of AHRI Std. 210/240-2008 (Table below).

$C$ = 0.77, a correction factor which tends to improve the agreement between calculated and measured building loads, dimensionless.

$D H R$ = Design Heating Requirement, W. This is the amount of heating required to maintain a given indoor temperature at a particular outdoor design temperature.

For a single speed heat pump with a fixed speed indoor fan installed, a constant-air-volume-rate indoor fan installed, or with no indoor fan installed, the minimum and maximum design heating requirements for each generalized climatic region can be calculated as follows:

$D H R_{min} = ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ \begin{matrix} \cdot Q h (8.33) \cdot (1.8) \cdot [\frac{18.33 - T_{O D}}{60.0}], f o r r e g i o n s I, I I, I I I, I V, a n d V I \cdot Q h (8.33), f o r r e g i o n s V \end{matrix} ⎫ ⎪ ⎪ ⎪ ⎬ ⎪ ⎪ ⎪ ⎭$

$X (T_{j}) = ⎧ ⎪ ⎪ ⎨ ⎪ ⎪ ⎩ \begin{matrix} B L (T_{J}) / \cdot Q_{h} (T_{j}) O r 1 \end{matrix} ⎫ ⎪ ⎪ ⎬ ⎪ ⎪ ⎭$

Both $D H R_{min}$ and $D H R_{max}$ above should be rounded to the nearest standardized DHR given in Table.

The intermediate term used in the calculations of HSPF can be calculated using the equations described below.

$\frac{e_{h} (T_{j})}{N} = \frac{X (T_{j}) \cdot \cdot E_{h} (T_{j}) \cdot δ (T_{j})}{P L F_{j}} \cdot \frac{n_{j}}{N}$

$\frac{R H (T_{j})}{N} = [B L (T_{j}) - (X (T_{j}) \cdot \cdot Q_{h} (T_{j}) \cdot δ (T_{j}))] \cdot \frac{n_{j}}{N}$

Where,

$j$

whichever is less; the heating mode load factor for temperature bin $j$ , dimensionless.

$\cdot Q_{h} (T_{j})$ = the space heating capacity of the heat pump when operating at outdoor temperature $T_{j}$ , W.

$\cdot E_{h} (T_{j})$ = the electrical power consumption of the heat pump when operating at outdoor temperature $T_{j}$ , W.

$δ (T_{j})$ = the heat pump low temperature cut-out factor, dimensionless.

$P L F_{j}$ = the part load factor, dimensionless.

Part load factor is determined as follows:

$P L F_{j} = 1 - C_{D^{h}} [1 - X (T_{j})]$

Where,

$C_{D^{h}}$ = Degradation coefficient which is defaulted to the value of 0.25

Low temperature cut-out factor $δ (T_{j})$ is determined as follows:

$δ (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ \begin{matrix} 0, i f T_{j} \leq T_{o f f} o r \frac{\cdot Q_{h} (T_{j})}{\cdot E_{h} (T_{j})} < 1 1 / 2, i f T_{o f f} < T_{j} \leq T_{o n} a n d \frac{\cdot Q_{h} (T_{j})}{\cdot E_{h} (T_{j})} \geq 1 1, i f T_{j} > T_{o n} a n d \frac{\cdot Q_{h} (T_{j})}{\cdot E_{h} (T_{j})} \geq 1 \end{matrix}$

where,

$T_{o f f}$ = the outdoor temperature when the compressor is automatically shut off, °C. (If no such temperature exists, $T_{j}$ is always greater than $T_{o f f}$ and $T_{o n}$ ).

$T_{o n}$ = the outdoor temperature when the compressor is automatically turned back on, if applicable, following an automatic shut-off, °C.

Calculate $\cdot Q_{h} (T_{j})$ and $\cdot E_{h} (T_{j})$ using,

$\cdot Q_{h} (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ \begin{matrix} \cdot Q_{h} (- 8.33) + \frac{[\cdot Q_{h} (8.33) - \cdot Q_{h} (- 8.33)] \cdot (T_{j} + 8.33)}{16.67}, i f T_{j} \geq {7.22}^{\circ} C o r T_{j} \leq - {8.33}^{\circ} C \cdot Q_{h} (- 8.33) + \frac{[\cdot Q_{h} (1.67) - \cdot Q_{h} (- 8.33)] \cdot (T_{j} + 8.33)}{10}, i f - {8.33}^{\circ} C < T_{j} < {7.22}^{\circ} C \end{matrix}$

$\cdot E_{h} (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ \begin{matrix} \cdot E_{h} (- 8.33) + \frac{[\cdot E_{h} (8.33) - \cdot E_{h} (- 8.33)] \cdot (T_{j} + 8.33)}{16.67}, i f T_{j} \geq {7.22}^{\circ} C o r T_{j} \leq - {8.33}^{\circ} C \cdot E_{h} (- 8.33) + \frac{[\cdot E_{h} (1.67) - \cdot E_{h} (- 8.33)] \cdot (T_{j} + 8.33)}{10}, i f - {8.33}^{\circ} C < T_{j} < {7.22}^{\circ} C \end{matrix}$

$\cdot Q_{h} (- 8.33), \cdot Q_{h} (1.67), \cdot Q_{h} (8.33), \cdot E_{h} (- 8.33), \cdot E_{h} (1.67)$ and $\cdot E_{h} (8.33)$ are the (net) values calculated using performance curves entered by the user and supply fan power per rated air flow rate.

Heating Mode Test Conditions for Units Having a Single-Speed Compressor and a Fixed-Speed Indoor Fan, a Constant Air Volume Rate Indoor Fan, or No Indoor Fan
Test description			Heating Air Volume Rate
(r)2-3 (l)4-5	Dry Bulb	Wet Bulb	Dry Bulb	Wet Bulb
H1 Test (Required, Steady)	21.11	15.56	8.33	6.11	Heating Full-load
H1C Test (Required, Steady)	21.11	15.56	8.33	6.11	(2)
H2 Test (Required, Steady)	21.11	15.56	1.67	0.56	Heating Full-load
H3 Test (Required, Steady)	21.11	15.56	-8.33	-9.44	Heating Full-load

Generalized Climatic Region Information
	I	II	III	IV	V	VI
	750	1250	1750	2250	2750	*2750
	2.78	-2.78	-8.33	-15	-23.33	-1.11
$j$	$T_{j}$ (°C)
1	16.67	.291	.215	.153	.132	.106	.113
2	13.89	0.239	0.189	0.142	0.111	0.092	0.206
3	11.11	0.194	0.163	0.138	0.103	0.086	0.215
4	8.33	0.129	0.143	0.137	0.093	0.076	0.204
5	5.56	0.081	0.112	0.135	0.1	0.078	0.141
6	2.78	0.041	0.088	0.118	0.109	0.087	0.076
7	0	0.019	0.056	0.092	0.126	0.102	0.034
8	-2.78	0.005	0.024	0.042	0.087	0.094	0.008
9	-5.56	0.001	0.008	0.021	0.055	0.074	0.003
10	-8.33	0	0.002	0.009	0.036	0.055	0
11	-11.11	0	0	0.005	0.026	0.047	0
12	-13.89	0	0	0.002	0.013	0.038	0
13	-16.67	0	0	0.001	0.006	0.029	0
14	-19.44	0	0	0	0.002	0.018	0
15	-22.22	0	0	0	0.001	0.01	0
16	-25	0	0	0	0	0.005	0
17	-27.78	0	0	0	0	0.002	0
18	-30.56	0	0	0	0	0.001	0

Standardized Design Heating Requirements (W)
Source: Table 18, Page 100, ANSI/AHRI Standard 210/240 -2008
1465.356	7326.78	14653.56	26376.41
2930.712	8792.136	17584.27	29307.12
4396.068	10257.49	20514.98	32237.83
5861.424	11722.85	23445.7	38099.26

Multi-Speed Electric Heat Pump DX Air Heating Coil[LINK]

Overview[LINK]

This model (object name Coil:Heating:DX:MultiSpeed:) simulates the performance of an air-to-air direct expansion (DX) heating system. The main difference compared to the other heating coil model (Coil:Heating:DX:SingleSpeed) is that this heating coil allows modeling of two to four discrete compressor speeds. Each speed has a set of corresponding performance information at rated conditions along with curve fits for variations in total capacity, energy input ratio and part-load fraction to determine the performance of the unit at part-load conditions (DOE 1982). The full load supply airflow rate is dependent on the speed number and is set by its parent object (Ref: AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed). The part-load impact on coil energy use is automatically applied to the lowest speed. A choice is provided to determine whether part-load impacts on coil energy use are applied when the coil is operating at speeds greater than Speed 1. Adjustment factors applied to total capacity and input power to account for frost formation on the outdoor coil are calculated at each speed.

This model simulates the thermal performance of the indoor DX heating coil, and the power consumption of the outdoor unit (multispeed compressor, fans, crankcase heaters and defrost heaters). The performance of the indoor supply air fan varies widely from system to system depending on control strategy (e.g., constant fan vs. AUTO fan), fan type, fan motor efficiency and pressure losses through the air distribution system. Therefore, this DX system model does not account for the thermal effects or electric power consumption of the indoor supply air fan. EnergyPlus contains separate models for simulating the performance of various indoor fan configurations, and these models can be easily linked with the DX system model described here to simulate the entire DX system being considered. For the time being, this coil model can only be called by the parent object AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed.

When the model determines performance at Speed 1 (the lowest speed) or cycling between OFF and Speed 1, its performance is almost the same as the performance for the Coil:Heating:DX:SingleSpeed model. However, the outlet conditions are calculated slightly differently. Therefore, the Coil:Heating:DX:SingleSpeed model may be considered as a subset of the model described here. When the multispeed coil model determines performance at higher speeds (above 1), the model linearly interpolates the performance at two consecutive speeds (n-1 and n) as needed to meet the heating load, with the fraction of time at each speed established by the speed ratio.

When single mode operation is specified at higher speeds (above 1), defined in the UnitarySystemPerformance:Multispeed object, a child object of AirLoopHVAC:UnitarySystem, its performance is almost the same as the performance for the Coil:Heating:DX:SingleSpeed model at different flow rate and capacity with given speed number. No linear interpolation is performed between two adjacent speeds. The highest speed which does not exceed the sensible load is used for each HVAC timestep. This operation does not apply to speed 1 operation and the case when the sensible load is above the highest sensible capacity.

Model Inputs[LINK]

The model inputs are also very similar to the inputs of the Coil:Heating:DX:SingleSpeed object. The main difference is that this multispeed model requires a set of fields at each speed, such as rated capacity, rated COP, two capacity modifiers, two energy input ratio modifiers, and part-load correction. The inputs also include waste heat fraction and modifier as a function of temperature to calculate recoverable waste heat for heat recovery, which are not available in the similar Coil:Heating:DX:SingleSpeed object.

Speed 1 Operation[LINK]

The calculation procedures in this model, including defrost and crankcase heater, are indentical to the Coil:Heating:DX:SingleSpeed object (Ref: Coil:Heating:DX:SingleSpeed) with one exception: outlet node condition calculation when the supply air fan operation mode is ContinuousFanWithCyclingCompressor. The following procedure provides the detailed description of the exception.

Total delivered heating capacity

The total delivered heating capacity for speed 1 operating at the cycling ratio needed to meet the requested heating load is:

$houtlet,average=hinlet−{{Q_{coil,cycling}}}/Qcoil,cycling⋅m{\mathop m\limits^\cdot }inlet$

where,

Q $_{c o i l, c y c l i n g}$ = delivered sensible heating capacity for Speed 1 operating at a specific cycling ratio [W]

$\cdot m S p e e d 1$ = air mass flow rate through heating coil at Speed 1 as set by the parent object [kg/s]

h $_{o u t l e t, f u l l}$ = specific enthalpy of the coil outlet air during full-load operation at Speed 1 (no cycling) [J/kg]

h $_{i n l e t}$ = specific enthalpy of the coil inlet air [J/kg]

CycRatio = cycling ratio at Speed 1, ratio of requested heating load to the full-load capacity of the coil at Speed 1 [dimensionless]

It is assumed that the coil provides no heating capacity when the coil is OFF, even if the supply air fan continues to operate.

Outlet air specific enthalpy

The average specific enthalpy of the coil outlet air is then calculated based on the delivered sensible heating capacity and the average air mass flow rate entering the coil:

$h_{o u t l e t, a v e r a g e} = h_{i n l e t} - \frac{Q_{c o i l, c y c l i n g}}{{˙ m}_{i n l e t}}$

where,

h $_{o u t l e t, a v e r a g e}$ = average specific enthalpy at the coil outlet [J/kg]

${˙ m}_{i n l e t}$ = mass flow rate at the inlet to the coil as established by the parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed, Mass Flow Rate Calculation). This flow rate is the average value determined by the parent object, accounting for the specified flow rate when the heating coil is ON and the specified flow rate when the heating coil is OFF for the time step being simulated.

Outlet air temperature

The main reason for using the above approach is that outlet air conditions are calculated in the same way for all operating speeds.

The crankcase heater defined for this DX heating coil is enabled during the time that the compressor is not running for either heating or cooling. The crankcase heater power use from either heating or cooling is reported in the heating coil.

Higher Speed Operation[LINK]

Total delivered heating capacity at Speed n-1 and Speed n

$T o t C a p_{n - 1} = R a t e d C a p_{n - 1} (T o t C a p T e m p M o d F a c_{n - 1}) (T o t C a p F l o w M o d F a c_{n - 1})$

$T o t C a p_{n} = R a t e d C a p_{n} (T o t C a p T e m p M o d F a c_{n}) (T o t C a p F l o w M o d F a c_{n})$

where,

TotCap $_{i}$ = total delivered heating capacity at given temperatures and flow rates at Speed i [W]

RatedCap $_{i}$ = heating capacity at the rated conditions at Speed i [W]

TotCapTempModFac $_{i}$ = total heating capacity modifier as a function of indoor and outdoor air dry-bulb temperature at Speed i

TotCapFlowModFac $_{i}$ = total heating capacity modifier as a function of the ratio of the actual flow rate across the heating coil to the rated airflow rate at Speed i

i = Speed n or Speed n-1

EIR at Speed n-1 and Speed n

$E I R_{n - 1} = R a t e d E I R_{n - 1} (E I R T e m p M o d F a c_{n - 1}) (E I R F l o w M o d F a c_{n - 1})$

$E I R_{n} = R a t e E I R_{n} (E I R T e m p M o d F a c_{n}) (E I R F l o w M o d F a c_{n})$

where,

EIR $_{i}$ = energy input ratio at given temperatures and flow rates at Speed i [W]

RatedEIR $_{i}$ = energy input ratio at the rated conditions at Speed i [W]

EIRTempModFac $_{i}$ = energy input ratio modifier as a function of indoor and outdoor air dry-bulb temperature at Speed i

EIRFlowModFac $_{i}$ = energy input ratio modifier as a function of the ratio of the actual flow rate across the heating coil to the rated airflow rate at Speed i

i = Speed n or Speed n-1

Full load outlet air specific enthalpy at Speed n-1 and Speed n

$h_{o u t l e t, f u l l_{S} p e e d, n} = h_{i n l e t} - \frac{T o t C a p_{n} \cdot H e a t i n g C a p a c i t y M u l t i p l i e r}{{˙ m}_{i n l e t}}$

$h_{o u t l e t, f u l l_{S} p e e d, n - 1} = h_{i n l e t} - \frac{T o t C a p_{n - 1} \cdot H e a t i n g C a p a c i t y M u l t i p l i e r}{{˙ m}_{i n l e t}}$

where,

HeatingCapacityMultiplier = frost adjustment factor for heating capacity (See Ref. Coil:Heating:DX:SingleSpeed)

h $_{o u t l e t, f u l l_S p e e d n}$ = specific enthalpy of the coil outlet air during full-load operation at Speed n (no cycling) [J/kg]

h $_{o u t l e t, f u l l_S p e e d n - 1}$ = specific enthalpy of the coil outlet air during full-load operation at Speed n-1 (no cycling) [J/kg]

Effective total heating capacity

$\begin{matrix} Q_{c o i l, S p e e d R a t i o} & = (S p e e d R a t i o) {˙ m}_{S p e e d n} (h_{i n l e t} - h_{o u t l e t, f u l l_{S} p e e d n}) + (1 - S p e e d R a t i o) {˙ m}_{S p e e d n - 1} (h_{i n l e t} - h_{o u t l e t, f u l l_{S} p e e d n - 1}) \end{matrix}$

where,

$Q_{c o i l, S p e e d R a t i o}$ = delivered sensible heating capacity at a given speed ratio between two consecutive speeds [W]

$\cdot m S p e e d n$ = air mass flow rate through heating coil at Speed n as set by the parent object [kg/s]

$\cdot m S p e e d n - 1$ = air mass flow rate through heating coil at Speed 1 as set by the parent object [kg/s]

Average outlet air enthalpy

$h_{o u t l e t, a v e r a g e} = h_{i n l e t} - \frac{Q_{c o i l, S p e e d R a t i o}}{{˙ m}_{i n l e t}}$

where,

h $_{o u t l e t, a v e r a g e}$ = average specific enthalpy at the coil outlet [J/kg]

h $_{i n l e t}$ = specific enthalpy of the coil inlet air [J/kg]

${˙ m}_{i n l e t}$ = Mass flow rate at the inlet to the coil as established by the parent object (Ref. AirLoopHVAC:UnitaryHeatPump:AirToAir:MultiSpeed, Mass Flow Rate Calculation). This flow rate is the average value determined by the parent object, accounting for the specified flow rate when the heating coil is at Speed n and the specified flow rate when the heating coil is at Speed n-1 for the time step being simulated.

Average outlet air temperature

Full load energy inputs at Speed n-1 and Speed n

$H e a t i n g P o w e r_{n} = T o t C a p_{n} (E I R_{n}) (H e a t i n g C a p a c i t y M u l t i p l i e r) (I n p u t P o w e r M u l t i p l i e r)$

$H e a t i n g P o w e r_{n - 1} = T o t C a p_{n - 1} (E I R_{n - 1}) (H e a t i n g C a p a c i t y M u l t i p l i e r) (I n p u t P o w e r M u l t i p l i e r)$

where,

InputPowerMultiplier = Frost adjustment factor for heating power calculation (Ref. Coil:Heating:DX:SingleSpeed)

Calculate combined energy input

When the input for the field ‘Apply Part Load Fraction to Speeds Greater than 1’ is No (equivalent to a single compressor), the combined energy output is calculated as follows:

$H e a t i n g P o w e r = H e a t i n g P o w e r_{n} (S p e e d R a t i o) + H e a t i n g P o w e r_{n - 1} (1.0 - S p e e d R a t i o)$

When the input for the field ’Apply Part Load Fraction to Speeds Greater than 1’ is Yes (equivalent to multiple compressors), the combined energy output is calculated as follows:

$H e a t i n g P o w e r = H e a t i n g P o w e r_{n} (R T F) + H e a t i n g P o w e r_{n - 1} (1.0 - R T F)$

where,

HeatingPower = Power used in Watt

RTF = Run time fraction (SpeedRatio/Part-load Fraction) at Speed n

Calculate defrost power

When the defrost strategy is resistive, the power calculation is the same as Speed 1 operation (Ref. Coil:Heating:DX:SingleSpeed). When the defrost strategy is reverse-cycle, the following calculations are performed:

$Qdefrost,n=0.01(tfrac,defrost)(7.222−Tdb,o)(Qtotal,rated,n1.01667)$

$Pdefrost,n−1=DefrostEIRTempModFac(Qtotal,rated,n−11.01667)(tfrac,defrost)$

$Pdefrost,n=DefrostEIRTempModFac(Qtotal,rated,n1.01667)(tfrac,defrost)$

where:

$Q_{d e f r o s t, n}$ = additional indoor heating load due to reverse-cycle defrost at Speed n (W)

$Q_{t o t a l, r a t e d, n}$ = total full-load heating capacity of the coil at rated conditions at Speed n (W)

$P_{d e f r o s t, n - 1}$ = full load defrost power for the simulation time step at Speed n-1 (W)

$P_{d e f r o s t, n}$ = full load defrost power for the simulation time step at Speed n (W)

$Q_{t o t a l, r a t e d, n - 1}$ = capacity of the resistive defrost heating element at Speed n-1 (W)

$Q_{t o t a l, r a t e d, n}$ = capacity of the resistive defrost heating element at Speed n (W)

DefrostEIRTempModFac = defrost energy input ratio (EIR) modifier curve (Ref. Coil:Heating:DX:SingleSpeed).

T $frac,defrost$ = fractional defrost time (Ref. Coil:Heating:DX:SingleSpeed)

When the input for the field ‘Apply Part Load Fraction to Speeds Greater than 1’ is No (equivalent to a single compressor), the average defrost power is calculated as follows:

$P_{d e f r o s t} = P_{d e f r o s t, n} (S p e e d R a t i o) + P_{d e f r o s t, n - 1} (1.0 - S p e e d R a t i o)$

When the input for the field ’Apply Part Load Fraction to Speeds Greater than 1’ is Yes (equivalent to multiple compressors), the combined defrost energy is calculated as follows:

$P_{d e f r o s t} = P_{d e f r o s t, n} (R T F) + P_{d e f r o s t, n - 1} (1.0 - R T F)$

where,

$P_{d e f r o s t}$ = average defrost power used in Watt

RTF = Run time fraction (SpeedRatio/Part-load Fraction) at Speed n

Crankcase heater

There is no power need at higher speed operation.

Higher Speed Operation with Single Mode Operation[LINK]

This section describes how higher speed operation is simulated with single mode operation. When the required sensible load is less than the full load sensible capacity at Speed n (Speed Number > 1), its performance is almost the same as the performance for the Coil:Heating:DX:SingleSpeed model at different flow rate and capacity with given Speed n. No linear interpolation is performed between two adjacent speeds. The highest speed which does not exceed the sensible load is used for each HVAC timestep. This operation does not apply to the case when the sensible load is above the highest sensible capacity.

Waste heat calculation[LINK]

The waste heat generated by this coil object is calculated as:

$Q_{W a s t e H e a t} = (F r a c t i o n) (T e m p M o d i f i e r) (H e a t i n g P o w e)$

where

Fraction = rated waste heat fraction of the energy input

TempModifier = waste heat modifier as a function of indoor and outdoor air dry-bulb temperature

Standard Rating of Multi-Speed DX Heating Coils[LINK]

For multi-speed direct expansion heating coils, the industry standard ratings of Standard Rating Heating Seasonal Performance Factor (HSPF) are calculated according to NSI/AHRI Standard 210/240 (AHRI 2008). These standard ratings can be calculated using the user-entered information for the Coil:Heating:DX:MultiSpeed object. According to Standard 210/240, the Standard Rating HSPF applies to air-to-air unitary heat pumps with rated heating capacities under 19,000 watts (65,000 Btu/h). The Heating Mode Test Conditions for Units Having two-capacity Compressor standard tests H0 $_{1}$ , H1 $_{1}$ , H2 $_{1}$ , H3 $_{1}$ , H1 $_{2}$ , H2 $_{2}$ , and H3 $_{2}$ are also provided in Table. The rated Energy Efficiency Ratio (EER) is not calculated for any DX Heating Coils as this time. The equations required to calculate the net heating capacity and HSPF values are described next.

Standard Ratings Heating Capacity[LINK]

The Standard Rating (Net) Heating Capacity calculation is calculated using the same procedure as the single speed DX heating coils except the capacities are evaluated at the maximum speed only. Refer to the standard ratings calculation for Coil:DX:Heating:SingelSpeed for description of the procedure.

The Rated Supply Fan Power Per Volume Flow Rate is a user-entered value, with a default of 773.3 W/(m $^{3}$ /s)) if the user leaves this input field blank. The default value is taken from ANSI/ASHRAE Standard 210/240-2008 where it is defined for systems which do not have an Indoor Coil (Evaporator) fan furnished as part of the system. The test conditions in ANSI/AHRI Standard 210/240 vary the external static pressure (i.e., pressure drop associated with ductwork and other devices external to the indoor fan/coil section) seen by the supply air fan based on the standard rating cooling capacity. New input field Rated Supply Fan Power Per Volume Flow Rate for each speed will be added to the IDD.

Heating seasonal performance factors (HSPF) for Multi-Speed DX Coil[LINK]

The following section describes the test condition for heating mode tests for Air-Source Heat Pumps, including heating-only heat pumps, and the equations required to calculate the HSPF. The HSPF for multi-speed compressor heat pumps in Btu/W-h is given by:

$H S P F = \frac{M \sum j \frac{n_{j}}{N} \cdot B L (T_{j})}{M \sum j \frac{e_{h} (T_{j})}{N} + M \sum j \frac{R H (T_{j})}{N}} \cdot F_{d e f}$

Where,

BL(T $_{j}$ ) = the building space conditioning load corresponding to an outdoor temperature of T $_{j}$ ; the heating season building load also depends on the generalized climatic region’s outdoor design temperature and the design heating requirement, Btu/h.

e $_{h}$ (T $_{j}$ )/N = the ratio of the electrical energy consumed by the heat pump during periods of the space heating season when the outdoor temperature fell within the range represented by bin temperature T $_{j}$ to the total number of hours in the heating season (N), W. For heat pumps having a heat comfort controller, this ratio may also include electrical energy used by resistive elements to maintain a minimum air delivery temperature.

RH(T $_{j}$ )/N = the ratio of the electrical energy used for resistive space heating during periods when the outdoor temperature fell within the range represented by bin temperature T $_{j}$ to the total number of hours in the heating season (N),W. Resistive space heating is modeled as being used to meet that portion of the building load that the heat pump does not meet because of insufficient capacity or because the heat pump automatically turns off at the lowest outdoor temperatures.

T $_{j}$ = the outdoor bin temperature, °C. Outdoor temperatures are “binned” such that calculations are only performed based one temperature within the bin. Bins of 2.78°C are used.

n $_{j}$ /N = Fractional bin hours for the heating season; the ratio of the number of hours during the heating season when the outdoor temperature fell within the range represented by bin temperature T $_{j}$ to the total number of hours in the heating season, dimensionless, given in Table.

j = the bin number, dimensionless.

M = for each generalized climatic region, the total number of temperature bins, dimensionless. Referring to Table, M is the highest bin number (j) having a nonzero entry for the fractional bin hours for the generalized climatic region of interest.

The building heating load is calculated as follows:

$B L (T_{j}) = \frac{18 .33 - T_{j}}{18 .33 - T_{O D}} \cdot C \cdot D H R$

The minimum and maximum design heating requirements for each generalized climate regions are given by:

$D H R_{m i n} = ⎧ ⎨ ⎩ \begin{matrix} 8.33 {˙ Q}_{h}^{k = 2} \frac{18.33 - T_{O D}}{33.33} & for regions I, II, III, IV, and VI 8.33 {˙ Q}_{h}^{k = 2} & for region V \end{matrix}$

and

$D H R_{m i n} = ⎧ ⎨ ⎩ \begin{matrix} 2 \cdot 8.33 {˙ Q}_{h}^{k = 2} \frac{18.33 - T_{O D}}{33.33} & for regions I, II, III, IV, and VI 2 \cdot 8.33 {˙ Q}_{h}^{k = 2} & for region V \end{matrix}$

Where,

C = 0.77, a correction factor which tends to improve the agreement between calculated and measured building loads dimensionless.

DHR = the design heating requirement for each generalized climate region, Btu/h.

${˙ Q}_{h}^{k = 2} (8.33)$ = the heating capacity determined from H1 $_{2}$ test, (W)

The heating capacity delivered and electric power input of the heat pump when operating at minimum (low) compressor capacity and outdoor air temperature of T $_{j}$ is given by:

${˙ Q}_{h}^{k = 1} (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ \begin{matrix} {˙ Q}_{h}^{k = 1} (8.33) + \frac{[{˙ Q}_{h}^{k = 1} (16 {.67) - ˙ Q}_{h}^{k = 1} (8 .33)] \cdot (T_{j} - 8.33)}{(16.67 - 8.33)} i f T_{j} \geq 4.44^{\circ} C {˙ Q}_{h}^{k = 1} (- 8.33) + \frac{[{˙ Q}_{h}^{k = 1} (1 {.67) - ˙ Q}_{h}^{k = 1} (- 8 .33)] \cdot (T_{j} + 8.33)}{(1.67 + 8.33)} i f - 8.33^{\circ} C \leq T_{j} < 4.44^{\circ} C {˙ Q}_{h}^{k = 1} (- 8.33) + \frac{[{˙ Q}_{h}^{k = 1} (8 {.33) - ˙ Q}_{h}^{k = 1} (- 8 .33)] \cdot (T_{j} + 8.33)}{(8.33 + 8.33)} i f T_{j} < - 8.33^{\circ} C \end{matrix} ⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭$

${˙ E}_{h}^{k = 1} (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎩ \begin{matrix} {˙ E}_{h}^{k = 1} (8.33) + \frac{[{˙ E}_{h}^{k = 1} (16 {.67) - ˙ E}_{h}^{k = 1} (8 .33)] \cdot (T_{j} - 8.33)}{(16.67 - 8.33)} i f T_{j} \geq 4.44^{\circ} C {˙ E}_{h}^{k = 1} (- 8.33) + \frac{[{˙ E}_{h}^{k = 1} (1 {.67) - ˙ E}_{h}^{k = 1} (- 8 .33)] \cdot (T_{j} + 8.33)}{(1.67 + 8.33)} i f - 8.33^{\circ} C \leq T_{j} < 4.44^{\circ} C {˙ E}_{h}^{k = 1} (- 8.33) + \frac{[{˙ E}_{h}^{k = 1} (8 {.33) - ˙ E}_{h}^{k = 1} (- 8 .33)] \cdot (T_{j} + 8.33)}{(8.33 + 8.33)} i f T_{j} < - 8.33^{\circ} C \end{matrix} ⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎬ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎭$

Where,

${˙ Q}_{h}^{k = 1} (16 .67), a n d {˙ E}_{h}^{k = 1} (16 .67)$ are determined from H0 $_{1}$ test and

${˙ Q}_{h}^{k = 1} (8 .33), a n d {˙ E}_{h}^{k = 1} (8 .33)$ are determined from H1 $_{1}$ test and

${˙ Q}_{h}^{k = 1} (1 .67), a n d {˙ E}_{h}^{k = 1} (1 .67)$ are determined from H2 $_{1}$ test, and

${˙ Q}_{h}^{k = 1} (- 8 .33), a n d {˙ E}_{h}^{k = 1} (- 8 .33)$ are determined from H3 $_{1}$ test.

An alternative to conducting or evaluating the performance at H2 $_{1}$ test is to approximate the capacity and electric power inputs as follows:

${˙ Q}_{h}^{k = 1} (1 .67) = 0.9 \cdot {{˙ Q}_{h}^{k = 1} (- 8 .33) + 0.6 [{˙ Q}_{h}^{k = 1} (8 .33) - {˙ Q}_{h}^{k = 1} (- 8 .33)]}$

${˙ E}_{h}^{k = 1} (1 .67) = 0.985 \cdot {{˙ E}_{h}^{k = 1} (- 8 .33) + 0.6 [{˙ E}_{h}^{k = 1} (8 .33) - {˙ E}_{h}^{k = 1} (- 8 .33)]}$

The heating capacity delivered and electric power input of the heat pump when operating at maximum (high) compressor speed and outdoor air temperature of T $_{j}$ is given by:

${˙ Q}_{h}^{k = 2} (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ \begin{matrix} {˙ Q}_{h}^{k = 2} (- 8.33) + \frac{[{˙ Q}_{h}^{k = 2} (8.33) - {˙ Q}_{h}^{k = 2} (- 8.33)] \cdot (T_{j} + 8.33)}{8.33 + 8.33}, I F T_{j} \geq {7.2}^{\circ} C o r T_{j} \leq - {8.33}^{\circ} C {˙ Q}_{h}^{k = 2} (- 8.33) + \frac{[{˙ Q}_{h}^{k = 2} (1.67) - {˙ Q}_{h}^{k = 2} (- 8.33)] \cdot (T_{j} + 8.33)}{1.67 + 8.33}, I F - {8.33}^{\circ} C < T_{j} < {7.2}^{\circ} C \end{matrix}$

${˙ E}_{h}^{k = 2} (T_{j}) = ⎧ ⎪ ⎪ ⎪ ⎨ ⎪ ⎪ ⎪ ⎩ \begin{matrix} {˙ E}_{h}^{k = 2} (- 8.33) + \frac{[{˙ E}_{h}^{k = 2} (8.33) - {˙ E}_{h}^{k = 2} (- 8.33)] \cdot (T_{j} + 8.33)}{8.33 + 8.33}, I F T_{j} \geq {7.2}^{\circ} C o r T_{j} \leq - {8.33}^{\circ} C {˙ E}_{h}^{k = 2} (- 8.33) + \frac{[{˙ E}_{h}^{k = 2} (1.67) - {˙ E}_{h}^{k = 2} (- 8.33)] \cdot (T_{j} + 8.33)}{1.67 + 8.33}, I F - {8.33}^{\circ} C < T_{j} < {7.2}^{\circ} C \end{matrix}$

Where,

${˙ Q}_{h}^{k = 2} (8 .33), a n d {˙ E}_{h}^{k = 2} (8 .33)$ are determined from H1 $_{2}$ test

${˙ Q}_{h}^{k = 2} (1 .67), a n d {˙ E}_{h}^{k = 2} (1 .67)$ are determined from H2 $_{2}$ test

${˙ Q}_{h}^{k = 2} (- 8 .33), a n d {˙ E}_{h}^{k = 2} (- 8 .33)$ are determined from H3 $_{2}$ test

The electric energy used by the heat pumps differs depending up on whether the heat pump would operate at low (minimum) capacity, cycle between successive lower and higher speed capacities, or operate at high speed capacity in responding to the building load. The procedure for each operation cases is described next.

Case 1: The steady state heating capacity when the unit is operating at minimum or low capacity, i.e., when the building heating load is less or equal to the low capacity:

$δ^{'} (T_{j}) = ⎧ ⎪ ⎨ ⎪ ⎩ \begin{matrix} 0, i f T_{j} \leq T_{o f f} 1 / 2, i f T_{o f f} < T_{j} \leq T_{o n} 1, i f T_{j} > T_{o n} \end{matrix} ⎫ ⎪ ⎬ ⎪ ⎭$

$\frac{e_{h} (T_{j})}{N} = \frac{X^{k = 1} (T_{j}) \cdot {˙ E}_{h}^{k = 1} (T_{j}) \cdot δ (T_{j})}{P L F_{j}} \cdot \frac{n_{j}}{N}$

${˙ Q}_{h}^{k = 1} (T_{j}) < B L (T_{j}) < {˙ Q}_{h}^{k = 2} (T_{j})$

$X^{k = 1} (T_{j}) = M I N {\begin{matrix} B L (T_{j}) {/ ˙ Q}_{h}^{k = 1} (T_{j}) 1.0 \end{matrix}}$

$δ^{'} (T_{j}) = ⎧ ⎪ ⎨ ⎪ ⎩ \begin{matrix} 0, i f T_{j} \leq T_{o f f} 1 / 2, i f T_{o f f} < T_{j} \leq T_{o n} 1, i f T_{j} > T_{o n} \end{matrix} ⎫ ⎪ ⎬ ⎪ ⎭$

$P L F_{j} = 1 - C_{D}^{h} \cdot [1 - X^{k = 1} (T_{j})]$

$C_{D}^{h} = M I N {\begin{matrix} C_{D}^{h}, c a l c u l a t e d \end{matrix}$

Engineering Reference — EnergyPlus 8.5

Engineering Reference

Coils[LINK]

Chilled-Water-Based Air Cooling Coil[LINK]

Heat Transfer and Energy Balance[LINK]

Model Description[LINK]

Intermediate calculated U-Factor Times Area Values: The Crux of the Model[LINK]

Design Block Calculations:[LINK]

Variable UA[LINK]

Operating Block Calculations:[LINK]

Effectiveness Equations:[LINK]

Coil Outlet Conditions:[LINK]

Coil Completely Dry Calculations: (operating block)[LINK]

Coil Completely Wet Calculations: (operating block)[LINK]

Coil Part Wet Part Dry Calculations: (operating block)[LINK]

References[LINK]

Chilled-Water-Based Detailed Geometry Air Cooling Coil[LINK]

Heat Transfer and Energy Balance[LINK]

Underlying Correlations, Properties, and Assumptions[LINK]

Solution Method of Model[LINK]

Application of Cooling Coil Model to Heating Coils[LINK]

Hot-Water-Based Air Heating Coil[LINK]

Overview[LINK]

Model Description[LINK]

References[LINK]

Single-Speed Electric DX Air Cooling Coil[LINK]

Overview[LINK]

Model Description[LINK]

Dry Coil Conditions[LINK]

SHR Calculation Using User Specified SHR Modifier Curves[LINK]

Condenser Options: AirCooled vs. EvaporativelyCooled[LINK]

Condenser Inlet Air Temperature[LINK]

Supply Air Fan Control: Cycling vs. Continuous[LINK]

Latent Capacity Degradation[LINK]

Basin Heater For Single-Speed DX Coil[LINK]

Special Calculations for Coil:Cooling:DX:TwoStageWithHumidityControlMode with CoilPerformance:DX:Cooling[LINK]

Bypass Operation[LINK]

Multi-Stage Operation[LINK]

Enhanced Dehumidification Mode[LINK]

Standard Rating of Single-Speed DX Cooling Coils[LINK]

ANSI/ASHRAE 127 - Standard Ratings of Single-Speed DX Cooling Coils[LINK]

Basin Heater For Two-Stage DX Coil[LINK]

References[LINK]

Multi-Speed Electric DX Air Cooling Coil[LINK]

Overview[LINK]

Model Inputs[LINK]

Speed 1 Operation[LINK]

Higher Speed Operation[LINK]

Higher Speed Operation with Single Mode Operation[LINK]

Waste heat calculation[LINK]

Basin Heater For Multi-Speed DX Coil[LINK]

Standard Rating of Multi-Speed DX Cooling Coils[LINK]

Standard Rating Cooling Capacity[LINK]

Seasonal Energy Efficiency Ratio (SEER) for Multi-Speed DX Coil[LINK]

References[LINK]

Two-Speed Electric DX Air Cooling Coil[LINK]

Overview[LINK]

Inputs and Data[LINK]

Calculation[LINK]

Simulation and Control[LINK]

Basin Heater For Two-Speed DX Coil[LINK]

Standard Ratings[LINK]

References[LINK]

Variable Speed DX Cooling Coil[LINK]

Overview[LINK]

Model Description[LINK]

Performance Curves:[LINK]

Lowest Speed Operation:[LINK]

Higher Speed Operation:[LINK]

Evaporatively Cooled Condenser:[LINK]

Electric Air Heating Coil[LINK]

Gas Air Heating Coil[LINK]

Field: Part Load Fraction Correlation (function of part load ratio)[LINK]

Multi-Stage Electric and Gas Air Heating Coil[LINK]

Overview[LINK]

Model Inputs[LINK]

Stage 1 Operation[LINK]

Higher Stage Operation[LINK]

Single-Speed Electric Heat Pump DX Air Heating Coil[LINK]

Overview[LINK]