[578] | 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 585 2008-07-28 23:31:29Z bicca $ |
---|
| 24 | *--------------------------------------------------------------------*# |
---|
| 25 | |
---|
| 26 | using "reactors/cstr"; |
---|
| 27 | |
---|
| 28 | FlowSheet cstr_PropyleneGlycol |
---|
| 29 | |
---|
| 30 | PARAMETERS |
---|
| 31 | PP as Plugin(Brief="Physical Properties", |
---|
| 32 | Type = "PP", |
---|
| 33 | Components = ["1,2-propylene oxide", "methanol","water","1,2-propylene glycol" ], |
---|
| 34 | LiquidModel = "PR", |
---|
| 35 | VapourModel = "PR" |
---|
| 36 | ); |
---|
| 37 | NComp as Integer; |
---|
| 38 | R as Real (Default=8.3145, Unit='J/mol/K'); |
---|
| 39 | |
---|
| 40 | DEVICES |
---|
| 41 | s1 as source; |
---|
| 42 | Reac as cstr; |
---|
| 43 | |
---|
| 44 | SET |
---|
| 45 | NComp = PP.NumberOfComponents; |
---|
| 46 | Reac.NReac = 1; |
---|
| 47 | Reac.stoic = [-1, 0, -1, 1]; |
---|
| 48 | |
---|
| 49 | s1.ValidPhases = "Vapour-Liquid"; |
---|
[585] | 50 | s1.CompositionBasis = "Molar"; |
---|
[578] | 51 | |
---|
| 52 | EQUATIONS |
---|
| 53 | |
---|
| 54 | Reac.Hr = -36400 * 'Btu/lbmol' - 7 * 'Btu/lbmol/degR' * (Reac.Outlet.T - 528 * 'degR'); |
---|
| 55 | Reac.r = 16.96e12 / 'h' * exp(- 32400 * 'Btu/lbmol'/R/Reac.Outlet.T) |
---|
| 56 | * Reac.C(1); |
---|
| 57 | |
---|
| 58 | SPECIFY |
---|
| 59 | s1.F = 52126 * 'mol/h'; |
---|
| 60 | s1.P = 1 * 'atm'; |
---|
| 61 | s1.T = 534.67* 'degR'; |
---|
| 62 | s1.Composition = [0.0468, 0.0783, 0.8749, 0]; |
---|
| 63 | |
---|
| 64 | Reac.Vr = 300 * 'gal'; |
---|
| 65 | Reac.Outlet.F = Reac.Inlet.F; |
---|
| 66 | |
---|
| 67 | Reac.q = 0 * 'J/s'; |
---|
| 68 | Reac.Outlet.P = 1 * 'atm'; |
---|
| 69 | |
---|
| 70 | CONNECTIONS |
---|
| 71 | |
---|
| 72 | s1.Outlet to Reac.Inlet; |
---|
| 73 | |
---|
| 74 | INITIAL |
---|
| 75 | |
---|
| 76 | Reac.Outlet.T = 533 * 'degR'; |
---|
| 77 | Reac.M = 100 * 'kmol'; |
---|
| 78 | Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874]; |
---|
| 79 | |
---|
| 80 | OPTIONS |
---|
| 81 | TimeStep = 0.1; |
---|
| 82 | TimeEnd = 5; |
---|
| 83 | |
---|
| 84 | end |
---|