temperature (°C) 0

Directions

The fugacities of water and carbon dioxide are calculated as a function of temperature for a closed container, which is a model of a can of soda. The concentrations of the two components are calculated in both the liquid and gas phases. As temperature increases, the pressure increases, and therefore the fugacities increase. Note that the CO2 concentration is much lower than the H2O concentration in the liquid phase, but the CO2 concentration is much higher than the H2O concentration in the gas phase. Because the gas phase is assumed to be ideal, the fugacities of CO2 and H2O in both phases are equal to their gas-phase partial pressures, and thus the CO2 fugacity is much higher than the H2O fugacity. As the temperature increases, the CO2 concentration in liquid water decreases. However, as the temperature increases the pressure increases, and a higher CO2 pressure increases the CO2 concentration in water, so the net effect is that CO2 concentration in the liquid phase does not change much as the temperature increases.

Details

The fugacity of CO2 dissolved in water was calculated using Henry's law and freezing point depression measurements:

\[ \Delta T_{fp} \equiv T_{\mathrm{H_2O}} - T_{sol} = k_{fp}\, m_{\mathrm{CO_2}}, \]
\[ k_{H} = k_{H}^{o}\,\exp\!\left[\dfrac{\delta\,\ln(k_{H})}{\delta\,T^{-1}}\right]\left(\dfrac{1}{T}-\dfrac{1}{273}\right), \]
\[ P_{\mathrm{CO_2}} = \dfrac{\Delta T_{fp}}{k_{H}\,k_{fp}}, \]
\[ f_{\mathrm{CO_2}} = P_{\mathrm{CO_2}}, \]

where \(\Delta T_{fp}\) is the difference between the freezing point of water and the freezing point of carbon dioxide in water (\(T_{\mathrm{H_2O}}\) and \(T_{sol}\)) (K), \(k_{fp}\) is the freezing point depression constant for water ([°C kg]/mol), \(m_{\mathrm{CO_2}}\) is the molality of CO2, \(k_{H}\) is Henry's law constant (kg/[mol bar]), \(k_{H}^{o}\) is Henry's constant at 273 K, \(\frac{\delta\,\ln(k_{H})}{\delta\,T^{-1}}=2400\) is a constant, \(T\) is temperature (K), \(P_{\mathrm{CO_2}}\) is the partial pressure of CO2 (bar), and \(f_{\mathrm{CO_2}}\) is the fugacity of CO2 (bar).

The fugacity of water was calculated from the saturation pressure of water \(P_{\mathrm{H_2O}}^{\,sat}\) using the Antoine equation:

\[ f_{\mathrm{H_2O}} = P_{\mathrm{H_2O}}^{\,sat} = 10^{\,8.07-\frac{1731}{T+233}}, \]

where \(f_{\mathrm{H_2O}}\) is the fugacity of water (bar).

Reference
[1] T. S. Kuntzleman and C. Richards, "Another Method for Determining the Pressure inside an Intact Carbonated Beverage Can (or Bottle)," Journal of Chemical Education, 87(9), 2010 p. 993. doi:10.1021/ed100255g.

About

This simulation was created in the Department of Chemical and Biological Engineering at University of Colorado Boulder for LearnChemE.com by John L. Falconer using Claude AI. It is a JavaScript/HTML5 implementation of a Mathematica simulation by Rachael L. Baumann. It was prepared with financial support from the National Science Foundation (DUE 2336987 and 2336988). Address any questions or comments to LearnChemE@gmail.com.

P = 3.12 bar
concentration (mol/L)
liquid vapor fugacity (bar)
CO2
H2O