source: branches/newlanguage/sample/reactors/fogler/chap6/hidrodesalkeletion.mso @ 171

Last change on this file since 171 was 171, checked in by gerson bicca, 16 years ago

some modifications in the models for the new language (reactors model)

  • Property svn:keywords set to Id
File size: 3.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* Hidrodesalkylation of mesitylene
17*----------------------------------------------------------------------
18* Solved problem from Fogler (1999)
19* Problem number: 6-6
20* Page: 273 (Brazilian version, 2002)
21*----------------------------------------------------------------------
22*
23*   Description:
24*               Sample to calculate of the molar concentration as functions of
25*       the resident time in a PFR with a fixed bed. There is the reaction
26*       of hidrodesalkylation of mesitylene:
27*                       mesitylene + H2 -> m-xylene + CH4
28*                       m-xylene   + H2 -> toluene  + CH4
29*       
30*   Assumptions:
31*               * change time in residence time
32*       * steady-state
33*       * isotermic and isobaric system
34*       * gaseous phase
35*
36*       Specify:
37*               * the inlet stream
38*               * the kinetic parameters
39*               * the parameters of reactor
40*
41*----------------------------------------------------------------------
42* Author: Christiano D. Wetzel Guerra and Rodolfo Rodrigues
43* $Id: hidrodesalkeletion.mso 171 2007-03-02 13:06:53Z bicca $
44*--------------------------------------------------------------------*#
45
46using "types";
47
48
49#*---------------------------------------------------------------------
50* Example 6-6: in a PFR
51*--------------------------------------------------------------------*#
52
53FlowSheet pfr
54        PARAMETERS
55        NComp           as Integer              (Brief="Number of components", Default=1);
56        NReac           as Integer              (Brief="Number of reactions");
57        stoic(NComp,NReac) as Real      (Brief="Stoichiometric coefficients");
58        k(NReac)        as Real                 (Brief="Specific velocity reaction", Unit='(ft^3/lbmol)^.5/h');
59        P                       as pressure     (Brief="Pressure of reactor");
60        T                       as temperature  (Brief="Temperature of reactor");
61        R                       as Real                 (Brief="Universal gas constant", Unit='atm*ft^3/degR/lbmol', Default=0.73);
62        yo(NComp)       as fraction     (Brief="Input molar fraction");
63        vo                      as flow_vol     (Brief="Input volumetric flow");
64       
65        VARIABLES
66        F(NComp)        as flow_mol     (Brief="Molar flow", Unit='lbmol/h');
67        C(NComp)        as conc_mol     (Brief="Molar concentration", Lower=0, Unit='lbmol/ft^3');
68        Co(NComp)       as conc_mol     (Brief="Input molar concentration", Lower=0, Unit='lbmol/ft^3');
69       
70        r(NComp,NReac) as reaction_mol (Brief="Relative rate of reaction", Unit='lbmol/h/ft^3');
71        rate(NComp)     as reaction_mol (Brief="Overall rate of reaction", Unit='lbmol/h/ft^3');
72        tau                     as time_h               (Brief="Residence time", Unit='h');
73        V                       as volume               (Brief="Reactor volume", Unit='ft^3');
74
75        EQUATIONS
76        "Change time in tau"
77        tau = time;
78       
79        "Molar balance"
80        diff(C) = rate;
81       
82        "Molar flow"
83        F = vo*C;
84       
85        "Inlet molar concentration"
86        Co = yo*P/R/T;
87       
88        "Residence time"
89        tau = V/vo;
90
91        "Relative rate of reaction 1"
92        r(:,1) = stoic(:,1)*k(1)*C(1)*sqrt(C(2));
93       
94        "Relative rate of reaction 2"
95        r(:,2) = stoic(:,2)*k(2)*C(3)*sqrt(C(2));
96       
97        "Overall rate of reaction"
98        rate = sumt(r);
99       
100        SET
101        NComp = 5;      # 1: mesitylene, 2: hydrogen, 3: m-xylene,
102                                # 4: methane and 5: toluene
103       
104        NReac = 2;      # 1: M + H -> X + Me,
105                                # 2: X + H -> T + Me
106       
107        stoic(:,1) = [-1, -1,  1, 1, 0]; # M + H -> X + Me
108        stoic(:,2) = [ 0, -1, -1, 1, 1]; # X + H -> T + Me
109       
110        vo= 476*'ft^3/h';
111        T = 1500*'degR';
112        P = 35*'atm';
113        yo= [1/3, 2/3, 0.0, 0.0, 0.0];
114        k = [55.2, 30.2]*'(ft^3/lbmol)^.5/h';
115       
116        INITIAL
117        "Molar concentration"
118        C = Co;
119       
120        OPTIONS
121        TimeStep =0.01;
122        TimeEnd  =0.5;
123        TimeUnit = 'h';
124end
Note: See TracBrowser for help on using the repository browser.