Chapter 13, Problem 13.47P

(a)

Interpretation Introduction

Interpretation:

The bubble point pressure for any one of the given binary systems in table $13.10$ is to be calculated using Wilson equation.

Concept Introduction:

Antoine equation is used to determine the vapor pressure of any substance at the given temperature by the equation:

  $\ln P^{\text{sat}}\left(\text{kPa}\right)=A-\frac{B}{\left(T^{\circ }\text{C}\right)+C}$ .......(1)

Here, $A,B\text{ and }C$ are the constants specific for a substance given in table B.2 of appendix B.

Equation $13.19$ to be used for Modified Raoult’s law is:

  $y_{i}P=x_i{\gamma }_i{P}_i{}^{\text{sat}}$ .......(2)

The Bubble point pressure for a binary system in vapor/liquid equilibrium is defined as the pressure where first bubble of vapor appears which is in equilibrium with the liquid present in the system. The equation which defines this pressure at this point is:

  $P_{\text{BUBL}}=x_{\text{1}}{\gamma }_1{P}_{\text{1}}^{\text{sat}}+x_{\text{2}}{\gamma }_2{P}_{\text{2}}^{\text{sat}}$ .......(3)

Wilson equations to be used are:

  $\begin{array}{c}\ln {\gamma }_1=-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\\ \ln {\gamma }_2=-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\end{array}$ .......(4)

Here, the parameters ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ are calculated by the formula:

  $\begin{array}{c}{\Lambda }_{12}=\frac{V_2}{V_1}\exp \frac{-a_{12}}{RT}\\ {\Lambda }_{21}=\frac{V_1}{V_2}\exp \frac{-a_{21}}{RT}\end{array}$ .......(5)

Where, $V_1,V_2,a_{12}\text{ and }{a}_{21}$ are the Wilson equation parameters given in table $13.10$ for binary systems.

Answer to Problem 13.47P

The bubble point pressure for $1-\text{Propanol(1)/Water(2)}$ using Wilson equation is:

  $P_{\text{BUBL}}=27.284\text{ kPa}$

Explanation of Solution

Given information:

The temperature at which the bubble point pressure is to be calculated is ${60}^{\circ }\text{C}$ . The composition of the liquid phase at this point is given as $x_1=0.3$ .

Wilson equation parameters are given in Table 13.10 as shown below:

The binary system for which the bubble point pressure will be calculated is $1-\text{Propanol(1)/Water(2)}$ at $t={60}^{\circ }\text{C}$ . The liquid phase composition is:

  $\begin{array}{c}{x}_1=0.3\\ x_2=0.7\end{array}$

From table B.2 of appendix B, the Antoine equation constants for $1-\text{Propanol(1) and Water(2)}$ are:

  $\begin{array}{c}{A}_1=16.1154\\ B_1=3483.67\\ C_1=205.807\\ A_2=16.3872\\ B_2=3885.70\\ C_2=230.170\end{array}$

Now, use equation (1) to calculate the vapor pressure of $1-\text{Propanol(1) and Water(2)}$ for temperature ${60}^{\circ }\text{C}$ as:

  $\begin{array}{c}\ln P_1{}^{\text{sat}}=16.1154-\frac{3483.67}{\left(60\right)+205.807}\\ P_1{}^{\text{sat}}=e^{\left(16.1154-\frac{3483.67}{\left(60\right)+205.807}\right)}\\ =20.2750\text{ kPa}\\ \ln P_2{}^{\text{sat}}=16.3872-\frac{3885.70}{\left(60\right)+230.170}\\ P_2{}^{\text{sat}}=e^{\left(16.3872-\frac{3885.70}{\left(60\right)+230.170}\right)}\\ =20.0070\text{ kPa}\end{array}$

From table $13.10$, the parameters needed for $1-\text{Propanol(1) and Water(2)}$ to be used in Wilson equation are:

  $\begin{array}{c}{V}_1=18.07{\text{ cm}}^3\text{/mol}\\ V_2=75.14{\text{ cm}}^3\text{/mol}\\ a_{12}=775.48\text{ cal/mol}\\ a_{21}=1351.90\text{ cal/mol}\end{array}$

The value of universal gas constant to be used is,

  $R=1.9859\text{ cal/}\left(\text{mol}\cdot \text{K}\right)$

