Changeset 175 for branches/newlanguage/eml

Timestamp:

Mar 2, 2007, 1:57:58 PM (17 years ago)

Author:

gerson bicca

Message:

updated sepComp model and samples of mixers/splitter

File:

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

branches/newlanguage/eml/mixers_splitters/sepComp.mso

@@ -52 +52 @@
 Model sepComp_n
         
+PARAMETERS
 
-        PARAMETERS
-ext PP   as CalcObject (Brief = "External Physical Properties");
-ext     NComp as Integer (Brief = "Number of chemical components", Lower = 1);
-        NOutlet as Integer (Brief = "Number of Outlet Streams", Lower = 1);
-        mainComp as Integer (Brief = "Component specified", Default = 1, Lower = 1);
+outer PP                        as Plugin       (Brief = "External Physical Properties");
+outer   NComp           as Integer      (Brief = "Number of chemical components", Lower = 1);
+                NOutlet                 as Integer      (Brief = "Number of Outlet Streams", Lower = 1);
+                mainComp        as Integer      (Brief = "Component specified", Default = 1, Lower = 1);
         
-        VARIABLES
-in      Inlet   as stream;
-out Outlet(NOutlet) as stream;
-        frac(NOutlet) as fraction (Brief = "Distribution of the Outlet streams");
-        recovery(NOutlet) as fraction (Brief = "Recovery of the component specified");
+VARIABLES
 
-        EQUATIONS
+in              Inlet                                   as stream;
+out     Outlet(NOutlet)         as stream;
+
+        frac(NOutlet)                   as fraction (Brief = "Distribution of the Outlet streams");
+        recovery(NOutlet)       as fraction (Brief = "Recovery of the component specified");
+
+EQUATIONS
         
-        "Flow"
+"Flow"
         sum(Outlet.F) = Inlet.F;
         
         
-        for i in [1:NOutlet-1]
-                
-                "Mol fraction normalisation"
-                sum(Outlet(i).z) = 1;
-        end
+for i in [1:NOutlet-1]
+
+"Mol fraction normalisation"
+        sum(Outlet(i).z) = 1;
+
+end
         
         
-        for i in [1:NComp]
+for i in [1:NComp]
         
-        "Composition"
-                sum(Outlet.F*Outlet.z(i)) = Inlet.F*Inlet.z(i);
+"Composition"
+        sum(Outlet.F*Outlet.z(i)) = Inlet.F*Inlet.z(i);
         
-        end     
+end     
         
         
-        for i in [1:NOutlet]
-                
-                "Flow"
-                Outlet(i).F = Inlet.F*frac(i);
-                
-                "Recovery"
-                recovery(i)*Inlet.z(mainComp) = frac(i)*Outlet(i).z(mainComp);
+for i in [1:NOutlet]
         
-                "Pressure"
-                Outlet(i).P = Inlet.P;
+"Flow"
+        Outlet(i).F = Inlet.F*frac(i);
         
-                "Enthalpy"
-                Outlet(i).h = (1-Outlet(i).v)*PP.LiquidEnthalpy(Outlet(i).T, Outlet(i).P, Outlet(i).z) + 
+"Recovery"
+        recovery(i)*Inlet.z(mainComp) = frac(i)*Outlet(i).z(mainComp);
+        
+"Pressure"
+        Outlet(i).P = Inlet.P;
+        
+"Enthalpy"
+        Outlet(i).h = (1-Outlet(i).v)*PP.LiquidEnthalpy(Outlet(i).T, Outlet(i).P, Outlet(i).z) + 
                                 Outlet(i).v*PP.VapourEnthalpy(Outlet(i).T, Outlet(i).P, Outlet(i).z);
         
-                "Temperature"   
-                Outlet(i).T = Inlet.T;
-         
-                "Vapourization Fraction"
-                Outlet(i).v = PP.VapourFraction(Outlet(i).T, Outlet(i).P, Outlet(i).z);
-        end
+"Temperature"   
+        Outlet(i).T = Inlet.T;
+
+"Vapourization Fraction"
+        Outlet(i).v = PP.VapourFraction(Outlet(i).T, Outlet(i).P, Outlet(i).z);
+        
+end
+
 end
 
@@ -112 +117 @@
 Model sepComp
         
-        PARAMETERS
-ext PP   as CalcObject (Brief = "External Physical Properties");
-ext     NComp as Integer (Brief = "Number of chemical components", Lower = 1); 
-        mainComp as Integer (Brief = "Component specified", Default = 1, Lower = 1);
+PARAMETERS
+
+outer PP                        as Plugin       (Brief = "External Physical Properties");
+outer   NComp           as Integer      (Brief = "Number of chemical components", Lower = 1); 
+                mainComp        as Integer      (Brief = "Component specified", Default = 1, Lower = 1);
         
-        VARIABLES
-in      Inlet   as stream;
-out Outlet1 as stream;
-out Outlet2 as stream;
-        frac as fraction (Brief = "Fraction to Outlet 1");
-        recovery as fraction (Brief = "Recovery of the component specified");
+VARIABLES
+
+in              Inlet           as stream;
+out     Outlet1         as stream;
+out     Outlet2         as stream;
+        
+        frac                    as fraction (Brief = "Fraction to Outlet 1");
+        recovery        as fraction (Brief = "Recovery of the component specified");
 
         EQUATIONS
         
-        "Flow"
+"Flow"
         Outlet1.F = Inlet.F * frac;
         Outlet1.F + Outlet2.F = Inlet.F;
@@ -134 +142 @@
         sum(Outlet1.z) = 1;
         
-        for i in [1:NComp]
-                
-        "Composition"
-                Outlet1.F*Outlet1.z(i) + Outlet2.F*Outlet2.z(i) = Inlet.F*Inlet.z(i);
-        end
+for i in [1:NComp]
         
-        "Pressure"
+ "Composition"
+        Outlet1.F*Outlet1.z(i) + Outlet2.F*Outlet2.z(i) = Inlet.F*Inlet.z(i);
+
+end
+        
+"Pressure"
         Outlet1.P = Inlet.P;
         Outlet2.P = Inlet.P;
         
-        "Enthalpy"
+"Enthalpy"
         Outlet1.h = (1-Outlet1.v)*PP.LiquidEnthalpy(Outlet1.T, Outlet1.P, Outlet1.z) + 
                                 Outlet1.v*PP.VapourEnthalpy(Outlet1.T, Outlet1.P, Outlet1.z);
@@ -150 +159 @@
                                 Outlet2.v*PP.VapourEnthalpy(Outlet2.T, Outlet2.P, Outlet2.z);
         
-        "Temperature"
+"Temperature"
         Outlet1.T = Inlet.T;
         Outlet2.T = Inlet.T;
         
-        "Vapourization Fraction"
+"Vapourization Fraction"
         Outlet1.v = PP.VapourFraction(Outlet1.T, Outlet1.P, Outlet1.z);
         Outlet2.v = PP.VapourFraction(Outlet2.T, Outlet2.P, Outlet2.z);
+
 end