| Name: |
Sixeq2c - Steady state solution
for reaction rate equations - 3rd 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 - 3rd 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 = 6.966 # |
| |
kr1 = 0.118 # |
| |
kr2 = 333.333 # |
| |
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 |
17.617 |
| |
x3 |
0.05 |
49.950 |
| |
x4 |
0.99 |
-25.117 |
| |
x5 |
0.05 |
-50.000 |
| |
x6 |
0 |
-0.04 |
| |
k1 |
17.721 |
|
| |
k2 |
6.966 |
|
| |
kr1 |
0.118 |
|
| |
kr2 |
333.333 |
|
| |
k3 |
505.051 |
|
| |
|
|
|
| Solution |
Variable |
Value |
f(x) |
| |
x1 |
0.9499356376871 |
2.776E-17 |
| |
x2 |
0.9666237584581 |
-6.939E-17 |
| |
x3 |
0.1001287246257 |
1.388E-17 |
| |
x4 |
0.9899863241315 |
-6.939E-18 |
| |
x5 |
0.0001001299968 |
1.388E-17 |
| |
x6 |
0.0099135458718 |
-3.990E-17 |
| |
|
|
|
| Additional information |
The solution here is different
only after the 4th significant |
| |
digit from the solution of
Sixeq2b. No infeasible solution |
| |
was found in this case |