Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Classical nucleation theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

A stress-sensitive precipitate nucleation model beyond classical nucleation theory

The dynamic evolution of precipitates and second phases dictates the strength and stability of most engineering alloys. By design, or as a consequence of thermo-mechanical aging, engineering metals and alloys often form precipitates of second phases when subjecting to diverse thermal and mechanical loads. Precipitation is governed by several factors, including the alloy’s composition, processing/operating temperature, and stresses — either as a result of external loads or from residual stresses. However, state-of-the-art models for precipitate nucleation (i.e., classical nucleation theory) typically lacks consistent method to capture the effects of externally applied and/or internal stresses on nucleation; thereby severely limiting the applicability of these models to complex materials systems and to representative loading scenarios. Here, in this work, we extend upon classical nucleation theory to account for the effect of stresses on precipitation kinetics and thermodynamics. This is achieved via the use of an Eshelbian micromechanics framework keeping track of (i) the stress build up resulting from second phase formation as a function of mechanical load and, (ii) the effects of dislocations on precipitate formation. This new model is applied to σ precipitate in Fe–Cr binary alloys and M 23 C 6 precipitate in 316H stainless steel (SS). Simulations demonstrate the important role of both the remotely applied loads and dislocation pile ups on precipitate nucleation.

36 MATERIALS SCIENCE

Data for The impact of aerosol mixing state on immersion freezing: Insights from classical nucleation theory and particle-resolved simulations

This dataset contains the values directly shown in the figures of the article "The impact of aerosol mixing state on immersion freezing: Insights from classical nucleation theory and particle-resolved simulations". This article is in preparation for submission to the journal Atmospheric Chemistry and Physics. The dataset consists of 12 NetCDF files processed from the raw output of the PartMC model. It does not include the theoretical values of frozen fraction, which can be computed using the equations provided in the paper. *New in V2: adding data for a newly included figure (INP_spectrum.nc), removing files that are no longer used in the revised manuscript figures (e.g., UNC_A_ratio=0.9_Dp=0.1.nc, UNC_A_ratio=0.9_Dp=10.0.nc, UNC_A_ratio=0.1_Dp=0.1.nc, and UNC_A_ratio=0.1_Dp=10.0.nc), and updating README.pdf accordingly.

Aerosol mixing state

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES

The mechanics of nucleation and growth and the surface tensions of active matter

Homogeneous nucleation, a textbook transition path for phase transitions, is typically understood on thermodynamic grounds through the prism of classical nucleation theory. However, recent studies have suggested the applicability of classical nucleation theory to systems far from equilibrium. Here, in this article, we formulate a purely mechanical perspective of homogeneous nucleation and growth, elucidating the criteria for the properties of a critical nucleus without appealing to equilibrium notions. Applying this theory to active fluids undergoing motility-induced phase separation, we find that nucleation proceeds in a qualitatively similar fashion to equilibrium systems, with concepts such as the Gibbs–Thomson effect and nucleation barriers remaining valid. We further demonstrate that the recovery of such concepts allows us to extend classical theories of nucleation rates and coarsening dynamics to active systems upon using the mechanically derived definitions of the nucleation barrier and surface tensions. Three distinct surface tensions—the mechanical, capillary, and Ostwald tensions—play a central role in our theory. While these three surface tensions are identical in equilibrium, our work highlights the distinctive role of each tension in the stability of active interfaces and the nucleation and growth of motility-induced phases.

Capillary waves

Deriving reliable nucleation rates from metadynamics simulations: Application to Yukawa fluids

In order to solidify the usefulness of metadynamics in studying nucleation of crystals from supercooled liquids, this work provides a specific procedure to calculate nucleation free energy barriers. After a pedagogical review of the important elements of classical nucleation theory and how metadynamics is used to find nucleation free energy barriers, we explain the benefits of local collective variables over more common global collective variables. We show how a metadynamics free energy barrier must be carefully post-processed so that classical nucleation theory can be applied to calculate nucleation rates. We apply our procedure to a Yukawa plasma and show that a particular physically motivated fit to metadynamics data reproduces low-temperature reference data, justifying the usefulness of metadynamics to predict nucleation rates.

42 ENGINEERING

How are Heterogeneous Nucleation Rate Observations Influenced by Instrument Resolution?

Experimental measurements of the heterogeneous nucleation rate rely on counting the number of nuclei with time. However, the size of a thermodynamically stable nucleus is often a few nanometers in diameter and is below the resolution of most (in situ) measurement techniques that provide a statistically valid sample. Due to the finite resolution of the instruments and analysis methods, it is challenging to capture the incipient nuclei and the subsequent evolution of nuclei density over time. In this work, we demonstrate the impact of instrument resolution on observed nuclei densities by comparing numerical modeling with experimental results. Further, to achieve this, we implemented heterogeneous nucleation within the pore-scale reactive transport modeling framework using classical nucleation theory (CNT). We compared the modeling results with nucleation rates measured using X-ray nanotomography (XnT) and evaluated how these impact the apparent values of the prefactor and interfacial energy based on CNT and the crystal growth rate. Specifically, we applied a resolution threshold (artificial resolution limit) in the model during nuclei counting to resemble an experimental resolution, ranging from 15 to 500 nm. The findings reveal that the instrument resolution significantly impacts the apparent prefactor and interfacial energy. Both apparent prefactor and interfacial energy decrease with a decrease in the instrument resolution. While deviation in the prefactor due to resolution is anticipated, those in the interfacial energy are unexpected. The approach described here allows one to correct apparent nucleation rates that depend on the instrument’s resolution to derive “intrinsic” CNT parameters for the prefactor and interfacial energy.

