source: branches/new_gui/my_folders/fogler/chap3/oxidation_of_so2.mso @ 896

Last change on this file since 896 was 896, checked in by gerson bicca, 13 years ago

added special folder for unconventional applications

File size: 3.3 KB
Line 
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 574 2008-07-25 14:18:50Z rafael $
42*--------------------------------------------------------------------*#
43
44using "types";
45
46
47#*---------------------------------------------------------------------
48* Model of a stream
49*--------------------------------------------------------------------*#
50
51Model stream
52        PARAMETERS
53outer 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");
58end
59
60
61#*---------------------------------------------------------------------
62* Example 3-7: Cj=hj(X)
63*--------------------------------------------------------------------*#
64
65FlowSheet 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");
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.z = [0.28, 0.1512, 0.0, 0.5688];
111        P = 1485*'kPa';
112        T = (227 + 273.15)*'K';
113       
114        OPTIONS
115        TimeStep = 0.005;
116        TimeEnd = 0.995;
117end
Note: See TracBrowser for help on using the repository browser.