1 | #*------------------------------------------------------------------- |
---|
2 | * EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC. |
---|
3 | * |
---|
4 | * This LIBRARY is free software; you can distribute it and/or modify |
---|
5 | * it under the therms of the ALSOC FREE LICENSE as available at |
---|
6 | * http://www.enq.ufrgs.br/alsoc. |
---|
7 | * |
---|
8 | * EMSO Copyright (C) 2004 - 2007 ALSOC, original code |
---|
9 | * from http://www.rps.eng.br Copyright (C) 2002-2004. |
---|
10 | * All rights reserved. |
---|
11 | * |
---|
12 | * EMSO is distributed under the therms of the ALSOC LICENSE as |
---|
13 | * available at http://www.enq.ufrgs.br/alsoc. |
---|
14 | * |
---|
15 | *--------------------------------------------------------------------- |
---|
16 | * Expressing Cj=hj(X) |
---|
17 | *---------------------------------------------------------------------- |
---|
18 | * Solved problem from Fogler (1999) |
---|
19 | * Problem number: 3-7 |
---|
20 | * Page: 87 (Brazilian edition, 2002) |
---|
21 | *---------------------------------------------------------------------- |
---|
22 | * |
---|
23 | * Description: |
---|
24 | * Expressing of the molar concentration as function of molar |
---|
25 | * conversion for a continuous reactor which occurs the oxidation |
---|
26 | * reaction: |
---|
27 | * 2SO2 + O2 -> 2SO3 |
---|
28 | * |
---|
29 | * Assumptions: |
---|
30 | * * steady-state |
---|
31 | * * isotermic and isobaric system |
---|
32 | * * gaseous phase |
---|
33 | * |
---|
34 | * Specify: |
---|
35 | * * the inlet stream (z,P,T) |
---|
36 | * * the parameters of reaction |
---|
37 | * * the outlet conversion |
---|
38 | * |
---|
39 | *---------------------------------------------------------------------- |
---|
40 | * Author: Christiano D. W. Guerra and Rodolfo Rodrigues |
---|
41 | * $Id: oxidation_of_so2.mso 188 2007-03-07 16:53:12Z rodolfo $ |
---|
42 | *--------------------------------------------------------------------*# |
---|
43 | |
---|
44 | using "types"; |
---|
45 | |
---|
46 | |
---|
47 | #*--------------------------------------------------------------------- |
---|
48 | * Model of a stream |
---|
49 | *--------------------------------------------------------------------*# |
---|
50 | |
---|
51 | Model stream |
---|
52 | PARAMETERS |
---|
53 | outer NComp as Integer (Brief="Number of chemical components", Lower=1); |
---|
54 | |
---|
55 | VARIABLES |
---|
56 | C(NComp)as conc_mol(Brief="Concentration", DisplayUnit='mol/l', Lower=0); |
---|
57 | z(NComp)as fraction(Brief="Molar fraction"); |
---|
58 | end |
---|
59 | |
---|
60 | |
---|
61 | #*--------------------------------------------------------------------- |
---|
62 | * Example 3-7: Cj=hj(X) |
---|
63 | *--------------------------------------------------------------------*# |
---|
64 | |
---|
65 | FlowSheet pfr |
---|
66 | PARAMETERS |
---|
67 | NComp as Integer; |
---|
68 | stoic(NComp) as Real(Brief="Stoichiometric coefficients"); |
---|
69 | k as Real (Brief="Specific rate of reaction", Unit='l/mol/s'); |
---|
70 | R as Real (Brief="Universal gas constant", Unit='atm*l/mol/K', Default=0.082); |
---|
71 | |
---|
72 | VARIABLES |
---|
73 | Inlet as stream; # Inlet stream |
---|
74 | Outlet as stream; # Outlet stream |
---|
75 | X as fraction (Brief="Molar conversion", Lower=0); |
---|
76 | r as reaction_mol (Brief="Rate of reaction of A", DisplayUnit='mol/l/s'); |
---|
77 | T as temperature (Brief="Temperature"); |
---|
78 | P as pressure (Brief="Pressure"); |
---|
79 | Theta(NComp)as Real (Brief="Parameter Theta"); |
---|
80 | epsilon as Real (Brief="Parameter epsilon"); |
---|
81 | |
---|
82 | EQUATIONS |
---|
83 | "Change time in X" |
---|
84 | X = time*'1/s'; |
---|
85 | |
---|
86 | "Outlet molar fraction" |
---|
87 | Outlet.C = Outlet.z*sum(Outlet.C); |
---|
88 | |
---|
89 | "Inlet concentration" |
---|
90 | Inlet.C = Inlet.z*P/(R*T); |
---|
91 | |
---|
92 | "Outlet concentration" |
---|
93 | Outlet.C = Inlet.C(1)*(Theta + stoic*X)/(1 + epsilon*X); |
---|
94 | |
---|
95 | "Parameter Theta" |
---|
96 | Theta = Inlet.z/Inlet.z(1); |
---|
97 | |
---|
98 | "Parameter epsilon" |
---|
99 | epsilon = Inlet.z(1)*sum(stoic); |
---|
100 | |
---|
101 | "Rate of reaction" |
---|
102 | (-r) = k*Outlet.C(1)*Outlet.C(2); |
---|
103 | |
---|
104 | SET |
---|
105 | NComp = 4; # A, B, C and I |
---|
106 | stoic = [-1.0, -0.5, 1.0, 0.0]; |
---|
107 | k = 200*'l/mol/s'; |
---|
108 | |
---|
109 | SPECIFY |
---|
110 | "Inlet molar fraction" |
---|
111 | Inlet.z = [0.28, 0.1512, 0.0, 0.5688]; |
---|
112 | "Inlet pressure" |
---|
113 | P = 1485*'kPa'; |
---|
114 | "Inlet temperature" |
---|
115 | T = (227 + 273.15)*'K'; |
---|
116 | |
---|
117 | OPTIONS |
---|
118 | TimeStep = 0.005; |
---|
119 | TimeEnd = 0.995; |
---|
120 | end |
---|