Gas-liquid nucleation at large metastability: unusual features and a new formalism

Nucleation at large metastability is still largely an unsolved problem, even though it is a problem of tremendous current interest, with wide-ranging practical value, from atmospheric research to materials science. It is now well accepted that the classical nucleation theory (CNT) fails to provide a qualitative picture and gives incorrect quantitative values for such quantities as activation-free energy barrier and supersaturation dependence of nucleation rate, especially at large metastability. In this paper, we present an alternative formalism to treat nucleation at large supersaturation by introducing an extended set of order parameters in terms of the kth largest liquid-like clusters, where k = 1 is the largest cluster in the system, k = 2 is the second largest cluster and so on. At low supersaturation, the size of the largest liquid-like cluster acts as a suitable order parameter. At large supersaturation, the free energy barrier for the largest liquid-like cluster disappears. We identify this supersaturation as the one at the onset of kinetic spinodal. The kinetic spinodal is system-size-dependent. Beyond kinetic spinodal many clusters grow simultaneously and competitively and hence the nucleation and growth become collective. In order to describe collective growth, we need to consider the full set of order parameters. We derive an analytic expression for the free energy of formation of the kth largest cluster. The expression predicts that, at large metastability (beyond kinetic spinodal), the barrier of growth for several largest liquid-like clusters disappears, and all these clusters grow simultaneously. The approach to the critical size occurs by barrierless diffusion in the cluster size space. The expression for the rate of barrier crossing predicts weaker supersaturation dependence than what is predicted by CNT at large metastability. Such a crossover behavior has indeed been observed in recent experiments (but eluded an explanation till now). In order to understand the large numerical discrepancy between simulation predictions and experimental results, we carried out a study of the dependence on the range of intermolecular interactions of both the surface tension of an equilibrium planar gas-liquid interface and the free energy barrier of nucleation. Both are found to depend significantly on the range of interaction for the Lennard-Jones potential, both in two and three dimensions. The value of surface tension and also the free energy difference between the gas and the liquid phase increase significantly and converge only when the range of interaction is extended beyond 6-7 molecular diameters. We find, with the full range of interaction potential, that the surface tension shows only a weak dependence on supersaturation, so the reason for the breakdown of CNT (with simulated values of surface tension and free energy gap) cannot be attributed to the supersaturation dependence of surface tension. This remains an unsettled issue at present because of the use of the value of surface tension obtained at coexistence.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

Thermodynamics of oxygen in liquid copper-tin alloys

The activity coefficients of oxygen in copper-tin alloys at 1 1 00°C have been measured by two different equilibrium methodsthe cell Pt, Ni + NiO I ZrOz solid electrolyte I O[Cu + Sn], cermet. Pt and the equilibrium between Cu + Sn alloys and SnO + SiO, slags established via SnO vapour. The results from both types of measurement confirm the work of Block and co-workers and show that other data are in error. The deoxidation equilibria for Sn in liquid copper, with solid SnO, as deoxidation product, have been evaluated at temperatures of interest in copper smelting.

Solubility and activity of oxygen in liquid gallium and gallium copper alloys

The solubility of oxygen in liquid gallium in the temperature range 775 –1125 °C and in liquid gallium-copper alloys at 1100 °C, in equilibrium with β-Ga2O3, has been measured by an isopiestic equilibrium technique. The solubility of oxygen in pure gallium is given by the equation log (at.% O) = −7380/T + 4.264 (±0.03) Using recently measured values for the standard free energy of formation of β-Ga2O3 and assuming that oxygen obeys Sievert's law up to the saturation limit, the standard free energy of solution of oxygen in liquid gallium may be calculated : View the MathML sourceΔ°298 = −52 680 + 6.53T (±200) cal where the standard state for dissolved oxygen is an infinitely dilute solution in which the activity is equal to atomic per cent. The effect of copper on the activity of oxygen dissolved in liquid gallium is found to be in good agreement with that predicted by a recent quasichemical model in which it was assumed that each oxygen is interstitially coordinated to four metal atoms and that the nearest neighbour metal atoms lose approximately half their metallic cohesive energies.

