O mesmo exemplo que o anterior, com outros valores de k's e com cáculo de selectividades e conversão
% Reactor descontínuo isotérmico com reacção em série paralelo. % A <--> B <-->C clc; clear; clf; % as concentrações iniciais são: Ao=1;Bo=0;Co=0; % mole/L % as constantes de velocidade são: global k
k=[0.1 0.05;0.1 0.05]; % Como são todas de 1ª ordem, têm unidade de 1/s % A->B A<-B % B->C B<-C % Construção e resolução dos balanços de massa feitos na subrotina % BMser.m % Serão resolvidos pelo metodo ode45 do Matlab entre t=0 e 40s com % condições iniciais C(0)=Co [t,C]=ode45(@BMser,[0 40],[Ao Bo Co]); % C é uma matriz com 3 colunas (uma por espécie) e N linhas, N determinado % pelo dt determinado pelo método ode45. subplot(2,1,1); plot(t,C); % Melhoramentos ao gráfico title('Reaccao A<->B<->C')
legend('A', 'B', 'C');
ylabel('C (M)');% porque as concentações estão em moles/L grid on xlabel('t (s.)');% porque as constantes k estão em s % Definir outro gráfico subplot(2,1,2); % S=A que reagiu S=C(:,2)+C(:,3); % X a conversão X=S/Ao*100; % Calculo das Selectividades a B e a C SB=C(:,2)'./S';SB=SB'*100; SC=C(:,3)'./S';SC=SC'*100; plot(t,X,'b',t,SB,'g',t,SC,'r');
% Melhoramentos ao gráfico legend('X','SB','SC');
ylabel('%')
grid on xlabel('t (s.)')
function dydt = BMser(t,y)
global k
dAdt = -k(1,1)*y(1)+k(1,2)*y(2);
dBdt=k(1,1)*y(1)-k(1,2)*y(2)-k(2,1)*y(2)+k(2,2)*y(3);
dCdt=k(2,1)*y(2)-k(2,2)*y(3);
dydt=[dAdt; dBdt; dCdt];
end
