Exercício 1 - Responda as questões teóricas a seguir embasando seus argumentos com

as teorias apresentadas no curso:

a) Em que situações a Lei de Raoult para o equilíbrio líquido-vapor pode ser aplicável?
Justifique.
A Lei de Raoult pode ser aplicada se as seguintes considerações forem
satisfeitas:
● Os componentes na fase vapor se comportarem como gases ideais, ou seja, a
baixas pressões f vi =1;
● Quando a pressão do sistema for aproximadamente a pressão de saturação o
fator de Poynting pode ser considerado igual a 1.
● Quando a pressão de saturação é baixa pode considerar que f sat
i =1
● Quando as interações intermoleculares na fase líquida forem parecidas e
permitindo considerar o comportamento da fase líquida como ideal.
b) Quais as diferenças das abordagens 𝛾-𝜑 e 𝜑 -𝜑? Quando cada estratégia costuma ser
utilizada?
● 𝛾-𝜑: pressão baixa; pela EoS(equação de estado) calcula-se o coeficiente de
fugacidade (𝜑) que representa a não idealidade na fase vapor e pela energia de
Gibbs de excesso calcula-se o coeficiente de atividade (𝛾) que representa a não
idealidade na fase líquida e é capaz de descrever o comportamento do sistema
independente do valor de pressão.
● 𝜑 -𝜑: pressão alta; deve haver uma EoS (equação de estado) válida para toda a
faixa de pressão que represente as fugacidades tanto na fase líquida quanto na de
vapor pelo coeficiente de fugacidade (𝜑) e as propriedades volumétricas em
função da temperatura, pressão e composições.
c) Sobre os azeótropos presentes em diagramas de fases, responda:
i) O que são azeótropos?
O nome azeótropo é proveniente do grego e significa ferver sem alteração. Eles são os
pontos de extremos, tanto máximo quanto mínimo, nos diagramas Pxy. Nesses pontos a
fração molar na fase liquida é igual fração molar na fase vapor para todos os
componentes.

ii) Cite uma situação prática em que o azeótropo é desfavorável e uma situação em que
ele pode ser útil, justificando suas escolhas.
Desfavorável: a mistura não pode ser separada por destilação simples
Favorável: cálculo dos coeficientes de atividade (𝛾):
y i∗P
γ i=
x i∗Psat
i

Para azeótropos: x i= y i, assim a equação de γ i simplifica-se para:

 P
γ i=
Psat
i

iii) Em uma mistura binária entre as substâncias A e B, o diagrama de ELV Pxy

apresenta um azeótropo mínimo. O que pode inferir sobre as interações intermoleculares
a-a, b-b e a-b?
As interações entre moléculas distintas (a-b) são mais fortes que as interações entre
moléculas iguais (a-a e b-b).
d) Em que situação prática o estado referência de Henry é comum de ser utilizado?
Mostre as equações que descrevem o equilíbrio nessa situação para ambas as
substâncias envolvidas, considerando uma mistura binária e o sistema mais
simplificado possível.
O estado referência de Henry é comum de ser utilizado para o equilíbrio líquido vapor
(ELV) de um gás dissolvido em um líquido.
Sistema mais simplificado possível:
● sistema binário
● O soluto de comportando como gás ideal
● Solubilidade do gás no líquido baixa
● O solvente se comporta como solução ideal no estado de referência de Lewis-
Randall
Assim, a Lei de Henry descreve as interações entre os compostos diferentes (a-b) pela
equação do equilíbrio.
Lei de Henry:
Pi=H i∗xi (1)
Onde:
Pi= y i∗P (2)
Então, temos que:
y i∗P=H i∗x i (3)
Como consideramos que a solubilidade do gás o líquido é baixa e, portanto, que o
solvente se comporta como solução ideal no estado de referência de Lewis-Randall,
podemos escrever para o líquido que:
Hl
γ ∞l = (4)
fl
f l=Psat
l (5)
γ ∞l =1 (6)
H l=Plsat (7)
Somando as duas pressões:

P=Pl + P g (8)
Ficamos que:
P=H g∗x g+ P satl ∗x l (9)
Isolando y i pela Equação 2, e substituindo Pi pela Equação 1 encontramos que:
H g∗x g
y g= (10)
H g∗x g+ P sat
l ∗x l
P sat
l ∗x l
y l= (11)
H g∗x g + Psat
l ∗x l
Onde H g é obtido experimentalmente. Escolhendo uma abordagem para o cálculo do
coeficiente de atividade, conseguimos descrever o equiíbrio nessa situação para ambos
os componentes pela Equação 12.
y i∗φvi ∗P=x i∗γ iHenry ' s∗H i

Exercício 2 - Julgue os itens abaixo como verdadeiros ou falsos, justificando sua

resposta.
a) A fase extrato, na extração líquido-líquido, corresponde à fase rica no solvente.
Falso. A fase extrato, na extração líquido-líquido, corresponde a fase rica em
soluto, ou seja, no componente de interesse da extração.
b) Na extração líquido-líquido as linhas de amarração são sempre retas paralelas ao eixo
das abscissas (eixo que corresponde à base do diagrama ternário).
Falso. As linhas de amarração no diagrama ternário não são necessariamente paralelas
ao eixo das abscissas.
c) O parâmetro seletividade indica a viabilidade de um processo de extração
líquido-líquido. Nesse sentido, valores menores que a unidade indicam que o processo
é aplicável.
Falso. Valores menores que uma unidade indicam que a concentração de soluto na fase
rafinado, ou seja, a de descarte/ sobra que é pobre no soluto, é maior que a concentração
de soluto na fase extrato, ou seja, a de interesse e que possui o soluto, indicando assim
que p processo não é aplicável.
d) É desejável que a fase rafinado seja rica no solvente.

Figura 1 – Esquema de separação.

Fonte: COSTA, Mariana Conceição. TERMODINÂMICA DE EQUILÍBRIO:

Equilíbrio Líquido-Líquido. Laboratório de Equilíbrio de Fases, Campinas, 2020.

Falso. É desejável que a fase extrato seja rica em solvente e rafinada pobre em soluto,
assim como ilustra a Figura 1.
e) O modelo de Wilson e o modelo Uniquac descrevem muito bem o equilíbrio
líquido-líquido.
Falso. A parte “o modelo Uniquac descrevem muito bem o equilíbrio líquido-líquido”
está correta, entretanto, o modelo de Wilson não consegue descrever o ELL, já que a
segunda derivada da energia de Gibbs de excesso nunca será menor ou igual a zero.
f) A temperatura consoluta superior é a menor temperatura a partir da qual acima dela a
fase líquida não se separa em duas fases.
Falso. Existem sistemas como o representado pela Figura 2 em que é observado além da
temperatura consoluta superior, a temperatura consoluta inferior. Abaixo da T ci a fase
líquida também não se separa em duas fases.
Figura 2 – Exemplo de sistema com da temperatura consoluta superior e inferior.

.
Fonte: Amélia M.G. do Val; Rosana Z. Domingues; Tullo Matencio- Unidade II,
Diagrama de fases. Química, modalidade a distância, UFMG, 2013.
g) O sistema é instável quando se encontra entre a curva binodal e a spinodal.
Falso. O sistema é metaestável quando se encontra entre a curva binodal e a spinodal e
instável dentro da curva spinodal.
h) Em um mesmo sistema é possível existir temperatura consoluta inferior e superior.
Verdadeiro. A Figura 1 ilustra um exemplo de sistema onde há Tci e Tcs.

Exercício 3 - Considere que o sistema metanol (1)/acetato de metila (2) tem suas
pressões de vapor dadas pelas equações de Antoine a seguir:

a) Considerando a fase líquida e vapor como ideal, determine a composição da fase

líquida e a pressão, quando a temperatura é 318,15 K e y1 = 0,6.
Nesse caso é possível aplicar a Lei de Raoult e a Lei de Dalton para o equilíbrio líquido-
vapor:
Pi= y i∗P (1)
sat
x i∗Pi
y i= (2)
P
P=P2sat + y 1∗P∗¿) (3)
Sabendo que y 1=0,6 e que T=318,15 K, é possível encontrar P pela Equação 3 e isolar
x i pela Equação 2:
y i∗P
x i= sat (4)
Pi
Encontrando o valor de Psat
i :
16,59158∗3643,31
Psat
1 =exp exp (318,15−33,424 )=44,51 kPa (5)
14,25326∗2665,54
Psat
2 =exp exp (318,15−53,424 )=65,64 kPa (6)
Encontrando o valor de P:
65,64 kPa
P−0,6∗P∗ 1− ( 44,51 kPa )=65,64 kPa
P=51,09 kPa
Substituindo os valores encontrados na Equação 4, encontramos x 1 e x 2:
0,6∗51,09
x 1= =0,69
44,51
0,4∗51,09
x 2= =0,31
65,64

