Changeset 270 for trunk/sample/stage_separators/sample_column.mso

Timestamp:

Jun 16, 2007, 4:18:47 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

Updated ATTRIBUTES section and some stage_separators icons

File:

• trunk/sample/stage_separators/sample_column.mso (modified) (2 diffs)

trunk/sample/stage_separators/sample_column.mso

@@ -281 +281 @@
 end
 
+FlowSheet Distillation_thermosyphon_subcooling_Test
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "isobutane", "n-pentane", "propylene",
+                "benzene", "isobutene" ],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+        NComp   as Integer;
+        
+        VARIABLES
+        Qc as heat_rate (Brief="Heat rate removed from condenser");
+        Qr as heat_rate (Brief="Heat rate supplied to reboiler");
+        Qttop as heat_rate (Brief="Heat rate removed from condenser");
+        Qtbottom as heat_rate (Brief="Heat rate supplied to reboiler");
+        
+        SET
+        NComp = PP.NumberOfComponents;
+
+        DEVICES
+        col as Distillation_thermosyphon_subcooling;
+        feed as source;
+        zero as stream;
+        
+        CONNECTIONS
+        feed.Outlet to col.trays(5).Inlet;
+        zero to col.trays([1:4]).Inlet;
+        zero to col.trays([6:col.NTrays]).Inlet;
+        Qttop to col.ttop.Q;
+        Qtbottom to col.tbottom.Q;
+        Qc to col.cond.Q;
+        Qr to col.reb.Q;
+        
+        SPECIFY
+        feed.Outlet.F = 113.4 * 'kmol/h';
+        feed.Outlet.T = 291 * 'K';
+        feed.Outlet.P = 168.3 * 'kPa';
+        feed.Outlet.z = 1/NComp;
+        
+        zero.F = 0 * 'kmol/h';
+        zero.T = 300 * 'K';
+        zero.P = 1 * 'atm';
+        zero.z = 1/NComp;
+        zero.v = 0;
+        zero.h = 0 * 'J/mol';
+        
+        col.sptop.Outlet2.F = 85 * 'kmol/h';
+        col.sptop.frac = 0.444445;
+        col.spbottom.Outlet1.F = 100 * 'kmol/h';
+        
+        Qr = 3.7743e6 * 'kJ/h';
+        Qc = -3.71e6 * 'kJ/h';
+        Qttop = 0 * 'kJ/h';
+        Qtbottom = 0 * 'kJ/h';
+        
+        col.pump1.dP = 16 * 'kPa';
+        col.trays.Emv = 1;
+        col.cond.DP = 0.5 * 'atm';
+        
+        SET
+        col.NTrays = 8;
+        col.trays.V = 4 * 'ft^3';
+        col.trays.Ah = 0.394 * 'ft^2';
+        col.trays.lw = 20.94 * 'in';
+        col.trays.hw = 0.125 * 'ft';
+        col.trays.Q = 0 * 'kW';
+        col.trays.beta = 0.6;
+        col.trays.alfa = 4;
+        col.trays.Ap = 3.94 * 'ft^2';
+        col.reb.DP = 0.5 * 'atm';
+        col.tbottom.Across = 2.20 * 'm^2';
+        
+        INITIAL
+        # ttop
+        col.ttop.Outlet.T = 260 *'K';
+        col.ttop.Level = 1 * 'm';
+        col.ttop.Outlet.z([1:4]) = [0.65, 0.05, 0.01, 0.01];
+
+        # tbottom
+        col.tbottom.Outlet.T = 330 *'K';
+        col.tbottom.Level = 1 * 'm';
+        col.tbottom.Outlet.z([1:4]) = [0.1, 0.7, 0.01, 0.01];
+
+        # column trays
+        col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
+        col.trays.Level = 1.2 * col.trays.hw;
+        col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01];
+
+        OPTIONS
+        TimeStep = 0.1;
+        TimeEnd = 50;
+        #time = [0:0.01:1, 2:50];
+        #GuessFile="Distillation_kettle_cond_Test.rlt";
+        #Dynamic = false;       
+end
+
 FlowSheet Column_ctrl
         PARAMETERS
@@ -493 +590 @@
 end
 
+