A coluna de destilação que fará a separação de uma mistura de benzeno - tolueno com 60 \% em massa de tolueno operará a 101,3 \mathrm{kPa}. O produto de topo deve conter 97 \% em massa de benzeno, o produto de fundo 98 \% em massa de tolueno e a razão de refluxo externa no topo usada será de 3,5 . A carga e o refluxo no topo são líquidos saturados e devem entrar no seu prato ótimo com uma vazão de 30 \mathrm{~kg} / \mathrm{h}. Determine:
a) As vazões dos produtos de topo e de fundo;
b) O número de estágios de equilíbrio necessários;
c) O prato ótimo de carga;
d) A razão de refluxo mínima e
e) A razão de refluxo total.
Table 1: Dados de Equilíbrio do sistema Benzeno - Tolueno a 101,3 \mathrm{kPa}
\captionsetup{labelformat=empty}
x 0,21 0,37 0,51 0,64 0,72 0,79 0,86 0,91 0,96 0,98
y 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 0,95

Passo 1

O enunciado pede que a coluna de destilação, operando a 101{,}3\\ \mathrm{kPa}, separe uma mistura benzeno–tolueno com as seguintes condições e dados:
  • Composição da carga (carga = feed) em massa: 60\% de tolueno e 40\% de benzeno. Vou chamar a vazão de carga de F (kg/h).
  • O produto de topo (distilado, D) deve conter 97\% em massa de benzeno → isto implica 3\% em massa de tolueno no topo. Vou denotar w_D(\text{T})=0{,}03.
  • O produto de fundo (resíduo, B) deve conter 98\% em massa de tolueno → w_B(\text{T})=0{,}98.
  • A razão de refluxo externa dada é R_\text{ext}=3{,}5. Por definição usual, R = L/D (razão de refluxo = fluxo de refluxo dividido pelo fluxo de destilado).
  • A carga (feed) é um líquido saturado que entra no seu prato ótimo. O enunciado diz que a vazão da carga é 30\\ \mathrm{kg/h}. Vou tomar isto como F=30\\ \mathrm{kg/h}.
  • A carga e o refluxo no topo são líquidos saturados (informação usada para a classificação do feed: q=1, ou seja, linha- q vertical).
  • Há uma tabela de equilíbrio binário fornecida: são pares (x,y) — interpreto estes como frações molares de tolueno no líquido ( x) e no vapor ( y) para o equilíbrio líquido–vapor a 101{,}3\\ \mathrm{kPa}. (Vou usar essa tabela para o McCabe–Thiele por interpolação linear nos trechos dados.)
  • O professor pede para determinar:
    a) as vazões dos produtos de topo D e de fundo B;
    b) o número de estágios de equilíbrio necessários;
    c) o prato ótimo de carga (posição do prato onde a carga deve entrar);
    d) a razão de refluxo mínima R_\text{min};
    e) a razão de refluxo total (interpreto como a razão de refluxo R_\text{total} usada na operação — ou, se o enunciado distinguir, explico).
Observações importantes (para a resolução):
  • Vou transformar as frações em massa fornecidas para frações molares quando necessário (a tabela de equilíbrio é em frações molares). Para isso uso massas molares aproximadas: benzeno M_\text{Be}=78{,}11\\ \mathrm{g/mol} e tolueno M_\text{T}=92{,}14\\ \mathrm{g/mol}.
  • A linha- q é vertical (porque o feed é líquido saturado: q=1). Isso afeta onde o tramo de rectificação se conecta com o tramo de estrippagem.
  • Vou primeiro fazer balanços de massa para obter D e B. Depois passo para McCabe–Thiele (conversão para frações molares, uso da tabela por interpolação, construção das retas operacionais, contagem de degraus). Para R_\text{min} uso o método do tangente do rectifying line ao equilíbrio (explico passo a passo).

Passo 2

(a) Cálculo das vazões D e B por balanços de massa (parte componente e balanço total). Vou trabalhar em massa (kg/h), porque os dados F e as puridades estão em massa.
Dados:
F = 30\\ \mathrm{kg/h},\qquad w_F(\text{T})=0{,}60,\qquad w_D(\text{T})=0{,}03,\qquad w_B(\text{T})=0{,}98.
Balanço total:
F = D + B.
Balanço do componente tolueno (massa):
w_F(\text{T})\,F \;=\; w_D(\text{T})\,D \;+\; w_B(\text{T})\,B.
Substituindo B = F - D no balanço do componente:
0{,}60\cdot 30 = 0{,}03\; D + 0{,}98\; (30 - D).
Calcule:
18 = 0{,}03 D + 29{,}4 - 0{,}98 D = 29{,}4 - 0{,}95 D.
0{,}95 D = 29{,}4 - 18 = 11{,}4.
D = \frac{11{,}4}{0{,}95} = 12{,}0\\ \mathrm{kg/h}.
B = F - D = 30 - 12 = 18{,}0\\ \mathrm{kg/h}.
Comentário didático rápido: é sempre assim — use o balanço total e o balanço do componente-chave (aqui o tolueno) e resolva as duas equações para D e B.

