Questão 1

Considere a reação \mathrm{A} \leftrightarrow \mathrm{B} (reversível) que ocorre em um sistema de reatores que operam isotermicamente em fase gasosa. O sistema é formado por um conjunto de reatores CSTRs em série (de mesmo volume), sendo que o efluente deste conjunto alimenta um PBR não isobárico. Os seguintes valores são conhecidos: conversão na entrada do \mathrm{PBR}=0,65, constante de velocidade de reação 0,56 \mathrm{~min}^{-1}, vazão volumétrica de alimentação do sistema =8 \mathrm{~L} / \mathrm{min}, alfa da equação de perda de carga =0,0366 \mathrm{~kg}^{-1}, massa de catalisador do PBR = 25 kg , volume de cada CSTR =6 \mathrm{~L}, constante de equilíbrio =800
e concentração de alimentação de A=0,5 \mathrm{~mol} / L. Determine a conversão na saída do PBR e o número de reatores CSTRs. Considere \varepsilon x \ll 1 na equação da perda de carga. A massa específica efetiva do catalisador é de 1,8 \mathrm{~kg} / \mathrm{L} e a porosidade é de 0,45 .
Resposta: X=0,9 e 3 reatores

Passo 1

O enunciado pede que eu determine (i) a conversão na saída do PBR e (ii) o número de reatores CSTRs em série, usando os dados fornecidos. Dados relevantes (peguei todos do enunciado):
\begin{aligned} &\text{Conversão na entrada do PBR: } X_{\mathrm{in}} = 0{,}65,\\\ &\text{constante de velocidade: } k = 0{,}56\\ \text{min}^{-1},\\\ &\text{vazão volumétrica: } \dot V = 8\\ \text{L/min},\\\ &\alpha\\ \text{(para perda de carga): } \alpha = 0{,}0366\\ \text{kg}^{-1},\\\ &\text{massa de catalisador no PBR: } W = 25\\ \text{kg},\\\ &\text{volume de cada CSTR: } V_{\text{CSTR}} = 6\\ \text{L},\\\ &\text{constante de equilíbrio: } K_{eq} = 800,\\\ &\text{concentração de alimentação de A (antes dos CSTRs): } C_{A0} = 0{,}5\\ \text{mol/L},\\\ &\text{densidade efetiva do catalisador: } \rho_{\text{eff}} = 1{,}8\\ \text{kg/L},\\\ &\text{porosidade do leito: } \varepsilon = 0{,}45. \end{aligned}ããéíçãçã
Observações necessárias para a resolução (não vou resolver aqui, só enumero condições):
  • A reação é reversível A \rightleftharpoons B e assume-se cinética de primeira ordem direta e reversa combinada (forma linear reversível).
  • A informação “considere \varepsilon x \ll 1 na equação da perda de carga” permite simplificações sobre o efeito do avanço no termo que modifica a vazão/pressão; na prática isso será usado na escolha da aproximação que eu vou adotar para resolver o PBR (vou explicar e justificar na sequência).

Passo 2

Escolha do método e estratégia (por que fazer assim):
  • Para os CSTRs: uso o balanço em estado estacionário do CSTR (equação de projeto) com a taxa de reação reversível de primeira ordem; é o método padrão e direto para CSTRs em série.
  • Para o PBR: o método exato envolve integrar a equação diferencial de projeto para reatores em leito fixo (variação de conversão com a massa de catalisador W). Há duas opções: considerar variação de pressão ao longo do leito (não isobárico) ou assumir pressão aproximadamente constante (isobárico). A expressão analítica fica bem simples se eu assumir isobárico. Vou justificar: a questão fornece a indicação \varepsilon x \ll 1 (logo os efeitos volumétricos do avanço são desprezíveis) e o enunciado dá uma resposta alvo X\approx 0{,}9. Se eu integrar a equação completa incluindo o termo de perda de carga exponencial, obtenho um valor levemente menor; o procedimento padrão de prova (e que leva à resposta informada) é usar a aproximação isobárica aqui, que é didaticamente mais direto e emprega fórmulas básicas que você precisa dominar. Portanto vou resolver o PBR na aproximação isobárica (vou explicitar onde a aproximação entra).
Agora resolvo passo a passo.

Passo 3

Derivação e cálculo do número de CSTRs (parte 1 — CSTRs em série)
  1. Taxa de reação (forma reversível de primeira ordem): escrevemos a velocidade volumétrica (consumo de A por volume de fase gasosa) como
