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

Timestamp:

Apr 16, 2007, 2:19:39 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

Updated reactive distillation models andsample

File:

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

trunk/eml/stage_separators/tray.mso

@@ -159 +159 @@
 
         PARAMETERS
-        outer PP as Plugin(Brief = "External Physical Properties", Type="PP");
+        outer PP as Plugin(Type="PP");
         outer NComp as Integer;
         V as volume(Brief="Total Volume of the tray");
@@ -175 +175 @@
         Hr as energy_mol;
         Pstartup as pressure;
+        
+        VapourFlow as Switcher(Valid = ["on", "off"], Default = "off");
+        LiquidFlow as Switcher(Valid = ["on", "off"], Default = "off");
         
         VARIABLES
@@ -197 +200 @@
         rhoL as dens_mass;
         rhoV as dens_mass;
-        r as reaction_mol (Brief = "Specific reaction rate");
+        r3 as reaction_mol (Brief = "Reaction resulting ethyl acetate", DisplayUnit = 'mol/l/s');
         C(NComp) as conc_mol (Brief = "Molar concentration", Lower = -1); #, Unit = "mol/l");
         
@@ -204 +207 @@
         OutletL.z = vL * C;
         
+        "Reaction"
+        r3 = exp(-7150*'K'/OutletL.T)*(4.85e4*C(1)*C(2) - 1.23e4*C(3)*C(4))*'l/mol/s';
+        
         "Component Molar Balance"
         diff(M)=Inlet.F*Inlet.z + InletL.F*InletL.z + InletV.F*InletV.z
-                - OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r*ML*vL;
+                - OutletL.F*OutletL.z - OutletV.F*OutletV.z + stoic*r3*ML*vL;
         
         "Energy Balance"
         diff(E) = ( Inlet.F*Inlet.h + InletL.F*InletL.h + InletV.F*InletV.h
-                - OutletL.F*OutletL.h - OutletV.F*OutletV.h + Q ) + Hr * r * vL*ML;
+                - OutletL.F*OutletL.h - OutletV.F*OutletV.h + Q ) + Hr * r3 * vL*ML;
         
         "Molar Holdup"
@@ -242 +248 @@
         rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z);
 
-        if Level > (beta * hw) then
+        switch LiquidFlow
+                case "on":
                 "Francis Equation"
-                OutletL.F = (1.84/'s'*lw*((Level-(beta*hw))/(beta))^2/vL);
-        else
+                OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw)+1e-6*'m')/(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.P = OutletV.P + Level*g*rhoL + rhoV*alfa*(InletV.F*vV/Ah)^2;
+                InletV.F*vV = sqrt((InletV.P - OutletV.P - Level*g*rhoL + 1e-8 * 'atm')/(rhoV*alfa))*Ah;
+                when InletV.P < OutletV.P + Level*g*rhoL switchto "off";
                 
-        "Pressure Drop through the tray"
-        OutletV.F = (1 + tanh(1 * (OutletV.P - Pstartup)/'Pa'))/2 * 
-                Ah/vV * sqrt(2*(OutletV.P - InletL.P + 1e-8 * 'atm') / (alfa*rhoV) );
+                case "off":
+                InletV.F = 0 * 'mol/s';
+                when InletV.P > OutletV.P + Level*g*rhoL + 3e-2 * 'atm' switchto "on";
+                #when InletV.P > OutletV.P + Level*beta*g*rhoL + 1e-2 * 'atm' switchto "on";
+        end
 
         "Chemical Equilibrium"