Changeset 346 for trunk/sample

Timestamp:

Aug 19, 2007, 10:32:18 PM (16 years ago)

Author:

Argimiro Resende Secchi

Message:

Adding entropy in streamPH, sources and sinks.

Location:

trunk/sample/stage_separators

Files:

• Distillation_kettle_cond_Test.rlt (added)
• sample_column.mso (modified) (5 diffs)
• sample_columnReact.mso (modified) (3 diffs)

trunk/sample/stage_separators/sample_column.mso

@@ -276 +276 @@
         TimeStep = 0.1;
         TimeEnd = 50;
-        #time = [0:0.01:1, 2:50];
         #GuessFile="Distillation_kettle_cond_Test.rlt";
         #Dynamic = false;       
 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 energy_source (Brief="Heat rate removed from condenser");
-        Qr as energy_source (Brief="Heat rate supplied to reboiler");
-        Qttop as energy_source (Brief="Heat rate removed from condenser");
-        Qtbottom as energy_source (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.OutletQ to col.ttop.InletQ;
-        Qtbottom.OutletQ to col.tbottom.InletQ;
-        Qc.OutletQ to col.cond.InletQ;
-        Qr.OutletQ to col.reb.InletQ;
-        
-        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.OutletQ.Q = 3.7743e6 * 'kJ/h';
-        Qc.OutletQ.Q = -3.71e6 * 'kJ/h';
-        Qttop.OutletQ.Q = 0 * 'kJ/h';
-        Qtbottom.OutletQ.Q = 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
@@ -467 +370 @@
         PC.Parameters.gain = 0.5;
         PC.Parameters.derivTime = 1*'s';
-        PC.Ports.setPoint = (2.0*'bar'-Pmin)/(Pmax-Pmin);
+        PC.Ports.setPoint = (4.0*'bar'-Pmin)/(Pmax-Pmin);
         PC.Ports.input = Pad;
         Pad = (col.cond.OutletV.P-Pmin)/(Pmax-Pmin);
@@ -519 +422 @@
         feed.Outlet.F = 113.4 * 'kmol/h';
         feed.Outlet.T = 291 * 'K';
-        feed.Outlet.P = 168.3 * 'kPa';
+        feed.Outlet.P = 5 * 'bar';
         feed.Outlet.z = 1/NComp;
         
@@ -529 +432 @@
         zero.h = 0 * 'J/mol';
         
-        Qr.OutletQ.Q = 3e6 * 'kJ/h';
+        Qr.OutletQ.Q = 4e6 * 'kJ/h';
         col.pump1.dP = 16 * 'kPa';
         col.trays.Emv = 1;
@@ -565 +468 @@
         Hmaxb = 2 * 'm';
         Pmin = 0.5 * 'bar';
-        Pmax = 4 * 'bar';
+        Pmax = 6 * 'bar';
         Qcmax = -5e5 * 'kJ/h';
         Qcmin = -5e6 * 'kJ/h';

trunk/sample/stage_separators/sample_columnReact.mso

@@ -105 +105 @@
         feed.Outlet.F = 1.076 * 'mol/s';
         feed.Outlet.T = 300 * 'K';
-        feed.Outlet.P = 1.0 * 'atm';
+        feed.Outlet.P = 3.0 * 'atm';
         feed.Outlet.z = [0.4962, 0.4808, 0, 0.0229];
 
@@ -115 +115 @@
         zero.h = 0 * 'J/mol';
 
-        col.p.dP = 0.95 * 'atm' - col.sp.Outlet2.P;
+        col.p.dP = 2 * 'atm' - col.sp.Outlet2.P;
         col.trays.Emv = 1;
         
@@ -241 +241 @@
         # condenser
         col.cond.OutletL.T = 300 *'K';
-        col.cond.Level = 0.1 * 'm';
+        col.cond.Level = 0.4 * 'm';
         col.cond.OutletL.z([1:3]) = [0.4962, 0.4808, 0];
 
         # reboiler
         col.reb.OutletL.T = 300 * 'K';
-        col.reb.Level = 0.1 * 'm';
+        col.reb.Level = 0.4 * 'm';
         col.reb.OutletL.z([1:3]) = [0.4962, 0.4808, 0];