Grey-box model for reversible chillers using a black-box to
simulate the refrigeration cycle
Deprecation Warning
This model is deprecated, switch to the ModularReversible
approach. Read
https://github.com/RWTH-EBC/AixLib/wiki/Switch-from-model-HeatPump-to-ModularReversible
for further information.
This generic grey-box chiller model uses empirical data to model
the refrigerant cycle. The modelling of system inertias and heat
losses allow the simulation of transient states.
Resulting in the choosen model structure, several configurations
are possible:
- Compressor type: on/off or inverter controlled
- Reversible chiller / only cooling
- Source/Sink: Any combination of mediums is possible
- Generik: Losses and inertias can be switched on or off.
Concept
Using a signal bus as a connector, this chiller model can be
easily combined within a chiller system model including several
control or safety blocks analogous to AixLib.Controls.HeatPump.
The relevant data is aggregated. The mode signal chooses the type
of the chiller operation. As a result, this model can also be used
as a heat pump:
- mode = true: Chilling
- mode = false: Heating
To model both on/off and inverter controlled chillers, the
compressor speed is normalizd to a relative value between 0 and
1.
Possible icing of the evaporator is modelled with an input value
between 0 and 1.
The model structure is as follows. To understand each submodel,
please have a look at the corresponding model information:
- InnerCycle
(Black Box): Here, the user can use between several input models or
just easily create his own, modular black box model. Please look at
the model description for more info.
- Inertia: A n-order element is used to model system inertias
(mass and thermal) of components inside the refrigerant cycle
(compressor, pipes, expansion valve)
-
HeatExchanger: This new model also enable modelling of thermal
interias and heat losses in a heat exchanger. Please look at the
model description for more info.
Parametrization
To simplify the parametrization of the evaporator and condenser
volumes and nominal mass flows there exists an option of automatic
estimation based on the nominal usable cooling power of the
Chiller. This function uses a linear correlation of these
parameters, which was established from the linear regression of
more than 20 data sets of water-to-water chillers from different
manufacturers (e.g. Carrier, Trane, Lennox) ranging from about 25kW
to 1MW nominal power. The linear regressions with coefficients of
determination above 91% give a good approximation of these
parameters. Nevertheless, estimates for machines outside the given
range should be checked for plausibility during simulation.
Assumptions
Several assumptions where made in order to model the chiller.
For a detailed description see the corresponding model.
-
Performance data 2D: In order to model inverter controlled
chillers, the compressor speed is scaled linearly
-
Performance data 2D: Reduced evaporator power as a result of
icing. The icing factor is multiplied with the evaporator
power.
- Inertia: The default value of the n-th order element is
set to 3. This follows comparisons with experimental data.
- Scaling factor: A scaling facor is implemented for
scaling of the chiller power and capacity. The factor scales the
parameters V, m_flow_nominal, C, GIns, GOut and dp_nominal. As a
result, the chiller can supply more heat with the COP staying
nearly constant. However, one has to make sure that the supplied
pressure difference or mass flow is also scaled with this factor,
as the nominal values do not increase said mass flow.
Known Limitations
- The n-th order element has a big influence on computational
time. Reducing the order or disabling it completly will decrease
computational time.
- Reversing the mode: A normal 4-way-exchange valve suffers from
heat losses and irreversibilities due to switching from one mode to
another. Theses losses are not taken into account.
Contents
| Name |
Description |
PerDataMainChi |
Performance data of a chiller in main operation mode |
PerDataRevChi |
Performance data of a chiller in reversible operation mode |
- May 22, 2019 by Julian Matthes:
First implementation (see issue #715)
Generated at 2026-04-15T18:18:36Z by OpenModelicaOpenModelica 1.26.3 using
GenerateDoc.mos