Example 9-2 Safety in Chemical Plants with Exothermic Reactions  |POLVER05_0
d(T)/d(t) = if (t<45) then (0) else ((Qg-Qr)/NCp) #
d(X)/d(t) =  (-ra)*V/Nao  #
NCp =  2504  #
V =  3.265+1.854  #
Nao =  9.0448  #
UA =  35.83  #
dH =  -590000  #
Nbo =  33  #
k =  .00017*exp(11273/(1.987)*(1/461-1/T))  #
Qr =  UA*(T-298)  #
Theata =  Nbo/Nao  #
ra =  -k*Nao^2*(1-X)*(Theata-2*X)/V^2  #
rate =  -ra  #
Qg =  ra*V*(dH)  #
t(0)=55
T(0)=467.992
X(0)=0.0423866
t(f)=121