0% found this document useful (0 votes)
226 views8 pages

VLE Exercise Solution

This document contains the solutions to an exercise involving vapor-liquid equilibrium calculations for a binary system using thermodynamic charts. 1) Calculations are shown to determine (a) vapor composition y1 and pressure P given liquid composition x1 and temperature T for bubble point, (b) liquid composition x1 and pressure P given vapor composition y1 and temperature T for dew point, (c) temperature T given liquid composition x1 and pressure P for bubble point, and (d) temperature T given vapor composition y1 and pressure P for dew point. 2) Using De Priester charts for a methane/ethylene/ethane system, calculations are shown to find (a) bubble point pressure given

Uploaded by

Hiba
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
226 views8 pages

VLE Exercise Solution

This document contains the solutions to an exercise involving vapor-liquid equilibrium calculations for a binary system using thermodynamic charts. 1) Calculations are shown to determine (a) vapor composition y1 and pressure P given liquid composition x1 and temperature T for bubble point, (b) liquid composition x1 and pressure P given vapor composition y1 and temperature T for dew point, (c) temperature T given liquid composition x1 and pressure P for bubble point, and (d) temperature T given vapor composition y1 and pressure P for dew point. 2) Using De Priester charts for a methane/ethylene/ethane system, calculations are shown to find (a) bubble point pressure given

Uploaded by

Hiba
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 8

UNIVERSITI MALAYSIA PERLIS

Pusat Pengajian Kejuruteraan Bioproses

ERT 206: Thermodynamics

Exercise 3 Solutions

1)
(a) Given x1 = 0.33 and T= 100 °C, find y1 and P.

BUBL P Calculation
Calculation:

Step 1: Find P1sat & P2sat using Antoine equation.

B
ln Pi sat / kPa = A −
T / °C + C

2726.81
ln P1sat / kPa = 13.7819 −
100°C + 217.572
= 5.195

sat
Thus, P1 = 180.45kPa

3056.96
ln P2sat / kPa = 13.932 −
100°C + 217.625

= 4.3076

sat
Thus, P2 = 74.2597kPa

Step 2: Find P.

P = x1 P1sat + x2 P2sat

= x1 P1sat + (1 − x1 )P2sat

= 0.33(180.45kPa ) + (1 − 0.33)(74.2597 kPa )


= 109.3025kPa
Step 3: Calculate y1 using Raoult’s law.

yi P = xi Pi sat
y1 = x1 P1sat P = 0.33(180.45kPa ) 109.3025kPa = 0.5448
(b) Given y1 = 0.33 and T= 100 °C, find x1 and P.

DEW P Calculation:

Step 1: Find P from Raoult’s Law assuming ∑xi i =1


1
P= sat
y1 / P1 + y 2 / P2sat

1
=
0.33 / 180.45 + (1 − 0.33) / 74.2597

= 92.156kPa

Step 2: Find x

yi P = xi Pi sat
x1 = y1 P P1sat
= (0.33)(92.156 ) (180.45)
= 0.1685

(c) Given x1 = 0.33 and P= 120 kPa, find y1 and T.

BUBL T Calculation:

Step 1: Start with α=1, find P2sat


P 120
P2sat = = = 120kPa
x1α + x 2 (0.33)(1) + (1 − 0.33)

Step 2: Find T using Antoine eq & substitute P2sat obtained in step 1


B
T= −C
A − ln P2sat

3056.96
= − 217.625
13.9320 − ln 120

= 116.67°C
Step 3: Find new α by substituting T

* find the expression for α first


2726.81
ln P1sat / kPa = 13.7819 −
T / °C + 217.572

3056.96
- ln P2sat / kPa = 13.932 −
T / °C + 217.625

2726.81 3056.96
( )
ln P1sat / P2sat = −0.1501 − +
T / °C + 217.572 T / °C + 217.625

2726.81 3056.96
ln α = −0.1501 − +
T / °C + 217.572 T / °C + 217.625

* substituting T=116.67̊C in the equation;


2726.81 3056.96
ln α = −0.1501 − +
116.67°C + 217.572 116.67°C + 217.625
= 0.8362
α = 2.3076
Step 4: Repeat step 1 by using new α until similar value of α is obtained

2nd iteration: 3rd iteration:


α = 2.3076 α = 2.3985
P2sat = 83.8277 kPa P2sat = 82.10714kPa
T = 104.0507°C T = 103.3502°C
α new = 2.3985 α new = 2.40387

