%Syed Imran Farogh
%Student number 336020
%For multi steady state values
%constant values
q = 0.01666667; %volumetric Flowrate (m^3/sec)
v = 0.1; %volume of CSTR (m^3)
rho = 1000; %density (Kg/m^3)
cp = 239; %Heat capacity (J/Kg.K)
dH = -5e4; %heat of reaction (J/mol)
EoverR = 8750; % E=activation energy (J/mol),R=Universal gas constant(J/mol.K)
ko = 1.2e8; %(1/sec)kinetic pre-exponential factor
UA = 833.3333; %(J/sec.K)or (W/sec.K), U=heat transfer coefficient(W/m^2-K), A=Area(m2)
caf = 1000; %fluid concentration (mol/m^3)
tf = 350; %temperature(K)
tcf = 300; %coolant temperature(K)
T = 300:10:1000;
n = length(T);
for i=1:n
ca(i) = caf/(1+(ko*v*exp(-EoverR/T(i)))/q);
rate(i) = ko*exp(-EoverR/T(i))*ca(i);
d(i) = q/v*(tf-T(i))+((-dH/(rho*cp))*rate(i))+((UA/(v*rho*cp))*(tcf-T(i)));
end
plot(T,d);
xlabel('temperature (T) K','FontSize',12)
ylabel ('d (objective fuction)','FontSize', 12)
评论3