The object simulates the performance of a chiller-heater which can receive pre-cooled or pre-heated water from the source loop, and provide cooling, heating, or simultaneous cooling-heating. The object needs to work with the Central Heat Pump System object to be controlled properly. This model does not simulate the thermal performance or the power consumption of associated pumps or cooling towers. The Central Heat Pump System object holds the input/output nodes connection of the chiller-heater and its control scheme, once the chiller-heater is properly referred.

## Model Description[LINK]

The model uses user-input performance information at design conditions along with three performance curves (curve objects) for cooling capacity and efficiency to determine chiller operation at off-design conditions. Three additional performance curves for heating capacity and efficiency are used when the chiller is operating in a heating-only mode or simultaneous cooling-heating mode.

### Cooling-only mode[LINK]

The following nomenclature is used in the cooling equations:

*CompMotorEffic* is the compressor motor efficiency

*CompPowerclg* is the compressor power (W)

*CompPower@PLRmin* is the compressor power at the minimum part-load ratio (W)

Cp,cw is the chilled water specific heat (J/kg-K)

*CyclingRatio* is the compressor cycling ratio which is equal to *PLRactual* / *PLRmin*

*EvapCapAvailclg* is the available full-load cooling capacity at current conditions (W)

*EvapCapFTclg* is the cooling capacity function of temperature curve

*EIRFTclg* is the electric input to cooling output factor for temperature function curve

*EIRFPLRclg* is the electric input to cooling output factor for part-load function curve

˙mcw is the chilled water mass flow rate (kg/s)

˙mavail,cw is the chilled water maximum available mass flow rate (kg/s)

*PLRclg* is the cooling part-load ratio which is equal to *CoolingLoad* / *EvapCapAvailclg*

*PLRactual* is the actual part-load ratio at current conditions

*PLRmin* is the minimum part-load ratio

QCondclg is the total condenser heat transfer energy (J)

˙QCondclg is the condenser heat transfer rate (W)

QEvapclg is the total evaporator heat transfer energy (J)

˙QEvapclg is the evaporator heat transfer rate (W)

˙Qfalseloading,clg is the false loading rate (W)

*RefCOPclg* is the reference coefficient of performance (W/W)

*RefEvapCapclg* is the reference evaporator capacity (W)

*FullLoadPwrclg* is the reference full load power which is equal to *EvapCapAvailclg* / *RefCOPclg* (W)

*Tcond* is either the entering or the leaving condenser water temperature depending on user input for the condenser water independent variable: *Tcond,l* if “LeavingCondenser” is chosen or *Tcond,e* if “EnteringCondenser” is chosen.

*Tcond,e* is the entering condenser water temperature (∘C)

*Tcond,l* is the leaving condenser water temperature (∘C)

*Tcw,e* is the entering chilled water temperature (W)

*Tcw,l* is the leaving chilled water temperature (W)

ΔTcw is the chilled water inlet and outlet temperature difference (∘C)

ΔTcw,max is the maximum chilled water inlet and outlet temperature difference (∘C).

The model sequentially calls each chiller-heater module in the order defined in the Central Heat Pump System object. It then determines cooling load that each chiller-heater needs to meet and water flow rates delivered to each chiller-heater. Once each chiller-heater is assumed to operate, it determines cooling capacity and efficiency using user-supplied performance information.

Three performance curves are used in the calculation of cooling capacity and efficiency as follows:

Cooling mode cooling capacity function of temperature curve (*EvapCapFTclg*)

Cooling mode electric input to cooling output ratio function of temperature curve (*EIRFTclg*)

Cooling mode electric input to cooling output ratio function of part load ratio curve (*EIRFPLRclg*)

