Lithium Bromide-Water Heat Pump Media and Components

Models/NHES/Fluid/ClosureModels/HeatTransfer/package.order

@@ -2,4 +2,4 @@ Examples
 alpha_Chen_TwoPhase
 alpha_DittusBoelter
 kandlikar_lowP
-Nu_DittusBoelter_Heat_Cool
+Nu_DittusBoelter_Heat_Cool

Models/NHES/Fluid/HeatExchangers/Generic_HXs/Generic_HXR_NoMass.mo

@@ -0,0 +1,127 @@
+within NHES.Fluid.HeatExchangers.Generic_HXs;
+model Generic_HXR_NoMass
+
+replaceable package Medium_h =
+    Media.TwoPhaseMixtures.LithiumBromideWater.LithiumBromideWater_pTX constrainedby
+    Modelica.Media.Interfaces.PartialMedium
+                                          annotation (choicesAllMatching);
+replaceable package Medium_c =
+    Media.TwoPhaseMixtures.LithiumBromideWater.LithiumBromideWater_pTX constrainedby
+    Modelica.Media.Interfaces.PartialMedium
+                                          annotation (choicesAllMatching);
+
+parameter Real epsilon(min=0, max=1)=0.9 "Heat exchanger effectiveness";
+Modelica.Units.SI.EntropyFlowRate UA(start=1603e3, displayUnit="kW/K")
+  "Overall HTC x HT area";
+
+//SI.EntropyFlowRate C_h "Mean hot side heat rate";
+//SI.EntropyFlowRate C_c "Mean cold side heat rate";
+Medium_h.Temperature T_hi(start=545.7) = Modelica.Fluid.Utilities.regStep(
+  port_hi.m_flow,
+  Medium_h.temperature(hi_state),
+  Medium_h.temperature(Medium_h.setState_phX(
+    port_hi.p,
+    port_hi.h_outflow,
+    port_hi.Xi_outflow)));
+
+Medium_h.Temperature T_ho(start=322.9) = Modelica.Fluid.Utilities.regStep(
+  port_ho.m_flow,
+  Medium_h.temperature(ho_state),
+  Medium_h.temperature(Medium_h.setState_phX(
+    port_ho.p,
+    port_ho.h_outflow,
+    port_ho.Xi_outflow)));
+
+Medium_c.Temperature T_ci(start=298.2) = Modelica.Fluid.Utilities.regStep(
+  port_ci.m_flow,
+  Medium_c.temperature(ci_state),
+  Medium_c.temperature(Medium_c.setState_phX(
+    port_ci.p,
+    port_ci.h_outflow,
+    port_ci.Xi_outflow)));
+
+Medium_c.Temperature T_co(start=485.7) = Modelica.Fluid.Utilities.regStep(
+  port_co.m_flow,
+  Medium_c.temperature(co_state),
+  Medium_c.temperature(Medium_c.setState_phX(
+    port_co.p,
+    port_co.h_outflow,
+    port_co.Xi_outflow)));
+
+ //Medium_h.SpecificEnthalpy h_ho(start=170.7e3) = inStream(port_ho.h_outflow);
+ //Medium_c.SpecificEnthalpy h_co(start=458.3e3) = inStream(port_co.h_outflow);
+
+SI.Power Q_HT(start=63799e3) "Heat transfer";
+SI.TemperatureDifference LMTD=Utilities.Functions.LMTD(
+      T_hi,
+      T_ho,
+      T_ci,
+      T_co,
+      true) "Log-mean temperature difference";
+
+outer Modelica.Fluid.System system;
+
+Modelica.Fluid.Interfaces.FluidPort_a port_ci(redeclare package Medium =
+      Medium_c)
+  annotation (Placement(transformation(extent={{-70,-110},{-50,-90}})));
+Modelica.Fluid.Interfaces.FluidPort_b port_co(redeclare package Medium =
+      Medium_c,
+      h_outflow(start=458.3e3))
+  annotation (Placement(transformation(extent={{-70,90},{-50,110}})));
+Modelica.Fluid.Interfaces.FluidPort_a port_hi(redeclare package Medium =
+      Medium_h)
+  annotation (Placement(transformation(extent={{50,90},{70,110}})));
+Modelica.Fluid.Interfaces.FluidPort_b port_ho(redeclare package Medium =
+      Medium_h,
+      h_outflow(start=170.7e3))
+  annotation (Placement(transformation(extent={{50,-110},{70,-90}})));
+
+Medium_h.ThermodynamicState hi_state;
+Medium_h.ThermodynamicState ho_state;
+Medium_c.ThermodynamicState ci_state;
+Medium_c.ThermodynamicState co_state;
+
+equation
+port_hi.m_flow + port_ho.m_flow = 0;
+port_ci.m_flow + port_co.m_flow = 0;
+
+port_hi.Xi_outflow = inStream(port_ho.Xi_outflow);
+port_ho.Xi_outflow = inStream(port_hi.Xi_outflow);
+
+port_ci.Xi_outflow = inStream(port_co.Xi_outflow);
+port_co.Xi_outflow = inStream(port_ci.Xi_outflow);
+
+port_hi.p = port_ho.p;
+port_ci.p = port_co.p;
+
+port_hi.h_outflow = inStream(port_ho.h_outflow);
+Q_HT = port_ho.m_flow*(port_ho.h_outflow - inStream(port_hi.h_outflow));
+
+port_ci.h_outflow = inStream(port_co.h_outflow);
+Q_HT = port_co.m_flow*(inStream(port_ci.h_outflow) - port_co.h_outflow);
+
+//C_h = port_hi.m_flow*Medium_h.specificHeatCapacityCp(hi_state);
+//C_c = port_ci.m_flow*Medium_c.specificHeatCapacityCp(ci_state);
+
+Q_HT = UA*LMTD;
+
+epsilon = Utilities.Functions.epsilon(T_hi, T_co, T_ho, T_ci);
+
+// Medium states for inflowing fluid
+hi_state = Medium_h.setState_phX(
+  port_hi.p,
+  inStream(port_hi.h_outflow),
+  inStream(port_hi.Xi_outflow));
+ho_state = Medium_h.setState_phX(
+  port_ho.p,
+  inStream(port_ho.h_outflow),
+  inStream(port_ho.Xi_outflow));
+ci_state = Medium_c.setState_phX(
+  port_ci.p,
+  inStream(port_ci.h_outflow),
+  inStream(port_ci.Xi_outflow));
+co_state = Medium_c.setState_phX(
+  port_co.p,
+  inStream(port_co.h_outflow),
+  inStream(port_co.Xi_outflow));
+end Generic_HXR_NoMass;

