Buildings.Fluid.Chillers

Package with chiller models

Information

This package contains component models for chillers.

Extends from Modelica.Icons.VariantsPackage (Icon for package containing variants).

Package Content

Name	Description
AbsorptionIndirectSteam	Indirect steam heated absorption chiller based on performance curves
Carnot_TEva	Chiller with prescribed evaporator leaving temperature and performance curve adjusted based on Carnot efficiency
Carnot_y	Chiller with performance curve adjusted based on Carnot efficiency
ElectricEIR	Electric chiller based on the DOE-2.1 model
ElectricReformulatedEIR	Electric chiller based on the DOE-2.1 model, but with performance as a function of condenser leaving instead of entering temperature
Data	Performance data for electric chillers
Examples	Collection of models that illustrate model use and test models
Validation	Collection of models that validate the chiller models
BaseClasses	Package with base classes for Buildings.Fluid.Chillers

Buildings.Fluid.Chillers.AbsorptionIndirectSteam

Indirect steam heated absorption chiller based on performance curves

Information

Model for an indirect steam heated absorption chiller based on performance curves. The model uses performance curves similar to the EnergyPlus model Chiller:Absorption:Indirect.

The model uses six functions to predict the chiller cooling capacity, power consumption for the chiller pump and the generator heat flow rate and the condenser heat flow. These functions use the performance data stored in the record per. The computations are as follows:

The capacity function of the evaporator is

capFun_eva = A₁ + A₂ T_eva,lvg + A₃ T²_eva,lvg + A₄ T³_eva,lvg.

The capacity function of the condenser is

capFun_con = B₁ + B₂ T_con,ent + B₃ T²_con,ent + B₄ T³_con,ent.

These capacity functions are used to compute the available cooling capacity of the evaporator as

Q̇_eva,ava = capFun_eva   capFun_con   Q̇_eva,0,

where Q̇_eva,0 is obtained from the performance data per.QEva_flow_nominal. Let Q̇_eva,set denote the heat required to meet the set point TSet. Then, the model computes the part load ratio as

PLR =min(Q̇_eva,set/Q̇_eva,ava, PLR_max).

Hence, the model ensures that the chiller capacity does not exceed the chiller capacity specified by the parameter per.PLRMax. The cycling ratio is computed as

CR = min(PLR/PLR_min, 1.0),

where PRL_min is obtained from the performance record per.PLRMin. This ratio expresses the fraction of time that a chiller would run if it were to cycle because its load is smaller than the minimal load at which it can operate. Note that this model continuously operates even if the part load ratio is below the minimum part load ratio. Its leaving evaporator and condenser temperature can therefore be considered as an average temperature between the modes when the compressor is off and on.

Using the part load ratio, the energy input ratio of the chiller pump is

EIRP = C₁ + C₂PLR+C₃PLR².

The generator heat input ratio is

genHIR = D₁ + D₂PLR+D₃PLR²+D₄PLR³.

Two additional curves modify the heat input requirement based on the condenser inlet water temperature and the evaporator outlet water temperature. Specifically, the generator heat modifier based on the condenser inlet water temperature is

genT_con = E₁ + E₂ T_con,ent + E₃ T²_con,ent + E₄ T³_con,ent,

and the generator heat modifier based on the evaporator inlet water temperature is

genT_eva= F₁ + F₂ T_eva,lvg + F₃ T²_eva,lvg + F₄ T³_eva,lvg.

The main outputs of the model that are to be used in energy analysis are the required generator heat QGen_flow and the electric power consumption of the chiller pump P. For example, if the chiller were to be regenerated with steam, then QGen_flow is the heat that must be provided by a steam loop. This model computes the required generator heat as

Q̇_gen = -Q̇_eva,ava genHIR genT_con genT_eva CR.

The pump power consumption is

P = EIRP CR P₀,

where P₀ is the pump nominal power obtained from the performance data per.P_nominal. The heat balance of the chiller is

Q̇_con = -Q̇_eva + Q̇_gen + P.

Performance data

The equipment performance data is obtained from the record per, which is an instance of Buildings.Fluid.Chillers.Data.AbsorptionIndirectSteam. Additional performance curves can be developed using two available techniques (Hydeman and Gillespie, 2002). The first technique is called the Least-squares Linear Regression method and is used when sufficient performance data exist to employ standard least-square linear regression techniques. The second technique is called Reference Curve Method and is used when insufficient performance data exist to apply linear regression techniques. A detailed description of both techniques can be found in Hydeman and Gillespie (2002).