Now, use equation (5) to calculate the values of ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ as:

  $\begin{array}{c}{\Lambda }_{12}=\frac{V_2}{V_1}\exp \frac{-a_{12}}{RT}\\ =\frac{18.07\frac{{\text{cm}}^3}{\text{mol}}}{75.14\frac{{\text{cm}}^3}{\text{mol}}}\exp \frac{-775.48\frac{\text{cal}}{\text{mol}}}{\left(1.9859\frac{\text{cal}}{\text{mol}\cdot \text{K}}\right)\left(333.15\text{ K}\right)}\\ =0.07448\\ {\Lambda }_{21}=\frac{V_1}{V_2}\exp \frac{-a_{21}}{RT}\\ =\frac{75.14\frac{{\text{cm}}^3}{\text{mol}}}{18.07\frac{{\text{cm}}^3}{\text{mol}}}\exp \frac{-1351.90\frac{\text{cal}}{\text{mol}}}{\left(1.9859\frac{\text{cal}}{\text{mol}\cdot \text{K}}\right)\left(333.15\text{ K}\right)}\\ =0.53887\end{array}$

Now, use these values of ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ along with the values of $x_1\text{ and }{x}_2$ and calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.3+0.7\left(0.07448\right)\right)+0.7\left(\frac{0.07448}{0.3+0.7\left(0.07448\right)}-\frac{0.53887}{0.7+0.3\left(0.53887\right)}\right)\right]\\ =2.12554\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.7+0.3\left(0.53887\right)\right)+0.3\left(\frac{0.07448}{0.3+0.7\left(0.07448\right)}-\frac{0.53887}{0.7+0.3\left(0.53887\right)}\right)\right]\\ =1.02504\end{array}$

Calculate the bubble point pressure of the system using equation (3) as:

  $\begin{array}{c}{P}_{\text{BUBL}}=x_{\text{1}}{\gamma }_1{P}_{\text{1}}^{\text{sat}}+x_{\text{2}}{\gamma }_2{P}_{\text{2}}^{\text{sat}}\\ =0.3\left(2.12554\right)\left(20.2750\right)+0.7\left(1.02504\right)\left(20.0070\right)\\ =27.284\text{ kPa}\end{array}$

(b)

Interpretation Introduction

Interpretation:

The dew point pressure for any one of the given binary systems in table $13.10$ is to be calculated using Wilson equation.

Concept Introduction:

Equation $13.19$ to be used for Modified Raoult’s law is:

  $y_{i}P=x_i{\gamma }_i{P}_i{}^{\text{sat}}$ .......(2)

Wilson equations to be used are:

  $\begin{array}{c}\ln {\gamma }_1=-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\\ \ln {\gamma }_2=-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\end{array}$ .......(4)

Here, the parameters ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ are calculated by the formula:

  $\begin{array}{c}{\Lambda }_{12}=\frac{V_2}{V_1}\exp \frac{-a_{12}}{RT}\\ {\Lambda }_{21}=\frac{V_1}{V_2}\exp \frac{-a_{21}}{RT}\end{array}$ .......(5)

Where, $V_1,V_2,a_{12}\text{ and }{a}_{21}$ are the Wilson equation parameters given in table $13.10$ for binary systems.

The Dew point pressure for a binary system in vapor/liquid equilibrium is defined as the pressure where first drop of liquid appears which is in equilibrium with the vapor present in the system at a particular temperature. The equation that defines this pressure at this point is:

  $P_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}$ .......(6)

Answer to Problem 13.47P

The dew point pressure for $1-\text{Propanol(1)/Water(2)}$ using Wilson equation is:

  $P_{\text{DEW}}=27.993\text{ kPa}$

Explanation of Solution

Given information:

The temperature at which the dew point pressure is to be calculated is ${60}^{\circ }\text{C}$ . The composition of the vapor phase at this point is given as $y_1=0.3$ .

Wilson equation parameters are given in Table 13.10 as shown below:

Use the values of $P_1{}^{\text{sat}},\text{ and }{P}_2{}^{\text{sat}}$ at ${60}^{\circ }\text{C}$ as calculated in part (a). The vapor phase composition is:

  $\begin{array}{c}{y}_1=0.3\\ y_2=0.7\end{array}$

From table $13.10$, the parameters needed for $1-\text{Propanol(1) and Water(2)}$ to be used in Wilson equation are:

  $\begin{array}{c}{V}_1=18.07{\text{ cm}}^3\text{/mol}\\ V_2=75.14{\text{ cm}}^3\text{/mol}\\ a_{12}=775.48\text{ cal/mol}\\ a_{21}=1351.90\text{ cal/mol}\end{array}$

