Library logo

Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression

Package with substations that have a vapor compression engine for heating or cooling

Information

This package contains models for district heating and cooling substations that provide heating or cooling using vapor compression engines.

Extends from Modelica.Icons.Package (Icon for standard packages).

Package Content

Name Description
Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.HeatingCoolingHotwaterTimeSeries_dT HeatingCoolingHotwaterTimeSeries_dT Substation for heating, cooling and domestic hot water with load as a time series
Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.Validation Validation Package with models for validation
Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.BaseClasses BaseClasses  

Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.HeatingCoolingHotwaterTimeSeries_dT Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.HeatingCoolingHotwaterTimeSeries_dT

Substation for heating, cooling and domestic hot water with load as a time series

Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.HeatingCoolingHotwaterTimeSeries_dT

Information

Model for a substation with space heating, space cooling and domestic hot water.

The model takes as parameters the temperature lift on the primary side and and then draws the required amount of water. The load is specified by a file that contains time series for the load profiles. This file needs to have the following format: 

#1
#Heating and Cooling Model loads for a SF large office
#First column: Seconds in the year (loads are hourly)
#Second column: cooling loads in Watts (as negative numbers).
#Third column: space heating loads in Watts. Heating is a combination of electric space heating, gas space heating
#Gas heaters in the model were 80% efficient where electric heaters were 100% efficient. Here the total watts = electric watts + .80*(gas watts)
#Fourth column: water heating = 0.8 * gas water watts
#Peak space cooling load = -383165.6989 Watts
#Peak space heating load = 893931.4335 Watts
#Peak water heating load = 19496.90012 Watts
double tab1(8760,4)
0,0,5972.314925,16
3600,0,4925.839944,1750.915684
7200,0,7470.393385,1750.971979
[to be continued]

Values at intermediate times are interpolated using cubic Hermite splines.

Implementation

The time series data are interpolated using Fritsch-Butland interpolation. This uses cubic Hermite splines such that y preserves the monotonicity and der(y) is continuous, also if extrapolated.

There is a control volume at each of the two fluid ports that are exposed by this model. These approximate the dynamics of the substation, and they also generally avoid nonlinear system of equations if multiple substations are connected to each other.

Parameters

TypeNameDefaultDescription
replaceable package MediumModelica.Media.Interfaces.Pa...Medium model for water
RealgaiCoo1Gain to scale cooling load
RealgaiHeagaiCooGain to scale heating load
RealgaiHotWatgaiHeaGain to scale hot water load
TemperatureTColMin273.15 + 8Minimum temperature of district cold water supply [K]
TemperatureTHotMax273.15 + 18Maximum temperature of district hot water supply [K]
TemperatureDifferencedTCooEva_nominal-4Temperature difference evaporator of the chiller [K]
StringfilNam Name of data file with heating and cooling load
TemperatureDifferencedTHotWatCon_nominal60 - 40Temperature difference condenser of hot water heat pump [K]
Design parameter
TemperatureDifferencedTCooCon_nominal4Temperature difference condenser of the chiller (positive) [K]
TemperatureDifferencedTHeaEva_nominal-4Temperature difference evaporator of the heat pump for space heating (negative) [K]
HeatFlowRateQCoo_flow_nominalgaiCoo*Buildings.Experimenta...Design heat flow rate [W]
HeatFlowRateQHea_flow_nominalgaiHea*Buildings.Experimenta...Design heat flow rate [W]
HeatFlowRateQHotWat_flow_nominalgaiHotWat*Buildings.Experime...Design heat flow rate for domestic hot water [W]
Pressuredp_nominal30000Pressure difference at nominal flow rate (for each flow leg) [Pa]
Nominal conditions
TemperatureTChiSup_nominal273.15 + 16Chilled water leaving temperature at the evaporator [K]
TemperatureTHeaSup_nominal273.15 + 30Supply temperature space heating system at TOut_nominal [K]
TemperatureTHeaRet_nominal273.15 + 25Return temperature space heating system at TOut_nominal [K]
TemperatureTOut_nominal Outside design temperature for heating [K]
Advanced
Diagnostics
Booleanshow_Tfalse= true, if actual temperature at port is computed
Dynamics
DynamicsmixingVolumeEnergyDynamicsModelica.Fluid.Types.Dynamic...Formulation of energy balance for mixing volume at inlet and outlet