b) Agora considere que a fase líquida não é ideal e os coeficientes de atividade são
fornecidos pela equação abaixo. Qual a pressão e a composição da fase líquida?
Discuta brevemente sobre o algoritmo utilizado.

Sabemos que:
y 1=0,6

y 2=0,4
Primeiro devemos encontrar a relação entre os coeficientes de atividade e as
composições na fase líquida para ambos os componentes pelas equações dadas no
exercício. Em seguida, encontramos a relação entre P, y 1 e y 2 pela Equação 3:
A' =2,771−0,00523∗318,15=1,107

γ 1=exp exp ( 1,107∗x 21 )

(1)
γ 2=exp exp ( 1,107∗x 22 )
(2)
y ∗P
γ i= i sat (3)
x i∗Pi
Substituindo a Equação 1 e 2 na Equação 3, chegamos em:
0,6∗P
exp exp ( 1,107∗x 22 )= (4)
x 1∗P1sat
 0,4∗P
exp exp ( 1,107∗x 21 )= (5)
x 2∗P2sat

Sabemos que a relação entre x 1 e x2 é:

x 1+ x2=1 (6)
Podemos isolar a pressão P nas Equações 4 e 5 e calcular o valor da composição
na fase líquida para um componente para depois pela Equação 6 obter a composição do
outro componente.

exp exp ( 1,107∗x 21 )∗x2∗Psat 2 sat

2 ∗0,6= ( 1,107∗x 2 )∗x 1∗P1

exp exp ( 1,107∗(1−x 2)2 )∗x 2∗Psat 2 sat

2 ∗0,6=( 1,107∗x2 )∗(1−x 2)∗P 1

Encontramos a partir dessa equação que x 2=0,18, portanto x 1=0,82. Substituindo esses
valores tanto na equação 4 quanto na Equação 5 encontramos o valor de P.
P=62,89 kPa
Exercício 4 - Sabe-se que o sistema etanol (a)-benzeno (b), à 25 ºC, apresenta um
azeótropo na pressão de 122,3 torr com composição x a = 0,28.
a) A partir dessa informação, determine os parâmetros da equação de g E de van
Larr.

Sabemos que:
x a=0,28

x b=0,72

Precisamos encontrar agora o valor dos coeficientes de atividade para os

componentes a e b.
Inicialmente encontramos a pressão do sistema em unidades mais usuais.
P=122,3 torr=16,31 kPa=0,1631 ¯¿
Assim, temos que:
y ∗P
γ i= i sat (1)
x i∗Pi
Tratando-se de um azeótropo sabemos que nesse ponto y i=x i e, por isso,
podemos simplificar a Equação acima para:
P
γ i= sat (2)
Pi
Agora, precisamos encontrar as pressões de saturação. Para isso, usaremos a
Equação de Antoine e as constantes fornecidas pelo Apêndice A.1 do Koretsky.
 Etanol:
A B C

Benzeno:
A B C

Dessa forma, as pressões de saturação calculadas foram:

Psat
a =7,88 kPa

Psat
b =12,69 kPa

Substituindo esses valores na equação 2, temos que os coeficientes de atividade

são:
γ a=2,07

γ b=1,29

Substituindo os valores de coeficiente de atividade e de composição na fase

líquida na equação dada no exercício chegamos em um sistema 2x2, onde as
incógnitas são A e B. Utilizando a ferramenta solver do Excel encontramos que:
A= 6 417,6
B=2 812,1
b) Determine a composição da fase líquida e a pressão quando o sistema com y a
=0,75 estiver em equilíbrio à 25 ºC.
No equilíbrio temos que:
y i∗P
γ i= (*)
x i∗Psat
i
Sabemos que:
y a=0,75

y b=0,25

x b=1−x a

γ a=2,07

γ b=1,29

Assim, chegamos em duas Equações com duas incógnitas P e x a . Utilizando a

