source: trunk/sample/reactors/sample_yield.mso @ 650

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

Fixing convergence problems!

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 of an yield reactor
17*----------------------------------------------------------------------
18*
19* This sample file needs VRTherm (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22*
23* There are three parallel reactions involved:
24*       (1) C4H6(1,3-butadiene) + H2 --> C4H8(1-butene)
25*       (2) C4H6(1,3-butadiene) + H2 --> C4H8(cis-2-butene)
26*       (3) C4H6(1,3-butadiene) + H2 --> C4H8(trans-2-butene)
27*
28*----------------------------------------------------------------------
29* Author: Rodolfo Rodrigues
30* $Id$
31*--------------------------------------------------------------------*#
32
33using "reactors/yield";
34
35
36FlowSheet sample_yield
37        PARAMETERS
38        PP      as Plugin(Brief="Physical properties",
39                Type="PP",
40                Components = ["1,3-butadiene","hydrogen","1-butene",
41                                                "cis-2-butene","trans-2-butene"],
42                LiquidModel = "PR",
43                VapourModel = "PR"
44        );
45        NComp   as Integer;
46       
47        DEVICES
48        Fin             as source;
49        R                       as yield_vap;
50
51        CONNECTIONS
52        Fin.Outlet      to R.Inlet;
53       
54        SET
55        NComp = PP.NumberOfComponents;
56        R.NReac = 3;
57        R.KComp = 2;
58       
59        SPECIFY
60        Fin.F = 1000*'kmol/h';
61        Fin.Composition = [0.25, 0.69, 0.03, 0.02, 0.01];
62        Fin.P = 10*'kgf/cm^2';
63        Fin.T = 500*'K';
64       
65        R.Outlet.F = 300*'kmol/h'*sqrt(R.Tank.Level/'m');
66        R.Outlet.P = 10*'kgf/cm^2';
67        R.Outlet.T = Fin.Outlet.T;
68        R.Q = 0*'W';   
69       
70        R.yield = [0.6, 1, 0.4, 0.5, 0.05];
71        R.Tank.Across = 50000*'cm^2';
72        R.Tank.L = 10000*'cm';
73
74        INITIAL
75        R.Outletm.T = 500*'K';
76        R.M = [75, 50, 10, 5, 2]*'kmol';
77       
78        OPTIONS
79#       Dynamic = false;
80        TimeStep = 50;
81        TimeEnd = 600;
82end
Note: See TracBrowser for help on using the repository browser.