47 OTHER INSTRUMENTATION

Prognostic simulations of mixed-phase clouds with model AC-1D v1.0: the impact of aerosol types and freezing parameterizations on ice crystal budgets

Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.

Arctic clouds

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE

Quantifying and Modeling the Impact of Phase State on the Ice Nucleation Abilities of 2-Methyltetrols as a Key Component of Secondary Organic Aerosol Derived from Isoprene Epoxydiols

Organic aerosols (OAs) may serve as ice-nucleating particles (INPs), impacting the formation and properties of cirrus clouds when their phase state and viscosity are in the semisolid to glassy range. However, there is a lack of direct parameterization between aerosol viscosity and their ice nucleation capabilities. In this study, we experimentally measured the ice nucleation rate of 2-methyltetrols (2-MT) aerosols, a key component of isoprene-epoxydiol-derived secondary organic aerosols (IEPOX-SOA), at different viscosities. These results demonstrate that the phase state has a significant impact on the ice nucleation abilities of OA under typical cirrus cloud conditions, with the ice nucleation rate increasing by 2 to 3 orders of magnitude when the phase state changes from liquid to semisolid. An innovative parametric model based on classical nucleation theory was developed to directly quantify the impact of viscosity on the heterogeneous nucleation rate. This model accurately represents our laboratory measurement and can be implemented into climate models due to its simple, equation-based form. Based on data collected from the ACRIDICON-CHUVA field campaign, our model predicts that the INP concentration from IEPOX-SOA can reach the magnitude of 1 to tens per liter in the cirrus cloud region impacted by the Amazon rainforest, consistent with recent field observations and estimations. This novel parameterization framework can also be applied in regional and global climate models to further improve representations of cirrus cloud formation and associated climate impacts.

2-methyltetrol

Thermodynamic Control of Interface Directs MnO 2 Nucleation Chemistry for Dense and Conformal Electrodeposition

Manganese dioxide (MnO 2 ) is widely recognized as a promising material for high-energy-density energy storage systems due to its broad applicability and facile electrodeposition. However, achieving uniform, thin, and high-mass-loading MnO 2 coatings on high-surface-area electrodes remains a significant challenge. Conventional electrodeposition methods typically yield nonuniform, thick layers with poor conductivity and limited material utilization, restricting their practical use. Here, we uncover a thermodynamically engineered vanadyl/pervanadyl (VO 2+ /VO 2 + ) interface that fundamentally reshapes MnO 2 electrodeposition chemistry, enabling highly uniform and dense coatings. Here, combining in situ AFM measurement, Classical Nucleation Theory, and Johnson–Mehl–Avrami–Kolmogorov modeling, we show that this interface reduces early-stage detectable MnO 2 island size by 35-fold and shifts the MnO 2 growth from diffusion-limited to reaction-limited progressive nucleation. This thermodynamically controlled interface yields highly dense and conformal MnO 2 films with record-high mass loading of 241 mg cm –2 (1607 mg cm –3 ) on 3D-printed graphene aerogels, without compromising porosity or inducing thickness gradient. As a prototype demonstration, the resulting MnO 2 electrodes deliver record-setting volumetric performance in both capacitors (106 F cm –3 ) and Zn//MnO 2 pouch cells (162 mAh cm –3 ). Beyond energy storage, our findings demonstrate the significance of thermodynamic interface control in MnO 2 nucleation chemistry for achieving dense and uniform coatings on various substrates, with implications for electrocatalysis, semiconductor processing, and advanced materials manufacturing.

Batteries

Interpreting dynamic-compression experiments to uncover the time dependence of freezing: Application to gallium

Using pulsed-power magnetic field sources to compress gallium to gigapascal pressures on nanosecond timescales, we report here experiments on shockless dynamic compression of a liquid metal. Time-resolved velocimetry data reveal signatures of rapid freezing from a metastable liquid state, and we demonstrate that the kinetics of this nonequilibrium solidification can be accurately simulated with a computational modeling framework we have developed in previous studies, where classical nucleation theory is coupled with hydrodynamics. Notably, velocity traces in some of our experiments show evidence of a phase transition, while others do not, even though other types of evidence suggest that solidification may be occurring in all of them. We explain how predictions made by our models regarding the presence or absence of these phase-transition signatures motivated additional experiments that later confirmed the theoretical predictions. Our analysis shows that due to the rapid, quasi-isentropic nature of the loading path, our experiments were able to compress liquid gallium to metastable states that are undercooled below the equilibrium melt temperature by more than 300 K and exhibit pressures that approach five times the equilibrium melt pressure. The understanding gained in this study should form the basis for future dynamic-compression experiments aimed at interrogating melt curves at high pressures.