The value of universal gas constant to be used is,

  $R=1.9859\text{ cal/}\left(\text{mol}\cdot \text{K}\right)$

Now, use equation (5) to calculate the values of ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ as:

  $\begin{array}{c}{\Lambda }_{12}=\frac{V_2}{V_1}\exp \frac{-a_{12}}{RT}\\ =\frac{18.07\frac{{\text{cm}}^3}{\text{mol}}}{75.14\frac{{\text{cm}}^3}{\text{mol}}}\exp \frac{-775.48\frac{\text{cal}}{\text{mol}}}{\left(1.9859\frac{\text{cal}}{\text{mol}\cdot \text{K}}\right)\left(333.15\text{ K}\right)}\\ =0.07448\\ {\Lambda }_{21}=\frac{V_1}{V_2}\exp \frac{-a_{21}}{RT}\\ =\frac{75.14\frac{{\text{cm}}^3}{\text{mol}}}{18.07\frac{{\text{cm}}^3}{\text{mol}}}\exp \frac{-1351.90\frac{\text{cal}}{\text{mol}}}{\left(1.9859\frac{\text{cal}}{\text{mol}\cdot \text{K}}\right)\left(333.15\text{ K}\right)}\\ =0.53887\end{array}$

1^st iteration:

Now, use these values of ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ along with the values of $x_1\text{ and }{x}_2$ guessed as $0.1\text{ and 0}\text{.9}$ respectively. Calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.1+\left(0.9\right)\left(0.07448\right)\right)+\left(0.9\right)\left(\frac{0.07448}{0.1+\left(0.9\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.9\right)+0.1\left(0.53887\right)}\right)\right]\\ =5.379\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(\left(0.9\right)+0.1\left(0.53887\right)\right)+0.1\left(\frac{0.07448}{0.1+\left(0.9\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.9\right)+0.1\left(0.53887\right)}\right)\right]\\ =1.036\end{array}$

Now, calculate the dew point pressure of the system using equation (6) as:

  $\begin{array}{c}{P}_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}\\ =\frac{1}{\frac{0.3}{\left(5.379\right)\left(20.2750\right)}+\frac{0.7}{\left(1.036\right)\left(20.0070\right)}}\\ =27.38\text{ kPa}\end{array}$

Apply Raoult’s law on both the components and use this dew point pressure to calculate $x_1\text{ and }{x}_2$ as:

  $\begin{array}{c}{y}_1{P}_{\text{DEW}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}\\ 0.3\left(27.38\right)=x_1\left(5.379\right)\left(20.2750\right)\\ x_1=0.075\\ y_2{P}_{\text{DEW}}=x_2{\gamma }_2{P}_2{}^{\text{sat}}\\ 0.7\left(27.38\right)=\left(x_2\right)\left(1.036\right)\left(20.0070\right)\\ x_2=0.925\end{array}$

Now, use this calculated value of $x_1\text{ and }{x}_2$ to again calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ and further calculate $P_{\text{DEW}}$ .

2^nd iteration:

Calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) and calculated values of $x_1\text{ and }{x}_2$ from 1^st iteration as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.075+\left(0.925\right)\left(0.07448\right)\right)+\left(0.925\right)\left(\frac{0.07448}{0.075+\left(0.925\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.925\right)+0.075\left(0.53887\right)}\right)\right]\\ =6.694\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(\left(0.925\right)+0.075\left(0.53887\right)\right)+0.075\left(\frac{0.07448}{0.075+\left(0.925\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.925\right)+0.075\left(0.53887\right)}\right)\right]\\ =1.033\end{array}$

Now, calculate the dew point pressure of the system using equation (6) as:

  $\begin{array}{c}{P}_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}\\ =\frac{1}{\frac{0.3}{\left(6.694\right)\left(20.2750\right)}+\frac{0.7}{\left(1.033\right)\left(20.0070\right)}}\\ =27.716\text{ kPa}\end{array}$

Apply Raoult’s law on both the components and use this dew point pressure to calculate $x_1\text{ and }{x}_2$ as:

  $\begin{array}{c}{y}_1{P}_{\text{DEW}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}\\ 0.3\left(27.716\right)=x_1\left(6.694\right)\left(20.2750\right)\\ x_1=0.061\\ y_2{P}_{\text{DEW}}=x_2{\gamma }_2{P}_2{}^{\text{sat}}\\ 0.7\left(27.716\right)=\left(x_2\right)\left(1.033\right)\left(20.0070\right)\\ x_2=0.939\end{array}$

