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