Termodinâmica para

Engenharia Química II
EQA5342
Turmas 06215 e 06216

Aula 11: Programas para Cálculos de

Propriedades Termodinâmicas de
Substâncias Puras e Misturas

Prof. Marcelo Lanza

Profa. Sandra R. S. Ferreira
 TIPOS DE CÁLCULOS

Basicamente, o objetivo dessa aula é mostrar três tipos de cálculos por

meio de programas:

1) CRIAR PROGRAMA PARA O CÁLCULO DA PRESSÃO DE

VAPOR DE UMA SUBSTÂNCIA PURA UTILIZANDO A EQUAÇÃO
DE ESTADO DE PENG-ROBINSON

2) CÁLCULO DE DADOS PVT E PROPRIEDADES

TERMODINÂMICAS DE SUBSTÂNCIAS PURAS UTILIZANDO
PROGRAMAS

3) CÁLCULO DO ELV DE MISTURAS (BINÁRIAS) POR MEIO DE

EQUAÇÕES DE ESTADO E/OU MODELOS DE GAMA
IMPLEMENTADAS EM PROGRAMAS OPEN SOURCE

https://www.sciencedirect.com/science/article/pii/S1749772811000091
 ATIVIDADE 1: PROBLEMA INDIVIDUAL

Exercício Individual: Uso de Programas Termodinâmicos

Realizar cálculos para o sistema contendo componente (1)

+ componente (2) na Temperatura T (°C) e Pressão P
(kPa) conforme será indicado na tabela.

Apresentar a resolução do exercício na forma de um

relatório incluindo prints, imagens e os arquivos eletrônicos
gerados nos cálculos.
 ATIVIDADE 2: TRABALHO EM GRUPO

GRUPOS: número de participantes definidos em cada turma.

O que fazer:

- Utilizar artigos com dados experimentais;

- Utilizar software para os cálculos de equilíbrio líquido-vapor (ELV) de
misturas binárias;
- Elaborar apresentação visual (por ex.: Power Point) com visão crítica
dos resultados;
- Apresentar seminário presencial:

Data da apresentações:

- Data: definidas em cada turma.

1) CRIAR PROGRAMA PARA O CÁLCULO
DA PRESSÃO DE VAPOR DE UM FLUIDO
UTILIZANDO A EQUAÇÃO DE ESTADO DE
PENG-ROBINSON

- Objetivo: Calcular a Pressão de Vapor (Pvap) para uma

substância pura em diversas temperaturas.
- Implementar o fluxograma da Fig. 7.5.1, pág. 306,
disponível no livro “Chemical, Biochemical, and
Engineering Thermodynamics, 4th, Sandler, 2006.
- Pode-se utilizar qualquer programa de planilhas como
Excel, Origin, Matlab, etc.
EQUAÇÕES:

𝑅𝑇 𝑎(𝑇)
𝑃= − (6.4-2)
V−𝑏 𝑉 𝑉+𝑏 +𝑏(𝑉−𝑏)

Com:
𝑅 2 𝑇𝑐2
𝑎 𝑇 = 0.45724 α 𝑇 (6.7-1)
𝑃𝑐

𝑅𝑇𝐶
𝑏 = 0.07780 (6.7-2)
𝑃𝑐

𝑇
𝛼 =1+κ 1− (6.7-3)
𝑇𝐶

Κ = 0.37464 + 1.54226 ∙ 𝜔 − 0.26992 ∙ 𝜔2 (6.7-4)

𝑍 3 − 1 − 𝐵 𝑍 2 + 𝐴 − 3𝐵2 − 2𝐵 𝑍 − 𝐴𝐵 − 𝐵2 − 𝐵3 = 0 (6.7-5)

𝑃𝑏 𝑎𝑃
Onde: 𝐵= e 𝐴=
𝑅𝑇 𝑅2 𝑇 2
EQUAÇÕES:

𝑉
𝑓𝑉 𝑉 𝑉
𝑏𝑃 𝑎 𝑍 + 1 + 2 𝑏𝑃Τ𝑅𝑇
𝑙𝑛 = 𝑍 − 1 − 𝑙𝑛 𝑍 − − 𝑙𝑛
𝑃 𝑅𝑇 2 2𝑏𝑅𝑇 𝑍 𝑉 + 1 − 2 𝑏𝑃Τ𝑅𝑇

𝐴 𝑍𝑉 + 1 + 2 𝐵
= 𝑍 𝑉 − 1 − 𝑙𝑛 𝑍 𝑉 − 𝐵 − 𝑙𝑛
2 2𝐵 𝑍𝑉 + 1 − 2 𝐵
(7.4-14a)

𝐿
𝑓𝐿 𝐿 𝐿
𝑏𝑃 𝑎 𝑍 + 1 + 2 𝑏𝑃Τ𝑅𝑇
𝑙𝑛 = 𝑍 − 1 − 𝑙𝑛 𝑍 − − 𝑙𝑛
𝑃 𝑅𝑇 2 2𝑏𝑅𝑇 𝑍𝐿 + 1 − 2 𝑏𝑃Τ𝑅𝑇

𝐴 𝑍𝐿 + 1 + 2 𝐵
= 𝑍𝐿 − 1 − 𝑙𝑛 𝑍𝐿 −𝐵 − 𝑙𝑛
2 2𝐵 𝑍𝐿 + 1 − 2 𝐵
(7.4-14b)
 Entrar com 𝑇𝐶 , 𝑃𝐶 e ω
Fig. 7.5-1 – Fluxograma de um programa
de computador para o cálculo da Pvap de
Entrar com T e chutar valor de P
um fluido usando a equação de estado de
Peng-Robinson. Pág. 306. Calcular a e b utilizando as
Eqs. 6.7-1, 6.7-2 e 6.7-3

Calcular A e B utilizando
𝑎𝑃 𝑃𝑏
𝐴= 2 2 𝑒𝐵=
𝑅 𝑇 𝑅𝑇

Resolver a Eq. 6.7-5 para 𝑍 𝐿 𝑒 𝑍 𝑉

Calcular 𝑓 𝑉 substituindo
𝑍 𝑉 na Eq. 7.4-14a

Fonte: Calcular 𝑓 𝐿 substituindo

𝑍 𝐿 na Eq. 7.4-14b
Chemical, Biochemical, and Engineering
Thermodynamics – 4th ed./2006. SANDLER, Stanley I. 𝑓𝐿
New York: J. Wiley, 1999. Número de chamada: 66.021- Se 𝑉 − 1 < 0.0001
𝑓
97 S217c 4. ed.
Não Sim
Chapter 7: Equilibrium and Stability in One-Component
Systems 𝑓𝐿 Imprimir a Pressão de
𝑃=𝑃 𝑉
𝑓 Equilíbrio (vapor)
Section 7.5: The Calculation of Pure Fluid-Phase
Equilibrium: The computation of Vapor Pressure from a Sair ou repetir o cálculo
Equation of State; pág. 305. para outra temperatura
PROGRAMA PARA O CÁLCULO DA PRESSÃO DE VAPOR DE UM FLUIDO
UTILIZANDO A EQUAÇÃO DE ESTADO DE PENG-ROBINSON
2) CÁLCULO DE DADOS PVT E
PROPRIEDADES TERMODINÂMICAS
DE SUBSTÂNCIAS PURAS
UTILIZANDO PROGRAMAS

- Utilizado para substâncias puras;

- Calcula dados PVT e prop. Termodinâmicas de
substâncias puras;
- Cria gráficos e tabelas com os valores obtidos.

Exemplos de programas : XSEOS, GPEC, PR1, Thermo

Solver,
Endereços eletrônicos de programas para cálculos de propriedades
termodinâmicas e equilíbrio de fases:

PHASE EQUILIBRIA (PE 2000):

https://www.tuhh.de/v8/links/pe-2000/download.html

