PFR Runaway Finding pm, p1and p2 versus T  |POLVER05_0
d(T)/d(t) =   2   #
Ca =   ((alpha*(T-Ta)*gam+alpha+beta)+((alpha*(Ta-T)*gam-alpha-beta)^2+4*gam*alpha^2*(Ta-T))^0.5)/2/gam/k1   #
Ca1 =   ((alpha*(T-Ta)*gam+alpha+beta)-((alpha*(Ta-T)*gam-alpha-beta)^2+4*gam*alpha^2*(Ta-T))^0.5)/2/gam/k1   #
alpha =      Ua/(-DH)      #
Ua =   82.7*2/.0125   #
DH =      -307000      #
beta =   k1*CaoCp/(-DH)   #
k =      exp((-a/T)+b)      #
a =      13636      #
b =      19.873      #
CaoCp =      .323      #
gam =      a/T/T      #
Ta =    625    #
p =   Ca*R*T   #
R =      .082      #
p1 =   Ca1*R*T   #
p1o =   Ca*29.48/1.293*1   #
p2o =   Ca1*29.48/1.293*1   #
P1 =     P11*C/(B*k)     #
P2 =     P12*C/(B*k)     #
C =   Ua/.323   #
A =   29.48*1/1.293*1300*0.208   #
B =   (-DH)*1300*0.208/.323   #
Cst =     T-Ta     #
Bst1 =     1+gam*(T-Ta)     #
Bst2 =     A*k/C     #
P11 =   ((Bst1+Bst2)+((Bst1+Bst2)^2-4*gam*Cst)^.5)/gam/2   #
P12 =   ((Bst1+Bst2)-((Bst1+Bst2)^2-4*gam*Cst)^.5)/gam/2   #
pm =    C*(T-Ta)/B/k    #
k1 =   k*1300*0.208*22.8   #
t(0)=0
T(0)=625
t(f)=50