Last change
on this file since 832 was
720,
checked in by gerson bicca, 14 years ago
|
fixed sample reconciliation
|
File size:
1.2 KB
|
Rev | Line | |
---|
[720] | 1 | using "stage_separators/flash"; |
---|
| 2 | |
---|
| 3 | # exemplo baseado nos dados do artigo do Gani 1985. |
---|
| 4 | FlowSheet flash_Test |
---|
| 5 | PARAMETERS |
---|
| 6 | PP as Plugin(Brief="Physical Properties", |
---|
| 7 | Type="PP", |
---|
| 8 | Components = ["1,3-butadiene", "isobutene", "n-pentane", |
---|
| 9 | "1-pentene", "1-hexene", "benzene"], |
---|
| 10 | LiquidModel = "PR", |
---|
| 11 | VapourModel = "PR" |
---|
| 12 | ); |
---|
| 13 | NComp as Integer; |
---|
| 14 | |
---|
| 15 | VARIABLES |
---|
| 16 | Q as energy_source (Brief="Heat supplied"); |
---|
| 17 | |
---|
| 18 | SET |
---|
| 19 | NComp = PP.NumberOfComponents; |
---|
| 20 | |
---|
| 21 | DEVICES |
---|
| 22 | fl as flash_steady; |
---|
| 23 | s1 as source; |
---|
| 24 | |
---|
| 25 | CONNECTIONS |
---|
| 26 | s1.Outlet to fl.Inlet; |
---|
| 27 | Q.OutletQ to fl.InletQ; |
---|
| 28 | |
---|
| 29 | SPECIFY |
---|
| 30 | s1.F = 496.3 * 'kmol/h'; |
---|
| 31 | s1.T = 338 * 'K'; |
---|
| 32 | s1.P = 507.1 * 'kPa'; |
---|
| 33 | s1.Composition = [0.2379,0.3082,0.09959,0.1373,0.08872,0.1283]; |
---|
| 34 | |
---|
| 35 | fl.OutletLiquid.P = 2.5 * 'atm'; |
---|
| 36 | |
---|
| 37 | #Q.OutletQ.Q = 0 * 'kJ/h'; |
---|
| 38 | fl.OutletLiquid.T = 315.06 * 'K'; |
---|
| 39 | |
---|
| 40 | OPTIONS |
---|
| 41 | Dynamic = false; |
---|
| 42 | end |
---|
| 43 | |
---|
| 44 | Reconciliation flash_rec as flash_Test |
---|
| 45 | |
---|
| 46 | RECONCILE |
---|
| 47 | s1.F; fl.OutletLiquid.F; fl.OutletVapour.F; |
---|
| 48 | |
---|
| 49 | FREE |
---|
| 50 | s1.F; |
---|
| 51 | |
---|
| 52 | EXPERIMENTS |
---|
| 53 | # FILE WEIGTH |
---|
| 54 | "flash.dat" 1; |
---|
| 55 | |
---|
| 56 | |
---|
| 57 | OPTIONS |
---|
| 58 | Filter = "mean"; |
---|
| 59 | Significance = 0.95; |
---|
| 60 | ObjectiveFunction = "Normal"; |
---|
| 61 | |
---|
| 62 | GrossErrorTests( |
---|
| 63 | Global = true, |
---|
| 64 | Nodal = true, |
---|
| 65 | Measurements = true |
---|
| 66 | ); |
---|
| 67 | |
---|
| 68 | NLPSolver( |
---|
| 69 | MaxIterations=1000, |
---|
| 70 | File = "complex" |
---|
| 71 | #File = "ipopt_emso" |
---|
| 72 | ); |
---|
| 73 | Dynamic = false; |
---|
| 74 | end |
---|
| 75 | |
---|
Note: See
TracBrowser
for help on using the repository browser.