source: mso/sample/miscellaneous/sample_feed_batch.mso @ 1

Last change on this file since 1 was 1, checked in by Rafael de Pelegrini Soares, 16 years ago

Initial import of the library

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 2.0 KB
Line 
1#*---------------------------------------------------------------------
2* This file is property of the author and cannot be used, copyed
3* or modified without permission.
4*
5* Copyright (C) 2002-2004  the author
6*----------------------------------------------------------------------
7* Author: Paula Bettio Staudt and Rafael de Pelegrini Soares
8* $Id: sample_feed_batch.mso 1 2006-06-20 17:33:53Z rafael $
9*----------------------------------------------------------------------
10*
11* Sample file for a feed-batch reactor.
12*
13*-------------------------------------------------------------------*#
14
15
16using "feed_batch";
17
18FlowSheet SampleBatch as feed_batch
19        PARAMETERS
20        NoComps as Integer(Brief="Number of components");
21        k1  as Real(Unit="(l/mol)^3.5/min", Default=18.7);
22        k2  as Real(Unit="l/mol/min", Default=0.58);
23        k3  as Real(Unit="(l/mol)^2/min", Default=0.09);
24        k4  as Real(Unit="(l/mol)^1.5/min", Default=0.42);
25        Ks  as inv_conc_mol(Default=115.83);
26        K   as Real(Default=34.4);
27        kla as flow_vol(Default=3.3);
28        H   as inv_conc_mol(Default=737);
29        pCO2;
30       
31        VARIABLES
32        C6 as conc_mol(Brief="Special component concentration");
33
34        SET
35        NoComps = 5;
36        NoReacs = 5;
37        stoic(1,:) = [-2,   1,  -1, -1,  0  ];
38        stoic(2,:) = [-0.5, 0,  0,  -1, -0.5];
39        stoic(3,:) = [1,   -1,  1,   0,   0 ];
40        stoic(4,:) = [0,   -1,  1,  -2,   0 ];
41        stoic(5,:) = [0,    1, -1,   0,   1 ];
42        pCO2 = 0.9;
43       
44        EQUATIONS
45        "Fin"
46        Fin.F = kla*(pCO2/H-C(2));
47        "Taxas"
48        r(1)=k1   * C(1)^4 * sqrt(C(2));
49        r(2)=k2   * C(3)   * C(4);
50        r(3)=k2/K * C(1)   * C(5);
51        r(4)=k3   * C(1)   * C(4)^2;
52        r(5)=k4   * C6^2   * sqrt(C(2));
53       
54        "C6 Equilibrium"
55        C6  = Ks * C(1) * C(4);
56       
57        SPECIFY
58        "Feed conditions"
59        Fin.T = 300*"K";
60        Fin.P = 1*"atm";
61        Fin.h = 0*"J/mol";
62        Fin.z = [0,1,0,0,0];
63
64        "Reactor volume"
65        Vr = 1 * "l";
66       
67        "Isothermic"
68        Hr = 0 * "J/mol";
69        q  = 0 * "W";
70        T  = Fin.T;
71        P  = Fin.P;     
72       
73        INITIAL
74        "Initial concentration"
75        C=[0.444, 0.00122, 0, 0.007, 0] * "mol/l";
76        "Initial energy"
77        h = 0 * "J/mol";
78
79        OPTIONS
80        integration = "index1";
81        outputLevel = "high";
82        time = [0:10:180] * "min";
83end
Note: See TracBrowser for help on using the repository browser.