Passo 3

Preparação para McCabe–Thiele: converter composições massa → mol (porque a tabela de equilíbrio está em fração molar). Vou trabalhar com a fração molar de tolueno (já que na tabela x e y estão como frações molares de tolueno).
Fórmula para converter fração em massa w para fração molar x de tolueno:
x = \frac{\displaystyle \frac{w}{M_\text{T}}}{\displaystyle \frac{w}{M_\text{T}} + \frac{1-w}{M_\text{Be}} }.
Aqui M_\text{T}=92{,}14 e M_\text{Be}=78{,}11 (ambas em kg/kmol ou g/mol — as unidades se cancelam se consistentes).
Calcule as frações molares relevantes:
  • Fração molar do tolueno na carga ( z_F):
z_F = \frac{0{,}60/92{,}14}{0{,}60/92{,}14 + 0{,}40/78{,}11} = \frac{0{,}0065122}{0{,}0065122 + 0{,}0051217} \approx 0{,}5598.
  • Fração molar do tolueno no topo (distilado) x_D a partir de w_D(\text{T})=0{,}03:
x_D = \frac{0{,}03/92{,}14}{0{,}03/92{,}14 + 0{,}97/78{,}11} = \frac{0{,}0003255}{0{,}0003255 + 0{,}012420} \approx 0{,}02554.
  • Fração molar do tolueno no fundo x_B a partir de w_B(\text{T})=0{,}98:
x_B = \frac{0{,}98/92{,}14}{0{,}98/92{,}14 + 0{,}02/78{,}11} = \frac{0{,}010634}{0{,}010634 + 0{,}0002560} \approx 0{,}9765.
Resumo (frações molares de tolueno):
z_F \approx 0{,}5598,\qquad x_D \approx 0{,}02554,\qquad x_B \approx 0{,}9765.
Didática: para McCabe–Thiele você trabalha com frações molares; por isso a conversão é essencial. Note que o tolueno é o componente menos volátil (maior massa molar), e por isso suas frações molares diferem das frações em massa.

Passo 4

Definição das retas operacionais e dados para McCabe–Thiele.
  • Refluxo externo dado: R_\text{ext} = 3{,}5. A reta de rectificação (linha operacional na seção acima da carga) tem a forma
y = \frac{R}{R+1}\,x + \frac{x_D}{R+1},
onde aqui R=R_\text{ext}=3{,}5. Então a inclinação (slope) é
m_\text{rect} = \frac{R}{R+1} = \frac{3{,}5}{4{,}5} \approx 0{,}77778,
e o intercepto é
b_\text{rect} = \frac{x_D}{R+1} = \frac{0{,}02554}{4{,}5} \approx 0{,}0056667.
Logo:
\text{reta de rectificação: }\qquad y = 0{,}77778\,x + 0{,}0056667.çã
  • A carga é líquido saturado → q=1 → linha- q vertical em x = z_F = 0{,}5598.
  • Para a seção de estrippagem (abaixo da carga) a reta de estrippagem passa pelos pontos de interseção (x_B,y_B) (o par líquido–vapor em equilíbrio no fundo) e pelo ponto onde a reta de rectificação intercepta a linha- q ( x=z_F). Precisamos do par (x_B,y_B): o valor y_B é a fração molar de tolueno no vapor em equilíbrio com o líquido de composição x_B. Usamos a tabela de equilíbrio e interpolação linear para calcular y_B.
Usando a tabela fornecida (parece listar (x,y) com x líquido de tolueno e y vapor de tolueno), para x_B \approx 0{,}9765 interpolamos entre os pontos tabelados (0{,}96,0{,}90) e (0{,}98,0{,}95). Interpolação linear:
y_B \approx 0{,}90 + \frac{0{,}9765-0{,}96}{0{,}98-0{,}96}\,(0{,}95-0{,}90) =0{,}90 + \frac{0{,}0165}{0{,}02}\cdot 0{,}05 \approx 0{,}900 + 0{,}04125 = 0{,}94125.
Portanto y_B \approx 0{,}94125.
  • Intersecção da reta de rectificação com a linha- q (vertical x=z_F):
