Exercícios Resolvidos de Difusão com Reação Química

Ver Teoria

Enunciado

Uma corrente gasosa contendo um reagente “A” entra em contato com um catalisador de geometria esférica de raio R. Esta partícula está dentro de um reator catalítico. Nas imediações da partícula catalítica, a concentração do reagente “A” é C A S (moles/cm3 ). A espécie “A” difunde através dos poros existentes no catalisador e converte no produto “B” através de uma reação irreversível e de primeira ordem no sítio ativo do mesmo. O produto “B” difunde no sentido contrário do reagente “A”. Determine o perfil de concentração do reagente “A” em função do raio da partícula considerando que o processo de transferência de massa ocorra em regime permanente e a temperatura e pressão constante.

Sendo:

a = área da superfície catalítica por unidade de volume do catalisador c m 2 / c m 3

k S = constante de velocidade c m / s

C A = concentração do reagente m o l / c m 3

Passo 1

Mais uma vez vamos trabalhar com a equação em coordenadas esféricas:

∂ C A ∂ t + 1 r 2 ⋅ ∂ r 2 N A ,   r ∂ r + 1 r ⋅ s e n   θ ⋅ ∂ s e n   θ ⋅ N A ,   θ ∂ θ + 1 r ⋅ s e n   θ ⋅ ∂ N A ,   φ ∂ φ = R A "

Fazendo nossas considerações…

Como estamos trabalhando em regime permanente, então: ∂ C A ∂ t = 0 ;

O fluxo é radial e unidirecional, então: ∇ ⋅ N A = ∂ r 2 N A ,   r ∂ r ;

E considerando que a reação acontece dentro do poro do catalisador, temos que: R A " =   - k S ⋅ a ⋅ C A

Sendo assim nossa equação, fica:

1 r 2 ⋅ ∂ r 2 N A ,   r ∂ r =   - k S ⋅ a ⋅ C A

O fluxo molar no interior da parte porosa, levando em conta o coeficiente de difusão efetiva, D e f , será:

N A ,   r =   - C D e f d y A   d r + y A ( N A ,   r + N B ,   r )

Sendo:

D e f = D A B ⋅ ε p τ

Onde ε p = porosidade do catalisador e τ = tortuosidade do catalisador.

E sendo N B ,   r =   - N A ,   r , nossa equação será:

N A ,   r =   - C D e f d y A   d r

N A ,   r =   - D e f d C A   d r

Agora substituindo isso na equação ali de cima, vamos ter:

1 r 2 ⋅ ∂ ∂ r r 2 ⋅ D e f d C A   d r =   k S ⋅ a ⋅ C A

Como o enunciado disse pra gente que estamos trabalhando com pressão e temperatura constantes, então sabemos que D e f também é constante, ficamos então com:

∂ ∂ r r 2 ⋅ d C A   d r = r 2 ⋅   k S ⋅ a ⋅ C A D e f

Chamando λ 2 = k S ⋅ a D e f , teremos:

∂ ∂ r r 2 ⋅ d C A   d r = r 2 ⋅ λ 2 ⋅ C A

Precisamos encontrar um jeito de melhor essa equação, porque esse C A vai nos atrapalhar nas integrações, ok? Vamos fazer o seguinte então:

r ⋅ C A = Ψ C A = Ψ r

Teremos então que:

d C A   d r = d Ψ d r r   -   Ψ d r d r r 2

d C A   d r = d Ψ d r r   -   Ψ r 2

Substituindo isso na nossa equação, ficaremos com:

d d r r 2 ⋅ d Ψ d r r   -   Ψ r 2 = r 2 ⋅ λ 2 ⋅ Ψ r

d d r d Ψ d r r   -   Ψ = r ⋅ λ 2 ⋅ Ψ

d 2 Ψ d r 2 r + d Ψ d r - d Ψ d r = r ⋅ λ 2 ⋅ Ψ

d 2 Ψ d r 2 = λ 2 ⋅ Ψ

Integrando isso aí, duas vezes, vamos ficar com:

Ψ = C 1 c o s h λ r + C 2 s e n h λ r

Agora voltando com o C A , vamos ter:

r ⋅ C A = C 1 c o s h λ r + C 2 s e n h λ r

Passo 2

Agora é a hora de estabelecer nossas condições de contorno, pra que a gente descubra quem são C 1 e C 2

A nossa primeira condição de contorno, se refere a entrada do poro, onde:

r = R   →   C A = C A S

Temos então:

R ⋅ C A S = C 1 c o s h λ R + C 2 s e n h λ R

E usando a segunda condição de contorno, que se refere ao centro da partícula catalítica, temos:

r = 0   →   C A = C A *

Isso quer dizer que temos:

0 ⋅ C A * = C 1 c o s h λ 0 + C 2 s e n h λ 0

0 = C 1 ⋅ 1 + C 2 ⋅ 0

O C 1 será então:

C 1 = 0

Agora vamos substituir C 1 lá na primeira condição de contorno, beleza?

R ⋅ C A S = 0 ⋅ c o s h λ R + C 2 s e n h λ R

R ⋅ C A S = C 2 s e n h λ R

Logo, o C 2 vai ser:

C 2 = R ⋅ C A S s e n h λ R

Agora substituindo C 1 e C 2 na equação que a gente tinha antes, ficaremos com:

r ⋅ C A = 0 ⋅ c o s h λ r + R ⋅ C A S s e n h λ R ⋅ s e n h λ r

Por fim, rearranjando essa galera, o perfil de concentração do reagente “A” em função do raio da partícula será:

C A C A S = R r ⋅ s e n h λ r s e n h λ R

Uhull! Conseguimos!

Resposta

O perfil de concentração do reagente “A” em função do raio da partícula:

C A C A S = R r ⋅ s e n h λ r s e n h λ R

Exercícios de Livros Relacionados

Para aumentar a superfície efetiva, e assim a taxa de reação química, as superfícies dos catalisadores frequentemente apresentam a forma de sólidos porosos. Um desses sólidos pode ser visualizado como
Um novo processo foi proposto para criar um tubo de compósito de paládio para ser usado como uma membrana para separação de hidrogênio na produção de hidrogênio com alta pureza. Para fabricar o tubo d
1. Refaça o exemplo (6.2) considerando reação de primeira ordem na superfície do grafite.
Emissões de óxido nítrico (NO) pela descarga de automóveis podem ser reduzidas com o uso de um conversor catalítico. A seguinte reação ocorre na superfície do catalisador: N O+C O→N 2+C O 2 A concentr
Carvão pulverizado, com partículas que podem ser aproximadas por esferas com raio r o = 1m m , é queimado em uma atmosfera de oxigênio puro a 1450 K e 1 atm. Oxigênio é transferido para a superfície d