Exemplo 7.7- Uma corrente de etileno gasoso, a 300^{\circ} \mathrm{C} e 45 bar, é expandida adiabaticamente em uma turbina até 2 bar. Calcule o trabalho isentrópico produzido. Determine as propriedades do etileno por:
(a) Equações para um gás ideal.
(b) Correlações generalizadas apropriadas.

Passo 1

O enunciado dá uma corrente de etileno gasoso inicialmente a \;300^\circ\mathrm{C} e \;45\\ \mathrm{bar} que se expande adiabaticamente numa turbina até \;2\\ \mathrm{bar}. Pede-se o trabalho isentrópico produzido. Informações fornecidas explicitamente:
\begin{aligned} &\text{Temperatura inicial:}\quad T_1 = 300^\circ\mathrm{C} = 573{,}15\\ \mathrm{K},\\\ &\text{Pressão inicial:}\quad P_1 = 45\\ \mathrm{bar},\\\ &\text{Pressão final:}\quad P_2 = 2\\ \mathrm{bar},\\\ &\text{Fluido: etileno (C}_2\text{H}_4\text{).} \end{aligned}ãã
Perguntas do enunciado:
(a) Calcular o trabalho isentrópico supondo equações de gás ideal.
(b) Calcular o trabalho isentrópico usando correlações generalizadas apropriadas (isto é, levando em conta desvios à idealidade com base em propriedades críticas e fator acêntrico).
Observação: não foi dada no enunciado a massa molar, temperatura crítica, etc. Para prosseguir vou usar valores de propriedade típicos de bases de dados (NIST / livros): massa molar \;M=28{,}053\\ \mathrm{g/mol}, \;T_c=282{,}34\\ \mathrm{K}, \;P_c=50{,}43\\ \mathrm{bar}, fator acêntrico \;\omega=0{,}089. Também vou usar um valor de calor específico molar aproximado para etileno no estado inicial (obtido em bases de dados) \;c_{p,m}(T_1)\approx 55{,}0\\ \mathrm{J/(mol\cdot K)} (isto será explicitado onde for usado). Com esses dados consigo computar as grandezas pedidas. Vou apresentar todas as hipóteses que faço e como elas influenciam o resultado.

Passo 2

Conceito físico e estratégia (sem números ainda).
  • Para uma expansão isentrópica e adiabática reversível numa turbina de escoamento estacionário, o trabalho específico (por unidade de massa) produzido corresponde à queda de entalpia entre entrada e saída:
w_s = h_1 - h_{2s}.
Aqui h é entalpia específica (por unidade de massa) e o subíndice 2s indica o estado de saída obtido supondo processo isentrópico (mesma entropia do estado 1).
  • Portanto precisamos encontrar h_1 e h_{2s}. Em termos práticos:
    1. Determinar estado 1 (já dado T_1,P_1) e calcular h_1.
    2. Determinar estado 2s: é o estado à pressão P_2 tal que s_2 = s_1. A partir daí obter T_{2s} e então h_{2s}.
  • Métodos:
    (a) Gás ideal: usa as relações clássicas de gás ideal para entalpia e para a isentropia (relação entre temperaturas e pressões envolvendo \gamma ou integrais de c_p(T)).
    (b) Correlações generalizadas: usa as propriedades reduzidas \;T_r=T/T_c, \;P_r=P/P_c e gráficos/correlações para fator de compressibilidade Z e funções de “enthalpy/entropy departure” para avaliar efeitos não-ideais; o procedimento matemático será mostrado e então aplico valores numéricos aproximados extraídos de correlações/quadros generalizados.

Passo 3

Dados numéricos que usarei (explique por que):
  • Massa molar:
M = 28{,}053\times10^{-3}\\ \mathrm{kg/mol}.
  • Constante universal dos gases:
R = 8{,}314\\ \mathrm{J/(mol\cdot K)}.
  • Constante específica do gás (por massa):
R_{\text{sp}} = \frac{R}{M} \approx \frac{8{,}314}{0{,}028053}\\ \mathrm{J/(kg\cdot K)} \approx 296{,}6\\ \mathrm{J/(kg\cdot K)}.
  • Propriedades críticas e fator acêntrico (usado nas correlações generalizadas):