Solubility and activity of oxygen in liquid indium and copper indium alloys

The solubility of oxygen in liquid indium in the temperature range 650–820 °C and in liquid copper-indium alloys at 1100 °C in equilibrium with indium sesquioxide has been measured by a phase equilibration technique. The solubility of oxygen in pure indium is given by the relation log(at.% O) = −4726/T + 3.73 (±0.08) Using the recently measured values for the standard free energy of formation of In2O3 and assuming that oxygen obeys Sievert's law up to saturation, the standard free energy of solution of molecular oxygen in liquid indium is calculated as View the MathML sourceΔG°= −51 440 + 8.07 T (±500) cal where the standard state for dissolved oxygen is an infinitely dilute solution in which activity is equal to atomic per cent. The effect of indium additions on the activity coefficient of oxygen dissolved in liquid copper was measured by a solid oxide galvanic cell. The interaction parameter ϵ0In is given by View the MathML source The experimentally determined variation of the activity coefficient of oxygen in dilute solution in Cu-In alloys is in fair agreement with that predicted by a quasichemical model in which each oxygen atom is assumed to be interstitially coordinated to four metal atoms and the nearest neighbour metal atoms are assumed to lose approximately half their metallic cohesive energies.

Solubility and activity of oxygen in liquid germanium and germanium copper alloys

The solubility of oxygen in liquid germanium in the temperature range 1233 to 1397 K, and in liquid germanium-copper alloys at 1373 K, in equilibrium with GeO2 has been measured by the phase equilibration technique. The solubility of oxygen in pure germanium is given by the relation R470 log(at. pct 0)=-6470/T+4.24 (±0.07). The standard free energy of solution of oxygen in liquid germanium is calculated from the saturation solubility, and recently measured values for the free energy of formation of GeO2, assuming that oxygen obeys Sievert’s law up to the saturation limit. For the reaction, 1/2 O2(g)→ OGe ΔG° =-39,000 + 3.21T (±500) cal = -163,200 + 13.43T (±2100) J. where the standard state for dissolved oxygen is that which makes the value of activity equal to the concentration (in at. pct), in the limit, as concentration approaches zero. The effect of copper on the activity of oxygen dissolved in liquid germanium is found to be in good agreement with that predicted by a quasichemical model in which each oxygen was assumed to be bonded to four metal atoms and the nearest neighbor metal atoms to an oxygen atom are assumed to lose approximately half of their metallic bonds.

Vapor pressure and dissociation energy of (In2O)

The vapor pressure of pure liquid indium, and the sum of pressures of (In) and (In2O) species over the condensed phase mixture {In} + , contained in a silica vessel, have been measured by Knudsen effusion and Langmuir free vaporization methods in the temperatue range 600 to 950°C. Mass spectrometric studies reported in the literature show that (In) and (In2O) are the important species in the vapor phase over the {In} + ; mixture. The vapor pressure of (In2O) corresponding to the reaction, deduced from the present measurements is given by the equation, The “apparent evaporation coefficient” for the condensed phase mixture is approximately 0.8. The energy for the dissociation (In2O) molecule into atoms calculated from the above equation is D°0 = 180.0 (± 1.0) kcal mol−1.

Solubility and activity of oxygen in liquid antimony

The solubility limit of oxygen in liquid antimony has been measured by a novel isopiestic technique in the temperature range 995--1175 deg K. The results can be expressed by the equation log c = --5500/T + 3.754 ( plus/minus 0.04) with c in at.% O and T in deg K. The oxygen potential over Sb + O alloys equilibrium with Sb2O3 has been measured by a solid state cell using a fully stabilized CaO.ZrO2 electrolyte. The cell was designed to contain the Sb + Sb2O3 mixture in a closed volume, that the vaporization of the oxide can be minimized and true equilibrium attained. The Gibbs free energy of the reaction 2 Sb(s) + 3/2 O2 = Sb2O3(s) is Delta G deg = --719560 + 274.51 T( plus/minus 500) and Sb(l) + 3/2 O2 = Sb2O3(l), Delta G deg = --704711 + ( plus/minus 500) ( Delta G deg in J/mole, T in deg K). The combination of these results with Sieverts' law yields the standard free energy of solution of oxygen in liquid antimony according to the reaction 1/2 O2 = \O\Sb,at.% as Delta G deg = --129620 + 14.23 T ( plus/minus 950). The standard enthalopy and entropy of the solution of oxygen in Sb are compared with values for other metal- oyxgen systems, and with the standard enthalpies of formation of corresponding oxides. The resulting correlations permit the estimation of the standard free energy of solution of oxygen in pure metals for which experimental information is lacking. 24 ref.--AA

