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

Timestamp:

Jan 26, 2009, 8:08:28 PM (14 years ago)

Author:

Argimiro Resende Secchi

Message:

Removing temporary variable from tray model.

File:

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

trunk/eml/stage_separators/tray.mso

@@ -144 +144 @@
         beta as fraction (Brief="Aeration fraction");
         alfa as fraction (Brief="Dry pressure drop coefficient");
-        w as Real (Brief="Feehery's correlation coefficient", Unit='1/m^4', Default=1);
-        btray as Real (Brief="Elgue's correlation coefficient", Unit='kg/m/mol^2', Default=1);
-        fw as Real      (Brief = "Olsen's correlation coefficient", Default=1);
+        w as Real (Brief="Feehery correlation coefficient", Unit='1/m^4', Default=1);
+        btray as Real (Brief="Elgue correlation coefficient", Unit='kg/m/mol^2', Default=1);
+        fw as Real      (Brief = "Olsen correlation coefficient", Default=1);
         Np as Real      (Brief = "Number of liquid passes in the tray", Default=1);
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
@@ -161 +161 @@
         rhoL as dens_mass;
         rhoV as dens_mass;
-
-        btemp as Real (Brief="Temporary variable of Roffel's liquid flow equation");
         
         EQUATIONS
@@ -187 +185 @@
                         
                                 case "Roffel_Fl":
-                                OutletL.F = 2/3*sqrt(2*g)*rhoL/sum(Mw*OutletL.z)*lw*(2*btemp-1)*(ML*sum(Mw*OutletL.z)/(Ap*1.3)/rhoL/(2*btemp-1))^1.5;
+                                OutletL.F = 2/3*rhoL/sum(Mw*OutletL.z)*lw*(ML*sum(Mw*OutletL.z)/(Ap*1.3)/rhoL)^1.5*sqrt(2*g/
+                                                        (2*(1 - 0.3593/'Pa^0.0888545'*(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709)-1));  
                         end
                 when Level < (beta * hw) switchto "off";
@@ -197 +196 @@
         end
         
-        btemp = 1 - 0.3593/'Pa^0.0888545'*(OutletV.F*sum(Mw*OutletV.z)/(Ap*1.3)/sqrt(rhoV))^0.177709; #/'(kg/m)^0.0888545/s^0.177709';
-
         switch VapourFlow
                 case "on":