References

• Hydeman, M. and K.L. Gillespie. 2002. Tools and Techniques to Calibrate Electric Chiller Component Models. ASHRAE Transactions, AC-02-9-1.

Extends from Buildings.Fluid.Interfaces.FourPortHeatMassExchanger (Model transporting two fluid streams between four ports with storing mass or energy).

Parameters

Connectors

Modelica definition

Buildings.Fluid.Chillers.Carnot_TEva

Chiller with prescribed evaporator leaving temperature and performance curve adjusted based on Carnot efficiency

Information

This is a model of a chiller whose coefficient of performance COP changes with temperatures in the same way as the Carnot efficiency changes. The control input is the setpoint of the evaporator leaving temperature, which is met exactly at steady state if the chiller has sufficient capacity.

Set use_eta_Carnot_nominal=true to specify directly the Carnot effectiveness η_Carnot,0, in which case the value of the parameter COP_nominal will not affect the simulation. If use_eta_Carnot_nominal=false, the model will use the value of the parameter COP_nominal together with the specified nominal temperatures to compute the Carnot effectiveness as

η_Carnot,0 = COP₀ ⁄ (T_eva,0 ⁄ (T_con,0 + T_app,con,0 - (T_eva,0-T_app,eva,0))),

where T_eva,0 is the evaporator temperature, T_con,0 is the condenser temperature, T_app,eva,0 is the evaporator approach temperature and T_app,con,0 is the condenser approach temperature.

The COP is computed as the product

COP = η_Carnot,0 COP_Carnot η_PL,

where COP_Carnot is the Carnot efficiency and η_PL is the part load efficiency, expressed using a polynomial. This polynomial has the form

η_PL = a₁ + a₂ y + a₃ y² + ...,

where y ∈ [0, 1] is the part load for cooling and the coefficients a_i are declared by the parameter a.

On the Dynamics tag, the model can be parametrized to compute a transient or steady-state response. The transient response of the model is computed using a first order differential equation for the evaporator and condenser fluid volumes. The chiller outlet temperatures are equal to the temperatures of these lumped volumes.

Typical use and important parameters

When using this component, make sure that the condenser has sufficient mass flow rate. Based on the evaporator mass flow rate, temperature difference and the efficiencies, the model computes how much heat will be added to the condenser. If the mass flow rate is too small, very high outlet temperatures can result.

The evaporator heat flow rate QEva_flow_nominal is used to assign the default value for the mass flow rates, which are used for the pressure drop calculations. It is also used to compute the part load efficiency. Hence, make sure that QEva_flow_nominal is set to a reasonable value.

The maximum cooling capacity is set by the parameter QEva_flow_min, which is by default set to negative infinity.

The coefficient of performance depends on the evaporator and condenser leaving temperature since otherwise the second law of thermodynamics may be violated.

Notes

For a similar model that can be used as a heat pump, see Buildings.Fluid.HeatPumps.Examples.Carnot_TCon.

Extends from Buildings.Fluid.Chillers.BaseClasses.PartialCarnot_T (Partial model for chiller with performance curve adjusted based on Carnot efficiency).

Parameters

Connectors

Modelica definition

Buildings.Fluid.Chillers.Carnot_y

Chiller with performance curve adjusted based on Carnot efficiency

Information

This is a model of a chiller whose coefficient of performance COP changes with temperatures in the same way as the Carnot efficiency changes. The input signal y is the control signal for the compressor.

Set use_eta_Carnot_nominal=true to specify directly the Carnot effectiveness η_Carnot,0, in which case the value of the parameter COP_nominal will not affect the simulation. If use_eta_Carnot_nominal=false, the model will use the value of the parameter COP_nominal together with the specified nominal temperatures to compute the Carnot effectiveness as

η_Carnot,0 = COP₀ ⁄ (T_eva,0 ⁄ (T_con,0 + T_app,con,0 - (T_eva,0-T_app,eva,0))),

where T_eva,0 is the evaporator temperature, T_con,0 is the condenser temperature, T_app,eva,0 is the evaporator approach temperature and T_app,con,0 is the condenser approach temperature.

The COP is computed as the product

COP = η_Carnot,0 COP_Carnot η_PL,

where COP_Carnot is the Carnot efficiency and η_PL is the part load efficiency, expressed using a polynomial. This polynomial has the form

