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