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

Last change on this file since 208 was 208, checked in by Argimiro Resende Secchi, 16 years ago

Fix more sample problems to new language.

  • Property svn:keywords set to Id
File size: 2.1 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 208 2007-03-15 05:18:39Z arge $
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","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, 1];
48       
49EQUATIONS
50
51        Reac.Hr = -36400 * 'Btu/lbmol' - 7 * 'Btu/lbmol/degR' * (Reac.Outlet.T - 528 * 'degR');
52        Reac.r = 16.96e12 / 'h' * exp(- 32400 * 'Btu/lbmol'/R/Reac.Outlet.T)
53                * Reac.C(1);
54
55SPECIFY
56        s1.Outlet.F = 52126 * 'mol/h';
57        s1.Outlet.P = 1     * 'atm';
58        s1.Outlet.T = 534.67* 'degR';
59        s1.Outlet.z = [0.0468, 0.0783, 0.8749, 0];
60
61        Reac.Vr                 = 300   * 'gal';
62        Reac.Outlet.F   = Reac.Inlet.F;
63
64        "Adiabatic"
65        Reac.q                  = 0     * 'J/s';
66        Reac.Outlet.P   = 1     * 'atm';
67       
68CONNECTIONS
69
70        s1.Outlet  to  Reac.Inlet;
71       
72INITIAL
73
74        Reac.Outlet.T = 530 * 'degR';
75        Reac.M = 100 * 'kmol';
76        Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874];
77
78OPTIONS
79        TimeStep        = 0.1;
80        TimeEnd         = 5;
81
82end
Note: See TracBrowser for help on using the repository browser.