| 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 245 2007-04-19 18:37:24Z 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 | EQUATIONS |
|---|
| 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 | |
|---|
| 55 | SPECIFY |
|---|
| 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 | |
|---|
| 68 | CONNECTIONS |
|---|
| 69 | |
|---|
| 70 | s1.Outlet to Reac.Inlet; |
|---|
| 71 | |
|---|
| 72 | INITIAL |
|---|
| 73 | |
|---|
| 74 | Reac.Outlet.T = 533 * 'degR'; |
|---|
| 75 | Reac.M = 100 * 'kmol'; |
|---|
| 76 | Reac.Outlet.z(1:3) = [0.0468, 0.0783, 0.874]; |
|---|
| 77 | |
|---|
| 78 | OPTIONS |
|---|
| 79 | TimeStep = 0.1; |
|---|
| 80 | TimeEnd = 5; |
|---|
| 81 | |
|---|
| 82 | end |
|---|
