Changeset 1018 for trunk

Timestamp:

01/29/08 12:57:27 (7 months ago)

Author:

hubertus

Message:

- Fixed usage of sqrt() to take into account new Dymola semantics
- Added min=0 to a few more units in SIunits in order to guarantee positive argument to some calls to sqrt()

Location:

trunk/Modelica

Files:

• Media/IdealGases/Common/package.mo (modified) (3 diffs)
• Media/Water/IF97_Utilities.mo (modified) (4 diffs)
• Media/package.mo (modified) (1 diff)
• SIunits.mo (modified) (1 diff)

trunk/Modelica/Media/IdealGases/Common/package.mo

@@ -353 +353 @@
     extends Modelica.Icons.Function;
   algorithm
-    a := sqrt(data.R*state.T*cp_T(data, state.T)/specificHeatCapacityCv(state));
+    a := sqrt(max(0,data.R*state.T*cp_T(data, state.T)/specificHeatCapacityCv(state)));
   end velocityOfSound;
 
@@ -598 +598 @@
     parameter Real Const2_SI=131.3/1000.0
       "Constant in formula for mur converted to SI units";
-    parameter Real mur=Const2_SI*mu/sqrt(Vc*Tc)
+    Real mur=Const2_SI*mu/sqrt(Vc*Tc)
       "Dimensionless dipole moment of gas molecule";
-    parameter Real Fc=1 - 0.2756*w + 0.059035*mur^4 + k
+    Real Fc=1 - 0.2756*w + 0.059035*mur^4 + k
       "Factor to account for molecular shape and polarities of gas";
     Real Tstar "Dimensionless temperature defined by equation below";
@@ -1025 +1025 @@
     input ThermodynamicState state "properties at upstream location";
   algorithm
-    a := sqrt(gasConstant(state)*state.T*specificHeatCapacityCp(state)/specificHeatCapacityCv(state));
+    a := sqrt(max(0,gasConstant(state)*state.T*specificHeatCapacityCp(state)/specificHeatCapacityCv(state)));
   end velocityOfSound;
 

trunk/Modelica/Media/Water/IF97_Utilities.mo

@@ -3416 +3416 @@
         Real[49] o "vector of auxiliary variables";
       algorithm
-        pi := p*data.IPSTAR;
+        pi := max(Modelica.Constants.small,p*data.IPSTAR);
         o[1] := pi^0.75;
         o[2] := 1/o[1];
@@ -6818 +6818 @@
   algorithm
     // dp/drho at constant s
-    v_sound := if aux.region == 3 then sqrt((aux.pd*aux.rho*aux.rho*aux.cv + aux.pt*aux.pt*aux.T)/(aux.rho*aux.rho*aux.cv)) else 
+    v_sound := if aux.region == 3 then sqrt(max(0,(aux.pd*aux.rho*aux.rho*aux.cv + aux.pt*aux.pt*aux.T)/(aux.rho*aux.rho*aux.cv))) else 
       if aux.region == 4 then 
-      sqrt(1/((aux.rho*(aux.rho*aux.cv/aux.dpT + 1.0)/(aux.dpT*aux.T)) - 1/aux.rho*aux.rho*aux.rho/(aux.dpT*aux.T))) else 
-           sqrt(-aux.cp/(aux.rho*aux.rho*(aux.vp*aux.cp+aux.vt*aux.vt*aux.T)));
+      sqrt(max(0,1/((aux.rho*(aux.rho*aux.cv/aux.dpT + 1.0)/(aux.dpT*aux.T)) - 1/aux.rho*aux.rho*aux.rho/(aux.dpT*aux.T)))) else 
+           sqrt(max(0,-aux.cp/(aux.rho*aux.rho*(aux.vp*aux.cp+aux.vt*aux.vt*aux.T))));
   end velocityOfSound_props_ph;
 
@@ -7234 +7234 @@
   algorithm
     // dp/drho at constant s
-    v_sound := if aux.region == 3 then sqrt((aux.pd*aux.rho*aux.rho*aux.cv + aux.pt*aux.pt*aux.T)/(aux.rho*aux.rho*aux.cv)) else 
-      sqrt(-aux.cp/(aux.rho*aux.rho*(aux.vp*aux.cp+aux.vt*aux.vt*aux.T)));
+    v_sound := if aux.region == 3 then sqrt(max(0,(aux.pd*aux.rho*aux.rho*aux.cv + aux.pt*aux.pt*aux.T)/(aux.rho*aux.rho*aux.cv))) else 
+      sqrt(max(0,-aux.cp/(aux.rho*aux.rho*(aux.vp*aux.cp+aux.vt*aux.vt*aux.T))));
   end velocityOfSound_props_pT;
 
@@ -7640 +7640 @@
   algorithm
     // dp/drho at constant s
-    v_sound := if aux.region == 3 then sqrt((aux.pd*aux.rho*aux.rho*aux.cv + aux.pt*aux.pt*aux.T)/(aux.rho*aux.rho*aux.cv)) else 
+    v_sound := if aux.region == 3 then sqrt(max(0,((aux.pd*aux.rho*aux.rho*aux.cv + aux.pt*aux.pt*aux.T)/(aux.rho*aux.rho*aux.cv)))) else 
       if aux.region == 4 then 
-      sqrt(1/((aux.rho*(aux.rho*aux.cv/aux.dpT + 1.0)/(aux.dpT*aux.T)) - 1/aux.rho*aux.rho*aux.rho/(aux.dpT*aux.T))) else 
-           sqrt(-aux.cp/(aux.rho*aux.rho*(aux.vp*aux.cp+aux.vt*aux.vt*aux.T)));
+      sqrt(max(0,1/((aux.rho*(aux.rho*aux.cv/aux.dpT + 1.0)/(aux.dpT*aux.T)) - 1/aux.rho*aux.rho*aux.rho/(aux.dpT*aux.T)))) else 
+           sqrt(max(0,-aux.cp/(aux.rho*aux.rho*(aux.vp*aux.cp+aux.vt*aux.vt*aux.T))));
   end velocityOfSound_props_dT;
 

trunk/Modelica/Media/package.mo

@@ -4963 +4963 @@
         extends Modelica.Icons.Function;
       algorithm
-        a := sqrt(1/(kappa_const*density(state) -beta_const*beta_const*state.T/cp_const));
+        a := sqrt(max(0,1/(kappa_const*density(state) -beta_const*beta_const*state.T/cp_const)));
       end velocityOfSound;
 

trunk/Modelica/SIunits.mo

@@ -1381 +1381 @@
       final unit="mol",
       min=0);
-  type MolarMass = Real (final quantity="MolarMass", final unit="kg/mol");
-  type MolarVolume = Real (final quantity="MolarVolume", final unit="m3/mol");
+  type MolarMass = Real (final quantity="MolarMass", final unit="kg/mol",min=0);
+  type MolarVolume = Real (final quantity="MolarVolume", final unit="m3/mol", min=0);
   type MolarInternalEnergy = Real (final quantity="MolarInternalEnergy", final unit
         =    "J/mol");