RWTH-EBC · PGorzalka · Apr 20, 2023 · Apr 20, 2023 · Apr 21, 2023 · Apr 26, 2023
diff --git a/AixLib/BoundaryConditions/GroundTemperature/DataSourceGroTem.mo b/AixLib/BoundaryConditions/GroundTemperature/DataSourceGroTem.mo
@@ -0,0 +1,20 @@
+within AixLib.BoundaryConditions.GroundTemperature;
+type DataSourceGroTem = enumeration(
+    Constant "Use constant value",
+    File "Use data from file with tabular data",
+    Kusuda "Use Kusuda model for undisturbed ground temperature with selectable parameters",
+    Sine "Use a sine model with selectable parameters")
+  "Enumeration to define data source for the ground temperature"
+      annotation(Documentation(info="<html>
+ <p>
+ Enumeration to define the data source used in AixLib.BoundaryConditions.GroundTemperature.Options.
+ </p>
+ </html>",
+        revisions="<html>
+ <ul>
+ <li>
+ April 21, 2023, by Philip Groesdonk:<br/>
+ First implementation.
+ </li>
+ </ul>
+ </html>"));
diff --git a/AixLib/BoundaryConditions/GroundTemperature/Options.mo b/AixLib/BoundaryConditions/GroundTemperature/Options.mo
@@ -0,0 +1,157 @@
+within AixLib.BoundaryConditions.GroundTemperature;
+model Options
+  "A model that provides different options for ground temperature output"
+  import Modelica.Constants.pi;
+
+  parameter AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem datSou
+    =AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Constant
+    "Choice by which option for ground temperature is to be provided"
+    annotation (Dialog(group="General"));
+  parameter Modelica.Units.SI.Temperature TMea "constant or mean ground temperature (for sine) or average air temperature over the year (for Kusuda)"
+    annotation (Dialog(group="General", enable=not datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.File));
+  parameter Real offTime=0
+    "Offset time from simulation time 0 until minimum ground temperature in s (for sine and Kusuda)"
+    annotation (Dialog(group="General", enable=datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Kusuda
+           or datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Sine));
+  parameter Modelica.Units.SI.TemperatureDifference ampTGro=0
+    "amplitude of ground temperature (for sine option) or of surface temperature [(maximum air temperature - minimum air temperature)/2] (for Kusuda)"
+    annotation (Dialog(group="General", enable=datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Kusuda
+           or datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Sine));
+  parameter String filDatSou="NoName"
+    "path to a (.mos) file that has the temperature stored in a table called TGround"
+    annotation (Dialog(group="General", enable=datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.File));
+  parameter Modelica.Units.SI.ThermalDiffusivity theDifKusPerDay=0.04
+    "Thermal diffusivity of the ground for Kusuda. Declare in m2/day!"
+    annotation (Dialog(group="Special Kusuda parameters", enable=datSou ==
+          AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Kusuda));
+  parameter Modelica.Units.SI.Distance depTGroKus=0
+    "Depth of ground temperature for Kusuda"
+    annotation (Dialog(group="Special Kusuda parameters", enable=datSou ==
+          AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Kusuda));
+
+  // add conditions
+  AixLib.BoundaryConditions.GroundTemperature.GroundTemperatureKusuda kusUnd(
+    T_mean=TMea,
+    T_amp=ampTGro,
+    D=depTGroKus,
+    alpha=theDifKusPerDay,
+    t_shift=offTime/86400.0) if datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Kusuda
+    annotation (Placement(transformation(extent={{-10,50},{10,70}})));
+  Modelica.Thermal.HeatTransfer.Sensors.TemperatureSensor TSenKus if datSou ==
+    AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Kusuda
+    "Sensor to show Kusuda ground temperature"
+    annotation (Placement(transformation(extent={{22,50},{42,70}})));
+  Modelica.Blocks.Interfaces.RealOutput TGroOut(
+    final quantity="ThermodynamicTemperature",
+    final unit="K",
+    displayUnit="degC") "Ground temperature"
+    annotation (Placement(transformation(extent={{100,-10},{120,10}})));
+  Modelica.Blocks.Sources.Sine sin(
+    amplitude=ampTGro,
+    f=1.0/31536000.0,
+    phase=-offTime/31536000.0*2*pi - pi/2,
+    offset=TMea,
+    startTime=-31536000.0) if datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Sine
+    "Source for sine function for temperature"
+    annotation (Placement(transformation(extent={{-10,-28},{10,-8}})));
+  Modelica.Blocks.Sources.Constant TGroConSou(final k=TMea) if datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.Constant
+    "Constant ground temperature"
+    annotation (Placement(transformation(extent={{10,-10},{-10,10}}, rotation=180,
+        origin={0,22})));
+  Modelica.Blocks.Sources.CombiTimeTable tabTGro(
+    tableOnFile=true,
+    tableName="TGround",
+    extrapolation=Modelica.Blocks.Types.Extrapolation.Periodic,
+    fileName=filDatSou) if datSou == AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.File
+    "Ground temperatures from table"
+    annotation (Placement(transformation(extent={{-10,-66},{10,-46}})));
+
+equation
+  connect(kusUnd.port, TSenKus.port) annotation (Line(points={{9.4,55},{15.7,55},
+          {15.7,60},{22,60}}, color={191,0,0}));
+  connect(TSenKus.T, TGroOut) annotation (Line(points={{43,60},{56,60},{56,0},{110,
+          0}}, color={0,0,127}));
+  connect(TGroConSou.y, TGroOut)
+    annotation (Line(points={{11,22},{56,22},{56,0},{110,0}},
+                                                            color={0,0,127}));
+  connect(sin.y, TGroOut) annotation (Line(points={{11,-18},{56,-18},{56,0},{
+          110,0}},
+               color={0,0,127}));
+  connect(tabTGro.y[1], TGroOut) annotation (Line(points={{11,-56},{56,-56},{56,
+          0},{110,0}}, color={0,0,127}));
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false), graphics={
+        Polygon(
+          points={{-100,20},{-46,14},{-10,20},{24,16},{60,18},{78,12},{100,20},{
+              100,-100},{-100,-100},{-100,20}},
+          lineColor={0,0,0},
+          fillColor={127,66,38},
+          fillPattern=FillPattern.Solid),
+        Polygon(
+          points={{-100,20},{-100,32},{-42,22},{-4,30},{32,22},{66,24},{82,20},{
+              100,28},{100,20},{78,12},{60,18},{24,16},{-10,20},{-46,14},{-100,20}},
+          lineColor={0,0,0},
+          fillColor={0,140,72},
+          fillPattern=FillPattern.Solid),
+        Polygon(
+          points={{-100,100},{100,100},{100,28},{82,20},{66,24},{32,22},{-4,30},
+              {-42,22},{-100,32},{-100,100}},
+          lineColor={0,0,0},
+          fillColor={0,157,255},
+          fillPattern=FillPattern.Solid),
+        Ellipse(
+          extent={{-70,84},{-40,54}},
+          lineColor={0,0,0},
+          fillColor={255,255,85},
+          fillPattern=FillPattern.Solid),
+        Text(
+          extent={{-139,-104},{161,-144}},
+          lineColor={0,0,255},
+          textString="%name")}), Diagram(
+        coordinateSystem(preserveAspectRatio=false)),
+    Documentation(revisions="<html><ul>
+    <li>April 23, 2023, by Philip Groesdonk:<br/>
+    First implementation (see <a href=
+    \"https://github.com/RWTH-EBC/AixLib/issues/1080\">issue 1080</a>).
+  </li>
+</ul>
+</html>", info="<html>
+This model provides a block to choose from different options to model ground temperature:
+<ul>
+<li>
+a constant value
+  </li>
+<li>
+a sine curve (yearly)
+  </li>
+  <li>
+  a table stored in a file
+  </li>
+  <li>a simple model (called Kusuda here) that was originally created for use with district heating grids or LTN in order to model
+  thermal pipe losses
+  </li>
+</ul>
+<h4>
+  Assumption and limitations
+</h4>
+<p>
+The model does not model any influence on the ground temperature. 
+The ground temperature is only dependent on the set parameters.
+</p>
+<h4>
+  Typical use and important parameters
+</h4>
+<p>
+The model is used as a boundary condition for ROM model floors outside surfaces.
+</p>
+<p>
+<b>Differently then stated in the parameter description below for Kusuda,
+  α needs to be declared in m2/day! (The Modelica SI unit diffusivity does not
+  support m2/day as a display unit.)</b>
+</p>
+<p>
+  A typical value for α is bewteen 0.03 and 0.05 m2/day.
+</p>
+
+</html>
+"));
+end Options;
diff --git a/AixLib/BoundaryConditions/GroundTemperature/package.order b/AixLib/BoundaryConditions/GroundTemperature/package.order
@@ -1,2 +1,4 @@
 GroundTemperatureKusuda
+Options
 Examples
+DataSourceGroTem
diff --git a/.../DataBase/ThermalZones/MultizoneInterzonalsFixedHeater/MorschenichSfhFixedHeater_attic.mo b/.../DataBase/ThermalZones/MultizoneInterzonalsFixedHeater/MorschenichSfhFixedHeater_attic.mo
@@ -0,0 +1,107 @@
+within AixLib.DataBase.ThermalZones.MultizoneInterzonalsFixedHeater;
+record MorschenichSfhFixedHeater_attic
+  "Record for attic of MorschenichSfhFixedHeater example"
+  extends AixLib.DataBase.ThermalZones.ZoneBaseRecord(
+    T_start = 293.15,
+    withAirCap = true,
+    VAir = 75.78795696372634,
+    AZone = 50.52377022336238,
+    hRad = 7.750349138760984,
+    lat = 0.88645272708792,
+    nOrientations = 2,
+    AWin = {0.14869610304800834, 0.0854780546614556},
+    ATransparent = {0.14869610304800834, 0.0854780546614556},
+    hConWin = 2.7000000000000006,
+    RWin = 1.5215845888330657,
+    gWin = 0.6,
+    UWin= 1.9004689468829974,
+    ratioWinConRad = 0.02,
+    AExt = {7.82825722254796, 8.06360011905509},
+    hConExt = 2.7000000000000006,
+    nExt = 1,
+    RExt = {0.004230796101479976},
+    RExtRem = 0.02032557242299069,
+    CExt = {6273497.6266179895},
+    AInt = 0.0,
+    hConInt = 0.0,
+    nInt = 1,
+    RInt = {0.0},
+    CInt = {0.0},
+    AFloor = 0.0,
+    hConFloor = 0.0,
+    nFloor = 1,
+    RFloor = {0.0},
+    RFloorRem =  0.0,
+    CFloor = {0.0},
+    ARoof = 75.72370770529446,
+    hConRoof = 4.999999999999999,
+    nRoof = 1,
+    RRoof = {1.9138990862524574e-05},
+    RRoofRem = 9.569495530074631e-05,
+    CRoof = {2362579.5220066374},
+    nOrientationsRoof = 2,
+    tiltRoof = {0.8377580409572781, 0.8552113334772214},
+    aziRoof = {1.239183768915974, -1.902408884673819},
+    wfRoof = {0.5020416925535248, 0.49795830744647523},
+    aRoof = 0.4,
+    aExt = 0.5999999999999999,
+    nIze = 1,
+    AIze = {50.52377022336238},
+    hConIze = {5.0},
+    nIzeRC = 1,
+    RIze = {{0.01989650625627703}},
+    RIzeRem = {0.003523280978428377},
+    CIze = {{7416783.851736545}},
+    othZoneInd = {3},
+    zoneInd = 1,
+    TSoil = 286.15,
+    TSoiDatSou = AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.File,
+    TSoiOffTim = 6004800,
+    TSoiAmp = 0,
+    TSoiFil = Modelica.Utilities.Files.loadResource("modelica://AixLib/Resources/LowOrder_ExampleData/t_soil_MorschenichSfhFixedHeater.txt"),
+    hConWallOut = 20.0,
+    hRadWall = 5.0,
+    hConWinOut = 20.0,
+    hConRoofOut = 19.999999999999996,
+    hRadRoof = 4.999999999999999,
+    tiltExtWalls = {1.5707963267948966, 1.5707963267948966},
+    aziExtWalls = {-0.47123889803846897, 2.670353755551324},
+    wfWall = {0.49259548800847114, 0.5074045119915288},
+    wfWin = {0.6349808386307645, 0.36501916136923546},
+    wfGro = 0.0,
+    specificPeople = 0.02,
+    internalGainsMoistureNoPeople = 0.5,
+    fixedHeatFlowRatePersons = 70,
+    activityDegree = 1.2,
+    ratioConvectiveHeatPeople = 0.5,
+    internalGainsMachinesSpecific = 2.0,
+    ratioConvectiveHeatMachines = 0.75,
+    lightingPowerSpecific = 7.0,
+    ratioConvectiveHeatLighting = 0.5,
+    useConstantACHrate = true,
+    baseACH = 10.0,
+    maxUserACH = 1.0,
+    maxOverheatingACH = {3.0, 2.0},
+    maxSummerACH = {1.0, 283.15, 290.15},
+    winterReduction = {0.2, 273.15, 283.15},
+    maxIrr = {10000.0, 10000.0},
+    shadingFactor = {0.0, 0.0},
+    withAHU = false,
+    minAHU = 0.3,
+    maxAHU = 0.6,
+    hHeat = 25655.980698795167,
+    lHeat = 0,
+    KRHeat = 10000,
+    TNHeat = 1,
+    HeaterOn = false,
+    hCool = 0,
+    lCool = -25655.980698795167,
+    heaLoadFacOut = 801.749396837349,
+    heaLoadFacGrd = 0.0,
+    KRCool = 10000,
+    TNCool = 1,
+    CoolerOn = false,
+    withIdealThresholds = false,
+    TThresholdHeater = 288.15,
+    TThresholdCooler = 295.15);
+end MorschenichSfhFixedHeater_attic;
diff --git a/...taBase/ThermalZones/MultizoneInterzonalsFixedHeater/MorschenichSfhFixedHeater_basement.mo b/...taBase/ThermalZones/MultizoneInterzonalsFixedHeater/MorschenichSfhFixedHeater_basement.mo
@@ -0,0 +1,107 @@
+within AixLib.DataBase.ThermalZones.MultizoneInterzonalsFixedHeater;
+record MorschenichSfhFixedHeater_basement
+  "Record for basement of MorschenichSfhFixedHeater"
+  extends AixLib.DataBase.ThermalZones.ZoneBaseRecord(
+    T_start = 293.15,
+    withAirCap = true,
+    VAir = 176.72156159931643,
+    AZone = 79.60430702671911,
+    hRad = 6.654164349463362,
+    lat = 0.88645272708792,
+    nOrientations = 4,
+    AWin = {0.0, 0.0001, 0.0, 0.0},
+    ATransparent = {0.0, 0.0001, 0.0, 0.0},
+    hConWin = 2.7,
+    RWin = 1871.4285714285713,
+    gWin = 0.75,
+    UWin= 2.8010185522008,
+    ratioWinConRad = 0.020000000000000004,
+    AExt = {6.418846336889531, 7.081167133233686, 6.419032155135023, 7.167126477220699},
+    hConExt = 2.6999999999999993,
+    nExt = 1,
+    RExt = {0.004229837490597628},
+    RExtRem = 0.02060924234483541,
+    CExt = {5576900.6389134005},
+    AInt = 81.6764,
+    hConInt = 2.7000000000000006,
+    nInt = 1,
+    RInt = {0.00022719717562130076},
+    CInt = {6212483.036120574},
+    AFloor = 131.85521768744456,
+    hConFloor = 0.8999999999999999,
+    nFloor = 1,
+    RFloor = {0.0008556017690475546},
+    RFloorRem =  0.0036349634514251787,
+    CFloor = {16730172.436092455},
+    ARoof = 0.0,
+    hConRoof = 0.0,
+    nRoof = 1,
+    RRoof = {0.0},
+    RRoofRem = 0.0,
+    CRoof = {0.0},
+    nOrientationsRoof = 1,
+    tiltRoof = {0},
+    aziRoof = {0},
+    wfRoof = {0},
+    aRoof = 0.0,
+    aExt = 0.5,
+    nIze = 1,
+    AIze = {79.60430702671911},
+    hConIze = {0.8999999999999999},
+    nIzeRC = 1,
+    RIze = {{0.0060208915205140715}},
+    RIzeRem = {0.0014172041906433724},
+    CIze = {{10100425.349640682}},
+    othZoneInd = {3},
+    zoneInd = 2,
+    TSoil = 286.15,
+    TSoiDatSou = AixLib.BoundaryConditions.GroundTemperature.DataSourceGroTem.File,
+    TSoiOffTim = 6004800,
+    TSoiAmp = 0,
+    TSoiFil = Modelica.Utilities.Files.loadResource("modelica://AixLib/Resources/LowOrder_ExampleData/t_soil_MorschenichSfhFixedHeater.txt"),
+    hConWallOut = 19.999999999999996,
+    hRadWall = 4.999999999999999,
+    hConWinOut = 20.0,
+    hConRoofOut = 0.0,
+    hRadRoof = 0.0,
+    tiltExtWalls = {1.5707963267948966, 1.5707963267948966, 1.5707963267948966, 1.5707963267948966},
+    aziExtWalls = {-0.47123889803846897, 1.0995574287564276, 2.670353755551324, -2.0420352248333655},
+    wfWall = {0.23697871787139957, 0.2614310765818997, 0.23698557813370572, 0.26460462741299495},
+    wfWin = {0.0, 1.0, 0.0, 0.0},
+    wfGro = 0.0,
+    specificPeople = 0.02,
+    internalGainsMoistureNoPeople = 0.5,
+    fixedHeatFlowRatePersons = 70,
+    activityDegree = 1.2,
+    ratioConvectiveHeatPeople = 0.5,
+    internalGainsMachinesSpecific = 2.0,
+    ratioConvectiveHeatMachines = 0.75,
+    lightingPowerSpecific = 7.0,
+    ratioConvectiveHeatLighting = 0.5,
+    useConstantACHrate = true,
+    baseACH = 1.0,
+    maxUserACH = 1.0,
+    maxOverheatingACH = {3.0, 2.0},
+    maxSummerACH = {1.0, 283.15, 290.15},
+    winterReduction = {0.2, 273.15, 283.15},
+    maxIrr = {9999.9, 10000.0, 9999.9, 9999.9},
+    shadingFactor = {1.0, 0.0, 1.0, 1.0},
+    withAHU = false,
+    minAHU = 0.3,
+    maxAHU = 0.6,
+    hHeat = 4206.03587151335,
+    lHeat = 0,
+    KRHeat = 10000,
+    TNHeat = 1,
+    HeaterOn = false,
+    hCool = 0,
+    lCool = -4206.03587151335,
+    heaLoadFacOut = 93.56640313904289,
+    heaLoadFacGrd = 173.13013872342532,
+    KRCool = 10000,
+    TNCool = 1,
+    CoolerOn = false,
+    withIdealThresholds = false,
+    TThresholdHeater = 288.15,
+    TThresholdCooler = 295.15);
+end MorschenichSfhFixedHeater_basement;