Changeset 339 for branches/tests/eml/streams.mso

Timestamp:

Aug 14, 2007, 2:46:50 PM (16 years ago)

Author:

gerson bicca

Message:

testing a new streamTherm

File:

• branches/tests/eml/streams.mso (modified) (14 diffs)

branches/tests/eml/streams.mso

@@ -35 +35 @@
 
         VARIABLES
-        F as flow_mol;
-        T as temperature;
-        P as pressure;
+        F as flow_mol(Brief = "Molar Flow");
+        T as temperature(Brief = "Temperature");
+        P as pressure(Brief = "Pressure");
         z(NComp) as fraction(Brief = "Overall Molar Fraction");
-        h as enth_mol;
+        h as enth_mol(Brief = "Molar Enthalpy");
         v as fraction(Brief = "Vapourization fraction");
+end
+
+Model streamTherm as stream
+        ATTRIBUTES
+        Pallete = false;
+        Brief = "General Material Stream";
+        Info =
+        "comments.";
+        
+        PARAMETERS
+        outer NComp as Integer (Brief = "Number of chemical components", Lower = 1); 
+
+        VARIABLES
+        x(NComp) as fraction(Brief = "Liquid Molar Fraction");
+        y(NComp) as fraction(Brief = "Vapour Molar Fraction");
+
 end
 
@@ -60 +76 @@
         "Liquid stream"
         v = 0;
+
 end
 
@@ -79 +96 @@
         "Vapour stream"
         v = 1;
-end
-
-Model streamPH as stream
+
+end
+
+Model streamPH as streamTherm
         ATTRIBUTES
         Brief = "Stream with built-in flash calculation";
@@ -99 +117 @@
         PARAMETERS
         outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
-        
-        VARIABLES
-        x(NComp) as fraction(Brief = "Liquid Molar Fraction");
-        y(NComp) as fraction(Brief = "Vapour Molar Fraction");
         
         EQUATIONS
@@ -150 +164 @@
 
         VARIABLES
-        out Outlet                      as stream                       (Brief = "Outlet stream", PosX=1, PosY=0.5256);
-        x(NComp)                        as fraction                     (Brief = "Liquid Molar Fraction");
-        y(NComp)                        as fraction                     (Brief = "Vapour Molar Fraction");
+        out Outlet                      as streamTherm                  (Brief = "Outlet stream", PosX=1, PosY=0.5256);
         hl                                      as enth_mol;
         hv                                      as enth_mol;
@@ -165 +177 @@
         EQUATIONS
         "Flash Calculation"
-        [Outlet.v, x, y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
+        [Outlet.v, Outlet.x, Outlet.y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
         
         "Overall Enthalpy"
-        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, x) +
-                Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
-        
-        hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, x);
-        hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
+        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, Outlet.x) +
+                Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, Outlet.y);
+        
+        hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, Outlet.x);
+        hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, Outlet.y);
         
         "Average Molecular Weight"
@@ -184 +196 @@
                 case "correlation":
         "Mass Density"
-                rho*((1-Outlet.v)/PP.LiquidDensity(Outlet.T,Outlet.P,x) + Outlet.v/PP.VapourDensity(Outlet.T,Outlet.P,y)) = 1;
+                rho*((1-Outlet.v)/PP.LiquidDensity(Outlet.T,Outlet.P,Outlet.x) + Outlet.v/PP.VapourDensity(Outlet.T,Outlet.P,Outlet.y)) = 1;
         end
         
@@ -194 +206 @@
 
         "Molar Volume"
-        vm = (1-Outlet.v)*PP.LiquidVolume(Outlet.T, Outlet.P, x) + Outlet.v*PP.VapourVolume(Outlet.T,Outlet.P,y);
+        vm = (1-Outlet.v)*PP.LiquidVolume(Outlet.T, Outlet.P, Outlet.x) + Outlet.v*PP.VapourVolume(Outlet.T,Outlet.P,Outlet.y);
         
         "Volumetric Flow"
@@ -226 +238 @@
         
         VARIABLES
-        out Outlet                      as stream               (Brief = "Outlet stream", PosX=1, PosY=0.5256);
-        x(NComp)                        as fraction             (Brief = "Liquid Molar Fraction");
-        y(NComp)                        as fraction             (Brief = "Vapour Molar Fraction");
+        out Outlet                      as streamTherm          (Brief = "Outlet stream", PosX=1, PosY=0.5256);
         hl                                      as enth_mol;
         hv                                      as enth_mol;
@@ -234 +244 @@
         EQUATIONS
         "Flash Calculation"
-        [Outlet.v, x, y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
+        [Outlet.v, Outlet.x, Outlet.y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
         
         "Overall Enthalpy"
-        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, x) +
-                Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
-
-        hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, x);
-        hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
+        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, Outlet.x) +
+                Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, Outlet.y);
+
+        hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, Outlet.x);
+        hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, Outlet.y);
 end
 
@@ -274 +284 @@
         
         VARIABLES
-        in Inlet                as stream               (Brief = "Inlet Stream", PosX=0, PosY=0.5308);
-        v                               as fraction             (Brief = "Vapourization fraction");
-        x(NComp)                as fraction             (Brief = "Liquid Molar Fraction");
-        y(NComp)                as fraction             (Brief = "Vapour Molar Fraction");
+        in Inlet                as streamTherm          (Brief = "Inlet Stream", PosX=0, PosY=0.5308);
         zmass(NComp)    as fraction             (Brief = "Mass Fraction");
         Mw                              as molweight    (Brief = "Average Mol Weight");
@@ -287 +294 @@
         
         EQUATIONS
-        "Flash Calculation"
-        [v, x, y] = PP.FlashPH(Inlet.P, Inlet.h, Inlet.z);
-        
         "Average Molecular Weight"
         Mw = sum(M*Inlet.z);
@@ -300 +304 @@
                 case "correlation":
         "Mass Density"
-                rho * ((1-v)/PP.LiquidDensity(Inlet.T,Inlet.P,x) + v/PP.VapourDensity(Inlet.T,Inlet.P,y)) = 1;
+                rho * ((1-Inlet.v)/PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.x) + Inlet.v/PP.VapourDensity(Inlet.T,Inlet.P,Inlet.y)) = 1;
         end
         
@@ -310 +314 @@
 
         "Molar Volume"
-        vm = (1-v)*PP.LiquidVolume(Inlet.T, Inlet.P, x) + v*PP.VapourVolume(Inlet.T,Inlet.P,y);
+        vm = (1-Inlet.v)*PP.LiquidVolume(Inlet.T, Inlet.P, Inlet.x) + Inlet.v*PP.VapourVolume(Inlet.T,Inlet.P,Inlet.y);
         
         "Volumetric Flow"