Solubility and activity of oxygen in liquid manganese

The oxygen concentration of liquid manganese in equilibrium with MnAl2+2xO4+3x and α−Al2O3 has been determined in the temperature range 1520 to 1875 K. The oxygen content of quenched samples, wrapped in oxygen-free nickel foil, was determined by an inert gas fusion technique. The results are combined with accurate data now available on the Gibbs energies of formation of MnO and Al2O3−saturated MnAl2+2xO4+3x to derive the oxygen content of liquid manganese in equilibrium with MnO and the Gibbs energy of solution of diatomic oxygen gas in liquid manganese. The enthalpy and entropy of solution of oxygen in manganese are compared with similar data on other metal-oxygen systems.

Current views on microscopic thermodynamics of liquid alloys

An attempt has been made to review current information on the microscopic thermodynamics of liquid alloys. For complex alloys, and for alloys of simple metals with strong "compound-forming" tendencies, the fluctuation approach developed by Bhatia and his co-workers provides a useful link between the fluctuation in concentration and number density of atoms in the mixture on the one hand, and macroscopic thermodynamic properties on the other. Some selected examples of the application of structural data of liquid alloys to estimating macroscopic thermodynamic properties such as the Gibbs free energy of mixing, coupled with the fluctuation approach are given. The relevant thermodynamic quantities such as vapor pressure and entropy are also discussed, to facilitate the understanding of the present status of the fundamental and powerful links between macroscopic and microscopic (atomic scale) structure of liquid alloys (Mg--Sn, Li--Pb, Hg--K). 63 ref.--AA

Measurement and modeling of Alloy-Spinel-Corundum equilibrium in the Ni-Mn-Al-0 system at 1873 K

The oxygen content of liquid Ni-Mn alloy equilibrated with spinel solid solution, (Ni,Mn)O. (1 +x)A12O3, and α-Al2O3 has been measured by suction sampling and inert gas fusion analysis. The corresponding oxygen potential of the three-phase system has been determined with a solid state cell incorporating (Y2O3)ThO2 as the solid electrolyte and Cr + Cr2O3 as the reference electrode. The equilibrium composition of the spinel phase formed at the interface of the alloy and alumina crucible was obtained using EPMA. The experimental data are compared with a thermodynamic model based on the free energies of formation of end-member spinels, free energy of solution of oxygen in liquid nickel, interaction parameters, and the activities in liquid Ni-Mn alloy and spinel solid solution. Mixing properties of the spinel solid solution are derived from a cation distribution model. The computational results agree with the experimental data on oxygen concentration, potential, and composition of the spinel phase.

Relation between Local Structure and Glass Forming Ability of Liquid Alloys

An attempt has been made to describe the glass forming ability (GFA) of liquid alloys, using the concepts of the short range order (SRO) and middle range order (MRO) characterizing the liquid structure.A new approach to obtain good GFA of liquid alloys is based on the following four main factors: (1) formation of new SRO and competitive correlation with two or more kinds of SROs for crystallization, (2) stabilization of dense random packing by interaction between different types of SRO, (3) formation of stable cluster (SC) or middle range order (MRO) by harmonious coupling of SROs, and (4) difference between SRO characterizing the liquid structure and the near-neighbor environment in the corresponding equilibrium crystalline phases. The atomic volume mismatch estimated from the cube of the atomic radius was found to be a close relation with the minimum solute concentration for glass formation. This empirical guideline enables us to provide the optimum solute concentration for good GFA in some ternary alloys. Model structures, denoted by Bernal type and the Chemical Order type, were again tested in the novel description for the glass structure as a function of solute concentration. We illustrated the related energetics of the completion between crystal embryo and different types of SRO. Recent systematic measurements also provide that thermal diffusivity of alloys in the liquid state may be a good indicator of their GFA.

