Example 13-9 RTD and Complex Reactions: A. Reactor with Asymmetric RTD |POLVER05_0
d(ca)/d(t) = ra #
d(cb)/d(t) = rb #
d(cc)/d(t) = rc #
d(cabar)/d(t) = ca*E #
d(cbbar)/d(t) = cb*E #
d(ccbar)/d(t) = cc*E #
d(cd)/d(t) = rd #
d(ce)/d(t) = re #
d(cdbar)/d(t) = cd*E #
d(cebar)/d(t) = ce*E #
k1 = 1 #
k2 = 1 #
k3 = 1 #
E1 = -2.104*t^4+4.167*t^3-1.596*t^2+0.353*t-0.004 #
E2 = -2.104*t^4+17.037*t^3-50.247*t^2+62.964*t-27.402 #
rc = k1*ca*cb #
re = k3*cb*cd #
ra = -k1*ca*cb-k2*ca #
rb = -k1*ca*cb-k3*cb*cd #
E = if(t<=1.26)then(E1)else(E2) #
rd = k2*ca-k3*cb*cd #
t(0)=0
ca(0)=1
cb(0)=1
cc(0)=0
cabar(0)=0
cbbar(0)=0
ccbar(0)=0
cd(0)=0
ce(0)=0
cdbar(0)=0
cebar(0)=0
t(f)=2.52