THERMO SOLVER:
http://bcs.wiley.com/he-
bcs/Books?action=resource&itemId=0471385867&bcsId=1793&resourceId=2674

XSEOS:
http://people.qatar.tamu.edu/marcelo.castier/index.htm

GPEC (Deve-se solicitar aos desenvolvedores/MOODLE):

E-mail: gpec@phasety.com
https://www.youtube.com/watch?v=b1VD4_2EpwA

PR1/UNIFAC/VLMU/CHEMEQ (Solicitar por e-mail/MOODLE):

E-mail: tmwalsh@udel.edu
http://www.che.udel.edu/thermo/basicprograms.htm
Endereços eletrônicos de programas para cálculos de propriedades
termodinâmicas e equilíbrio de fases:

TA – Thermodynamics Application:
http://mypages.iit.edu/~abbasian/thermo-app.html

G&P Engineering Software (Pago):

http://gpengineeringsoft.com/pages/products.html

SOPE:
http://www.eq.uc.pt/~abel/simul.htm

CHEMSEP:
http://www.chemsep.org/downloads/index.html
https://www.youtube.com/watch?v=mwH_wBQUSiM

STEAM TABLES ONLINE:

https://www.steamtablesonline.com/steam97web.aspx

NIST WEBBOOK:
https://webbook.nist.gov/chemistry/
Outros programas para cálculos de propriedades termodinâmicas e
equilíbrio de fases disponíveis no MOODLE:

GPEC e PR1: Mostrados anteriormente.

ELV 1.0: Calcula o ELV de misturas.

PROPITERM – LASEFI: Calcula dados PVT para substâncias puras usando a

equação de Peng-Robinson.

UNIFAC 3.0 e UNIFAC/LASEFI: Calculam os coeficientes de atividade em

misturas.

TERMOFIS 1.1: Calcula propriedades termodinâmicas de substâncias puras.

PCO2: Calcula propriedades termodinâmicas do CO2.

OCTOPUS 1.0: Suplemento para uso no Excel que calcula propriedades

termodinâmicas como Pvap, densidade, etc.
APLICATIVOS PARA SMARTPHONE

DWSIM CALCULATOR:
http://dwsim.inforside.com.br/wiki/index.php?title=Calculator

https://play.google.com/store/apps/details?id=com.danielmedeiros.dwsim_calculator&hl=enc

PHASEDIAGRAM
https://play.google.com/store/apps/details?id=com.andrew.liashuk.phasediagram&hl=pt_BR

STEAM PROPERTY - LITE

https://play.google.com/store/apps/details?id=com.ADLS.steamproperty&hl=pt_BR

STEAM CALCULATOR
https://play.google.com/store/apps/details?id=ru.lextop.steamcalculator&hl=pt_BR

INTERNATIONAL STEAM TABLES

https://play.google.com/store/apps/details?id=mk.internationalsteamtables.app&hl=pt_BR

STEAM TABLES
https://play.google.com/store/apps/details?id=chemical.applications.steam.tables&hl=pt_BR
APLICATIVOS PARA SMARTPHONE

VAPOUR PRESSURES
https://play.google.com/store/apps/details?id=chemical.applications.pressures.vapour&hl=pt_
BR

STEAMTABLE
https://play.google.com/store/apps/details?id=in.augsys.steamtable&hl=pt_BR

STEAM TABLE
https://play.google.com/store/apps/details?id=com.aswdc_steamtable&hl=pt_BR

WATER VAPOUR PRESSURE CALCULATOR

https://play.google.com/store/apps/details?id=com.kuldemontorenes.pressurecalculator&hl=pt
_BR

CALC PSYCHROMÉTRICO
https://play.google.com/store/apps/details?id=com.deltatsoftware.pyschrometriccalc&hl=pt_B
R

PSYCHROMETRIC CALCULATOR
https://play.google.com/store/apps/details?id=com.swastic.psychrometricairpropertycalculator
&hl=pt_BR
PESQUISA DE PROPRIEDADES TERMODINÂMICAS PARA
SUBSTÂNCIAS PURAS: NIST WEBBOOK