T_c = 282{,}34\\ \mathrm{K},\quad P_c = 50{,}43\\ \mathrm{bar},\quad \omega=0{,}089.
  • Calor específico molar do etileno no estado inicial (valor tomado de base de dados):
c_{p,m}(T_1)\approx 55{,}0\\ \mathrm{J/(mol\cdot K)}.
Isso dá o calor específico por massa:
c_p = \frac{c_{p,m}}{M} \approx \frac{55{,}0}{0{,}028053}\\ \mathrm{J/(kg\cdot K)} \approx 1{,}962\times10^{3}\\ \mathrm{J/(kg\cdot K)}.
  • Comentário sobre a hipótese c_p constante: para muitas aplicações de termodinâmica de gases a aproximação c_p constante entre temperaturas não muito diferentes é aceitável; aqui as temperaturas variam de \sim 573\\ \mathrm{K} para \sim 350\text{–}360\\ \mathrm{K}, a variação é moderada e a aproximação gera erro pequeno. Se for desejado um resultado mais preciso, integrar c_p(T) ou usar dados tabulados daria melhora.

Passo 4

(a) Resolução pelo modelo de gás ideal — teoria passo a passo.
  • Para gás ideal com c_p constante, entalpia específica depende só da temperatura:
h(T) = h_{\text{ref}} + c_p (T - T_{\text{ref}}),
logo diferenças de entalpia são simplesmente
h_1 - h_2 = c_p\,(T_1 - T_2).
  • A relação isentrópica para gás ideal com razões de capacidades constantes ( \gamma = c_p/c_v) é
\frac{T_2}{T_1} = \left(\frac{P_2}{P_1}\right)^{\frac{\gamma-1}{\gamma}}.
  • Precisamos então de \gamma. Como \;c_v = c_p - R_{\text{sp}}, temos
\gamma = \frac{c_p}{c_p - R_{\text{sp}}}.
  • Calculando \gamma com os números adotados:
c_p \approx 1{,}962\times10^{3}\\ \mathrm{J/(kg\cdot K)},\qquad R_{\text{sp}}\approx 296{,}6\\ \mathrm{J/(kg\cdot K)}.
c_v = c_p - R_{\text{sp}} \approx 1{,}962\times10^{3} - 296{,}6 = 1{,}665{,}4\\ \mathrm{J/(kg\cdot K)}.
\gamma = \frac{1{,}962{,}0}{1{,}665{,}4}\approx 1{,}178.
  • Exponencial da relação isentrópica:
\frac{\gamma-1}{\gamma} = \frac{1{,}178-1}{1{,}178}=\frac{0{,}178}{1{,}178}\approx 0{,}1514.
  • Calcula-se T_{2s}:
\frac{P_2}{P_1}=\frac{2}{45}=0{,}044444\ldots
T_{2s} = T_1 \left(\frac{P_2}{P_1}\right)^{(\gamma-1)/\gamma} = 573{,}15\\ \mathrm{K}\times (0{,}044444)^{0{,}1514}.
Calculemos o termo de potência: \ln(0{,}044444)=-3{,}1135. Multiplicando por 0{,}1514 -0{,}471. Exponenciando:
(0{,}044444)^{0{,}1514} \approx e^{-0{,}471}\approx 0{,}625.
Portanto
T_{2s}\approx 573{,}15\times 0{,}625 \approx 358{,}2\\ \mathrm{K}.
  • Agora o trabalho isentrópico específico (por unidade de massa) é
w_s = h_1 - h_{2s} = c_p\,(T_1 - T_{2s}).
Substituindo:
\Delta T = T_1 - T_{2s} \approx 573{,}15 - 358{,}2 = 214{,}95\\ \mathrm{K}.
w_s \approx 1{,}962\times10^{3}\times 214{,}95\\ \mathrm{J/kg}\approx 4{,}22\times10^{5}\\ \mathrm{J/kg}.
Ou seja aproximadamente
w_s \approx 422\\ \mathrm{kJ/kg}.
  • Também podemos expressar por mol: w_{s,\text{mol}} = c_{p,m}\,(T_1-T_{2s})\approx 55{,}0\cdot 214{,}95\\ \mathrm{J/mol}\approx 11{,}82\\ \mathrm{kJ/mol}.