3^rd iteration:

Calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) and calculated values of $x_1\text{ and }{x}_2$ from 2^nditeration as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.061+\left(0.939\right)\left(0.07448\right)\right)+\left(0.939\right)\left(\frac{0.07448}{0.061+\left(0.939\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.939\right)+0.061\left(0.53887\right)}\right)\right]\\ =7.741\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(\left(0.939\right)+0.061\left(0.53887\right)\right)+0.061\left(\frac{0.07448}{0.061+\left(0.939\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.939\right)+0.061\left(0.53887\right)}\right)\right]\\ =1.030\end{array}$

Now, calculate the dew point pressure of the system using equation (6) as:

  $\begin{array}{c}{P}_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}\\ =\frac{1}{\frac{0.3}{\left(7.741\right)\left(20.2750\right)}+\frac{0.7}{\left(1.030\right)\left(20.0070\right)}}\\ =27.871\text{ kPa}\end{array}$

Apply Raoult’s law on both the components and use this dew point pressure to calculate $x_1\text{ and }{x}_2$ as:

  $\begin{array}{c}{y}_1{P}_{\text{DEW}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}\\ 0.3\left(27.871\right)=x_1\left(7.741\right)\left(20.2750\right)\\ x_1=0.053\\ y_2{P}_{\text{DEW}}=x_2{\gamma }_2{P}_2{}^{\text{sat}}\\ 0.7\left(27.871\right)=\left(x_2\right)\left(1.030\right)\left(20.0070\right)\\ x_2=0.947\end{array}$

4^thiteration:

Calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) and calculated values of $x_1\text{ and }{x}_2$ from 3^rd iteration as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.053+\left(0.947\right)\left(0.07448\right)\right)+\left(0.947\right)\left(\frac{0.07448}{0.053+\left(0.947\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.947\right)+0.053\left(0.53887\right)}\right)\right]\\ =8.492\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(\left(0.947\right)+0.053\left(0.53887\right)\right)+0.053\left(\frac{0.07448}{0.053+\left(0.947\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.947\right)+0.053\left(0.53887\right)}\right)\right]\\ =1.028\end{array}$

Now, calculate the dew point pressure of the system using equation (6) as:

  $\begin{array}{c}{P}_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}\\ =\frac{1}{\frac{0.3}{\left(8.492\right)\left(20.2750\right)}+\frac{0.7}{\left(1.028\right)\left(20.0070\right)}}\\ =27.951\text{ kPa}\end{array}$

Apply Raoult’s law on both the components and use this dew point pressure to calculate $x_1\text{ and }{x}_2$ as:

  $\begin{array}{c}{y}_1{P}_{\text{DEW}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}\\ 0.3\left(27.951\right)=x_1\left(8.492\right)\left(20.2750\right)\\ x_1=0.049\\ y_2{P}_{\text{DEW}}=x_2{\gamma }_2{P}_2{}^{\text{sat}}\\ 0.7\left(27.951\right)=\left(x_2\right)\left(1.028\right)\left(20.0070\right)\\ x_2=0.951\end{array}$

5^th iteration:

Calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) and calculated values of $x_1\text{ and }{x}_2$ from 4^th iteration as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.049+\left(0.951\right)\left(0.07448\right)\right)+\left(0.951\right)\left(\frac{0.07448}{0.049+\left(0.951\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.951\right)+0.049\left(0.53887\right)}\right)\right]\\ =8.921\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(\left(0.951\right)+0.049\left(0.53887\right)\right)+0.049\left(\frac{0.07448}{0.049+\left(0.951\right)\left(0.07448\right)}-\frac{0.53887}{\left(0.951\right)+0.049\left(0.53887\right)}\right)\right]\\ =1.027\end{array}$

Now, calculate the dew point pressure of the system using equation (6) as:

  $\begin{array}{c}{P}_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}\\ =\frac{1}{\frac{0.3}{\left(8.942\right)\left(20.2750\right)}+\frac{0.7}{\left(1.027\right)\left(20.0070\right)}}\\ =27.993\text{ kPa}\end{array}$

