PFR Runaway Reactions FB |POLVER05_3
d(p)/d(t) = -rate #
d(T)/d(t) = (Qg-Qr)/3600 #
Qg =     B*p*k     #
Qr =     C*(T-Tw)     #
Ratio =     A*p*k/(C*(T-Tw)-B*k*p)     #
Tw =     625     #
k =     exp(-(a/T)+b)     #
a =     13636     #
b =     19.873     #
A =     6165     #
B =     2.57*10^8     #
C =     43195     #
rate =     A*p*k/3600     #
z =     t     #
Tdif =     T-600     #
Ca =     p/(.082*T)     #
Cm =     Ua*(T-Tw)/(-DH*k)/1.987/T     #
Ua =     87/.0125     #
DH =     -307000     #
pm =     Qr/B/k     #
Deriv =    (Qg-Qr)/3600    #
Bst1 =  1+gam*(T-Tw)  #
Bst2 =  A*k/C  #
gam =  a/T/T  #
Cst =  T-Tw  #
P11 =  ((Bst1+Bst2)+((Bst1+Bst2)^2-4*gam*Cst)^.5)/gam/2  #
P12 =  ((Bst1+Bst2)-((Bst1+Bst2)^2-4*gam*Cst)^.5)/gam/2  #
P1 =  P11*C/B/k  #
P2 =  P12*C/B/k  #
t(0)=0
p(0)=.01825
T(0)=625
t(f)=1