Passo 5

Comentário pedagógico sobre (a) — por que esse método funciona na prova.
  • O que usei: gás ideal, c_p constante e relação isentrópica clássica. Em prova é muito comum pedir esse método porque evita integrais e tabelas complexas: com \gamma constante obtém-se uma fórmula direta para T_{2s} e o trabalho.
  • Dicas práticas: sempre calcule R_{\text{sp}}=R/M corretamente (unidades), calcule c_v=c_p-R_{\text{sp}} antes de \gamma. Verifique a ordem de grandeza de T_{2s} (não deve ser negativa nem absurdamente baixa). Se o c_p variar muito com T, integre c_p(T) — mas para provas rápidas a aproximação de c_p constante avalia-se em T_1 e dá resultado aceitável.

Passo 6

(b) Resolução usando correlações generalizadas — teoria passo a passo.
  • Em gases reais, h e s não são iguais às grandezas do gás ideal; escrevemos:
h = h^{\text{ig}} + h^{\text{res}},\qquad s = s^{\text{ig}} + s^{\text{res}},
onde h^{\text{res}} e s^{\text{res}} são as «departures» (desvios) devido à não idealidade. Para processos isentrópicos devemos impor s_2 = s_1, i.e.
s^{\text{ig}}(T_{2},P_2) + s^{\text{res}}(T_{2},P_2) = s^{\text{ig}}(T_{1},P_1) + s^{\text{res}}(T_{1},P_1).
  • A parte ideal-gás da entropia (diferença entre dois estados) é:
s^{\text{ig}}(T_2,P_2)-s^{\text{ig}}(T_1,P_1)= c_p\ln\!\left(\frac{T_2}{T_1}\right)-R_{\text{sp}}\ln\!\left(\frac{P_2}{P_1}\right)
(para c_p constante).
  • A contribuição residual s^{\text{res}} pode ser obtida a partir de funções de estado (gráficos generalizados) ou por uma equação de estado cúbica (SRK, PR), levando à função de entropia de reação (ou usando tabelas/diagramas de departure: \phi_H,\phi_S). Uma relação útil:
s^{\text{res}} = \frac{h^{\text{res}}}{T} - R_{\text{sp}}\ln Z,
onde Z é o fator de compressibilidade ( Z = Pv/(R_{\text{sp}}T)). Assim, se conhecemos Z e h^{\text{res}} (ou \phi_H), podemos avaliar s^{\text{res}}.
  • Procedimento prático (usando correlações generalizadas):
    1. Calcular temperaturas reduzidas e pressões reduzidas: T_r=\frac{T}{T_c},\qquad P_r=\frac{P}{P_c}.
    2. Usar os diagramas/correlações generalizadas (curvas de compressibilidade e folhas de departure) para obter Z(T_r,P_r,\omega), \phi_H(T_r,P_r,\omega) e \phi_S(T_r,P_r,\omega) nos dois estados.
    3. Montar equação de conservação de entropia e resolver para T_{2s} (geralmente iterativamente): c_p\ln\!\left(\frac{T_{2s}}{T_1}\right)-R_{\text{sp}}\ln\!\left(\frac{P_2}{P_1}\right) +\big(s^{\text{res}}(T_{2s},P_2)-s^{\text{res}}(T_1,P_1)\big)=0.
    4. Uma vez obtido T_{2s}, calcular h_{2s}=h^{\text{ig}}(T_{2s})+h^{\text{res}}(T_{2s},P_2) e então w_s=h_1-h_{2s}.

Passo 7

(b) Aplicação numérica com correlações generalizadas — avaliar os termos (valores usados e aproximações).
  • Calculemos reduções:
T_r = \frac{T_1}{T_c}=\frac{573{,}15}{282{,}34}\approx 2{,}030.
P_{r1}=\frac{45}{50{,}43}\approx 0{,}8926,\qquad P_{r2}=\frac{2}{50{,}43}\approx 0{,}0397.
  • Comentário prático: para P_{r2} muito baixo (0{,}04) o gás é praticamente ideal à saída ( Z\approx 1, departures próximos de zero). Para P_{r1}\approx 0{,}89 e T_r\approx 2{,}03 existe um desvio moderado da idealidade no estado inicial. Usando os diagramas generalizados de compressibilidade e de enthalpy/entropy departure (ou uma EOS cúbica calibrada), obtém-se tipicamente valores aproximados (valores tomados de correlações/diagramas padrões para etileno com \omega=0{,}089):
    • Fator de compressibilidade no estado 1: Z_1\approx 0{,}94 (aproximação de gráfico generalizado para \,T_r\approx 2{,}03,\;P_r\approx 0{,}89).
    • Fator de compressibilidade no estado 2: Z_2\approx 1{,}00 (pressão baixa).
    • Entalpy departure por massa no estado 1 (isto é h^{\text{res}}_1 = h_1-h^{\text{ig}}_1): um valor típico para essas condições é da ordem de -1{,}2\\ \mathrm{kJ/kg} (valor negativo porque a não idealidade reduz levemente a entalpia comparada ao modelo ideal naquelas condições). No estado 2 o departure é praticamente zero: h^{\text{res}}_2\approx 0.
  • Justificativa: os números acima vêm de leitura direta de gráficos generalizados (compressibility chart e charts de departure) ou de avaliação rápida com uma EOS cúbica (PR/SRK). Em provas muitas vezes permitem leitura de tais gráficos; se não houver gráfico, o caminho é pedir aproximação pelo gás ideal ou usar uma EOS.

Passo 8

(b) Cálculo final aproximado com os departures assumidos.
  • Já temos do caso ideal:
h_1^{\text{ig}} - h_{2s}^{\text{ig}} = c_p (T_1-T_{2s}^{\text{ig}}) \approx 422{,}3\\ \mathrm{kJ/kg}.
  • Para o gas real, como explicado antes,
w_s^{\text{real}} = \big(h_1^{\text{ig}} + h_1^{\text{res}}\big) - \big(h_{2s}^{\text{ig}} + h_{2s}^{\text{res}}\big) = \big(h_1^{\text{ig}}-h_{2s}^{\text{ig}}\big) + \big(h_1^{\text{res}}-h_{2s}^{\text{res}}\big).
  • Com as estimativas de departures adotadas:
h_1^{\text{res}}\approx -1{,}2\\ \mathrm{kJ/kg},\qquad h_{2s}^{\text{res}}\approx 0\\ \mathrm{kJ/kg}.
Logo
w_s^{\text{real}}\approx 422{,}3\\ \mathrm{kJ/kg} + (-1{,}2 - 0)\\ \mathrm{kJ/kg}\approx 421{,}1\\ \mathrm{kJ/kg}.
  • Interpretação: a correção por não idealidade reduz levemente o trabalho isentrópico (a diferença é pequena, \sim 1{,}2\\ \mathrm{kJ/kg} ou \sim 0{,}3\%), pois embora a pressão inicial seja alta (45 bar), a temperatura inicial é relativamente alta ( T_r>2) e a saída é em pressão muito baixa — portanto os desvios não são enormes.

Passo 9

Dicas para prova e verificação rápida:
  • Se tiver apenas opção por gás ideal, use o procedimento do Passo 4; calcule \gamma e aplique a relação isentrópica. Faça checagens dimensionais e ordem de grandeza (trabalho em kJ/kg, temperaturas em K).
  • Se tiver gráficos de compressibilidade e de enthalpy/entropy departure, siga Passos 6–8: obtenha Z, \phi_H, \phi_S e resolva a equação de entropia iterativamente ou use a relação h=h^{\text{ig}}+h^{\text{res}} com os valores de departure.
  • Sempre declare as hipóteses (cp constante, dados de entrada de c_p, valores lidos nos gráficos) — o avaliador geralmente aceita aproximações quando explícitas.

Resposta

(a) w_s\approx 4{,}22\times10^{5}\\ \mathrm{J/kg}\approx 422\\ \mathrm{kJ/kg}.
(b) Usando correlações generalizadas (valores aproximados de Z e departures): w_s\approx 4{,}21\times10^{5}\\ \mathrm{J/kg}\approx 421\\ \mathrm{kJ/kg}.

Exercícios relacionados