%************************************************************************** % Programa para mostrar a evolução da conversão de estado estacionário com a temperatura %de um reactor CSTR com uma reacção elementar equilibrada (A <-> B) endo- ou exotérmica. % Notar que o sistema é de 1ª ordem e por isso independente de CAo. % Notar que como o nº de moles se mantem, o sistema é válido para fase líquida e para gases % ideais. % (c)FGF 2008 %************************************************************************** clear; clc; clf; %limpar tudo ti=300; % temperatura ,mais baixa
Xe=zeros(1,1000-ti); xr=zeros(1,1000-ti); %préalocamento das matrizes para poupar tempo de cpu R=1.9872; % V=10;Q=1; %Volume (l), caudal (l t-1) dHr=-10000; %entalpia da reacção (<0 exotérmica) cal/mol/K Ead=20000; %Energia de activação da reacção directa (cal/mol/K) Eai=-dHr+20000; %Energia de activação da reacção inversa, levando em conta a entalpia da reacção for T=ti+1:1:1000;% para todas as temperaturas calcule kD=0.1*exp(-Ead/R*(1/T-1/500)); % constante de velocidade directa ki=0.1*exp(-Eai/R*(1/T-1/500)); %% constante de velocidade inversa %********************************************************************* % As constantes de velocidade estão definidas de modo a que 500K ambas % sejam 0.1 (t-1) %********************************************************************** K=kD/ki; % constante de equilíbrio Xe(T-ti)=K/(K+1);%conversão de equilibrio xr(T-ti)=1/(1/K+Q/(V*kD)+1);%conversão do CSTR (balanço de massa resolvido em ordem a x). end %fazer gráfico T=ti+1:1:1000; %vector de temperaturas plot(T,xr,T,Xe); legend('xr','Xe');
if dHr>0, title (' Conversão de uma reacção Endotérmica'),elseif title (' Conversão de uma reacção Exotérmica'),end xlabel('T (K)');
ylabel('conversão');