The Cooling Capacity Function of Temperature Curve (*EvapCapFTclg*) represents the fraction of the cooling capacity of the chiller-heater as it varies by temperature. The curve should have a value of 1.0 at the reference conditions. The output of a bi-quadratic curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EvapCapFTclg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The Cooling Mode Electric Input to Cooling Output Ratio Function of Temperature (EIRFTclg) curve represents the fraction of electricity to the chiller-heater at full load as it varies by temperature. The output of a bi-quadratic curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EIRFTclg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The Cooling Mode Electric Input to Cooling Output Ratio Function of Part Load Ratio (EIRFPLRclg) curve represents the fraction of electricity to the chiller-heater as the load on the chiller varies at a given set of operating temperatures. The curve is normalized so that at full load the value of the curve should be 1.0. Note that the bi-cubic formulation below is generally only valid when LeavingCondenser variable is chosen for the field of Cooling Mode Condenser Water Temperature Curve Input Variable whereas the quadratic curve can be used for both choices, i.e., LeavingCondenser and EnteringCondenser. Bi-cubic may also be used when the chiller-heater uses a variable-speed compressor motor drive. The output of this curve can be determined by one of the following three performance curves:

EIRFPLRclg=a+bTCond,l+cTCond,l2+dPLR+ePLR2+fTCond,lPLR+g⋅0+hPLR3+i⋅0+j⋅0

EIRFPLRclg=a+bPLR+cPLR2

EIRFPLRclg=a+bPLR+cPLR2+dPLR3

The full-load cooling capacity at specific temperature operating conditions (i.e., at temperatures different from the design temperatures) is then computed as follows:

EvapCapAvailclg=RefEvapCapclg⋅EvapCapFTclg

The model then determines current chiller-heater’s evaporator heat transfer rate based on the total cooling load required a central heat pump system to meet and the maximum available chiller-heater cooling capacity. The maximum evaporator temperature difference between the entering chilled water temperature (*Tcw,e*) and the leaving chilled water temperature (*Tcw,l*) obtained from the plant loop setpoint temperature can also be determined. It then calculates mass flow rate for variable flow control chiller-heaters and the temperature difference for constant flow control chiller-heaters, setting the cooling load each chiller-heater needs to meet equal to the evaporator heat transfer rate.

As for variable flow control chiller-heaters, the chilled water mass flow rate is computed as follows:

˙mcw=˙QEvapclgCp,cwΔTcw,max

The chilled water mass flow rate calculated is then compared to the maximum available mass flow rate for individual chiller-heaters. If the calculated one is bigger than the maximum, the model sets the chilled water mass flow rate equal to the maximum. It then adjusts the temperature difference based on the evaporator heat transfer rate and the maximum mass flow rate. If the adjusted temperature difference also exceeds the maximum, the model finally adjusts the evaporator heat transfer rate at the maximum temperature difference and mass flow rate as follows:

˙QEvapclg=˙mavail,cwCp,cwΔTcw,max

As for constant flow control chiller-heaters, the model calculates chilled water temperature difference as follows:

ΔTcw=˙QEvapclg˙mavail,cwCp,cw

The temperature difference calculated is then compared to the maximum temperature difference allowed. If the calculated one is bigger than the maximum, the model sets the chilled water temperature difference equal the maximum, and then adjusts the evaporator heat transfer rate at the given conditions as follows:

˙QEvapclg=˙mavail,cwCp,cwΔTcw,max

The model then calculates the part-load ratio as the ratio of the evaporator heat transfer rate to the available chiller-heater capacity as follows:

PLRclg=˙QEvapclg˙QEvapAvailclg

The part-load ratio calculated is set to be between the maximum of 1.0 and the minimum of 0.0 when it is out of the range. Once the part-load ratio is calculated the cycling ratio and false loading rate can be obtained as follows:

CyclingRatio=PLRactualPLRmin

˙Qfalseloading,clg=EvapCapAvailclg⋅PLRclg⋅CyclingRatio−˙QEvapclg

The compressor power demand is then computed by:

CompPowerclg=RefFullLoadPowerclg⋅EIRFTclg⋅EIRFPLRclg⋅CyclingRatio

The heat transfer rate for the chiller-heater condenser can then be computed as follows:

˙QCondclg=ComPowerclg⋅CompMotorEffic+˙QEvapclg+˙Qfalseloading,clg

The total heat transfer energy by the evaporator and condenser can be calculated as follows:

QEvapclg=˙QEvapclg⋅TimeStepSys⋅3600