ferramenta solver do Excel encontramos que:
x a=0,93

P=9,88 kPa
Sabendo que x b=1−x a=0,07
Exercício 5 - Considere uma mistura líquida parcialmente miscível formada por hexano
(1) e acetona (2) a 15 ºC e 300 bar. Sabe-se que a fase α contém um total de 20 mols e
x α1 = 0,2, enquanto a fase β contém um total de 10 mols e x β1 = 0,8. Além disso, tem-se as
seguintes informações:

a) A partir das informações dadas, as interações entre as moléculas iguais são mais
fortes ou mais fracas que as interações entre as moléculas distintas? Justifique.
As interações entre as moléculas iguais são mais fortes que as interações entre
moléculas distintas já que ocorre a separação em 2 fases.

b) Calcule o valor da fugacidade do hexano puro, utilizando a referência de Lewis

Randall.
Sabendo que a pressão de saturação do hexano puro é bem menor que a pressão do
sistema podemos considerar que o coeficiente de fugacidade na saturação é igual a 1.
Podemos considerar ainda que a variação do volume molar não varia consideravelmente
com a pressão. Assim, utilizamos a seguinte equação para calcular a fugacidade do
hexano puro para o estado de referência de Lewis Randall :

l sat sat v li sat

f i=φ i ∗Pi ∗exp ⁡( ∗( P−Pi ) ) (1)
RT

φ sat
i =1

cm3
l
f 1=0,0127 MPa∗exp ⁡(
130,5 ( )
mol
∗(30 MPa−0,0127 MPa ) )
8,314 ( cm3∗MPa∗K −1∗mol−1 )∗288,15 K

f l1=0,065 MPa=65 kPa

c) Determine o parâmetro A da equação de Margules com um parâmetro para esse

sistema, evidenciando as simplificações necessárias.
Podemos escrever a seguinte equação para o ELL:

x αi ∗γ αi =x iβ∗γ iβ (1)

A equação de Margules de um parâmetro é dada por:

A∗x 2i =RT∗ln γ j (2)

Podemos isolar o coeficiente de atividade da Equação 2:

A∗x 2i
γ j =exp ⁡( ) (3)
RT
Substituindo a Equação 3 na Equação 1 temos:

x α1∗exp ⁡¿ (4)

x α2∗exp ⁡¿ (5)

Aplicando o logaritmo neperiano em ambos os lados das Equações 4 e 5, chegamos nas

Equações 6 e 7, respectivamente:

ln ⁡( x α1 )∗¿

x α1
ln ⁡( β )=A∗¿ ¿ (6)
x1

x α2
ln ⁡( β
)=A∗¿ ¿ (7)
x2

Do enunciado: x α1 = 0,2 e x β1 = 0,8. Dessa forma, inferimos que x α2 = 0,8 e x β2 = 0,2.

Isolando o A, chegamos em:

x α2
( ) x β2
∗RT

¿¿
A=5535,19 J /mol
d) Determine a temperatura que o sistema precisa ser elevado para que ele apresente
apenas uma fase, enunciando as suas hipóteses.

Figura 3 – Curva binodal e Espinodal.

Fonte: COSTA, Mariana Conceição. TERMODINÂMICA DE EQUILÍBRIO:

Equilíbrio Líquido-Líquido. Laboratório de Equilíbrio de Fases, Campinas, 2020.
Onde as Equações 5, 6 e 7 são:
Assumindo que a constante A é independente da temperatura, podemos escrever que o
sistema precisa ser elevado a temperatura dada pela equação (*):
RT
=2∗A (*)
x a∗x b
Isolando T temos que:
2∗A∗( x a∗x b)
T= (**)
R
Sabemos pelo enunciado que: a fase α contém um total de 20 mols e x α1 = 0,2, enquanto a
fase β contém um total de 10 mols e x β1 = 0,8.
Dessa forma, para que o sistema apresente apenas uma fase, a quantidade de
matéria total deverá ser de 30 mols (soma da fase α com a fase β). Precisamos
saber agora a concentração de cada componente. Para isso calcularemos para o
componente 1 e usaremos a relação x b=1−x apara encontrar a concentração do
componente 2.
0,2∗20+ 0,8∗10
x 1= =0,4
30
x 2=1−0,4=0,6