Apply Raoult’s law on both the components and use this dew point pressure to calculate $x_1\text{ and }{x}_2$ as:

  $\begin{array}{c}{y}_1{P}_{\text{DEW}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}\\ 0.3\left(27.993\right)=x_1\left(8.942\right)\left(20.2750\right)\\ x_1=0.046\\ y_2{P}_{\text{DEW}}=x_2{\gamma }_2{P}_2{}^{\text{sat}}\\ 0.7\left(27.993\right)=\left(x_2\right)\left(1.027\right)\left(20.0070\right)\\ x_2=0.954\end{array}$

Since, ${\left(x_1\right)}_{5^{\text{th}}}-{\left(x_1\right)}_{4^{\text{th}}}=0.003$, this value of $P_{\text{DEW}}$ calculated in 5^th iteration is well accepted.

Therefore,

  $P_{\text{DEW}}=27.993\text{ kPa}$

(c)

Interpretation Introduction

Interpretation:

$P,T-$ flash is to be performed for any one of the given binary systems in table $13.10$ using Wilson equation.

Concept Introduction:

Equation $13.19$ to be used for Modified Raoult’s law is:

  $y_{i}P=x_i{\gamma }_i{P}_i{}^{\text{sat}}$ .......(2)

The Bubble point pressure for a binary system in vapor/liquid equilibrium is defined as the pressure where first bubble of vapor appears which is in equilibrium with the liquid present in the system. The equation which defines this pressure at this point is:

  $P_{\text{BUBL}}=x_{\text{1}}{\gamma }_1{P}_{\text{1}}^{\text{sat}}+x_{\text{2}}{\gamma }_2{P}_{\text{2}}^{\text{sat}}$ .......(3)

The Dew point pressure for a binary system in vapor/liquid equilibrium is defined as the pressure where first drop of liquid appears which is in equilibrium with the vapor present in the system at a particular temperature. The equation that defines this pressure at this point is:

  $P_{\text{DEW}}=\frac{1}{\frac{y_{\text{1}}}{{\gamma }_1{P}_{\text{1}}^{\text{sat}}}+\frac{y_{\text{2}}}{{\gamma }_2{P}_{\text{2}}^{\text{sat}}}}$ .......(6)

  $P,T-$ flash calculations are done to find the values of $\mathcal{V}\text{, }\mathcal{L},y_i,\text{ and }{x}_i$ for a system which flashes to produce a two-phase system of vapor and liquid in equilibrium.

The equation for equilibrium ratio, $K_i$ also known as K -value is:

  $K_i\equiv \frac{y_i}{x_i}$ .......(7)

Here, $y_i$ is the vapor phase mole fraction and $x_i$ is the liquid phase mole fraction of species $i$ in the binary system in VLE.

The equations for flash calculations to be used are:

  $\mathcal{L}+\mathcal{V}=1$ .......(8)

Here, $\mathcal{V}$ is the total moles of vapor and $\mathcal{L}$ is the total moles of liquid.

In terms of $K_i$, the equation relating $y_i\text{ and }{z}_i$ is:

  $y_i=\frac{z_i{K}_i}{1+\mathcal{V}\left(K_i-1\right)}$ .......(9)

Here, $z_i$ is the total mole fraction of component $i$ in the system.

Answer to Problem 13.47P

The result of the $P,T-$ flash calculations is:

  $\begin{array}{c}\mathcal{V}=0.837\\ \mathcal{L}=0.163\\ y_1=0.338\\ y_2=0.662\\ x_1=0.215\\ x_2=0.895\end{array}$

Explanation of Solution

Given information:

The flash temperature at which the $P,T-$ flash is to be calculated is ${60}^{\circ }\text{C}$ . The total composition of component 1 at this point is given as $z_1=0.3$ .

The condition for the flash pressure for this system is,

  $P=\frac{1}{2}\left(P_{\text{BUBL}}+P_{\text{DEW}}\right)$

Wilson equation parameters are given in Table 13.10 as shown below:

Use the values of $P_1{}^{\text{sat}},\text{ and }{P}_2{}^{\text{sat}}$ at ${60}^{\circ }\text{C}$ as calculated in part (a).

To perform $P,T-$ flash calculations, first find the bubble point and dew point pressure for the binary system $1-\text{Propanol(1) and Water(2)}$ at temperature ${60}^{\circ }\text{C}$ .

To calculate bubble point pressure, let

  $\begin{array}{c}{z}_1=x_1=0.3\\ z_2=x_2=0.7\end{array}$

