Search

A new method for predicting homogeneous bubble nucleation rates of pure compounds from vapor-liquid equilibrium (VLE) data is presented. It combines molecular dynamics simulation on the one side with density gradient theory using an equation of state (EOS) on the other. The new method is applied here to predict bubble nucleation rates in metastable liquid carbon dioxide (CO2). The molecular model of CO2 is taken from previous work of our group. PC-SAFT is used as an EOS. The consistency between the molecular model and the EOS is achieved by adjusting the PC-SAFT parameters to VLE data obtained from the molecular model. The influence parameter of density gradient theory is fitted to the surface tension of the molecular model. Massively parallel molecular dynamics simulations are performed close to the spinodal to compute bubble nucleation rates. From these simulations, the kinetic prefactor of the hybrid nucleation theory is estimated, whereas the nucleation barrier is calculated from density gradient theory. This enables the extrapolation of molecular simulation data to the whole metastable range including technically relevant densities. The results are tested against available experimental data and found to be in good agreement. The new method does not suffer from typical deficiencies of classical nucleation theory concerning the thermodynamic barrier at the spinodal and the bubble size dependence of surface tension, which is typically neglected in classical nucleation theory. In addition, the density in the center of critical bubbles and their surface tension is determined as a function of their radius. The usual linear Tolman correction to the capillarity approximation is found to be invalid. [ABSTRACT FROM AUTHOR]

Abstract: Vapour–liquid equilibria (VLE) and the influence of an inert carrier gas on homogeneous vapour to liquid nucleation are investigated by molecular simulation for quaternary mixtures of carbon dioxide, nitrogen, oxygen, and argon. Canonical ensemble molecular dynamics simulation using the Yasuoka–Matsumoto method is applied to nucleation in supersaturated vapours that contain more carbon dioxide than in the saturated state at the dew line. Established molecular models are employed that are known to accurately reproduce the VLE of the pure fluids as well as their binary and ternary mixtures. On the basis of these models, also the quaternary VLE properties of the bulk fluid are determined with the Grand Equilibrium method. Simulation results for the carrier gas influence on the nucleation rate are compared with the classical nucleation theory (CNT) considering the ‘pressure effect’ [Phys. Rev. Lett. 101: 125703 (2008)]. It is found that the presence of air as a carrier gas decreases the nucleation rate only slightly and, in particular, to a significantly lower extent than predicted by CNT. The nucleation rate of carbon dioxide is generally underestimated by CNT, leading to a deviation between one and two orders of magnitude for pure carbon dioxide in the vicinity of the spinodal line and up to three orders of magnitude in the presence of air as a carrier gas. Furthermore, CNT predicts a temperature dependence of the nucleation rate in the spinodal limit, which cannot be confirmed by molecular simulation. [Copyright &y& Elsevier]

An abstract of the article "Molecular Dynamics Simulation of Nucleation in CO2," by M. T. Horsch and colleagues is presented.

An abstract of the article "Molecular Dynamics Simulation of Nucleation in CO2," by M. T. Horsch and colleagues is presented.