Substituindo na Equação (**) temos:

J
2∗5535,19( )∗(0,4∗0,6)
mol
T= =319,6 K
J
8,314( )
K∗mol
e) Determine a constante de Henry para o hexano a 15 ºC e 300 bar.
A Lei de Henry aplica-se para a condição de diluição infinita:
∞ H1
γ1 =
f1

A∗x 2i
A f 1 foi obtida no item a) e o coeficiente de atividade é dado por γ j =exp exp ( RT )
Assim:
J
( mol )

( )
5535,19 ∗1
∞
γ =exp exp
1 =10,08
J
8,314 (
K∗mol )
∗288,15 K

γ ∞1 ∗f 1 =H 1

H 1=10,08∗0,065 MPa=655,2 kPa

Exercício 6 - Com base no diagrama binário formado por metanol + água na pressão de
1atm confirme se as seguintes informações estão corretas:

a) Existem 3 fases homogêneas neste diagrama formadas pela mistura de ambos os

componentes do sistema.
Verdadeiro.
b) Existem 3 regiões de equilíbrio sólido-líquido e uma de equilíbrio líquido-líquido
neste diagrama de fases.
Falso. Existem 3 regiões de equilíbrio sólido-líquido mas não existe uma de equilíbrio
líquido-líquido neste diagrama de fases, se tivesse teria que ter uma curva binodal.
c) Para a composição de 30 % de metanol na temperatura de -80 °C a fase sólida é
composta só por água (gelo) e a fase líquida tem 52% de água.
Verdadeiro.
d) A 90 °C coexistem em equilíbrio uma fase vapor cuja
composição é y CH3OH ≈ 0,65 e uma fase líquida cuja
composição y H2O ≈ 0,70.
Falso.
e) A composição do ponto eutético desse sistema é 75% de água e do ponto peritético
42% de água.
Falso. A composição do ponto eutético desse sistema é 75% de metanol e do ponto
peritético, está correto e é de 42% de água.
Exercício 7 - Uma corrente líquida comprimida contendo uma mistura equimolar de n-
butano e isobutano flui para dentro de uma unidade de flash a uma vazão F. No estado
estacionário, 40 % da vazão de alimentação é vaporizada e sai do tanque como uma
corrente de vapor a uma vazão V. O resto sai como líquido a uma vazão L. Se a pressão
do tanque flash é 1 bar, qual a temperatura requerida para o processo? Qual a
composição da fase líquida e da fase vapor? Se fizer alguma consideração, deixe-a
expressa.
Chamando o n-butano de 1 e o isobutano de 2, temos:
-Balanço de massa em estado estacionário:
1. Por componente:
x 1 ,f ∗F−x 1 ,v∗V −x 1 , L∗L=0 (1)
x 2 ,f ∗F−x 2 ,v ∗V −x 2 , L∗L=0 (2)
Pelo enunciado sabemos que:
x 1 ,f =x 2 ,f =0,5

V =0,4∗F
2. Global:
F=L+V
L=0,6∗F
Substituindo esses valores nas equações 1 e 2, temos:
0,5∗F−x 1, v∗0,4∗F−x 1 , L∗0,6∗F=0

0,5∗F−x 2, v∗0,4∗F−x 2 , L∗0,6∗F=0

Podemos dividir os dois lados das duas equações por F para obtermos:
0,5−x1 , v∗0,4−x 1 , L∗0,6=0 (3)
0,5−x 2, v∗0,4−x 2 , L∗0,6=0 (4)
Sabemos também que:
x 1 ,v + x 2 ,v =1

x 1 , L + x 2 ,L =1

Assim, precisamos de uma relação entre a fase vapor e a fase líquida para poder resolver
o sistema de equações.

Consideraremos que a solução apresenta um comportamento ideal na fase líquida já que

o n-butano e o isobutano são compostos semelhantes. Consideraremos também que na
fase vapor o comportamento desses compostos é de gás ideal já que a pressão do
sistema é baixa. Dessa forma, a Lei de Raoult pode ser aplicada.

x L Psat
i
x v=
P
Por isso, precisamos calcular a pressão de saturação pela equação de Antoine para
ambos os componentes com as constantes fornecidas pelo Apêndice A.1 do Koretsky.

n-butano
A B C

