Issue1576 Heat pump model integration

IBPSA/Fluid/Chillers/BaseClasses/BlackBoxVapourCompressionCycle.mo

@@ -0,0 +1,192 @@
+within IBPSA.Fluid.Chillers.BaseClasses;
+model BlackBoxVapourCompressionCycle
+  "Blackbox model of refrigerant cycle of a chiller"
+  extends
+    IBPSA.Fluid.HeatPumps.BaseClasses.PartialBlackBoxVapourCompressionCycle;
+  // Setting all values to zero avoids errors when checking this model.
+  // The values are correctly propagated by the heat pump / chiller model anyway
+  replaceable model BlackBoxChillerCooling =
+      IBPSA.Fluid.Chillers.BlackBoxData.BaseClasses.PartialChillerBlackBox (
+      QUse_flow_nominal=0,
+      QUseBlaBox_flow_nominal=0,
+      scaFac=0,
+      TCon_nominal=0,
+      TEva_nominal=0,
+      dTCon_nominal=0,
+      dTEva_nominal=0,
+      mCon_flow_nominal=0,
+      mEva_flow_nominal=0,
+      y_nominal=0)
+    constrainedby
+    IBPSA.Fluid.Chillers.BlackBoxData.BaseClasses.PartialChillerBlackBox
+    "Replaceable model for black-box model of a chiller in main operation mode"
+    annotation (choicesAllMatching=true);
+
+  replaceable model BlackBoxChillerHeating =
+      IBPSA.Fluid.HeatPumps.BlackBoxData.BaseClasses.NoHeating
+    constrainedby
+    IBPSA.Fluid.HeatPumps.BlackBoxData.BaseClasses.PartialHeatPumpBlackBox
+    "Replaceable model for black-box model of a chiller in reversible operation mode"
+    annotation (Dialog(enable=use_rev),choicesAllMatching=true);
+
+  BlackBoxChillerCooling blaBoxChiCoo
+    annotation (Placement(transformation(extent={{21,40},{60,80}},rotation=0)));
+  BlackBoxChillerHeating blaBoxChiHea
+    annotation (Placement(transformation(extent={{-80,38},{-39,80}}, rotation=0)));
+
+  Modelica.Blocks.Math.Gain gainEva(final k=-1)
+    "Negate QEva to match definition of heat flow direction" annotation (Placement(
+        transformation(extent={{-10,-10},{10,10}},
+                                               rotation=180,
+        origin={-10,12})));
+  Modelica.Blocks.Math.Gain gainCon(final k=-1) if use_rev
+    "Negate QCon to match definition of heat flow direction" annotation (
+      Placement(transformation(
+        extent={{-10,-10},{10,10}},
+        rotation=0,
+        origin={30,-8})));
+protected
+  IBPSA.Utilities.IO.Strings.StringPassThrough strPasThr
+    "String pass through to enable conditional string data";
+  IBPSA.Utilities.IO.Strings.ConstStringSource conStrSou(
+    final k=blaBoxChiCoo.datSou) if not use_rev
+   "Constant String data source";
+initial equation
+  assert(
+    strPasThr.y == blaBoxChiCoo.datSou,
+    "Data sources for reversible operation are not equal! 
+    Only continue if this is intended",
+    AssertionLevel.warning);
+equation
+  connect(conStrSou.y, strPasThr.u);
+  connect(blaBoxChiHea.datSouOut, strPasThr.u);
+  connect(blaBoxChiCoo.Pel, switchPel.u1) annotation (Line(points={{40.5,38},{
+          40,38},{40,8},{8,8},{8,-58}},
+                                 color={0,0,127}));
+  connect(blaBoxChiHea.Pel, switchPel.u3) annotation (Line(
+      points={{-59.5,35.9},{-58,35.9},{-58,26},{-24,26},{-24,-58},{-8,-58}},
+      color={0,0,127},
+      pattern=LinePattern.Dash));
+  connect(constZero.y, switchPel.u3) annotation (Line(points={{-59,-70},{-34,-70},
+          {-34,-58},{-8,-58}},
+      color={0,0,127},
+      pattern=LinePattern.Dash));
+  connect(constZero.y, switchQEva.u3) annotation (Line(
+      points={{-59,-70},{-52,-70},{-52,-8},{-58,-8}},
+      color={0,0,127},
+      pattern=LinePattern.Dash));
+  connect(constZero.y, switchQCon.u3) annotation (Line(points={{-59,-70},{-52,-70},
+          {-52,-38},{58,-38},{58,-8}},       color={0,0,127},
+      pattern=LinePattern.Dash));
+  connect(sigBus, blaBoxChiHea.sigBus) annotation (Line(
+      points={{0,102},{0,92},{-59.295,92},{-59.295,80.84}},
+      color={255,204,51},
+      thickness=0.5));
+  connect(sigBus, blaBoxChiCoo.sigBus) annotation (Line(
+      points={{0,102},{0,92},{40.695,92},{40.695,80.8}},
+      color={255,204,51},
+      thickness=0.5));
+
+  connect(blaBoxChiCoo.QEva_flow, gainEva.u) annotation (Line(points={{56.1,38},
+          {56.1,26},{14,26},{14,12},{2,12}},             color={0,0,127}));
+  connect(gainEva.y, switchQEva.u1) annotation (Line(points={{-21,12},{-50,12},
+          {-50,8},{-58,8}},                   color={0,0,127}));
+  connect(blaBoxChiHea.QEva_flow, switchQEva.u3) annotation (Line(
+      points={{-43.1,35.9},{-42,35.9},{-42,-10},{-52,-10},{-52,-8},{-58,-8}},
+      color={0,0,127},
+      pattern=LinePattern.Dash));
+  connect(blaBoxChiHea.QCon_flow, gainCon.u) annotation (Line(
+      points={{-75.9,35.9},{-74,35.9},{-74,20},{-34,20},{-34,-8},{18,-8}},
+      color={0,0,127},
+      pattern=LinePattern.Dash));
+  connect(gainCon.y, switchQCon.u3)
+    annotation (Line(points={{41,-8},{58,-8}},     color={0,0,127},
+        pattern=LinePattern.Dash));
+  connect(blaBoxChiCoo.QCon_flow, switchQCon.u1) annotation (Line(points={{24.9,38},
+          {24.9,36},{24,36},{24,24},{50,24},{50,8},{58,8}},
+                                                color={0,0,127}));
+  annotation (Icon(coordinateSystem(preserveAspectRatio=false), graphics={
+        Rectangle(
+          extent={{-100,100},{100,-100}},
+          lineColor={238,46,47},
+          fillColor={0,0,0},
+          fillPattern=FillPattern.Solid),
+        Ellipse(
+          extent={{-24,88},{22,44}},
+          lineColor={28,108,200},
+          lineThickness=0.5),
+        Line(
+          points={{-16,82},{20,74}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{-18,52},{20,58}},
+          color={28,108,200},
+          thickness=0.5),
+        Rectangle(
+          extent={{-98,40},{-60,-28}},
+          lineColor={28,108,200},
+          lineThickness=0.5),
+        Line(
+          points={{-20,-60},{-20,-70},{-20,-80},{20,-60},{20,-80},{-20,-60}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{-122,34},{-66,34},{-82,10},{-66,-22},{-120,-22}},
+          color={28,108,200},
+          thickness=0.5),
+        Rectangle(
+          extent={{60,40},{98,-28}},
+          lineColor={28,108,200},
+          lineThickness=0.5),
+        Line(
+          points={{120,34},{64,34},{80,10},{64,-22},{118,-22}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{-80,40},{-80,68},{-24,68}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{22,66},{80,66},{80,40}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{78,-28},{78,-70}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{78,-70},{62,-70},{20,-70}},
+          color={28,108,200},
+          thickness=0.5),
+        Line(
+          points={{-80,-26},{-80,-68},{-20,-68}},
+          color={28,108,200},
+          thickness=0.5),
+        Text(
+          extent={{-151,147},{149,107}},
+          textColor={0,0,255},
+          fillPattern=FillPattern.HorizontalCylinder,
+          fillColor={0,127,255},
+          textString="%name")}), Diagram(coordinateSystem(preserveAspectRatio=false)),
+    Documentation(revisions="<html><ul>
+  <li>
+    <i>May 22, 2019&#160;</i> by Julian Matthes:<br/>
+    First implementation (see issue <a href=
+    \"https://github.com/RWTH-EBC/AixLib/issues/715\">AixLib #715</a>)
+  </li>
+</ul>
+</html>", info="<html>
+<p>This black-box model represents the refrigerant cycle of a chiller. 
+Used in <a href=\"modelica://IBPSA.Fluid.Chillers.Chiller\">IBPSA.Fluid.Chillers.Chiller</a>, 
+this model serves the simulation of a 
+reversible chiller. Thus, data both of chillers and heat pumps can be used to
+calculate the three relevant values <code>P_el</code>, <code>QCon</code> and 
+<code>QEva</code>. The <code>mode</code> of the chiller is used to switch
+between the black-box model of the chiller and the heat pump. </p>
+ <p>The user can choose between different types of black-box model or implement 
+ a new black-box model by extending from 
+ <a href=\"modelica://IBPSA.Fluid.Chillers.BlackBoxData.BaseClasses.PartialChillerBlackBox\">
+ IBPSA.Fluid.Chillers.BlackBoxData.BaseClasses.PartialChillerBlackBox</a>. </p>
+</html>"));
+end BlackBoxVapourCompressionCycle;

IBPSA/Fluid/Chillers/BaseClasses/Carnot.mo

@@ -13,11 +13,11 @@ partial model Carnot
   parameter Modelica.Units.SI.HeatFlowRate QCon_flow_nominal(min=0)
     "Nominal heating flow rate" annotation (Dialog(group="Nominal condition"));
 
-  parameter Modelica.Units.SI.TemperatureDifference dTEva_nominal(final max=0)
-     = -10 "Temperature difference evaporator outlet-inlet"
+  parameter Modelica.Units.SI.TemperatureDifference dTEva_nominal(final max=0)=
+       -10 "Temperature difference evaporator outlet-inlet"
     annotation (Dialog(group="Nominal condition"));
-  parameter Modelica.Units.SI.TemperatureDifference dTCon_nominal(final min=0)
-     = 10 "Temperature difference condenser outlet-inlet"
+  parameter Modelica.Units.SI.TemperatureDifference dTCon_nominal(final min=0)=
+       10 "Temperature difference condenser outlet-inlet"
     annotation (Dialog(group="Nominal condition"));
 
   // Efficiency
@@ -52,12 +52,12 @@ partial model Carnot
     "Pressure difference over evaporator"
     annotation (Dialog(group="Nominal condition"));
 
-  parameter Modelica.Units.SI.TemperatureDifference TAppCon_nominal(min=0) =
+  parameter Modelica.Units.SI.TemperatureDifference TAppCon_nominal(min=0)=
     if cp1_default < 1500 then 5 else 2
     "Temperature difference between refrigerant and working fluid outlet in condenser"
     annotation (Dialog(group="Efficiency"));
 
-  parameter Modelica.Units.SI.TemperatureDifference TAppEva_nominal(min=0) =
+  parameter Modelica.Units.SI.TemperatureDifference TAppEva_nominal(min=0)=
     if cp2_default < 1500 then 5 else 2
     "Temperature difference between refrigerant and working fluid outlet in evaporator"
     annotation (Dialog(group="Efficiency"));
@@ -137,12 +137,12 @@ partial model Carnot
     x2=TConAct - TEvaAct,
     deltaX=0.25) "Carnot efficiency";
 
-  Modelica.Units.SI.Temperature TConAct(start=TCon_nominal + TAppCon_nominal)
-     = Medium1.temperature(staB1) + QCon_flow/QCon_flow_nominal*TAppCon_nominal
+  Modelica.Units.SI.Temperature TConAct(start=TCon_nominal + TAppCon_nominal)=
+       Medium1.temperature(staB1) + QCon_flow/QCon_flow_nominal*TAppCon_nominal
     "Condenser temperature used to compute efficiency, taking into account pinch temperature between fluid and refrigerant";
 
-  Modelica.Units.SI.Temperature TEvaAct(start=TEva_nominal - TAppEva_nominal)
-     = Medium2.temperature(staB2) - QEva_flow/QEva_flow_nominal*TAppEva_nominal
+  Modelica.Units.SI.Temperature TEvaAct(start=TEva_nominal - TAppEva_nominal)=
+       Medium2.temperature(staB2) - QEva_flow/QEva_flow_nominal*TAppEva_nominal
     "Evaporator temperature used to compute efficiency, taking into account pinch temperature between fluid and refrigerant";
 
 protected

IBPSA/Fluid/Chillers/BaseClasses/PartialCarnot_y.mo

@@ -36,9 +36,9 @@ partial model PartialCarnot_y
     annotation (Placement(transformation(extent={{-140,70},{-100,110}})));
 
 protected
-  Modelica.Units.SI.HeatFlowRate QCon_flow_internal(start=QCon_flow_nominal) =
+  Modelica.Units.SI.HeatFlowRate QCon_flow_internal(start=QCon_flow_nominal)=
     P - QEva_flow_internal "Condenser heat input";
-  Modelica.Units.SI.HeatFlowRate QEva_flow_internal(start=QEva_flow_nominal) =
+  Modelica.Units.SI.HeatFlowRate QEva_flow_internal(start=QEva_flow_nominal)=
     if COP_is_for_cooling then -COP*P else (1 - COP)*P "Evaporator heat input";
 
   Modelica.Blocks.Sources.RealExpression yEva_flow_in(

IBPSA/Fluid/Chillers/BaseClasses/package.order

@@ -1,3 +1,4 @@
+BlackBoxVapourCompressionCycle
 Carnot
 PartialCarnot_T
 PartialCarnot_y

IBPSA/Fluid/Chillers/BlackBoxData/BaseClasses/NoCooling.mo

@@ -0,0 +1,41 @@
+within IBPSA.Fluid.Chillers.BlackBoxData.BaseClasses;
+model NoCooling
+  "No Cooling. Used to avoid warnings about partial model must be replaced"
+  extends PartialChillerBlackBox(
+    redeclare final IBPSA.Fluid.HeatPumps.BlackBoxData.Frosting.NoFrosting
+      iceFacCalc,
+    datSou="",
+    QUseBlaBox_flow_nominal=0,
+    scaFac=0,
+    y_nominal=0,
+    mEva_flow_nominal=1,
+    mCon_flow_nominal=1,
+    dTEva_nominal=0,
+    dTCon_nominal=0,
+    TEva_nominal=273.15,
+    TCon_nominal=273.15,
+    QUse_flow_nominal=0);
+  Modelica.Blocks.Sources.Constant const(final k=0)
+    annotation (Placement(transformation(extent={{-88,16},{-68,36}})));
+equation
+  connect(const.y, feedbackHeatFlowEvaporator.u1) annotation (Line(points={{-67,
+          26},{-62,26},{-62,4},{-94,4},{-94,-10},{-78,-10}}, color={0,0,127}));
+  connect(const.y, feedbackHeatFlowEvaporator.u2) annotation (Line(points={{-67,
+          26},{-62,26},{-62,4},{-94,4},{-94,-10},{-84,-10},{-84,-24},{-70,-24},{
+          -70,-18}}, color={0,0,127}));
+  connect(const.y, calcRedQCon.u2) annotation (Line(points={{-67,26},{-62,26},{-62,
+          4},{64,4},{64,-58}}, color={0,0,127}));
+  connect(const.y, Pel)
+    annotation (Line(points={{-67,26},{0,26},{0,-110}}, color={0,0,127}));
+  annotation (Documentation(info="<html>
+  <p>Model to use to avoid warnings about partial model must be replaced. 
+  Using this model, the chiller will always be off.</p>
+</html>", revisions="<html><ul>
+  <li>
+    <i>October 2, 2022</i> by Fabian Wuellhorst:<br/>
+    First implementation (see issue <a href=
+    \"https://github.com/ibpsa/modelica-ibpsa/issues/1576\">#1576</a>)
+  </li>
+</ul>
+</html>"));
+end NoCooling;

IBPSA/Fluid/Chillers/BlackBoxData/BaseClasses/PartialChillerBlackBox.mo

@@ -0,0 +1,30 @@
+within IBPSA.Fluid.Chillers.BlackBoxData.BaseClasses;
+partial model PartialChillerBlackBox
+  "Partial black-box model of vapour compression cycles used for chiller applications"
+  extends IBPSA.Fluid.HeatPumps.BlackBoxData.BaseClasses.PartialBlackBox;
+
+equation
+  connect(iceFacCalc.iceFac, sigBus.iceFacChiMea) annotation (Line(points={{
+          -79,-42},{-72,-42},{-72,-26},{-102,-26},{-102,104},{1,104}}, color=
+          {0,0,127}), Text(
+      string="%second",
+      index=1,
+      extent={{6,3},{6,3}},
+      horizontalAlignment=TextAlignment.Left));
+  annotation (Icon(graphics={Rectangle(
+        extent={{-100,-100},{100,100}},
+        lineColor={0,0,255},
+        fillColor={255,255,255},
+        fillPattern=FillPattern.Solid),   Text(
+          extent={{-57.5,-35},{57.5,35}},
+          lineColor={0,0,255},
+          pattern=LinePattern.Dash,
+          textString="%name", origin={-3.5,-15},
+          rotation=180)}),
+      Documentation(
+  info="<html>
+<p>Partial model for calculation of electrical power <code>P_el</code>, 
+condenser heat flow <code>QCon</code> and evaporator heat flow <code>QEva</code> 
+based on the values in the sigBus for a chiller.</p>
+</html>"));
+end PartialChillerBlackBox;

IBPSA/Fluid/Chillers/BlackBoxData/BaseClasses/package.mo

@@ -0,0 +1,31 @@
+within IBPSA.Fluid.Chillers.BlackBoxData;
+package BaseClasses "Package with partial classes of performance Data"
+annotation (Icon(graphics={
+        Rectangle(
+          lineColor={200,200,200},
+          fillColor={248,248,248},
+          fillPattern=FillPattern.HorizontalCylinder,
+          extent={{-100.0,-100.0},{100.0,100.0}},
+          radius=25.0),
+        Rectangle(
+          lineColor={128,128,128},
+          extent={{-100.0,-100.0},{100.0,100.0}},
+          radius=25.0),
+        Ellipse(
+          extent={{-30.0,-30.0},{30.0,30.0}},
+          lineColor={128,128,128},
+          fillColor={255,255,255},
+          fillPattern=FillPattern.Solid)}), Documentation(revisions="<html><ul>
+  <li>
+    <i>May 22, 2019&#160;</i> by Julian Matthes:<br/>
+    First implementation (see issue <a href=
+    \"https://github.com/RWTH-EBC/AixLib/issues/715\">AixLib #715</a>)
+  </li>
+</ul>
+</html>", info="<html>
+<p>
+  This package contains base classes for the package <a href=
+  \"IBPSA.Fluid.Chillers.BlackBoxData\">IBPSA.Fluid.Chillers.BlackBoxData</a>.
+</p>
+</html>"));
+end BaseClasses;

IBPSA/Fluid/Chillers/BlackBoxData/BaseClasses/package.order

@@ -0,0 +1,2 @@
+NoCooling
+PartialChillerBlackBox