Models/NHES/Fluid/HeatExchangers/Generic_HXs/package.mo

@@ -3,4 +3,129 @@ package Generic_HXs "Generic heat exchangers"
 
 
 
+  model recuperator_simple
+
+  replaceable package Medium_h =
+      Media.TwoPhaseMixtures.LithiumBromideWater.LithiumBromideWater_pTX constrainedby
+    Modelica.Media.Interfaces.PartialMedium annotation (choicesAllMatching);
+  replaceable package Medium_c =
+      Media.TwoPhaseMixtures.LithiumBromideWater.LithiumBromideWater_pTX constrainedby
+    Modelica.Media.Interfaces.PartialMedium annotation (choicesAllMatching);
+
+  parameter Real epsilon(min=0, max=1)=0.9 "Heat exchanger effectiveness";
+  Modelica.Units.SI.EntropyFlowRate UA(start=1603e3, displayUnit="kW/K")
+    "Overall HTC x HT area";
+
+  //SI.EntropyFlowRate C_h "Mean hot side heat rate";
+  //SI.EntropyFlowRate C_c "Mean cold side heat rate";
+  Medium_h.Temperature T_hi(start=545.7) = Modelica.Fluid.Utilities.regStep(
+    port_hi.m_flow,
+    Medium_h.temperature(hi_state),
+    Medium_h.temperature(Medium_h.setState_phX(
+      port_hi.p,
+      port_hi.h_outflow,
+      port_hi.Xi_outflow)));
+
+  Medium_h.Temperature T_ho(start=322.9) = Modelica.Fluid.Utilities.regStep(
+    port_ho.m_flow,
+    Medium_h.temperature(ho_state),
+    Medium_h.temperature(Medium_h.setState_phX(
+      port_ho.p,
+      port_ho.h_outflow,
+      port_ho.Xi_outflow)));
+
+  Medium_c.Temperature T_ci(start=298.2) = Modelica.Fluid.Utilities.regStep(
+    port_ci.m_flow,
+    Medium_c.temperature(ci_state),
+    Medium_c.temperature(Medium_c.setState_phX(
+      port_ci.p,
+      port_ci.h_outflow,
+      port_ci.Xi_outflow)));
+
+  Medium_c.Temperature T_co(start=485.7) = Modelica.Fluid.Utilities.regStep(
+    port_co.m_flow,
+    Medium_c.temperature(co_state),
+    Medium_c.temperature(Medium_c.setState_phX(
+      port_co.p,
+      port_co.h_outflow,
+      port_co.Xi_outflow)));
+
+   //Medium_h.SpecificEnthalpy h_ho(start=170.7e3) = inStream(port_ho.h_outflow);
+   //Medium_c.SpecificEnthalpy h_co(start=458.3e3) = inStream(port_co.h_outflow);
+
+  SI.Power Q_HT(start=63799e3) "Heat transfer";
+  SI.TemperatureDifference LMTD=Utilities.Functions.LMTD(
+        T_hi,
+        T_ho,
+        T_ci,
+        T_co,
+        true) "Log-mean temperature difference";
+
+  outer Modelica.Fluid.System system;
+
+  Modelica.Fluid.Interfaces.FluidPort_a port_ci(redeclare package Medium =
+        Medium_c)
+    annotation (Placement(transformation(extent={{-70,-110},{-50,-90}})));
+  Modelica.Fluid.Interfaces.FluidPort_b port_co(redeclare package Medium =
+        Medium_c,
+        h_outflow(start=458.3e3))
+    annotation (Placement(transformation(extent={{-70,90},{-50,110}})));
+  Modelica.Fluid.Interfaces.FluidPort_a port_hi(redeclare package Medium =
+        Medium_h)
+    annotation (Placement(transformation(extent={{50,90},{70,110}})));
+  Modelica.Fluid.Interfaces.FluidPort_b port_ho(redeclare package Medium =
+        Medium_h,
+        h_outflow(start=170.7e3))
+    annotation (Placement(transformation(extent={{50,-110},{70,-90}})));
+
+  Medium_h.ThermodynamicState hi_state;
+  Medium_h.ThermodynamicState ho_state;
+  Medium_c.ThermodynamicState ci_state;
+  Medium_c.ThermodynamicState co_state;
+
+  equation
+  port_hi.m_flow + port_ho.m_flow = 0;
+  port_ci.m_flow + port_co.m_flow = 0;
+
+  port_hi.Xi_outflow = inStream(port_ho.Xi_outflow);
+  port_ho.Xi_outflow = inStream(port_hi.Xi_outflow);
+
+  port_ci.Xi_outflow = inStream(port_co.Xi_outflow);
+  port_co.Xi_outflow = inStream(port_ci.Xi_outflow);
+
+  port_hi.p = port_ho.p;
+  port_ci.p = port_co.p;
+
+  port_hi.h_outflow = inStream(port_ho.h_outflow);
+  Q_HT = port_ho.m_flow*(port_ho.h_outflow - inStream(port_hi.h_outflow));
+
+  port_ci.h_outflow = inStream(port_co.h_outflow);
+  Q_HT = port_co.m_flow*(inStream(port_ci.h_outflow) - port_co.h_outflow);
+
+  //C_h = port_hi.m_flow*Medium_h.specificHeatCapacityCp(hi_state);
+  //C_c = port_ci.m_flow*Medium_c.specificHeatCapacityCp(ci_state);
+
+  Q_HT = UA*LMTD;
+
+  epsilon = Utilities.Functions.epsilon(T_hi, T_co, T_ho, T_ci);
+
+  // Medium states for inflowing fluid
+  hi_state = Medium_h.setState_phX(
+    port_hi.p,
+    inStream(port_hi.h_outflow),
+    inStream(port_hi.Xi_outflow));
+  ho_state = Medium_h.setState_phX(
+    port_ho.p,
+    inStream(port_ho.h_outflow),
+    inStream(port_ho.Xi_outflow));
+  ci_state = Medium_c.setState_phX(
+    port_ci.p,
+    inStream(port_ci.h_outflow),
+    inStream(port_ci.Xi_outflow));
+  co_state = Medium_c.setState_phX(
+    port_co.p,
+    inStream(port_co.h_outflow),
+    inStream(port_co.Xi_outflow));
+  end recuperator_simple;
+
 end Generic_HXs;