Since, the given conditions are same as in part (a), the calculated value of $P_{\text{BUBL}}$ as in part (a) is:

  $P_{\text{BUBL}}=27.284\text{ kPa}$

To calculate dew point pressure, let

  $\begin{array}{c}{z}_1=y_1=0.3\\ z_2=y_2=0.7\end{array}$

Since, the given conditions are same as in part (a), the calculated value of $P_{\text{DEW}}$ as in part (a) is:

  $P_{\text{DEW}}=27.993\text{ kPa}$

From the given condition of the flash pressure, it is calculated as:

  $\begin{array}{c}P=\frac{1}{2}\left(P_{\text{BUBL}}+P_{\text{DEW}}\right)\\ =\frac{1}{2}\left(27.284+27.993\right)\\ =27.639\text{ kPa}\end{array}$

Now, using the modified Raoult’s law, calculate the values of equilibrium ratio of component 1 and 2 using equations (2) and (7) as:

  $\begin{array}{c}{K}_i=\frac{y_i}{x_i}\\ =\frac{{\gamma }_i{P}_i{}^{\text{sat}}}{P}\\ K_1=\frac{{\gamma }_1{P}_1{}^{\text{sat}}}{P}=\frac{\left(2.12554\right)\left(20.2750\right)}{27.369}=1.57\\ K_2=\frac{{\gamma }_2{P}_2{}^{\text{sat}}}{P}=\frac{\left(1.02504\right)\left(20.0070\right)}{27.639}=0.74\end{array}$

Now, use equation (9) and write it for both the component, 1 and 2 as shown below:

  $\begin{array}{c}{y}_1=\frac{0.3\left(1.57\right)}{1+\mathcal{V}\left(1.57-1\right)}\\ y_2=\frac{0.7\left(0.74\right)}{1+\mathcal{V}\left(0.74-1\right)}\end{array}$ .......(10)

Since, $y_1+y_2=1$, the above equations are solved to calculate $\mathcal{V}$ as:

  $\begin{array}{c}\frac{0.3\left(1.57\right)}{1+\mathcal{V}\left(1.57-1\right)}+\frac{0.7\left(0.74\right)}{1+\mathcal{V}\left(0.74-1\right)}=1\\ \mathcal{V}=0.837\end{array}$

Now, use equation (8) to calculate the value of $\mathcal{L}$ as:

  $\begin{array}{c}\mathcal{L}+\mathcal{V}=1\\ \mathcal{L}+0.837=1\\ \mathcal{L}=0.163\end{array}$

Also, use the calculated value of $\mathcal{V}$ in equation set (10) and find the values of $y_1\text{ and }{y}_2$ as:

  $\begin{array}{c}{y}_1=\frac{0.3\left(1.57\right)}{1+\mathcal{V}\left(1.57-1\right)}=\frac{0.3\left(1.57\right)}{1+0.837\left(1.57-1\right)}=0.338\\ y_2=\frac{0.7\left(0.74\right)}{1+\mathcal{V}\left(0.74-1\right)}=\frac{0.7\left(0.74\right)}{1+0.837\left(0.74-1\right)}=0.662\end{array}$

Using these values and the calculated values of $K_1\text{ and }{K}_2$, find the values of $x_1\text{ and }{x}_2$ by equation (7) as:

  $\begin{array}{c}{K}_i=\frac{y_i}{x_i}\\ x_i=\frac{y_i}{K_i}\\ x_1=\frac{y_1}{K_1}=\frac{0.338}{1.57}=0.215\\ x_2=\frac{y_2}{K_2}=\frac{0.662}{0.74}=0.895\end{array}$

The result of the $P,T-$ flash calculations is:

  $\begin{array}{c}\mathcal{V}=0.837\\ \mathcal{L}=0.163\\ y_1=0.338\\ y_2=0.662\\ x_1=0.215\\ x_2=0.895\end{array}$

(d)

Interpretation Introduction

Interpretation:

The values of the azeotropic pressure and composition of the system is to be calculated if it exists for the given binary system.

Concept Introduction:

Equation $13.19$ to be used for Modified Raoult’s law is:

  $y_{i}P=x_i{\gamma }_i{P}_i{}^{\text{sat}}$ .......(2)

Wilson equations to be used are:

  $\begin{array}{c}\ln {\gamma }_1=-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\\ \ln {\gamma }_2=-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\end{array}$ .......(4)

