model Case600 "Case 600FF, but with dual-setpoint for heating and cooling" extends Case600FF( redeclare Buildings.ThermalZones.Detailed.Validation.BESTEST.Data.StandardResults staRes( annualHea(Min=4.296*3.6e9, Max=5.709*3.6e9, Mean=5.090*3.6e9), annualCoo(Min=-6.137*3.6e9, Max=-7.964*3.6e9, Mean=-6.832*3.6e9), peakHea(Min=3.437*1000, Max=4.354*1000, Mean=4.000*1000), peakCoo(Min=-5.965*1000, Max=-6.827*1000, Mean=-6.461*1000))); Buildings.Controls.OBC.CDL.Continuous.LimPID conHea( k=0.1, Ti=300, controllerType=Buildings.Controls.OBC.CDL.Types.SimpleController.PI) "Controller for heating"; Buildings.Controls.OBC.CDL.Continuous.LimPID conCoo( k=0.1, Ti=300, reverseAction=true, controllerType=Buildings.Controls.OBC.CDL.Types.SimpleController.PI) "Controller for cooling"; Buildings.Controls.OBC.CDL.Continuous.Gain gaiHea(k=1E6) "Gain for heating"; Buildings.Controls.OBC.CDL.Continuous.Gain gaiCoo(k=-1E6) "Gain for cooling"; Modelica.Blocks.Math.Sum sumHeaCoo(nin=2) "Sum of heating and cooling heat flow rate"; Modelica.Blocks.Routing.Multiplex2 multiplex2; Modelica.Thermal.HeatTransfer.Sources.PrescribedHeatFlow preHea "Prescribed heat flow for heating and cooling"; Modelica.Blocks.Continuous.Integrator EHea( k=1, initType=Modelica.Blocks.Types.Init.InitialState, y_start=0, u(unit="W"), y(unit="J")) "Heating energy in Joules"; Modelica.Blocks.Continuous.Integrator ECoo( k=1, initType=Modelica.Blocks.Types.Init.InitialState, y_start=0, u(unit="W"), y(unit="J")) "Cooling energy in Joules"; BaseClasses.DaySchedule TSetHea(table=[0.0,273.15 + 20]) "Heating setpoint"; BaseClasses.DaySchedule TSetCoo(table=[0.0,273.15 + 27]) "Cooling setpoint"; Buildings.Controls.OBC.CDL.Continuous.MovingMean PHea(delta=3600) "Hourly averaged heating power"; Buildings.Controls.OBC.CDL.Continuous.MovingMean PCoo(delta=3600) "Hourly averaged cooling power"; equation connect(TRooAir.T,conHea. u_m); connect(conCoo.u_m, TRooAir.T); connect(conHea.y,gaiHea. u); connect(conCoo.y,gaiCoo. u); connect(gaiHea.y,multiplex2. u1[1]); connect(gaiCoo.y,multiplex2. u2[1]); connect(multiplex2.y, sumHeaCoo.u); connect(sumHeaCoo.y, preHea.Q_flow); connect(EHea.u,gaiHea. y); connect(ECoo.u,gaiCoo. y); connect(TSetHea.y[1],conHea. u_s); connect(TSetCoo.y[1],conCoo. u_s); connect(PCoo.u,gaiCoo. y); connect(PHea.u,gaiHea. y); connect(preHea.port, roo.heaPorAir); end Case600;

