| Name: |
Twoeq6 - Conversion in a chemical
reactor |
| Source: |
Shacham, M., 1986, Int.
J. numerical Meth. Engng, 23, 1455 |
| Reference/s |
Shacham, M. 1989, Chem.
Engng. Sci., 44(7), 1495 |
| |
|
|
|
| Model: |
2 implicit equation, indep.
variable names: x1, x2 |
| |
Higher difficulty level |
| |
Constraints: 0<=x1<=1 |
| |
Discontinuities: undefined
for x1=1, x2=0, (x1>1 & x2>0) |
| |
Initial estimates: 1. x1=0.9,
x2=0.5; 2. x1=0.5, x2=0.5; 3.x1=0.4, x2=0.05; |
| |
4. x1=0.6, x2=0.1; |
| |
|
|
|
| Solved by |
Shacham, M., POLYMATH 5.1,
build 19, April 2, 2001 |
| |
|
|
|
| Model Eqs. |
Conversion in a Chemical Reactor
|POLVER05_3 |
|
EXCEL
FILE |
f(x1) = x1/(1-x1)-5*ln(0.4*(1-x1)/x2)+4.45977
# |
|
TEXT
FILE |
f(x2) = x2-(0.4-0.5*x1) # |
|
POLYMATH
FILE |
x1(0)=0.9 |
| |
x2(0)=0.5 |
| |
|
|
|
| Variable/function values |
Variable |
Value |
f(x) |
| |
x1 |
0.9 |
2.6088E+01 |
| |
x2 |
0.5 |
5.5000E-01 |
| |
|
|
|
| Solution |
Variable |
Value |
f(x) |
| |
x1 |
0.7573962468236 |
6.2172E-15 |
| |
x2 |
0.0213018765882 |
-2.0817E-17 |
| |
|
|
|
| Infeasible Solution |
Variable |
Value |
f(x) |
| |
x1 |
1.0989839337750 |
4.8850E-13 |
| |
x2 |
-0.1494919668876 |
0.0000E+00 |
| |
|
|
|
| Additional information |
Modified version of problem
Oneeq7. Most programs won't |
| |
converge from the 1st initial
guess and will converge to the |
| |
infeasible solution from the
4rt initial guess |