Menghitung Temperatur RATB

% Program untuk Menghitung Temperatur RATB yang dilengkapi jaket pendingin

………….program start…………

clc, clear
T=535:650;%Rentang Temperatur yang dicoba dengan interval 1 R
Ta1=545;%Asumsi suhu pendingin masuk (derajat Reamur)
Ta2=565;%Asumsi suhu pendingin keluar (derajat Reamur)
XEB=(403.3.*(T-535)+92.9.*((Ta1-Ta2)./log((T-Ta1)./(T-Ta2))))./(36400+7.*(T-528));
XMB=2.084e12.*exp(-16306./T)./(1+2.084e12.*exp(-16306./T));
plot(T,XEB,’r’)
hold on
plot(T,XMB,’b’)
hold off
legend(‘XEB’,’XMB’)
grid on
xlabel(‘Temperatur(R)’)
ylabel(‘Konversi’)
f=inline(‘(403.3.*(T-535)+92.9.*((545-565)./log((T-545)./(T-565))))./(36400+7.*(T-528))-2.084e12.*exp(-16306./T)./(1+2.084e12.*exp(-16306./T))’);
t=fzero(f,650);
disp([‘Temperatur reaktor=’,num2str(t), ‘  derajat Reamur’])

……………program end……….

Command Window

Temperatur reaktor=643.1438 derajat Reamur