ideal gas limit, related to #11

thorade · Jan 28, 2013 · 610bea8 · 610bea8
1 parent f9554f1
commit 610bea8
Show file tree

Hide file tree

Showing 2 changed files with 36 additions and 0 deletions.
diff --git a/Examples/Validation/idealGasLimit.mo b/Examples/Validation/idealGasLimit.mo
@@ -0,0 +1,35 @@
+within HelmholtzMedia.Examples.Validation;
+model idealGasLimit
+  package Medium = HelmholtzFluids.Helium;
+
+  parameter Medium.Temperature T=298.15;
+  Medium.Density d;
+
+  Medium.EoS.HelmholtzDerivs f = Medium.EoS.setHelmholtzDerivsThird(d=d, T=T, phase=1);
+  Medium.ThermodynamicState state = Medium.setState_dT(d=d, T=T, phase=1);
+
+equation
+  d = 100/exp(time);
+
+  annotation (Documentation(info="
+<html>
+<style type=\"text/css\">
+  code{background-color:#EEEEEE; padding:2px; margin:2px;}
+</style>
+<body>
+This model is used to validate the behavior of the EoS in the ideal gas limit, that is at <code> lim d -> 0 </code>.
+<br />
+How to use:
+<ol>
+  <li>Simulate with any <code>T_trip &lt; T &lt; T_max</code></li>
+  <li>Plot <code>f.rd</code> and <code>B</code> </li>
+  <li>Make <code>d</code> the independent variable of the plot </li>
+  <li>Set the horizontal axis <code>(d)</code> to logarithmic scale</li>
+  <li>Set line style to dotted and marker style to cross</li>
+  <li>Look at the graphs </li>
+</ol>
+
+</body>
+</html>"),experiment(StopTime=1000, __Dymola_NumberOfIntervals=10000),
+    __Dymola_experimentSetupOutput);
+end idealGasLimit;
diff --git a/Examples/Validation/package.order b/Examples/Validation/package.order
@@ -8,3 +8,4 @@ CriticalState
 acentricFactor
 ThermoSurface
 MaxwellLoop
+idealGasLimit