Scribd
Carregar
9 visualizações8 páginas

Cinética

O documento analisa a conversão em quatro tipos de operação de um reator de leito fixo (PBR) com troca térmica: adiabático, temperatura constante do fluido refrigerante, escoamento cocorrente e contracorrente. A melhor conversão, de aproximadamente 55%, foi obtida com a operação de temperatura constante do fluido refrigerante.

Enviado por

Poliana Goes
Direitos autorais
© © All Rights Reserved
Levamos muito a sério os direitos de conteúdo. Se você suspeita que este conteúdo é seu, reivindique-o aqui.
Formatos disponíveis
Baixe no formato PDF, TXT ou leia on-line no Scribd
9 visualizações8 páginas

Cinética

O documento analisa a conversão em quatro tipos de operação de um reator de leito fixo (PBR) com troca térmica: adiabático, temperatura constante do fluido refrigerante, escoamento cocorrente e contracorrente. A melhor conversão, de aproximadamente 55%, foi obtida com a operação de temperatura constante do fluido refrigerante.

Enviado por

Poliana Goes
Direitos autorais
© © All Rights Reserved
Levamos muito a sério os direitos de conteúdo. Se você suspeita que este conteúdo é seu, reivindique-o aqui.
Formatos disponíveis
Baixe no formato PDF, TXT ou leia on-line no Scribd
Você está na página 1/ 8

Prova 3 - Exercício 1

–– Reação em um PBR com Troca Térmica ––

Gustavo Junqueira Valias Meira Filho | 11748


Poliana Larissa de Souza Goes | 119780

Enunciado
A reação elementar em fase líquida A ↔ 2B é conduzida em um reator de leito fixo. A puro entra no reator a
305 K a uma vazão de 14,67 kmol/h e uma concentração de 1,86 mol/dm 3. O PFR cercado por um trocador
de calor com fluido refrigerante a 315 K. Compare a conversão atingida para os quatro tipos de operação do
trocador de calor: adiabático, Ta constante, escoamento cocorrente e escoamento em contracorrente.

% Variáveis
KC = 3.03; % kmol-1 (T = 333,15 K)
k = 31.1; % (T = 360 K) h-1
Cp_c = 28; % kJ/kmol·K
Ua= 5000; % kJ/m³·h·K
DHRx = -34500; % kJ/kmol
T_R = 298; % K
CP_A = 141; % kJ/kmol·K
CP_B = 159; % kJ/kmol·K
E = 65700; % kJ/kmol
m_c = 500; % kmol/h
T0 = 305; % K
FA0 = 14.67; % kmol/k
CA0 = 1.86; % kmol/m³
Ta0 = 315; % K
R = 8.314; % kJ/kmol·K

Razão da Taxa:

A ↔ 2B

Balanço Molar:

Equação da Taxa e Estequiometria:

1
Conversão de Equilíbrio:

Balanço de Energia:

(a) Transferência em Cocorrente


Balanço para Fluido de Troca Térmica:

F0 = [Ta0; 0; T0];
[V,F] = ode89(@Cocorrente,[0 5],F0);
figure(1);
plot(V,F(:,1),'LineWidth',1.3,'Color',[0.3010 0.7450 0.9330]); hold on;
plot(V,F(:,3),'LineWidth',1.3,'Color',[0.6350 0.0780 0.1840]); hold off;
grid on; grid minor;
title('Temperaturas Cocorrentes');
legend('Ta [K]','T [K]');
xlabel('V [m3]'); ylabel('T [K]'); axis([0 5 290 380]);

figure(2);
plot(V,F(:,2),'LineWidth',1.3,'Color',[0.4660 0.6740 0.1880]);
grid on; grid minor;
title('Conversão Cocorrentes');
legend('X','Location','southeast');
xlabel('V [m3]'); ylabel('X');

2
(b) Contracorrente
Balanço para Fluido de Troca Térmica:

Para esses casos, devemos supor um valor inicial de em . Se no final da reação consiguirmos

atingir , então acertamos o valor. Chutaremos e verificaremos se em .

F0 = [333; 0; T0];
[V,F] = ode89(@Contracorrente,[0 5],F0);
figure(3);
plot(V,F(:,1),'LineWidth',1.3,'Color',[0.3010 0.7450 0.9330]); hold on;
plot(V,linspace(315,315,length(V)),'k--','LineWidth',1);
plot(V,F(:,3),'LineWidth',1.3,'Color',[0.6350 0.0780 0.1840]); hold off;
grid on; grid minor;
title('Temperaturas Contracorrentes');
legend('Ta [K]','Ta0 = 315 K','T [K]');
xlabel('V [m3]'); ylabel('T [K]'); axis([0 5 300 400]);

3
figure(4);
plot(V,F(:,2),'LineWidth',1.3,'Color',[0.4660 0.6740 0.1880]);
grid on; grid minor;
title('Conversão Contracorrente');
legend('X','Location','southeast');
xlabel('V [m3]'); ylabel('X');

(c) Temperatura Constante


Balanço para Fluido de Troca Térmica:

F0 = [Ta0; 0; T0];
[V,F] = ode89(@T_Constante,[0 5],F0);
figure(3);
plot(V,F(:,1),'LineWidth',1.3,'Color',[0.3010 0.7450 0.9330]); hold on;

4
plot(V,F(:,3),'LineWidth',1.3,'Color',[0.6350 0.0780 0.1840]); hold off;
grid on; grid minor;
title('Temperaturas Ta Constante');
legend('Ta [K]','T [K]');
xlabel('V [m3]'); ylabel('T [K]'); axis([0 5 300 380]);

