Changeset 210 for branches/newlanguage/eml/pressure_changers/turbine.mso

Timestamp:

Mar 15, 2007, 9:52:28 AM (16 years ago)

Author:

Argimiro Resende Secchi

Message:

Remove convergence problems of some samples.

File:

• branches/newlanguage/eml/pressure_changers/turbine.mso (modified) (5 diffs)

branches/newlanguage/eml/pressure_changers/turbine.mso

@@ -37 +37 @@
 outer PP                as Plugin                       (Brief = "External Physical Properties", Type="PP");
         Mw(NComp)       as molweight            (Brief = "Molar Weight");
-        Eff     as positive                     (Default = 0.72, Brief = "Pump efficiency");
-        Meff    as positive                     (Default = 1.0, Brief = "Brake efficiency");
-        Beta    as positive                     (Default = 0, Brief = "Volumetric expansivity", Unit = '1/K');
-        
         
         VARIABLES
+        Eff     as efficiency           (Brief = "Pump efficiency");
+        Meff    as efficiency           (Brief = "Brake efficiency");
+        Beta    as positive             (Brief = "Volumetric expansivity", Unit = '1/K');
         Head    as head                         (Brief = "Head Developed");
         FPower  as power                        (Brief = "Fluid Power");
@@ -59 +58 @@
         #Mixtures Properties
         "Calculate Mwm for Inlet Mixture"
-        Mwm = sum(Mw([1:NComp])*Inlet.z([1:NComp]));
-        
+        Mwm = sum(Mw*Inlet.z);
+
         "Calculate rho using a External Physical Properties Routine"
         rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
@@ -78 +77 @@
         
         "Calculate Fluid Power"
-        FPower = Pdiff * Inlet.F / (rho/Mwm);
+        FPower * rho = Pdiff * Inlet.F * Mwm;
         
         "Calculate Brake Power"
@@ -87 +86 @@
         
         "Calculate Outlet Temperature"
-        (Outlet.T - Inlet.T)*Cp = (Outlet.h - Inlet.h) -  Pdiff /(rho/Mwm) * (1-Beta*Inlet.T);
+        (Outlet.T - Inlet.T) * rho * Cp = (Outlet.h - Inlet.h) * rho
+        -  Pdiff * Mwm * (1-Beta*Inlet.T);
         
         "Calculate Outlet Enthalpy"
-        Outlet.h - Inlet.h =  Pdiff / (rho/Mwm);
+        (Outlet.h - Inlet.h) * rho =  Pdiff * Mwm;
         
         "Molar Balance"
@@ -99 +99 @@
 
         "Calculate Head"
-        Head = Pdiff/rho;
-
+        Head * rho = Pdiff;
 end