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