1 | #*--------------------------------------------------------------------- |
---|
2 | * This file is property of the author and cannot be used, copyed |
---|
3 | * or modified without permission. |
---|
4 | * |
---|
5 | * Copyright (C) 2002-2004 the author |
---|
6 | *---------------------------------------------------------------------- |
---|
7 | * Author: Paula Bettio Staudt |
---|
8 | * Rafael de Pelegrini Soares |
---|
9 | * $Id: sample_batch.mso 8 2006-08-26 18:05:18Z rafael $ |
---|
10 | *---------------------------------------------------------------------- |
---|
11 | * |
---|
12 | * Sample file for a isotermic batch reactor. |
---|
13 | * |
---|
14 | *-------------------------------------------------------------------*# |
---|
15 | |
---|
16 | using "reactors/batch"; |
---|
17 | |
---|
18 | FlowSheet SampleBatch |
---|
19 | PARAMETERS |
---|
20 | NComp as Integer(Brief="Number of components"); |
---|
21 | |
---|
22 | DEVICES |
---|
23 | feed as stream; |
---|
24 | reac as batch_isotermic; |
---|
25 | |
---|
26 | CONNECTIONS |
---|
27 | feed to reac.Inlet; |
---|
28 | |
---|
29 | PARAMETERS |
---|
30 | k1 as Real(Unit="(l/mol)^3.5/min", Default=18.7); |
---|
31 | k2 as Real(Unit="l/mol/min", Default=0.58); |
---|
32 | k3 as Real(Unit="(l/mol)^2/min", Default=0.09); |
---|
33 | k4 as Real(Unit="(l/mol)^1.5/min", Default=0.42); |
---|
34 | Ks as inv_conc_mol(Default=115.83); |
---|
35 | K as Real(Default=34.4); |
---|
36 | kla as flow_vol(Default=3.3); |
---|
37 | H as inv_conc_mol(Default=737); |
---|
38 | pCO2; |
---|
39 | |
---|
40 | VARIABLES |
---|
41 | C6 as conc_mol(Brief="Special component concentration"); |
---|
42 | |
---|
43 | SET |
---|
44 | NComp = 5; |
---|
45 | reac.NReac = 5; |
---|
46 | pCO2 = 0.9; |
---|
47 | |
---|
48 | reac.stoic(1,:) = [-2, 1, -1, -1, 0 ]; |
---|
49 | reac.stoic(2,:) = [-0.5, 0, 0, -1, -0.5]; |
---|
50 | reac.stoic(3,:) = [1, -1, 1, 0, 0 ]; |
---|
51 | reac.stoic(4,:) = [0, -1, 1, -2, 0 ]; |
---|
52 | reac.stoic(5,:) = [0, 1, -1, 0, 1 ]; |
---|
53 | |
---|
54 | EQUATIONS |
---|
55 | "Inlet Flow" |
---|
56 | feed.F = kla*(pCO2/H-reac.C(2)); |
---|
57 | |
---|
58 | "Reaction rates" |
---|
59 | reac.r(1)=k1 * reac.C(1)^4 * sqrt(reac.C(2)); |
---|
60 | reac.r(2)=k2 * reac.C(3) * reac.C(4); |
---|
61 | reac.r(3)=k2/K * reac.C(1) * reac.C(5); |
---|
62 | reac.r(4)=k3 * reac.C(1) * reac.C(4)^2; |
---|
63 | reac.r(5)=k4 * C6^2 * sqrt(reac.C(2)); |
---|
64 | |
---|
65 | "C6 Equilibrium" |
---|
66 | C6 = Ks * reac.C(1) * reac.C(4); |
---|
67 | |
---|
68 | SPECIFY |
---|
69 | "Feed conditions" |
---|
70 | feed.T = 300*"K"; |
---|
71 | feed.P = 1*"atm"; |
---|
72 | feed.h = 0*"J/mol"; |
---|
73 | feed.v = 0; |
---|
74 | feed.z = [0,1,0,0,0]; |
---|
75 | |
---|
76 | "Reactor volume" |
---|
77 | reac.Vr = 1 * "l"; |
---|
78 | |
---|
79 | INITIAL |
---|
80 | "Initial concentration" |
---|
81 | reac.C=[0.444, 0.00122, 0, 0.007, 0] * "mol/l"; |
---|
82 | |
---|
83 | OPTIONS |
---|
84 | outputLevel = "high"; |
---|
85 | time = [0:100:180*60]; |
---|
86 | end |
---|