Models/NHES/Fluid/HeatExchangers/Generic_HXs/package.order

@@ -1,3 +1,5 @@
 Examples
 Generic_HX
 Generic_STHX
+Generic_HXR_NoMass
+recuperator_simple

Models/NHES/Fluid/HeatExchangers/Utilities/Functions/LMTD.mo

@@ -26,7 +26,7 @@ algorithm
                    else (DT_1 - DT_2)/log(max(0.000001,DT_1/DT_2)));
 
   annotation (Documentation(info="<html>
-<p>The Peclet numver is defined to be the ratio of the rate of&nbsp;advection&nbsp;of a physical quantity by the flow to the rate of&nbsp;diffusion&nbsp;of the same quantity driven by an appropriate gradient.</p>
+<p>The Peclet number is defined to be the ratio of the rate of&nbsp;advection&nbsp;of a physical quantity by the flow to the rate of&nbsp;diffusion&nbsp;of the same quantity driven by an appropriate gradient.</p>
 <p> In the context of species or&nbsp;mass transfer, the P&eacute;clet number is the product of the&nbsp;Reynolds number&nbsp;and the&nbsp;Schmidt number. </p>
 <ul>
 <li>Pe = Re*Sc</li>

Models/NHES/Fluid/HeatExchangers/Utilities/Functions/package.mo

@@ -2,4 +2,21 @@ within NHES.Fluid.HeatExchangers.Utilities;
 package Functions
 
 
+function epsilon
+  extends Modelica.Icons.Function;
+  input Modelica.Units.SI.Temperature High_1;
+  input Modelica.Units.SI.Temperature Low_1;
+  input Modelica.Units.SI.Temperature High_2;
+  input Modelica.Units.SI.Temperature Low_2;
+  output Real eps;
+protected
+  Modelica.Units.SI.Temperature DELTAT_den=High_1 - Low_2;
+algorithm
+  if
+    (DELTAT_den<=0) then
+    eps :=0.8;
+  else
+    eps := if (High_1 - High_2)/DELTAT_den > 1 then 0.8 else (High_1 - High_2)/DELTAT_den;
+  end if;
+end epsilon;
 end Functions;

Models/NHES/Fluid/HeatExchangers/Utilities/Functions/package.order

@@ -1,2 +1,3 @@
 LMTD
 UA
+epsilon