source: trunk/sample/reactors/sample_equil.mso @ 988

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

Checking new EML.

File size: 2.4 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 equilibrium reactor
17*----------------------------------------------------------------------
18*
19* This sample file needs VRTherm (www.vrtech.com.br) to run.
20*
21*----------------------------------------------------------------------
22*
23*   This is a sample based on a problem of the book "Introduction to
24* Chemical Engineering Thermodynamics" by Smith, Van Ness and Abbott,
25* 5th ed., 1996.
26*
27*----------------------------------------------------------------------
28* Author: Rodolfo Rodrigues
29* $Id$
30*--------------------------------------------------------------------*#
31
32using "reactors/equil";
33
34
35#*---------------------------------------------------------------------
36* Resolved example 15.5 - p.527 (Brazilian edition)
37*
38* System in gas fase contents CO, H2O, CO2, and H2.
39*
40*       Co + H2O -> CO2 + H2
41*
42* Initial conditions
43*   1bar, 1100K, 1mol of CO, and 1mol of H2O
44*
45* Specify
46*       Outlet.T
47*       Outlet.P
48*
49* Return
50*       Outlet.z
51----------------------------------------------------------------------*#
52
53
54FlowSheet sample_equil_vap
55        PARAMETERS
56        PP              as Plugin(Brief="External physical properties", Type="PP",
57                        Components = ["carbon monoxide","water","carbon dioxide","hydrogen"],
58                        LiquidModel = "IdealLiquid",
59                        VapourModel = "Ideal");
60        NComp   as Integer;
61        NReac   as Integer;
62       
63        DEVICES
64        Fin     as source;
65        R               as equil_vap;
66
67        CONNECTIONS
68        Fin.Outlet      to R.Inlet;
69       
70        SET
71        NComp = PP.NumberOfComponents;
72       
73        R.NReac = 1;
74        R.stoic(:,1) = [-1.0, -1.0, 1.0, 1.0];
75       
76        SPECIFY
77        Fin.F = 2*'mol/s';
78        Fin.Composition= [0.5, 0.5, 0.0, 0.0];
79        Fin.P = 1*'atm';
80        Fin.T = 1100*'K';
81       
82        R.Outlet.P = 1*'bar';
83        R.Outlet.T = Fin.Outlet.T;
84        R.Outlet.F = 2*'mol/s'*sqrt(R.Tank.Level/'m');
85
86        R.Tank.Across = 5*'m^2';
87        R.Tank.L = 1.5*'m';
88       
89        INITIAL
90    R.Outletm.T = 1100*'K';
91        R.M = [1.0, 1.0, 0.0, 0.0]*'mol';
92
93        OPTIONS
94        Dynamic = false;
95end
Note: See TracBrowser for help on using the repository browser.