source: branches/newlanguage/sample/reactors/fogler/chap4/molarflow_pfr.mso @ 188

Last change on this file since 188 was 188, checked in by Rodolfo Rodrigues, 15 years ago

Updated reactor samples for the new language

  • Property svn:keywords set to Id
File size: 2.8 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 188 2007-03-07 16:53:12Z rodolfo $
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        "Input total molar concentration"
91        Cto = 0.1*'mol/l';
92       
93        "Input molar flow"
94        Fo = [10, 0]*'mol/min';
95       
96        INITIAL
97        "Molar flow of A"
98        F(1) = 10*'mol/min';
99       
100        OPTIONS
101        TimeStep = 1;
102        TimeEnd = 100;
103end
Note: See TracBrowser for help on using the repository browser.