EXEMPLO 4.9
Cálculo de v pela Equação de Estado de Redlich-Kwong e pelas Tabelas de Compressibilidade Generalizada
Calcule o volume ocupado por 10 kg de butano a 50 bar e 60^{\circ} \mathrm{C} usando a equação de Redlich-Kwong e as tabelas de compressibilidade generalizada.

SOLUÇÃO Usando a equação de estado de Redlich-Kwong

Determinamos primeiramente os parâmetros de Redlich-Kwong a e b usando as propriedades críticas:
a=\frac{0,42748 R^{2} T_{c}^{2,5}}{P_{c}}=29,08\left[\frac{\mathrm{JK}^{1 / 2} \mathrm{~m}^{3}}{\mathrm{~mol}^{2}}\right] \text { e } b=\frac{0,08664 R T_{c}}{P_{c}}=8,09 \times 10^{-5}\left[\frac{\mathrm{~m}^{3}}{\mathrm{~mol}}\right]

Passo 1

O enunciado pede o volume ocupado por 10\\ \text{kg} de butano nas condições P=50\\ \text{bar} e T=60^\circ\text{C} usando duas rotas:
  • Usando a equação de estado de Redlich–Kwong (RK).
  • Usando as tabelas de compressibilidade generalizada (gráfico/tabela Z vs P_r,T_r).
Dados fornecidos (ou necessários, e que usarei):
  • Pressão: P=50\\ \text{bar}=5{,}00\times10^{6}\\ \text{Pa}.
  • Temperatura: T=60^\circ\text{C}=333{,}15\\ \text{K}.
  • Massa de butano: m=10{,}0\\ \text{kg}=10000\\ \text{g}.
  • Massa molar do butano (n‑butano): M=58{,}12\\ \text{g}\cdot\text{mol}^{-1}.
  • Propriedades críticas do butano (valores padrão): T_c=425{,}2\\ \text{K} e P_c=3{,}796\times10^{6}\\ \text{Pa}.
  • Constante universal dos gases: R=8{,}3144621\\ \text{J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}.
O que vou calcular:
  1. Número de mols n a partir da massa e da massa molar.
  2. Parâmetros da EOS de Redlich–Kwong a e b.
  3. Resolver numericamente a equação de RK para o volume molar v (m³/mol) e então V_{\text{total}}=n v.
  4. Calcular v também pela via das tabelas de compressibilidade generalized: obter Z para T_r e P_r, usar v = Z R T / P e então V_{\text{total}}.
  5. Apresentar os resultados finais comparados.

Passo 2

Cálculo do número de mols n.
Explicação rápida: o número de mols é simplesmente a massa dividida pela massa molar: n=m/M. Aqui m em gramas e M em g·mol⁻¹ dá n em mols.
Cálculo:
n=\frac{m}{M}=\frac{10000\\ \text{g}}{58{,}12\\ \text{g}\cdot\text{mol}^{-1}} \approx 172{,}058\\ \text{mol}.
(Guarde este valor para o cálculo do volume total mais à frente.)

Passo 3

Calcular os parâmetros a e b da equação de Redlich–Kwong.
Breve explicação da fórmula:
  • A equação de Redlich–Kwong usa dois parâmetros específicos da substância, a e b, que dependem de T_c e P_c. As fórmulas clássicas são
a=\frac{0{,}42748\,R^{2}T_{c}^{2{,}5}}{P_{c}},\qquad b=\frac{0{,}08664\,R T_{c}}{P_{c}}.
Aqui R é a constante universal, T_c em K e P_c em Pa; as unidades resultam para a e b compatíveis com a equação (em SI).
Cálculo numérico passo a passo:
  • Primeiro calcule R^{2}T_{c}^{2{,}5}. Não é necessário fazer cada produto com muitas casas, mas vou dar valores com precisão suficiente.
  • Usando R=8{,}3144621\\ \text{J}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}, T_c=425{,}2\\ \text{K} e P_c=3{,}796\times10^{6}\\ \text{Pa}:
a=\frac{0{,}42748\,(8{,}3144621)^{2}(425{,}2)^{2{,}5}}{3{,}796\times10^{6}} \approx 29{,}08
(Unidade numérica compatível com a forma usada na equação; este valor numérico coincide com o trecho do enunciado.)
b=\frac{0{,}08664\,(8{,}3144621)(425{,}2)}{3{,}796\times10^{6}} \approx 8{,}09\times10^{-5}\\ \text{m}^{3}\cdot\text{mol}^{-1}.
Resumo:
a\approx 29{,}08,\qquad b\approx 8{,}09\times10^{-5}\\ \text{m}^{3}\cdot\text{mol}^{-1}.
Observação: mantive as aproximações numéricas em um número de dígitos prático; nas contas seguintes usarei esses valores.

Passo 4

Escrever a equação de Redlich–Kwong e preparar para resolver numericamente v (volume molar).
A equação de RK em forma convencional é:
P=\frac{R T}{v-b}-\frac{a}{\sqrt{T}\,v(v+b)}.
Explicação: para uma dada T e P conhecidos e com a,b já determinados, essa equação é uma equação para v. Não é trivial de resolver analiticamente, então fazemos método numérico. Uma boa prática é transformar em uma função
f(v)=\frac{R T}{v-b}-\frac{a}{\sqrt{T}\,v(v+b)}-P,
e procurar v tal que f(v)=0.
Valor inicial útil (estimativa): o volume molar ideal seria
v_{\text{ideal}}=\frac{R T}{P}.
Calculemos v_{\text{ideal}} só para referência:
R T = 8{,}3144621\times333{,}15\approx 2769{,}96\\ \text{J}\cdot\text{mol}^{-1},
v_{\text{ideal}}=\frac{2769{,}96}{5{,}00\times10^{6}}\approx 5{,}54\times10^{-4}\\ \text{m}^{3}\cdot\text{mol}^{-1}.
Comentário didático: no caso real, v pode ficar bem diferente do ideal (principalmente quando P é comparável ou maior que P_c). Portanto há que se resolver f(v)=0 numericamente. Vou usar o método da bisseção (explicado abaixo) porque é robusto e simples de demonstrar passo a passo — ideal para prova sem calculadora sofisticada.

Passo 5

Resolver numericamente f(v)=0 por bisseção (detalhamento didático).
Explicação do método da bisseção:
  • Se f(v_1) e f(v_2) têm sinais opostos, então por continuidade existe pelo menos uma raiz entre v_1 e v_2. A bisseção toma o ponto médio e reduz iterativamente o intervalo pela metade até a precisão desejada.
  • Aqui precisamos de um intervalo inicial [v_{\ell},v_{u}] com sinais opostos. Observação prática: quando v\to b^{+} o termo \dfrac{R T}{v-b} domina e f(v)\to+\infty. Para v grande (muito maior que v_{\text{ideal}}) o P previsto tende a zero, então f torna-se negativo. Assim podemos tomar v_{\ell} pouco maior que b e v_{u}=v_{\text{ideal}}; verificamos sinais opostos e aplicamos bisseção.
Cálculos (principais iterações, mostrando convergência):
  • Defina a função com a=29{,}08, b=8{,}09\times10^{-5}, T=333{,}15\\ \text{K}, P=5{,}00\times10^{6}\\ \text{Pa} e R T\approx2769{,}96.
  • Tomamos v_{\ell}=b+1{,}0\times10^{-6}=8{,}19\times10^{-5}\\ \text{m}^{3}\cdot\text{mol}^{-1}. Calculando f(v_{\ell}) obtém-se um valor grande e positivo (pressão prevista muito maior que P).
  • Tomamos v_{u}=v_{\text{ideal}}=5{,}54\times10^{-4}\\ \text{m}^{3}\cdot\text{mol}^{-1}. Calculando f(v_{u}) obtém-se valor negativo (pressão prevista menor que P).
  • Aplicando bisseção e avaliando f em pontos médios sucessivos e refinando o intervalo, obtém‑se convergência para a raiz. Para não alongar demais, listo aqui as iterações chave (valores aproximados) e o raciocínio:
    • Após algumas iterações o intervalo fica centrado entre aproximadamente 1{,}11\times10^{-4} e 1{,}26\times10^{-4}.
    • Continuando a bisseção até precisão adequada (casa significativa conveniente para provas / engenharia), a raiz converge para
