Changeset 513 for trunk

Timestamp:

May 13, 2008, 4:07:40 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

More hydraulic correlations implemented in tray model

File:

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

trunk/eml/stage_separators/tray.mso

@@ -121 +121 @@
 * the liquid height (Level) OR the liquid flow OutletL.F
 * (NoComps - 1) OutletL compositions
+
+== Options ==
+You can choose the equation for the liquid outlet flow and the vapour
+inlet flow calculation through the VapourFlowModel and LiquidFlowModel
+switchers.
+
+== References ==
+* ELGUE, S.; PRAT, L.; CABASSUD, M.; LANN, J. L.; CéZERAC, J. Dynamic models for start-up operations of batch distillation columns with experimental validation. Computers and Chemical Engineering, v. 28, p. 2735-2747, 2004.
+* FEEHERY, W. F. Dynamic Optimization with Path Constraints. Tese (Doutorado) - Massachusetts Institute of Technology, June 1998.
+* KLINGBERG, A. Modeling and Optimization of Batch Distillation. Dissertação (Mestrado) - Department of Automatic Control, Lund Institute of Technology, Lund, Sweden, fev. 2000.
+* OLSEN, I.; ENDRESTOL, G. O.; SIRA, T. A rigorous and efficient distillation column model for engineering and training simulators. Computers and Chemical Engineering,v. 21, n. Suppl, p. S193-S198, 1997.
+* REEPMEYER, F.; REPKE, J.-U.; WOZNY, G. Analysis of the start-up process for reactive distillation. Chemical Engineering Technology, v. 26, p. 81-86, 2003.
+* ROFFEL, B.; BETLEM, B.; RUIJTER, J. de. First principles dynamic modeling and multivariable control of a cryogenic distillation column process. Computers and Chemical Engineering, v. 24, p. 111-123, 2000.
+* WANG, L.; LI, P.; WOZNY, G.; WANG, S. A start-up model for simulation of batch distillation starting from a cold state. Computers and Chemical Engineering, v. 27, p.1485-1497, 2003.
 ";      
 
@@ -131 +145 @@
         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);
+        Np as Real      (Brief = "Number of liquid passes in the tray", Default=1);
         Mw(NComp)       as molweight    (Brief = "Component Mol Weight");
         
         VapourFlow as Switcher(Valid = ["on", "off"], Default = "on");
         LiquidFlow as Switcher(Valid = ["on", "off"], Default = "on");
-        VapourFlowModel as Switcher(Valid = ["Reepmeyer", "Feehery", "Roffel"], Default = "Reepmeyer");
-#       LiquidFlow as Switcher(Valid = ["on", "off"], Default = "on");
+        VapourFlowModel as Switcher(Valid = ["Reepmeyer", "Feehery_Fv", "Roffel_Fv", "Klingberg", "Wang_Fv", "Elgue"], Default = "Reepmeyer");
+        LiquidFlowModel as Switcher(Valid = ["default", "Wang_Fl", "Olsen", "Feehery_Fl", "Roffel_Fl"], Default = "default");
 
         SET
@@ -145 +162 @@
         rhoV as dens_mass;
 
+        btemp as Real (Brief="Temporary variable of Roffel's liquid flow equation");
+        
         EQUATIONS
         "Liquid Density"
@@ -153 +172 @@
         switch LiquidFlow
                 case "on":
-                "Francis Equation"
-#               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;
+                        switch LiquidFlowModel
+                                case "default":
+                                "Francis Equation"
+                                OutletL.F*vL = 1.84*'1/s'*lw*((Level-(beta*hw))/(beta))^2;
+                        
+                                case "Wang_Fl":
+                                OutletL.F*vL = 1.84*'m^0.5/s'*lw*((Level-(beta*hw))/(beta))^1.5;
+                        
+                                case "Olsen":
+                                OutletL.F / 'mol/s'= lw*Np*rhoL/sum(Mw*OutletV.z)/(0.665*fw)^1.5 * ((ML*sum(Mw*OutletL.z)/rhoL/Ap)-hw)^1.5 * 'm^0.5/mol';
+                        
+                                case "Feehery_Fl":
+                                OutletL.F = lw*rhoL/sum(Mw*OutletL.z) * ((Level-hw)/750/'mm')^1.5 * 'm^2/s';
+                        
+                                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;
+                        end
                 when Level < (beta * hw) switchto "off";
                 
@@ -163 +196 @@
                 when Level > (beta * hw) + 1e-6*'m' switchto "on";
         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
@@ -170 +205 @@
                                 InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;
                         
-                                case "Feehery":
+                                case "Feehery_Fv":
                                 InletV.F = rhoV/Ap/w/sum(Mw*OutletV.z) * sqrt(((InletV.P - OutletV.P)-(rhoV*g*ML*vL/Ap))/rhoV);
                         
-                                case "Roffel":
-                                InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92' = (InletV.P - OutletV.P) - (beta*sum(M*Mw)/(Ap*1)*g) * (rhoV*Ah/sum(Mw*OutletV.z))^1.08 * 'm^1.08/mol^1.08';
+                                case "Roffel_Fv":
+                                InletV.F^1.08 * 0.0013 * 'kg/m/mol^1.08/s^0.92*1e5' = (InletV.P - OutletV.P)*1e5 - (beta*sum(M*Mw)/(Ap*1.3)*g*1e5) * (rhoV*Ah/sum(Mw*OutletV.z))^1.08 * 'm^1.08/mol^1.08';
+                        
+                                case "Klingberg":
+                                InletV.F * vV = Ap * sqrt(((InletV.P - OutletV.P)-rhoL*g*Level)/rhoV);
+                        
+                                case "Wang_Fv":
+                                InletV.F * vV = Ap * sqrt(((InletV.P - OutletV.P)-rhoL*g*Level)/rhoV*alfa);
+                                
+                                case "Elgue":
+                                InletV.F  = sqrt((InletV.P - OutletV.P)/btray);
                         end
                 when InletV.F < 1e-6 * 'kmol/h' switchto "off";