4th iteration: 5th iteration:


α = 2.40387 α = 2.4042
P2sat = 82.00765kPa P2sat = 82.00kPa
T = 103.3094°C T = 103.307°C
α new = 2.4042 α new = 2.4042

Thus, T = 103.307°C

Step 5: Find P1sat & find y1 using Raoult’s law

P1sat = αP2sat y1 = x1 P1sat P


= 2.4042(82.00kPa ) = 0.33(197.144kPa ) 120kPa
= 197.144kPa = 0.542
(d) Given y1 = 0.33 and P= 120 kPa, find x1 and T.

DEW T Calculation:

Step 1: Start with α=1, find P1sat

P1sat = P( y1 + y2α )
= 120[0.33 + 0.67(1)]
= 120kPa

Step 2: Find T using Antoine eq & substitute P1sat obtained in step 1


B
T= −C
A − ln P1sat

2726.81
= − 217.572
13.7819 − ln 120

= 85.595°C

Step 3: Find new α by substituting T


2726.81 3056.96
ln α = −0.1501 − +
T / °C + 217.572 T / °C + 217.625
2726.81 3056.96
ln α = −0.1501 − +
85.595 / °C + 217.572 85.595 / °C + 217.625
α = 2.5527
Step 4: Repeat step 1 by using new α until similar value of α is obtained

2nd iteration: 3rd iteration:


α = 2.5527 α = 2.34215
P1sat = 244.8347kPa P1sat = 227.9091kPa
T = 111.7005°C T = 108.8766°C
α new = 2.34215 α new = 2.36250

4th iteration: 5th iteration:


α = 2.36250 α = 2.3604
P1sat = 229.5446 Pa P1sat = 229.3807kPa
T = 109.1563°C T = 109.1284°C
α new = 2.3604 α new = 2.3606
Step 5: Find x1 using Raoult’s law

y1 = x1 P1sat P

x1 = y1 P P1sat

= 0.33(120kPa ) 229.3807kPa
= 0.1726

2) Assuming the validity of the De Priester charts, make the following VLE calculations
for methane(1)/ethylene(2)/ethane(3) system:
(a) BUBL P, given x1 = 0.10, x2 = 0.50, and T = -60 °F

P = 200 psia P = 250 psia P = 215 psia


Component xi
Ki yi=Kixi Ki yi=Kixi Ki yi=Kixi
Methane 0.10 5.600 0.560 4.600 0.460 5.150 0.515
Ethylene 0.50 0.700 0.350 0.575 0.288 0.650 0.325
Ethane 0.40 0.445 0.178 0.380 0.152 0.420 0.168

SUM = 1.088 SUM = 0.900 SUM = 1.008


(close enough)

Thus, BUBL P = 215 psia

(b) DEW P, given y1 = 0.50, y2 = 0.25, and T = -60 °F

P = 190 psia P = 200 psia


Component yi
Ki xi=yi/Ki Ki xi=yi/Ki
Methane 0.50 5.900 0.085 5.600 0.089
Ethylene 0.25 0.730 0.342 0.700 0.357
Ethane 0.25 0.460 0.543 0.445 0.562
SUM = 1.008
SUM = 0.971
(close enough)

Thus, DEW P = 210 psia


(c) BUBL T, given x1 = 0.12, x2 = 0.40, and P = 250 psia

T = -50 F T = -60 F T = -57 F


Component xi
Ki yi=Kixi Ki yi=Kixi Ki yi=Kixi
Methane 0.12 4.900 0.588 4.600 0.552 4.700 0.564
Ethylene 0.40 0.680 0.272 0.570 0.228 0.615 0.246
Ethane 0.48 0.450 0.216 0.380 0.182 0.405 0.194

SUM = 1.076 SUM = 0.962 SUM = 1.004


(close enough)

Thus, BUBL T = -57 F

(d) DEW T, given y1 = 0.43, y2 = 0.36, and P = 250 psia


T = -40 F T = -50 F T = -45 F
Component yi
Ki xi= yi /Ki Ki xi= yi /Ki Ki xi= yi /Ki
Methane 0.43 5.200 0.083 4.900 0.088 5.050 0.085
Ethylene 0.36 0.800 0.450 0.680 0.529 0.740 0.486
Ethane 0.21 0.520 0.404 0.450 0.467 0.485 0.433

SUM = 0.937 SUM = 1.084 SUM = 1.005


(close enough)
Thus, DEW T = -45 F

You might also like