η_PL = a₁ + a₂ y + a₃ y² + ...,

where y ∈ [0, 1] is the part load for cooling and the coefficients a_i are declared by the parameter a.

On the Dynamics tag, the model can be parametrized to compute a transient or steady-state response. The transient response of the model is computed using a first order differential equation for the evaporator and condenser fluid volumes. The chiller outlet temperatures are equal to the temperatures of these lumped volumes.

Typical use and important parameters

When using this component, make sure that the evaporator and the condenser have sufficient mass flow rate. Based on the mass flow rates, the compressor power, temperature difference and the efficiencies, the model computes how much heat will be added to the condenser and removed at the evaporator. If the mass flow rates are too small, very high temperature differences can result.

The evaporator heat flow rate QEva_flow_nominal is used to assign the default value for the mass flow rates, which are used for the pressure drop calculations. It is also used to compute the part load efficiency. Hence, make sure that QEva_flow_nominal is set to a reasonable value.

The maximum cooling capacity is set by the parameter QEva_flow_min, which is by default set to negative infinity.

The coefficient of performance depends on the evaporator and condenser leaving temperature since otherwise the second law of thermodynamics may be violated.

Notes

For a similar model that can be used as a heat pump, see Buildings.Fluid.HeatPumps.Carnot_y.

Extends from Buildings.Fluid.Chillers.BaseClasses.PartialCarnot_y (Partial chiller model with performance curve adjusted based on Carnot efficiency).

Parameters

Connectors

Modelica definition

Buildings.Fluid.Chillers.ElectricEIR

Electric chiller based on the DOE-2.1 model

Information

Model of an electric chiller, based on the DOE-2.1 chiller model and the EnergyPlus chiller model Chiller:Electric:EIR.

This model uses three functions to predict capacity and power consumption:

These functions take the same form as documented in EnergyPlus v22.1.0 Engineering Reference section 14.3.9.2 (equations 14.231 through 14.233). These curves are stored in the data record per and are available from Buildings.Fluid.Chillers.Data.ElectricEIR. How they are used to compute the adjusted capacity and compressor power can be found in the documentation of Buildings.Fluid.Chillers.BaseClasses.PartialElectric. Additional performance curves can be developed using two available techniques (Hydeman and Gillespie, 2002). The first technique is called the Least-squares Linear Regression method and is used when sufficient performance data exist to employ standard least-square linear regression techniques. The second technique is called Reference Curve Method and is used when insufficient performance data exist to apply linear regression techniques. A detailed description of both techniques can be found in Hydeman and Gillespie (2002).

The model takes as an input the set point for the leaving chilled water temperature, which is met if the chiller has sufficient capacity. Thus, the model has a built-in, ideal temperature control. The model has three tests on the part load ratio and the cycling ratio:

1. The test

     PLR1 =min(QEva_flow_set/QEva_flow_ava, per.PLRMax);
   

   ensures that the chiller capacity does not exceed the chiller capacity specified by the parameter per.PLRMax.
2. The test

     CR = min(PLR1/per.PRLMin, 1.0);
   

   computes a cycling ratio. This ratio expresses the fraction of time that a chiller would run if it were to cycle because its load is smaller than the minimal load at which it can operate. Note that this model continuously operates even if the part load ratio is below the minimum part load ratio. Its leaving evaporator and condenser temperature can therefore be considered as an average temperature between the modes where the compressor is off and on.
3. The test

     PLR2 = max(per.PLRMinUnl, PLR1);
   

   computes the part load ratio of the compressor. The assumption is that for a part load ratio below per.PLRMinUnl, the chiller uses hot gas bypass to reduce the capacity, while the compressor power draw does not change.

The electric power only contains the power for the compressor, but not any power for pumps or fans.

The model can be parametrized to compute a transient or steady-state response. The transient response of the chiller is computed using a first order differential equation for the evaporator and condenser fluid volumes. The chiller outlet temperatures are equal to the temperatures of these lumped volumes.

Optionally, the model can be configured to represent heat recovery chillers with a switchover option by setting the parameter have_switchover to true. In that case an additional Boolean input connector coo is used. The chiller is tracking a chilled water supply temperature setpoint at the outlet of the evaporator barrel if coo is true. Otherwise, if coo is false, the chiller is tracking a hot water supply temperature setpoint at the outlet of the condenser barrel. See Buildings.Fluid.Chillers.Examples.ElectricEIR_HeatRecovery for an example with a chiller operating in heating mode.

