Buildings.Experimental.DHC.Plants.Combined.Subsystems
Package of models for subsystems of combined plants
Information
This package contains models of components used in combined heating and cooling plants.
Extends from Modelica.Icons.VariantsPackage (Icon for package containing variants).
Package Content
| Name | Description |
|---|---|
 ChillerGroup
|
Model of multiple identical chillers in parallel |
 ChillerHeatRecoveryGroup
|
Model of multiple identical heat recovery chillers in parallel |
 CoolingTowerGroup
|
Model of multiple identical cooling towers in parallel |
 HeatPumpGroup
|
Model of multiple identical air-source heat pumps in parallel |
 MultiplePumpsDp
|
Model of multiple identical pumps in parallel with dp-controlled pump model |
| MultiplePumpsFlow
|
Model of multiple identical pumps in parallel with flow-controlled pump model |
| MultiplePumpsSpeed
|
Model of multiple identical pumps in parallel with speed-controlled pump model |
 Validation
|
Package with validation models |
 BaseClasses
|
Package with base classes |
Buildings.Experimental.DHC.Plants.Combined.Subsystems.ChillerGroup
Model of multiple identical chillers in parallel
Information
This model represents a set of identical water-cooled compression chillers that are piped in parallel. Modulating isolation valves are included on condenser and evaporator side.
Control points
The following input and output points are available.
-
On/Off command
y1: DO signal dedicated to each unit, with a dimensionality of one -
CHW supply temperature setpoint
TSet: AO signal common to all units, with a dimensionality of zero -
Condenser and evaporator isolation valve commanded position
yVal(Con|Eva): AO signal dedicated to each unit, with a dimensionality of one -
Condenser and evaporator leaving temperature
T(Con|Eva)Lvg: AI signal dedicated to each unit, with a dimensionality of one -
Condenser and evaporator mass flow rate
m(Con|Eva)_flow: AI signal dedicated to each unit, with a dimensionality of one
Details
Modeling approach
In a parallel arrangement, all operating units have the same operating point, provided that the isolation valves are commanded to the same position. This allows modeling the heat transfer through the condenser and evaporator barrel with a single instance of Buildings.Fluid.Chillers.ElectricReformulatedEIR. Hydronics are resolved with mass flow rate multiplier components in conjunction with instances of Buildings.Experimental.DHC.Plants.Combined.Subsystems.BaseClasses.MultipleValves which represent the parallel network of valves and fixed resistances.
Actuators
By default, linear valve models are used. Those are configured with a pressure drop varying linearly with the flow rate, as opposed to the quadratic dependency usually considered for a turbulent flow regime. This is because the whole plant model contains large nonlinear systems of equations and this configuration limits the risk of solver failure while reducing the time to solution. This has no significant impact on the operating point of the circulation pumps when a control loop is used to modulate the valve opening and maintain the flow rate or the leaving temperature at setpoint. Then, whatever the modeling assumptions for the valve, the control loop ensures that the valve creates the adequate pressure drop and flow, which will simply be reached at a different valve opening with the above simplification.
Extends from Buildings.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 | |
| Integer | nUni | Number of units operating at design conditions | |
| Generic | dat | redeclare parameter Fluid.Ch... | Chiller parameters (each unit) |
| Nominal condition | |||
| MassFlowRate | m1_flow_nominal | mConWat_flow_nominal | Nominal mass flow rate [kg/s] |
| MassFlowRate | m2_flow_nominal | mChiWat_flow_nominal | Nominal mass flow rate [kg/s] |
| PressureDifference | dpEva_nominal | Design chiller evaporator pressure drop (each unit) [Pa] | |
| PressureDifference | dpCon_nominal | Chiller condenser design pressure drop (each unit) [Pa] | |
| 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 |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Filtered opening | |||
| Boolean | use_inputFilter | energyDynamics <> Modelica.F... | = true, if opening is filtered with a 2nd order CriticalDamping filter |
| Time | riseTime | 120 | Rise time of the filter (time to reach 99.6 % of an opening step) [s] |
| Init | init | Modelica.Blocks.Types.Init.I... | Type of initialization (no init/steady state/initial state/initial output) |
| Real | y_start | 1 | Initial position of actuator |
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 BooleanInput | y1[nUni] | Chiller On/Off command |
| input RealInput | TSet | Supply temperature setpoint [K] |
| output RealOutput | P | Power drawn [W] |
| input RealInput | yValCon[nUni] | Chiller condenser isolation valve commanded position [1] |
| input RealInput | yValEva[nUni] | Chiller evaporator isolation valve commanded position [1] |
| output RealOutput | TConLvg[nUni] | Chiller condenser leaving temperature [K] |
| output RealOutput | mCon_flow[nUni] | Chiller condenser barrel mass flow rate [kg/s] |
| output RealOutput | TEvaLvg[nUni] | Chiller evaporator leaving temperature [K] |
| output RealOutput | mEva_flow[nUni] | Chiller evaporator barrel mass flow rate [kg/s] |
Modelica definition
model ChillerGroup
"Model of multiple identical chillers in parallel"
extends Buildings.Fluid.Interfaces.PartialFourPortInterface(
final m1_flow_nominal = mConWat_flow_nominal,
final m2_flow_nominal = mChiWat_flow_nominal);
parameter Integer nUni(final min=1, start=1)
"Number of units operating at design conditions";
final parameter Modelica.Units.SI.Temperature TChiWatSup_nominal=
dat.TEvaLvg_nominal
"Design (minimum) CHW supply temperature";
final parameter Modelica.Units.SI.Temperature TConWatLvg_nominal=
dat.TConLvg_nominal
"Design (maximum) CW leaving temperature";
final parameter Modelica.Units.SI.HeatFlowRate QChiWatUni_flow_nominal=
dat.QEva_flow_nominal
"Design cooling heat flow rate (each unit, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QConWatUni_flow_nominal=
-dat.QEva_flow_nominal * (1 + 1 / dat.COP_nominal * dat.etaMotor)
"Design CW heat flow rate (each unit, >0)";
final parameter Modelica.Units.SI.HeatFlowRate QChiWat_flow_nominal=
nUni * QChiWatUni_flow_nominal
"Design cooling heat flow rate (all units, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QConWat_flow_nominal=
nUni * QConWatUni_flow_nominal
"Design CW heat flow rate (all units, >0)";
final parameter Modelica.Units.SI.MassFlowRate mChiWatUni_flow_nominal(
final min=0)=dat.mEva_flow_nominal
"Design chiller CHW mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mConWatUni_flow_nominal(
final min=0)=dat.mCon_flow_nominal
"Design chiller CW mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mChiWat_flow_nominal(
final min=0)=nUni * mChiWatUni_flow_nominal
"Design CHW mass flow rate (all units)";
final parameter Modelica.Units.SI.MassFlowRate mConWat_flow_nominal(
final min=0)=nUni * mConWatUni_flow_nominal
"Design CW mass flow rate (all units)";
parameter Modelica.Units.SI.PressureDifference dpEva_nominal(displayUnit="Pa")
"Design chiller evaporator pressure drop (each unit)";
parameter Modelica.Units.SI.PressureDifference dpCon_nominal(displayUnit="Pa")
"Chiller condenser design pressure drop (each unit)";
replaceable parameter Fluid.Chillers.Data.ElectricReformulatedEIR.Generic dat
"Chiller parameters (each unit)";
// Assumptions
parameter Modelica.Fluid.Types.Dynamics energyDynamics=
Modelica.Fluid.Types.Dynamics.DynamicFreeInitial
"Type of energy balance: dynamic (3 initialization options) or steady state";
parameter Boolean use_inputFilter=energyDynamics<>Modelica.Fluid.Types.Dynamics.SteadyState
"= true, if opening is filtered with a 2nd order CriticalDamping filter";
parameter Modelica.Units.SI.Time riseTime=120
"Rise time of the filter (time to reach 99.6 % of an opening step)";
parameter Modelica.Blocks.Types.Init init=Modelica.Blocks.Types.Init.InitialOutput
"Type of initialization (no init/steady state/initial state/initial output)";
parameter Real y_start=1 "Initial position of actuator";
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput y1[nUni]
"Chiller On/Off command";
Buildings.Controls.OBC.CDL.Interfaces.RealInput TSet(
final unit="K", displayUnit="degC")
"Supply temperature setpoint";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput P(
final unit="W")
"Power drawn";
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValCon[nUni](
each final unit="1", each final min=0, each final max=1)
"Chiller condenser isolation valve commanded position";
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValEva[nUni](
each final unit="1", each final min=0, each final max=1)
"Chiller evaporator isolation valve commanded position";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput TConLvg[nUni](each final
unit="K", each displayUnit="degC")
"Chiller condenser leaving temperature";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput mCon_flow[nUni](
each final unit="kg/s")
"Chiller condenser barrel mass flow rate";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput TEvaLvg[nUni](each final
unit="K", each displayUnit="degC")
"Chiller evaporator leaving temperature";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput mEva_flow[nUni](
each final unit="kg/s")
"Chiller evaporator barrel mass flow rate";
Fluid.Chillers.ElectricReformulatedEIR chi(
PLR1(start=0),
final per=dat,
redeclare final package Medium1=Medium1,
redeclare final package Medium2=Medium2,
final have_switchover=false,
final dp1_nominal=0,
final dp2_nominal=0,
final allowFlowReversal1=allowFlowReversal1,
final allowFlowReversal2=allowFlowReversal2,
final energyDynamics=energyDynamics,
final show_T=show_T) "Chiller";
Fluid.BaseClasses.MassFlowRateMultiplier mulConInl(
redeclare final package Medium = Medium1,
final allowFlowReversal=allowFlowReversal1,
final use_input=true)
"Flow rate multiplier";
Fluid.BaseClasses.MassFlowRateMultiplier mulConOut(
redeclare final package Medium = Medium1,
final allowFlowReversal=allowFlowReversal1,
final use_input=true)
"Flow rate multiplier";
Fluid.BaseClasses.MassFlowRateMultiplier mulEvaInl(
redeclare final package Medium = Medium2,
final allowFlowReversal=allowFlowReversal2,
final use_input=true)
"Flow rate multiplier";
Fluid.BaseClasses.MassFlowRateMultiplier mulEvaOut(
redeclare final package Medium = Medium2,
final allowFlowReversal=allowFlowReversal2,
final use_input=true)
"Flow rate multiplier";
BaseClasses.MultipleCommands com(
final nUni=nUni)
"Convert command signals";
BaseClasses.MultipleValves valEva(
redeclare final package Medium = Medium2,
redeclare final Buildings.Fluid.Actuators.Valves.TwoWayLinear val,
linearized=true,
final nUni=nUni,
final mUni_flow_nominal=mChiWatUni_flow_nominal,
final dpFixed_nominal=dpEva_nominal,
dpValve_nominal=1E3,
final allowFlowReversal=allowFlowReversal2,
final energyDynamics=energyDynamics,
final use_inputFilter=use_inputFilter,
final riseTime=riseTime,
final init=init,
final y_start=y_start,
final show_T=show_T)
"Chiller evaporator isolation valves";
BaseClasses.MultipleValves valCon(
redeclare final package Medium = Medium1,
redeclare final Buildings.Fluid.Actuators.Valves.TwoWayLinear val,
linearized=true,
final nUni=nUni,
final mUni_flow_nominal=mConWatUni_flow_nominal,
final dpFixed_nominal=dpCon_nominal,
dpValve_nominal=1E3,
final allowFlowReversal=allowFlowReversal1,
final energyDynamics=energyDynamics,
final use_inputFilter=use_inputFilter,
final riseTime=riseTime,
final init=init,
final y_start=y_start,
final show_T=show_T)
"Chiller condenser isolation valves";
Buildings.Controls.OBC.CDL.Reals.Multiply mulP "Scale power";
Fluid.Sensors.TemperatureTwoPort temConLvg(
redeclare final package Medium = Medium1,
final m_flow_nominal=mConWatUni_flow_nominal,
final allowFlowReversal=allowFlowReversal1)
"Chiller condenser leaving temperature";
Fluid.Sensors.TemperatureTwoPort temEvaLvg(
redeclare final package Medium = Medium2,
final m_flow_nominal=mChiWatUni_flow_nominal,
final allowFlowReversal=allowFlowReversal2)
"Chiller evaporator leaving temperature";
Fluid.Sensors.MassFlowRate floCon(
redeclare final package Medium =Medium1,
final allowFlowReversal=allowFlowReversal1)
"Chiller condenser barrel mass flow rate";
Fluid.Sensors.MassFlowRate floEva(
redeclare final package Medium =Medium2,
final allowFlowReversal=allowFlowReversal2)
"Chiller evaporator barrel mass flow rate";
Buildings.Controls.OBC.CDL.Routing.RealScalarReplicator rep(nout=nUni)
"Replicate";
Buildings.Controls.OBC.CDL.Routing.RealScalarReplicator rep1(nout=nUni)
"Replicate";
Buildings.Controls.OBC.CDL.Routing.RealScalarReplicator rep2(nout=nUni)
"Replicate";
Buildings.Controls.OBC.CDL.Routing.RealScalarReplicator rep3(nout=nUni)
"Replicate";
equation
connect(mulEvaInl.port_b,chi. port_a2);
connect(TSet,chi. TSet);
connect(mulConOut.uInv, mulConInl.u);
connect(mulEvaOut.uInv, mulEvaInl.u);
connect(port_a1, mulConInl.port_a);
connect(port_a2, mulEvaInl.port_a);
connect(port_b2, valEva.port_b);
connect(mulEvaOut.port_b, valEva.port_a);
connect(y1, com.y1);
connect(chi.P, mulP.u1);
connect(com.nUniOn, mulP.u2);
connect(yValCon, valCon.y);
connect(yValEva, valEva.y);
connect(mulConInl.port_b,chi. port_a1);
connect(mulP.y, P);
connect(com.nUniOnBou, mulConOut.u);
connect(mulConInl.uInv, mulEvaOut.u);
connect(com.y1One,chi. on);
connect(valCon.port_b, port_b1);
connect(mulConOut.port_b, valCon.port_a);
connect(temEvaLvg.port_a, chi.port_b2);
connect(chi.port_b1, temConLvg.port_a);
connect(temConLvg.port_b, floCon.port_a);
connect(floCon.port_b, mulConOut.port_a);
connect(temEvaLvg.port_b, floEva.port_a);
connect(floEva.port_b, mulEvaOut.port_a);
connect(rep.y, TConLvg);
connect(rep1.y, mCon_flow);
connect(floCon.m_flow, rep1.u);
connect(temConLvg.T, rep.u);
connect(rep3.y, TEvaLvg);
connect(rep2.y, mEva_flow);
connect(floEva.m_flow, rep2.u);
connect(temEvaLvg.T, rep3.u);
end ChillerGroup;
Buildings.Experimental.DHC.Plants.Combined.Subsystems.ChillerHeatRecoveryGroup
Model of multiple identical heat recovery chillers in parallel
Information
This model represents a set of identical heat recovery chillers that are piped in parallel. Modulating isolation valves and modulating switchover valves are included on condenser and evaporator side. The switchover valves allow indexing the condenser (resp. the evaporator) either to the CWC loop or to the HW loop (resp. to the CWE loop or to the CHW loop). Modulating valves are used to allow for sequences of operation that bleed CWE into the HW return flow to modulate the condenser entering temperature.
Control points
The following input and output points are available.
-
On/Off command
y1: DO signal dedicated to each unit, with a dimensionality of one -
Cooling switchover command
y1Coo: DO signal dedicated to each unit, with a dimensionality of one -
Supply temperature setpoint
TSet: AO signal dedicated to each unit, with a dimensionality of one The signal corresponds either to the HW supply temperature setpoint when the unit operates in heating mode, or to the CHW supply temperature setpoint when the unit operates in cooling mode. -
Condenser and evaporator isolation valve commanded position
yVal(Con|Eva): AO signal dedicated to each unit, with a dimensionality of one -
Condenser and evaporator switchover valve commanded position
yVal(Con|Eva)Swi: AO signal dedicated to each unit, with a dimensionality of one -
Condenser and evaporator leaving temperature
T(Con|Eva)Lvg: AI signal dedicated to each unit, with a dimensionality of one -
Condenser entering temperature
TConEnt: AI signal dedicated to each unit, with a dimensionality of one -
Condenser and evaporator mass flow rate
m(Con|Eva)_flow: AI signal dedicated to each unit, with a dimensionality of one
Details
HRC performance data
The performance data should cover the HRC lift envelope,
that is when the HRC is operating in direct heat recovery mode,
producing CHW and HW at their setpoint value at full load.
In this case, and to allow for cascading heat recovery where
a third fluid circuit is used to generate a cascade of thermodynamic cycles,
two additional parameters TCasEntCoo_nominal and
TCasEntHea_nominal are exposed to specify the
entering temperature of the third fluid circuit when
the HRC is operating in cooling mode and in heating mode,
respectively.
In cooling mode the third fluid circuit is connected to the
condenser barrel.
In heating mode, the third fluid circuit is connected to the
evaporator barrel.
The parameters TCasEnt*_nominal are then used to assess the
design capacity in heating and cooling mode, respectively.
Actuators
By default, linear valve models are used. Those are configured with a pressure drop varying linearly with the flow rate, as opposed to the quadratic dependency usually considered for a turbulent flow regime. This is because the whole plant model contains large nonlinear systems of equations and this configuration limits the risk of solver failure while reducing the time to solution. This has no significant impact on the operating point of the circulation pumps when a control loop is used to modulate the valve opening and maintain the flow rate or the leaving temperature at setpoint. Then, whatever the modeling assumptions for the valve, the control loop ensures that the valve creates the adequate pressure drop and flow, which will simply be reached at a different valve opening with the above simplification.
Extends from Buildings.Fluid.Interfaces.PartialEightPortInterface (Partial model transporting fluid between eight 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 | |
| replaceable package Medium3 | PartialMedium | Medium 3 in the component | |
| replaceable package Medium4 | PartialMedium | Medium 4 in the component | |
| replaceable package Medium | Buildings.Media.Water | Medium model for all four fluid circuits | |
| Integer | nUni | Number of units operating at design conditions | |
| Temperature | TCasHeaEnt_nominal | 298.15 | Design evaporator entering temperature in cascading heating mode [K] |
| Temperature | TCasCooEnt_nominal | 288.15 | Design condenser entering temperature in cascading cooling mode [K] |
| Generic | dat | redeclare parameter Fluid.Ch... | Chiller parameters (each unit) |
| Nominal condition | |||
| MassFlowRate | m1_flow_nominal | mConWat_flow_nominal | Nominal mass flow rate [kg/s] |
| MassFlowRate | m2_flow_nominal | mChiWat_flow_nominal | Nominal mass flow rate [kg/s] |
| MassFlowRate | m3_flow_nominal | mConWat_flow_nominal | Nominal mass flow rate [kg/s] |
| MassFlowRate | m4_flow_nominal | mChiWat_flow_nominal | Nominal mass flow rate [kg/s] |
| PressureDifference | dpEva_nominal | Design evaporator pressure drop (each unit) [Pa] | |
| PressureDifference | dpCon_nominal | Design condenser pressure drop (each unit) [Pa] | |
| Assumptions | |||
| Boolean | allowFlowReversal1 | allowFlowReversal | = true to allow flow reversal in medium 1, false restricts to design direction (port_a -> port_b) |
| Boolean | allowFlowReversal2 | allowFlowReversal | = true to allow flow reversal in medium 2, false restricts to design direction (port_a -> port_b) |
| Boolean | allowFlowReversal3 | allowFlowReversal | = true to allow flow reversal in medium 3, false restricts to design direction (port_a -> port_b) |
| Boolean | allowFlowReversal4 | allowFlowReversal | = true to allow flow reversal in medium 4, false restricts to design direction (port_a -> port_b) |
| Boolean | allowFlowReversal | true | = false to simplify equations, assuming, but not enforcing, no flow reversal |
| Advanced | |||
| Initialization | |||
| SpecificEnthalpy | h_outflow_a1_start | Medium1.h_default | Start value for enthalpy flowing out of port a1 [J/kg] |
| SpecificEnthalpy | h_outflow_b1_start | Medium1.h_default | Start value for enthalpy flowing out of port b1 [J/kg] |
| SpecificEnthalpy | h_outflow_a2_start | Medium2.h_default | Start value for enthalpy flowing out of port a2 [J/kg] |
| SpecificEnthalpy | h_outflow_b2_start | Medium2.h_default | Start value for enthalpy flowing out of port b2 [J/kg] |
| SpecificEnthalpy | h_outflow_a3_start | Medium3.h_default | Start value for enthalpy flowing out of port a1 [J/kg] |
| SpecificEnthalpy | h_outflow_b3_start | Medium3.h_default | Start value for enthalpy flowing out of port b1 [J/kg] |
| SpecificEnthalpy | h_outflow_a4_start | Medium4.h_default | Start value for enthalpy flowing out of port a1 [J/kg] |
| SpecificEnthalpy | h_outflow_b4_start | Medium4.h_default | Start value for enthalpy flowing out of port b1 [J/kg] |
| 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] |
| MassFlowRate | m3_flow_small | 1E-4*abs(m3_flow_nominal) | Small mass flow rate for regularization of zero flow [kg/s] |
| MassFlowRate | m4_flow_small | 1E-4*abs(m4_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 |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Filtered opening | |||
| Boolean | use_inputFilter | energyDynamics <> Modelica.F... | = true, if opening is filtered with a 2nd order CriticalDamping filter |
| Time | riseTime | 120 | Rise time of the filter (time to reach 99.6 % of an opening step) [s] |
| Init | init | Modelica.Blocks.Types.Init.I... | Type of initialization (no init/steady state/initial state/initial output) |
| Real | y_start | 1 | Initial position of actuator |
Connectors
| Type | Name | Description |
|---|---|---|
| replaceable package Medium1 | Medium 1 in the component | |
| replaceable package Medium2 | Medium 2 in the component | |
| replaceable package Medium3 | Medium 3 in the component | |
| replaceable package Medium4 | Medium 4 in the component | |
| 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) |
| FluidPort_a | port_a3 | Fluid connector a1 (positive design flow direction is from port_a3 to port_b3) |
| FluidPort_b | port_b3 | Fluid connector b2 (positive design flow direction is from port_a3 to port_b3) |
| FluidPort_a | port_a4 | Fluid connector a1 (positive design flow direction is from port_a4 to port_b4) |
| FluidPort_b | port_b4 | Fluid connector b2 (positive design flow direction is from port_a4 to port_b4) |
| replaceable package Medium | Medium model for all four fluid circuits | |
| input BooleanInput | y1[nUni] | Chiller On/Off command |
| input BooleanInput | y1Coo[nUni] | Cooling switchover command: true for cooling, false for heating |
| input RealInput | TSet[nUni] | Supply temperature setpoint [K] |
| input RealInput | yValConSwi[nUni] | Chiller condenser switchover valve commanded position [1] |
| input RealInput | yValEvaSwi[nUni] | Chiller evaporator switchover valve commanded position [1] |
| output RealOutput | P | Power drawn [W] |
| input RealInput | yValCon[nUni] | Chiller condenser isolation valve commanded position [1] |
| input RealInput | yValEva[nUni] | Chiller evaporator isolation valve commanded position [1] |
| output RealOutput | mEva_flow[nUni] | Chiller evaporator barrel mass flow rate [kg/s] |
| output RealOutput | TEvaLvg[nUni] | Chiller evaporator leaving temperature [K] |
| output RealOutput | mCon_flow[nUni] | Chiller condenser barrel mass flow rate [kg/s] |
| output RealOutput | TConLvg[nUni] | Chiller condenser leaving temperature [K] |
| output RealOutput | TConEnt[nUni] | Chiller condenser entering temperature [K] |
Modelica definition
model ChillerHeatRecoveryGroup
"Model of multiple identical heat recovery chillers in parallel"
extends Buildings.Fluid.Interfaces.PartialEightPortInterface(
redeclare final package Medium1=Medium,
redeclare final package Medium2=Medium,
redeclare final package Medium3=Medium,
redeclare final package Medium4=Medium,
final m1_flow_nominal = mConWat_flow_nominal,
final m2_flow_nominal = mChiWat_flow_nominal,
final m3_flow_nominal = mConWat_flow_nominal,
final m4_flow_nominal = mChiWat_flow_nominal,
final allowFlowReversal1=allowFlowReversal,
final allowFlowReversal2=allowFlowReversal,
final allowFlowReversal3=allowFlowReversal,
final allowFlowReversal4=allowFlowReversal);
replaceable package Medium=Buildings.Media.Water
constrainedby Modelica.Media.Interfaces.PartialMedium
"Medium model for all four fluid circuits";
parameter Integer nUni(final min=1, start=1)
"Number of units operating at design conditions";
parameter Modelica.Units.SI.Temperature TCasHeaEnt_nominal=298.15
"Design evaporator entering temperature in cascading heating mode";
parameter Modelica.Units.SI.Temperature TCasCooEnt_nominal=288.15
"Design condenser entering temperature in cascading cooling mode";
final parameter Modelica.Units.SI.Temperature TChiWatSup_nominal=
dat.TEvaLvg_nominal
"Design (minimum) CHW supply temperature";
final parameter Modelica.Units.SI.Temperature THeaWatSup_nominal=
dat.TConLvg_nominal
"Design (maximum) HW supply temperature";
final parameter Modelica.Units.SI.HeatFlowRate QChiWatUni_flow_nominal=
dat.QEva_flow_nominal
"Design cooling heat flow rate (each unit, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QHeaWatUni_flow_nominal=
-dat.QEva_flow_nominal * (1 + 1 / dat.COP_nominal * dat.etaMotor)
"Design heating heat flow rate in direct heat recovery mode (each unit, >0)";
final parameter Modelica.Units.SI.Efficiency COPCasHea_nominal(fixed=false)
"Coefficient of performance in cascading heating mode";
final parameter Modelica.Units.SI.Efficiency COPCasCoo_nominal(fixed=false)
"Coefficient of performance in cascading cooling mode";
final parameter Modelica.Units.SI.Temperature TCasHeaLvg_nominal(fixed=false)
"Design value of evaporator leaving temperature in cascading heating mode";
final parameter Modelica.Units.SI.Temperature TCasCooLvg_nominal(fixed=false)
"Design value of condenser leaving temperature in cascading cooling mode";
final parameter Modelica.Units.SI.HeatFlowRate QEvaCasHeaUni_flow_nominal(
fixed=false,
start=QChiWatUni_flow_nominal * (1 + 2E-2 * (TCasHeaEnt_nominal -
(dat.TEvaLvg_nominal - QChiWatUni_flow_nominal / mChiWatUni_flow_nominal / cpCas))))
"Design evaporator heat flow rate in cascading heating mode (each unit, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QChiWatCasCooUni_flow_nominal(
fixed=false,
start=QChiWatUni_flow_nominal * (1 - 2E-2 * (TCasCooEnt_nominal -
(dat.TConLvg_nominal - QHeaWatUni_flow_nominal / mConWatUni_flow_nominal / cpCas))))
"Design cooling heat flow rate in cascading cooling mode (each unit, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QHeaWatCasHeaUni_flow_nominal=
-QEvaCasHeaUni_flow_nominal * (1 + 1 / COPCasHea_nominal * dat.etaMotor)
"Design heating heat flow rate in cascading heating mode (each unit, >0)";
final parameter Modelica.Units.SI.HeatFlowRate QConCasCooUni_flow_nominal=
-QChiWatCasCooUni_flow_nominal * (1 + 1 / COPCasCoo_nominal * dat.etaMotor)
"Design condenser heat flow rate in cascading cooling mode (each unit, >0)";
final parameter Modelica.Units.SI.HeatFlowRate QChiWat_flow_nominal=
nUni * QChiWatUni_flow_nominal
"Design cooling heat flow rate (all units, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QHeaWat_flow_nominal=
nUni * QHeaWatUni_flow_nominal
"Design heating heat flow rate (all units, >0)";
final parameter Modelica.Units.SI.HeatFlowRate QEvaCasHea_flow_nominal=
nUni * QEvaCasHeaUni_flow_nominal
"Design evaporator heat flow rate in cascading heating mode (all units, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QChiWatCasCoo_flow_nominal=
nUni * QChiWatCasCooUni_flow_nominal
"Design cooling heat flow rate in cascading cooling mode (all units, <0)";
final parameter Modelica.Units.SI.HeatFlowRate QHeaWatCasHea_flow_nominal=
nUni * QHeaWatCasHeaUni_flow_nominal
"Design heating heat flow rate in cascading heating mode (all units, >0)";
final parameter Modelica.Units.SI.HeatFlowRate QConCasCoo_flow_nominal=
nUni * QConCasCooUni_flow_nominal
"Design condenser heat flow rate in cascading cooling mode (all units, >0)";
final parameter Modelica.Units.SI.MassFlowRate mChiWatUni_flow_nominal(
final min=0)=dat.mEva_flow_nominal
"Design CHW mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mConWatUni_flow_nominal(
final min=0)=dat.mCon_flow_nominal
"Design CW mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mChiWat_flow_nominal(
final min=0)=nUni * mChiWatUni_flow_nominal
"Design CHW mass flow rate (all units)";
final parameter Modelica.Units.SI.MassFlowRate mConWat_flow_nominal(
final min=0)=nUni * mConWatUni_flow_nominal
"Design CW mass flow rate (all units)";
parameter Modelica.Units.SI.PressureDifference dpEva_nominal(displayUnit="Pa")
"Design evaporator pressure drop (each unit)";
parameter Modelica.Units.SI.PressureDifference dpCon_nominal(displayUnit="Pa")
"Design condenser pressure drop (each unit)";
replaceable parameter Fluid.Chillers.Data.ElectricReformulatedEIR.Generic dat
"Chiller parameters (each unit)";
// Assumptions
parameter Boolean allowFlowReversal = true
"= false to simplify equations, assuming, but not enforcing, no flow reversal";
parameter Modelica.Fluid.Types.Dynamics energyDynamics=
Modelica.Fluid.Types.Dynamics.DynamicFreeInitial
"Type of energy balance: dynamic (3 initialization options) or steady state";
parameter Boolean use_inputFilter=energyDynamics<>Modelica.Fluid.Types.Dynamics.SteadyState
"= true, if opening is filtered with a 2nd order CriticalDamping filter";
parameter Modelica.Units.SI.Time riseTime=120
"Rise time of the filter (time to reach 99.6 % of an opening step)";
parameter Modelica.Blocks.Types.Init init=Modelica.Blocks.Types.Init.InitialOutput
"Type of initialization (no init/steady state/initial state/initial output)";
parameter Real y_start=1 "Initial position of actuator";
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput y1[nUni]
"Chiller On/Off command";
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput y1Coo[nUni]
"Cooling switchover command: true for cooling, false for heating";
Buildings.Controls.OBC.CDL.Interfaces.RealInput TSet[nUni](
each final unit="K", each displayUnit="degC")
"Supply temperature setpoint";
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValConSwi[nUni](
each final unit="1",
each final min=0,
each final max=1) "Chiller condenser switchover valve commanded position";
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValEvaSwi[nUni](
each final unit="1",
each final min=0,
each final max=1) "Chiller evaporator switchover valve commanded position";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput P(
final unit="W") "Power drawn";
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValCon[nUni](
each final unit="1", each final min=0, each final max=1)
"Chiller condenser isolation valve commanded position";
Buildings.Controls.OBC.CDL.Interfaces.RealInput yValEva[nUni](
each final unit="1", each final min=0, each final max=1)
"Chiller evaporator isolation valve commanded position";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput mEva_flow[nUni](each final
unit="kg/s") "Chiller evaporator barrel mass flow rate";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput TEvaLvg[nUni](
each final unit="K", each displayUnit="degC")
"Chiller evaporator leaving temperature";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput mCon_flow[nUni](
each final unit="kg/s")
"Chiller condenser barrel mass flow rate";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput TConLvg[nUni](
each final unit="K", each displayUnit="degC")
"Chiller condenser leaving temperature";
Fluid.Chillers.ElectricReformulatedEIR chi[nUni](
PLR1(each start=0),
each final per=dat,
redeclare each final package Medium1 = Medium,
redeclare final package Medium2 = Medium,
each final have_switchover=true,
each final dp1_nominal=0,
each final dp2_nominal=0,
each final allowFlowReversal1=allowFlowReversal,
each final allowFlowReversal2=allowFlowReversal,
each final energyDynamics=energyDynamics,
each final show_T=show_T) "Chiller";
Fluid.Actuators.Valves.TwoWayLinear valCon[nUni](
redeclare each final package Medium = Medium,
each linearized=true,
each final CvData=Buildings.Fluid.Types.CvTypes.OpPoint,
each final m_flow_nominal=mConWatUni_flow_nominal,
each dpValve_nominal=1E3,
each final dpFixed_nominal=dpCon_nominal,
each final allowFlowReversal=allowFlowReversal,
each final use_inputFilter=use_inputFilter,
each final riseTime=riseTime,
each final init=init,
each final y_start=y_start) "Condenser isolation valve";
Fluid.Actuators.Valves.TwoWayLinear valEva[nUni](
redeclare each final package Medium = Medium,
each linearized=true,
each final CvData=Buildings.Fluid.Types.CvTypes.OpPoint,
each final m_flow_nominal=mChiWatUni_flow_nominal,
each dpValve_nominal=1E3,
each final dpFixed_nominal=dpEva_nominal,
each final allowFlowReversal=allowFlowReversal,
each final use_inputFilter=use_inputFilter,
each final riseTime=riseTime,
each final init=init,
each final y_start=y_start) "Evaporator isolation valve";
Fluid.FixedResistances.Junction junConWatEvaOut[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mChiWatUni_flow_nominal*{1,-1,1},
each final dp_nominal=fill(1E3, 3),
each final energyDynamics=energyDynamics,
each final portFlowDirection_1=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_2=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_3=if allowFlowReversal then Modelica.Fluid.Types.PortFlowDirection.Bidirectional
else Modelica.Fluid.Types.PortFlowDirection.Entering) "Fluid junction";
Fluid.FixedResistances.Junction junConWatEvaInl[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mChiWatUni_flow_nominal*{1,-1,-1},
each final dp_nominal=fill(1E3, 3),
each final energyDynamics=energyDynamics,
each final portFlowDirection_1=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_2=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_3=if allowFlowReversal then Modelica.Fluid.Types.PortFlowDirection.Bidirectional
else Modelica.Fluid.Types.PortFlowDirection.Leaving) "Fluid junction";
Fluid.Sensors.TemperatureTwoPort temEvaLvg[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mChiWatUni_flow_nominal,
each final allowFlowReversal=allowFlowReversal)
"Chiller evaporator leaving temperature";
Fluid.FixedResistances.Junction junHeaWatConInl[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mConWatUni_flow_nominal*{1,-1,-1},
each final dp_nominal=fill(1E3, 3),
each final energyDynamics=energyDynamics,
each final portFlowDirection_1=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_2=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_3=if allowFlowReversal then Modelica.Fluid.Types.PortFlowDirection.Bidirectional
else Modelica.Fluid.Types.PortFlowDirection.Leaving) "Fluid junction";
Fluid.FixedResistances.Junction junHeaWatConOut[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mConWatUni_flow_nominal*{1,-1,1},
each final dp_nominal=fill(1E3, 3),
each final energyDynamics=energyDynamics,
each final portFlowDirection_1=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_2=Modelica.Fluid.Types.PortFlowDirection.Bidirectional,
each final portFlowDirection_3=if allowFlowReversal then Modelica.Fluid.Types.PortFlowDirection.Bidirectional
else Modelica.Fluid.Types.PortFlowDirection.Entering) "Fluid junction";
Fluid.Sensors.TemperatureTwoPort temConEnt[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mConWatUni_flow_nominal,
each final allowFlowReversal=allowFlowReversal)
"Chiller condenser entering temperature";
Fluid.Sensors.TemperatureTwoPort temConLvg[nUni](
redeclare each final package Medium = Medium,
each final m_flow_nominal=mConWatUni_flow_nominal,
each final allowFlowReversal=allowFlowReversal)
"Chiller condenser leaving temperature";
Buildings.Controls.OBC.CDL.Reals.MultiSum mulSum(nin=nUni)
"Sum up power of all units";
Fluid.Sensors.MassFlowRate floEva[nUni](
redeclare each final package Medium =Medium,
each final allowFlowReversal=allowFlowReversal)
"Chiller evaporator barrel mass flow rate";
Fluid.Sensors.MassFlowRate floCon[nUni](
redeclare each final package Medium =Medium,
each final allowFlowReversal=allowFlowReversal)
"Chiller condenser barrel mass flow rate";
Fluid.Actuators.Valves.TwoWayLinear valConSwi[nUni](
redeclare each final package Medium = Medium,
each linearized=true,
each final CvData=Buildings.Fluid.Types.CvTypes.OpPoint,
each final m_flow_nominal=mConWatUni_flow_nominal,
each dpValve_nominal=1E3,
each final dpFixed_nominal=0,
each final allowFlowReversal=true,
each final use_inputFilter=use_inputFilter,
each final riseTime=riseTime,
each final init=init,
each final y_start=y_start) "Condenser switchover valve";
Fluid.Actuators.Valves.TwoWayLinear valEvaSwi[nUni](
redeclare each final package Medium = Medium,
each linearized=true,
each final CvData=Buildings.Fluid.Types.CvTypes.OpPoint,
each final m_flow_nominal=mChiWatUni_flow_nominal,
each dpValve_nominal=1E3,
each final dpFixed_nominal=0,
each final allowFlowReversal=true,
each final use_inputFilter=use_inputFilter,
each final riseTime=riseTime,
each final init=init,
each final y_start=y_start) "Evaporator switchover valve";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput TConEnt[nUni](each final
unit="K", each displayUnit="degC")
"Chiller condenser entering temperature";
protected
parameter Medium.ThermodynamicState staCas=Medium.setState_pTX(
T=(TCasHeaEnt_nominal + TCasCooEnt_nominal) / 2,
p=Medium.p_default,
X=Medium.X_default)
"State of source medium in cascading mode (at mean source temperature)";
parameter Modelica.Units.SI.SpecificHeatCapacity cpCas=
Medium.specificHeatCapacityCp(state=staCas)
"Heat capacity of source medium in cascading mode";
initial equation
TCasHeaLvg_nominal = TCasHeaEnt_nominal +
QEvaCasHeaUni_flow_nominal / cpCas / mConWatUni_flow_nominal;
TCasCooLvg_nominal = TCasCooEnt_nominal +
QConCasCooUni_flow_nominal / cpCas / mConWatUni_flow_nominal;
QEvaCasHeaUni_flow_nominal = dat.QEva_flow_nominal *
Buildings.Utilities.Math.Functions.biquadratic(
a=dat.capFunT,
x1=Modelica.Units.Conversions.to_degC(TCasHeaLvg_nominal),
x2=Modelica.Units.Conversions.to_degC(THeaWatSup_nominal));
QChiWatCasCooUni_flow_nominal = dat.QEva_flow_nominal *
Buildings.Utilities.Math.Functions.biquadratic(
a=dat.capFunT,
x1=Modelica.Units.Conversions.to_degC(TChiWatSup_nominal),
x2=Modelica.Units.Conversions.to_degC(TCasCooLvg_nominal));
COPCasHea_nominal = dat.COP_nominal /
Buildings.Utilities.Math.Functions.biquadratic(
a=dat.EIRFunT,
x1=Modelica.Units.Conversions.to_degC(TCasCooLvg_nominal),
x2=Modelica.Units.Conversions.to_degC(THeaWatSup_nominal));
COPCasCoo_nominal = dat.COP_nominal /
Buildings.Utilities.Math.Functions.biquadratic(
a=dat.EIRFunT,
x1=Modelica.Units.Conversions.to_degC(TChiWatSup_nominal),
x2=Modelica.Units.Conversions.to_degC(TCasCooLvg_nominal));
equation
if nUni > 1 then
connect(junHeaWatConOut[1:nUni - 1].port_1, junHeaWatConOut[2:nUni].port_2);
connect(valEvaSwi[1:nUni - 1].port_a, junConWatEvaInl[2:nUni].port_2);
connect(valConSwi[1:nUni - 1].port_a, junHeaWatConInl[2:nUni].port_2);
connect(junConWatEvaOut[1:nUni - 1].port_1, junConWatEvaOut[2:nUni].port_2);
end if;
connect(port_a4, valEvaSwi[nUni].port_a);
connect(yValConSwi, valConSwi.y);
connect(y1, chi.on);
connect(TSet, chi.TSet);
connect(y1Coo, chi.coo);
connect(yValCon, valCon.y);
connect(chi.port_b2, temEvaLvg.port_a);
connect(yValEva, valEva.y);
connect(temConEnt.port_b, chi.port_a1);
connect(chi.port_b1, temConLvg.port_a);
connect(valCon.port_b, junHeaWatConOut.port_3);
connect(valEva.port_b,junConWatEvaOut. port_3);
connect(chi.P, mulSum.u);
connect(temEvaLvg.port_b, floEva.port_a);
connect(floEva.port_b, valEva.port_a);
connect(mulSum.y, P);
connect(floEva.m_flow, mEva_flow);
connect(temEvaLvg.T, TEvaLvg);
connect(temConLvg.port_b, floCon.port_a);
connect(floCon.port_b, valCon.port_a);
connect(temConLvg.T, TConLvg);
connect(floCon.m_flow, mCon_flow);
connect(junConWatEvaInl.port_3, chi.port_a2);
connect(yValEvaSwi, valEvaSwi.y);
connect(junHeaWatConInl.port_3, temConEnt.port_a);
connect(port_a1, junHeaWatConInl[1].port_2);
connect(valConSwi.port_b, junHeaWatConInl.port_1);
connect(junHeaWatConOut[1].port_2, port_b1);
connect(valEvaSwi.port_b, junConWatEvaInl.port_1);
connect(port_a2, junConWatEvaInl[1].port_2);
connect(junHeaWatConOut[nUni].port_1, port_b3);
connect(junConWatEvaOut[1].port_2, port_b2);
connect(junConWatEvaOut[nUni].port_1, port_b4);
connect(valConSwi[nUni].port_a, port_a3);
connect(temConEnt.T, TConEnt);
end ChillerHeatRecoveryGroup;
Buildings.Experimental.DHC.Plants.Combined.Subsystems.CoolingTowerGroup
Model of multiple identical cooling towers in parallel
Information
This model represents a set of identical cooling towers that are piped in parallel. No actuated isolation valves are included.
Control points
The following input and output points are available.
-
Start command (VFD Run)
y1: DO signal dedicated to each unit, with a dimensionality of one -
Speed command
y: AO signal common to all units, with a dimensionality of zero -
CW supply temperature
TConWatSup: AI signal common to all units, with a dimensionality of zero
Details
Modeling approach
In a parallel arrangement, all operating units have the same operating point. This allows modeling the heat transfer from outdoor air to condenser water with a single instance of Buildings.Fluid.HeatExchangers.CoolingTowers.Merkel. Hydronics are resolved with mass flow rate multiplier components. Note that the elevation head of open-circuit cooling towers is not modeled, which is an inherent limitation of the cooling tower component model.
The fan cycling On and Off is implicitly modeled in the cooling tower component which uses a low limit of the control signal to switch to a free convection regime at zero fan power.
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 | |
| Integer | nUni | Number of units operating at design conditions | |
| Real | yFan_min | 0.1 | CT fan minimum speed [1] |
| Nominal condition | |||
| MassFlowRate | mConWatUni_flow_nominal | CW design mass flow rate (each unit) [kg/s] | |
| PressureDifference | dpConWatFriUni_nominal | CW flow-friction losses through tower and piping only (without elevation head or valve) [Pa] | |
| MassFlowRate | mAirUni_flow_nominal | Air mass flow rate (each unit) [kg/s] | |
| Temperature | TWetBulEnt_nominal | Entering air wetbulb temperature [K] | |
| Temperature | TConWatRet_nominal | CW return temperature (tower entering) [K] | |
| Temperature | TConWatSup_nominal | CW supply temperature (tower leaving) [K] | |
| Power | PFanUni_nominal | Fan power (each unit) [W] | |
| 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 |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Nominal condition | |||
| Time | tau | 30 | Time constant of fluid volume for nominal flow, used if energy or mass balance is dynamic [s] |
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) |
| input BooleanInput | y1[nUni] | Cooling tower Start command |
| input RealInput | y | Cooling tower speed (common to all tower cells) [1] |
| Bus | weaBus | Bus with weather data |
| output RealOutput | P | Power drawn by tower fans [W] |
| output RealOutput | TConWatSup | CW supply temperature (tower leaving) [K] |
Modelica definition
model CoolingTowerGroup "Model of multiple identical cooling towers in parallel"
extends Buildings.Fluid.Interfaces.PartialTwoPortInterface(
final m_flow_nominal(final min=Modelica.Constants.small)=mConWat_flow_nominal);
parameter Integer nUni(final min=1, start=1)
"Number of units operating at design conditions";
parameter Modelica.Units.SI.MassFlowRate mConWatUni_flow_nominal(
final min=0)
"CW design mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mConWat_flow_nominal(
final min=0)=nUni * mConWatUni_flow_nominal
"CW design mass flow rate (all units)";
parameter Modelica.Units.SI.PressureDifference dpConWatFriUni_nominal(
displayUnit="Pa")
"CW flow-friction losses through tower and piping only (without elevation head or valve)";
parameter Modelica.Units.SI.MassFlowRate mAirUni_flow_nominal(
final min=0,
start=mConWatUni_flow_nominal / 1.45)
"Air mass flow rate (each unit)";
parameter Modelica.Units.SI.Temperature TWetBulEnt_nominal(
final min=273.15)
"Entering air wetbulb temperature";
parameter Modelica.Units.SI.Temperature TConWatRet_nominal(
final min=273.15)
"CW return temperature (tower entering)";
parameter Modelica.Units.SI.Temperature TConWatSup_nominal(
final min=273.15)
"CW supply temperature (tower leaving)";
parameter Modelica.Units.SI.Power PFanUni_nominal(
final min=0,
start=340 * mConWatUni_flow_nominal)
"Fan power (each unit)";
parameter Real yFan_min(unit="1")=0.1
"CT fan minimum speed";
// Assumptions
parameter Modelica.Fluid.Types.Dynamics energyDynamics=
Modelica.Fluid.Types.Dynamics.DynamicFreeInitial
"Type of energy balance: dynamic (3 initialization options) or steady state";
parameter Modelica.Units.SI.Time tau=30
"Time constant of fluid volume for nominal flow, used if energy or mass balance is dynamic";
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput y1[nUni]
"Cooling tower Start command";
Buildings.Controls.OBC.CDL.Interfaces.RealInput y(
final unit="1", final min=0)
"Cooling tower speed (common to all tower cells)";
BoundaryConditions.WeatherData.Bus weaBus
"Bus with weather data";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput P(final unit="W")
"Power drawn by tower fans";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput TConWatSup(
final unit="K", displayUnit="degC")
"CW supply temperature (tower leaving)";
Fluid.BaseClasses.MassFlowRateMultiplier mulInl(
redeclare final package Medium = Medium,
final allowFlowReversal=allowFlowReversal,
final use_input=false,
final k=1/nUni) "Flow rate multiplier";
Fluid.BaseClasses.MassFlowRateMultiplier mulOut(
redeclare final package Medium = Medium,
final allowFlowReversal=allowFlowReversal,
final use_input=false,
final k=nUni) "Flow rate multiplier";
BaseClasses.MultipleCommands com(final nUni=nUni)
"Convert command signals";
Buildings.Controls.OBC.CDL.Reals.Multiply mulP "Scale power";
Fluid.HeatExchangers.CoolingTowers.Merkel coo(
redeclare final package Medium = Medium,
final m_flow_nominal=mConWatUni_flow_nominal,
final dp_nominal=dpConWatFriUni_nominal,
final yMin=yFan_min,
final ratWatAir_nominal=mConWatUni_flow_nominal/mAirUni_flow_nominal,
final TAirInWB_nominal=TWetBulEnt_nominal,
final TWatIn_nominal=TConWatRet_nominal,
final TWatOut_nominal=TConWatSup_nominal,
final PFan_nominal=PFanUni_nominal,
final tau=tau,
final allowFlowReversal=allowFlowReversal,
final show_T=show_T,
final energyDynamics=energyDynamics)
"Cooling tower";
Buildings.Controls.OBC.CDL.Logical.Pre preY1[nUni]
"Left limit of signal avoiding direct feedback of status to controller";
Buildings.Controls.OBC.CDL.Conversions.BooleanToReal booToRea
"Convert to real";
Buildings.Controls.OBC.CDL.Reals.Multiply inp
"Compute pump input signal";
equation
connect(com.nUniOn, mulP.u2);
connect(mulP.y, P);
connect(mulOut.port_b, port_b);
connect(mulInl.port_b, coo.port_a);
connect(coo.port_b, mulOut.port_a);
connect(weaBus.TWetBul, coo.TAir);
connect(coo.PFan, mulP.u1);
connect(coo.TLvg, TConWatSup);
connect(y1, preY1.u);
connect(preY1.y, com.y1);
connect(com.y1One, booToRea.u);
connect(inp.y, coo.y);
connect(y, inp.u1);
connect(booToRea.y, inp.u2);
connect(port_a, mulInl.port_a);
end CoolingTowerGroup;
Buildings.Experimental.DHC.Plants.Combined.Subsystems.HeatPumpGroup
Model of multiple identical air-source heat pumps in parallel
Information
This model represents a set of identical air-to-water heat pumps that are piped in parallel. Dedicated constant-speed condenser pumps are included.
Control points
The following input and output points are available.
-
On/Off command
y1: DO signal dedicated to each unit, with a dimensionality of one -
Supply temperature setpoint
TSet: AO signal common to all units, with a dimensionality of zero -
CW supply temperature
TConWatSup: AI signal common to all units, with a dimensionality of zero
Details
Modeling approach
In a parallel arrangement, all operating units have the same operating point. This allows modeling the heat transfer from outdoor air to condenser water with a single instance of Buildings.Fluid.HeatPumps.EquationFitReversible. Hydronics are resolved with mass flow rate multiplier components.
The model Buildings.Fluid.HeatPumps.EquationFitReversible does not capture the sensitivity of the heat pump performance to the HW supply temperature setpoint. This means that a varying HW supply temperature setpoint has no impact on the heat pump COP (all other variables such as the HW return temperature being kept invariant). This is a limitation of the model.
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 package MediumAir | Buildings.Media.Air | Air medium | |
| Integer | nUni | Number of units operating at design conditions | |
| Generic | dat | redeclare parameter Building... | Heat pump parameters (each unit) |
| Generic | datPum | redeclare parameter Fluid.Mo... | Condenser pump parameters (each unit) |
| Nominal condition | |||
| PressureDifference | dpPumHeaWatUni_nominal | dat.dpHeaLoa_nominal | HW pump design head (each unit) [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 |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Nominal condition | |||
| Time | tau | 30 | Time constant of fluid volume for nominal HW flow, used if energy or mass balance is dynamic [s] |
| Filtered pump speed | |||
| Boolean | use_inputFilter | energyDynamics <> Modelica.F... | = true, if signal is filtered with a 2nd order CriticalDamping filter |
| Time | riseTime | 30 | Rise time of the filter (time to reach 99.6 % of the speed) [s] |
| Init | init | Modelica.Blocks.Types.Init.I... | Type of initialization (no init/steady state/initial state/initial output) |
| Real | y_start | 1 | Initial position of actuator |
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 package MediumAir | Air medium | |
| input BooleanInput | y1[nUni] | Heat pump On/Off command |
| input RealInput | TSet | Supply temperature setpoint [K] |
| Bus | weaBus | Bus with weather data |
| output RealOutput | PPum | Power drawn by HW pumps [W] |
| output RealOutput | P | Power drawn by heat pumps [W] |
Modelica definition
model HeatPumpGroup
"Model of multiple identical air-source heat pumps in parallel"
extends Buildings.Fluid.Interfaces.PartialTwoPortInterface(
final m_flow_nominal(final min=Modelica.Constants.small)=mHeaWat_flow_nominal);
replaceable package MediumAir=Buildings.Media.Air
"Air medium";
parameter Integer nUni(final min=1, start=1)
"Number of units operating at design conditions";
final parameter Modelica.Units.SI.HeatFlowRate capUni_nominal(
final min=0)=dat.hea.Q_flow
"Heat pump design capacity (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mHeaWatUni_flow_nominal(
final min=0)=dat.hea.mLoa_flow
"HW design mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mAirUni_flow_nominal(
final min=0)=dat.hea.mSou_flow
"Air design mass flow rate (each unit)";
final parameter Modelica.Units.SI.MassFlowRate mHeaWat_flow_nominal(
final min=0)=nUni * mHeaWatUni_flow_nominal
"HW design mass flow rate (all units)";
parameter Modelica.Units.SI.PressureDifference dpPumHeaWatUni_nominal(
final min=0,
displayUnit="Pa")=dat.dpHeaLoa_nominal
"HW pump design head (each unit)";
replaceable parameter Buildings.Fluid.HeatPumps.Data.EquationFitReversible.Generic
dat
"Heat pump parameters (each unit)";
replaceable parameter Fluid.Movers.Data.Generic datPum
constrainedby Fluid.Movers.Data.Generic(
pressure(V_flow={0,1,2}*mHeaWatUni_flow_nominal/rho_nominal,
dp={1.14,1,0.42}*dpPumHeaWatUni_nominal),
motorCooledByFluid=false)
"Condenser pump parameters (each unit)";
// Assumptions
parameter Modelica.Fluid.Types.Dynamics energyDynamics=
Modelica.Fluid.Types.Dynamics.DynamicFreeInitial
"Type of energy balance: dynamic (3 initialization options) or steady state";
parameter Modelica.Units.SI.Time tau=30
"Time constant of fluid volume for nominal HW flow, used if energy or mass balance is dynamic";
// Pump speed filter parameters
parameter Boolean use_inputFilter=energyDynamics<>Modelica.Fluid.Types.Dynamics.SteadyState
"= true, if signal is filtered with a 2nd order CriticalDamping filter";
parameter Modelica.Units.SI.Time riseTime=30
"Rise time of the filter (time to reach 99.6 % of the speed)";
parameter Modelica.Blocks.Types.Init init=Modelica.Blocks.Types.Init.InitialOutput
"Type of initialization (no init/steady state/initial state/initial output)";
parameter Real y_start=1 "Initial position of actuator";
Buildings.Controls.OBC.CDL.Interfaces.BooleanInput y1[nUni]
"Heat pump On/Off command";
Buildings.Controls.OBC.CDL.Interfaces.RealInput TSet(
final unit="K", displayUnit="degC")
"Supply temperature setpoint";
BoundaryConditions.WeatherData.Bus weaBus
"Bus with weather data";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput PPum(
final unit="W")
"Power drawn by HW pumps";
Buildings.Controls.OBC.CDL.Interfaces.RealOutput P(
final unit="W")
"Power drawn by heat pumps";
Fluid.BaseClasses.MassFlowRateMultiplier mulConInl(
redeclare final package Medium = Medium,
final allowFlowReversal=allowFlowReversal,
final use_input=true)
"Flow rate multiplier";
Fluid.BaseClasses.MassFlowRateMultiplier mulConOut(
redeclare final package Medium = Medium,
final allowFlowReversal=allowFlowReversal,
final use_input=true)
"Flow rate multiplier";
BaseClasses.MultipleCommands com(final nUni=nUni)
"Convert command signals";
Buildings.Controls.OBC.CDL.Reals.Multiply mulP "Scale power";
Fluid.HeatPumps.EquationFitReversible heaPum(
redeclare final package Medium1=Medium,
redeclare final package Medium2=MediumAir,
final per=dat,
final tau1=tau,
final show_T=show_T,
final energyDynamics=energyDynamics,
final allowFlowReversal1=allowFlowReversal,
final allowFlowReversal2=false)
"Heat pump";
Fluid.Movers.SpeedControlled_y pum(
redeclare final package Medium = Medium,
inputType=Buildings.Fluid.Types.InputType.Continuous,
final tau=tau,
final show_T=show_T,
final allowFlowReversal=allowFlowReversal,
final energyDynamics=energyDynamics,
final use_inputFilter=use_inputFilter,
final riseTime=riseTime,
final init=init,
final per=datPum,
addPowerToMedium=false) "HW pump";
Fluid.Sources.MassFlowSource_WeatherData airSou(
redeclare final package Medium=MediumAir,
final use_m_flow_in=true,
final nPorts=1)
"Air flow source";
Fluid.Sources.Boundary_pT airSin(
redeclare final package Medium=MediumAir,
final nPorts=1)
"Air flow sink";
Buildings.Controls.OBC.CDL.Reals.MultiplyByParameter comFan(
final k=mAirUni_flow_nominal)
"Convert On/Off command to air flow setpoint";
Buildings.Controls.OBC.CDL.Reals.Multiply mulP1 "Scale power";
Buildings.Controls.OBC.CDL.Logical.Pre preY1[nUni]
"Left limit of signal avoiding direct feedback of status to controller";
Buildings.Controls.OBC.CDL.Conversions.BooleanToReal booToRea
"Convert to real";
Buildings.Controls.OBC.CDL.Reals.Multiply inp
"Compute pump input signal";
Buildings.Controls.OBC.CDL.Reals.Sources.Constant cst(final k=1)
"Constant";
Buildings.Controls.OBC.CDL.Conversions.RealToInteger reaToInt
"Convert to integer";
protected
parameter Medium.ThermodynamicState sta_nominal=Medium.setState_pTX(
T=dat.hea.TRefLoa,
p=3E5,
X=Medium.X_default)
"State of the medium at the medium default properties";
parameter Modelica.Units.SI.Density rho_nominal=Medium.density(sta_nominal)
"Density at the medium default properties";
equation
connect(mulConOut.uInv, mulConInl.u);
connect(com.nUniOn, mulP.u2);
connect(mulP.y, P);
connect(com.nUniOnBou, mulConOut.u);
connect(pum.port_b, heaPum.port_a1);
connect(heaPum.port_b1, mulConOut.port_a);
connect(mulConInl.port_b, pum.port_a);
connect(mulConOut.port_b, port_b);
connect(heaPum.P, mulP.u1);
connect(port_a, mulConInl.port_a);
connect(airSou.ports[1], heaPum.port_a2);
connect(airSin.ports[1], heaPum.port_b2);
connect(TSet, heaPum.TSet);
connect(comFan.y, airSou.m_flow_in);
connect(weaBus, airSou.weaBus);
connect(mulP1.y, PPum);
connect(pum.P, mulP1.u1);
connect(com.nUniOn, mulP1.u2);
connect(y1, preY1.u);
connect(preY1.y, com.y1);
connect(com.y1One, booToRea.u);
connect(inp.y, pum.y);
connect(booToRea.y, comFan.u);
connect(booToRea.y, inp.u2);
connect(cst.y, inp.u1);
connect(booToRea.y, reaToInt.u);
connect(reaToInt.y, heaPum.uMod);
end HeatPumpGroup;
Buildings.Experimental.DHC.Plants.Combined.Subsystems.MultiplePumpsDp
Model of multiple identical pumps in parallel with dp-controlled pump model
Information
This model represents a set of identical dp-controlled variable speed pumps that are piped in parallel. An optional check valve in series with each pump is included.
Control points
The following input and output points are available.
-
Start command (VFD Run)
y1: DO signal dedicated to each unit, with a dimensionality of one -
(Optionally if
have_varSetistrue) Differential pressure setpointdp_in: AO signal common to all units, with a dimensionality of zero -
Pump status
y1_actual: DI signal dedicated to each unit, with a dimensionality of one
Details
See the base class Buildings.Experimental.DHC.Plants.Combined.Subsystems.BaseClasses.PartialMultiplePumps. for a description of the modeling approach.
Extends from BaseClasses.PartialMultiplePumps (Base class for modeling multiple identical pumps in parallel).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | PartialMedium | Medium in the component | |
| Integer | nPum | Number of pumps | |
| Boolean | have_var | have_varSet | Set to true for variable speed pumps, false for constant speed |
| Boolean | have_valve | true | Set to true for inline check valve |
| Generic | per | redeclare parameter Fluid.Mo... | Pump parameters |
| Boolean | have_varSet | true | Set to true for variable setpoint, false for constant setpoint (design value) |
| Nominal condition | |||
| MassFlowRate | mPum_flow_nominal | Design mass flow rate (each pump) [kg/s] | |
| PressureDifference | dpPum_nominal | Design head (each pump) [Pa] | |
| PressureDifference | dpValve_nominal | 10000 | Pressure drop of check valve fully open [Pa] |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Nominal condition | |||
| Time | tau | 1 | Time constant of fluid volume for nominal flow, used if energy or mass balance is dynamic [s] |
| Filtered speed | |||
| Boolean | use_inputFilter | true | = true, if speed is filtered with a 2nd order CriticalDamping filter |
| Time | riseTime | 30 | Rise time of the filter (time to reach 99.6 % of the speed) [s] |
| Init | init | Modelica.Blocks.Types.Init.I... | Type of initialization (no init/steady state/initial state/initial output) |
| Initialization | |||
| AbsolutePressure | p_start | Medium.p_default | Start value of pressure [Pa] |
| Temperature | T_start | Medium.T_default | Start value of temperature [K] |
| MassFraction | X_start[Medium.nX] | Medium.X_default | Start value of mass fractions m_i/m [kg/kg] |
| ExtraProperty | C_start[Medium.nC] | fill(0, Medium.nC) | Start value of trace substances |
| ExtraProperty | C_nominal[Medium.nC] | fill(1E-2, Medium.nC) | Nominal value of trace substances. (Set to typical order of magnitude.) |
| 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) |
| input BooleanInput | y1[nPum] | Start signal (VFD Run or motor starter contact) |
| output RealOutput | P | Total power (all pumps) [W] |
| output BooleanOutput | y1_actual[nPum] | Pump status |
| input RealInput | dp_in | Differential pressure setpoint [Pa] |
Modelica definition
model MultiplePumpsDp
"Model of multiple identical pumps in parallel with dp-controlled pump model"
extends BaseClasses.PartialMultiplePumps(
final have_var=have_varSet,
redeclare final Buildings.Fluid.Movers.FlowControlled_dp pum(
final m_flow_nominal=mPum_flow_nominal,
final dp_nominal=dpPum_nominal),
cst(final k=dpPum_nominal));
parameter Boolean have_varSet = true
"Set to true for variable setpoint, false for constant setpoint (design value)";
Buildings.Controls.OBC.CDL.Interfaces.RealInput dp_in(
final unit="Pa",
final min=0)
if have_var
"Differential pressure setpoint";
equation
connect(dp_in, inp.u1);
connect(inp.y, pum.dp_in);
end MultiplePumpsDp;
Buildings.Experimental.DHC.Plants.Combined.Subsystems.MultiplePumpsFlow
Model of multiple identical pumps in parallel with flow-controlled pump model
Information
This model represents a set of identical flow-controlled variable speed pumps that are piped in parallel. An optional check valve in series with each pump is included.
Control points
The following input and output points are available.
-
Start command (VFD Run)
y1: DO signal dedicated to each unit, with a dimensionality of one -
(Optionally if
have_varSetistrue) Mass flow rate setpointm_flow_in: AO signal with a dimensionality of zero The setpoint corresponds to the total flow rate of all units. -
Pump status
y1_actual: DI signal dedicated to each unit, with a dimensionality of one
Details
See the base class Buildings.Experimental.DHC.Plants.Combined.Subsystems.BaseClasses.PartialMultiplePumps. for a description of the modeling approach.
Extends from BaseClasses.PartialMultiplePumps (Base class for modeling multiple identical pumps in parallel).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | PartialMedium | Medium in the component | |
| Integer | nPum | Number of pumps | |
| Boolean | have_var | have_varSet | Set to true for variable speed pumps, false for constant speed |
| Boolean | have_valve | true | Set to true for inline check valve |
| Generic | per | redeclare parameter Fluid.Mo... | Pump parameters |
| Boolean | have_varSet | true | Set to true for variable setpoint, false for constant setpoint (design value) |
| Nominal condition | |||
| MassFlowRate | mPum_flow_nominal | Design mass flow rate (each pump) [kg/s] | |
| PressureDifference | dpPum_nominal | Design head (each pump) [Pa] | |
| PressureDifference | dpValve_nominal | 10000 | Pressure drop of check valve fully open [Pa] |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Nominal condition | |||
| Time | tau | 1 | Time constant of fluid volume for nominal flow, used if energy or mass balance is dynamic [s] |
| Filtered speed | |||
| Boolean | use_inputFilter | true | = true, if speed is filtered with a 2nd order CriticalDamping filter |
| Time | riseTime | 30 | Rise time of the filter (time to reach 99.6 % of the speed) [s] |
| Init | init | Modelica.Blocks.Types.Init.I... | Type of initialization (no init/steady state/initial state/initial output) |
| Initialization | |||
| AbsolutePressure | p_start | Medium.p_default | Start value of pressure [Pa] |
| Temperature | T_start | Medium.T_default | Start value of temperature [K] |
| MassFraction | X_start[Medium.nX] | Medium.X_default | Start value of mass fractions m_i/m [kg/kg] |
| ExtraProperty | C_start[Medium.nC] | fill(0, Medium.nC) | Start value of trace substances |
| ExtraProperty | C_nominal[Medium.nC] | fill(1E-2, Medium.nC) | Nominal value of trace substances. (Set to typical order of magnitude.) |
| 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) |
| input BooleanInput | y1[nPum] | Start signal (VFD Run or motor starter contact) |
| output RealOutput | P | Total power (all pumps) [W] |
| output BooleanOutput | y1_actual[nPum] | Pump status |
| input RealInput | m_flow_in | Mass flow rate setpoint (total over all pumps) [kg/s] |
Modelica definition
model MultiplePumpsFlow
"Model of multiple identical pumps in parallel with flow-controlled pump model"
extends BaseClasses.PartialMultiplePumps(
final have_var=have_varSet,
redeclare final Buildings.Fluid.Movers.FlowControlled_m_flow pum(
final m_flow_nominal=mPum_flow_nominal,
final dp_nominal=dpPum_nominal),
cst(final k=mPum_flow_nominal));
parameter Boolean have_varSet = true
"Set to true for variable setpoint, false for constant setpoint (design value)";
Buildings.Controls.OBC.CDL.Interfaces.RealInput m_flow_in(
final unit="kg/s",
final min=0)
if have_var
"Mass flow rate setpoint (total over all pumps)";
Buildings.Controls.OBC.CDL.Reals.Divide div1
"Divide total flow setpoint by number of pumps commanded On";
equation
connect(m_flow_in, div1.u1);
connect(com.nUniOnBou, div1.u2);
connect(div1.y, inp.u1);
connect(inp.y, pum.m_flow_in);
end MultiplePumpsFlow;
Buildings.Experimental.DHC.Plants.Combined.Subsystems.MultiplePumpsSpeed
Model of multiple identical pumps in parallel with speed-controlled pump model
Information
This model represents a set of identical speed-controlled pumps that are piped in parallel. The model may be configured to represent either constant speed pumps or variable speed pumps. An optional check valve in series with each pump is included.
Control points
The following input and output points are available.
-
Start command (VFD Run for variable speed pumps or Starter contact for constant speed pumps)
y1: DO signal dedicated to each unit, with a dimensionality of one -
(Optionally if
have_varistrue) Speed commandy: AO signal common to all units, with a dimensionality of zero -
Pump status
y1_actual: DI signal dedicated to each unit, with a dimensionality of one
Details
See the base class Buildings.Experimental.DHC.Plants.Combined.Subsystems.BaseClasses.PartialMultiplePumps. for a description of the modeling approach.
Extends from Buildings.Experimental.DHC.Plants.Combined.Subsystems.BaseClasses.PartialMultiplePumps (Base class for modeling multiple identical pumps in parallel).
Parameters
| Type | Name | Default | Description |
|---|---|---|---|
| replaceable package Medium | PartialMedium | Medium in the component | |
| Integer | nPum | Number of pumps | |
| Boolean | have_var | true | Set to true for variable speed pumps, false for constant speed |
| Boolean | have_valve | true | Set to true for inline check valve |
| Generic | per | redeclare parameter Fluid.Mo... | Pump parameters |
| Nominal condition | |||
| MassFlowRate | mPum_flow_nominal | Design mass flow rate (each pump) [kg/s] | |
| PressureDifference | dpPum_nominal | Design head (each pump) [Pa] | |
| PressureDifference | dpValve_nominal | 10000 | Pressure drop of check valve fully open [Pa] |
| Dynamics | |||
| Conservation equations | |||
| Dynamics | energyDynamics | Modelica.Fluid.Types.Dynamic... | Type of energy balance: dynamic (3 initialization options) or steady state |
| Nominal condition | |||
| Time | tau | 1 | Time constant of fluid volume for nominal flow, used if energy or mass balance is dynamic [s] |
| Filtered speed | |||
| Boolean | use_inputFilter | true | = true, if speed is filtered with a 2nd order CriticalDamping filter |
| Time | riseTime | 30 | Rise time of the filter (time to reach 99.6 % of the speed) [s] |
| Init | init | Modelica.Blocks.Types.Init.I... | Type of initialization (no init/steady state/initial state/initial output) |
| Initialization | |||
| AbsolutePressure | p_start | Medium.p_default | Start value of pressure [Pa] |
| Temperature | T_start | Medium.T_default | Start value of temperature [K] |
| MassFraction | X_start[Medium.nX] | Medium.X_default | Start value of mass fractions m_i/m [kg/kg] |
| ExtraProperty | C_start[Medium.nC] | fill(0, Medium.nC) | Start value of trace substances |
| ExtraProperty | C_nominal[Medium.nC] | fill(1E-2, Medium.nC) | Nominal value of trace substances. (Set to typical order of magnitude.) |
| 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) |
| input BooleanInput | y1[nPum] | Start signal (VFD Run or motor starter contact) |
| output RealOutput | P | Total power (all pumps) [W] |
| output BooleanOutput | y1_actual[nPum] | Pump status |
| input RealInput | y | Pump speed (common to all pumps) [1] |
Modelica definition
model MultiplePumpsSpeed
"Model of multiple identical pumps in parallel with speed-controlled pump model"
extends Buildings.Experimental.DHC.Plants.Combined.Subsystems.BaseClasses.PartialMultiplePumps
(
redeclare final Buildings.Fluid.Movers.SpeedControlled_y pum,
cst(final k=1));
Buildings.Controls.OBC.CDL.Interfaces.RealInput y(
final unit="1",
final min=0,
final max=1)
if have_var
"Pump speed (common to all pumps)";
equation
connect(y, inp.u1);
connect(inp.y, pum.y);
end MultiplePumpsSpeed;