source: trunk/sample/reactors/sample_stoic.mso @ 1006

Last change on this file since 1006 was 913, checked in by Argimiro Resende Secchi, 13 years ago

Checking new EML.

File size: 3.6 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* Samples of a stoichiometric 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/stoic";
34
35
36#*----------------------------------------------------------------------
37* only vapour phase and
38* the extents of reactions are known
39*---------------------------------------------------------------------*#
40FlowSheet sample_stoic_extent_vap
41        PARAMETERS
42        PP              as Plugin(Brief="External physical properties", Type="PP",
43        Components = ["1,3-butadiene", "hydrogen", "1-butene",
44                                        "cis-2-butene", "trans-2-butene"],
45        LiquidModel = "IdealLiquid",
46        VapourModel = "Ideal");
47        NComp   as Integer;
48       
49        DEVICES
50        Fin                     as source;
51        R                       as stoic_extent_vap;
52
53        CONNECTIONS
54        Fin.Outlet      to R.Inlet;
55       
56        SET
57        NComp = PP.NumberOfComponents;
58       
59        R.NReac = 3;
60       
61        R.stoic(:,1) = [-1, -1, 1, 0, 0]; # (1)
62        R.stoic(:,2) = [-1, -1, 0, 1, 0]; # (2)
63        R.stoic(:,3) = [-1, -1, 0, 0, 1]; # (3)
64       
65        SPECIFY
66        Fin.F = 1000*'kmol/h';
67        Fin.Composition = [0.25, 0.75, 0, 0, 0];
68        Fin.P = 10*'kgf/cm^2';
69        Fin.T = 500*'K';
70
71        R.Outlet.F = 1000*'kmol/h'*sqrt(R.Tank.Level/'m');
72        R.Outlet.P = 10*'kgf/cm^2';
73        R.Q = 0*'W';
74       
75        R.extent = [85, 65, 75]*'kmol/h';
76        R.Tank.Across = 150*'cm^2';
77        R.Tank.L = 350*'cm';
78       
79        INITIAL
80    R.Outletm.T = 500*'K';
81        R.M = [25, 75, 0, 0, 0]*'kmol';
82       
83        OPTIONS
84        TimeStep = 5;
85        TimeEnd = 150;
86end
87
88
89#*----------------------------------------------------------------------
90* only vapour phase and
91* the conversion of a key component is known
92*---------------------------------------------------------------------*#
93FlowSheet sample_stoic_conv_vap
94        PARAMETERS
95        PP              as Plugin(Brief="External physical properties", Type="PP",
96        Components = ["1,3-butadiene", "hydrogen", "1-butene",
97                                                "cis-2-butene", "trans-2-butene"],
98        LiquidModel = "IdealLiquid",
99        VapourModel = "Ideal");
100        NComp   as Integer;
101       
102        DEVICES
103        Fin                     as source;
104        R                       as stoic_conv_vap;
105
106        CONNECTIONS
107        Fin.Outlet      to R.Inlet;
108
109        SET
110        NComp = PP.NumberOfComponents;
111       
112        R.NReac = 3;
113        R.KComp = 1; # 1,3-butadiene
114
115        R.stoic(:,1) = [-1, -1, 1, 0, 0]; # (1)
116        R.stoic(:,2) = [-1, -1, 0, 1, 0]; # (2)
117        R.stoic(:,3) = [-1, -1, 0, 0, 1]; # (3)
118       
119        SPECIFY
120        Fin.F = 1000*'kmol/h';
121        Fin.Composition = [0.25, 0.75, 0, 0, 0];
122        Fin.P = 10*'kgf/cm^2';
123        Fin.T = 500*'K';
124       
125        R.Outlet.F = 1000*'kmol/h'*sqrt(R.Tank.Level/'m');
126        R.Outlet.P = 10*'kgf/cm^2';
127        R.Q = 0*'W';
128
129        R.kconv = 0.9;
130        R.Tank.Across = 150*'cm^2';
131        R.Tank.L = 350*'cm';
132
133        OPTIONS
134        Dynamic = false;
135end
Note: See TracBrowser for help on using the repository browser.