|
Last change
on this file since 390 was
390,
checked in by Argimiro Resende Secchi, 15 years ago
|
|
Adding AUTO DAE examples.
|
|
File size:
1.1 KB
|
| Rev | Line | |
|---|
| [390] | 1 | using "types"; |
|---|
| 2 | |
|---|
| 3 | FlowSheet cstr_dae |
|---|
| 4 | |
|---|
| 5 | PARAMETERS |
|---|
| 6 | fi as Real; |
|---|
| 7 | beta as Real; |
|---|
| 8 | delta as Real; |
|---|
| 9 | q as Real; |
|---|
| 10 | alfa as Real; |
|---|
| 11 | S as Real; |
|---|
| 12 | csi as Real; |
|---|
| 13 | delta1 as Real; |
|---|
| 14 | delta2 as Real; |
|---|
| 15 | gama as Real; |
|---|
| 16 | x1f as Real; |
|---|
| 17 | x2f as Real; |
|---|
| 18 | x3f as Real; |
|---|
| 19 | x4f as Real; |
|---|
| 20 | qc as Real; |
|---|
| 21 | |
|---|
| 22 | VARIABLES |
|---|
| 23 | u1 as Real; |
|---|
| 24 | u2 as Real; |
|---|
| 25 | u3 as Real; |
|---|
| 26 | u4 as Real; |
|---|
| 27 | u5 as Real; |
|---|
| 28 | u6 as Real; |
|---|
| 29 | |
|---|
| 30 | SET |
|---|
| 31 | fi = 0.04; |
|---|
| 32 | beta = 8; |
|---|
| 33 | delta = 0.78; |
|---|
| 34 | q = 1; |
|---|
| 35 | alfa = 0.19; |
|---|
| 36 | S = 1.015; |
|---|
| 37 | csi = 0.32; |
|---|
| 38 | delta1 = 10; |
|---|
| 39 | delta2 = 0.952; |
|---|
| 40 | gama = 27.85; |
|---|
| 41 | x1f = 1; |
|---|
| 42 | x2f = 0; |
|---|
| 43 | x3f = 0; |
|---|
| 44 | x4f = -1; |
|---|
| 45 | qc = 0; |
|---|
| 46 | |
|---|
| 47 | EQUATIONS |
|---|
| 48 | |
|---|
| 49 | q*(x1f-u1)-u1*u5*fi - diff(u1)*'s' = 0; |
|---|
| 50 | q*(x2f-u2)-u2*fi*S*u6+u1*fi*u5 - diff(u2)*'s' = 0; |
|---|
| 51 | q*(x3f-u3)+delta*(u4-u3)+beta*fi*(u1*u5+alfa*u2*u6*S) - diff(u3)*'s' = 0; |
|---|
| 52 | delta1*(qc*(x4f-u4)+delta*delta2*(u3-u4)) - diff(u4)*'s' = 0; |
|---|
| 53 | -u5 + exp(u3/(1+(u3/gama))) = 0; |
|---|
| 54 | -u6 + exp(csi*u3/(1+(u3/gama))) = 0; |
|---|
| 55 | |
|---|
| 56 | GUESS |
|---|
| 57 | u1 = 0.6419; |
|---|
| 58 | u2 = 0.3270; |
|---|
| 59 | u3 = 2.9101; |
|---|
| 60 | u4 = 2.9101; |
|---|
| 61 | u5 = 13.9403; |
|---|
| 62 | u6 = 2.3236; |
|---|
| 63 | |
|---|
| 64 | INITIAL |
|---|
| 65 | u1 = 0.6; |
|---|
| 66 | u2 = 0.3; |
|---|
| 67 | u3 = 3; |
|---|
| 68 | u4 = 3; |
|---|
| 69 | |
|---|
| 70 | OPTIONS |
|---|
| 71 | TimeEnd = 10; |
|---|
| 72 | TimeStep = 0.2; |
|---|
| 73 | Dynamic = false; |
|---|
| 74 | end |
|---|
Note: See
TracBrowser
for help on using the repository browser.