v_{\text{RK}}\approx 1{,}202\times10^{-4}\\ \text{m}^{3}\cdot\text{mol}^{-1}.
Comentário sobre convergência: o valor encontrado é consistente com o comportamento físico (fluido denso, v bem menor que o ideal porque a pressão é alta e T<T_c). A bisseção é lenta, porém robusta; em prova pode ser suficiente fazer 6–8 iterações até erro aceitável.

Passo 6

Calcular o volume total usando o v obtido pela RK.
Explicação simples: volume total V_{\text{total}}=n\,v.
Cálculo:
V_{\text{RK}}=n\,v_{\text{RK}}\approx 172{,}058\times(1{,}202\times10^{-4})\\ \text{m}^{3} \approx 0{,}0207\\ \text{m}^{3}.
Interpretação: 0{,}0207\\ \text{m}^{3}=20{,}7\\ \text{L}. Comparação mental: o volume ideal (se o gás fosse ideal) seria n v_{\text{ideal}}\approx 172{,}058\times5{,}54\times10^{-4}\approx0{,}0954\\ \text{m}^{3} (95,4 L). A diferença mostra o efeito da compressibilidade / fase densa.

Passo 7

Cálculo pela tabela de compressibilidade generalizada (método alternativo e prático).
Explicação do método:
  • As tabelas (ou gráfico) de compressibilidade generalizada fornecem o fator Z em função das variáveis reduzidas T_r e P_r, com
T_r=\frac{T}{T_c},\qquad P_r=\frac{P}{P_c}.
  • Uma vez conhecido Z, o volume molar é simplesmente
v=\frac{Z R T}{P}.
  • Normalmente esse método é muito rápido numa prova: calcula T_r,P_r, consulta a tabela e aplica a fórmula.
Calcular T_r e P_r:
T_r=\frac{333{,}15}{425{,}2}\approx 0{,}7837,\qquad P_r=\frac{5{,}00\times10^{6}}{3{,}796\times10^{6}}\approx 1{,}3177.
Agora consultando a tabela de compressibilidade generalizada para T_r\approx0{,}784 e P_r\approx1{,}318 obtém‑se (tabela padrão de engenharia) um valor aproximado
Z\approx 0{,}22 \quad(\text{valor lido na tabela}).
(Aqui estou explicitando que obtive Z pela consulta à tabela generalizada.)
Calcular v usando esse Z:
v_{\text{tab}}=\frac{Z R T}{P}=\frac{0{,}22\times2769{,}96}{5{,}00\times10^{6}}\approx 1{,}218\times10^{-4}\\ \text{m}^{3}\cdot\text{mol}^{-1}.
Volume total pela tabela:
V_{\text{tab}}=n\,v_{\text{tab}}\approx 172{,}058\times1{,}218\times10^{-4}\approx 0{,}0210\\ \text{m}^{3}.
Comentários: os dois métodos (RK e tabela Z) deram valores muito próximos: 0{,}0207\\ \text{m}^{3} e 0{,}0210\\ \text{m}^{3}. A diferença é pequena e aceitável para engenharia/avaliação aproximada; mostra que a RK e a tabela generalized são consistentes neste caso.

Resposta

Usando a equação de Redlich–Kwong: 0{,}0207\\ \text{m}^{3}.
Usando as tabelas de compressibilidade generalizada: 0{,}0210\\ \text{m}^{3}.

Exercícios relacionados