iso-butano
A B C
 B
ln ln P sat =A−
T+C

¿ A−ln ln Psat

¿ ( A−ln ln Psat )∗(T +C)

B
¿
( A−ln ln P sat )
B
¿ −C
( A−ln ln P sat )
Onde
x v∗P
Psat
i =
xL

Assim chegamos em um sistema 2x2 que por acaso é independente do sistema original
(Equações 3 e 4):
B
¿ −C
x ∗P
( A−ln ln 1 , v
x 1 ,L )
B
¿ −C
(1−x 1 , v )∗P
( A−ln ln
(1−x 1, L ) )
Igualando essas duas Equações e usando o solver do Excel, chegamos em:
x 1 ,v =0,48 x 1 , L =0,59

Assim, encontramos:
x 2 ,v =0,52 x 2 , L=0,41

¿ 267,2 K
Exercício 8 - A corrente de alimentação de um tanque flash contém uma mistura de 40
% de n-butano, 30 % de n-pentano e 30 % n-hexano. A temperatura do flash é de 290 K
e a pressão do flash é de 0,6 bar. Qual a razão entre a vazão de vapor e a vazão de
alimentação (V/F)? Qual a composição do líquido (x i ) e do vapor (y i )? Descreva
como você chegou na equação final, bem como a forma de resolvê-la.

Chamando o n-butano de 1, o n-pentano de 2 e o n-hexano de 3, temos:

-Balanço de massa em estado estacionário:
1. Por componente:
x 1 ,f ∗F−x 1 ,v∗V −x 1 , L∗L=0 (1)
x 2 ,f ∗F−x 2 ,v ∗V −x 2 , L∗L=0 (2)
x 3 ,f ∗F−x 3 , v∗V −x 3 , L∗L=0 (3)
Pelo enunciado sabemos que:
x 1 ,f =0,4

x 2 ,f =0,3

x 3 ,f =0,3

2. Global:
F=L+V
Substituindo esses valores nas equações 1 e 2, temos:
0,4∗F−x 1 ,v ∗V −x 1 , L∗L=0

0,3∗F−x 2, v∗V −x 2, L∗L=0

0,3∗F−x 3 ,v∗V −x 3 , L∗L=0

Dividindo ambos os lados das três equações por F chegamos em:

0,4−x 1 ,v ∗V /F− x1 , L∗L/ F=0

0,3−x 2, v∗V /F−x 2 , L∗L/F=0

0,3−x3 , v∗V /F−x 3 , L∗L/ F =0

Sabemos também que:

x 1 ,v + x 2 ,v + x 3 , v =1

x 1 , L + x 2 ,L + x 3 ,L =1

São necessárias mais equações para viabilizar o cálculo das composições nas correntes
L e V.
Faremos a mesma consideração realizada no exercício anterior, ou seja, fase líquida e
vapor comportando-se idealmente já que a pressão do sistema é baixa e os compostos
possuem naturezas semelhantes. Assim, a Lei de Raoult é aplicável:
sat x v∗P
Pi =
xL

Calculando as pressões de saturação por Antoine é possível chegar em três equações

relacionando a composição da fase líquida e da fase vapor.
n-butano (1)

n-pentano (2)

n-hexano (3)

2154,9
Psat ( 290−34,42
1 =exp exp 9,058− )=1,87123 ¯¿ 187,123 kPa
2477,07
=exp exp ( 9,2131−
290−39,94 )
Psat
2 =0,50017 ¯¿ 50,017 kPa

2697,55
=exp exp ( 9,2164−
290−48,78 )
Psat
3 =0,13995 ¯¿ 13,995 kPa

x 1 ,v ∗60
187,123=
x1 , L
x 2 , v∗60
50,017=
x2 , L
x 3 ,v ∗60
13,995=
x3 , L
x 1 ,v + x 2 ,v + x 3 , v =1

x 1 , L + x 2 ,L + x 3 ,L =1

0,4−x 1 ,v ∗V /F− x1 , L∗L/ F=0

0,3−x 2, v∗V /F−x 2 , L∗L/F=0

0,3−x3 , v∗V /F−x 3 , L∗L/ F =0

Agora possuímos 8 equações e 8 incógnitas possibilitando a resolução do sistema.

