source: branches/new_gui/my_folders/fogler/chap4/molarflow_pfr.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: 2.7 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* Working with molar flow in a PFR
17*----------------------------------------------------------------------
18* Solved problem from Fogler (1999)
19* Problem number: 4-9
20* Page: 160 (Brazilian version, 2002)
21*----------------------------------------------------------------------
22*
23*   Description:
24*               Sample to calculate of the molar flows as function of the
25*       volume in a PFR with the reaction:
26*                       A <-> 2B
27*       
28*   Assumptions:
29*               * change time in volume
30*               * elementary reaction
31*       * steady-state
32*       * isotermic and isobaric system
33*       * gaseous phase
34*
35*       Specify:
36*               * the inlet stream
37*               * the kinetic parameters
38*
39*----------------------------------------------------------------------
40* Author: Christiano D. W. Guerra and Rodolfo Rodrigues
41* $Id: molarflow_pfr.mso 574 2008-07-25 14:18:50Z rafael $
42*--------------------------------------------------------------------*#
43
44using "types";
45
46
47FlowSheet MF_PFR
48        PARAMETERS
49        NComp   as Integer;
50        ka              as Real                 (Brief="Specific rate of reaction", Unit='1/min');
51        Kc              as conc_mol     (Brief="Equilibrium constant");
52       
53        VARIABLES
54    F(NComp)as flow_mol         (Brief="Molar flow", DisplayUnit='mol/min');
55        C(NComp)as conc_mol     (Brief="Molar concentration", DisplayUnit='mol/l');
56        r(NComp)as reaction_mol (Brief="Reaction rate", DisplayUnit='mol/min/l');
57        Fo(NComp) as flow_mol   (Brief="Input molar flow of A", DisplayUnit='mol/min');
58    Ft          as flow_mol     (Brief="Total molar flow", DisplayUnit='mol/min');
59        Cto     as conc_mol     (Brief="Initial concentration", DisplayUnit='mol/l');
60        V               as volume               (Brief="Reactor volume", DisplayUnit='l');
61       
62        EQUATIONS
63        "Change time in V"
64        V = time*'l/s';
65
66        "Molar balance"
67        diff(F(1))= r(1)*'l/s';
68       
69        "Reaction rate of A"
70        (-r(1)) = ka*(C(1) - C(2)^2/Kc);
71       
72        "Reaction rate of B"
73        r(2) = 2*(-r(1));
74
75        "Molar flow of B"
76        F(2) = 2*(Fo(1) - F(1));
77       
78        "Total molar flow"
79        Ft = sum(F);
80       
81        "Molar concentration"
82        C = Cto*F/Ft;
83       
84        SET
85        NComp = 2; # components A and B
86        ka = 2.7*'1/min';
87        Kc = 1.2*'mol/l';
88       
89    SPECIFY
90        Cto = 0.1*'mol/l';
91       
92        Fo = [10, 0]*'mol/min';
93       
94        INITIAL
95        "Molar flow of A"
96        F(1) = 10*'mol/min';
97       
98        OPTIONS
99        TimeStep = 1;
100        TimeEnd = 100;
101end
Note: See TracBrowser for help on using the repository browser.