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

VLE Thermodynamics Exercise Solutions

1. The document contains solutions to thermodynamics exercises involving vapor-liquid equilibrium calculations using Raoult's law and the Antoine equation for various systems. 2. For a binary system at 100°C, calculations determine that given x1=0.33, y1=0.5448 and P=109.3 kPa. Conversely, given y1=0.33, x1 is calculated to be 0.1685 and P=92.156 kPa. 3. For a ternary methane/ethylene/ethane system, vapor-liquid calculations using De Priester charts determine that given inputs of temperature, pressure and compositions, the equilibrium pressures and temperatures can be

Uploaded by

Zidni Ilma
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)
854 views8 pages

VLE Thermodynamics Exercise Solutions

1. The document contains solutions to thermodynamics exercises involving vapor-liquid equilibrium calculations using Raoult's law and the Antoine equation for various systems. 2. For a binary system at 100°C, calculations determine that given x1=0.33, y1=0.5448 and P=109.3 kPa. Conversely, given y1=0.33, x1 is calculated to be 0.1685 and P=92.156 kPa. 3. For a ternary methane/ethylene/ethane system, vapor-liquid calculations using De Priester charts determine that given inputs of temperature, pressure and compositions, the equilibrium pressures and temperatures can be

Uploaded by

Zidni Ilma
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
  • Exercise 3 Solutions
  • Appendix B: Properties of Pure Species
  • VLE from K-Value Correlations

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

Common questions

Powered by AI

Relative volatility (α) is a measure of a component's tendency to vaporize as compared to another component in mixture calculations. In BUBL T and DEW T calculations, α is used to correlate temperature with equilibrium compositions. Initially, α is assumed, and saturation temperatures are calculated. Iterative adjustments are made as new values of α are substituted back into the calculations until convergence is reached, reflecting accurate equilibrium states. This iterative adjustment accounts for non-ideal behavior and ensures precise equilibrium predictions through recalibration of assumptions based on results .

To determine DEW T, start with the given gaseous mole fractions and assume an initial guess for the temperature . Calculate each component's saturation pressure using the Antoine equation and convert gaseous mole fractions to liquid using the K-values (equilibrium constants) derived. Iterate by updating the temperature guess based on the liquid phase sum approaching unity (1) across iterations, calculating new K-values each time by reapplying the Antoine parameters until convergence . This process ensures the temperature where the first droplet condenses is accurately identified.

The Antoine equation is preferred over empirical tables for calculating saturation pressures because it provides a continuous mathematical model that can be applied at any temperature within the range of coefficients. This adaptability ensures precision and minimizes interpolation errors that could arise with discrete data points from tables. Furthermore, it aligns with algorithms in simulation software, enabling automated calculations and integration into dynamic temperature-sensitive processes, enhancing the accuracy and efficiency in engineering applications .

BUBL P calculation aims to determine the pressure at which the first bubble of vapor forms from a liquid mixture at a given temperature. It sums partial pressures using Raoult’s law and calculates total pressure from known mole fractions in the liquid phase . DEW P calculation, however, determines the pressure where the first drop of liquid appears from a vapor mixture. It begins with gaseous mole fractions, converts them using K-values to liquid mole fractions, and sums to find total pressure . Methodologically, BUBL P uses the liquid phase mole fractions, while DEW P starts with vapor phase mole fractions .

The iterative process in a BUBL T calculation involves beginning with an initial guess for the relative volatility (α). The saturation temperature for each component is found using the Antoine equation, and the total moles are recalculated. This value is used to adjust α and repeat the calculation until convergence is achieved, where the relative volatilities do not change significantly between iterations. This is necessary because the relationship between temperature and pressure is non-linear, and assumptions made initially might not hold after re-calculation, requiring adjustments for precision .

In a BUBL process, Raoult’s law is applied to each component in a liquid mixture to determine the partial pressure contribution to total pressure, using the equation P_i = x_i * P_sat_i, where P_i is the partial pressure, x_i is the mole fraction, and P_sat_i is the saturation pressure . The vapor-phase mole fraction y_i is then calculated by dividing the component's partial pressure by the total pressure P, i.e., y_i = (x_i * P_sat_i) / P, facilitating the transition from liquid to vapor composition at equilibrium .

To calculate the BUBL P for a binary mixture: Step 1 involves finding the saturation pressure of each component at the given temperature using the Antoine equation . Step 2 involves applying Raoult's law, where the total pressure P is calculated as the sum of the products of mole fraction (x_i) and saturation pressure (P_sat_i) of each component, i.e., P = Σ(x_i * P_sat_i). Finally, Step 3 determines the vapor-phase mole fractions y_i by dividing each component's product of x_i and P_sat_i by the total pressure P .

The De Priester chart provides empirical K-values for components in VLE systems, which are used to convert between vapor and liquid phase compositions. For a given system like methane/ethylene/ethane, the chart helps determine distribution coefficients at specific pressures and temperatures without conducting experimental measurements. These coefficients simplify calculations like BUBL P, DEW P, BUBL T, and DEW T by providing a direct way to obtain the equilibrium constant K, necessary for relating mole fractions in equilibrium calculations .

The Antoine equation is used to calculate the saturation pressure (P_sat) of a component based on its specific coefficients A, B, and C, and the temperature T. The equation is of the form log(P_sat) = A - (B / (T + C)), where the temperature T is in degrees Celsius. By substituting the given temperature into the equation, P_sat can be determined .

Achieving convergence in iterative calculations like BUBL T and DEW T is crucial as it indicates that the solution has reached a stable temperature or pressure that satisfies phase equilibrium relationships and balances the material under specified conditions. Without convergence, results could be arbitrary or divergent, leading to incorrect interpretations. Convergence ensures consistency and reliability, which is critical for designing and optimizing industrial processes where precise understanding of phase behavior is required .

You might also like