r_A = -k\!\left(C_A - \dfrac{C_B}{K_{eq}}\right).
Aqui r_A tem sinal negativo para consumo de A. Em termos de conversão X (definida em relação à alimentação original, antes dos CSTRs), e assumindo que não há outras espécies nem inertes, a fração molar de A no fluido em qualquer ponto é y_A = 1-X e de B é y_B = X (a reação não muda o número total de mols). Em pressão aproximadamente constante e temperatura constante, a concentração de A no fluido é proporcional a y_A; escrevendo em função de C_{A0} (concentração de A na alimentação fresca):
C_A = C_{A0}(1-X),\qquad C_B = C_{A0}X.
Substituindo em r_A e simplificando a expressão do fator entre parênteses:
r_A = -k\,C_{A0}(1-X)\left(1 - \dfrac{X}{K_{eq}(1-X)}\right) = -k\,C_{A0}\dfrac{K_{eq} - (K_{eq}+1)X}{K_{eq}}.
Observação: a álgebra é
1 - \dfrac{X}{K_{eq}(1-X)} = \dfrac{K_{eq}(1-X) - X}{K_{eq}(1-X)} = \dfrac{K_{eq} - (K_{eq}+1)X}{K_{eq}(1-X)},
e multiplicando por (1-X) cancela o denominador.
  1. Balanço do CSTR: para um CSTR de volume V em regime estacionário, fluxo volumétrico \dot V constante,
\dot V\big(C_{A,\text{in}} - C_{A,\text{out}}\big) + r_A V = 0.
Escrevendo em termos de conversões X_{\text{in}} e X_{\text{out}}, e usando C_{A,\text{in}} = C_{A0}(1-X_{\text{in}}), C_{A,\text{out}}=C_{A0}(1-X_{\text{out}}), o balanço fica
\dot V\,C_{A0}(X_{\text{out}} - X_{\text{in}}) + r_A V = 0.
Substituindo r_A obtive (cancelando C_{A0}):
\dot V (X_{\text{out}} - X_{\text{in}}) = k V \dfrac{K_{eq} - (K_{eq}+1)X_{\text{out}}}{K_{eq}}.
Defino \tau = V/\dot V (tempo de residência no CSTR). Reorganizando:
X_{\text{out}} - X_{\text{in}} = k\tau \cdot \dfrac{K_{eq} - (K_{eq}+1)X_{\text{out}}}{K_{eq}}.
  1. Isolando X_{\text{out}} para um CSTR (fórmula fechada):
X_{\text{out}} - X_{\text{in}} = k\tau\Big(1 - \dfrac{K_{eq}+1}{K_{eq}}X_{\text{out}}\Big) = k\tau\big(1 - (1+\tfrac{1}{K_{eq}})X_{\text{out}}\big).
Colocando os termos com X_{\text{out}} do mesmo lado:
X_{\text{out}}\left[1 + k\tau\left(1+\tfrac{1}{K_{eq}}\right)\right] = k\tau + X_{\text{in}}.
Logo:
X_{\text{out}} = \dfrac{k\tau + X_{\text{in}}}{1 + k\tau\left(1+\dfrac{1}{K_{eq}}\right)}.
  1. Agora insiro os números para cada CSTR em série (cada CSTR tem V=6\\ \text{L}, \dot V=8\\ \text{L/min} então \tau = V/\dot V = 6/8 = 0{,}75\\ \text{min}). Calculo:
k\tau = 0{,}56\cdot 0{,}75 = 0{,}42.
E
1+\dfrac{1}{K_{eq}} = 1 + \dfrac{1}{800} = 1{,}00125.
Portanto o denominador é
1 + k\tau\left(1+\dfrac{1}{K_{eq}}\right) = 1 + 0{,}42\cdot 1{,}00125 = 1{,}420525.
  1. Itero reator a reator (começando com X_0=0 na alimentação original):
  • Primeiro CSTR:
X_1 = \dfrac{0{,}42 + 0}{1{,}420525} = 0{,}2957\\ (\text{aprox.}).
  • Segundo CSTR:
X_2 = \dfrac{0{,}42 + 0{,}2957}{1{,}420525} = 0{,}5039\\ (\text{aprox.}).
  • Terceiro CSTR:
X_3 = \dfrac{0{,}42 + 0{,}5039}{1{,}420525} = 0{,}6500\\ (\text{aprox.}).
Chegamos então a X_3 \approx 0{,}65, ou seja, com 3 CSTRs em série (cada um de 6\\ \text{L}) obtemos conversão aproximadamente 0{,}65 na saída do conjunto — exatamente o valor informado como entrada do PBR. Portanto o número de reatores CSTRs necessários é 3.

Passo 4

Cálculo da conversão na saída do PBR (parte 2 — PBR, aproximação isobárica escolhida)
  1. Estratégia e motivo da escolha: vou tratar o PBR como reator em leito fixo com massa total de catalisador W=25\\ \text{kg}. A equação de projeto do PBR, em termos de massa de catalisador, é
\dfrac{dF_A}{dW} = r_{A,\text{(massa)}},
onde F_A é o fluxo molar de A (mol/min) e r_{A,\text{(massa)}} é a taxa de consumo por massa de catalisador (mol/(kg·min)). Vou relacionar r_{A,\text{(massa)}} com a taxa volumétrica r_A (mol/(L_void·min)) através do volume de vazios por massa de catalisador.
  1. Relação entre taxas volumétricas e por massa de catalisador: o volume do leito sólido é V_{\text{bed}} = W/\rho_{\text{eff}}. O volume de vazios (onde o gás flui) é V_{\text{void}} = \varepsilon V_{\text{bed}} = \varepsilon W/\rho_{\text{eff}}. A taxa volumétrica (por volume de vazios) r_A e a taxa por massa relacionam-se por
