Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

Modeling and Simulation of Pore Formation in a Bainitic Steel During Creep

In the field of power engineering, where materials are subjected to high pressures at elevated temperatures for many decades, creep-resistant steels are put to work. Their service life is still, however, finite, as the many changes in their microstructure can merely be mitigated and not avoided. Creep cavitation is one of those changes and, in many cases, ultimately causes failure by rupture. In this work, a model is proposed to simulate the nucleation and growth of cavities during creep. This exclusively physics-based model uses modified forms of Classical Nucleation Theory and the Onsager Extremum Principle in a newly developed Kampmann–Wagner framework. The model is validated on P23 steel which underwent creep rupture experiments at 600 °C and stresses of 50, 70, 80, 90 and 100 MPa for creep times up to 46000 hours. The model predicts qualitatively the shape and prevalence of cavities at different sites in the microstructure, and quantitatively the number density, size of cavities and their phase fraction contributing to a reduction in density. Finally, we find good agreement between the simulation and the experimental results especially at low stresses and longer creep times.

Nucleation Parameters of SPC/E And TIP4P/2005 Water Vapor Measured In NPT Molecular Dynamics Simulations

Nucleation rates for droplet formation in water vapor are measured in molecular dynamics simulations of SPC/E and TIP4P/2005 water by monitoring individual nucleation events. The nucleation process is simulated in the NPT ensemble to evaluate the steady-state nucleation rate in accordance with the assumptions of classical nucleation theory (CNT). Nucleation rates measured between 300 K and 425 K for the SPC/E model, and between 325 K and 475 K for the TIP4P/2005 model, agree with the CNT predictions roughly within the standard deviation of the MD measurements of the nucleation rates.

Application of the Classical Nucleation Theory for Quasi-real Systems Differ in Dipole Moment Value – the Role of Diffusion

In this paper, we examine the crystallization tendency for two quasi-real systems, which differ exclusively in the dipole moment's value. The main advantage of the studied system is the fact that despite that their structures are entirely identical, they exhibit different physical properties. Hence, the results obtained for one of the proposed model systems cannot be scaled to reproduce the results for another corresponding system, as it can be done for simple model systems, where structural differences are modeled by the different parameters of the intermolecular interactions. Our results show that both examined systems exhibit similar stability behavior below the melting temperature. This finding is contrary to the classical nucleation theory predictions, which differ significantly for them. On the basis of the performed studies, we suggest that a kinetic aspect of the classical nucleation theory seems to be a reason for reported discrepancies.

Dynamic Surface Tension Enhances the Stability of Nanobubbles in Xylem Sap

Air seeded nanobubbles have recently been observed within tree sap under negative pressure. They are stabilized by an as yet unidentified process, although some embolize their vessels in extreme circumstances. Current literature suggests that a varying surface tension helps bubbles survive, but few direct measurements of this quantity have been made. Here, we present calculations of dynamic surface tension for two biologically relevant lipids using molecular dynamics simulations. We find that glycolipid monolayers resist expansion proportionally to the rate of expansion. Their surface tension increases with the tension applied, in a similar way to the viscosity of a non-Newtonian fluid. In contrast, a prototypical phospholipid was equally resistant to all applied tensions, suggesting that the fate of a given nanobubble is dependent on its surface composition. By incorporating our results into a Classical Nucleation Theory (CNT) framework, we predict nanobubble stability with respect to embolism. We find that the metastable radius of glycolipid coated nanobubbles is approximately 35 nm, and that embolism is in this case unlikely when the external pressure is less negative than –1.5 MPa.

A lab-on-a-chip approach integrating in-situ characterization and reactive transport modelling diagnostics to unravel (Ba,Sr)SO4 oscillatory zoning