−187,123∗x1 , L + 0∗x 2 , L +0∗x 3 , L +60∗x 1 ,v +0∗x 2, v +0∗x3 , v =0

0∗x 1 ,L −50,017∗x 2 ,L +0∗x3 , L +0∗x1 , v + 60∗x 2 , v + 0∗x 3 , v =0

0∗x 1 ,L +0∗x2 , L −13,995∗x 3 , L + 0∗x 1 , v + 0∗x 2 , v +60∗x 3 ,v =0

0∗x 1 ,L +0∗x2 , L + 0∗x 3 , L + x 1 , v + x 2 , v + x 3 , v =1

x 1 , L + x 2 ,L + x 3 ,L +0∗x1 , v +0∗x2 , v + 0∗x 3 , v =1

L x 1 , v∗V
∗x 1 , L +0∗x 2 , L +0∗x 3 ,L + + 0∗x 2 , v + 0∗x 3 , v =0,4
F F
 L V
0∗x 1 ,L + ∗x 2 , L +0∗x 3 ,L +0∗x1 , v + ∗x 2 , v + 0∗x 3 , v =0,3
F F
L V
0∗x 1 ,L +0∗x2 , L + ∗x 3 ,L +0∗x1 , v + 0∗x 2 , v + ∗x 3 , v =0,3
F F
Utilizamos solver Excel para resolver esse problema e multiplicação de matrizes.
Criamos uma matriz 6x6 (matriz A) com as seis primeiras equações com os termos
dependentes e uma matriz 6x1(matriz b) com os termos independentes dessas equações.
Invertemos a matriz 6x6 e multiplicamos pela matriz dos termos independentes para
L
encontrar a composição. Em seguida, criamos duas células para variar os valores de e
F
V
e passamos os termos independentes para o lado esquerdo subtraindo os valores de
F
ambos os lados da igualdade. Assim, pedimos para o solver zerar a equação 7 com a
restrição de que a equação 8 também deveria ser zero. A metodologia e os resultados
estão apresentados na Figura 4.
Figura 4 – Resolução do sistema.

Resposta final:
V
=0,482
F
L
=0,528
F
x 1 , L=0,198

x 2 , L=0,326

x 3 , L=0,476

x 1 ,v =0,617

x 2 ,v =0,272
 x 3 ,v =0,111

OBS= Denominamos y ide x i ,v .

Exercício 9 - Calcule a composição dos pontos indicados no diagrama ternário.

Valores aproximados:
Fase α: rica em B
Fase β: rica em A
Ponto A: Global: - 60% de C
- 70% de B - 40% de B
- 10% de A - 20% de A Ponto E:
- 20% de C - 40% de C Global:
Ponto B:
Binodal: Ponto D: - 0% de B
- 40% de B Global: - 10% de A
- 11% de A - 20% de B - 90% de C
- 49% de C - 60% de A
Ponto C: - 20% de C Fase β
Fase β - 0% de B
- 20% de B Fase β - 0% de A
- 70% de A - 15% de B - 100% de C
- 10% de C - 75% de A Fase α
Fase α - 10% de C - 0% de B
- 45% de B Fase α - 90% de A
- 10% de A - 32% de B - 10% de C
- 45% de C - 8% de A

Exercício 10 - Um processo industrial tem por objetivo produzir um determinado

produto A.
Sabe-se que na corrente de saída do reator existe o produto A misturado com o excesso
de reagente B que não foi consumido. O produto A precisa ser separado desta corrente
para ser vendido comercialmente. Sabe-se que a melhor maneira de manter as
propriedades e a qualidade do produto A apropriadas é separá-lo através de um processo
de extração líquido-líquido. Na indústria em questão, existem dois solventes (S1 e S2)
que poderiam ser aplicados nesta extração. Baseado nos diagramas ternários de ELL de
cada sistema avalie a possibilidade de se separar o produto A do reagente B utilizando
um dos solventes disponíveis. Na sua análise, considere que se deseja obter o produto A
na maior pureza possível, e que na corrente de saída do reator existem apenas o produto
A e o reagente B
Deve-se analisar o coeficiente de partição para os dois casos:
Analisar o ponto nos diagramas em que a fase extrato possua a maior quantidade de
soluto, ou seja, do componente A.
Solvente 1:
Fase Rafinado:
-11% de A
-89% de S1
-0% de B
Fase Extrato:
-70% de A
-11% de S1
-19% de B
Solvente 2:
Fase Rafinado:
-67% de A
-12% de S1
-21% de B
Fase Extrato:
-52% de A
-1% de S1
-47% de B
Cálculo do coeficiente de partição do soluto A:

