Changeset 506 for branches/rate

Timestamp:

Apr 28, 2008, 3:46:20 PM (15 years ago)

Author:

Gabriela Sporleder Straatmann

Message:

Adding a example for trayRate.

File:

• branches/rate/sample/stage_separators/sample_tray.mso (modified) (1 diff)

branches/rate/sample/stage_separators/sample_tray.mso

@@ -215 +215 @@
         TimeEnd = 100;
 end
+
+FlowSheet trayRate_Test
+        PARAMETERS
+        PP      as Plugin(Brief="Physical Properties",
+                Type="PP",
+                Components = [ "n-pentane", "benzene"],
+                LiquidModel = "PR",
+                VapourModel = "PR"
+        );
+        NComp   as Integer;
+
+        SET
+        NComp = PP.NumberOfComponents;
+        
+        DEVICES
+        t1 as trayRate;
+        feed as source;
+        feedV as source;
+        inL as liquid_stream;
+        inV as vapour_stream;
+        
+        CONNECTIONS
+        feed.Outlet to t1.Inlet;
+        feedV.Outlet to t1.InletFV;
+        inL to t1.InletL;
+        inV to t1.InletV;
+        
+        SPECIFY
+        feed.Outlet.F = 113.4 * 'kmol/h';
+        feed.Outlet.T = 291 * 'K';
+        feed.Outlet.P = 1.66 * 'atm';
+        feed.Outlet.z = [0.5, 0.5];
+        
+        feedV.Outlet.F = 0 * 'kmol/h';
+        feedV.Outlet.T = 291 * 'K';
+        feedV.Outlet.P = 1.66 * 'atm';
+        feedV.Outlet.z = [0.5, 0.5];
+        
+        inL.P = 165 * 'kPa';
+        inL.T = 310 * 'K';
+        inL.F = 61.99 * 'kmol/h';
+        inL.z = [0.1641, 0.8359];
+
+        inV.P = 150 * 'kPa';
+        inV.T = 321 * 'K';
+        inV.z = [0.0584, 0.9416];
+        inV.F = 61.99 * 'kmol/h';
+        
+        #t1.Emv = 1;
+        t1.OutletV.F = 147.1 * 'kmol/h';
+        t1.OutletL.F = 147.1 * 'kmol/h';
+        t1.interf.a=0.01*'m^2';
+        t1.interf.htL=0.4*'kW/m^2/K';
+        t1.interf.htV=0.4*'kW/m^2/K';
+        t1.interf.kL=0.001*'m/s';
+        t1.interf.kV=0.001*'m/s';
+        t1.ML=3*'mol';
+        t1.MV=3*'mol';
+
+        SET
+        t1.V = 4 * 'ft^3';
+        #*t1.Ah = 0.394 * 'ft^2';
+        t1.lw = 20.94 * 'in';
+        t1.hw = 0.125 * 'ft';
+        t1.beta = 0.6;
+        t1.alfa = 4;*#
+        
+        t1.Ap = 3.94 * 'ft^2';
+        t1.Q = 0 * 'kW';
+        
+        
+        INITIAL
+        t1.OutletL.T = 290 *'K';
+        t1.OutletV.T = 290 *'K';
+        #t1.Level = 0.9*'m' ;
+        t1.OutletL.z(1) = 0.5;
+        t1.OutletV.z(1) = 0.5;
+        #t1.MV = 0.01*'mol' ;
+        
+        OPTIONS
+        TimeEnd = 100;
+end