Steady state operation of an exothermic CSTR   |POLVER05_3
f(ca) = F*(2.88-ca)/V-k*ca*ca #
f(t1) = F*(66-t1)/V-dhr*k*ca*ca/(rho*cp)-u*a*(t1-tc)/(V*rho*cp) #
f(tc) = u*a*(t1-tc)/(1.82*1000*4184)-fc*(tc-27)/1.82 #
f(b) = ((t1-80)/20-b)/20 #
f(y) = (m-y)/taui #
k=0.0744*exp(-1.182e7/(8314.39*(t1+273.16)))  #
fc = 0.02*50^(-m)  #
m = y+kc*(10/20-b) #
F=0.0075  #
V=7.08  #
dhr=-9.86e7  #
rho=19.2  #
cp=1.815e5  #
u=3550  #
a=5.4  #
taui=600  #
kc=1  #
ca(0)=1
t1(0)=100
tc(0)=50
b(0)=0.4
y(0)=0.25