Changeset 842 for branches/gui/eml

Timestamp:

Sep 1, 2009, 3:00:55 PM (13 years ago)

Author:

gerson bicca

Message:

updated packed column

File:

• branches/gui/eml/stage_separators/column.mso (modified) (3 diffs)

branches/gui/eml/stage_separators/column.mso

@@ -2426 +2426 @@
         HeatSupply                      as heat_rate    (Brief="Rate of heat supply"); 
         ColumnDiameter          as length               (Brief="Column diameter");
-        VoidFraction            as Real                 (Brief="Void fraction of packing, m^3/m^3");
+        VoidFraction            as Real                 (Brief="Void fraction of packing, (m^3 void space/m^3 packed bed)");
         ResistanceCoeff         as positive     (Brief="Resistance coefficient on the liquid load", Default=1);
         AreaPerPackingVol       as Real                 (Brief="surface area per packing volume", Unit='m^2/m^3');
@@ -2558 +2558 @@
 end
 
-Model Packed_Section_Column as Packed_Section_ColumnBasic
+Model Packed_Section_Column     as Packed_Section_ColumnBasic
 
         ATTRIBUTES
@@ -2632 +2632 @@
 "Inlet Tray -  Enthalpy"
     0*'J/mol' = STAGES(i).Inlet.h;
+
+end
+
+end
+
+end
+
+Model Packed_Section_Column2    as Packed_Section_ColumnBasic
+
+        ATTRIBUTES
+        Pallete         = true;
+        Icon            = "icon/PackedSectionColumn"; 
+        Brief           = "Model of a packed column section.";
+        Info            =
+"== Model of a packed column section containing ==
+* NStages theoretical stages.
+        
+== Specify ==
+* the feed stream of each tray (Inlet);
+* the InletLiquid stream of the top tray;
+* the InletVapour stream of the bottom tray;
+* the total pressure drop (dP) of the section.
+        
+== Initial Conditions ==
+* the stages temperature (OutletLiquid.T);
+* the stages liquid holdup;
+* (NoComps - 1) OutletLiquid (OR OutletVapour) compositions for each tray.
+";
+
+SET
+        NumberOfFeeds = 2;
+
+VARIABLES
+        
+        in      FeedStageOne    as stream               (Brief="Feed stream", PosX=0, PosY=0.45);
+        in      FeedStageTwo    as stream               (Brief="Feed stream", PosX=0, PosY=0.65);
+
+EQUATIONS
+#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+# Equating Feed Stage Variables to Stages Variables
+#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
+
+for i in 1:NumberOfStages do
+
+if   i equal FeedStageLocation(1) then    
+
+"FeedTrayTop Inlet Flow"        
+        FeedStageOne.F= STAGES(i).Inlet.F;
+
+"FeedTrayTop Inlet Temperature"
+        FeedStageOne.T = STAGES(i).Inlet.T;
+
+"FeedTrayTop Inlet Pressure"
+        FeedStageOne.P = STAGES(i).Inlet.P;
+
+"FeedTrayTop Inlet Composition"
+        FeedStageOne.z = STAGES(i).Inlet.z;
+
+"FeedTrayTop Inlet Vapour Fraction"
+        FeedStageOne.v = STAGES(i).Inlet.v;
+
+"FeedTrayTop Inlet Enthalpy"
+        FeedStageOne.h = STAGES(i).Inlet.h;
+
+else if   i equal FeedStageLocation(2) then    
+
+"FeedTrayBottom Inlet Flow"     
+        FeedStageTwo.F= STAGES(i).Inlet.F;
+
+"FeedTrayBottom Inlet Temperature"
+        FeedStageTwo.T = STAGES(i).Inlet.T;
+
+"FeedTrayBottom Inlet Pressure"
+        FeedStageTwo.P = STAGES(i).Inlet.P;
+
+"FeedTrayBottom Inlet Composition"
+        FeedStageTwo.z = STAGES(i).Inlet.z;
+
+"FeedTrayBottom Inlet Vapour Fraction"
+        FeedStageTwo.v = STAGES(i).Inlet.v;
+
+"FeedTrayBottom Inlet Enthalpy"
+        FeedStageTwo.h = STAGES(i).Inlet.h;
+
+else
+
+"Inlet Tray - Flow Sealed"
+    0*'mol/h'= STAGES(i).Inlet.F;
+        
+"Inlet Tray -  Temperature"
+    300*'K' = STAGES(i).Inlet.T;
+        
+"Inlet Tray -  Pressure"
+    1*'atm' = STAGES(i).Inlet.P;
+        
+"Inlet Tray -  Composition"
+    0.1 = STAGES(i).Inlet.z;
+        
+"Inlet Tray -  Vapour Fraction"
+    0 = STAGES(i).Inlet.v;
+        
+"Inlet Tray -  Enthalpy"
+    0*'J/mol' = STAGES(i).Inlet.h;
+
+end
 
 end