No Title  |POLVER05_0 |1
  #Example 12-5 Parallel Reaction in a PFR with Heat Effects  |POLVER05_0 |1
d(Fa)/d(V) = r1a+r2a #
d(Fb)/d(V) = -r1a #
d(Fc)/d(V) = -r2a/2 #
d(T)/d(V) = (Qg-Qr)/(90*Fa+90*Fb+180*Fc) #
Qr = 4000*(T-373)#
Qg = (-r1a)*(-deltaH1)+(-r2a)*(-deltaH2) #
k1a = 10*exp(4000*(1/300-1/T)) #
k2a = 0.09*exp(9000*(1/300-1/T)) #
Cto = 0.1 #
deltaH1 = -20000 #
deltaH2 = -60000 #
Ft = Fa+Fb+Fc #
To = 423 #
Ca = Cto*(Fa/Ft)*(To/T) #
Cb = Cto*(Fb/Ft)*(To/T) #
Cc = Cto*(Fc/Ft)*(To/T) #
r1a = -k1a*Ca #
r2a = -k2a*Ca^2 #
V(0)=0
Fa(0)=100
Fb(0)=0
Fc(0)=0
T(0)=423
V(f)=1