QCondclg=˙QCondclg⋅TimeStepSys⋅3600

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:

Vwater=˙mair(ωevapcond,out−ωevapcond,in)ρwater(PLR)(TimeStepSys)(3600)

where:

Vwater is the DX cooling coil evap condenser water consumption (m3)

˙mair is the evaporative condenser air mass flow rate (kg/s)

ωevapcond,in is the humidity ratio of outdoor air entering the evap condenser (kg/kg)

ωevapcond,out is the humidity ratio of air leaving the evap condenser (kg/kg)

ρwater is the density of water at the outdoor dry-bulb temperature (kg/m3)

*TimeStepSys* is the HVAC system simulation time step (hr).

### Heating-only mode and Simultaneous cooling-heating mode[LINK]

The following nomenclature is used in the heating equations:

*CompMotorEffic* is the compressor motor efficiency

*CompPowerhtg* is the compressor power demand (W)

*CompPower@PLRmin* is the compressor power at the minimum part-load ratio (W)

Cp,evap is the evaporator water specific heat (J/kg-K)

Cp,hw is the hot water specific heat (J/kg-K)

*CyclingRatio* is the compressor cycling ratio which is equal to *PLRactual* / *PLRmin*

*EvapCapAvailhtg* is the available full-load cooling capacity at current conditions (W)

*EvapCapFThtg* is the heating mode cooling capacity function of temperature curve

*EIRFThtg* is the electric input to cooling output factor for temperature function curve

*EIRFPLRhtg* is the electric input to cooling output factor for part-load function curve

˙mavail,evap is the evaporator water maximum available mass flow rate (kg/s)

˙mavail,hw is the condenser water maximum available mass flow rate (kg/s)

˙mhw is the hot water mass flow rate (kg/s)

*PLRhtg* is the cooling part-load ratio which is equal to *RefCap* / *EvapCapAvailhtg*

*PLRmax* is the maximum part-load ratio at current conditions

*PLRmin* is the minimum part-load ratio

QCondhtg is the total condenser heat transfer energy (J)

˙QCondAvailhtg is the available full-load heating capacity at current conditions (W)

˙QCondhtg is the condenser heat transfer rate (W)

QEvaphtg is the total evaporator heat transfer energy (J)

˙QEvaphtg is the evaporator heat transfer rate (W)

˙Qfalseloading,htg is the false loading rate (W)

*RefCOPhtg* is the reference coefficient of performance (W/W)

*RefEvapCaphtg* is the reference evaporator capacity (W)

*FullLoadPwrhtg* is the reference full load power which is equal to *EvapCapAvailhtg*/ *RefCOPhtg* (W)

*Tcond* is either entering or leaving condenser water temperature depending on user input for condenser water independent variable. *Tcond,l* if “LeavingCondenser” is chosen or *Tcond,e* if “EnteringCondenser” is chosen.

*Tcond,e* is the entering condenser water temperature (∘C)

*Tcond,l* is the leaving condenser water temperature (∘C)

*Tcw,l* is the leaving chilled water temperature (∘C)

*Thw,e* is the entering hot water temperature (∘C)

*Thw,l* is the leaving hot water temperature (∘C)

ΔTevap is the evaporator inlet and outlet water temperature difference (∘C)

ΔThw is the hot water inlet and outlet temperature difference (∘C)

ΔThw,max is the maximum hot water inlet and outlet temperature difference (∘C).

The calculations for the evaporator side are similar to the cooling-only mode calculations. The evaporator capacity and efficiency is determined by a different set of three performance curves read in the cooling-only mode, and the performance curve set is used for both heating-only mode and simultaneous cooling-heating mode. During these modes, the evaporator side is not connected to the chilled water loop, but source water loop. The model thus assumes that each chiller-heater does not meet the plant loop chilled water setpoint temperature while the evaporator operates at the full load capacity to produce heating at a constant water flow rate.

The model sequentially calls each chiller-heater module in the order of the definition in the central heat pump system. It then determines heating load that each chiller-heater needs to meet and water flow rates delivered to each chiller-heater. Once each chiller-heater is assumed to operate, it determines heating capacity and efficiency using the following performance curves:

Heating mode cooling capacity function of temperature curve (*EvapCapFThtg*)

Heating mode electric input to cooling output ratio function of temperature curve (*EIRFThtg*)

Heating mode electric input to cooling output ratio function of part load ratio curve (*EIRFPLRhtg*)

The output of a Heating Mode Cooling Capacity Function of Temperature curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EvapCapFThtg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The output of a Heating Mode Cooling Output Ratio Function of Temperature curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EIRFThtg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The output of Heating Mode Cooling Output Ratio Function of Part Load Ratio curve can be determined by one of the following three performance curves as follows:

EIRFPLRhtg=a+bTCond,l+cTCond,l2+dPLR+ePLR2+fTCond,lPLR+g⋅0+hPLR3+i⋅0+j⋅0

EIRFPLRhtg=a+bPLR+cPLR2

EIRFPLRhtg=a+bPLR+cPLR2+dPLR3

The full-load evaporator capacity at specific temperature operating conditions is then given by:

EvapCapAvailhtg=RefEvapCaphtg⋅EvapCapFThtg

The part-load ratio is set to be between zero and the maximum, and the evaporator heat transfer rate is computed by:

˙QEvaphtg=EvapCapAvailhtgPLRhtg

The evaporator inlet and outlet temperature difference is then given by:

ΔTevap=˙QEvaphtg˙mavail,evapCp,evap

Once the part-load ratio is calculated the cycling ratio and false loading rate are computed by:

CyclingRatio=PLRactualPLRmin

˙Qfalseloading,htg=EvapCapAvailhtg⋅PLRhtg⋅CyclingRatio−˙QEvaphtg

The compressor power demand is then computed by:

CompPowerhtg=FullLoadPowerhtg⋅EIRFThtg⋅EIRFPLRhtg⋅CyclingRatio

The heat transfer rate of the chiller-heater condenser is then computed as follows:

˙QCondAvailhtg=ComPowerhtg⋅CompMotorEffic+˙QEvaphtg+˙Qfalseloading,htg

Once condenser available heating capacity is determined, the model calculates current chiller-heater’s condenser heat transfer rate based on the total heating load required a central heat pump system to meet as well as available heating capacity of the chiller-heater. The maximum condenser temperature difference between the entering hot water temperature (*Thw,e*) and the leaving hot water temperature (*Thw,l*) obtained from the plant loop setpoint temperature can also be obtained. It then calculates condenser water mass flow rate for variable flow control chiller-heaters and the hot water temperature difference for constant flow control chiller-heaters, setting the cooling load that each chiller-heater needs to meet equal the evaporator heat transfer rate.

As for variable flow control chiller-heaters, the condenser water mass flow rate is computed as follows:

˙mhw=˙QCondhtgCp,hwΔThw,max

The condenser water mass flow rate calculated is then compared to the maximum available mass flow rate for individual chiller-heaters. If the calculated one is bigger than the maximum, the model sets the condenser water mass flow rate equal the maximum. It then adjusts the hot water temperature difference at the maximum mass flow rate. If the adjusted temperature difference also exceeds the maximum, the model finally adjusts the condenser heat transfer rate at the maximum allowable conditions as follows:

˙QCondhtg=˙mavail,hwCp,hwΔThw,max

As for constant flow control chiller-heaters, the model calculates condenser temperature difference as follows:

ΔThw=˙QCondhtg˙mavail,hwCp,hw

The temperature difference calculated is then compared to maximum hot water temperature difference. If the calculated one is bigger than the maximum, the model sets the hot water temperature difference equal the maximum, and then adjusts the condenser heat transfer rate at the given conditions as follows:

˙QCondhtg=˙mavail,hwCp,hwΔThw,max

Finally, the total heat transfer energy by the evaporator and condenser can then be calculated as follows:

QEvaphtg=˙QEvaphtg⋅TimeStepSys⋅3600

QCondhtg=˙QCondhtg⋅TimeStepSys⋅3600

Central Geothermal Systems, Applications Engineering Manual, Trane Company, April 2010, SYS-APM009-EN.

