Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems
Package of models for subsystems of fifth generation DHC ETS
Information
This package contains models for subsystems composing energy transfer stations in fifth generation district heating and cooling systems.
Extends from Modelica.Icons.VariantsPackage (Icon for package containing variants).
Package Content
| Name | Description |
|---|---|
 Borefield
|
Base subsystem with geothermal borefield |
 Chiller
|
Base subsystem with heat recovery chiller |
 HeatExchanger
|
Base subsystem with district heat exchanger |
 HeatPump
|
Base subsystem with water-to-water heat pump |
 SwitchBox
|
Model for mass flow rate redirection with three-port two-position directional valves |
 Validation
|
Collection of validation models |
Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems.Borefield
Base subsystem with geothermal borefield
Information
This is a model for a borefield system with a variable speed pump and a mixing valve modulated to maintain a maximum inlet temperature.
The system is controlled based on the logic described in Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.Borefield. The pump flow rate is considered proportional to the pump speed under the assumption of a constant flow resistance in the borefield loop. This assumption is justified by the connection of the loop to the buffer tanks, and by the additional assumption that the bypass branch of the mixing valve is balanced with the direct branch.
(The parameter from_dp of the valve model is set to false to
simplify the system of algebraic equations, which, in this specific case,
alleviates non-convergence issues.)
Extends from Buildings.Fluid.Interfaces.PartialTwoPortInterface (Partial model transporting fluid between two ports without storing mass or energy).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | PartialMedium | Medium in the component | |
| replaceable model BoreFieldType | Fluid.Geothermal.Borefields.... | Wall heat transfer | |
| Example | datBorFie | datBorFie(conDat=Buildings.F... | Borefield parameters |
| Generic | perPum | perPum(motorCooledByFluid=fa... | Record with performance data for borefield pump |
| Pressure | dpValBorFie_nominal | dp_nominal/2 | Nominal pressure drop of control valve [Pa] |
| Temperature | TBorWatEntMax | Maximum value of borefield water entering temperature [K] | |
| Real | spePumBorMin | 0.1 | Borefield pump minimum speed |
| Nominal condition | |||
| MassFlowRate | m_flow_nominal | datBorFie.conDat.mBorFie_flo... | Nominal mass flow rate [kg/s] |
| Pressure | dp_nominal | Pressure losses for the entire borefield (control valve excluded) [Pa] | |
| Assumptions | |||
| Boolean | allowFlowReversal | true | = false to simplify equations, assuming, but not enforcing, no flow reversal |
| Advanced | |||
| MassFlowRate | m_flow_small | 1E-4*abs(m_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| Diagnostics | |||
| Boolean | show_T | false | = true, if actual temperature at port is computed |
Connectors
| Type | Name | Description |
|---|---|---|
| FluidPort_a | port_a | Fluid connector a (positive design flow direction is from port_a to port_b) |
| FluidPort_b | port_b | Fluid connector b (positive design flow direction is from port_a to port_b) |
| replaceable model BoreFieldType | Wall heat transfer | |
| input RealInput | yValIso_actual[2] | Isolation valves return position (fractional) |
| input RealInput | u | Control signal from supervisory |
| output RealOutput | PPum | Pump power [W] |
Modelica definition
model Borefield
"Base subsystem with geothermal borefield"
extends Buildings.Fluid.Interfaces.PartialTwoPortInterface(
final m_flow_nominal=datBorFie.conDat.mBorFie_flow_nominal);
replaceable model BoreFieldType=Fluid.Geothermal.Borefields.OneUTube
constrainedby Fluid.Geothermal.Borefields.BaseClasses.PartialBorefield(
redeclare package Medium=Medium,
allowFlowReversal=allowFlowReversal,
borFieDat=datBorFie,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial)
"Wall heat transfer";
replaceable parameter Fluid.Geothermal.Borefields.Data.Borefield.Example datBorFie(
conDat=Buildings.Fluid.Geothermal.Borefields.Data.Configuration.Example())
constrainedby Fluid.Geothermal.Borefields.Data.Borefield.Template(
conDat(
dp_nominal=0))
"Borefield parameters";
replaceable parameter Buildings.Fluid.Movers.Data.Generic perPum(
motorCooledByFluid=false)
constrainedby Buildings.Fluid.Movers.Data.Generic
"Record with performance data for borefield pump";
parameter Modelica.SIunits.Pressure dp_nominal(
displayUnit="Pa")
"Pressure losses for the entire borefield (control valve excluded)";
parameter Modelica.SIunits.Pressure dpValBorFie_nominal=dp_nominal/2
"Nominal pressure drop of control valve";
parameter Modelica.SIunits.Temperature TBorWatEntMax(
displayUnit="degC")
"Maximum value of borefield water entering temperature";
parameter Real spePumBorMin=0.1
"Borefield pump minimum speed";
// IO VARIABLES
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValIso_actual[2]
"Isolation valves return position (fractional)";
Buildings.Controls.OBC.CDL.Interfaces.RealInput u
"Control signal from supervisory";
// COMPONENTS
Fluid.Actuators.Valves.ThreeWayEqualPercentageLinear val(
redeclare final package Medium=Medium,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
from_dp=false,
use_inputFilter=false,
final m_flow_nominal=m_flow_nominal,
final dpValve_nominal=dpValBorFie_nominal,
final dpFixed_nominal=fill(
dp_nominal,
2))
"Mixing valve controlling entering temperature";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pum(
redeclare final package Medium=Medium,
final per=perPum,
final m_flow_nominal=m_flow_nominal,
final dp_nominal=dpValBorFie_nominal+dp_nominal)
"Pump with prescribed mass flow rate";
Fluid.Sensors.TemperatureTwoPort senTEnt(
final tau=
if allowFlowReversal then
1
else
0,
redeclare final package Medium=Medium,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=m_flow_nominal)
"Entering temperature";
BoreFieldType borFie
"Geothermal borefield";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Junction spl(
redeclare final package Medium=Medium,
final m_flow_nominal=m_flow_nominal .* {1,-1,-1})
"Flow splitter";
Fluid.Sensors.TemperatureTwoPort senTLvg(
redeclare final package Medium=Medium,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=m_flow_nominal)
"Leaving temperature";
Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.Borefield
con(
final TBorWatEntMax=TBorWatEntMax,
final spePumBorMin=spePumBorMin)
"Controller";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PPum(
final unit="W")
"Pump power";
Buildings.Controls.OBC.CDL.Continuous.Gain gai1(
final k=m_flow_nominal)
"Scale to nominal mass flow rate";
initial equation
assert(
abs(
datBorFie.conDat.dp_nominal) < Modelica.Constants.eps,
"In "+getInstanceName()+": dp_nominal in the parameter record should be set to zero as the nominal
pressure drop is lumped in the valve model. Use the exposed parameter
dp_nominal instead.",
level=AssertionLevel.warning);
equation
connect(pum.port_b,senTEnt.port_a);
connect(senTEnt.port_b,borFie.port_a);
connect(port_b,spl.port_2);
connect(spl.port_1,senTLvg.port_b);
connect(borFie.port_b,senTLvg.port_a);
connect(spl.port_3,val.port_3);
connect(u,con.u);
connect(yValIso_actual,con.yValIso_actual);
connect(con.yValMix,val.y);
connect(port_a,val.port_1);
connect(val.port_2,pum.port_a);
connect(senTEnt.T,con.TBorWatEnt);
connect(pum.P,PPum);
connect(con.yPum,gai1.u);
connect(gai1.y,pum.m_flow_in);
end Borefield;
Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems.Chiller
Base subsystem with heat recovery chiller
Information
This is a model for a chiller system with constant speed evaporator and condenser pumps, and mixing valves modulated to maintain a minimum condenser inlet temperature (resp. maximum evaporator inlet temperature).
The system is controlled based on the logic described in Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.Chiller. The pump flow rate is considered proportional to the pump speed under the assumption of a constant flow resistance for both the condenser and the evaporator loops. This assumption is justified by the connection of the loops to the buffer tanks, and the additional assumption that the bypass branch of the mixing valves is balanced with the direct branch.
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | Modelica.Media.Interfaces.Pa... | Medium model | |
| Generic | dat | redeclare parameter Building... | Chiller performance data |
| Generic | perPumCon | perPumCon(motorCooledByFluid... | Record with performance data for condenser pump |
| Generic | perPumEva | perPumEva(motorCooledByFluid... | Record with performance data for evaporator pump |
| Nominal condition | |||
| PressureDifference | dpCon_nominal | Nominal pressure drop accross condenser [Pa] | |
| PressureDifference | dpEva_nominal | Nominal pressure drop accross evaporator [Pa] | |
| Pressure | dpValCon_nominal | dpCon_nominal/2 | Nominal pressure drop accross control valve on condenser side [Pa] |
| Pressure | dpValEva_nominal | dpEva_nominal/2 | Nominal pressure drop accross control valve on evaporator side [Pa] |
| Controls | |||
| Temperature | TConWatEntMin | dat.TConEntMin | Minimum value of condenser water entering temperature [K] |
| Temperature | TEvaWatEntMax | dat.TEvaLvgMax - dat.QEva_fl... | Maximum value of evaporator water entering temperature [K] |
| Assumptions | |||
| Boolean | allowFlowReversal | false | = true to allow flow reversal, false restricts to design direction (port_a -> port_b) |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package Medium | Medium model | |
| input BooleanInput | uHea | Heating enable signal |
| input BooleanInput | uCoo | Cooling enable signal |
| input RealInput | TChiWatSupSet | Chilled water supply temperature set point (may be reset down) [K] |
| FluidPort_a | port_aChiWat | Fluid port for chilled water return |
| FluidPort_b | port_bChiWat | Fluid port for chilled water supply |
| FluidPort_a | port_aHeaWat | Fluid port for heating water return |
| FluidPort_b | port_bHeaWat | Fluid port for heating water supply |
| output RealOutput | PChi | Chiller power [W] |
| output RealOutput | PPum | Pump power [W] |
Modelica definition
model Chiller
"Base subsystem with heat recovery chiller"
replaceable package Medium=Modelica.Media.Interfaces.PartialMedium
"Medium model";
parameter Boolean allowFlowReversal=false
"= true to allow flow reversal, false restricts to design direction (port_a -> port_b)";
replaceable parameter Buildings.Fluid.Chillers.Data.ElectricEIR.Generic dat
"Chiller performance data";
replaceable parameter Buildings.Fluid.Movers.Data.Generic perPumCon(
motorCooledByFluid=false)
constrainedby Buildings.Fluid.Movers.Data.Generic
"Record with performance data for condenser pump";
replaceable parameter Buildings.Fluid.Movers.Data.Generic perPumEva(
motorCooledByFluid=false)
constrainedby Buildings.Fluid.Movers.Data.Generic
"Record with performance data for evaporator pump";
parameter Modelica.SIunits.PressureDifference dpCon_nominal(
displayUnit="Pa")
"Nominal pressure drop accross condenser";
parameter Modelica.SIunits.PressureDifference dpEva_nominal(
displayUnit="Pa")
"Nominal pressure drop accross evaporator";
parameter Modelica.SIunits.Pressure dpValCon_nominal=dpCon_nominal/2
"Nominal pressure drop accross control valve on condenser side";
parameter Modelica.SIunits.Pressure dpValEva_nominal=dpEva_nominal/2
"Nominal pressure drop accross control valve on evaporator side";
parameter Modelica.SIunits.Temperature TConWatEntMin(
displayUnit="degC")=dat.TConEntMin
"Minimum value of condenser water entering temperature";
parameter Modelica.SIunits.Temperature TEvaWatEntMax(
displayUnit="degC")=dat.TEvaLvgMax-dat.QEva_flow_nominal/cp_default/dat.mEva_flow_nominal
"Maximum value of evaporator water entering temperature";
// IO CONNECTORS
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput uHea
"Heating enable signal";
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput uCoo
"Cooling enable signal";
Buildings.Controls.OBC.CDL.Interfaces.RealInput TChiWatSupSet(
final unit="K",
displayUnit="degC")
"Chilled water supply temperature set point (may be reset down)";
Modelica.Fluid.Interfaces.FluidPort_a port_aChiWat(
redeclare final package Medium=Medium,
m_flow(
min=
if allowFlowReversal then
-Modelica.Constants.inf
else
0),
h_outflow(
start=Medium.h_default,
nominal=Medium.h_default))
"Fluid port for chilled water return";
Modelica.Fluid.Interfaces.FluidPort_b port_bChiWat(
redeclare final package Medium=Medium,
m_flow(
max=
if allowFlowReversal then
+Modelica.Constants.inf
else
0),
h_outflow(
start=Medium.h_default,
nominal=Medium.h_default))
"Fluid port for chilled water supply";
Modelica.Fluid.Interfaces.FluidPort_a port_aHeaWat(
redeclare final package Medium=Medium,
m_flow(
min=
if allowFlowReversal then
-Modelica.Constants.inf
else
0),
h_outflow(
start=Medium.h_default,
nominal=Medium.h_default))
"Fluid port for heating water return";
Modelica.Fluid.Interfaces.FluidPort_b port_bHeaWat(
redeclare final package Medium=Medium,
m_flow(
max=
if allowFlowReversal then
+Modelica.Constants.inf
else
0),
h_outflow(
start=Medium.h_default,
nominal=Medium.h_default))
"Fluid port for heating water supply";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PChi(
final unit="W")
"Chiller power";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PPum(
final unit="W")
"Pump power";
// COMPONENTS
Fluid.Chillers.ElectricEIR chi(
redeclare final package Medium1=Medium,
redeclare final package Medium2=Medium,
final allowFlowReversal1=allowFlowReversal,
final allowFlowReversal2=allowFlowReversal,
final dp1_nominal=0,
final dp2_nominal=0,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
final per=dat)
"Water cooled chiller (ports indexed 1 are on condenser side)";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pumCon(
redeclare final package Medium=Medium,
final per=perPumCon,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=dat.mCon_flow_nominal,
final dp_nominal=dpCon_nominal+dpValCon_nominal)
"Condenser pump";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pumEva(
redeclare final package Medium=Medium,
final per=perPumEva,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=dat.mEva_flow_nominal,
final dp_nominal=dpEva_nominal+dpValEva_nominal)
"Evaporator pump";
Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.Chiller
con(
final TConWatEntMin=TConWatEntMin,
final TEvaWatEntMax=TEvaWatEntMax)
"Controller";
Buildings.Fluid.Sensors.TemperatureTwoPort senTConLvg(
redeclare final package Medium=Medium,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=dat.mCon_flow_nominal)
"Condenser water leaving temperature";
Buildings.Fluid.Sensors.TemperatureTwoPort senTConEnt(
redeclare final package Medium=Medium,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=dat.mCon_flow_nominal)
"Condenser water entering temperature";
Buildings.Fluid.Sensors.TemperatureTwoPort senTEvaEnt(
redeclare final package Medium=Medium,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=dat.mEva_flow_nominal)
"Evaporator water entering temperature";
Buildings.Fluid.Sensors.TemperatureTwoPort senTEvaLvg(
redeclare final package Medium=Medium,
final allowFlowReversal=allowFlowReversal,
final m_flow_nominal=dat.mEva_flow_nominal)
"Evaporator water leaving temperature";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Junction splEva(
redeclare final package Medium=Medium,
final m_flow_nominal=dat.mEva_flow_nominal .* {1,-1,-1})
"Flow splitter for the evaporator water circuit";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Junction splConMix(
redeclare final package Medium=Medium,
final m_flow_nominal=dat.mCon_flow_nominal .* {1,-1,-1})
"Flow splitter";
Fluid.Actuators.Valves.ThreeWayEqualPercentageLinear valEva(
redeclare final package Medium=Medium,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
from_dp=false,
use_inputFilter=false,
final m_flow_nominal=dat.mEva_flow_nominal,
final dpValve_nominal=dpValEva_nominal,
final dpFixed_nominal=fill(
dpEva_nominal,
2))
"Control valve for maximum evaporator water entering temperature";
Fluid.Actuators.Valves.ThreeWayEqualPercentageLinear valCon(
redeclare final package Medium=Medium,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
from_dp=false,
use_inputFilter=false,
final m_flow_nominal=dat.mCon_flow_nominal,
final dpValve_nominal=dpValCon_nominal,
final dpFixed_nominal=fill(
dpCon_nominal,
2))
"Control valve for minimum condenser water entering temperature";
Buildings.Controls.OBC.CDL.Conversions.BooleanToReal booToRea
"Constant speed primary pumps control signal";
Buildings.Controls.OBC.CDL.Continuous.Add add2;
Buildings.Controls.OBC.CDL.Continuous.Gain gai1(
final k=dat.mCon_flow_nominal)
"Scale to nominal mass flow rate";
Buildings.Controls.OBC.CDL.Continuous.Gain gai2(
final k=dat.mEva_flow_nominal)
"Scale to nominal mass flow rate";
protected
final parameter Medium.ThermodynamicState sta_default=Medium.setState_pTX(
T=Medium.T_default,
p=Medium.p_default,
X=Medium.X_default[1:Medium.nXi])
"Medium state at default properties";
final parameter Modelica.SIunits.SpecificHeatCapacity cp_default=Medium.specificHeatCapacityCp(
sta_default)
"Specific heat capacity of the fluid";
equation
connect(splConMix.port_3,valCon.port_3);
connect(valCon.port_2,pumCon.port_a);
connect(pumEva.port_b,splEva.port_1);
connect(splEva.port_3,valEva.port_3);
connect(con.yValEva,valEva.y);
connect(con.yValCon,valCon.y);
connect(con.yChi,chi.on);
connect(uHea,con.uHea);
connect(uCoo,con.uCoo);
connect(senTConEnt.T,con.TConWatEnt);
connect(senTEvaEnt.T,con.TEvaWatEnt);
connect(splConMix.port_2,port_bHeaWat);
connect(splEva.port_2,port_bChiWat);
connect(port_aHeaWat,valCon.port_1);
connect(port_aChiWat,valEva.port_1);
connect(valEva.port_2,senTEvaEnt.port_a);
connect(senTEvaLvg.port_b,pumEva.port_a);
connect(senTEvaLvg.port_a,chi.port_b2);
connect(senTEvaEnt.port_b,chi.port_a2);
connect(chi.port_b1,senTConLvg.port_a);
connect(senTConLvg.port_b,splConMix.port_1);
connect(pumCon.port_b,senTConEnt.port_a);
connect(senTConEnt.port_b,chi.port_a1);
connect(chi.P,PChi);
connect(add2.y,PPum);
connect(pumEva.P,add2.u2);
connect(pumCon.P,add2.u1);
connect(con.yChi,booToRea.u);
connect(booToRea.y,gai2.u);
connect(gai2.y,pumEva.m_flow_in);
connect(gai1.y,pumCon.m_flow_in);
connect(booToRea.y,gai1.u);
connect(TChiWatSupSet,chi.TSet);
end Chiller;
Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems.HeatExchanger
Base subsystem with district heat exchanger
Information
This is a model for a district heat exchanger system with a variable speed pump on the secondary side, and a variable speed pump (in case of a passive network) or a two-way modulating valve (in case of an active network) on the primary side.
The system is controlled based on the logic described in Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.HeatExchanger. The pump flow rate is considered proportional to the pump speed under the assumption of a constant flow resistance in both the primary and the secondary loops.
Extends from Fluid.Interfaces.PartialFourPortInterface (Partial model transporting fluid between two ports without storing mass or energy).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium1 | PartialMedium | Medium 1 in the component | |
| replaceable package Medium2 | PartialMedium | Medium 2 in the component | |
| ConnectionConfiguration | conCon | District connection configuration | |
| Generic | perPum1 | perPum1(motorCooledByFluid=f... | Record with performance data for primary pump |
| Generic | perPum2 | perPum2(motorCooledByFluid=f... | Record with performance data for secondary pump |
| Nominal condition | |||
| MassFlowRate | m1_flow_nominal | abs(QHex_flow_nominal/4200/(... | Nominal mass flow rate [kg/s] |
| MassFlowRate | m2_flow_nominal | abs(QHex_flow_nominal/4200/(... | Nominal mass flow rate [kg/s] |
| PressureDifference | dp1Hex_nominal | Nominal pressure drop across heat exchanger on district side [Pa] | |
| PressureDifference | dp2Hex_nominal | Nominal pressure drop across heat exchanger on building side [Pa] | |
| PressureDifference | dpVal1Hex_nominal | dp1Hex_nominal/2 | Nominal pressure drop of primary control valve [Pa] |
| PressureDifference | dpVal2Hex_nominal | dp2Hex_nominal/2 | Nominal pressure drop of secondary control valve [Pa] |
| HeatFlowRate | QHex_flow_nominal | Nominal heat flow rate (from district to building) [W] | |
| Temperature | T_a1Hex_nominal | Nominal water inlet temperature on district side [K] | |
| Temperature | T_b1Hex_nominal | Nominal water outlet temperature on district side [K] | |
| Temperature | T_a2Hex_nominal | Nominal water inlet temperature on building side [K] | |
| Temperature | T_b2Hex_nominal | Nominal water outlet temperature on building side [K] | |
| Controls | |||
| Real | spePum1HexMin | 0.1 | Heat exchanger primary pump minimum speed (fractional) [1] |
| Real | yVal1HexMin | 0.1 | Minimum valve opening for temperature measurement (fractional) [1] |
| Real | spePum2HexMin | 0.1 | Heat exchanger secondary pump minimum speed (fractional) [1] |
| TemperatureDifference | dT1HexSet[2] | Primary side deltaT set point schedule (index 1 for heat rejection) [K] | |
| Real | k[2] | {0.05,0.1} | Gain schedule for controller (index 1 for heat rejection) |
| Time | Ti | 120 | Time constant of integrator block [s] |
| Assumptions | |||
| Boolean | allowFlowReversal1 | true | = false to simplify equations, assuming, but not enforcing, no flow reversal for medium 1 |
| Boolean | allowFlowReversal2 | true | = false to simplify equations, assuming, but not enforcing, no flow reversal for medium 2 |
| Advanced | |||
| MassFlowRate | m1_flow_small | 1E-4*abs(m1_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| MassFlowRate | m2_flow_small | 1E-4*abs(m2_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| Diagnostics | |||
| Boolean | show_T | false | = true, if actual temperature at port is computed |
Connectors
| Type | Name | Description |
|---|---|---|
| FluidPort_a | port_a1 | Fluid connector a1 (positive design flow direction is from port_a1 to port_b1) |
| FluidPort_b | port_b1 | Fluid connector b1 (positive design flow direction is from port_a1 to port_b1) |
| FluidPort_a | port_a2 | Fluid connector a2 (positive design flow direction is from port_a2 to port_b2) |
| FluidPort_b | port_b2 | Fluid connector b2 (positive design flow direction is from port_a2 to port_b2) |
| input RealInput | yValIso_actual[2] | Isolation valves return position (index 1 for condenser) |
| input RealInput | u | Control signal for secondary side (from supervisory) |
| output RealOutput | PPum | Power drawn by pump motors [W] |
Modelica definition
model HeatExchanger
"Base subsystem with district heat exchanger"
extends Fluid.Interfaces.PartialFourPortInterface(
final m1_flow_nominal=abs(
QHex_flow_nominal/4200/(T_b1Hex_nominal-T_a1Hex_nominal)),
final m2_flow_nominal=abs(
QHex_flow_nominal/4200/(T_b2Hex_nominal-T_a2Hex_nominal)));
parameter Buildings.Experimental.DHC.EnergyTransferStations.Types.ConnectionConfiguration
conCon
"District connection configuration";
replaceable parameter Buildings.Fluid.Movers.Data.Generic perPum1(
motorCooledByFluid=false)
constrainedby Buildings.Fluid.Movers.Data.Generic
"Record with performance data for primary pump";
replaceable parameter Buildings.Fluid.Movers.Data.Generic perPum2(
motorCooledByFluid=false)
constrainedby Buildings.Fluid.Movers.Data.Generic
"Record with performance data for secondary pump";
parameter Modelica.SIunits.PressureDifference dp1Hex_nominal(
displayUnit="Pa")
"Nominal pressure drop across heat exchanger on district side";
parameter Modelica.SIunits.PressureDifference dp2Hex_nominal(
displayUnit="Pa")
"Nominal pressure drop across heat exchanger on building side";
parameter Modelica.SIunits.PressureDifference dpVal1Hex_nominal(
displayUnit="Pa")=dp1Hex_nominal/2
"Nominal pressure drop of primary control valve";
parameter Modelica.SIunits.PressureDifference dpVal2Hex_nominal(
displayUnit="Pa")=dp2Hex_nominal/2
"Nominal pressure drop of secondary control valve";
parameter Modelica.SIunits.HeatFlowRate QHex_flow_nominal
"Nominal heat flow rate (from district to building)";
parameter Modelica.SIunits.Temperature T_a1Hex_nominal
"Nominal water inlet temperature on district side";
parameter Modelica.SIunits.Temperature T_b1Hex_nominal
"Nominal water outlet temperature on district side";
parameter Modelica.SIunits.Temperature T_a2Hex_nominal
"Nominal water inlet temperature on building side";
parameter Modelica.SIunits.Temperature T_b2Hex_nominal
"Nominal water outlet temperature on building side";
parameter Real spePum1HexMin(
final unit="1",
min=0)=0.1
"Heat exchanger primary pump minimum speed (fractional)";
parameter Real yVal1HexMin(
final unit="1",
min=0.01)=0.1
"Minimum valve opening for temperature measurement (fractional)";
parameter Real spePum2HexMin(
final unit="1",
min=0.01)=0.1
"Heat exchanger secondary pump minimum speed (fractional)";
parameter Modelica.SIunits.TemperatureDifference dT1HexSet[2]
"Primary side deltaT set point schedule (index 1 for heat rejection)";
parameter Real k[2]={0.05,0.1}
"Gain schedule for controller (index 1 for heat rejection)";
parameter Modelica.SIunits.Time Ti=120
"Time constant of integrator block";
// IO CONNECTORS
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValIso_actual[2]
"Isolation valves return position (index 1 for condenser)";
Buildings.Controls.OBC.CDL.Interfaces.RealInput u
"Control signal for secondary side (from supervisory)";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PPum(
final unit="W")
"Power drawn by pump motors";
// COMPONENTS
Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.HeatExchanger
con(
final conCon=conCon,
final spePum1HexMin=spePum1HexMin,
final spePum2HexMin=spePum2HexMin,
final dT1HexSet=dT1HexSet,
final k=k,
final Ti=Ti)
"District heat exchanger loop controller";
Fluid.HeatExchangers.PlateHeatExchangerEffectivenessNTU hex(
redeclare final package Medium1=Medium1,
redeclare final package Medium2=Medium2,
final use_Q_flow_nominal=true,
configuration=Buildings.Fluid.Types.HeatExchangerConfiguration.CounterFlow,
final allowFlowReversal1=allowFlowReversal1,
final allowFlowReversal2=allowFlowReversal2,
final dp1_nominal=
if have_val1Hex then
0
else
dp1Hex_nominal,
final dp2_nominal=0,
final m1_flow_nominal=m1_flow_nominal,
final m2_flow_nominal=m2_flow_nominal,
final Q_flow_nominal=QHex_flow_nominal,
final T_a1_nominal=T_a1Hex_nominal,
final T_a2_nominal=T_a2Hex_nominal)
"Heat exchanger";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pum1Hex(
redeclare final package Medium=Medium1,
final per=perPum1,
final m_flow_nominal=m1_flow_nominal,
final dp_nominal=dp1Hex_nominal,
final allowFlowReversal=allowFlowReversal1) if not have_val1Hex
"District heat exchanger primary pump";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pum2Hex(
redeclare final package Medium=Medium2,
final per=perPum2,
final m_flow_nominal=m2_flow_nominal,
final dp_nominal=dp2Hex_nominal+dpVal2Hex_nominal,
final allowFlowReversal=allowFlowReversal2)
"Secondary pump";
Fluid.Sensors.TemperatureTwoPort senT2HexWatEnt(
redeclare final package Medium=Medium2,
final m_flow_nominal=m2_flow_nominal,
final allowFlowReversal=allowFlowReversal2)
"Heat exchanger secondary water entering temperature";
Fluid.Sensors.TemperatureTwoPort senT2HexWatLvg(
redeclare final package Medium=Medium2,
final m_flow_nominal=m2_flow_nominal,
final allowFlowReversal=allowFlowReversal2)
"Heat exchanger secondary water leaving temperature";
Buildings.Controls.OBC.CDL.Continuous.Gain gai2(
final k=m2_flow_nominal)
"Scale to nominal mass flow rate";
Fluid.Actuators.Valves.TwoWayEqualPercentage val1Hex(
redeclare final package Medium=Medium1,
final m_flow_nominal=m1_flow_nominal,
from_dp=true,
final dpValve_nominal=dpVal1Hex_nominal,
use_inputFilter=false,
final dpFixed_nominal=dp1Hex_nominal) if have_val1Hex
"Heat exchanger primary control valve";
Fluid.Sensors.TemperatureTwoPort senT1HexWatEnt(
redeclare final package Medium=Medium1,
final m_flow_nominal=m1_flow_nominal,
final allowFlowReversal=allowFlowReversal1)
"Heat exchanger primary water entering temperature";
Fluid.Sensors.TemperatureTwoPort senT1HexWatLvg(
redeclare final package Medium=Medium1,
final m_flow_nominal=m1_flow_nominal,
final allowFlowReversal=allowFlowReversal1)
"Heat exchanger primary water leaving temperature";
Buildings.Controls.OBC.CDL.Continuous.Gain gai1(
final k=m1_flow_nominal) if not have_val1Hex
"Scale to nominal mass flow rate";
Buildings.Controls.OBC.CDL.Continuous.MultiSum totPPum(
final nin=
if have_val1Hex then
1
else
2)
"Total pump power";
Fluid.Actuators.Valves.ThreeWayEqualPercentageLinear val2Hex(
redeclare final package Medium=Medium2,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial,
from_dp=false,
use_inputFilter=false,
final m_flow_nominal=m2_flow_nominal,
final dpValve_nominal=dpVal2Hex_nominal,
final dpFixed_nominal=fill(
dp2Hex_nominal,
2))
"Control valve";
Buildings.Experimental.DHC.EnergyTransferStations.BaseClasses.Junction spl(
redeclare final package Medium=Medium2,
final m_flow_nominal=m2_flow_nominal .* {1,-1,-1})
"Flow splitter";
protected
parameter Boolean have_val1Hex=conCon == Buildings.Experimental.DHC.EnergyTransferStations.Types.ConnectionConfiguration.TwoWayValve
"True in case of control valve on district side, false in case of a pump";
equation
if not have_val1Hex then
connect(senT1HexWatLvg.port_b,port_b1);
else
connect(port_a1,senT1HexWatEnt.port_a);
end if;
connect(gai2.y,pum2Hex.m_flow_in);
connect(port_a1,pum1Hex.port_a);
connect(val1Hex.port_b,port_b1);
connect(hex.port_b1,senT1HexWatLvg.port_a);
connect(pum2Hex.port_b,senT2HexWatEnt.port_a);
connect(senT2HexWatEnt.port_b,hex.port_a2);
connect(senT1HexWatEnt.port_b,hex.port_a1);
connect(pum1Hex.port_b,senT1HexWatEnt.port_a);
connect(val1Hex.port_a,senT1HexWatLvg.port_b);
connect(con.y1Hex,val1Hex.y);
connect(con.y1Hex,gai1.u);
connect(gai1.y,pum1Hex.m_flow_in);
connect(pum1Hex.P,totPPum.u[2]);
connect(pum2Hex.P,totPPum.u[1]);
connect(totPPum.y,PPum);
connect(yValIso_actual,con.yValIso);
connect(con.yPum2Hex,gai2.u);
connect(u,con.u);
connect(hex.port_b2,senT2HexWatLvg.port_a);
connect(val2Hex.port_2,pum2Hex.port_a);
connect(port_a2,val2Hex.port_1);
connect(spl.port_1,senT2HexWatLvg.port_b);
connect(spl.port_2,port_b2);
connect(spl.port_3,val2Hex.port_3);
connect(con.yVal2Hex,val2Hex.y);
connect(senT1HexWatLvg.T,con.T1HexWatLvg);
connect(senT1HexWatEnt.T,con.T1HexWatEnt);
end HeatExchanger;
Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems.HeatPump
Base subsystem with water-to-water heat pump
Information
This model represents a water-to-water heat pump, an evaporator water pump,
and an optional condenser water pump if have_pumCon is set to
true.
The heat pump model is described in
Buildings.Fluid.HeatPumps.Carnot_TCon.
By default variable speed pumps are considered.
Constant speed pumps may also be represented by setting have_varFloEva
and have_varFloCon to false.
Controls
The system is enabled when the input control signal uEna switches to
true.
When enabled,
- the evaporator and optionally the condenser water pumps are commanded on and supply either the mass flow rate set point provided as an input in the case of variable speed pumps, or the nominal mass flow rate in the case of constant speed pumps,
- the heat pump is commanded on when the evaporator and optionally the condenser water pump are proven on. When enabled, the heat pump controller—idealized in this model—tracks the supply temperature set point at the condenser outlet.
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium1 | Modelica.Media.Interfaces.Pa... | Medium model on condenser side | |
| replaceable package Medium2 | Modelica.Media.Interfaces.Pa... | Medium model on evaporator side | |
| Boolean | have_pumCon | true | Set to true to include a condenser pump (false for external pump) |
| Boolean | have_varFloCon | true | Set to true for a variable condenser flow |
| Boolean | have_varFloEva | true | Set to true for a variable evaporator flow |
| Nominal condition | |||
| Real | COP_nominal | Heat pump COP [1] | |
| Temperature | TCon_nominal | Condenser outlet temperature used to compute COP_nominal [K] | |
| Temperature | TEva_nominal | Evaporator outlet temperature used to compute COP_nominal [K] | |
| HeatFlowRate | Q1_flow_nominal | Heating heat flow rate [W] | |
| TemperatureDifference | dT1_nominal | 5 | Temperature difference condenser outlet-inlet [K] |
| TemperatureDifference | dT2_nominal | -5 | Temperature difference evaporator outlet-inlet [K] |
| Pressure | dp1_nominal | Pressure difference over condenser [Pa] | |
| Pressure | dp2_nominal | Pressure difference over evaporator [Pa] | |
| Assumptions | |||
| Boolean | allowFlowReversal1 | false | Set to true to allow flow reversal on condenser side |
| Boolean | allowFlowReversal2 | false | Set to true to allow flow reversal on evaporator side |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package Medium1 | Medium model on condenser side | |
| replaceable package Medium2 | Medium model on evaporator side | |
| input BooleanInput | uEna | Enable signal |
| input RealInput | TSupSet | Supply temperature set point [K] |
| input RealInput | m1_flow | Condenser mass flow rate [kg/s] |
| input RealInput | m2_flow | Evaporator mass flow rate [kg/s] |
| FluidPort_a | port_a2 | Fluid port for entering evaporator water |
| FluidPort_b | port_b2 | Fluid port for leaving evaporator water |
| FluidPort_a | port_a1 | Fluid port for entering condenser water |
| FluidPort_b | port_b1 | Fluid port for leaving condenser water |
| output RealOutput | PHea | Heat pump power [W] |
| output RealOutput | PPum | Pump power [W] |
| output RealOutput | mEva_flow | Evaporator water mass flow rate [kg/s] |
Modelica definition
model HeatPump "Base subsystem with water-to-water heat pump"
replaceable package Medium1=Modelica.Media.Interfaces.PartialMedium
"Medium model on condenser side";
replaceable package Medium2=Modelica.Media.Interfaces.PartialMedium
"Medium model on evaporator side";
parameter Boolean have_pumCon = true
"Set to true to include a condenser pump (false for external pump)";
parameter Boolean have_varFloCon = true
"Set to true for a variable condenser flow";
parameter Boolean have_varFloEva = true
"Set to true for a variable evaporator flow";
parameter Real COP_nominal(final unit="1")
"Heat pump COP";
parameter Modelica.SIunits.Temperature TCon_nominal
"Condenser outlet temperature used to compute COP_nominal";
parameter Modelica.SIunits.Temperature TEva_nominal
"Evaporator outlet temperature used to compute COP_nominal";
parameter Modelica.SIunits.HeatFlowRate Q1_flow_nominal(min=0)
"Heating heat flow rate";
parameter Modelica.SIunits.TemperatureDifference dT1_nominal(
final min=0) = 5 "Temperature difference condenser outlet-inlet";
parameter Modelica.SIunits.TemperatureDifference dT2_nominal(
final max=0) = -5 "Temperature difference evaporator outlet-inlet";
parameter Modelica.SIunits.Pressure dp1_nominal(displayUnit="Pa")
"Pressure difference over condenser";
parameter Modelica.SIunits.Pressure dp2_nominal(displayUnit="Pa")
"Pressure difference over evaporator";
parameter Boolean allowFlowReversal1=false
"Set to true to allow flow reversal on condenser side";
parameter Boolean allowFlowReversal2=false
"Set to true to allow flow reversal on evaporator side";
final parameter Modelica.SIunits.MassFlowRate m1_flow_nominal(min=0)=
heaPum.m1_flow_nominal
"Mass flow rate on condenser side";
final parameter Modelica.SIunits.MassFlowRate m2_flow_nominal(min=0)=
heaPum.m2_flow_nominal
"Mass flow rate on evaporator side";
// IO CONNECTORS
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput uEna(start=false)
"Enable signal";
Buildings.Controls.OBC.CDL.Interfaces.RealInput TSupSet(
final unit="K",
displayUnit="degC")
"Supply temperature set point";
Buildings.Controls.OBC.CDL.Interfaces.RealInput m1_flow(
final unit="kg/s") if have_varFloCon and have_pumCon
"Condenser mass flow rate";
Buildings.Controls.OBC.CDL.Interfaces.RealInput m2_flow(
final unit="kg/s") if have_varFloEva
"Evaporator mass flow rate";
Modelica.Fluid.Interfaces.FluidPort_a port_a2(
redeclare final package Medium = Medium2,
m_flow(min=if allowFlowReversal2 then -Modelica.Constants.inf else 0),
h_outflow(start=Medium2.h_default, nominal=Medium2.h_default))
"Fluid port for entering evaporator water";
Modelica.Fluid.Interfaces.FluidPort_b port_b2(
redeclare final package Medium = Medium2,
m_flow(max=if allowFlowReversal2 then +Modelica.Constants.inf else 0),
h_outflow(start=Medium2.h_default, nominal=Medium2.h_default))
"Fluid port for leaving evaporator water";
Modelica.Fluid.Interfaces.FluidPort_a port_a1(
redeclare final package Medium = Medium1,
m_flow(min=if allowFlowReversal1 then -Modelica.Constants.inf else 0),
h_outflow(start=Medium1.h_default, nominal=Medium1.h_default))
"Fluid port for entering condenser water";
Modelica.Fluid.Interfaces.FluidPort_b port_b1(
redeclare final package Medium = Medium1,
m_flow(max=if allowFlowReversal1 then +Modelica.Constants.inf else 0),
h_outflow(start=Medium1.h_default, nominal=Medium1.h_default))
"Fluid port for leaving condenser water";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PHea(
final unit="W") "Heat pump power";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PPum(
final unit="W") "Pump power";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput mEva_flow(final unit="kg/s")
"Evaporator water mass flow rate";
// COMPONENTS
Fluid.HeatPumps.Carnot_TCon heaPum(
redeclare final package Medium1 = Medium1,
redeclare final package Medium2 = Medium2,
final dTEva_nominal=dT2_nominal,
final dTCon_nominal=dT1_nominal,
final TCon_nominal=TCon_nominal,
final TEva_nominal=TEva_nominal,
final allowFlowReversal1=allowFlowReversal1,
final allowFlowReversal2=allowFlowReversal2,
final use_eta_Carnot_nominal=false,
final COP_nominal=COP_nominal,
final QCon_flow_nominal=Q1_flow_nominal,
final dp1_nominal=dp1_nominal,
final dp2_nominal=dp2_nominal)
"Heat pump (index 1 for condenser side)";
DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pumEva(
redeclare final package Medium = Medium2,
final m_flow_nominal=m2_flow_nominal,
final allowFlowReversal=allowFlowReversal2)
"Heat pump evaporator water pump";
DHC.EnergyTransferStations.BaseClasses.Pump_m_flow pumCon(
redeclare final package Medium = Medium1,
final m_flow_nominal=m1_flow_nominal,
final allowFlowReversal=allowFlowReversal1) if have_pumCon
"Heat pump condenser water pump";
Buildings.Controls.OBC.CDL.Continuous.Sources.Constant floConNom(
final k=m1_flow_nominal) if not have_varFloCon
"Nominal flow rate";
Buildings.Controls.OBC.CDL.Conversions.BooleanToReal booToRea;
Buildings.Controls.OBC.CDL.Continuous.Sources.Constant floEvaNom(
final k=m2_flow_nominal) if not have_varFloEva
"Nominal flow rate";
Fluid.Sensors.TemperatureTwoPort senTConLvg(
redeclare final package Medium = Medium1,
final allowFlowReversal=allowFlowReversal1,
final m_flow_nominal=m1_flow_nominal)
"Condenser water leaving temperature";
Fluid.Sensors.TemperatureTwoPort senTConEnt(
redeclare final package Medium = Medium1,
final allowFlowReversal=allowFlowReversal1,
final m_flow_nominal=m1_flow_nominal)
"Condenser water entering temperature";
Buildings.Controls.OBC.CDL.Logical.Switch enaHeaPum(
u2(start=false))
"Enable heat pump by switching to actual set point";
Buildings.Controls.OBC.CDL.Continuous.Add add2 "Adder";
Modelica.Blocks.Sources.Constant zer(final k=0) if not have_pumCon
"Replacement variable";
Buildings.Controls.OBC.CDL.Continuous.GreaterThreshold staPum[2](
y(each start=false),
t=1e-2 .* {m1_flow_nominal,m2_flow_nominal},
h=0.5e-2 .* {m1_flow_nominal, m2_flow_nominal})
"Pump return status";
Buildings.Controls.OBC.CDL.Logical.And ena
"Enable heat pump if pump return status on";
Modelica.Blocks.Sources.Constant one(final k=1) if not have_pumCon
"Replacement variable";
Buildings.Controls.OBC.CDL.Continuous.Product floCon if have_pumCon
"Zero flow rate if not enabled";
Buildings.Controls.OBC.CDL.Continuous.Product floEva
"Zero flow rate if not enabled";
protected
Fluid.FixedResistances.LosslessPipe pip(
redeclare package Medium = Medium1,
m_flow_nominal=m1_flow_nominal) if not have_pumCon
"Dummy connection used if model is configured to have no pump";
equation
connect(pumEva.port_b,heaPum. port_a2);
connect(heaPum.port_b1,senTConLvg. port_a);
connect(senTConEnt.port_b,heaPum. port_a1);
connect(senTConEnt.T,enaHeaPum. u3);
connect(enaHeaPum.y,heaPum. TSet);
connect(uEna, booToRea.u);
connect(TSupSet, enaHeaPum.u1);
connect(heaPum.port_b2, port_b2);
connect(senTConLvg.port_b, port_b1);
connect(pumEva.m_flow_actual, mEva_flow);
connect(port_a2, pumEva.port_a);
connect(port_a1, pumCon.port_a);
connect(add2.y, PPum);
connect(heaPum.P, PHea);
connect(pumCon.P, add2.u2);
connect(pumEva.P, add2.u1);
connect(pumCon.port_b, senTConEnt.port_a);
connect(zer.y, add2.u2);
connect(pumCon.m_flow_actual, staPum[1].u);
connect(pumEva.m_flow_actual, staPum[2].u);
connect(staPum[1].y, ena.u1);
connect(staPum[2].y, ena.u2);
connect(ena.y, enaHeaPum.u2);
connect(one.y, staPum[1].u);
connect(booToRea.y, floCon.u1);
connect(m1_flow, floCon.u2);
connect(booToRea.y, floEva.u1);
connect(floConNom.y, floCon.u2);
connect(m2_flow, floEva.u2);
connect(floEvaNom.y, floEva.u2);
connect(floEva.y, pumEva.m_flow_in);
connect(floCon.y, pumCon.m_flow_in);
connect(port_a1, pip.port_a);
connect(pip.port_b, senTConEnt.port_a);
end HeatPump;
Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems.SwitchBox
Model for mass flow rate redirection with three-port two-position directional valves
Information
This model represents a hydronic arrangement avoid flow reversal in the service line, for instance when connecting an energy transfer station such as the one modeled in Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.HeatPumpHeatExchanger. For that intent, two three-port two-position directional valves are used. The valves are actuated based on the logic described in Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Controls.SwitchBox.
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | Modelica.Media.Interfaces.Pa... | Medium model | |
| Real | trueHoldDuration | 60 | true hold duration [s] |
| Real | falseHoldDuration | trueHoldDuration | false hold duration [s] |
| Nominal condition | |||
| MassFlowRate | m_flow_nominal | Nominal mass flow rate [kg/s] | |
| PressureDifference | dpValve_nominal | 5000 | Valve pressure drop at nominal conditions [Pa] |
| Dynamics | |||
| Equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance (except for the pump always modeled in steady state) |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package Medium | Medium model | |
| FluidPort_b | port_bSup | Supply line outlet port |
| FluidPort_b | port_bRet | Return line outlet port |
| FluidPort_a | port_aSup | Supply line inlet port |
| FluidPort_a | port_aRet | Return line inlet port |
| input RealInput | mRev_flow | Service water mass flow rate in reverse direction [kg/s] |
| input RealInput | mPos_flow | Service water mass flow rate in positive direction [kg/s] |
Modelica definition
model SwitchBox
"Model for mass flow rate redirection with three-port two-position directional valves"
replaceable package Medium = Modelica.Media.Interfaces.PartialMedium
"Medium model";
parameter Modelica.SIunits.MassFlowRate m_flow_nominal
"Nominal mass flow rate";
parameter Modelica.SIunits.PressureDifference dpValve_nominal(
min=0, displayUnit="Pa") = 5000
"Valve pressure drop at nominal conditions";
parameter Real trueHoldDuration(
final unit="s") = 60
"true hold duration";
parameter Real falseHoldDuration(
final unit="s") = trueHoldDuration
"false hold duration";
parameter Modelica.Fluid.Types.Dynamics energyDynamics=
Modelica.Fluid.Types.Dynamics.FixedInitial
"Type of energy balance (except for the pump always modeled in steady state)";
// IO CONECTORS
Modelica.Fluid.Interfaces.FluidPort_b port_bSup(
redeclare package Medium = Medium)
"Supply line outlet port";
Modelica.Fluid.Interfaces.FluidPort_b port_bRet(
redeclare final package Medium = Medium)
"Return line outlet port";
Modelica.Fluid.Interfaces.FluidPort_a port_aSup(
redeclare final package Medium = Medium)
"Supply line inlet port";
Modelica.Fluid.Interfaces.FluidPort_a port_aRet(
redeclare final package Medium = Medium)
"Return line inlet port";
Buildings.Controls.OBC.CDL.Interfaces.RealInput mRev_flow(final unit="kg/s")
"Service water mass flow rate in reverse direction";
Buildings.Controls.OBC.CDL.Interfaces.RealInput mPos_flow(final unit="kg/s")
"Service water mass flow rate in positive direction";
// COMPONENTS
DHC.EnergyTransferStations.BaseClasses.Junction splSup(
redeclare final package Medium = Medium,
m_flow_nominal={1,1,1}*m_flow_nominal)
"Flow splitter";
DHC.EnergyTransferStations.BaseClasses.Junction splRet(
redeclare final package Medium = Medium,
m_flow_nominal={1,1,1}*m_flow_nominal)
"Flow splitter";
DHC.EnergyTransferStations.Combined.Generation5.Controls.SwitchBox con(
final trueHoldDuration=trueHoldDuration,
final falseHoldDuration=falseHoldDuration)
"Switch box controller";
Fluid.Actuators.Valves.ThreeWayEqualPercentageLinear valSup(
redeclare package Medium = Medium,
dpValve_nominal=dpValve_nominal,
use_inputFilter=false,
m_flow_nominal=m_flow_nominal,
linearized={true,true},
energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState)
"Directional valve";
Fluid.Actuators.Valves.ThreeWayEqualPercentageLinear valRet(
redeclare package Medium = Medium,
dpValve_nominal=dpValve_nominal,
use_inputFilter=false,
m_flow_nominal=m_flow_nominal,
linearized={true,true},
energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState)
"Directional valve";
equation
connect(port_bSup, splSup.port_2);
connect(mRev_flow, con.mRev_flow);
connect(splRet.port_1, port_aRet);
connect(valSup.port_1, splSup.port_1);
connect(valSup.port_3, splRet.port_3);
connect(splRet.port_2, valRet.port_1);
connect(splSup.port_3, valRet.port_3);
connect(valRet.port_2, port_bRet);
connect(valSup.port_2, port_aSup);
connect(con.y, valRet.y);
connect(mPos_flow, con.mPos_flow);
connect(con.y, valSup.y);
end SwitchBox;
Buildings.Experimental.DHC.EnergyTransferStations.Combined.Generation5.Subsystems.Borefield.BoreFieldType
Wall heat transfer
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | PartialMedium | Medium in the component | |
| Time | tLoaAgg | 300 | Time resolution of load aggregation [s] |
| Integer | nCel | 5 | Number of cells per aggregation level |
| Integer | nSeg | 10 | Number of segments to use in vertical discretization of the boreholes |
| Template | borFieDat | Borefield data | |
| Assumptions | |||
| Boolean | allowFlowReversal | true | = false to simplify equations, assuming, but not enforcing, no flow reversal |
| Advanced | |||
| MassFlowRate | m_flow_small | 1E-4*abs(m_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| Boolean | forceGFunCalc | false | Set to true to force the thermal response to be calculated at the start instead of checking whether this has been pre-computed |
| Diagnostics | |||
| Boolean | show_T | false | = true, if actual temperature at port is computed |
| Flow resistance | |||
| Boolean | from_dp | false | = true, use m_flow = f(dp) else dp = f(m_flow) |
| Boolean | linearizeFlowResistance | false | = true, use linear relation between m_flow and dp for any flow rate |
| Real | deltaM | 0.1 | Fraction of nominal flow rate where flow transitions to laminar |
| Dynamics | |||
| Equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Boolean | dynFil | true | Set to false to remove the dynamics of the filling material. |
| Initialization | |||
| AbsolutePressure | p_start | Medium.p_default | Start value of pressure [Pa] |
| Temperature | TFlu_start[nSeg] | TGro_start | Start value of fluid temperature [K] |
| Soil | |||
| Temperature | TExt0_start | 283.15 | Initial far field temperature [K] |
| Temperature | TExt_start[nSeg] | {if z[i] >= z0 then TExt0_st... | Temperature of the undisturbed ground [K] |
| Filling material | |||
| Temperature | TGro_start[nSeg] | TExt_start | Start value of grout temperature [K] |
| Temperature profile | |||
| Height | z0 | 10 | Depth below which the temperature gradient starts [m] |
| Real | dT_dz | 0.01 | Vertical temperature gradient of the undisturbed soil for h below z0 [K/m] |
Connectors
| Type | Name | Description |
|---|---|---|
| FluidPort_a | port_a | Fluid connector a (positive design flow direction is from port_a to port_b) |
| FluidPort_b | port_b | Fluid connector b (positive design flow direction is from port_a to port_b) |
| output RealOutput | TBorAve | Average borehole wall temperature in the borefield [K] |
Modelica definition
replaceable model BoreFieldType=Fluid.Geothermal.Borefields.OneUTube
constrainedby Fluid.Geothermal.Borefields.BaseClasses.PartialBorefield(
redeclare package Medium=Medium,
allowFlowReversal=allowFlowReversal,
borFieDat=datBorFie,
energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial)
"Wall heat transfer";