Added NOVEC649 incompressible fluid model, tested fluid system with it

985b5990 · i.elizarov · ef6b1944 · 985b5990 · 985b5990 · 985b5990
Commit 985b5990 authored Dec 10, 2021 by i.elizarov
--- a/Novec649.mo
+++ b/Novec649.mo
+package Novec649 "NOVEC649 incompressible model"
+  extends Modelica.Media.Incompressible.TableBased(
+    mediumName="NOVEC649",
+    T_min = Modelica.SIunits.Conversions.from_degC(-50), T_max = Modelica.SIunits.Conversions.from_degC(30),
+    TinK = false, T0=273.15,
+    tableDensity=
+      [-50, 1813; 0, 1676; 30, 1587.8],
+       
+    tableHeatCapacity=
+       [-50, 1099; 0, 1091; 30, 1105.4],
+    
+    // Replace later
+    tableConductivity=
+       [-50, 0.073; 0, 0.0634; 30, 0.058],
+
+    tableViscosity=
+       [-50, 2.175; 0, 0.87152; 30, 0.63455],
+       
+    tableVaporPressure=
+       [-50, 407.49; 0, 12365.28; 30, 49354]);
+    annotation (Documentation(info="<html>
+
+</html>"));
+end Novec649;
--- a/Piping_local.mo
+++ b/Piping_local.mo
 model Piping_local "Sections of channels"

-package Medium =
+replaceable package Medium =
      Modelica.Media.Water.ConstantPropertyLiquidWater;

 parameter Modelica.SIunits.Height channel_height = 40e-3 "Height of the channels (CH)";
 parameter Modelica.SIunits.Length channel_cross_section_width = 5e-3 "Width of cross-section of the channels";
-parameter Modelica.SIunits.Length channel_horizontal_lenght_middle = 42e-3 "Horizontal length of the middle channel (HLM)";
-parameter Modelica.SIunits.Length channel_horizontal_lenght_right_left = 29e-3 "Horizontal length of the right and left channel (HLRL)";
 parameter Modelica.SIunits.Height channel_cross_section_height = 1e-3 "Height of cross-section of the channels";
 parameter Modelica.SIunits.Height channel_roughness = 0.015e-3 "Rougness of wet surface";
-parameter Modelica.SIunits.Density novec_density = 1600 "NOVEC density";
 parameter Integer channel_number = 22 "Vertical array of channels";
 parameter Real epsilon_left_right(min = 0)= 0.5 "Local hydraulic resistance coefficient";
 parameter Real epsilon_left_right_middle(min = 0)= 0.5 "Local hydraulic resistance coefficient";
@@ -129,6 +126,6 @@ equation
    Line(points = {{80, 100}, {80, 40}, {50, 40}, {50, 10}}));
  annotation(
    uses(Modelica(version = "3.2.3")),
-    Icon(graphics = {Rectangle(fillColor = {255, 255, 255}, fillPattern = FillPattern.Solid, extent = {{-100, 100}, {100, -100}})}),
+    Icon(graphics = {Rectangle(fillColor = {255, 255, 255}, fillPattern = FillPattern.Solid, extent = {{-100, 100}, {100, -100}}), Line(origin = {-82, 1}, points = {{0, 43}, {0, -43}}), Line(origin = {-68, 2}, points = {{0, 44}, {0, -44}}), Line(origin = {-26, 5}, points = {{0, 45}, {0, -45}}), Line(origin = {-12, 6}, points = {{0, 46}, {0, -46}}), Line(origin = {30, 9}, points = {{0, 49}, {0, -49}}), Line(origin = {44, 10}, points = {{0, 50}, {0, -50}}), Line(origin = {76, 12}, points = {{0, 52}, {0, -52}}), Line(origin = {88, 13}, points = {{0, 53}, {0, -53}})}),
  Diagram(graphics = {Text(origin = {-77, 0}, extent = {{-17, 10}, {17, -10}}, textString = "N elements")}));
 end Piping_local;
--- a/banana_simple.mo
+++ b/banana_simple.mo
-model banana_simple
+model banana_local

  package Medium =
      Modelica.Media.Water.ConstantPropertyLiquidWater;
  
-  Piping_simple piping(channel_number = 22)  annotation(
+  Piping_local piping(channel_number = 22)  annotation(
    Placement(visible = true, transformation(origin = {-6, -4}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
  inner Modelica.Fluid.System system(m_flow_start = 0.008)  annotation(
    Placement(visible = true, transformation(origin = {-70, 70}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
@@ -30,5 +30,5 @@ equation
    Line(points = {{2, -14}, {2, -34}, {-34, -34}}, color = {0, 127, 255}));
  annotation(
    uses(Modelica(version = "3.2.3")),
-    Diagram(graphics = {Text(origin = {41, 70}, extent = {{-39, 12}, {39, -12}}, textString = "Piping_simple")}));
-end banana_simple;
+    Diagram(graphics = {Text(origin = {41, 70}, extent = {{-39, 12}, {39, -12}}, textString = "Piping_local")}));
+end banana_local;
--- a/test_novec.mo
+++ b/test_novec.mo
+model test_novec
+  package Medium = Novec649;
+  Medium.BaseProperties medium;
+
+  Medium.DynamicViscosity eta=Medium.dynamicViscosity(medium.state);
+  Medium.ThermalConductivity lambda=Medium.thermalConductivity(medium.state);
+  Medium.SpecificEntropy s=Medium.specificEntropy(medium.state);
+  Medium.SpecificHeatCapacity cv=Medium.specificHeatCapacityCv(medium.state);
+  Medium.SpecificInternalEnergy u=Medium.specificInternalEnergy(medium.state);
+  Medium.SpecificInternalEnergy h=Medium.specificEnthalpy(medium.state);
+  Medium.Density d=Medium.density(medium.state);
+  protected
+  constant Modelica.SIunits.Time timeUnit = 1;
+  constant Modelica.SIunits.Temperature Ta = 1;
+equation
+  medium.p = 1.013e5;
+  medium.T = Medium.T_min + time/timeUnit*Ta;
+    annotation (experiment(StopTime=1.01));
+
+
+
+end test_novec;
--- a/waffle.mo
+++ b/waffle.mo
--- a/waffle_validation.mo
+++ b/waffle_validation.mo
+model waffle_validation
+
+  package Medium =
+      Novec649;
+  inner Modelica.Fluid.System system(m_flow_start = 0.008)  annotation(
+    Placement(visible = true, transformation(origin = {-70, 70}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+  Modelica.Fluid.Sources.FixedBoundary boundary(redeclare package Medium = Medium, T = 273.15 - 40, nPorts = 1, p = 2e5)  annotation(
+    Placement(visible = true, transformation(origin = {-44, 32}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+  
+  Modelica.Fluid.Sources.FixedBoundary boundary1(redeclare package Medium = Medium, nPorts = 1, p = 1e5) annotation(
+    Placement(visible = true, transformation(origin = {-44, -34}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+  waffle waffle_test(redeclare package Medium = Medium, channel_cross_section_height = 3e-3, channel_cross_section_width = 3e-3) annotation(
+    Placement(visible = true, transformation(origin = {0, 4}, extent = {{-10, -10}, {10, 10}}, rotation = 0)));
+equation
+  connect(boundary.ports[1], waffle_test.port_a) annotation(
+    Line(points = {{-34, 32}, {0, 32}, {0, 14}}, color = {0, 127, 255}));
+  connect(waffle_test.port_b, boundary1.ports[1]) annotation(
+    Line(points = {{0, -6}, {0, -34}, {-34, -34}}, color = {0, 127, 255}));
+  annotation(
+    uses(Modelica(version = "3.2.3")),
+    Diagram);
+end waffle_validation;