One of the easiest ways to begin to explore fugacities in nonideal solutions is to model solubilities of crystalline solids dissolved in high pressure gases. In this case, the crystalline solids remain as a pure phase in equilibrium with a vapor mixture, and the fugacity of the “solid” component must be the same in the crystalline phase as in the vapor phase. Consider biphenyl dissolved in carbon dioxide, using kij = 0.100. The molar volume of crystalline biphenyl is 156 cm3/mol.
(a) Calculate the fugacity (in MPa) of pure crystalline biphenyl at 310 K and 330 K and 0.1, 1, 10, 15, and 20 MPa.
(b) Calculate and plot the biphenyl solubility for the isotherm over the pressure range. Compare the solubility to the ideal gas solubility of biphenyl where the Poynting correction is included, but the gas phase nonidealities are ignored.