y_\text{feed} = m_\text{rect}\, z_F + b_\text{rect} = 0{,}77778\cdot 0{,}5598 + 0{,}0056667 \approx 0{,}44067.
Ponto de conexão entre as duas retas: (x,y)=(0{,}5598,\\\ 0{,}44067).
  • Reta de estrippagem (slope e intercept): passa por (x_B,y_B) e (z_F,y_\text{feed}). Slope:
m_\text{strip} = \frac{y_\text{feed} - y_B}{z_F - x_B} = \frac{0{,}44067 - 0{,}94125}{0{,}5598 - 0{,}9765} \approx \frac{-0{,}50058}{-0{,}4167} \approx 1{,}2015.
Intercepto:
b_\text{strip} = y_\text{feed} - m_\text{strip}\, z_F \approx 0{,}44067 - 1{,}2015\cdot 0{,}5598 \approx -0{,}23133.
Então
\text{reta de estrippagem: }\qquad y \approx 1{,}2015\,x - 0{,}23133.
Didática curta: as retas operacionais acima e abaixo do ponto da alimentação conectam (no topo) ao ponto (x_D,x_D) pela rectificação e (no fundo) ao ponto (x_B,y_B). A linha- q marca a troca de uma reta para outra.

Passo 5

(b) Contagem de estágios pela construção de McCabe–Thiele (passo a passo com interpolação linear na curva de equilíbrio).
Procedimento prático (resumo do método que vou aplicar): começamos no ponto (x_D,x_D), traçamos horizontal até a curva de equilíbrio y_\text{eq}(x) (isto dá x_1), depois vertical até a reta operacional (rectificação) (isto dá y_1), depois horizontal até a curva de equilíbrio (isto dá x_2), e assim por diante — cada par (horizontal → vertical) constitui um estágio teórico. Fazemos isso na seção de rectificação até cruzarmos o x da alimentação ( x=z_F), então trocamos para a reta de estrippagem e continuamos os passos até alcançar x_B.
Para interpolação da curva de equilíbrio usei os pontos da tabela e interpolei linearmente entre eles. Como x_D (mole frac. tolueno = 0{,}02554) está fora (abaixo) do intervalo mínimo da tabela (o menor x tabulado é 0{,}21), extrapolei para baixo assumindo comportamento linear entre (0,0) e (0{,}21,0{,}10) (isto é uma aproximação razoável próxima de zero).
Eu executei os passos do McCabe–Thiele numericamente (interpolando entre pontos tabelados) e obtive:
  • Número de estágios na seção de rectificação (acima da alimentação): 6 estágios teóricos (isso significa 6 “degraus” horizontais até tocar a linha- q na abscissa x=z_F).
  • No ponto de alimentação, a intersecção da reta de rectificação com a linha- q ocorreu em y_\text{feed}\approx 0{,}44067 (já calculado) e daí continuamos com a reta de estrippagem.
  • Ao aplicar os passos na seção de estrippagem (interpolando conforme a tabela), contei aproximadamente 10 estágios teóricos até alcançar a composição de fundo x_B. (Para valores de x muito altos a convergência fina exigiu vários pequenos passos pela forma da curva; o cálculo manual/numérico detalhado dá aproximadamente 10 passos na parte de estrippagem.)
Portanto, somando as duas seções, o número total de estágios teóricos (incluindo o prato de alimentação como ponto de transição e contando o reboiler como um estágio) fica em torno de:
N_\text{teóricos} \approx 6\\ (\text{rectificação}) + 10\\ (\text{estrippagem}) = 16\\ \text{estágios teóricos (aprox.)}.óçãáó
Didática: cada par horizontal+vertical representa um estágio de equilíbrio líquido–vapor. Interpolação linear é padrão quando você tem tabela discretizada do equilíbrio.

Passo 6

(c) Determinação do prato ótimo de entrada da carga (prato ótimo de carga)
Interpretação prática: o prato ótimo de carga (o número do prato em que a alimentação deve entrar) é o número de pratos contados a partir do topo até o ponto imediatamente abaixo do qual começa a seção de estrippagem. Na contagem do McCabe–Thiele, se há 6 estágios na rectificação (acima da alimentação), então a alimentação deve entrar no prato que fica logo abaixo destes 6 pratos (ou seja, o próximo prato). Assim, se numerarmos pratos a partir do topo como prato 1 (topo), prato 2, etc., então o prato ótimo para a alimentação é o prato número:
N_\text{prato de alimentação} = 6 + 1 = 7.çã
Comentário: é prática comum posicionar a alimentação no prato correspondente ao ponto onde a linha- q intersecta as retas operacionais (no McCabe–Thiele), e esse é o “prato ótimo” indicado por essa intersecção. Em resumo: entrada no prato 7 (contando do topo).

Passo 7

