No Title  |POLVER05_0 |1
#Example 13-5 Multiple Reactions in a Semibatch Reactor  |POLVER05_0 |1
d(Ca)/d(t) = ra+(Cao-Ca)*vo/V #
d(Cb)/d(t) = rb-Cb*vo/V #
d(Cc)/d(t) = rc-Cc*vo/V #
d(T)/d(t) = (Qg-Qr)/((Ca*30+Cb*60+Cc*20)*V+100*35) #
UA = 35000 #
dHrx1 = -6500 #
dHrx2 = 8000 #
Qg = (dHrx1*(-k1a*Ca)+dHrx2*(-k2b*Cb))*V #
Qr = UA*(T-298) + Cao*vo*30*(T-305) #
Cao = 4 #
vo = 240 #
k1a = 1.25*exp((9500/1.987)*(1/320-1/T)) #
k2b = 0.08*exp((7000/1.987)*(1/300-1/T)) #
ra = -k1a*Ca #
V = 100+vo*t #
rc = 3*k2b*Cb #
rb = k1a*Ca/2-k2b*Cb #
t(0)=0
Ca(0)=1
Cb(0)=0
Cc(0)=0
T(0)=290
t(f)=1.5