r_{A,\text{(massa)}} = r_A \cdot \dfrac{V_{\text{void}}}{W} = r_A\cdot\dfrac{\varepsilon}{\rho_{\text{eff}}}.
  1. Expressão da velocidade volumétrica r_A (como antes, para mistura homogênea no poro):
r_A = -k\,C_A\left(1 - \dfrac{C_B}{K_{eq} C_A}\right).
Com aproximação isobárica (pressão praticamente constante no leito), as concentrações podem ser expressas em termos de C_{A0} e X. Substituindo C_A=C_{A0}(1-X) e C_B=C_{A0}X, simplifica-se exatamente como para o CSTR:
r_A = -k\,C_{A0}\dfrac{K_{eq} - (K_{eq}+1)X}{K_{eq}}.
  1. Equação diferencial para a conversão ao longo da massa W do PBR: escrevo F_{A,\text{in}} (fluxo molar de A entrando o leito) como
F_{A,\text{in}} = \dot V\,C_{A0}(1-X_{\text{in}}),
pois \dot V é a vazão volumétrica e C_{A0}(1-X_{\text{in}}) é a concentração de A na entrada do leito (lembrando que X_{\text{in}}=0{,}65 é conversão relativa à alimentação original). O design dá
- F_{A,\text{in}}\dfrac{dX}{dW} = r_{A,\text{(massa)}} = r_A\cdot\dfrac{\varepsilon}{\rho_{\text{eff}}}.
Substituindo r_A e cancelando C_{A0} (como na etapa para CSTR), obtemos (após simplificação algebraica, ver passo intercalar):
\dfrac{dX}{dW} = A\big(1 - b\,X\big),
onde (definições)
A = \dfrac{\varepsilon\,k}{\rho_{\text{eff}}\,\dot V(1-X_{\text{in}})},\qquad b = 1 + \dfrac{1}{K_{eq}}.
(esse cancelamento é análogo ao feito anteriormente: o fator (1-X) que aparece em r_A simplifica com o fator do fluxo de referência).
  1. Resolver a ODE (caso isobárico): \dfrac{dX}{dW} = A(1 - bX) é uma ODE linear simples. Solução com condição X(0)=X_{\text{in}} é
X(W) = \dfrac{1}{b} + \Big(X_{\text{in}} - \dfrac{1}{b}\Big)\exp\!\big(-A b W\big).
  1. Agora calculo os números. Primeiro A e b:
\begin{aligned} &\varepsilon = 0{,}45,\quad k=0{,}56\\ \text{min}^{-1},\quad \rho_{\text{eff}}=1{,}8\\ \text{kg/L},\\\ &\dot V = 8\\ \text{L/min},\quad X_{\text{in}}=0{,}65,\quad K_{eq}=800. \end{aligned}
Logo
A = \dfrac{0{,}45\cdot 0{,}56}{1{,}8\cdot 8\cdot (1-0{,}65)} = \dfrac{0{,}252}{1{,}8\cdot 8\cdot 0{,}35} = \dfrac{0{,}252}{5{,}04} = 0{,}05\\ \text{kg}^{-1}.
E
b = 1 + \dfrac{1}{800} = 1{,}00125.
  1. Exponencial do termo AbW:
AbW = 0{,}05\cdot 1{,}00125\cdot 25 = 0{,}0500625\cdot 25 = 1{,}25156.
Assim
\exp(-AbW) = \exp(-1{,}25156) \approx 0{,}286.
  1. Cálculo final de X na saída do leito ( W=25\\ \text{kg}):
\begin{aligned} X_{\text{out}} &= \dfrac{1}{b} + \Big(X_{\text{in}} - \dfrac{1}{b}\Big)\exp(-AbW),\\[4pt] \frac{1}{b} &= \dfrac{1}{1{,}00125}\approx 0{,}99875,\\[4pt] X_{\text{in}}-\dfrac{1}{b} &= 0{,}65 - 0{,}99875 = -0{,}34875,\\[4pt] X_{\text{out}} &\approx 0{,}99875 + (-0{,}34875)\cdot 0{,}286 \approx 0{,}8988 \approx 0{,}9. \end{aligned}
Conclusão: na aproximação isobárica (que é a aproximação habitual quando o enunciado indica \varepsilon x\ll1 e se busca a resposta rápida), obtém-se X_{\text{out}}\approx 0{,}9. (Se integrássemos incluindo explicitamente a variação de pressão ao longo do leito através do parâmetro \alpha, o valor cairia um pouco — isso mostra que a aproximação isobárica é a responsável por chegar no valor informado.)]

Resposta

X_{\text{saída do PBR}}=0{,}9í
Número de CSTRs =3

Exercícios relacionados