model Case600FF "Basic test with light-weight construction and free floating temperature" extends Modelica.Icons.Example; package MediumA = Buildings.Media.Air "Medium model"; parameter Modelica.SIunits.Angle S_= Buildings.Types.Azimuth.S "Azimuth for south walls"; parameter Modelica.SIunits.Angle E_= Buildings.Types.Azimuth.E "Azimuth for east walls"; parameter Modelica.SIunits.Angle W_= Buildings.Types.Azimuth.W "Azimuth for west walls"; parameter Modelica.SIunits.Angle N_= Buildings.Types.Azimuth.N "Azimuth for north walls"; parameter Modelica.SIunits.Angle C_= Buildings.Types.Tilt.Ceiling "Tilt for ceiling"; parameter Modelica.SIunits.Angle F_= Buildings.Types.Tilt.Floor "Tilt for floor"; parameter Modelica.SIunits.Angle Z_= Buildings.Types.Tilt.Wall "Tilt for wall"; parameter Integer nConExtWin = 1 "Number of constructions with a window"; parameter Integer nConBou = 1 "Number of surface that are connected to constructions that are modeled inside the room"; parameter Buildings.HeatTransfer.Data.OpaqueConstructions.Generic matExtWal( nLay=3, absIR_a=0.9, absIR_b=0.9, absSol_a=0.6, absSol_b=0.6, material={Buildings.HeatTransfer.Data.Solids.Generic( x=0.009, k=0.140, c=900, d=530, nStaRef=nStaRef),Buildings.HeatTransfer.Data.Solids.Generic( x=0.066, k=0.040, c=840, d=12, nStaRef=nStaRef),Buildings.HeatTransfer.Data.Solids.Generic( x=0.012, k=0.160, c=840, d=950, nStaRef=nStaRef)}) "Exterior wall"; parameter Buildings.HeatTransfer.Data.OpaqueConstructions.Generic matFlo( final nLay=2, absIR_a=0.9, absIR_b=0.9, absSol_a=0.6, absSol_b=0.6, material={Buildings.HeatTransfer.Data.Solids.Generic( x=1.003, k=0.040, c=0, d=0, nStaRef=nStaRef),Buildings.HeatTransfer.Data.Solids.Generic( x=0.025, k=0.140, c=1200, d=650, nStaRef=nStaRef)}) "Floor"; parameter Buildings.HeatTransfer.Data.Solids.Generic soil( x=2, k=1.3, c=800, d=1500) "Soil properties"; Buildings.ThermalZones.Detailed.MixedAir roo( redeclare package Medium = MediumA, hRoo=2.7, nConExtWin=nConExtWin, nConBou=1, nPorts=3, energyDynamics=Modelica.Fluid.Types.Dynamics.FixedInitial, AFlo=48, datConBou( layers={matFlo}, each A=48, each til=F_), datConExt( layers={roof,matExtWal,matExtWal,matExtWal}, A={48,6*2.7,6*2.7,8*2.7}, til={C_,Z_,Z_,Z_}, azi={S_,W_,E_,N_}), nConExt=4, nConPar=0, nSurBou=0, datConExtWin( layers={matExtWal}, A={8*2.7}, glaSys={window600}, wWin={2*3}, hWin={2}, fFra={0.001}, til={Z_}, azi={S_}), lat=weaDat.lat, massDynamics=Modelica.Fluid.Types.Dynamics.SteadyState) "Room model"; Modelica.Blocks.Sources.Constant qConGai_flow(k=80/48) "Convective heat gain"; Modelica.Blocks.Sources.Constant qRadGai_flow(k=120/48) "Radiative heat gain"; Modelica.Blocks.Routing.Multiplex3 multiplex3_1 "Multiplex for internal gains"; Modelica.Blocks.Sources.Constant qLatGai_flow(k=0) "Latent heat gain"; Buildings.BoundaryConditions.WeatherData.ReaderTMY3 weaDat(filNam= Modelica.Utilities.Files.loadResource("modelica://Buildings/Resources/weatherdata/DRYCOLD.mos"), computeWetBulbTemperature=false); Modelica.Blocks.Sources.Constant uSha(k=0) "Control signal for the shading device"; Modelica.Blocks.Routing.Replicator replicator(nout=max(1,nConExtWin)); Modelica.Thermal.HeatTransfer.Sources.FixedTemperature TSoi[nConBou](each T= 283.15) "Boundary condition for construction"; parameter Buildings.HeatTransfer.Data.OpaqueConstructions.Generic roof( nLay=3, absIR_a=0.9, absIR_b=0.9, absSol_a=0.6, absSol_b=0.6, material={Buildings.HeatTransfer.Data.Solids.Generic( x=0.019, k=0.140, c=900, d=530, nStaRef=nStaRef),Buildings.HeatTransfer.Data.Solids.Generic( x=0.1118, k=0.040, c=840, d=12, nStaRef=nStaRef),Buildings.HeatTransfer.Data.Solids.Generic( x=0.010, k=0.160, c=840, d=950, nStaRef=nStaRef)}) "Roof"; parameter Buildings.ThermalZones.Detailed.Validation.BESTEST.Data.Win600 window600( UFra=3, haveExteriorShade=false, haveInteriorShade=false) "Window"; Buildings.HeatTransfer.Conduction.SingleLayer soi( A=48, material=soil, steadyStateInitial=true, stateAtSurface_a=false, stateAtSurface_b=true, T_a_start=283.15, T_b_start=283.75) "2m deep soil (per definition on p.4 of ASHRAE 140-2007)"; Buildings.Fluid.Sources.MassFlowSource_T sinInf( redeclare package Medium = MediumA, m_flow=1, use_m_flow_in=true, use_T_in=false, use_X_in=false, use_C_in=false, nPorts=1) "Sink model for air infiltration"; Buildings.Fluid.Sources.Outside souInf(redeclare package Medium = MediumA, nPorts=1) "Source model for air infiltration"; Modelica.Blocks.Sources.Constant InfiltrationRate(k=-48*2.7*0.5/3600) "0.41 ACH adjusted for the altitude (0.5 at sea level)"; Modelica.Blocks.Math.Product product "Product to compute infiltration mass flow rate"; Buildings.Fluid.Sensors.Density density( redeclare package Medium = MediumA) "Air density inside the building"; Buildings.BoundaryConditions.WeatherData.Bus weaBus "Weather data bus"; Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensor TRooAir "Room air temperature"; replaceable parameter Buildings.ThermalZones.Detailed.Validation.BESTEST.Data.StandardResultsFreeFloating staRes( minT( Min=-18.8+273.15, Max=-15.6+273.15, Mean=-17.6+273.15), maxT( Min=64.9+273.15, Max=69.5+273.15, Mean=66.2+273.15), meanT(Min=24.2+273.15, Max=25.9+273.15, Mean=25.1+273.15)) constrainedby Modelica.Icons.Record "Reference results from ASHRAE/ANSI Standard 140"; Modelica.Blocks.Math.MultiSum multiSum(nu=1) "Multi sum for infiltration air flow rate"; Controls.OBC.CDL.Continuous.MovingMean TRooHou(delta=3600) "Hourly averaged room air temperature"; Controls.OBC.CDL.Continuous.MovingMean TRooAnn(delta=86400*365) "Annual averaged room air temperature"; equation connect(qRadGai_flow.y,multiplex3_1. u1[1]); connect(qLatGai_flow.y,multiplex3_1. u3[1]); connect(multiplex3_1.y, roo.qGai_flow); connect(roo.uSha, replicator.y); connect(qConGai_flow.y, multiplex3_1.u2[1]); connect(weaDat.weaBus, roo.weaBus); connect(uSha.y, replicator.u); connect(product.y, sinInf.m_flow_in); connect(density.port, roo.ports[1]); connect(density.d, product.u2); connect(weaDat.weaBus, weaBus); connect(TSoi[1].port, soi.port_a); connect(soi.port_b, roo.surf_conBou[1]); connect(weaBus, souInf.weaBus); connect(roo.heaPorAir, TRooAir.port); connect(sinInf.ports[1], roo.ports[2]); connect(multiSum.y, product.u1); connect(InfiltrationRate.y, multiSum.u[1]); connect(TRooAir.T, TRooHou.u); connect(TRooAir.T, TRooAnn.u); connect(souInf.ports[1], roo.ports[3]); end Case600FF;