70
K 1= =6,36
11
52
K 2= =0,78
67
Percebe-se que o maior coeficiente de partição é o do solvente 1 e isso significa que há
mais soluto na fase extrato que na fase rafinada e dessa forma, ele deve ser o solvente
indicado para a realização dessa extração.
Exercício 11 - Bismuto (1) e cádmio (2) formam um eutético a 144 °C e x1 = 0,45. Eles
são completamente imiscíveis na fase sólida. Os pontos de fusão para o bismuto e
cádmio são 271°C e 321°C, respectivamente. As entalpias de fusão de seus elementos
são 210,46 kJ/mol e 26,1 kJ/mol, respectivamente. Estime os parâmetros do modelo de
Margules (2 parâmetros) baseado apenas nestes dados fornecidos. Faça as considerações
que julgar necessárias, desde que enunciadas na resolução feita por você.

Sabendo que o bismuto e o cádmio são completamente imiscíveis na fase sólida, temos
que:

f S1 =x1 y 1 f l1 (1)
f S2 =x2 y 2 f l2 (2)

Desconsiderando ∆ c SL
p , chegamos nas Equações 3 e 4:

∆ hfus ,Tm 1 1
ln ( γ a ) =
R [ −
T Tm ]
−ln ( x a ) (3)
∆h 1 1
ln ( γ b ) = fus ,Tm
R [ −
T Tm ]
−ln ( x b ) (4)

As expressões para os coeficientes de atividade pelo modelo de Margules (2

parâmetros) são dadas por:
Figura 1

Igualando as expressões da Figura 1 com as equações 3 e 4 e substituindo as variáveis

com os dados fornecidos pelo exercício
-Ponto eutético a 144 °C e x1 = 0,45
- Os pontos de fusão para o bismuto e cádmio são 271°C e 321°C, respectivamente. As
entalpias de fusão de seus elementos são 210,46 kJ/mol e 26,1 kJ/mol, respectivamente
Chegamos nas Equações 5 e 6.

( A +3 B )( 0,55 )2−4 B(0,55)3 210460 1 1

(8,314 )(417,15)
=
8,314 417,15 [
−
544,15 ]
−ln (0,45)
(5)
2 3
( A−3 B ) ( 0,45 ) + 4 B(0,45) 26100 1 1
(8,314 )(417,15)
=
[ −
8,314 417,15 594,15 ] −ln(0,55)
(6)

Usando o solver do Excel nós obtemos:

A=122374,57 J/mol
B=61449,08 J/mol

Exercício 12 - Para cada um dos pontos do diagrama de ELLV abaixo, especifique

quais fases existem em cada ponto da mistura e qual a composição de cada uma dessas
fases.

M1: vapor
Co=0,25 de X
Fação de X=1

M2: α - vapor

Vapor:
0,3 de x
α:
0,12 de x
Peso:
S 0,11
Wα= = =0,61
R +S 0,11+0,07
R 0,07
Wβ= = =0,39
R +S 0,11+0,07

M3: β vapor
Vapor:
0,6 de x
β:
0,95 de x
Peso:
S 3,1
Wα= = =0,98
R +S 3,1+ 0,05
R 0,05
Wβ= = =0,02
R +S 3,1+ 0,05

M4: ponto eutético (α líquido + β líquido + vapor)

Vapor – 0,4
x
Líquido α = 0,13
Líquido β =0,61
y
Líquido α = 0,87
Líquido β = 0,39
Peso:
S 0,2
Wα= = =0,44
R +S 0,2+ 0,25
R 0,25
Wβ= = =0,55
R +S 0,2+ 0,25

M5: α
Co=0,05 de X
Fação de X=1
M6: α + β (líquido)
Vapor:
0,73 de x
α:
0,03 de x
β:
0,89 de x
Peso:
S 0,1
Wα= = =0,11
R +S 0,1+ 0,8
R 0,8
Wβ= = =0,89
R +S 0,1+ 0,8