Changeset 398 for trunk/eml/stage_separators/tray.mso

Timestamp:

Oct 26, 2007, 12:44:56 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

Updated tray and column models and samples

File:

• trunk/eml/stage_separators/tray.mso (modified) (2 diffs)

trunk/eml/stage_separators/tray.mso

@@ -131 +131 @@
         alfa as fraction (Brief="Dry pressure drop coefficient");
         
+        VapourFlow as Switcher(Valid = ["on", "off"], Default = "on");
+        LiquidFlow as Switcher(Valid = ["on", "off"], Default = "on");
+        
         VARIABLES
         rhoL as dens_mass;
@@ -141 +144 @@
         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
 
-        if Level > (beta * hw) then
+        switch LiquidFlow
+                case "on":
                 "Francis Equation"
-                OutletL.F = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2/vL;
-        else
+#               OutletL.F*vL = 1.84*'m^0.5/s'*lw*((Level-(beta*hw))/(beta))^1.5;
+                OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2;
+                when Level < (beta * hw) switchto "off";
+                
+                case "off":
                 "Low level"
                 OutletL.F = 0 * 'mol/h';
+                when Level > (beta * hw) + 1e-6*'m' switchto "on";
+        end
+
+        switch VapourFlow
+                case "on":
+                InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;
+                when InletV.F < 1e-6 * 'kmol/h' switchto "off";
+                
+                case "off":
+                InletV.F = 0 * 'mol/s';
+                when InletV.P > OutletV.P + Level*g*rhoL + 1e-1 * 'atm' switchto "on";
         end