https://webbook.nist.gov/chemistry/

https://webbook.nist.gov/chemistry/name-ser/
EXEMPLO - CÁLCULO :

Utilizar os programas a seguir para calcular a Pressão de

Saturação ou Vapor (Psat ou Pvap) para o Etanol a 298,15 K (25
°C);

1) PHASE EQUILIBRIA (PE 2000)

2) THERMO SOLVER
3) PR1
4) TA – Thermodynamics Application
5) PROPITERM – LASEFI
6) G&P Engineering Software (14 days free trial/pago)

Calcular propriedades termodinâmicas para o CO2:

1) PCO2
2) NIST WEBBOOK
EXEMPLO:

Utilizar os programas a seguir para determinar dados PVT e

construir um diagrama PV (com várias isotermas) para uma
substância pura;

1) PHASE EQUILIBRIA (PE 2000)

https://www.tuhh.de/v8/links/pe-2000/download.html

2) XSEOS
http://people.qatar.tamu.edu/marcelo.castier/index.htm

3) PROPITERM – LASEFI
Disponível no Moodle
1) DIAGRAMA PVT - PHASE EQUILIBRIA (PE 2000)
Exemplo: Amonia (Nitreto de Hidrogênio, NH3)

Propriedades NH3
M (g/mol) 17,031
Tc (K) 405,7
Pc (bar) 112,8
ω 0,253
Tn (K) 239,7
R (cm^3*bar/mol*K) 83,14
 Propriedades NH3
2) PropiTerm - LASEFI M (g/mol)
Tc (K)
17,031
405,7
Pc (bar) 112,8
ω 0,253
Exemplo: Amonia (Nitreto de Hidrogênio, NH3) Tn (K) 239,7
R (cm^3*bar/mol*K) 83,14
A) Pvap foi calculada com o PE2000 (Vapor Pressure
Curve). b) Entra com dados usando o PropiTerm:
CURVA DE SATURAÇÃO
T (K) Pvap (bar) T (K) Pvap (bar)
195,4 0,059 305,4 12,380
200,4 0,087 310,4 14,300
205,4 0,126 315,4 16,438
210,4 0,179 320,4 18,810
215,4 0,249 325,4 21,430
220,4 0,341 330,4 24,317
225,4 0,459 335,4 27,486
230,4 0,609 340,4 30,955
235,4 0,797 345,4 34,742
240,4 1,030 350,4 38,865 c) Calcula as isotermas:
245,4 1,315 355,4 43,343
250,4 1,661 360,4 48,196
255,4 2,076 365,4 53,442
260,4 2,569 370,4 59,105
265,4 3,151 375,4 65,203
270,4 3,832 380,4 71,759
275,4 4,623 385,4 78,795
280,4 5,537 390,4 86,337
285,4 6,584 395,4 94,407
290,4 7,779 400,4 103,033
295,4 9,134 403,4 107,808
300,4 10,663 405,4 112,197
3) XSEOS/OCTOPUS
Exemplo: Amonia
(Nitreto de
Hidrogênio, NH3)

6 Diagrama PV amônia Líquido saturado

5 Vapor saturado
4 Ponto crítico
3 T=265.4 líq
Região líquido
ln P (bar)

2 subresfriado T=265.4 vap

1

0 Região líquido-vapor
Região vapor
-1 superaquesido

-2

-3
0 2 4 6 8 10 12 14

ln Volume (mol/cm^3)
3) CÁLCULO DO ELV DE MISTURAS
BINÁRIAS UTILIZANDO SOFTWARE
OPEN SOURCE

PHASE EQUILIBRIA (PE 2000):

https://www.tuhh.de/v8/links/pe-2000/download.html

THERMO SOLVER:
http://bcs.wiley.com/he-
bcs/Books?action=resource&itemId=0471385867&bcsId=1793&resourceId=2674

XSEOS:
http://people.qatar.tamu.edu/marcelo.castier/index.htm

GPEC (Deve-se solicitar aos desenvolvedores/MOODLE):

