| Name: |
Sixeq2b - Steady state solution
for reaction rate equations - 2nd parameter set |
| Source: |
Shacham, M., pp. 891-923 in
A. W. Westerberg and H. H. Chien (Eds), |
| |
Proc. of FOCAPD 2, CACHE Publications,
Ann Arbor, Michigan, 1984. |
| Reference/s |
Bullard, L. G. and Biegler,
L. T., Computers Chem. Engng. 15(4), 239 (1991). |
| |
|
|
|
| Model: |
6 implicit equations, indep.
variables x1 to x6 |
| |
Average difficulty level |
| |
Constraints: xi are nonnegative
for I=1,2,…6 |
| |
Discontinuities: none |
| |
Initial estimates: 1. (0.99,
0.05, 0.05, 0.99, 0.05, 0); 2.(0.05, 0.99, 0.05, 0.05, 0.99,
0); |
| |
|
|
|
| Solved by |
Shacham, M., POLYMATH 5.1,
build 19, April 9, 2001 |
| |
|
|
|
| Model Eqs. |
Steady state solution for
reaction rate equations - 2nd parameter set |POLVER05_1 |
|
EXCEL
FILE |
f(x1) = 1-x1-k1*x1*x6+kr1*x4
# |
|
TEXT
FILE |
f(x2) = 1-x2-k2*x2*x6+kr2*x5
# |
|
POLYMATH
FILE |
f(x3) = -x3+2*k3*x4*x5 # |
| |
f(x4) = k1*x1*x6-kr1*x4-k3*x4*x5
# |
| |
f(x5) = 1.5*(k2*x2*x6-kr2*x5)-k3*x4*x5
# |
| |
f(x6) = 1-x4-x5-x6 # |
| |
k1 = 17.721 # |
| |
k2 = 3.483 # |
| |
kr1 = 0.118 # |
| |
kr2 = 0.033 # |
| |
k3 = 505.051 # |
| |
x1(0)=.99 |
| |
x2(0)=.05 |
| |
x3(0)=.05 |
| |
x4(0)=.99 |
| |
x5(0)=.05 |
| |
x6(0)=0 |
| |
|
|
|
| Variable/function values |
Variable |
Value |
f(x) |
| |
x1 |
0.99 |
0.127 |
| |
x2 |
0.05 |
0.952 |
| |
x3 |
0.05 |
49.950 |
| |
x4 |
0.99 |
-25.117 |
| |
x5 |
0.05 |
-25.002 |
| |
x6 |
0 |
-0.04 |
| |
k1 |
17.721 |
|
| |
k2 |
3.483 |
|
| |
kr1 |
0.118 |
|
| |
kr2 |
0.033 |
|
| |
k3 |
505.051 |
|
| |
|
|
|
| Solution |
Variable |
Value |
f(x) |
| |
x1 |
0.9499424500947 |
0 |
| |
x2 |
0.9666283000631 |
-4.423E-17 |
| |
x3 |
0.1001150998106 |
-1.388E-17 |
| |
x4 |
0.9899868097778 |
0 |
| |
x5 |
0.0001001163356 |
6.939E-18 |
| |
x6 |
0.0099130738866 |
2.776E-17 |
| |
|
|
|
| Infeasible solution |
Variable |
Value |
f(x) |
| |
x1 |
1.0698344865210 |
1.554E-17 |
| |
x2 |
1.0465563243480 |
6.939E-18 |
| |
x3 |
-0.1396689730426 |
0 |
| |
x4 |
-0.0001377456645 |
0 |
| |
x5 |
1.0038221407600 |
0 |
| |
x6 |
-0.0036843950958 |
4.597E-17 |
| |
|
|
|
| Additional information |
Most algorithms converge to
the infeasible solution from |
| |
the second initial estimate |