Changeset 235 for trunk

Timestamp:

Apr 11, 2007, 10:16:05 PM (16 years ago)

Author:

Argimiro Resende Secchi

Message:

Added orientation type (vertical and horizontal) into the flash model. "Across" is now a variable and "diameter" is a new parameter.

Location:

trunk

Files:

• eml/stage_separators/flash.mso (modified) (3 diffs)
• sample/stage_separators/sample_flash.mso (modified) (2 diffs)

trunk/eml/stage_separators/flash.mso

@@ -41 +41 @@
 Model flash
         PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
+        outer PP as Plugin (Brief = "External Physical Properties", Type="PP");
         outer NComp as Integer (Brief = "Number of chemical components", Lower = 1); 
-        V as volume(Brief="Total Volume of the flash");
+        V as volume (Brief="Total Volume of the flash");
         Mw(NComp) as molweight;
-        Across as area  (Brief="Flash Cross section area");
-        
+        orientation as Switcher (Valid=["vertical","horizontal"],Default="vertical");
+        diameter as length (Brief="Vessel diameter");
+
         SET
         Mw=PP.MolecularWeight();
-        
+
         VARIABLES
         in      Inlet as stream(Brief="Feed Stream");
@@ -63 +64 @@
         vV as volume_mol (Brief="Vapour Molar volume");
         Level as length (Brief="liquid height");
-        
+        Across as area (Brief="Flash Cross section area");
+
         EQUATIONS
         "Component Molar Balance"
@@ -98 +100 @@
         
         "Geometry Constraint"
-        V = ML* vL + MV*vV;
+        V = ML * vL + MV * vV;
+
+        switch orientation
+        case "vertical":
+        "Cross Section Area"
+                Across = 0.5 * asin(1) * diameter^2;
+        
         "Liquid Level"
-        ML* vL = Across * Level;
+                ML * vL = Across * Level;
+
+        case "horizontal":
+        "Cylindrical Side Area"
+                Across = 0.25*diameter^2 * (asin(1) - asin((diameter - 2*Level)/diameter)) + 
+                                (Level - 0.5*diameter)*sqrt(Level*(diameter - Level));
+
+        "Liquid Level"
+                0.5 * asin(1) * diameter^2 * ML* vL = Across * V;
+        end
 end
 

trunk/sample/stage_separators/sample_flash.mso

@@ -65 +65 @@
         
         SET
-        fl.V = 2000 * 'm^3';
-        fl.Across = 1 * 'm^2';
-        
+        fl.V = 50 * 'm^3';
+        fl.diameter = 2 * 'm';
+        fl.orientation = "vertical";
+#       fl.orientation = "horizontal";
+
         INITIAL
         fl.OutletL.T = 338 * 'K';
-        fl.Level = 1 * 'm';
+        fl.Level = 0.4 * 'm';
 
         fl.OutletL.z(1) = 0.1;
@@ -83 +85 @@
         Dynamic = true;
         TimeStep = 0.1;
-        TimeEnd = 20;
+        TimeEnd = 5;
         TimeUnit = 'h';
 end