Here, the parameters ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ are calculated by the formula:

  $\begin{array}{c}{\Lambda }_{12}=\frac{V_2}{V_1}\exp \frac{-a_{12}}{RT}\\ {\Lambda }_{21}=\frac{V_1}{V_2}\exp \frac{-a_{21}}{RT}\end{array}$ .......(5)

Where, $V_1,V_2,a_{12}\text{ and }{a}_{21}$ are the Wilson equation parameters given in table $13.10$ for binary systems.

Relative volatility is defined by,

  ${\alpha }_{12}\equiv \frac{y_1/x_1}{y_2/x_2}$ .......(11)

When ${\left({\alpha }_{12}\right)}_{x_1=0}>1\text{ and }{\left({\alpha }_{12}\right)}_{x_1=1}<1$, there exists an azeotrope in the given system.

At the azeotropic point, ${\alpha }_{12}=1\text{ and }{x}_1=y_1$ .

Answer to Problem 13.47P

The azeotropic values of pressure and composition for the binary system is calculated as:

  $\begin{array}{c}{P}^{\text{az}}=31.44\text{ kPa}\\ x_1=y_1=0.437\end{array}$

Explanation of Solution

Given information:

The temperature at which the azeotrope of the system may exists is ${60}^{\circ }\text{C}$ .

Wilson equation parameters are given in Table 13.10 as shown below:

Use the values of $P_1{}^{\text{sat}},P_2{}^{\text{sat}},{\gamma }_1,\text{ and }{\gamma }_2$ at ${60}^{\circ }\text{C}$ as calculated in part (a).

From table $13.10$, the parameters needed for $1-\text{Propanol(1) and Water(2)}$ to be used in Wilson equation are:

  $\begin{array}{c}{V}_1=18.07{\text{ cm}}^3\text{/mol}\\ V_2=75.14{\text{ cm}}^3\text{/mol}\\ a_{12}=775.48\text{ cal/mol}\\ a_{21}=1351.90\text{ cal/mol}\end{array}$

The value of universal gas constant to be used is,

  $R=1.9859\text{ cal/}\left(\text{mol}\cdot \text{K}\right)$

Now, use equation (5) to calculate the values of ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ as:

  $\begin{array}{c}{\Lambda }_{12}=\frac{V_2}{V_1}\exp \frac{-a_{12}}{RT}\\ =\frac{18.07\frac{{\text{cm}}^3}{\text{mol}}}{75.14\frac{{\text{cm}}^3}{\text{mol}}}\exp \frac{-775.48\frac{\text{cal}}{\text{mol}}}{\left(1.9859\frac{\text{cal}}{\text{mol}\cdot \text{K}}\right)\left(333.15\text{ K}\right)}\\ =0.07448\\ {\Lambda }_{21}=\frac{V_1}{V_2}\exp \frac{-a_{21}}{RT}\\ =\frac{75.14\frac{{\text{cm}}^3}{\text{mol}}}{18.07\frac{{\text{cm}}^3}{\text{mol}}}\exp \frac{-1351.90\frac{\text{cal}}{\text{mol}}}{\left(1.9859\frac{\text{cal}}{\text{mol}\cdot \text{K}}\right)\left(333.15\text{ K}\right)}\\ =0.53887\end{array}$

Now, use these values of ${\Lambda }_{12}\text{ and }{\Lambda }_{21}$ along with the values of $x_1=0\text{ and }{x}_2=1$ and calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0+1\left(0.07448\right)\right)+1\left(\frac{0.07448}{0+1\left(0.07448\right)}-\frac{0.53887}{1+0\left(0.53887\right)}\right)\right]\\ =21.2925\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(1+0\left(0.53887\right)\right)+0\left(\frac{0.07448}{0+1\left(0.07448\right)}-\frac{0.53887}{1+0\left(0.53887\right)}\right)\right]\\ =1.00\end{array}$

Using equation (11) along with the modified Raoult’s law, calculate the value of relative volatility at $x_1=0$ as:

  $\begin{array}{c}{\left({\alpha }_{12}\right)}_{x_1=0}\equiv \frac{y_1/x_1}{y_2/x_2}\\ =\frac{{\gamma }_1{P}_1{}^{\text{sat}}/P}{{\gamma }_2{P}_2{}^{\text{sat}}/P}\\ =\frac{{\gamma }_1{P}_1{}^{\text{sat}}}{{\gamma }_2{P}_2{}^{\text{sat}}}\\ =\frac{\left(21.2925\right)\left(20.2750\right)}{\left(1.00\right)\left(20.0070\right)}\\ =21.578\end{array}$