gallium

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics

The importance of accounting for the Tolman correction to surface tension for nucleation and growth modeling of Fe clusters

Gibbs free energies of clusters are required for predictive modeling of cluster growth during condensation of a cooling vapor. Here, we present a straightforward method of calculating free energies of cluster formation using the data from molecular dynamics (MD) simulations. We apply this method to iron clusters having from 2 to 100 atoms. The energies obtained are verified by comparing to an MD-simulated equilibrium cluster size distribution in a sub-saturated vapor. We show that these free energies differ significantly from those obtained with a commonly used spherical cluster approximation, which relies on a surface tension coefficient of a flat surface, as it is used in the classical nucleation theory (CNT). We show that the spherical cluster approximation in CNT can be improved by using a cluster-size-dependent Tolman correction for the surface tension. The Tolman length and effective surface tension values were derived for iron clusters, and they significantly differ from the commonly used experimentally measured values. This improved approximation does not account for geometric magic number effects responsible for spikes and troughs in densities of neighbor cluster sizes. Nonetheless, it allows to more accurately model cluster formation from a cooling vapor. It better reproduces the condensation timeline, overall shape of the cluster size distribution, average cluster size, and the distribution width. In contrast, using a constant surface tension coefficient (as done in CNT) resulted in incorrect condensation dynamics and cluster size distributions. The analytical expression for cluster nucleation rate from CNT was updated to account for the size-dependence of cluster surface tension.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science

Barrierless nucleation in glassy precursors drives zeolite formation

Zeolites are crystalline, microporous silicates widely used in catalysis and separations, yet the molecular mechanisms of their formation remain unresolved. Experiments indicate that hydrothermal synthesis of silica zeolites from clear solution proceeds through amorphous nanoaggregates that gradually develop zeolite order in an apparently continuous amorphous to crystal transformation. Here, we combine molecular simulations with advanced algorithms that identify zeolite order and computer vision to elucidate the pathway from clear solution to zeolite nanocrystal. We show that at conditions of hydrothermal synthesis of silica zeolites, the transformation of precursor aggregate into zeolite is not limited by nucleation barriers but by the slow dynamics of reorganization in the glassy precursor matrix. The negligible nucleation barriers result in spinodal-like crystallization that leads to a gradual formation of a mosaic of small crystallites that explain the seemingly continuous character of zeolite crystallization and the catalytic activity of X-ray amorphous, protozeolites and embryonic zeolites. We find that zeolite-like porosity and short-range order emerge early within glassy precursors, well before crystallinity is detected in transmission electron microscopy (TEM) images or X-ray diffraction. The nanoaggregate’s temperature-size phase diagram reveals a convergence of the zeolite–amorphous equilibrium and maximum crystallization rate at ~3 nm diameter nanoparticle diameters and ~200 °C. This convergence signals the termination of the first-order amorphous-to-zeolite transition. Our results provide a unifying framework for understanding nucleation of silica zeolites from solution and suggest that barrierless nucleation may govern the formation of other nanoparticle systems, including minerals and oxides synthesized far below their bulk melting points.

36 MATERIALS SCIENCE

Revealing and Engineering Assembly Pathways of 3D DNA Origami Crystals

Recent developments in nanomaterial self-assembly demonstrate the capability to create tailored nanostructures by engineering both the binding coordination and specificity of interactions between material subunits. DNA origami frames allow for the design and fabrication of a broad variety of ordered 3D nanoscale architectures through self-assembly, facilitated by frame-to-frame bonds with designable strength and specificity. While the bond design is critical to lattice formation, the assembly process itself is often dependent on a thermal pathway. Highly ordered nanoscale frameworks, assembled from DNA frames, are predominantly crystallized through thermal annealing pathways that typically follow a “slow” cooling approach, with experiments on the time scale of days yielding DNA origami crystals in the range of 1−10 μm. This extended assembly time scale hinders the study of crystal formation pathways, necessitating a deeper understanding of factors governing successful annealing. Lack of insight into time scale also presents a practical limitation for material fabrication. Here, we investigate key factors affecting lattice assembly pathways and demonstrate that precise engineering of assembly conditions greatly reduces assembly times by up to nearly 2 orders of magnitude. We evaluate the nucleation and growth of crystals via optical and electron microscopy, and small-angle X-ray scattering techniques, mapping the time−temperature-transformation of superlattices from the melt through single-crystal optical tracking. The results show that origami frame assembly can be described by classical nucleation and growth theory, which can, in turn, be used to prescribe the growth of the crystals. Lastly, these findings are applied to demonstrate thermal pathway-dependent assembly, forming distinct assemblies based on different thermal annealing profiles.

36 MATERIALS SCIENCE