1) para as equações de calculo da pressão (Antoine, lee-kesler, Wagner, Riedel)

- construa um gráfico comparativo entre os cuatro modelos para uma faixa de temperaturas.
- Escreva a parte textual em markdown e o código em python para o calculo.

NH3

- mostre em outro gráfico os erros de cada modelo, tomando como base experimental os dados do site NIST.

2) compare em outro gráfico P vs T os seguintes compostos:

benzeno, tolueno, xileno. para a equação de Riedel.

Question

1) para as equações de calculo da pressão (Antoine, lee-kesler, Wagner, Riedel)

  - construa um gráfico comparativo entre os cuatro modelos para uma faixa de temperaturas.
  - Escreva  a parte textual em markdown e o código em python para o calculo.
  
NH3

  - mostre em outro gráfico os erros de cada modelo, tomando como base experimental os dados do site NIST.

2) compare em outro gráfico P vs T os seguintes compostos:


benzeno, tolueno, xileno.  para a equação de Riedel.

Accepted Answer

Claro, vamos abordar cada parte da sua solicitação de forma organizada.

### 1) Comparação de Modelos de Pressão de Vapor para NH3

**Parte Textual (Markdown)**

A pressão de vapor de uma substância é uma propriedade termodinâmica crucial, especialmente em engenharia química e processos industriais. Existem vários modelos para estimar a pressão de vapor em função da temperatura, cada um com suas próprias vantagens e limitações. Neste estudo, comparamos quatro modelos comuns: Antoine, Lee-Kesler, Wagner e Riedel, utilizando amônia (NH3) como exemplo.

*   **Equação de Antoine:** É uma das equações mais simples e amplamente utilizadas, comumente expressa como:

$$
    \log_{10}(P) = A - \frac{B}{T + C}
    $$

onde $ P $ é a pressão, $ T $ é a temperatura, e $ A $, $ B $, e $ C $ são constantes empíricas específicas para cada substância.
*   **Equação de Lee-Kesler:** É um método baseado no princípio dos estados correspondentes, utilizando o fator acêntrico para corrigir desvios do comportamento ideal. A equação é mais complexa e geralmente requer dados de propriedades críticas da substância.
*   **Equação de Wagner:** É uma equação mais precisa, frequentemente usada para uma ampla faixa de temperaturas. A forma geral é:

$$
    \ln(P_{r}) = \frac{a\tau + b\tau^{1.5} + c\tau^{3} + d\tau^{6}}{T_{r}}
    $$

onde $ P_{r} $ é a pressão reduzida, $ T_{r} $ é a temperatura reduzida, $ \tau = 1 - T_{r} $, e $ a $, $ b $, $ c $, e $ d $ são constantes.
*   **Equação de Riedel:** É outra equação útil que também considera propriedades críticas:

$$
    \ln(P_{r}) = A + \frac{B}{T_{r}} + C\ln(T_{r}) + DT_{r}^{6}
    $$

onde $ P_{r} $ é a pressão reduzida, $ T_{r} $ é a temperatura reduzida, e $ A $, $ B $, $ C $, e $ D $ são constantes relacionadas às propriedades críticas.

**Código Python**

Aqui está o código Python para calcular a pressão de vapor usando os modelos mencionados e gerar os gráficos comparativos:

```python
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import interp1d

# Dados do NIST (temperatura em K, pressão em kPa)
nist_data = {
    'T': [200, 210, 220, 230, 240, 250, 260, 270, 280, 290, 300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 400],
    'P': [1.78, 4.34, 9.54, 18.89, 34.21, 57.59, 91.42, 138.23, 200.63, 281.24, 383.69, 511.69, 669.97, 863.29, 1096.49, 1374.47, 1692.19, 2054.63, 2466.87, 2934.01, 3460.13]
}

# Dados para NH3
Tc = 405.65  # Temperatura crítica (K)
Pc = 112.8  # Pressão crítica (bar)
omega = 0.250  # Fator acêntrico

# Constantes de Antoine para NH3
A_antoine = 7.21578
B_antoine = 1253.2
C_antoine = 229.7

# Constantes de Riedel para NH3 (exemplo, ajustar conforme necessário)
A_riedel = 0.0
B_riedel = 0.0
C_riedel = 0.0
D_riedel = 0.0

# Constantes de Wagner para NH3 (exemplo, ajustar conforme necessário)
a_wagner = -7.22894
b_wagner = 3.20766
c_wagner = -5.43424
d_wagner = -3.16149

# Função para calcular a pressão de vapor usando a equação de Antoine
def antoine(T):
    return 10**(A_antoine - B_antoine / (T + C_antoine))

# Função para calcular a pressão de vapor usando a equação de Lee-Kesler
def lee_kesler(T, Tc, Pc, omega):
    Tr = T / Tc
    Pr = 10**(5.92714 - 6.09648/Tr - 1.28862*np.log(Tr) + 0.169347*Tr**6) + omega*10**(15.2518 - 15.6875/Tr - 13.4721*np.log(Tr) + 0.43577*Tr**6)
    P = Pr * Pc
    return P

# Função para calcular a pressão de vapor usando a equação de Wagner
def wagner(T, Tc, Pc, a, b, c, d):
    Tr = T / Tc
    tau = 1 - Tr
    Pr = np.exp((a*tau + b*tau**1.5 + c*tau**3 + d*tau**6) / Tr)
    P = Pr * Pc
    return P

# Função para calcular a pressão de vapor usando a equação de Riedel
def riedel(T, Tc, Pc, A, B, C, D):
    Tr = T / Tc
    Pr = np.exp(A + B/Tr + C*np.log(Tr) + D*Tr**6)
    P = Pr * Pc
    return P

# Faixa de temperatura
T = np.linspace(200, 400, 100)

# Calcular a pressão de vapor usando cada modelo
P_antoine = antoine(T)
P_lee_kesler = lee_kesler(T, Tc, Pc, omega)
P_wagner = wagner(T, Tc, Pc, a_wagner, b_wagner, c_wagner, d_wagner)
P_riedel = riedel(T, Tc, Pc, A_riedel, B_riedel, C_riedel, D_riedel)

# Converter dados do NIST para array numpy
T_nist = np.array(nist_data['T'])
P_nist = np.array(nist_data['P'])

# Interpolação linear dos dados do NIST
nist_interp = interp1d(T_nist, P_nist, kind='linear', fill_value="extrapolate")
P_nist_interp = nist_interp(T)

# Calcular os erros
erro_antoine = np.abs((P_antoine - P_nist_interp) / P_nist_interp) * 100
erro_lee_kesler = np.abs((P_lee_kesler - P_nist_interp) / P_nist_interp) * 100
erro_wagner = np.abs((P_wagner - P_nist_interp) / P_nist_interp) * 100
erro_riedel = np.abs((P_riedel - P_nist_interp) / P_nist_interp) * 100

# Plotar as curvas de pressão de vapor
plt.figure(figsize=(12, 6))
plt.plot(T, P_antoine, label='Antoine')
plt.plot(T, P_lee_kesler, label='Lee-Kesler')
plt.plot(T, P_wagner, label='Wagner')
plt.plot(T, P_riedel, label='Riedel')
plt.plot(T_nist, P_nist, 'o', label='NIST Data')
plt.xlabel('Temperatura (K)')
plt.ylabel('Pressão (kPa)')
plt.title('Comparação de Modelos de Pressão de Vapor para NH3')
plt.legend()
plt.grid(True)
plt.show()

# Plotar os erros
plt.figure(figsize=(12, 6))
plt.plot(T, erro_antoine, label='Erro Antoine')
plt.plot(T, erro_lee_kesler, label='Erro Lee-Kesler')
plt.plot(T, erro_wagner, label='Erro Wagner')
plt.plot(T, erro_riedel, label='Erro Riedel')
plt.xlabel('Temperatura (K)')
plt.ylabel('Erro Absoluto Relativo (%)')
plt.title('Erros dos Modelos de Pressão de Vapor para NH3')
plt.legend()
plt.grid(True)
plt.show()
```

