Changeset 147 for branches/newlanguage/eml/streams.mso

Timestamp:

Jan 31, 2007, 5:02:46 PM (17 years ago)

Author:

gerson bicca

Message:

testing the new language

File:

• branches/newlanguage/eml/streams.mso (modified) (4 diffs)

branches/newlanguage/eml/streams.mso

@@ -107 +107 @@
 
         PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        outer NComp as Integer (Brief = "Number of chemical components", Lower = 1); 
-        
+        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
+        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1); 
+                        M(NComp)  as molweight  (Brief="Component Mol Weight");
+        
+        SET
+
+        M   = PP.MolecularWeight();
+
         VARIABLES
         out Outlet as stream;
@@ -116 +121 @@
         hl as enth_mol;
         hv as enth_mol;
+        zmass(NComp)    as fraction                     (Brief = "Mass Fraction");
+        Mw                                      as molweight            (Brief="Average Mol Weight");
+        vm                                      as volume_mol   (Brief="Molar Volume"); 
+        rho                                     as dens_mass            (Brief="Stream Density");
+        Fw                                      as flow_mass            (Brief="Stream Mass Flow");
+        Fvol                            as flow_vol         (Brief = "Volumetric Flow");
         
         EQUATIONS
         "Flash Calculation"
         [Outlet.v, x, y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
+        
         "Overall Enthalpy"
         Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, x) +
@@ -126 +138 @@
         hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, x);
         hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
+        
+        "Average Molecular Weight"
+        Mw = sum(M*Outlet.z);
+
+        "Mass Density"
+        rho =   (1-Outlet.v)*PP.LiquidDensity(Outlet.T,Outlet.P,x) + Outlet.v*PP.VapourDensity(Outlet.T,Outlet.P,y);
+
+        "Flow Mass"
+        Fw      =  Mw*Outlet.F;
+
+        "Molar Volume"
+        vm = (1-Outlet.v)*PP.LiquidVolume(Outlet.T, Outlet.P, x) + Outlet.v*PP.VapourVolume(Outlet.T,Outlet.P,y);
+        
+        "Volumetric Flow"
+        Fvol = Outlet.F*vm ;
+        
+"Mass Fraction"
+        zmass = M*Outlet.z / Mw;        
+        
 end
 
@@ -135 +166 @@
 
         PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        outer NComp as Integer (Brief = "Number of chemical components", Lower = 1); 
+        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
+        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1); 
+                        M(NComp)  as molweight  (Brief="Component Mol Weight");
+        
+        SET
+
+        M   = PP.MolecularWeight();
         
         VARIABLES
         in Inlet as stream;
         v as fraction;
-        x(NComp) as fraction(Brief = "Liquid Molar Fraction");
-        y(NComp) as fraction(Brief = "Vapour Molar Fraction");
+        x(NComp)        as fraction                     (Brief = "Liquid Molar Fraction");
+        y(NComp)        as fraction                     (Brief = "Vapour Molar Fraction");
+        zmass(NComp)    as fraction     (Brief = "Mass Fraction");
+        Mw                      as molweight            (Brief="Average Mol Weight");
+        vm                      as volume_mol   (Brief="Molar Volume"); 
+        rho                     as dens_mass            (Brief="Stream Density");
+        Fw                      as flow_mass            (Brief="Stream Mass Flow");
+        Fvol          as flow_vol          (Brief = "Volumetric Flow");
         
         EQUATIONS
         "Flash Calculation"
         [v, x, y] = PP.FlashPH(Inlet.P, Inlet.h, Inlet.z);
-end
+        
+        "Average Molecular Weight"
+        Mw = sum(M*Inlet.z);
+
+        "Mass Density"
+        rho =   (1-v)*PP.LiquidDensity(Inlet.T,Inlet.P,x) + v*PP.VapourDensity(Inlet.T,Inlet.P,y);
+
+        "Flow Mass"
+        Fw      =  Mw*Inlet.F;
+
+        "Molar Volume"
+        vm = (1-v)*PP.LiquidVolume(Inlet.T, Inlet.P, x) + v*PP.VapourVolume(Inlet.T,Inlet.P,y);
+        
+        "Volumetric Flow"
+        Fvol = Inlet.F*vm ;
+        
+        "Mass Fraction"
+        zmass = M*Inlet.z / Mw;
+        
+end