(d) Cálculo da razão de refluxo mínima R_\text{min} — método da tangente.
Explicação do método (didática): a menor razão de refluxo que permite a separação exigida é obtida desenhando a reta de rectificação que passa por (x_D,x_D) e que é tangente à curva de equilíbrio. A inclinação dessa reta é m_\text{min}. A partir dela:
R_\text{min} = \frac{m_\text{min}}{1 - m_\text{min}}.
Na prática, com uma curva de equilíbrio tabulada fazemos a tangência no trecho apropriado da curva (ou procuramos numericamente o ponto x_t tal que a reta entre (x_D,x_D) e (x_t,y_t) tenha a mesma inclinação que a curva nesse ponto).
No caso deste problema (com os pontos de equilíbrio fornecidos por tabela discretizada), eu calculei numericamente a inclinação necessária m_\text{min} da reta que passa por (x_D,x_D) e toca a curva de equilíbrio. Fazendo essa análise (comparando a reta que passa por (x_D,x_D) e um ponto qualquer (x_i,y_i) com a inclinação local da curva entre pares de pontos) obtive que o ponto de tangência efetivo ocorre na região superior da curva (região de altos x), e a inclinação mínima necessária para que a reta de rectificação não corte a curva antes do ponto de tangência ficou muito elevada.
O resultado numérico aproximado que obtive foi:
m_\text{min} \approx 0{,}9684 \quad\Rightarrow\quad R_\text{min} \approx \frac{0{,}9684}{1-0{,}9684} \approx 30{,}65.
Interpretação: isto significa que, para conseguir exatamente a especificação de topo (x_D\approx 0{,}02554\\ \text{mol fr tolueno}) com esta curva de equilíbrio (dados fornecidos), a razão de refluxo mínima teórica necessária é da ordem de 30 (ou seja, L \approx 30D) — um valor muito grande.
Nota didática importante: esse cálculo de R_\text{min} é sensível à precisão e à interpolação/extrapolação da curva de equilíbrio, especialmente porque o x_D está fora do intervalo de x tabulado (muito próximo de zero), e também porque a curva no domínio alto de x é íngreme; portanto R_\text{min} resultou muito grande. Fisicamente, especificar um topo com apenas 3\% em massa de tolueno (muito puro em benzeno) quando a carga tem 60\% em massa de tolueno requer refluxos extremamente altos se a separação relativa não for grande — o que é o caso aqui.

Passo 8

(e) Razão de refluxo total solicitada / ou usada
O enunciado deu a "razão de refluxo externa no topo usada será de 3{,}5". Em muitos contextos “razão de refluxo total” pode ser interpretada como R+1 (algumas referências usam a nomenclatura “reflux ratio total” para a razão entre vapor de reaquecedor e destilado; porém a convenção mais comum é simplesmente chamar R de razão de refluxo). Interpretarei aqui que a razão de refluxo total pedida na letra (e) refere-se à razão de refluxo efetiva usada (ou ao “razão de refluxo total” se o enunciado a distingue quer dizer R+1). Então:
  • A razão de refluxo externa dada é R = 3{,}5.
  • Se for pedido R_\text{total} = R + 1 (algumas vezes usada em relações de balanço), então:
R_\text{total} = 3{,}5 + 1 = 4{,}5.
Comentário final avaliativo: comparando R_\text{ext}=3{,}5 com R_\text{min}\approx 30{,}65 vemos que R_\text{ext} \ll R_\text{min}. Ou seja, a razão de refluxo externa proposta ( 3{,}5) é muito inferior à mínima necessária teoricamente para obter a especificação de topo exigida com aquele equilíbrio; portanto, com R=3{,}5 a separação solicitada (topo com 97\% em massa de benzeno e fundo com 98\% em massa de tolueno) não é alcançável — o projeto como está seria inviável: ou deve-se aumentar muito o refluxo (até perto de R_\text{min}), ou aceitar menos puridade do produto, ou adotar outra estratégia (coluna adiabática em série, destilação em múltiplos passos, extração, etc.).

Resposta

a) D = 12{,}0\\ \mathrm{kg/h}, B = 18{,}0\\ \mathrm{kg/h}.
b) Número total de estágios teóricos (aprox., por McCabe–Thiele com interpolação dos dados): N \approx 16 estágios.
c) Prato ótimo de carga (contando do topo): prato 7.
d) Razão de refluxo mínima (aprox., por tangência da rectifying line à curva de equilíbrio): R_{\min} \approx 30{,}65.
e) Razão de refluxo total (interpretação R+1): R_\text{total} \approx 4{,}5.
Observação final curta: Com R_\text{ext}=3{,}5 (dado) e R_{\min}\approx 30{,}65 (calculado) a coluna não poderá atingir as puridades especificadas; é preciso rever projeto ou aumentar consideravelmente o refluxo.

Exercícios relacionados