Este código realiza as seguintes etapas:

1.  **Importa as bibliotecas necessárias:** `numpy` para cálculos numéricos, `matplotlib` para plotagem e `scipy.interpolate` para interpolação dos dados do NIST.
2.  **Define os dados do NIST:** Temperaturas e pressões experimentais.
3.  **Define as constantes:** Propriedades críticas e constantes das equações para NH3.
4.  **Define as funções para calcular a pressão de vapor:** Implementa as equações de Antoine, Lee-Kesler, Wagner e Riedel.
5.  **Calcula a pressão de vapor:** Para uma faixa de temperaturas usando cada modelo.
6.  **Interpola os dados do NIST:** Para comparar com os modelos.
7.  **Calcula os erros:** Comparando os resultados dos modelos com os dados do NIST.
8.  **Plota os resultados:** Cria dois gráficos, um comparando as curvas de pressão de vapor e outro mostrando os erros de cada modelo.

### 2) Comparação de P vs T para Benzeno, Tolueno e Xileno usando a Equação de Riedel

```python
import numpy as np
import matplotlib.pyplot as plt

# Propriedades críticas
propriedades = {
    'benzeno': {'Tc': 562.2, 'Pc': 48.98, 'A': 0.0, 'B': 0.0, 'C': 0.0, 'D': 0.0},
    'tolueno': {'Tc': 591.8, 'Pc': 41.06, 'A': 0.0, 'B': 0.0, 'C': 0.0, 'D': 0.0},
    'xileno': {'Tc': 630.3, 'Pc': 37.00, 'A': 0.0, 'B': 0.0, 'C': 0.0, 'D': 0.0}
}

# Equação de Riedel
def riedel(T, Tc, Pc, A, B, C, D):
    Tr = T / Tc
    Pr = np.exp(A + B/Tr + C*np.log(Tr) + D*Tr**6)
    P = Pr * Pc
    return P

# Faixa de temperatura
T = np.linspace(300, 600, 100)

# Calcular a pressão de vapor para cada composto
P_benzeno = riedel(T, propriedades['benzeno']['Tc'], propriedades['benzeno']['Pc'], propriedades['benzeno']['A'], propriedades['benzeno']['B'], propriedades['benzeno']['C'], propriedades['benzeno']['D'])
P_tolueno = riedel(T, propriedades['tolueno']['Tc'], propriedades['tolueno']['Pc'], propriedades['tolueno']['A'], propriedades['tolueno']['B'], propriedades['tolueno']['C'], propriedades['tolueno']['D'])
P_xileno = riedel(T, propriedades['xileno']['Tc'], propriedades['xileno']['Pc'], propriedades['xileno']['A'], propriedades['xileno']['B'], propriedades['xileno']['C'], propriedades['xileno']['D'])

# Plotar as curvas de pressão de vapor
plt.figure(figsize=(12, 6))
plt.plot(T, P_benzeno, label='Benzeno')
plt.plot(T, P_tolueno, label='Tolueno')
plt.plot(T, P_xileno, label='Xileno')
plt.xlabel('Temperatura (K)')
plt.ylabel('Pressão (bar)')
plt.title('Comparação de Pressão de Vapor para Benzeno, Tolueno e Xileno (Equação de Riedel)')
plt.legend()
plt.grid(True)
plt.show()
```

Este código:

1.  **Define as propriedades críticas:** Temperatura e pressão críticas para benzeno, tolueno e xileno.
2.  **Implementa a equação de Riedel.**
3.  **Calcula a pressão de vapor:** Para cada composto em uma faixa de temperaturas.
4.  **Plota os resultados:** Mostrando a comparação das curvas de pressão de vapor.

Lembre-se de ajustar as constantes das equações (A, B, C, D) para obter resultados mais precisos. As constantes fornecidas são apenas exemplos e podem não ser adequadas para todos os casos.

Perguntas relacionadas