source: branches/newlanguage/sample/reactors/sample_cstr.mso @ 169

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

updated cstr reactor model for the new language

  • 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 169 2007-03-01 18:09:21Z bicca $
24*--------------------------------------------------------------------*#
25
26using "reactors/cstr";
27
28FlowSheet cstr_PropyleneGlycol
29
30PARAMETERS
31        PP              as Plugin       (Brief="Physical Properties",File="vrpp");
32        NComp   as Integer;
33        R                       as Real         (Default=8.3145, Unit='J/mol/K');
34
35DEVICES
36        s1      as source;
37        Reac    as cstr;
38
39SET
40
41        PP.Components   = ["1,2-propylene oxide", "methanol","water","sulfuric acid", "1,2-propylene glycol" ];
42        PP.LiquidModel  = "PR";
43        PP.VapourModel  = "PR";
44
45        NComp                   = PP.NumberOfComponents;
46        Reac.NReac      = 1;
47        Reac.stoic              = [-1, 0, -1, 0, 1];
48       
49EQUATIONS
50
51#       Reac.Hr = -36400 * "Btu/lbmol" - 7 * "Btu/lbmol/degR" * (Reac.Fout.T - 528 * "degR");
52#       Reac.r = 16.96e12 / "h" * exp(- 32400 * "Btu/lbmol"/R/Reac.Fout.T)
53#               * Reac.C(1);
54        Reac.Hr = 0 * 'Btu/lbmol';
55        Reac.r  = 0 * 'mol/l/s';
56
57        SPECIFY
58        s1.Outlet.F = 52126 * 'mol/h';
59        s1.Outlet.P = 1     * 'atm';
60        s1.Outlet.T = 534.67* 'degR';
61        s1.Outlet.z = [0.0468, 0.0783, 0.874, 0.0009, 0];
62        s1.Outlet.v = 0.0;
63
64        Reac.Vr                 = 300   * 'gal';
65        Reac.Outlet.F   = Reac.Inlet.F;
66
67"Adiabatic"
68        Reac.q                  = 0     * 'J/s';
69        Reac.Outlet.P   = 1     * 'atm';
70       
71CONNECTIONS
72
73        s1.Outlet  to  Reac.Inlet;
74       
75INITIAL
76
77        Reac.Outlet.T = Reac.Inlet.T;
78        Reac.Outlet.z = Reac.Inlet.z;
79
80OPTIONS
81        TimeStep        = 0.1;
82        TimeEnd         = 5;
83
84end
Note: See TracBrowser for help on using the repository browser.