For $x_1=1\text{ and }{x}_2=0$ and calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(1+0\left(0.07448\right)\right)+0\left(\frac{0.07448}{1+0\left(0.07448\right)}-\frac{0.53887}{0+1\left(0.53887\right)}\right)\right]\\ =1.00\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0+1\left(0.53887\right)\right)+0\left(\frac{0.07448}{1+0\left(0.07448\right)}-\frac{0.53887}{0+1\left(0.53887\right)}\right)\right]\\ =1.856\end{array}$

Using equation (11) along with the modified Raoult’s law, calculate the value of relative volatility at $x_1=1$ as:

  $\begin{array}{c}{\left({\alpha }_{12}\right)}_{x_1=1}\equiv \frac{y_1/x_1}{y_2/x_2}\\ =\frac{{\gamma }_1{P}_1{}^{\text{sat}}/P}{{\gamma }_2{P}_2{}^{\text{sat}}/P}\\ =\frac{{\gamma }_1{P}_1{}^{\text{sat}}}{{\gamma }_2{P}_2{}^{\text{sat}}}\\ =\frac{\left(1.00\right)\left(20.2750\right)}{\left(1.856\right)\left(20.0070\right)}\\ =0.546\end{array}$

Since ${\left({\alpha }_{12}\right)}_{x_1=0}>1\text{ and }{\left({\alpha }_{12}\right)}_{x_1=1}<1$, there exists an azeotrope for the binary system.

To calculate the azeotropic pressure, consider the condition ${\alpha }_{12}=1\text{ and }{x}_1=y_1$ .

1^st iteration:

Now, use the values of $x_1\text{ and }{x}_2$ guessed as $0.3\text{ and 0}\text{.7}$, and the values of $y_1\text{ and }{y}_2$ guessed as $0.3\text{ and 0}\text{.7}$, respectively. Calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ using equations set (4) as:

  $\begin{array}{c}{\gamma }_1=\exp \left[-\ln \left(x_1+x_2{\Lambda }_{12}\right)+x_2\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.3+0.7\left(0.07448\right)\right)+0.7\left(\frac{0.07448}{0.3+0.7\left(0.07448\right)}-\frac{0.53887}{0.7+0.3\left(0.53887\right)}\right)\right]\\ =2.12554\\ {\gamma }_2=\exp \left[-\ln \left(x_2+x_1{\Lambda }_{21}\right)+x_1\left(\frac{{\Lambda }_{12}}{x_1+x_2{\Lambda }_{12}}-\frac{{\Lambda }_{21}}{x_2+x_1{\Lambda }_{21}}\right)\right]\\ =\exp \left[-\ln \left(0.7+0.3\left(0.53887\right)\right)+0.3\left(\frac{0.07448}{0.3+0.7\left(0.07448\right)}-\frac{0.53887}{0.7+0.3\left(0.53887\right)}\right)\right]\\ =1.02504\end{array}$

Now, calculate the azeotropic pressure of the system as:

  $\begin{array}{c}{P}^{\text{az}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}+x_2{\gamma }_2{P}_2{}^{\text{sat}}\\ P^{\text{az}}=0.3\left(2.12554\right)\left(20.2750\right)+0.7\left(1.02504\right)\left(20.0070\right)\\ P^{\text{az}}=27.284\text{ kPa}\end{array}$

Apply Raoult’s law on both the components and use this pressure to calculate $y_1\text{ and }{y}_2$ as:

  $\begin{array}{c}{y}_1{P}^{\text{az}}=x_1{\gamma }_1{P}_1{}^{\text{sat}}\\ y_1\left(27.284\right)=0.3\left(1.12554\right)\left(20.2750\right)\\ y_1=0.251\\ y_2=1-y_1\\ y_2=0.749\end{array}$

Now, use this calculated value of $y_1\text{ and }{y}_2$ and substitute it as $x_1\text{ and }{x}_2$ . Again, calculate the values of ${\gamma }_1\text{ and }{\gamma }_2$ and further calculate $P^{\text{az}}$ . After 10 iterations, the values of $P^{\text{az}},x_1,\text{ and }{y}_1$ obtained are:

  $\begin{array}{c}{P}^{\text{az}}=31.44\text{ kPa}\\ x_1=y_1=0.437\end{array}$