[82] | 1 | #*------------------------------------------------------------------- |
---|
| 2 | * EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC. |
---|
[1] | 3 | * |
---|
[82] | 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 a isotermic batch reactor (model batch_isotermic). |
---|
[1] | 17 | *---------------------------------------------------------------------- |
---|
| 18 | * Author: Paula Bettio Staudt |
---|
| 19 | * Rafael de Pelegrini Soares |
---|
| 20 | * $Id: sample_batch.mso 574 2008-07-25 14:18:50Z rafael $ |
---|
| 21 | *-------------------------------------------------------------------*# |
---|
| 22 | |
---|
[8] | 23 | using "reactors/batch"; |
---|
| 24 | |
---|
[1] | 25 | FlowSheet SampleBatch |
---|
| 26 | PARAMETERS |
---|
[8] | 27 | NComp as Integer(Brief="Number of components"); |
---|
[1] | 28 | |
---|
| 29 | DEVICES |
---|
| 30 | feed as stream; |
---|
[313] | 31 | reac as batch_isothermic; |
---|
[1] | 32 | |
---|
| 33 | CONNECTIONS |
---|
| 34 | feed to reac.Inlet; |
---|
| 35 | |
---|
| 36 | PARAMETERS |
---|
[172] | 37 | k1 as Real(Unit='(l/mol)^3.5/min', Default=18.7); |
---|
| 38 | k2 as Real(Unit='l/mol/min', Default=0.58); |
---|
| 39 | k3 as Real(Unit='(l/mol)^2/min', Default=0.09); |
---|
| 40 | k4 as Real(Unit='(l/mol)^1.5/min', Default=0.42); |
---|
[1] | 41 | Ks as inv_conc_mol(Default=115.83); |
---|
| 42 | K as Real(Default=34.4); |
---|
| 43 | kla as flow_vol(Default=3.3); |
---|
| 44 | H as inv_conc_mol(Default=737); |
---|
| 45 | pCO2; |
---|
| 46 | |
---|
| 47 | VARIABLES |
---|
| 48 | C6 as conc_mol(Brief="Special component concentration"); |
---|
| 49 | |
---|
| 50 | SET |
---|
[8] | 51 | NComp = 5; |
---|
| 52 | reac.NReac = 5; |
---|
[1] | 53 | pCO2 = 0.9; |
---|
| 54 | |
---|
| 55 | reac.stoic(1,:) = [-2, 1, -1, -1, 0 ]; |
---|
| 56 | reac.stoic(2,:) = [-0.5, 0, 0, -1, -0.5]; |
---|
| 57 | reac.stoic(3,:) = [1, -1, 1, 0, 0 ]; |
---|
| 58 | reac.stoic(4,:) = [0, -1, 1, -2, 0 ]; |
---|
| 59 | reac.stoic(5,:) = [0, 1, -1, 0, 1 ]; |
---|
| 60 | |
---|
| 61 | EQUATIONS |
---|
| 62 | "Inlet Flow" |
---|
| 63 | feed.F = kla*(pCO2/H-reac.C(2)); |
---|
| 64 | |
---|
| 65 | "Reaction rates" |
---|
| 66 | reac.r(1)=k1 * reac.C(1)^4 * sqrt(reac.C(2)); |
---|
| 67 | reac.r(2)=k2 * reac.C(3) * reac.C(4); |
---|
| 68 | reac.r(3)=k2/K * reac.C(1) * reac.C(5); |
---|
| 69 | reac.r(4)=k3 * reac.C(1) * reac.C(4)^2; |
---|
| 70 | reac.r(5)=k4 * C6^2 * sqrt(reac.C(2)); |
---|
| 71 | |
---|
| 72 | "C6 Equilibrium" |
---|
| 73 | C6 = Ks * reac.C(1) * reac.C(4); |
---|
| 74 | |
---|
| 75 | SPECIFY |
---|
[172] | 76 | feed.T = 300*'K'; |
---|
| 77 | feed.P = 1*'atm'; |
---|
| 78 | feed.h = 0*'J/mol'; |
---|
[1] | 79 | feed.v = 0; |
---|
| 80 | feed.z = [0,1,0,0,0]; |
---|
| 81 | |
---|
[172] | 82 | reac.Vr = 1 * 'l'; |
---|
[1] | 83 | |
---|
| 84 | INITIAL |
---|
| 85 | "Initial concentration" |
---|
[172] | 86 | reac.C=[0.444, 0.00122, 0, 0.007, 0] * 'mol/l'; |
---|
[1] | 87 | |
---|
| 88 | OPTIONS |
---|
[172] | 89 | TimeStep = 100; |
---|
| 90 | TimeEnd = 180*60; |
---|
[1] | 91 | end |
---|