Buildings.Fluid.HeatPumps.Examples
Collection of models that illustrate model use and test models
Information
This package contains examples for the use of models that can be found in Buildings.Fluid.HeatPumps.
Extends from Modelica.Icons.ExamplesPackage (Icon for packages containing runnable examples).
Package Content
| Name | Description |
|---|---|
 Carnot_TCon
|
Test model for heat pump based on Carnot efficiency and condenser outlet temperature control signal |
| Carnot_y
|
Test model for heat pump based on Carnot efficiency |
| ScrollWaterToWater_OneRoomRadiator
|
Heat pump with scroll compressor connected to a simple room model with radiator |
Buildings.Fluid.HeatPumps.Examples.Carnot_TCon
Test model for heat pump based on Carnot efficiency and condenser outlet temperature control signal
Information
Example that simulates a chiller whose efficiency is scaled based on the
Carnot cycle.
The chiller takes as an input the evaporator leaving water temperature.
The condenser mass flow rate is computed in such a way that it has
a temperature difference equal to dTEva_nominal.
Extends from Modelica.Icons.Example (Icon for runnable examples).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| TemperatureDifference | dTEva_nominal | -5 | Temperature difference evaporator inlet-outlet [K] |
| TemperatureDifference | dTCon_nominal | 10 | Temperature difference condenser outlet-inlet [K] |
| HeatFlowRate | QCon_flow_nominal | 100E3 | Evaporator heat flow rate [W] |
| MassFlowRate | m1_flow_nominal | QCon_flow_nominal/dTCon_nomi... | Nominal mass flow rate at condenser [kg/s] |
Modelica definition
model Carnot_TCon
"Test model for heat pump based on Carnot efficiency and condenser outlet temperature control signal"
extends Modelica.Icons.Example;
package Medium1 = Buildings.Media.Water "Medium model";
package Medium2 = Buildings.Media.Water "Medium model";
parameter Modelica.SIunits.TemperatureDifference dTEva_nominal=-5
"Temperature difference evaporator inlet-outlet";
parameter Modelica.SIunits.TemperatureDifference dTCon_nominal=10
"Temperature difference condenser outlet-inlet";
parameter Modelica.SIunits.HeatFlowRate QCon_flow_nominal=100E3
"Evaporator heat flow rate";
parameter Modelica.SIunits.MassFlowRate m1_flow_nominal=QCon_flow_nominal/
dTCon_nominal/4200 "Nominal mass flow rate at condenser";
Buildings.Fluid.HeatPumps.Carnot_TCon heaPum(
redeclare package Medium1 = Medium1,
redeclare package Medium2 = Medium2,
dTEva_nominal=dTEva_nominal,
dTCon_nominal=dTCon_nominal,
m1_flow_nominal=m1_flow_nominal,
show_T=true,
allowFlowReversal1=false,
allowFlowReversal2=false,
use_eta_Carnot_nominal=true,
etaCarnot_nominal=0.3,
QCon_flow_nominal=QCon_flow_nominal,
dp1_nominal=6000,
dp2_nominal=6000) "Heat pump";
Buildings.Fluid.Sources.MassFlowSource_T sou1(
nPorts=1,
redeclare package Medium = Medium1,
m_flow=m1_flow_nominal,
T=293.15);
Buildings.Fluid.Sources.MassFlowSource_T sou2(
nPorts=1,
redeclare package Medium = Medium2,
use_T_in=false,
use_m_flow_in=true,
T=288.15);
Buildings.Fluid.Sources.FixedBoundary sin1(redeclare package Medium = Medium1,
nPorts=1);
Buildings.Fluid.Sources.FixedBoundary sin2(nPorts=1, redeclare package Medium
= Medium2);
Modelica.Blocks.Sources.Ramp TConLvg(
duration=60,
startTime=1800,
height=15,
offset=273.15 + 35) "Control signal for condenser leaving temperature";
Modelica.Blocks.Math.Gain mEva_flow(k=-1/cp2_default/dTEva_nominal)
"Evaporator mass flow rate";
Modelica.Blocks.Math.Add QEva_flow(k2=-1) "Evaporator heat flow rate";
final parameter Modelica.SIunits.SpecificHeatCapacity cp2_default=
Medium2.specificHeatCapacityCp(Medium2.setState_pTX(
Medium2.p_default,
Medium2.T_default,
Medium2.X_default))
"Specific heat capacity of medium 2 at default medium state";
equation
connect(sou1.ports[1], heaPum.port_a1);
connect(sou2.ports[1], heaPum.port_a2);
connect(sin2.ports[1], heaPum.port_b2);
connect(QEva_flow.y, mEva_flow.u);
connect(heaPum.port_b1, sin1.ports[1]);
connect(TConLvg.y, heaPum.TSet);
connect(mEva_flow.y, sou2.m_flow_in);
connect(QEva_flow.u1, heaPum.QCon_flow);
connect(QEva_flow.u2, heaPum.P);
end Carnot_TCon;
Buildings.Fluid.HeatPumps.Examples.Carnot_y
Test model for heat pump based on Carnot efficiency
Information
Example that simulates a heat pump whose efficiency is scaled based on the Carnot cycle. The control signal of the heat pump is the compressor speed.
Extends from Modelica.Icons.Example (Icon for runnable examples).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| Real | COP_nominal | 6 | Nominal COP |
| Power | P_nominal | 10E3 | Nominal compressor power (at y=1) [W] |
| TemperatureDifference | dTEva_nominal | -10 | Temperature difference evaporator outlet-inlet [K] |
| TemperatureDifference | dTCon_nominal | 10 | Temperature difference condenser outlet-inlet [K] |
| MassFlowRate | m2_flow_nominal | -P_nominal*(COP_nominal - 1)... | Nominal mass flow rate at chilled water side [kg/s] |
| MassFlowRate | m1_flow_nominal | P_nominal*COP_nominal/cp1_de... | Nominal mass flow rate at condenser water wide [kg/s] |
Modelica definition
model Carnot_y "Test model for heat pump based on Carnot efficiency"
extends Modelica.Icons.Example;
package Medium1 = Buildings.Media.Water "Medium model";
package Medium2 = Buildings.Media.Water "Medium model";
parameter Real COP_nominal=6 "Nominal COP";
parameter Modelica.SIunits.Power P_nominal=10E3
"Nominal compressor power (at y=1)";
parameter Modelica.SIunits.TemperatureDifference dTEva_nominal=-10
"Temperature difference evaporator outlet-inlet";
parameter Modelica.SIunits.TemperatureDifference dTCon_nominal=10
"Temperature difference condenser outlet-inlet";
parameter Modelica.SIunits.MassFlowRate m2_flow_nominal=-P_nominal*(
COP_nominal - 1)/cp2_default/dTEva_nominal
"Nominal mass flow rate at chilled water side";
parameter Modelica.SIunits.MassFlowRate m1_flow_nominal=P_nominal*COP_nominal/
cp1_default/dTCon_nominal
"Nominal mass flow rate at condenser water wide";
Buildings.Fluid.HeatPumps.Carnot_y heaPum(
redeclare package Medium1 = Medium1,
redeclare package Medium2 = Medium2,
P_nominal=P_nominal,
dTEva_nominal=dTEva_nominal,
dTCon_nominal=dTCon_nominal,
dp1_nominal=6000,
dp2_nominal=6000,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
show_T=true,
use_eta_Carnot_nominal=false,
COP_nominal=COP_nominal,
TCon_nominal=303.15,
TEva_nominal=278.15) "Heat pump model";
Buildings.Fluid.Sources.MassFlowSource_T sou1(
nPorts=1,
redeclare package Medium = Medium1,
use_T_in=true,
m_flow=m1_flow_nominal,
T=298.15);
Buildings.Fluid.Sources.MassFlowSource_T sou2(
nPorts=1,
redeclare package Medium = Medium2,
use_T_in=true,
m_flow=m2_flow_nominal,
T=291.15);
Buildings.Fluid.Sources.FixedBoundary sin1(nPorts=1, redeclare package Medium
= Medium1);
Buildings.Fluid.Sources.FixedBoundary sin2(nPorts=1, redeclare package Medium
= Medium2);
Modelica.Blocks.Sources.Ramp uCom(
height=-1,
duration=60,
offset=1,
startTime=1800) "Compressor control signal";
Modelica.Blocks.Sources.Ramp TCon_in(
height=10,
duration=60,
offset=273.15 + 20,
startTime=60) "Condenser inlet temperature";
Modelica.Blocks.Sources.Ramp TEva_in(
height=10,
duration=60,
startTime=900,
offset=273.15 + 15) "Evaporator inlet temperature";
final parameter Modelica.SIunits.SpecificHeatCapacity cp1_default=
Medium1.specificHeatCapacityCp(Medium1.setState_pTX(
Medium1.p_default,
Medium1.T_default,
Medium1.X_default))
"Specific heat capacity of medium 2 at default medium state";
final parameter Modelica.SIunits.SpecificHeatCapacity cp2_default=
Medium2.specificHeatCapacityCp(Medium2.setState_pTX(
Medium2.p_default,
Medium2.T_default,
Medium2.X_default))
"Specific heat capacity of medium 2 at default medium state";
equation
connect(sou1.ports[1], heaPum.port_a1);
connect(sou2.ports[1], heaPum.port_a2);
connect(heaPum.port_b1, sin1.ports[1]);
connect(sin2.ports[1], heaPum.port_b2);
connect(TCon_in.y, sou1.T_in);
connect(TEva_in.y, sou2.T_in);
connect(uCom.y, heaPum.y);
end Carnot_y;
Buildings.Fluid.HeatPumps.Examples.ScrollWaterToWater_OneRoomRadiator
Heat pump with scroll compressor connected to a simple room model with radiator
Information
Example that simulates one room equipped with a radiator. Hot water is produced by a 24 kW nominal capacity heat pump. The source side water temperature to the heat pump is constant at 10°C.
The heat pump is turned on when the room temperature falls below 19°C and turned off when the room temperature rises above 21°C.
Extends from Modelica.Icons.Example (Icon for runnable examples).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package MediumA | Buildings.Media.Air | Medium model for air | |
| replaceable package MediumW | Buildings.Media.Water | Medium model for water | |
| HeatFlowRate | Q_flow_nominal | 20000 | Nominal heat flow rate of radiator [W] |
| Temperature | TRadSup_nominal | 273.15 + 50 | Radiator nominal supply water temperature [K] |
| Temperature | TRadRet_nominal | 273.15 + 45 | Radiator nominal return water temperature [K] |
| MassFlowRate | mHeaPum_flow_nominal | Q_flow_nominal/4200/5 | Heat pump nominal mass flow rate [kg/s] |
| Volume | V | 6*10*3 | Room volume [m3] |
| MassFlowRate | mA_flow_nominal | V*6/3600 | Nominal mass flow rate [kg/s] |
| HeatFlowRate | QRooInt_flow | 4000 | Internal heat gains of the room [W] |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package MediumA | Medium model for air | |
| replaceable package MediumW | Medium model for water | |
| Bus | weaBus | Weather data bus |
Modelica definition
model ScrollWaterToWater_OneRoomRadiator
"Heat pump with scroll compressor connected to a simple room model with radiator"
extends Modelica.Icons.Example;
replaceable package MediumA = Buildings.Media.Air "Medium model for air";
replaceable package MediumW = Buildings.Media.Water "Medium model for water";
parameter Modelica.SIunits.HeatFlowRate Q_flow_nominal=20000
"Nominal heat flow rate of radiator";
parameter Modelica.SIunits.Temperature TRadSup_nominal=273.15 + 50
"Radiator nominal supply water temperature";
parameter Modelica.SIunits.Temperature TRadRet_nominal=273.15 + 45
"Radiator nominal return water temperature";
parameter Modelica.SIunits.MassFlowRate mHeaPum_flow_nominal=Q_flow_nominal/4200
/5 "Heat pump nominal mass flow rate";
//------------------------------------------------------------------------------//
Buildings.Fluid.MixingVolumes.MixingVolume vol(
redeclare package Medium = MediumA,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
m_flow_nominal=mA_flow_nominal,
V=V);
Modelica.Thermal.HeatTransfer.Components.ThermalConductor theCon(G=20000/40)
"Thermal conductance with the ambient";
parameter Modelica.SIunits.Volume V=6*10*3 "Room volume";
parameter Modelica.SIunits.MassFlowRate mA_flow_nominal=V*6/3600
"Nominal mass flow rate";
parameter Modelica.SIunits.HeatFlowRate QRooInt_flow=4000
"Internal heat gains of the room";
Modelica.Thermal.HeatTransfer.Sources.PrescribedHeatFlow preHea
"Prescribed heat flow";
Modelica.Thermal.HeatTransfer.Components.HeatCapacitor heaCap(C=2*V*1.2*1006)
"Heat capacity for furniture and walls";
Modelica.Blocks.Sources.CombiTimeTable timTab(
extrapolation=Modelica.Blocks.Types.Extrapolation.Periodic,
smoothness=Modelica.Blocks.Types.Smoothness.ConstantSegments,
table=[-6*3600, 0; 8*3600, QRooInt_flow; 18*3600, 0])
"Time table for internal heat gain";
Buildings.Fluid.HeatExchangers.Radiators.RadiatorEN442_2 rad(
redeclare package Medium = MediumW,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
Q_flow_nominal=Q_flow_nominal,
T_a_nominal=TRadSup_nominal,
T_b_nominal=TRadRet_nominal,
m_flow_nominal=mHeaPum_flow_nominal,
T_start=TRadSup_nominal) "Radiator";
Buildings.Fluid.Sensors.TemperatureTwoPort temSup(
redeclare package Medium = MediumW,
m_flow_nominal=mHeaPum_flow_nominal,
T_start=TRadSup_nominal) "Supply water temperature";
Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensor temRoo
"Room temperature";
//----------------------------------------------------------------------------//
Buildings.Fluid.Movers.FlowControlled_m_flow pumHeaPum(
redeclare package Medium = MediumW,
m_flow_nominal=mHeaPum_flow_nominal,
y_start=1,
m_flow_start=0.85,
T_start=TRadSup_nominal,
nominalValuesDefineDefaultPressureCurve=true,
use_inputFilter=false,
energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState)
"Pump for radiator side";
//----------------------------------------------------------------------------//
Buildings.Fluid.Sources.FixedBoundary preSou(
redeclare package Medium = MediumW,
nPorts=1,
T=TRadSup_nominal)
"Source for pressure and to account for thermal expansion of water";
Buildings.Fluid.Sensors.TemperatureTwoPort temRet(
redeclare package Medium = MediumW,
m_flow_nominal=mHeaPum_flow_nominal,
T_start=TRadSup_nominal) "Return water temperature";
//------------------------------------------------------------------------------------//
BoundaryConditions.WeatherData.ReaderTMY3 weaDat(filNam="modelica://Buildings/Resources/weatherdata/USA_IL_Chicago-OHare.Intl.AP.725300_TMY3.mos")
"Weather data reader";
BoundaryConditions.WeatherData.Bus weaBus "Weather data bus";
Modelica.Thermal.HeatTransfer.Sources.PrescribedTemperature TOut
"Outside temperature";
//--------------------------------------------------------------------------------------//
ScrollWaterToWater heaPum(
redeclare package Medium1 = MediumW,
redeclare package Medium2 = MediumW,
redeclare package ref = Media.Refrigerants.R410A,
dp1_nominal=2000,
dp2_nominal=2000,
tau1=15,
tau2=15,
show_T=true,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
m1_flow_nominal=mHeaPum_flow_nominal,
m2_flow_nominal=mHeaPum_flow_nominal,
datHeaPum=Data.ScrollWaterToWater.Heating.Daikin_WRA072_24kW_4_30COP_R410A(),
T1_start=TRadSup_nominal) "Heat pump";
Sources.FixedBoundary sou(
redeclare package Medium = MediumW,
use_T=true,
nPorts=1,
T=281.15) "Fluid source on source side";
Sources.FixedBoundary sin(
redeclare package Medium = MediumW,
use_T=true,
nPorts=1,
T=283.15) "Fluid sink on source side";
Buildings.Fluid.Movers.FlowControlled_m_flow pumHeaPumSou(
redeclare package Medium = MediumW,
y_start=1,
m_flow_start=0.85,
m_flow_nominal=mHeaPum_flow_nominal,
nominalValuesDefineDefaultPressureCurve=true,
use_inputFilter=false,
energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState)
"Pump for heat pump source side";
Modelica.Blocks.Logical.Hysteresis hysteresis(uLow=273.15 + 19, uHigh=273.15 +
21) "Hysteresis controller";
Modelica.Blocks.Logical.Not not2 "Negate output of hysteresis";
Modelica.Blocks.Math.BooleanToReal booToReaPum(realTrue=1, y(start=0))
"Pump signal";
Modelica.Blocks.Logical.And and1;
Modelica.Blocks.Logical.And and2;
Modelica.Blocks.Math.BooleanToReal booToReaPum1(realTrue=1, y(start=0))
"Pump signal";
Modelica.Blocks.Logical.Hysteresis tesConHea(uHigh=0.25*mHeaPum_flow_nominal,
uLow=0.20*mHeaPum_flow_nominal) "Test for flow rate of condenser pump";
Modelica.Blocks.Logical.Hysteresis tesEvaPum(uLow=0.20*mHeaPum_flow_nominal,
uHigh=0.25*mHeaPum_flow_nominal) "Test for flow rate of evaporator pump";
equation
connect(theCon.port_b, vol.heatPort);
connect(preHea.port, vol.heatPort);
connect(heaCap.port, vol.heatPort);
connect(timTab.y[1], preHea.Q_flow);
connect(temSup.port_b, rad.port_a);
connect(temRoo.port, vol.heatPort);
connect(rad.heatPortCon, vol.heatPort);
connect(rad.heatPortRad, vol.heatPort);
connect(weaDat.weaBus, weaBus);
connect(weaBus.TDryBul, TOut.T);
connect(TOut.port, theCon.port_a);
connect(temRet.port_b, heaPum.port_a1);
connect(preSou.ports[1], temRet.port_b);
connect(sou.ports[1], pumHeaPumSou.port_a);
connect(pumHeaPumSou.port_b, heaPum.port_a2);
connect(sin.ports[1], heaPum.port_b2);
connect(hysteresis.y, not2.u);
connect(temRoo.T, hysteresis.u);
connect(pumHeaPum.port_b, temSup.port_a);
connect(temRet.port_a, rad.port_b);
connect(and2.u2, and1.y);
connect(booToReaPum1.y, heaPum.y);
connect(booToReaPum1.u, and2.y);
connect(tesConHea.y, and1.u1);
connect(tesEvaPum.y, and1.u2);
connect(pumHeaPum.m_flow_actual, tesConHea.u);
connect(pumHeaPumSou.m_flow_actual, tesEvaPum.u);
connect(heaPum.port_b1, pumHeaPum.port_a);
connect(booToReaPum.y, pumHeaPum.m_flow_in);
connect(booToReaPum.y, pumHeaPumSou.m_flow_in);
connect(not2.y, and2.u1);
connect(not2.y, booToReaPum.u);
end ScrollWaterToWater_OneRoomRadiator;