Changeset 636 for branches

Timestamp:

Sep 26, 2008, 6:59:05 PM (15 years ago)

Author:

gerson bicca

Message:

added initialization mode for column model (temperature and composition)

Location:

branches/gui/eml/stage_separators

Files:

• columnTeste.mso (modified) (4 diffs)
• reboilerTeste.mso (added)

branches/gui/eml/stage_separators/columnTeste.mso

@@ -32 +32 @@
 
 using "trayTeste";
-using "reboiler";
+using "reboilerTeste";
 using "condenser";
 using "mixers_splitters/splitter";
@@ -68 +68 @@
         outer PP                                                                                as Plugin                       (Brief="External Physical Properties", Type="PP");
         outer NComp                                                             as Integer                      (Brief="Number of components");
-        NumberOfTrays                                                           as Integer                      (Brief="Number of trays", Default=2);
+        NumberOfTrays                                                           as Integer                      (Brief="Number of trays", Default=8);
         SpecialTrayIndex(NumberOfTrays)         as Integer                      (Brief="Number of trays", Default=0,Hidden=true);
         FeedTrayLocation                                                        as Integer                      (Brief="Feed tray Location", Default=2);
@@ -88 +88 @@
         OlsenCoeff                      as Real                         (Brief="Olsens correlation coefficient", Default=1);
         TrayLiquidPasses        as Real                         (Brief="Number of liquid passes in the tray", Default=1);
+        
+        TopTemperature          as temperature;
+        TbottomTemperature      as temperature;
+        TopComposition(NComp)                   as fraction     (Brief = "Component Molar Fraction at Top");
+        BottomComposition(NComp)                as fraction     (Brief = "Component Molar Fraction at Bottom");
         
         V                                       as volume               (Brief="Total Volume of the tray",Hidden=true);
@@ -256 +261 @@
 end
 
-end
-
+INITIAL
+
+# The initial temperature of the trays
+for i in 1:NumberOfTrays do
+        
+"The initial temperature of the trays"
+        trays(i).OutletL.T = TopTemperature+(TbottomTemperature-TopTemperature)*((i-1)/(NumberOfTrays-1));
+
+end
+
+for i in 1:NComp-1 do
+        
+for j in 1:NumberOfTrays do
+
+"The initial composition of the trays"
+        trays(j).OutletL.z(i) = TopComposition(i) +(BottomComposition(i)-TopComposition(i) )*((j-1)/(NumberOfTrays-1));
+
+end
+
+end
+
+end
 
 #*----------------------------------------------------------------------