## ChillerHeaterPerformance:Electric:EIR[LINK]

## Overview[LINK]

The object simulates the performance of a chiller-heater which can receive pre-cooled or pre-heated water from the source loop, and provide cooling, heating, or simultaneous cooling-heating. The object needs to work with the Central Heat Pump System object to be controlled properly. This model does not simulate the thermal performance or the power consumption of associated pumps or cooling towers. The Central Heat Pump System object holds the input/output nodes connection of the chiller-heater and its control scheme, once the chiller-heater is properly referred.

## Model Description[LINK]

The model uses user-input performance information at design conditions along with three performance curves (curve objects) for cooling capacity and efficiency to determine chiller operation at off-design conditions. Three additional performance curves for heating capacity and efficiency are used when the chiller is operating in a heating-only mode or simultaneous cooling-heating mode.

## Cooling-only mode[LINK]

The following nomenclature is used in the cooling equations:

CompMotorEfficis the compressor motor efficiencyCompPowerclgis the compressor power (W)CompPower@PLRminis the compressor power at the minimum part-load ratio (W)Cp,cw is the chilled water specific heat (J/kg-K)

CyclingRatiois the compressor cycling ratio which is equal toPLRactual/PLRminEvapCapAvailclgis the available full-load cooling capacity at current conditions (W)EvapCapFTclgis the cooling capacity function of temperature curveEIRFTclgis the electric input to cooling output factor for temperature function curveEIRFPLRclgis the electric input to cooling output factor for part-load function curve˙mcw is the chilled water mass flow rate (kg/s)

˙mavail,cw is the chilled water maximum available mass flow rate (kg/s)

PLRclgis the cooling part-load ratio which is equal toCoolingLoad/EvapCapAvailclgPLRactualis the actual part-load ratio at current conditionsPLRminis the minimum part-load ratioQCondclg is the total condenser heat transfer energy (J)

˙QCondclg is the condenser heat transfer rate (W)

QEvapclg is the total evaporator heat transfer energy (J)

˙QEvapclg is the evaporator heat transfer rate (W)

˙Qfalseloading,clg is the false loading rate (W)

RefCOPclgis the reference coefficient of performance (W/W)RefEvapCapclgis the reference evaporator capacity (W)FullLoadPwrclgis the reference full load power which is equal toEvapCapAvailclg/RefCOPclg(W)Tcondis either the entering or the leaving condenser water temperature depending on user input for the condenser water independent variable:Tcond,lif “LeavingCondenser” is chosen orTcond,eif “EnteringCondenser” is chosen.Tcond,eis the entering condenser water temperature (∘C)Tcond,lis the leaving condenser water temperature (∘C)Tcw,eis the entering chilled water temperature (W)Tcw,lis the leaving chilled water temperature (W)ΔTcw is the chilled water inlet and outlet temperature difference (∘C)

ΔTcw,max is the maximum chilled water inlet and outlet temperature difference (∘C).

The model sequentially calls each chiller-heater module in the order defined in the Central Heat Pump System object. It then determines cooling load that each chiller-heater needs to meet and water flow rates delivered to each chiller-heater. Once each chiller-heater is assumed to operate, it determines cooling capacity and efficiency using user-supplied performance information.

Three performance curves are used in the calculation of cooling capacity and efficiency as follows:

Cooling mode cooling capacity function of temperature curve (

EvapCapFTclg)Cooling mode electric input to cooling output ratio function of temperature curve (

EIRFTclg)Cooling mode electric input to cooling output ratio function of part load ratio curve (

EIRFPLRclg)The Cooling Capacity Function of Temperature Curve (

EvapCapFTclg) represents the fraction of the cooling capacity of the chiller-heater as it varies by temperature. The curve should have a value of 1.0 at the reference conditions. The output of a bi-quadratic curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:EvapCapFTclg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The Cooling Mode Electric Input to Cooling Output Ratio Function of Temperature (EIRFTclg) curve represents the fraction of electricity to the chiller-heater at full load as it varies by temperature. The output of a bi-quadratic curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EIRFTclg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The Cooling Mode Electric Input to Cooling Output Ratio Function of Part Load Ratio (EIRFPLRclg) curve represents the fraction of electricity to the chiller-heater as the load on the chiller varies at a given set of operating temperatures. The curve is normalized so that at full load the value of the curve should be 1.0. Note that the bi-cubic formulation below is generally only valid when LeavingCondenser variable is chosen for the field of Cooling Mode Condenser Water Temperature Curve Input Variable whereas the quadratic curve can be used for both choices, i.e., LeavingCondenser and EnteringCondenser. Bi-cubic may also be used when the chiller-heater uses a variable-speed compressor motor drive. The output of this curve can be determined by one of the following three performance curves:

EIRFPLRclg=a+bTCond,l+cTCond,l2+dPLR+ePLR2+fTCond,lPLR+g⋅0+hPLR3+i⋅0+j⋅0

EIRFPLRclg=a+bPLR+cPLR2

EIRFPLRclg=a+bPLR+cPLR2+dPLR3

The full-load cooling capacity at specific temperature operating conditions (i.e., at temperatures different from the design temperatures) is then computed as follows:

EvapCapAvailclg=RefEvapCapclg⋅EvapCapFTclg

The model then determines current chiller-heater’s evaporator heat transfer rate based on the total cooling load required a central heat pump system to meet and the maximum available chiller-heater cooling capacity. The maximum evaporator temperature difference between the entering chilled water temperature (

Tcw,e) and the leaving chilled water temperature (Tcw,l) obtained from the plant loop setpoint temperature can also be determined. It then calculates mass flow rate for variable flow control chiller-heaters and the temperature difference for constant flow control chiller-heaters, setting the cooling load each chiller-heater needs to meet equal to the evaporator heat transfer rate.As for variable flow control chiller-heaters, the chilled water mass flow rate is computed as follows:

˙mcw=˙QEvapclgCp,cwΔTcw,max

The chilled water mass flow rate calculated is then compared to the maximum available mass flow rate for individual chiller-heaters. If the calculated one is bigger than the maximum, the model sets the chilled water mass flow rate equal to the maximum. It then adjusts the temperature difference based on the evaporator heat transfer rate and the maximum mass flow rate. If the adjusted temperature difference also exceeds the maximum, the model finally adjusts the evaporator heat transfer rate at the maximum temperature difference and mass flow rate as follows:

˙QEvapclg=˙mavail,cwCp,cwΔTcw,max

As for constant flow control chiller-heaters, the model calculates chilled water temperature difference as follows:

ΔTcw=˙QEvapclg˙mavail,cwCp,cw

The temperature difference calculated is then compared to the maximum temperature difference allowed. If the calculated one is bigger than the maximum, the model sets the chilled water temperature difference equal the maximum, and then adjusts the evaporator heat transfer rate at the given conditions as follows:

˙QEvapclg=˙mavail,cwCp,cwΔTcw,max

The model then calculates the part-load ratio as the ratio of the evaporator heat transfer rate to the available chiller-heater capacity as follows:

PLRclg=˙QEvapclg˙QEvapAvailclg

The part-load ratio calculated is set to be between the maximum of 1.0 and the minimum of 0.0 when it is out of the range. Once the part-load ratio is calculated the cycling ratio and false loading rate can be obtained as follows:

CyclingRatio=PLRactualPLRmin

˙Qfalseloading,clg=EvapCapAvailclg⋅PLRclg⋅CyclingRatio−˙QEvapclg

The compressor power demand is then computed by:

CompPowerclg=RefFullLoadPowerclg⋅EIRFTclg⋅EIRFPLRclg⋅CyclingRatio

The heat transfer rate for the chiller-heater condenser can then be computed as follows:

˙QCondclg=ComPowerclg⋅CompMotorEffic+˙QEvapclg+˙Qfalseloading,clg

The total heat transfer energy by the evaporator and condenser can be calculated as follows:

QEvapclg=˙QEvapclg⋅TimeStepSys⋅3600

QCondclg=˙QCondclg⋅TimeStepSys⋅3600

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:

Vwater=˙mair(ωevapcond,out−ωevapcond,in)ρwater(PLR)(TimeStepSys)(3600)

where:

Vwater is the DX cooling coil evap condenser water consumption (m3)

˙mair is the evaporative condenser air mass flow rate (kg/s)

ωevapcond,in is the humidity ratio of outdoor air entering the evap condenser (kg/kg)

ωevapcond,out is the humidity ratio of air leaving the evap condenser (kg/kg)

ρwater is the density of water at the outdoor dry-bulb temperature (kg/m3)

TimeStepSysis the HVAC system simulation time step (hr).## Heating-only mode and Simultaneous cooling-heating mode[LINK]

The following nomenclature is used in the heating equations:

CompMotorEfficis the compressor motor efficiencyCompPowerhtgis the compressor power demand (W)CompPower@PLRminis the compressor power at the minimum part-load ratio (W)Cp,evap is the evaporator water specific heat (J/kg-K)

Cp,hw is the hot water specific heat (J/kg-K)

CyclingRatiois the compressor cycling ratio which is equal toPLRactual/PLRminEvapCapAvailhtgis the available full-load cooling capacity at current conditions (W)EvapCapFThtgis the heating mode cooling capacity function of temperature curveEIRFThtgis the electric input to cooling output factor for temperature function curveEIRFPLRhtgis the electric input to cooling output factor for part-load function curve˙mavail,evap is the evaporator water maximum available mass flow rate (kg/s)

˙mavail,hw is the condenser water maximum available mass flow rate (kg/s)

˙mhw is the hot water mass flow rate (kg/s)

PLRhtgis the cooling part-load ratio which is equal toRefCap/EvapCapAvailhtgPLRmaxis the maximum part-load ratio at current conditionsPLRminis the minimum part-load ratioQCondhtg is the total condenser heat transfer energy (J)

˙QCondAvailhtg is the available full-load heating capacity at current conditions (W)

˙QCondhtg is the condenser heat transfer rate (W)

QEvaphtg is the total evaporator heat transfer energy (J)

˙QEvaphtg is the evaporator heat transfer rate (W)

˙Qfalseloading,htg is the false loading rate (W)

RefCOPhtgis the reference coefficient of performance (W/W)RefEvapCaphtgis the reference evaporator capacity (W)FullLoadPwrhtgis the reference full load power which is equal toEvapCapAvailhtg/RefCOPhtg(W)Tcondis either entering or leaving condenser water temperature depending on user input for condenser water independent variable.Tcond,lif “LeavingCondenser” is chosen orTcond,eif “EnteringCondenser” is chosen.Tcond,eis the entering condenser water temperature (∘C)Tcond,lis the leaving condenser water temperature (∘C)Tcw,lis the leaving chilled water temperature (∘C)Thw,eis the entering hot water temperature (∘C)Thw,lis the leaving hot water temperature (∘C)ΔTevap is the evaporator inlet and outlet water temperature difference (∘C)

ΔThw is the hot water inlet and outlet temperature difference (∘C)

ΔThw,max is the maximum hot water inlet and outlet temperature difference (∘C).

The calculations for the evaporator side are similar to the cooling-only mode calculations. The evaporator capacity and efficiency is determined by a different set of three performance curves read in the cooling-only mode, and the performance curve set is used for both heating-only mode and simultaneous cooling-heating mode. During these modes, the evaporator side is not connected to the chilled water loop, but source water loop. The model thus assumes that each chiller-heater does not meet the plant loop chilled water setpoint temperature while the evaporator operates at the full load capacity to produce heating at a constant water flow rate.

The model sequentially calls each chiller-heater module in the order of the definition in the central heat pump system. It then determines heating load that each chiller-heater needs to meet and water flow rates delivered to each chiller-heater. Once each chiller-heater is assumed to operate, it determines heating capacity and efficiency using the following performance curves:

Heating mode cooling capacity function of temperature curve (

EvapCapFThtg)Heating mode electric input to cooling output ratio function of temperature curve (

EIRFThtg)Heating mode electric input to cooling output ratio function of part load ratio curve (

EIRFPLRhtg)The output of a Heating Mode Cooling Capacity Function of Temperature curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EvapCapFThtg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The output of a Heating Mode Cooling Output Ratio Function of Temperature curve with the input variables being the leaving chilled water temperature and either the entering or leaving condenser water temperature is given by:

EIRFThtg=a+bTcw,l+cTcw,l2+dTCond+eTCond2+fTcw,lTCond

The output of Heating Mode Cooling Output Ratio Function of Part Load Ratio curve can be determined by one of the following three performance curves as follows:

EIRFPLRhtg=a+bTCond,l+cTCond,l2+dPLR+ePLR2+fTCond,lPLR+g⋅0+hPLR3+i⋅0+j⋅0

EIRFPLRhtg=a+bPLR+cPLR2

EIRFPLRhtg=a+bPLR+cPLR2+dPLR3

The full-load evaporator capacity at specific temperature operating conditions is then given by:

EvapCapAvailhtg=RefEvapCaphtg⋅EvapCapFThtg

The part-load ratio is set to be between zero and the maximum, and the evaporator heat transfer rate is computed by:

˙QEvaphtg=EvapCapAvailhtgPLRhtg

The evaporator inlet and outlet temperature difference is then given by:

ΔTevap=˙QEvaphtg˙mavail,evapCp,evap

Once the part-load ratio is calculated the cycling ratio and false loading rate are computed by:

CyclingRatio=PLRactualPLRmin

˙Qfalseloading,htg=EvapCapAvailhtg⋅PLRhtg⋅CyclingRatio−˙QEvaphtg

The compressor power demand is then computed by:

CompPowerhtg=FullLoadPowerhtg⋅EIRFThtg⋅EIRFPLRhtg⋅CyclingRatio

The heat transfer rate of the chiller-heater condenser is then computed as follows:

˙QCondAvailhtg=ComPowerhtg⋅CompMotorEffic+˙QEvaphtg+˙Qfalseloading,htg

Once condenser available heating capacity is determined, the model calculates current chiller-heater’s condenser heat transfer rate based on the total heating load required a central heat pump system to meet as well as available heating capacity of the chiller-heater. The maximum condenser temperature difference between the entering hot water temperature (

Thw,e) and the leaving hot water temperature (Thw,l) obtained from the plant loop setpoint temperature can also be obtained. It then calculates condenser water mass flow rate for variable flow control chiller-heaters and the hot water temperature difference for constant flow control chiller-heaters, setting the cooling load that each chiller-heater needs to meet equal the evaporator heat transfer rate.As for variable flow control chiller-heaters, the condenser water mass flow rate is computed as follows:

˙mhw=˙QCondhtgCp,hwΔThw,max

The condenser water mass flow rate calculated is then compared to the maximum available mass flow rate for individual chiller-heaters. If the calculated one is bigger than the maximum, the model sets the condenser water mass flow rate equal the maximum. It then adjusts the hot water temperature difference at the maximum mass flow rate. If the adjusted temperature difference also exceeds the maximum, the model finally adjusts the condenser heat transfer rate at the maximum allowable conditions as follows:

˙QCondhtg=˙mavail,hwCp,hwΔThw,max

As for constant flow control chiller-heaters, the model calculates condenser temperature difference as follows:

ΔThw=˙QCondhtg˙mavail,hwCp,hw

The temperature difference calculated is then compared to maximum hot water temperature difference. If the calculated one is bigger than the maximum, the model sets the hot water temperature difference equal the maximum, and then adjusts the condenser heat transfer rate at the given conditions as follows:

˙QCondhtg=˙mavail,hwCp,hwΔThw,max

Finally, the total heat transfer energy by the evaporator and condenser can then be calculated as follows:

QEvaphtg=˙QEvaphtg⋅TimeStepSys⋅3600

QCondhtg=˙QCondhtg⋅TimeStepSys⋅3600

## References[LINK]

Central Geothermal Systems, Applications Engineering Manual, Trane Company, April 2010, SYS-APM009-EN.

Documentation content copyright © 1996-2021 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.

This documentation is made available under the EnergyPlus Open Source License v1.0.