Connectors

TypeNameDescription
replaceable package MediumMedium model for water
output RealOutputPHeaElectrical power consumed for space heating [W]
output RealOutputPHotWatElectrical power consumed for hot water [W]
output RealOutputPCooElectrical power consumed for space cooling [W]
output RealOutputQHea_flowSpace heating thermal load [W]
output RealOutputQHotWat_flowHot water thermal load [W]
output RealOutputQCoo_flowSpace cooling thermal load [W]
BusweaBus 
FluidPort_aport_aFluid connector a
FluidPort_bport_bFluid connector b
HeatPort_aheatPort_aHeat port for sensible heat input into volume a
HeatPort_aheatPort_bHeat port for sensible heat input into volume b

Modelica definition

model HeatingCoolingHotwaterTimeSeries_dT "Substation for heating, cooling and domestic hot water with load as a time series" replaceable package Medium = Modelica.Media.Interfaces.PartialMedium "Medium model for water"; parameter Real gaiCoo(min=0) = 1 "Gain to scale cooling load"; parameter Real gaiHea(min=0) = gaiCoo "Gain to scale heating load"; parameter Real gaiHotWat(min=0) = gaiHea "Gain to scale hot water load"; parameter Modelica.SIunits.Temperature TColMin = 273.15+8 "Minimum temperature of district cold water supply"; parameter Modelica.SIunits.Temperature THotMax = 273.15+18 "Maximum temperature of district hot water supply"; parameter Modelica.SIunits.TemperatureDifference dTCooCon_nominal( min=0.5, displayUnit="K") = 4 "Temperature difference condenser of the chiller (positive)"; parameter Modelica.SIunits.TemperatureDifference dTHeaEva_nominal( max=-0.5, displayUnit="K") = -4 "Temperature difference evaporator of the heat pump for space heating (negative)"; parameter Modelica.SIunits.TemperatureDifference dTCooEva_nominal=-4 "Temperature difference evaporator of the chiller"; parameter String filNam "Name of data file with heating and cooling load"; parameter Modelica.SIunits.HeatFlowRate QCoo_flow_nominal(max=-Modelica.Constants.eps)= gaiCoo * Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.BaseClasses.getPeakLoad( string="#Peak space cooling load", filNam=filNam) "Design heat flow rate"; parameter Modelica.SIunits.HeatFlowRate QHea_flow_nominal(min=Modelica.Constants.eps)= gaiHea * Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.BaseClasses.getPeakLoad( string="#Peak space heating load", filNam=filNam) "Design heat flow rate"; parameter Modelica.SIunits.HeatFlowRate QHotWat_flow_nominal(min=Modelica.Constants.eps)= gaiHotWat * Buildings.Experimental.DistrictHeatingCooling.SubStations.VaporCompression.BaseClasses.getPeakLoad( string="#Peak water heating load", filNam=filNam) "Design heat flow rate for domestic hot water"; parameter Modelica.SIunits.Temperature TChiSup_nominal = 273.15 + 16 "Chilled water leaving temperature at the evaporator"; parameter Modelica.SIunits.Temperature THeaSup_nominal = 273.15+30 "Supply temperature space heating system at TOut_nominal"; parameter Modelica.SIunits.Temperature THeaRet_nominal = 273.15+25 "Return temperature space heating system at TOut_nominal"; parameter Modelica.SIunits.Temperature TOut_nominal "Outside design temperature for heating"; parameter Modelica.SIunits.TemperatureDifference dTHotWatCon_nominal(min=0)=60-40 "Temperature difference condenser of hot water heat pump"; parameter Modelica.SIunits.Pressure dp_nominal(displayUnit="Pa")=30000 "Pressure difference at nominal flow rate (for each flow leg)"; final parameter Modelica.SIunits.MassFlowRate mCooCon_flow_nominal(min=0)= -QCoo_flow_nominal/cp_default/dTCooCon_nominal "Design mass flow rate for cooling at district side"; final parameter Modelica.SIunits.MassFlowRate mHeaEva_flow_nominal(min=0)= -QHea_flow_nominal/cp_default/dTHeaEva_nominal "Design mass flow rate for space heating at district side"; final parameter Modelica.SIunits.MassFlowRate mHotWatEva_flow_nominal(min=0)= QHotWat_flow_nominal/cp_default/dTHotWatCon_nominal "Design mass flow rate for domestic hot water at district side"; // Diagnostics parameter Boolean show_T = false "= true, if actual temperature at port is computed"; parameter Modelica.Fluid.Types.Dynamics mixingVolumeEnergyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial "Formulation of energy balance for mixing volume at inlet and outlet"; Modelica.Blocks.Interfaces.RealOutput PHea(unit="W") "Electrical power consumed for space heating"; Modelica.Blocks.Interfaces.RealOutput PHotWat(unit="W") "Electrical power consumed for hot water"; Modelica.Blocks.Interfaces.RealOutput PCoo(unit="W") "Electrical power consumed for space cooling"; Modelica.Blocks.Interfaces.RealOutput QHea_flow(unit="W") "Space heating thermal load"; Modelica.Blocks.Interfaces.RealOutput QHotWat_flow(unit="W") "Hot water thermal load"; Modelica.Blocks.Interfaces.RealOutput QCoo_flow(unit="W") "Space cooling thermal load"; BoundaryConditions.WeatherData.Bus weaBus; Modelica.Fluid.Interfaces.FluidPort_a port_a( redeclare final package Medium = Medium, h_outflow(start=Medium.h_default)) "Fluid connector a"; Modelica.Fluid.Interfaces.FluidPort_b port_b( redeclare final package Medium = Medium, h_outflow(start=Medium.h_default)) "Fluid connector b"; Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a heatPort_a "Heat port for sensible heat input into volume a"; Modelica.Thermal.HeatTransfer.Interfaces.HeatPort_a heatPort_b "Heat port for sensible heat input into volume b"; Medium.ThermodynamicState staHea_a2= Medium.setState_phX(heaPum.port_a2.p, noEvent(actualStream(heaPum.port_a2.h_outflow)), noEvent(actualStream(heaPum.port_a2.Xi_outflow))) if show_T "Medium properties in port_a2 of space heating heat pump intake"; Medium.ThermodynamicState staHea_b2= Medium.setState_phX(heaPum.port_b2.p, noEvent(actualStream(heaPum.port_b2.h_outflow)), noEvent(actualStream(heaPum.port_b2.Xi_outflow))) if show_T "Medium properties in port_b2 of space heating heat pump outlet"; Medium.ThermodynamicState staHotWat_a2= Medium.setState_phX(heaPumHotWat.port_a2.p, noEvent(actualStream(heaPumHotWat.port_a2.h_outflow)), noEvent(actualStream(heaPumHotWat.port_a2.Xi_outflow))) if show_T "Medium properties in port_a2 of hot water heat pump intake"; Medium.ThermodynamicState staHotWat_b2= Medium.setState_phX(heaPumHotWat.port_b2.p, noEvent(actualStream(heaPumHotWat.port_b2.h_outflow)), noEvent(actualStream(heaPumHotWat.port_b2.Xi_outflow))) if show_T "Medium properties in port_b2 of hot water heat pump outlet"; Medium.ThermodynamicState staCoo_a1= Medium.setState_phX(chi.port_a1.p, noEvent(actualStream(chi.port_a1.h_outflow)), noEvent(actualStream(chi.port_a1.Xi_outflow))) if show_T "Medium properties in port_a1 of chiller intake"; Medium.ThermodynamicState staCoo_b1= Medium.setState_phX(chi.port_b1.p, noEvent(actualStream(chi.port_b1.h_outflow)), noEvent(actualStream(chi.port_b1.Xi_outflow))) if show_T "Medium properties in port_b1 of chiller outlet"; constant Modelica.SIunits.SpecificHeatCapacity cp_default=4184 "Specific heat capacity of the fluid"; Buildings.Fluid.HeatPumps.Carnot_TCon heaPum( redeclare package Medium1 = Medium, redeclare package Medium2 = Medium, dTEva_nominal=dTHeaEva_nominal, dTCon_nominal=THeaSup_nominal-THeaRet_nominal, allowFlowReversal1=false, allowFlowReversal2=allowFlowReversal, use_eta_Carnot_nominal=true, etaCarnot_nominal=0.3, QCon_flow_nominal=QHea_flow_nominal, dp1_nominal=dp_nominal, dp2_nominal=dp_nominal) "Heat pump"; Buildings.Fluid.Movers.FlowControlled_m_flow pumHea( redeclare package Medium = Medium, energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState, allowFlowReversal=allowFlowReversal, m_flow_nominal=mHeaEva_flow_nominal, inputType=Buildings.Fluid.Types.InputType.Continuous, use_inputFilter=false, addPowerToMedium=false) "Pump for space heating heat pump"; Buildings.Controls.SetPoints.HotWaterTemperatureReset TSupHeaSet( TSup_nominal=THeaSup_nominal, TRet_nominal=THeaRet_nominal, TOut_nominal=TOut_nominal) "Set points for heating supply temperature"; Buildings.Fluid.Movers.FlowControlled_m_flow pumHotWat( redeclare package Medium = Medium, energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState, m_flow_nominal=mHotWatEva_flow_nominal, allowFlowReversal=allowFlowReversal, inputType=Buildings.Fluid.Types.InputType.Continuous, use_inputFilter=false, addPowerToMedium=false) "Pump"; Buildings.Fluid.HeatPumps.Carnot_TCon heaPumHotWat( redeclare package Medium1 = Medium, redeclare package Medium2 = Medium, dTEva_nominal=dTHeaEva_nominal, allowFlowReversal1=false, allowFlowReversal2=allowFlowReversal, use_eta_Carnot_nominal=true, etaCarnot_nominal=0.3, dp1_nominal=dp_nominal, dp2_nominal=dp_nominal, QCon_flow_nominal=QHotWat_flow_nominal, dTCon_nominal=dTHotWatCon_nominal) "Heat pump"; Buildings.Fluid.Chillers.Carnot_TEva chi( redeclare package Medium1 = Medium, redeclare package Medium2 = Medium, use_eta_Carnot_nominal=true, etaCarnot_nominal=0.3, dp1_nominal=dp_nominal, dp2_nominal=dp_nominal, QEva_flow_nominal=QCoo_flow_nominal, dTEva_nominal=dTCooEva_nominal, dTCon_nominal=dTCooCon_nominal, allowFlowReversal1=allowFlowReversal, allowFlowReversal2=false) "Chiller"; Buildings.Fluid.Movers.FlowControlled_m_flow pumChi( redeclare package Medium = Medium, energyDynamics=Modelica.Fluid.Types.Dynamics.SteadyState, allowFlowReversal=allowFlowReversal, inputType=Buildings.Fluid.Types.InputType.Continuous, use_inputFilter=false, m_flow_nominal=mCooCon_flow_nominal, addPowerToMedium=false) "Pump"; protected constant Boolean allowFlowReversal = false "= true to allow flow reversal, false restricts to design direction (port_a -> port_b)"; 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_check= Medium.specificHeatCapacityCp(sta_default) "Specific heat capacity of the fluid"; Modelica.Blocks.Sources.CombiTimeTable loa( tableOnFile=true, tableName="tab1", fileName=Modelica.Utilities.Files.loadResource(filNam), extrapolation=Modelica.Blocks.Types.Extrapolation.Periodic, y(each unit="W"), offset={0,0,0}, columns={2,3,4}, smoothness=Modelica.Blocks.Types.Smoothness.MonotoneContinuousDerivative1) "Loads"; Modelica.Blocks.Routing.DeMultiplex3 deMul "De multiplex"; Buildings.Fluid.Sources.MassFlowSource_T sou( redeclare package Medium = Medium, nPorts=1, use_m_flow_in=true, use_T_in=true) "Mass flow source"; Buildings.Fluid.Sources.FixedBoundary sin1( redeclare package Medium = Medium, nPorts=1); Modelica.Blocks.Math.Add dTHeaAct(k2=-1) "Temperature difference for heating"; Modelica.Blocks.Math.Division mConFlow "Mass flow rate for condenser"; Modelica.Blocks.Math.Gain gainH(k=cp_default) "Gain for heating"; Modelica.Blocks.Math.Add QEva_flow(k1=-1) "Heat flow rate at evaporator"; Modelica.Blocks.Math.Add PHeaAct "Power consumption for heating"; Modelica.Blocks.Math.Add QEvaHotWat_flow(k1=-1) "Heat flow rate at evaporator"; Modelica.Blocks.Math.Add PHotWatAct "Power consumption for heating"; Buildings.Fluid.Sources.FixedBoundary sinHotWat( redeclare package Medium = Medium, nPorts=1) "Pressure source"; Buildings.Fluid.Sources.MassFlowSource_T souHotWat( redeclare package Medium = Medium, nPorts=1, use_m_flow_in=true, use_T_in=false, T=313.15) "Mass flow source"; Modelica.Blocks.Math.Division mConHotWatFlow "Mass flow rate for condenser"; Modelica.Blocks.Sources.Constant gainHotWat(k=dTHotWatCon_nominal*cp_default) "Gain for hot water"; Buildings.Fluid.Sources.FixedBoundary sin2( redeclare package Medium = Medium, nPorts=1); Buildings.Fluid.Sources.MassFlowSource_T sou1( redeclare package Medium = Medium, nPorts=1, use_m_flow_in=true, use_T_in=true) "Mass flow source"; Modelica.Blocks.Math.Gain mPumCoo_flow(k=1/(cp_default*dTCooCon_nominal)) "Mass flow rate for cooling loop"; Modelica.Blocks.Math.Add QCon_flow(k1=-1) "Heat flow rate at condenser"; Modelica.Blocks.Math.Add PCooAct "Power consumption for cooling"; Modelica.Blocks.Sources.Constant TChiSup(k=TChiSup_nominal) "Supply water temperature set point for chilled water loop"; Buildings.Utilities.Math.SmoothMax smoothMax( deltaX=(THeaSup_nominal - THeaRet_nominal)*0.005); Modelica.Blocks.Sources.Constant dTHeaPumConMin(k=4) "Minimum temperature difference over condenser"; Modelica.Blocks.Sources.Constant TSetHotWat(k=273.15 + 60) "Set point for hot water temperature"; Modelica.Blocks.Math.Add THeaPumSup "Set point for heat pump outlet temperature. This block avoids TSup=TRet"; Modelica.Blocks.Sources.Constant dTChiEvaMin(k=-dTCooEva_nominal) "Temperature difference over evaporator"; Modelica.Blocks.Math.Add TEvaIn "Evaporator inlet temperature"; Buildings.Fluid.Delays.DelayFirstOrder del_a( redeclare final package Medium = Medium, nPorts=4, final m_flow_nominal=(mHeaEva_flow_nominal + mCooCon_flow_nominal + mHotWatEva_flow_nominal)/2, final allowFlowReversal=true, tau=600, final energyDynamics=mixingVolumeEnergyDynamics) "Mixing volume to break algebraic loops and to emulate the delay of the substation"; Buildings.Fluid.Delays.DelayFirstOrder del_b( redeclare final package Medium = Medium, nPorts=4, final m_flow_nominal=(mHeaEva_flow_nominal + mCooCon_flow_nominal + mHotWatEva_flow_nominal)/2, final allowFlowReversal=true, tau=600, final energyDynamics=mixingVolumeEnergyDynamics) "Mixing volume to break algebraic loops and to emulate the delay of the substation"; Modelica.Blocks.Math.Gain gaiLoa[3](final k={gaiCoo, gaiHea, gaiHotWat}) "Gain that can be used to scale the individual loads up or down. Components are cooling, heating and hot water"; Modelica.Blocks.Math.Gain mPumHotWat_flow( final k=1/(cp_default*dTHeaEva_nominal)) "Mass flow rate for hot water loop"; Modelica.Blocks.Math.Gain mPumHea_flow( final k=1/(cp_default*dTHeaEva_nominal)) "Mass flow rate for heating loop"; Modelica.Blocks.Math.Gain mPumCooEva_flow( final k=-1/(cp_default*dTCooCon_nominal)) "Mass flow rate for chiller evaporator water loop"; initial equation assert(abs((cp_default-cp_default_check)/cp_default) < 0.1, "Wrong cp_default value. Check cp_default constant."); assert(QCoo_flow_nominal < 0, "Nominal cooling rate must be strictly negative. Obtained QCoo_flow_nominal = " + String(QCoo_flow_nominal)); assert(QHea_flow_nominal > 0, "Nominal heating rate must be strictly positive. Obtained QHea_flow_nominal = " + String(QHea_flow_nominal)); assert(QHotWat_flow_nominal > 0, "Nominal hot water heating rate must be strictly positive. Obtained QHotWat_flow_nominal = " + String(QHotWat_flow_nominal)); equation connect(pumHea.port_b, heaPum.port_a2); connect(dTHeaAct.u1, TSupHeaSet.TSup); connect(TSupHeaSet.TRet, dTHeaAct.u2); connect(mConFlow.u1, deMul.y2[1]); connect(gainH.y, mConFlow.u2); connect(mConFlow.y, sou.m_flow_in); connect(sou.ports[1], heaPum.port_a1); connect(heaPum.port_b1, sin1.ports[1]); connect(sou.T_in, TSupHeaSet.TRet); connect(QEva_flow.u1, deMul.y2[1]); connect(QEva_flow.u2, heaPum.P); connect(PHeaAct.y, PHea); connect(PHeaAct.u1, pumHea.P); connect(PHeaAct.u2, heaPum.P); connect(QEvaHotWat_flow.u2, heaPumHotWat.P); connect(PHotWatAct.u2, heaPumHotWat.P); connect(PHotWatAct.u1, pumHotWat.P); connect(heaPumHotWat.port_b1, sinHotWat.ports[1]); connect(souHotWat.ports[1], heaPumHotWat.port_a1); connect(mConHotWatFlow.y, souHotWat.m_flow_in); connect(gainHotWat.y, mConHotWatFlow.u2); connect(QEvaHotWat_flow.u1, deMul.y3[1]); connect(PHotWatAct.y, PHotWat); connect(pumHotWat.port_b, heaPumHotWat.port_a2); connect(mConHotWatFlow.u1, deMul.y3[1]); connect(chi.port_b2, sin2.ports[1]); connect(chi.port_a2, sou1.ports[1]); connect(mPumCoo_flow.y, pumChi.m_flow_in); connect(mPumCoo_flow.u, QCon_flow.y); connect(pumChi.port_b, chi.port_a1); connect(QCon_flow.u1, deMul.y1[1]); connect(QCon_flow.u2, chi.P); connect(PCooAct.y, PCoo); connect(PCooAct.u1, pumChi.P); connect(PCooAct.u2, chi.P); connect(TChiSup.y, chi.TSet); connect(dTHeaAct.y, smoothMax.u1); connect(dTHeaPumConMin.y, smoothMax.u2); connect(TSetHotWat.y, heaPumHotWat.TSet); connect(gainH.u, smoothMax.y); connect(THeaPumSup.u1, TSupHeaSet.TRet); connect(THeaPumSup.u2, smoothMax.y); connect(THeaPumSup.y, heaPum.TSet); connect(TChiSup.y, TEvaIn.u1); connect(dTChiEvaMin.y, TEvaIn.u2); connect(TEvaIn.y, sou1.T_in); connect(del_a.ports[1], port_a); connect(pumHea.port_a, del_a.ports[2]); connect(pumHotWat.port_a, del_a.ports[3]); connect(chi.port_b1, del_a.ports[4]); connect(heaPum.port_b2, del_b.ports[1]); connect(heaPumHotWat.port_b2, del_b.ports[2]); connect(pumChi.port_a, del_b.ports[3]); connect(port_b, del_b.ports[4]); connect(weaBus.TDryBul, TSupHeaSet.TOut); connect(deMul.y2[1], QHea_flow); connect(deMul.y3[1], QHotWat_flow); connect(deMul.y1[1], QCoo_flow); connect(loa.y, gaiLoa.u); connect(gaiLoa.y, deMul.u); connect(del_a.heatPort, heatPort_a); connect(del_b.heatPort, heatPort_b); connect(QEvaHotWat_flow.y, mPumHotWat_flow.u); connect(mPumHotWat_flow.y, pumHotWat.m_flow_in); connect(QEva_flow.y, mPumHea_flow.u); connect(mPumHea_flow.y, pumHea.m_flow_in); connect(deMul.y1[1], mPumCooEva_flow.u); connect(mPumCooEva_flow.y, sou1.m_flow_in); end HeatingCoolingHotwaterTimeSeries_dT;