Vapor-Liquid Equilibrium Diagram for Non-Ideal Mixture

​
diagram
P-x-y
T-x-y
temperature (°C)
72
​
This Demonstration shows
P-x-y
and
T-x-y
diagrams for vapor-liquid equilibrium (VLE) of a benzene/ethanol mixture. This liquid mixture is non-ideal, and the system has an azeotrope (a condition where liquid and vapor have the same composition). The blue curve represents the liquid-phase boundary (bubble point), and the green curve represents the vapor-phase boundary (dew point). You can vary the temperature for the
P-x-y
diagram. Click and drag the black dot to change the benzene mole fraction and the temperature or pressure. The bar chart shows the moles of liquid (blue) and vapor (green) in equilibrium and the mole fraction of benzene in each phase (
x
B
for liquid phase,
y
B
for vapor phase). The system contains a total of 1 mol. The activity coefficients used in the modified Raoult's law are calculated using the two-parameter Margules model.

Details

The saturation pressure of component
i=(1,2)
is calculated using the Antoine equation:
sat
P
i
=
A
i
-
B
i
T+
C
i
10
,
where
A
i
,
B
i
, and
C
i
are Antoine constants,
T
is temperature (°C), and
sat
P
i
is in bar.
The two-parameter Margules model is used to calculate the activity coefficients for a non-ideal liquid mixture of benzene
1
and ethanol
2
. This model fits the excess Gibbs free energy:
E
G
RT
=
x
1
x
2
(
A
21
x
1
+
A
12
x
2
)
,
where
E
G
is excess Gibbs energy, and
R
is the ideal gas constant.
The activity coefficients
γ
1
,
γ
2
are given by:
ln
γ
1
=
2
x
2
(
A
12
+2(
A
21
-
A
12
)
x
1
)
,
ln
γ
2
=
2
x
1
(
A
21
+2(
A
12
-
A
21
)
x
2
)
,
where
x
1
and
x
2
are the liquid mole fractions of benzene and ethanol and
x
1
+
x
2
=1
, and
A
21
and
A
12
are the Margules parameters.
The modified Raoult's law is used to calculate the bubble-point and dew-point pressures using the
K
factors:
K
i
=
y
i
x
i
=
γ
i
sat
P
i
P
,
where
y
i
is the vapor mole fraction and
y
1
+
y
2
=1
, and
P
is the total pressure (bar).
Bubble-point pressure calculation:
P=
x
1
γ
1
sat
P
1
+
x
2
γ
2
sat
P
2
.
Dew-point pressure calculation:
P=
-1
y
1
γ
1
sat
P
1
+
y
2
γ
2
sat
P
2
.
The screencast video at[2] shows how to use this Demonstration.

References

[1] J. R. Elliott and C. T. Lira, Introductory Chemical Engineering Thermodynamics, New York: Prentice Hall, 2012 pp. 372–377, 430.
[2] Vapor-Liquid Equilibrium Diagram for Non-Ideal Mixture. sites.google.com/a/learncheme.com/learncheme/simulations/thermodynamics/thermo-2/vle-diagram-for-non-ideal-mixtures.

External Links

Non-Ideal Vapor-Liquid Equilibrium (VLE) Modeled by the Margules Equation
P-x-y and T-x-y Diagrams for Vapor-Liquid Equilibrium (VLE)
Mass Balances for Binary Vapor-Liquid Equilibrium (VLE)

Permanent Citation

Rachael L. Baumann, John L. Falconer, Nick Bongiardina
​
​"Vapor-Liquid Equilibrium Diagram for Non-Ideal Mixture"​
​http://demonstrations.wolfram.com/VaporLiquidEquilibriumDiagramForNonIdealMixture/​
​Wolfram Demonstrations Project​
​Published: June 17, 2014