References

• Hydeman, M. and K.L. Gillespie. 2002. Tools and Techniques to Calibrate Electric Chiller Component Models. ASHRAE Transactions, AC-02-9-1.
• EnergyPlus v22.1.0 Engineering Reference

Extends from Buildings.Fluid.Chillers.BaseClasses.PartialElectric (Partial model for electric chiller based on the model in DOE-2, CoolTools and EnergyPlus).

Parameters

Connectors

Modelica definition

Buildings.Fluid.Chillers.ElectricReformulatedEIR

Electric chiller based on the DOE-2.1 model, but with performance as a function of condenser leaving instead of entering temperature

Information

Model of an electric chiller, based on the model by Hydeman et al. (2002) that has been developed in the CoolTools project and that is implemented in EnergyPlus as the model Chiller:Electric:ReformulatedEIR. This empirical model is similar to Buildings.Fluid.Chillers.ElectricEIR. The difference is that to compute the performance, this model uses the condenser leaving temperature instead of the entering temperature, and it uses a bicubic polynomial to compute the part load performance.

This model uses three functions to predict capacity and power consumption:

These curves are stored in the data record per and are available from Buildings.Fluid.Chillers.Data.ElectricReformulatedEIR. How they are used to compute the adjusted capacity and compressor power can be found in the documentation of Buildings.Fluid.Chillers.BaseClasses.PartialElectric. Additional performance curves can be developed using two available techniques (Hydeman and Gillespie, 2002). The first technique is called the Least-squares Linear Regression method and is used when sufficient performance data exist to employ standard least-square linear regression techniques. The second technique is called Reference Curve Method and is used when insufficient performance data exist to apply linear regression techniques. A detailed description of both techniques can be found in Hydeman and Gillespie (2002).

The model takes as an input the set point for the leaving chilled water temperature, which is met if the chiller has sufficient capacity. Thus, the model has a built-in, ideal temperature control. The model has three tests on the part load ratio and the cycling ratio:

1. The test

     PLR1 =min(QEva_flow_set/QEva_flow_ava, per.PLRMax);
   

   ensures that the chiller capacity does not exceed the chiller capacity specified by the parameter per.PLRMax.
2. The test

     CR = min(PLR1/per.PRLMin, 1.0);
   

   computes a cycling ratio. This ratio expresses the fraction of time that a chiller would run if it were to cycle because its load is smaller than the minimal load at which it can operate. Note that this model continuously operates even if the part load ratio is below the minimum part load ratio. Its leaving evaporator and condenser temperature can therefore be considered as an average temperature between the modes where the compressor is off and on.
3. The test

     PLR2 = max(per.PLRMinUnl, PLR1);
   

   computes the part load ratio of the compressor. The assumption is that for a part load ratio below per.PLRMinUnl, the chiller uses hot gas bypass to reduce the capacity, while the compressor power draw does not change.

The electric power only contains the power for the compressor, but not any power for pumps or fans.

The model can be parametrized to compute a transient or steady-state response. The transient response of the chiller is computed using a first order differential equation for the evaporator and condenser fluid volumes. The chiller outlet temperatures are equal to the temperatures of these lumped volumes.

Optionally, the model can be configured to represent heat recovery chillers with a switchover option by setting the parameter have_switchover to true. In that case an additional Boolean input connector coo is used. The chiller is tracking a chilled water supply temperature setpoint at the outlet of the evaporator barrel if coo is true. Otherwise, if coo is false, the chiller is tracking a hot water supply temperature setpoint at the outlet of the condenser barrel. See Buildings.Fluid.Chillers.Examples.ElectricEIR_HeatRecovery for an example with a chiller operating in heating mode.

References

• Hydeman, M., N. Webb, P. Sreedharan, and S. Blanc. 2002. Development and Testing of a Reformulated Regression-Based Electric Chiller Model. ASHRAE Transactions, HI-02-18-2.
• Hydeman, M. and K.L. Gillespie. 2002. Tools and Techniques to Calibrate Electric Chiller Component Models. ASHRAE Transactions, AC-02-9-1.

Extends from Buildings.Fluid.Chillers.BaseClasses.PartialElectric (Partial model for electric chiller based on the model in DOE-2, CoolTools and EnergyPlus).

Parameters

Connectors

Modelica definition