11eq1

Name: 11eq1 - Stirred pot reactor process
Source: Henley, E. J. and Rosen, E. M., Material and Energy Balance
  Computations, Wiley: New York, 1969
Reference/s      
Model: 11 implicit equations, indep. variables x1, x2, x3, V, L, s11, s12, s13, s61, s62 and s63
  Lower difficulty level
  Constraints: All the variables are nonnegative
  Discontinuities: none
  Initial estimates: 1.(0.1, 0.2, 0, 500, 500, 970, 30, 0, 970, 30, 0) ;
  2. (0.1, 0.4, 0.5, 600, 1500, 970, 300, 0, 500, 900, 800)
  3. (0.3, 0.3, 0.4, 100, 100, 100, 100, 100, 100, 100, 100, 100, 100)
       
Solved by Shacham, M., POLYMATH 5.1, build 19, April 18, 2001
       
Model Eqs. Stirred pot reactor process |POLVER05_1

EXCEL FILE 

f(x1)=s61-x1*(L+V*2.523) #

TEXT FILE 

f(x2)=s62-x2*(L+V*1.57) #

POLYMATH FILE 

f(x3)=s63-x3*(L+V*0.0329) #
  f(V)=2.523*x1+1.57*x2+0.0329*x3-1 #
  f(L)=x1+x2+x3-1 #
  f(s11)=s11-970-s51 #
  f(s12)=s12-30-s52 #
  f(s13)=s13-s53 #
  f(s61)=s81+s101+s31-s61 #
  f(s62)=s82+s102+s32-s62 #
  f(s63)=s83+s103+s33-s63 #
  s102=0 #
  s21=s11/(1+0.211*1.5) #
  s71=L*x1 #
  s72=L*x2 #
  s73=L*x3 #
  s41=V*2.523*x1 #
  s42=V*1.57*x2 #
  s43=V*0.0329*x3 #
  s51=0.75*s41 #
  s52=0.75*s42 #
  s53=0.75*s43 #
  s81=0.1*s71 #
  s101=0.9*0.1*s71 #
  s82=0.5*s72 #
  s83=0.2*s73 #
  s103=0.8*0.4*s73 #
  s22=s12/(1+0.101*1.5)+0.211*s21*1.5/(1+0.101*1.5) #
  s31=s21/(1+0.44*2) #
  s23=s13+0.101*1.5*s22 #
  s32=s22/(1+0.219*2)+0.44*s31*2/(1+0.219*2) #
  s33=s23+0.219*2*s32 #
  x1(0)=0.1
  x2(0)=0.2
  x3(0)=0
  V(0)=500
  L(0)=500
  s11(0)=970
  s12(0)=30
  s13(0)=0
  s61(0)=970
  s62(0)=30
  s63(0)=0
       
Variable/function values Variable Value f(x)
  x1 0.1 793.85
  x2 0.2 -227
  x3 0 0
  V 500 -0.434
  L 500 -0.700
  s11 970 -94.613
  s12 30 -117.750
  s13 0 0
  s61 970 -568.584
  s62 30 418.787
  s63 0 209.297
  s102 0  
  s21 736.802  
  s71 50  
  s72 100  
  s73 0  
  s41 126.150  
  s42 157  
  s43 0  
  s51 94.613  
  s52 117.75  
  s53 0  
  s81 5  
  s101 4.5  
  s82 50  
  s83 0  
  s103 0  
  s22 228.570  
  s31 391.916  
  s23 34.628  
  s32 398.787  
  s33 209.297  
       
Solution Variable Value f(x)
  x1 0.16547930319775 0
  x2 0.36109556119140 0
  x3 0.4734251356109 -1.137E-13
  V 630.763337421609 0
  L 1554.563311908010 0
  s11 1167.509795711420 -1.990E-13
  s12 298.194278136889 0
  s13 7.368429217898 -1.776E-15
  s61 520.594447913209 -1.137E-13
  s62 918.938282370162 -1.137E-13
  s63 745.793919046246 0
  s102 0  
  s21 886.828557319726  
  s71 257.248053631316  
  s72 561.345911520977  
  s73 735.969346755715  
  s41 263.346394281892  
  s42 357.592370849185  
  s43 9.824572290531  
  s51 197.509795711419  
  s52 268.194278136889  
  s53 7.368429217898  
  s81 25.724805363132  
  s101 23.152324826819  
  s82 280.672955760488  
  s83 147.193869351143  
  s103 235.510190961829  
  s22 502.714300068243  
  s31 471.717317723258  
  s23 83.529645678237  
  s32 638.265326609674  
  s33 363.089858733274