model Case610 "Case 600 with south shading" extends Buildings.ThermalZones.Detailed.Validation.BESTEST.Cases6xx.Case600( roo( datConExtWin( ove( wR={0.5}, wL={0.5}, dep={1}, gap={0.5}))), staRes( annualHea(Min=4.355*3.6e9, Max=5.786*3.6e9, Mean=5.146*3.6e9), annualCoo(Min=-3.915*3.6e9, Max=-5.778*3.6e9, Mean=-4.964*3.6e9), peakHea(Min=3.437*1000, Max=4.354*1000, Mean=3.998*1000), peakCoo(Min=-5.689*1000, Max=-6.371*1000, Mean=-5.988*1000))); end Case610;

model Case620 "Case 600, but with windows on East and West side walls" extends Buildings.ThermalZones.Detailed.Validation.BESTEST.Cases6xx.Case600( roo( nConExtWin=2, datConExtWin( layers={matExtWal, matExtWal}, each A=6*2.7, glaSys={window600, window600}, each wWin=3, each hWin=2, each fFra=0.001, each til=Z_, azi={W_,E_}), nConExt=3, datConExt( layers={roof,matExtWal,matExtWal}, A={48,8*2.7,8*2.7}, til={C_,Z_,Z_}, azi={S_,S_,N_})), staRes( annualHea(Min=4.613*3.6e9, Max=5.944*3.6e9, Mean=5.407*3.6e9), annualCoo(Min=-3.417*3.6e9, Max=-5.004*3.6e9, Mean=-4.218*3.6e9), peakHea(Min=3.591*1000, Max=4.379*1000, Mean=4.062*1000), peakCoo(Min=-3.634*1000, Max=-5.096*1000, Mean=-4.343*1000))); end Case620;

