| 1 | #*--------------------------------------------------------------------- |
|---|
| 2 | * This file is property of the author and cannot be used, copyed |
|---|
| 3 | * or modified without permission. |
|---|
| 4 | * |
|---|
| 5 | * Copyright (C) 2004 the author |
|---|
| 6 | *---------------------------------------------------------------------- |
|---|
| 7 | * Author: Paula B. Staudt |
|---|
| 8 | * $Id: sample_cstr_PG.mso 6 2005-01-11 17:06:06Z paula $ |
|---|
| 9 | *--------------------------------------------------------------------*# |
|---|
| 10 | |
|---|
| 11 | using "reactors/cstr"; |
|---|
| 12 | |
|---|
| 13 | FlowSheet cstr_PropyleneGlycol |
|---|
| 14 | PARAMETERS |
|---|
| 15 | PP as CalcObject(Brief="Physical Properties",File="vrpp"); |
|---|
| 16 | NComp as Integer; |
|---|
| 17 | R as Real(Default=8.3145, Unit="J/mol/K"); |
|---|
| 18 | |
|---|
| 19 | DEVICES |
|---|
| 20 | s1 as stream_therm; |
|---|
| 21 | Reac as cstr; |
|---|
| 22 | |
|---|
| 23 | SET |
|---|
| 24 | PP.Components = ["1,2-propylene oxide", "methanol", |
|---|
| 25 | "water", "sulfuric acid", "1,2-propylene glycol" ]; |
|---|
| 26 | PP.LiquidModel = "PR"; |
|---|
| 27 | PP.VapourModel = "PR"; |
|---|
| 28 | |
|---|
| 29 | NComp = PP.NumberOfComponents; |
|---|
| 30 | Reac.NReac = 1; |
|---|
| 31 | Reac.stoic = [-1, 0, -1, 0, 1]; |
|---|
| 32 | |
|---|
| 33 | EQUATIONS |
|---|
| 34 | # Reac.Hr = -36400 * "Btu/lbmol" - 7 * "Btu/lbmol/degR" * (Reac.Fout.T - 528 * "degR"); |
|---|
| 35 | # Reac.r = 16.96e12 / "h" * exp(- 32400 * "Btu/lbmol"/R/Reac.Fout.T) |
|---|
| 36 | # * Reac.C(1); |
|---|
| 37 | Reac.Hr = 0 * "Btu/lbmol"; |
|---|
| 38 | Reac.r = 0 * "mol/l/s"; |
|---|
| 39 | |
|---|
| 40 | SPECIFY |
|---|
| 41 | s1.F = 52126 * "mol/h"; |
|---|
| 42 | s1.P = 1 * "atm"; |
|---|
| 43 | s1.T = 534.67* "degR"; |
|---|
| 44 | s1.z = [0.0468, 0.0783, 0.874, 0.0009, 0]; |
|---|
| 45 | s1.v = 0.0; |
|---|
| 46 | |
|---|
| 47 | Reac.Vr = 300 * "gal"; |
|---|
| 48 | Reac.Outlet.F= Reac.Inlet.F; |
|---|
| 49 | "Adiabatic" |
|---|
| 50 | Reac.q = 0 * "J/s"; |
|---|
| 51 | Reac.Outlet.P= 1 * "atm"; |
|---|
| 52 | |
|---|
| 53 | CONNECTIONS |
|---|
| 54 | s1 to Reac.Inlet; |
|---|
| 55 | |
|---|
| 56 | INITIAL |
|---|
| 57 | Reac.Outlet.T = Reac.Inlet.T; |
|---|
| 58 | Reac.Outlet.z = Reac.Inlet.z; |
|---|
| 59 | |
|---|
| 60 | OPTIONS |
|---|
| 61 | time = [0:0.1:5]*"h"; |
|---|
| 62 | end |
|---|
