Gradient theory computation of the radius-dependent surface tension and nucleation rate for n-nonane clusters

Three-dimensional phase field modeling of inhomogeneous gas-liquid systems using the PeTS equation of state

Recently, an equation of state (EoS) for the Lennard-Jones truncated and shifted (LJTS) fluid has become available. As it describes metastable and unstable states well, it is suited for predicting density profiles in vapor-liquid interfaces in combination with density gradient theory (DGT). DGT is usually applied to describe interfaces in Cartesian one-dimensional scenarios. In the present work, the perturbed LJ truncated and shifted (PeTS) EoS is implemented into a three-dimensional phase field (PF) model which can be used for studying inhomogeneous gas-liquid systems in a more general way. The results are compared with the results from molecular dynamics simulations for the LJTS fluid that are carried out in the present work and good agreement is observed. The PF model can therefore be used to overcome the scale limit of molecular simulations. A finite element approach is applied for the implementation of the PF model. This requires the first and second derivatives of the PeTS EoS which are calculated using hyper-dual numbers. Several tests and examples of applications of the new PeTS PF model are discussed.

Predictions of homogeneous nucleation rates for n-alkanes accounting for the diffuse phase interface and capillary waves

Homogeneous droplet nucleation has been studied for almost a century but has not yet been fully understood. In this work, we used the density gradient theory (DGT) and considered the influence of capillary waves (CWs) on the predicted size-dependent surface tensions and nucleation rates for selected n-alkanes. The DGT model was completed by an equation of state (EoS) based on the perturbed-chain statistical associating fluid theory and compared to the classical nucleation theory and the Peng-Robinson EoS. It was found that the critical clusters are practically free of CWs because they are so small that even the smallest wavelengths of CWs do not fit into their finite dimensions. The CWs contribute to the entropy of the system and thus decrease the surface tension. A correction for the effect of CWs on the surface tension is presented. The effect of the different EoSs is relatively small because by a fortuitous coincidence their predictions are similar in the relevant range of critical cluster sizes. The difference of the DGT predictions to the classical nucleation theory computations is important but not decisive. Of the effects investigated, the most pronounced is the suppression of CWs which causes a sizable decrease of the predicted nucleation rates. The major difference between experimental nucleation rate data and theoretical predictions remains in the temperature dependence. For normal alkanes, this discrepancy is much stronger than observed, e.g., for water. Theoretical corrections developed here have a minor influence on the temperature dependency. We provide empirical equations correcting the predicted nucleation rates to values comparable with experiments.

Curvature Dependence of the Liquid-Vapor Surface Tension beyond the Tolman Approximation

Surface tension is a macroscopic manifestation of the cohesion of matter, and its value σ_{∞} is readily measured for a flat liquid-vapor interface. For interfaces with a small radius of curvature R, the surface tension might differ from σ_{∞}. The Tolman equation, σ(R)=σ_{∞}/(1+2δ/R), with δ a constant length, is commonly used to describe nanoscale phenomena such as nucleation. Here we report experiments on nucleation of bubbles in ethanol and n-heptane, and their analysis in combination with their counterparts for the nucleation of droplets in supersaturated vapors, and with water data. We show that neither a constant surface tension nor the Tolman equation can consistently describe the data. We also investigate a model including 1/R and 1/R^{2} terms in σ(R). We describe a general procedure to obtain the coefficients of these terms from detailed nucleation experiments. This work explains the conflicting values obtained for the Tolman length in previous analyses, and suggests directions for future work.

Extending the nonequilibrium square-gradient model with temperature-dependent influence parameters

Nonequilibrium interface phenomena play a key role in crystallization, hydrate formation, pipeline depressurization, and a multitude of other examples. Square gradient theory extended to the nonequilibrium domain is a powerful tool for understanding these processes. The theory gives an accurate prediction of surface tension at equilibrium, only with temperature-dependent influence parameters. We extend in this work the nonequilibrium square gradient model to have temperature-dependent influence parameters. The extension leads to thermodynamic quantities which depend on temperature gradients. Remarkably the Gibbs relation proposed in earlier work is still valid. Also for the extended framework, the "Gibbs surface" described by excess variables is found to be in local equilibrium. The temperature-dependent influence parameters give significantly different interface resistivities (∼9%-50%), due to changed density gradients and additional terms in the enthalpy. The presented framework facilitates a more accurate description of transport across interfaces with square gradient theory.

Theoretical study of vapor-liquid homogeneous nucleation using stability analysis of a macroscopic phase

Stability analysis is generally used to verify that the solution to phase equilibrium calculations corresponds to a stable state (minimum of the free energy). In this work, tangent plane distance analysis for stability of macroscopic mixtures is also used for analyzing the nucleation process, reconciling thus this analysis with classical nucleation theories. In the context of the revised nucleation theory, the driving force and the nucleation work are expressed as a function of the Lagrange multiplier corresponding to the mole fraction constraint from the minimization problem of stability analysis. Using a van der Waals fluid applied to a ternary mixture, Lagrange multiplier properties are illustrated. In particular, it is shown how the Lagrange multiplier value is equal to one on the binodal and spinodal curves at the same time as the driving force of nucleation vanishes on these curves. Finally, it is shown that, on the spinodal curve, the nucleation work from the revised and generalized nucleation theories are characterized by two different local minima from stability analysis, irrespective of any interfacial tension models.