CSTR Code
CSTR Code
function f = cstr(x,y)
FV=1;
Tj=298;
Ti=298;
UAV=150;
H=5960;
rhcp=500;
k0=9705*3600;
ca0=10;
E=11840;
R=1.986;
DT =@(T) (FV*(Ti-T)
-(UAV/(rhcp))*(T-Tj)+(k0*exp(-E/(R*T))/(rhcp))*((H*FV*ca0)/(FV+(k0*exp(-E/
(R*T))))));
roots = fzero(DT,[298]);
ca=((FV*ca0)/(FV+k0*exp(-E/(R*roots))));
disp('roots')
disp(roots)
disp('ca')
disp(ca)
>> cstr
roots
311.0045
ca
8.5817