% Reactor descontínuo isotérmico com reacção em paralelo. % A --> B % \-->C clc; clear; clf; % as concentrações iniciais são: Ao=1;Bo=0;Co=0; % as constantes de velocidade são: global k
k=[0.01;0.02];
% A->B % A->C % Construção e resolução dos balanços de massa feitos na subrotina % BMpar.m % Serão resolvidos pelo metodo ode45 do Matlab entre t=0 e 1000 com % condições iniciais C(0)=Co [t,C]=ode45(@BMpar,[0 200],[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 A->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.)')
