source: mso/sample/miscellaneous/sample_batch.mso @ 8

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

Fixed some samples and removed duplicated ones

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 2.1 KB
RevLine 
[1]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
8*         Rafael de Pelegrini Soares
9* $Id: sample_batch.mso 8 2006-08-26 18:05:18Z rafael $
10*----------------------------------------------------------------------
11*
12* Sample file for a isotermic batch reactor.
13*
14*-------------------------------------------------------------------*#
15
[8]16using "reactors/batch";
17 
[1]18FlowSheet SampleBatch
19        PARAMETERS
[8]20        NComp as Integer(Brief="Number of components");
[1]21
22        DEVICES
23        feed as stream;
24        reac as batch_isotermic;
25       
26        CONNECTIONS
27        feed to reac.Inlet;
28       
29        PARAMETERS
30        k1 as Real(Unit="(l/mol)^3.5/min", Default=18.7);
31        k2 as Real(Unit="l/mol/min", Default=0.58);
32        k3 as Real(Unit="(l/mol)^2/min", Default=0.09);
33        k4 as Real(Unit="(l/mol)^1.5/min", Default=0.42);
34        Ks as inv_conc_mol(Default=115.83);
35        K  as Real(Default=34.4);
36        kla as flow_vol(Default=3.3);
37        H as inv_conc_mol(Default=737);
38        pCO2;
39       
40        VARIABLES
41        C6 as conc_mol(Brief="Special component concentration");
42
43        SET
[8]44        NComp = 5;
45        reac.NReac = 5;
[1]46        pCO2 = 0.9;
47       
48        reac.stoic(1,:) = [-2,   1,  -1, -1,  0  ];
49        reac.stoic(2,:) = [-0.5, 0,  0,  -1, -0.5];
50        reac.stoic(3,:) = [1,   -1,  1,   0,   0 ];
51        reac.stoic(4,:) = [0,   -1,  1,  -2,   0 ];
52        reac.stoic(5,:) = [0,    1, -1,   0,   1 ];
53
54        EQUATIONS
55        "Inlet Flow"
56        feed.F = kla*(pCO2/H-reac.C(2));
57       
58        "Reaction rates"
59        reac.r(1)=k1   * reac.C(1)^4 * sqrt(reac.C(2));
60        reac.r(2)=k2   * reac.C(3)   * reac.C(4);
61        reac.r(3)=k2/K * reac.C(1)   * reac.C(5);
62        reac.r(4)=k3   * reac.C(1)   * reac.C(4)^2;
63        reac.r(5)=k4   * C6^2   * sqrt(reac.C(2));
64       
65        "C6 Equilibrium"
66        C6  = Ks * reac.C(1) * reac.C(4);
67       
68        SPECIFY
69        "Feed conditions"
70        feed.T = 300*"K";
71        feed.P = 1*"atm";
72        feed.h = 0*"J/mol";
73        feed.v = 0;
74        feed.z = [0,1,0,0,0];
75
76        "Reactor volume"
77        reac.Vr = 1 * "l";
78       
79        INITIAL
80        "Initial concentration"
81        reac.C=[0.444, 0.00122, 0, 0.007, 0] * "mol/l";
82
83        OPTIONS
84        outputLevel = "high";
85        time = [0:100:180*60];
86end
Note: See TracBrowser for help on using the repository browser.