figure(4);
plot(V,F(:,2),'LineWidth',1.3,'Color',[0.4660 0.6740 0.1880]);
grid on; grid minor;
title('Conversão Ta Constante');
legend('X','Location','southeast');
xlabel('V [m3]'); ylabel('X');

(d) Adiabático
Neste tipo de operação, temos .

F0 = [Ta0; 0; T0];

5
[V,F] = ode89(@Adiabatico,[0 5],F0);
figure(3);
plot(V,F(:,3),'LineWidth',1.3,'Color',[0.6350 0.0780 0.1840]);
grid on; grid minor;
title('Temperaturas Adiabático');
legend('T [K]','Location','southeast');
xlabel('V [m3]'); ylabel('T [K]'); axis([0 5 300 400]);

figure(4);
plot(V,F(:,2),'LineWidth',1.3,'Color',[0.4660 0.6740 0.1880]);
grid on; grid minor;
title('Conversão Adiabático');
legend('X','Location','southeast');
xlabel('V [m3]'); ylabel('X');

Tendo esses resultados em vista, a melhor conversão foi atingida no trocador com Ta constante, com
aproximadamente X = 0,55.

6
Funções Puxadas

% F(1) = Ta
% F(2) = X
% F(3) = T

function dFdV = Cocorrente(V,F)


% Variáveis do Enunciado
Cp_c = 28; % kJ/kmol·K
Ua= 5000; % kJ/m³·h·K
DHRx = -34500; % kJ/kmol
E = 65700; % kJ/kmol
m_c = 500; % kmol/h
FA0 = 14.67; % kmol/k
CA0 = 1.86; % kmol/m³
R = 8.314; % kJ/kmol·K
CP_A = 141; % kJ/kmol·K

% Equações Dependentes
k = 31.1.*exp(E./(R).*((1./360)-(1./F(3))));
KC = 3.03.*exp((DHRx./R).*((1/333.15)-(1./F(3))));
ra = -k.*CA0.*((1-F(2))) + k.*(((2.*CA0.*F(2)).^2)./KC);
% Equações de Taxa
dFdV = [(Ua.*(F(3)-F(1)))./(m_c.*Cp_c);
-ra./FA0;
((ra.*DHRx)-(Ua.*(F(3)-F(1))))./(FA0.*CP_A)];
end

% ------------------ %

function dFdV = Contracorrente(V,F)


% Variáveis do Enunciado
Cp_c = 28; % kJ/kmol·K
Ua= 5000; % kJ/m³·h·K
DHRx = -34500; % kJ/kmol
E = 65700; % kJ/kmol
m_c = 500; % kmol/h
FA0 = 14.67; % kmol/k
CA0 = 1.86; % kmol/m³
R = 8.314; % kJ/kmol·K
CP_A = 141; % kJ/kmol·K

% Equações Dependentes
k = 31.1.*exp(E./(R).*((1./360)-(1./F(3))));
KC = 3.03.*exp((DHRx./R).*((1/333.15)-(1./F(3))));
ra = -k.*CA0.*((1-F(2))) + k.*(((2.*CA0.*F(2)).^2)./KC);
% Equações de Taxa
dFdV = [-(Ua.*(F(3)-F(1)))./(m_c.*Cp_c);
-ra./FA0;

7
((ra.*DHRx)-(Ua.*(F(3)-F(1))))./(FA0.*CP_A)];
end

% ------------------ %

function dFdV = T_Constante(V,F)


% Variáveis do Enunciado
Cp_c = 28; % kJ/kmol·K
Ua= 5000; % kJ/m³·h·K
DHRx = -34500; % kJ/kmol
E = 65700; % kJ/kmol
m_c = 500; % kmol/h
FA0 = 14.67; % kmol/k
CA0 = 1.86; % kmol/m³
R = 8.314; % kJ/kmol·K
CP_A = 141; % kJ/kmol·K

% Equações Dependentes
k = 31.1.*exp(E./(R).*((1./360)-(1./F(3))));
KC = 3.03.*exp((DHRx./R).*((1/333.15)-(1./F(3))));
ra = -k.*CA0.*((1-F(2))) + k.*(((2.*CA0.*F(2)).^2)./KC);
% Equações de Taxa
dFdV = [(Ua.*(F(3)-F(1)))./(m_c.*Cp_c)*0;
-ra./FA0;
((ra.*DHRx)-(Ua.*(F(3)-F(1))))./(FA0.*CP_A)];
end

% ------------------ %

function dFdV = Adiabatico(V,F)


% Variáveis do Enunciado
Cp_c = 28; % kJ/kmol·K
Ua= 0; % kJ/m³·h·K
DHRx = -34500; % kJ/kmol
E = 65700; % kJ/kmol
m_c = 500; % kmol/h
FA0 = 14.67; % kmol/k
CA0 = 1.86; % kmol/m³
R = 8.314; % kJ/kmol·K
CP_A = 141; % kJ/kmol·K

% Equações Dependentes
k = 31.1.*exp(E./(R).*((1./360)-(1./F(3))));
KC = 3.03.*exp((DHRx./R).*((1/333.15)-(1./F(3))));
ra = -k.*CA0.*((1-F(2))) + k.*(((2.*CA0.*F(2)).^2)./KC);
% Equações de Taxa
dFdV = [(Ua.*(F(3)-F(1)))./(m_c.*Cp_c);
-ra./FA0;
((ra.*DHRx)-(Ua.*(F(3)-F(1))))./(FA0.*CP_A)];
end

Você também pode gostar

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy