Changeset 511

Timestamp:

May 12, 2008, 2:25:33 PM (15 years ago)

Author:

Paula Bettio Staudt

Message:

More hydraulic correlations in tray model (part 1)

Location:

trunk/eml/stage_separators

Files:

• column.mso (modified) (2 diffs)
• icon/PackedDistillationKettleCond.png (added)
• icon/PackedDistillationKettleCond.svg (added)
• icon/PackedSectionColumn.png (added)
• icon/PackedSectionColumn.svg (added)
• icon/PackedStage.png (added)
• icon/PackedStage.svg (added)
• tray.mso (modified) (3 diffs)

trunk/eml/stage_separators/column.mso

@@ -1270 +1270 @@
         ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/SectionColumn"; 
+        Icon            = "icon/PackedSectionColumn"; 
         Brief           = "Model of a packed column section.";
         Info            =
@@ -1326 +1326 @@
         ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/DistillationKettleCond"; 
+        Icon            = "icon/PackedDistillationKettleCond"; 
         Brief           = "Model of a distillation column with dynamic condenser and dynamic reboiler.";
         Info            =

trunk/eml/stage_separators/tray.mso

@@ -130 +130 @@
         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);
+        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");
+
+        SET
+        Mw = PP.MolecularWeight();
         
         VARIABLES
@@ -159 +166 @@
         switch VapourFlow
                 case "on":
-                InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;
+                        switch VapourFlowModel
+                                case "Reepmeyer":
+                                InletV.F*vV = sqrt((InletV.P - OutletV.P)/(rhoV*alfa))*Ah;
+                        
+                                case "Feehery":
+                                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';
+                        end
                 when InletV.F < 1e-6 * 'kmol/h' switchto "off";
                 
@@ -330 +346 @@
         ATTRIBUTES
         Pallete         = false;
-        Icon            = "icon/Tray"; 
+        Icon            = "icon/PackedStage"; 
         Brief           = "Complete model of a packed column stage.";
         Info            =