source: branches/gui/sample/reactors/sample_cstr.mso

Last change on this file was 585, checked in by gerson bicca, 15 years ago

fixed

  • Property svn:keywords set to Id
File size: 2.2 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* Sample file for model cstr
17*--------------------------------------------------------------------
18*
19* This sample file needs VRTherm (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22* Author: Paula B. Staudt
23* $Id: sample_cstr.mso 585 2008-07-28 23:31:29Z bicca $
24*--------------------------------------------------------------------*#
25
26using "reactors/cstr";
27
28FlowSheet cstr_PropyleneGlycol
29
30PARAMETERS
31        PP              as Plugin(Brief="Physical Properties",
32                Type = "PP",
33                Components  = ["1,2-propylene oxide", "methanol","water","1,2-propylene glycol" ],
34                LiquidModel = "PR",
35                VapourModel = "PR"
36        );
37        NComp   as Integer;
38        R                       as Real         (Default=8.3145, Unit='J/mol/K');
39
40DEVICES
41        s1      as source;
42        Reac    as cstr;
43
44SET
45        NComp                   = PP.NumberOfComponents;
46        Reac.NReac      = 1;
47        Reac.stoic              = [-1, 0, -1, 1];
48       
49        s1.ValidPhases = "Vapour-Liquid";
50        s1.CompositionBasis = "Molar";
51       
52EQUATIONS
53
54        Reac.Hr = -36400 * 'Btu/lbmol' - 7 * 'Btu/lbmol/degR' * (Reac.Outlet.T - 528 * 'degR');
55        Reac.r = 16.96e12 / 'h' * exp(- 32400 * 'Btu/lbmol'/R/Reac.Outlet.T)
56                * Reac.C(1);
57
58SPECIFY
59        s1.F = 52126 * 'mol/h';
60        s1.P = 1     * 'atm';
61        s1.T = 534.67* 'degR';
62        s1.Composition = [0.0468, 0.0783, 0.8749, 0];
63
64        Reac.Vr                 = 300   * 'gal';
65        Reac.Outlet.F   = Reac.Inlet.F;
66
67        Reac.q                  = 0     * 'J/s';
68        Reac.Outlet.P   = 1     * 'atm';
69       
70CONNECTIONS
71
72        s1.Outlet  to  Reac.Inlet;
73       
74INITIAL
75
76        Reac.Outlet.T = 533 * 'degR';
77        Reac.M = 100 * 'kmol';
78        Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874];
79
80OPTIONS
81        TimeStep        = 0.1;
82        TimeEnd         = 5;
83
84end
Note: See TracBrowser for help on using the repository browser.