Ticket #124: ColunaDissertação.mso

File ColunaDissertação.mso, 2.7 KB (added by Paula Bettio Staudt, 13 years ago)
Line 
1
2using "stage_separators/column";
3
4# column section with 8 trays
5FlowSheet Distillation_kettle_cond_Test
6        PARAMETERS
7        PP      as Plugin(Brief="External Physical Properties",
8                Type="PP",
9                Components = [ "o-xylene", "m-xylene", "p-xylene","n-propylbenzene"],
10                LiquidModel = "PR",
11                VapourModel = "PR"
12        );
13        NComp   as Integer;
14       
15        VARIABLES
16        Qc as energy_source (Brief="Heat rate removed from condenser");
17        Qr as energy_source (Brief="Heat rate supplied to reboiler");
18       
19        SET
20        NComp = PP.NumberOfComponents;
21
22        DEVICES
23        col as Distillation_kettle_cond;
24        feed as source;
25        zero as stream;
26       
27        CONNECTIONS
28        feed.Outlet to col.trays(65).Inlet;
29        zero to col.reb.Inlet;
30        zero to col.trays([1:64]).Inlet;
31        zero to col.trays([66:col.NTrays]).Inlet;
32        Qc.OutletQ to col.cond.InletQ;
33        Qr.OutletQ to col.reb.InletQ;
34       
35        SPECIFY
36        feed.Fw = 71500 * 'kg/h';               #fornecida pelo artigo de Fabio, 71500 kg/h - deve ser colocada em mol/s - Calcular novo valor
37        feed.T = (172+273.15) * 'K';    #fornecida pelo artigo de Fabio
38        feed.P = 2.4 * 'kgf/cm^2';              #fornecida pelo artigo de Fabio
39        feed.Composition =[0.242, 0.170, 0.585, 0.003]; #estimada por mim com base em informações do artigo
40
41       
42        zero.F = 0 * 'kmol/h';                  #corrente necessária para conexão
43        zero.T = 300 * 'K';
44        zero.P = 1 * 'atm';
45        zero.z = 1/NComp;
46        zero.v = 0;
47        zero.h = 0 * 'J/mol';
48       
49        #INSERIR MAIS DADOS DA DISSERTAÇÃO
50       
51        col.sptop.Outlet2.F = 280000 * 'kmol/h';
52        col.reb.OutletL.F = 6000 * 'kmol/h';
53        col.sptop.frac = 0.2;
54        col.cond.OutletV.F = 0 * 'kmol/h';
55        Qr.OutletQ.Q = 3.7743e6 * 'kJ/h';
56        Qc.OutletQ.Q = -3.71e6 * 'kJ/h';
57        col.pump1.dP = 16 * 'kPa';
58        col.trays.Emv = 1;
59        #col.cond.OutletL.z(1)=0.98;            #fornecida pelo artigo
60       
61        col.alfaTopo = 2;
62       
63        SET
64        col.NTrays = 130;
65        col.cond.V = 2 * 'm^3';
66        col.cond.Across = 1 * 'm^2';
67        col.trays.V = 4 * 'ft^3';
68        col.trays.Ah = 0.394 * 'ft^2';
69        col.trays.lw = 20.94 * 'in';
70        col.trays.hw = 0.125 * 'ft';
71        col.trays.Q = 0 * 'kW';
72        col.trays.beta = 0.8;
73        col.trays.alfa = 4;
74        col.trays.Ap = 3.94 * 'ft^2';
75        col.reb.V = 2 * 'm^3';
76        col.reb.Across = 1 * 'm^2';
77       
78        INITIAL
79        # condenser
80        col.cond.OutletL.T = 260 *'K';
81        col.cond.Level = 1 * 'm';
82        col.cond.OutletL.z([1:3]) = [0.98, 0.01, 0.01]; #0.98 do artigo
83
84        # reboiler
85        col.reb.OutletL.T = 350 *'K';
86        col.reb.Level = 1 * 'm';
87        col.reb.OutletL.z([1:3]) = [0.1, 0.7, 0.15];
88
89        # column trays
90        col.trays.OutletL.T = [290:(330-290)/(col.NTrays-1):330] * 'K';
91        col.trays.Level = 1.2 * col.trays.hw;
92        col.trays.OutletL.z([1:3]) = [0.5, 0.05, 0.01];
93
94        OPTIONS
95        TimeStep = 0.1;
96        TimeEnd = 2;
97        NLASolver(File="nlasolver");
98        #GuessFile="Distillation_kettle_cond_Test.rlt";
99        #Dynamic = false;       
100end