Problem D9 - Reversible Exothermic Gas Phase Reaction in PBR    |POLVER05_0
d(x)/d(W) =   -rA/FA0   #
d(y)/d(W) =   -0.015*(1-.5*x)*(T/450)/(2*y)   #
d(T)/d(W) =   (.8*(Ta-T)+rA*delH)/(CPA*FA0)   #
rA =   -k*(CA^2-CC/Kc)   #
k =   .5*exp(5032*(1/450-1/T))   #
Kc =   25000*exp(delH/8.314*(1/450-1/T))   #
CA =   .271*(1-x)*(450/T)/(1-.5*x)*y   #
CC =   .271*.5*x*(450/T)/(1-.5*x)*y   #
FA0 =   5   #
Ta =   500   #
delH =   -40000   #
CPA =   40   #
W(0)=0
x(0)=0
y(0)=1
T(0)=450
W(f)=20