The co-precipitation of sulphate minerals such as celestine and barite is widely studied because their formation is ubiquitous in natural and anthropogenic systems. Co-precipitation in porous media results in crystallization of solid solutions yielding characteristics such as oscillatory zoning that are rarely observed in bulk solution or in batch experiments. In the past, the precipitation of compositionally-zoned (Ba,Sr)SO4 crystals was observed post-mortem in macroscopic silica gel counter-diffusion experiments. Their formation was originally explained by the difference in the solubility products of the end-members combined with diffusion-limited transport of solutes to the mineral-fluid interface, while a later study favored the idea of kinetically controlled reactions. With recent advances combining in-operando microfluidic experiments and reactive transport modelling, it is now possible to verify hypotheses on the driving forces of transport-coupled geochemical processes. We developed a “lab on a chip” experiment that enabled the systematic study of the nucleation and growth of oscillatory-zoned (Ba,Sr)SO4 crystals in a microfluidic reactor. The compositions of the solid solutions were determined by in-situ Raman spectroscopy. Our investigation shows (1) that the composition of the nucleating phases can be approximated using classical nucleation theory, (2) that the oscillatory zoning is not solely controlled by the limited diffusional transport of solutes, and (3) that nucleation kinetics plays a major role in the switch between different stoichiometric compositions. The zoning phenomena is governed by the complex interplay between the diffusion of reactants and the crystallization kinetics as well as other factors, e.g. surface tension and lattice mismatch.

Freezing of few nanometers water droplets

Water-ice transformation of few nm nanodroplets plays a critical role in nature including climate change, microphysics of clouds, survival mechanism of animals in cold environments, and a broad spectrum of technologies. In most of these scenarios, water-ice transformation occurs in a heterogenous mode where nanodroplets are in contact with another medium. Despite computational efforts, experimental probing of this transformation at few nm scales remains unresolved. Here, we report direct probing of water-ice transformation down to 2 nm scale and the length-scale dependence of transformation temperature through two independent metrologies. The transformation temperature shows a sharp length dependence in nanodroplets smaller than 10 nm and for 2 nm droplet, this temperature falls below the homogenous bulk nucleation limit. Contrary to nucleation on curved rigid solid surfaces, ice formation on soft interfaces (omnipresent in nature) can deform the interface leading to suppression of ice nucleation. For soft interfaces, ice nucleation temperature depends on surface modulus. Considering the interfacial deformation, the findings are in good agreement with predictions of classical nucleation theory. This understanding contributes to a greater knowledge of natural phenomena and rational design of anti-icing systems for aviation, wind energy and infrastructures and even cryopreservation systems.

Water cavitation from ambient to high temperatures

Predicting cavitation has proved a formidable task, particularly for water. Despite the experimental difficulty of controlling the sample purity, there is nowadays substantial consensus on the remarkable tensile strength of water, on the order of −120 MPa at ambient conditions. Recent progress significantly advanced our predictive capability which, however, still considerably depends on elaborate fitting procedures based on the input of external data. Here a self-contained model is discussed which is shown able to accurately reproduce cavitation data for water over the most extended range of temperatures for which accurate experiments are available. The computations are based on a diffuse interface model which, as only inputs, requires a reliable equation of state for the bulk free energy and the interfacial tension. A rare event technique, namely the string method, is used to evaluate the free-energy barrier as the base for determining the nucleation rate and the cavitation pressure. The data allow discussing the role of the Tolman length in determining the nucleation barrier, confirming that, when the size of the cavitation nuclei exceed the thickness of the interfacial layer, the Tolman correction effectively improves the predictions of the plain Classical Nucleation Theory.

Formation of ice particles through nucleation in the mesosphere

Observations of polar mesospheric clouds have revealed the presence of solid ice particles in the upper mesosphere at high latitudes; however, their formation mechanism remains uncertain. In this study, we investigated the formation process of ice particles through nucleation from small amounts of water vapor at low temperatures. Previous studies that used classical nucleation theory have shown that amorphous solid water particles can nucleate homogeneously at conditions that are present in the mesosphere. However, the rate predictions for water in classical nucleation theory disagree with experimental measurements by several orders of magnitude. We adopted a semi-phenomenological model for the nucleation process, which corrects the evaluation of the molecular cluster formation energy using the second virial coefficient, which agrees with both experiments and molecular dynamics simulations. To calculate the nucleation process, we applied atmospheric conditions for the temperature, pressure, numerical density of dust grains, and cooling rate. The results indicate that homogeneous water nucleation is extremely unlikely to occur in the mesosphere, while heterogeneous nucleation occurs effectively. Dust grains generated by meteor ablation can serve as nuclei for heterogeneous nucleation. We also showed that the ice can form directly in a crystalline state, rather than an amorphous state.