Effect of liquid precursor pyrolysis on phase selection in the MgO-MgAl2O4 system

Aqueous solutions of Al and Mg nitrates have been spray pyrolysed at 673 K to synthesize powders with compositions varying between MgO and MgAl2O4. This has been carried out with the aim of studying phase selection and phase evolution in this system. The powders have been subsequently heat treated and the sequence of phases characterised by X-ray diffraction and transmission electron microscopy. Metastable extensions of the different phase fields have been calculated based on functions which predict the equilibrium phase diagram accurately. The appearance of phases is closely related to the temperature and to the non-stoichiometry in different compositional ranges of the system. The sequence of phase evolution has been correlated to the thermodynamics of nucleation in the system.

Vaporization and collision modeling of liquid fuel sprays in a co-axial fuel and air pre-mixer

Droplet collision occurs frequently in regions where the droplet number density is high. Even for Lean Premixed and Pre-vaporized (LPP) liquid sprays, the collision effects can be very high on the droplet size distributions, which will in turn affect the droplet vaporization process. Hence, in conjunction with vaporization modeling, collision modeling for such spray systems is also essential. The standard O'Rourke's collision model, usually implemented in CFD codes, tends to generate unphysical numerical artifact when simulations are performed on Cartesian grid and the results are not grid independent. Thus, a new collision modeling approach based on no-time-counter method (NTC) proposed by Schmidt and Rutland is implemented to replace O'Rourke's collision algorithm to solve a spray injection problem in a cylindrical coflow premixer. The so called ``four-leaf clover'' numerical artifacts are eliminated by the new collision algorithm and results from a diesel spray show very good grid independence. Next, the dispersion and vaporization processes for liquid fuel sprays are simulated in a coflow premixer. Two liquid fuels under investigation are jet-A and Rapeseed Methyl Esters (RME). Results show very good grid independence in terms of SMD distribution, droplet number distribution and fuel vapor mass flow rate. A baseline test is first established with a spray cone angle of 90 degrees and injection velocity of 3 m/s and jet-A achieves much better vaporization performance than RME due to its higher vapor pressure. To improve the vaporization performance for both fuels, a series of simulations have been done at several different combinations of spray cone angle and injection velocity. At relatively low spray cone angle and injection velocity, the collision effect on the average droplet size and the vaporization performance are very high due to relatively high coalescence rate induced by droplet collisions. Thus, at higher spray cone angle and injection velocity, the results expectedly show improvement in fuel vaporization performance since smaller droplet has a higher vaporization rate. The vaporization performance and the level of homogeneity of fuel-air mixture can be significantly improved when the dispersion level is high, which can be achieved by increasing the spray cone angle and injection velocity. (C) 2012 Elsevier Ltd. All rights reserved.

Experimental analysis of shape deformation of evaporating droplet using Legendre polynomials

Experiments involving heating of liquid droplets which are acoustically levitated, reveal specific modes of oscillations. For a given radiation flux, certain fluid droplets undergo distortion leading to catastrophic bag type breakup. The voltage of the acoustic levitator has been kept constant to operate at a nominal acoustic pressure intensity, throughout the experiments. Thus the droplet shape instabilities are primarily a consequence of droplet heating through vapor pressure, surface tension and viscosity. A novel approach is used by employing Legendre polynomials for the mode shape approximation to describe the thermally induced instabilities. The two dominant Legendre modes essentially reflect (a) the droplet size reduction due to evaporation, and (b) the deformation around the equilibrium shape. Dissipation and inter-coupling of modal energy lead to stable droplet shape while accumulation of the same ultimately results in droplet breakup. (C) 2013 Elsevier B.V. All rights reserved.