E-mail: gpec@phasety.com
https://www.youtube.com/watch?v=b1VD4_2EpwA
EXEMPLO - CÁLCULO de PHI e GAMA:

Utilizar os programas a seguir para calcular a o coeficiente de

atividade (γ) e o coeficiente de fugacidade (φ) para uma mistura
de tolueno (1) + clorobenzeno (2) a 343 K e x1 = 0,4.

1) THERMO SOLVER
2) UNIFAC Calculator
3) XSEOS

Tolueno (C7H8) Clorobenzeno (C6H5Cl)

M (g/mol) 92,14 M (g/mol) 112,56
Tc (K) 591,8 Tc (K) 632,4
Pc (bar) 41,06 Pc (bar) 45,20
ω 0,262 ω 0,250
Tn (K) 383,8 Tn (K) 404,9
 ELV DE MISTURAS BINÁRIAS

1) Procurar artigos com dados experimentais de ELV utilizando

bases de dados

2) Verificar o tipo de cálculo e abordagem (gama-phi, phi-phi, ideal)

foram assumidas no artigo

3) Utilizar os softwares listados abaixo e efetuar os cálculos de

equilíbrio líquido-vapor (ELV) de misturas binárias:

PHASE EQUILIBRIA
THERMO SOLVER
GPEC
XSEOS
ELV 1.0
BASE DE DADOS PARA PESQUISA DE ARTIGOS DE EQUILÍBRIO
DE FASES:

Periódicos CAPES:
http://www.periodicos.capes.gov.br.ez46.periodicos.capes.gov.br/

Web of Science: www.isiknowledge.com/

American Chemical Society: http://pubs.acs.org/

Science Direct: http://www.sciencedirect.com/

Scielo: http://www.scielo.org/php/index.php

NIST: http://webbook.nist.gov/chemistry/
Fig. 1. P–x1–y1 diagram for tetrahydrofuran (1) + 1-chloropentane (2): (□,■) experimental data
at 298.15 K; (●,○) experimental data at 313.15 K; (▲) experimental data at 328.15 K; (—)
Wilson equation
 Tetrahuidrofurano (C4H8O)
PHASE EQUILIBRIA M (g/mol) 72.10
Tc (K) 540.2
Pc (bar) 51.9
ω 0.217
Tn (K) 338.7

1-cloropentano (C5H11Cl)
M (g/mol) 106.59
Tc (K) 568.5
Pc (bar) 33.22*
Ω 0.331*
Tn (K) 320.1
* Estimados
THERMO SOLVER
 TUTORIAIS PARA USO DO XSEOS:

Diagrama PVT: https://www.youtube.com/watch?v=ci3LIBPOkxI

Canal no YouTube Prof. André Zuber (UTFPR):

Instalando o XSEOS: https://www.youtube.com/watch?v=AkRJ5Ic8PTQ

Calcular a fugacidade: https://www.youtube.com/watch?v=jCLg6zPVYUA

Canal YouThermo do Prof. Marcelo Castier:

https://www.youtube.com/channel/UCYNDgYUJKT6JcSzAoMitTEA/videos

Resolução de Problemas de ELV usando o XSEOS (Parte 1, 2 e 3):

https://www.youtube.com/watch?v=DUO2VhpBmrc
https://www.youtube.com/watch?v=VFB-Zxm4Jdw
https://www.youtube.com/watch?v=gPmLNrP5ccQ
 REFERÊNCIAS:

MARTÍN, A.; BERMEJO, M. D.; MATO, F. A.; COCERO, M. J.

Teaching advanced equations of state in applied thermodynamics
courses using open source programs, v. 6, p. e114-e121, 2011.
https://www.sciencedirect.com/science/article/pii/S17497728110000
91

CASTIER, M.; AMER, M. M. XSEOS: An evolving tool for teaching

chemical engineering thermodynamics, v. 6, p. e62-e70, 2011.
https://www.sciencedirect.com/science/article/pii/S1749772811000030