model Case630 "Case 620, but with added overhang and sidefins" extends Buildings.ThermalZones.Detailed.Validation.BESTEST.Cases6xx.Case620( roo( datConExtWin( ove( each wR=0.0, each wL=0.0, each dep=1.0, each gap=0.5), sidFin( each h=0.5, each dep=1.0, each gap=0.0))), staRes( annualHea(Min=5.050*3.6e9, Max=6.469*3.6e9, Mean=5.783*3.6e9), annualCoo(Min=-2.129*3.6e9, Max=-3.701*3.6e9, Mean=-2.832*3.6e9), peakHea(Min=3.592*1000, Max=4.280*1000, Mean=4.006*1000), peakCoo(Min=-3.072*1000, Max=-4.116*1000, Mean=-3.626*1000))); end Case630;

model Case640 "Case 600, but with heating schedule" extends Case600( TSetHea(table=[ 0, 273.15 + 10; 7*3600, 273.15 + 10; 7*3600, 273.15 + 20; 23*3600, 273.15 + 20; 23*3600,273.15 + 10; 24*3600,273.15 + 10]), staRes( annualHea(Min=2.751*3.6e9, Max=3.803*3.6e9, Mean=3.207*3.6e9), annualCoo(Min=-5.952*3.6e9, Max=-7.811*3.6e9, Mean=-6.592*3.6e9), peakHea(Min = 5.232*1000, Max = 6.954*1000, Mean = 6.031 * 1000), peakCoo(Min= -5.892*1000, Max = -6.776*1000, Mean= -6.410 * 1000))); end Case640;

model Case650 "Case 600, but cooling based on schedule, night venting, and no heating" extends Case600( TSetHea(table=[0, 273.15 -200]), TSetCoo(table=[ 0, 273.15+100; 7*3600, 273.15+100; 7*3600, 273.15+27; 18*3600, 273.15+27; 18*3600, 273.15+100; 24*3600, 273.15+100]), staRes( annualHea(Min=0*3.6e9, Max=0*3.6e9, Mean=0*3.6e9), annualCoo(Min=-4.816*3.6e9, Max=-6.545*3.6e9, Mean=-5.482*3.6e9), peakHea(Min = 0*1000, Max = 0*1000, Mean = 0 * 1000), peakCoo(Min= -5.831*1000, Max = -6.679*1000, Mean= -6.321*1000)), gaiHea(k=0), multiSum(nu=2)); BaseClasses.DaySchedule vent(table=[ 0, -1703.16/3600; 7*3600, -1703.16/3600; 7*3600, 0; 18*3600, 0; 18*3600, -1703.16/3600; 24*3600, -1703.16/3600]) "Ventilation air flow rate"; equation connect(multiSum.u[2], vent.y[1]); end Case650;

model Case650FF "Case 600, no heating, no cooling, and ventilation as in case 650" extends Case600FF( redeclare Buildings.ThermalZones.Detailed.Validation.BESTEST.Data.StandardResultsFreeFloating staRes( minT( Min=-23.0+273.15, Max=-21.6+273.15, Mean=-22.7+273.15), maxT( Min=63.2+273.15, Max=68.2+273.15, Mean=64.7+273.15), meanT(Min=18.0+273.15, Max=19.6+273.15, Mean=18.7+273.15)), multiSum(nu=2)); BaseClasses.DaySchedule vent(table=[ 0, -1703.16/3600; 7*3600, -1703.16/3600; 7*3600, 0; 18*3600, 0; 18*3600, -1703.16/3600; 24*3600, -1703.16/3600]) "Ventilation air flow rate"; equation connect(multiSum.u[2],vent.y[1]); end Case650FF;