Search NASASearch

SEARCH · Search NASA

Results for “gas effluent data analysis”

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.

1,302 records · Page 68

ClassNMSW- a real-time classification approach for non-recycled municipal solid waste using hyperspectral imaging

Real-time classification of non-recycled municipal solid waste (NMSW) is essential for efficient valorization. This study introduces ClassNMSW, a comprehensive framework for classifying 22 NMSW subclasses under industrial constraints by using hyperspectral imaging (HSI). A primary innovation of this work is the development of a variance-controlled spectral extraction algorithm. Unlike traditional methods that rely on simple averaging, this approach systematically investigates the extent of pixel extraction to minimize the loss of critical chemical information while maximizing data reduction thus ensuring high spectral fidelity with low computational cost. The approach developed in this work integrates automated, computer-vision-based background removal, eliminating the need for the manual thresholding common in current literature. To resolve ambiguities among chemically similar subclasses, a tiered classification and multi-camera fusion strategy (NIR17 and NIR22) is implemented. Results demonstrate that ClassNMSW achieves an object-wise weighted accuracy of 98.70% for single-sensor configurations and 100% under sensor fusion. A novel rolling-window strategy satisfies desired end-to-end latency of <2 s, satisfying the strict deterministic requirements of high-speed industrial sorting environments. The ClassNMSW framework provides a scalable foundation for advancing circularity and resource recovery in large-scale waste valorization operations.

99 - GENERAL AND MISCELLANEOUS

Enabling Reliable, Fault-Tolerant Autonomous Lunar Habitats with High-Performance Spaceflight Computing

The lunar surface presents unfavorable constraints and harsh living conditions. To address these challenges, autonomous habitats will require complex integrated systems that combine advanced software, high-performance hardware, and cutting-edge sensors to ensure sustainability, safety, and operational efficiency. Consequently, maintaining a sustainable presence on the Moon requires reliable infrastructure and efficient development, precise monitoring, and utilization of resources within a lunar installation. These elements are essential not only to ensure that lunar settlement can be long-term, self-sustaining, and resource-efficient, but also to serve as a foundation for future missions and eventual human habitation on Mars. Humans are not native to the Moon; therefore, our survival and ability to thrive will depend on autonomous systems that can foster safety and resilience through high-availability architectures, graceful degradation, and highly fault-tolerant spaceflight hardware capable of continuing operation during failures. This requires advanced human-rated distributed systems architectures with specialized electronics, scalable capabilities, and an integrated design approach. Unlike current practices focused on short-term missions and regularly maintained components, permanent lunar compute systems must be designed for extended operations beyond mission durations. This paper explores the necessity of transitioning toward fault- tolerant, highly autonomous hardware systems designed for multi-year missions. It also identifies critical subsystems that require high levels of autonomy, supported by radiation-hardened processors and extreme thermal loads, which are essential to mitigate long-term degradation and ensure sustainable lunar habitation. Finally, the paper aligns with NASA’s identified Civil Space Shortfalls, particularly in high-performance onboard computing, advanced data acquisition, extreme-environment avionics, radiation monitoring and countermeasures, and autonomous health management. It proposes NASA’s new High-Performance Spaceflight Computing (HPSC) processor as a turnkey solution, delivering 100 times the performance-per-watt of legacy rad-hard CPUs and enabling onboard AI, edge computing, and fault-tolerant features essential for sustained lunar autonomy and beyond.

Sarkis S Mikaelian

Semi-Dynamic Leach Testing of Densified Silicon-based Iodine Waste Forms

Iodine waste forms (IWF) require a conceptual corrosion release model (CCRM) to provide iodine (I) release rates for performance modeling nuclear waste disposal repositories. To develop a CCRM, an understanding of the corrosion mechanisms of the IWF are required along with data from consistent test methods to parameterize the model. The present study has advanced both areas by providing and assessing the corrosion resistance of IWF types based on their processing history in minor variations of semi-dynamic leach tests. The test included a series of semi-dynamic leach tests using monolithic IWFs in deionized water (leachant). Several experiments were conducted under alternate test conditions with changes to temperature, leachant replacement, leachant pH, leachant volume, masking, and surface finish to elucidate if varying these conditions impacted IWF corrosion behavior. Tests were conducted on two classes of IWFs: (1) I-bearing silver-mordenite (AgZ) materials processed by hot isostatic pressing (HIP) at different temperatures, pressures, sizes, and times; and (2) I-bearing silver-functionalized silica aerogels (SFA) processed by either HIP or spark plasma sintering (SPS). The corrosion susceptibility of AgZ samples was influenced by HIP temperature and pressure. The SPS SFAs retained I far better than HIP SFAs. Additional findings in this study include: (1) The iodine dissolution rate decreased with decreasing temperature, (2) a common ion effect may occur and slow dissolution of the host phase if the leachant is not regularly replaced, (3) pH controls the dissolution rate, and (4) the iodine dissolution rate slows with extended test time (up to 224 days). Based on this work, these parameters should thus be represented when developing a CCRM.

corrosion

Derivation of low-energy Hamiltonians for heavy-fermion materials

Here, by utilizing a multiorbital periodic Anderson model with parameters obtained from ab initio band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both nonmagnetic 4⁢𝑓 0 and magnetic 4⁢𝑓 2 virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn 3 . The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of 4⁢𝑓 materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Decoupling the capacity fade contributions in polymer electrolyte-based high-voltage solid-state batteries

Polymer electrolyte (PE)-based solid-state batteries (PE-SSBs) made with high-voltage cathodes are known to suffer from severe capacity fade, stemming primarily from the poor oxidative stability of most PEs under high-voltage cycling conditions. PEs also suffer from greater ion-transport limitations compared to liquid or solid electrolytes. However, often, these limitations are collectively stated to be responsible for the observed capacity fade, and it is challenging to decouple the contributions of different factors. Herein, a tunable cell fabrication platform was developed to systematically investigate and decouple the two primary capacity fade drivers (cell impedance growth and kinetic limitations), while keeping the other cell parameters constant. Three PE types with distinct transport characteristics were compared. By utilizing a voltage profile analysis method, the contribution of the cell's internal impedance growth was quantitatively decoupled from the kinetic limitations stemming from the high concentration gradient in the polymer catholyte and slow charge transfer reactions. We demonstrate that the high interfacial impedance did not necessarily correlate with the high capacity fade rate. Kinetic limitations that are not reflected by impedance measurements can play a dominant role in causing cumulative capacity decay.

Ock, Ji-young [Oak Ridge National Laboratory (ORNL

Growth of deuterium supersaturated surface layer with increasing ion flux and fluence in plasma-exposed tungsten

Deuterium supersaturated surface layer (DSSL) in tungsten, a few nm thick layer exhibiting extremely high deuterium content (> 5%), has been studied as a function of deuterium plasma ion flux and fluence. Tungsten samples were exposed at 400 K and deuterium ion energy of ∼ 60 eV. The deuterium ion flux spanned over an order of magnitude, being 7.3 × 10 20 , 4.2 × 10 21 , and 3.8 × 10 22 D/m 2 . The fluence ranged from 3.4 × 1024 to 3.4 × 1025 D/m2. The samples and their deuterium content were analyzed by nuclear reaction analysis (NRA), scanning transmission electron microscopy (STEM), and thermal desorption spectroscopy (TDS). The largest thickness of DSSL was found to be around 10.5 nm and was observed in the case of the highest exposure flux and fluence, whereas the layer was only about 3.8 nm thick in the case of the lowest value of the two parameters. The thickness of the DSSL was found to monotonically increase with both deuterium flux and fluence. Finally, similar to the thickness, the estimated D concentration seems to follow the same trend, increasing from the lowest value of around 3 at.% to the highest value of around 5 at.%.

Deuterium

Dark energy survey year 3 results: likelihood-free, simulation-based w CDM inference with neural compression of weak-lensing map statistics

We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock dark energy survey lensing data sets. For wCDM inference, for which we allow –1 < w < –$\frac{1}{3}$⁠, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives Ω m = 0.283$^{+0.020}_{–0.027}$⁠, S 8 = 0.804$^{+0.025}_{–0.017⁠}$, and w < –0.80 (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter (Ω⁠ DE , w⁠) precision.

79 ASTRONOMY AND ASTROPHYSICS

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer

Photosynthetic capacity is reduced by warming but unaffected by elevated CO2 in seedlings of five boreal tree species

Abstract Increasing atmospheric CO2 concentrations fuel global warming, with boreal regions warming at a faster rate than many other areas. Boreal forests are an important component of the global carbon cycle, yet we have little data on photosynthetic responses of boreal trees to elevated CO2 (EC) and warming. We grew seedlings of 5 widespread North American boreal tree species (from Betula, Larix, Picea, and Pinus) under current (410 ppm) or elevated (750 ppm) CO2 and either ambient (+0 °C) or increased (+4 °C or +8 °C) temperature, then measured photosynthetic traits over a range of leaf temperatures. Our results were generally consistent across species: photosynthetic capacity (maximum rates of Rubisco carboxylation, Vcmax, and electron transport, Jmax) was unaffected by EC but decreased under +8 °C warming. Accordingly, net photosynthesis measured at the growth CO2 concentration (Agrowth) was reduced under warming and increased under EC. The thermal optimum for Agrowth (ToptA) increased by ∼1.8 °C with EC but increased with warming in only two species. In contrast, the activation energies and thermal optima for Vcmax and Jmax, which are used to estimate photosynthesis in Earth System Models, were unaffected by growth environment. There were a few interactions between growth, CO2, and warming. These results suggest increased photosynthesis of widespread boreal tree species under EC may be offset by future reductions in photosynthetic capacity related to warming. We also show that the temperature sensitivities of parameters used to estimate global photosynthesis in large-scale models are generally unaffected by simulated climate change in these species.

Plant Sciences

Giant Exfoliation Induced Magnetic Coercivity in Fe 3 GaTe 2

Permanent magnets with strong anisotropy and high coercivity underpin modern information and energy technologies, yet rare-earth-free alternatives remain limited. Here, we show that thickness engineering via mechanical exfoliation induces hard magnetic behavior in the van der Waals ferromagnet Fe 3 GaTe 2 . Bulk crystals exhibit Curie temperatures above 350 K but negligible room-temperature coercivity. When thinned below ∼100 nm, the coercive field is dramatically enhanced, reaching nearly 1 T at room temperature for in-plane fields—comparable to conventional hard magnets. Micromagnetic analysis reveals a crossover in magnetization reversal from domain-mediated processes in bulk samples to quasi-coherent rotation in thin flakes, driven by increased effective anisotropy and suppressed domain formation. This thickness-dependent transition enables tuning of magnetic hardness without chemical modification. Combined with high saturation magnetization and robust room-temperature performance, Fe 3 GaTe 2 emerges as a promising rare-earth-free material for spintronic applications. Its layered structure further allows integration into van der Waals heterostructures, where large in-plane coercivity can stabilize magnetic states against perturbations and interlayer coupling, offering potential for high-density nonvolatile memory and domain-wall-based devices.

36 MATERIALS SCIENCE

Examination of Regional Trends in Low Level Cloud Properties Found in the Aqua-MODIS Satellite Record

Clouds have a pronounced influence on the Earth?s climate. Relative to cloud free conditions, they cool the planet by increasing the amount of solar radiation reflected back to space and reducing the amount of sunlight reaching the surface, but they warm the planet by decreasing the amount of thermal infrared radiation escaping to space and increasing the amount reaching the surface (a greenhouse effect). The global mean net cloud radiative effect (CRE) is about -20 W/m2, a cooling effect at both the top-ofatmosphere and surface. Given the magnitude of CRE?s, it is expected that changes in cloud properties could be a significant factor in climate change due to anthropogenic forcing?s, yet cloud feedbacks are not well known and remain one of the largest uncertainties in climate prediction. This paper explores relationships between coincident observations of atmospheric aerosols, clouds and radiation derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) and from the Clouds and the Earth?s Radiant Energy System (CERES) instruments on the Aqua satellite. We investigate several interesting regional trends that have emerged in the nearly 18-year satellite record that suggest correlation between changes in low-level cloud properties and changes in aerosol optical depth that may be associated with changes in pollution emissions and possibly with other factors. MERRA reanalysis of meteorological conditions and aerosol particulate species are investigated to help better understand the potential mechanisms responsible for the observed cloud property trends. Finally, we analyze a new CERES flux by cloud type dataset in order to try and isolate the associated trends in low-level cloud radiative effects. It is anticipated that this study using long-term observations of clouds, aerosols and radiative fluxes combined with model reanalysis data will contribute to an improved understanding of cloud climate feedbacks.

William L Smith Jr.

Phase-field modeling of stored-energy-driven grain growth with intra-granular variation in dislocation density

Abstract We present a phase-field (PF) model to simulate the microstructure evolution occurring in polycrystalline materials with a variation in the intra-granular dislocation density. The model accounts for two mechanisms that lead to the grain boundary migration: the driving force due to capillarity and that due to the stored energy arising from a spatially varying dislocation density. In addition to the order parameters that distinguish regions occupied by different grains, we introduce dislocation density fields that describe spatial variation of the dislocation density. We assume that the dislocation density decays as a function of the distance the grain boundary has migrated. To demonstrate and parameterize the model, we simulate microstructure evolution in two dimensions, for which the initial microstructure is based on real-time experimental data. Additionally, we applied the model to study the effect of a cyclic heat treatment (CHT) on the microstructure evolution. Specifically, we simulated stored-energy-driven grain growth during three thermal cycles, as well as grain growth without stored energy that serves as a baseline for comparison. We showed that the microstructure evolution proceeded much faster when the stored energy was considered. A non-self-similar evolution was observed in this case, while a nearly self-similar evolution was found when the microstructure evolution is driven solely by capillarity. These results suggest a possible mechanism for the initiation of abnormal grain growth during CHT. Finally, we demonstrate an integrated experimental-computational workflow that utilizes the experimental measurements to inform the PF model and its parameterization, which provides a foundation for the development of future simulation tools capable of quantitative prediction of microstructure evolution during non-isothermal heat treatment.

Materials Science

Emissivity measurements of CuCrZr alloy

The Facility for Rare Isotope Beams (FRIB) heavy-ion accelerator, in user operation since 2022, produces rare isotope beams via interactions of high-intensity stable ion beams with a graphite production target. Approximately 20–40% of the primary beam power is deposited in the target, while the remaining 60–80% is absorbed by the beam dump. The minichannel beam dump (MCBD), currently operated at 20 kW and designed for operation up to 50 kW, uses CuCrZr alloy absorber plates. Thermal validation and thermal cycling tests of the MCBD were conducted at the Applied Research Laboratory (ARL) at Pennsylvania State University. Temperature measurements were obtained from an infrared (IR) camera. Since accurate temperature determination requires reliable emissivity values, the emissivity of CuCrZr was measured using the IR camera validated against thermocouple reference temperatures up to approximately 650 °C. The measurements were conducted under a vacuum level of approximately 10 -5 torr to minimize emissivity variations due to surface oxidation. The emissivity of CuCrZr was determined to be 0.057 ± 0.009 using a constant fit to the measured data over the surface temperature range from 100–650 °C.

Accelerator Subsystems and Technologies

Spatiotemporal Studies of Soluble Inorganic Nanostructures with X‐rays and Neutrons

This Review addresses the use of X-ray and neutron scattering as well as X-ray absorption to describe how inorganic nanostructured materials assemble, evolve, and function in solution. We first provide an overview of techniques and instrumentation (both large user facilities and benchtop). We review recent studies of soluble inorganic nanostructure assembly, covering the disciplines of materials synthesis, processes in nature, nuclear materials, and the widely applicable fundamental processes of hydrophobic interactions and ion pairing. Reviewed studies cover size regimes and length scales ranging from sub-Ångström (coordination chemistry and ion pairing) to several nanometers (molecular clusters, i.e. polyoxometalates, polyoxocations, and metal-organic polyhedra), to the mesoscale (supramolecular assembly processes). Reviewed studies predominantly exploit 1) SAXS/WAXS/SANS (small- and wide-angle X-ray or neutron scattering), 2) PDF (pair-distribution function analysis of X-ray total scattering), and 3) XANES and EXAFS (X-ray absorption near-edge structure and extended X-ray absorption fine structure, respectively). While the scattering techniques provide structural information, X-ray absorption yields the oxidation state in addition to the local coordination. Our goal for this Review is to provide information and inspiration for the inorganic/materials science communities that may benefit from elucidating the role of solution speciation in natural and synthetic processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Numerical Investigation of Heat Transfer and Fluid Flow within Electrochemical Hydrogen Peroxide Generation Unit

Long-term manned space missions require the onboard production of disinfectants essential for maintaining crew health and supporting life systems. Currently, disinfection aboard the International Space Station (ISS) relies on disposable wetted wipes, which are regularly resupplied from Earth. This approach imposes a significant burden on resupply logistics, storage, and waste management. To address these challenges and support future missions, efforts are underway to develop an in-situ solution that electrochemically generates hydrogen peroxide disinfectant using onboard resources. In collaboration with NASA, Faraday Technology, Inc. has advanced this concept through a series of Small Business Innovation Research (SBIR) projects, resulting in the development of a Peroxide Generation Unit (PGU). The PGU can produce up to 3 wt.% hydrogen peroxide on-demand at a rate of 1 liter per day, providing a sustainable alternative to Earth-dependent supplies. The resulting aqueous hydrogen peroxide (H₂O₂) is an effective disinfectant, safe for crew use, compatible with spacecraft systems, and free from volatiles, off-gassing, or residues. This innovation offers a reliable, efficient solution for onboard disinfection, reducing dependence on Earth-based resupply while ensuring the health and safety of space crews. Generating hydrogen peroxide at the required rate needs high voltages and currents, exceeding 20V and 2A respectively, which leads to significant heat generation from Joule heating. This temperature rise poses a risk to sensitive system components, especially critical and expensive membranes that can degrade under thermal stress. To mitigate this risk, the thermal, fluid, and electrical flows within the system are modeled computationally using the commercial software COMSOL. The numerical simulations are validated against experimental data from both sub-scale and alpha-scale systems. Once verified, the model is employed to identify thermal hotspots, investigate their underlying causes, and explore solutions to prevent them.

Life Support Systems (LSS)

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust

Heat Transfer Fluids as Co‐Diluents in Localized High‐Concentration Electrolytes for High‐Rate Lithium Metal Batteries With Enhanced Safety

Localized high-concentration electrolytes (LHCEs) have been identified as promising electrolyte formulations for lithium metal batteries, due to their effective interphase formation and promotion of compact Li deposition, yet their practical implementation is often limited by reduced ion transport kinetics. In this study, two industrially established fluorinated ethers are identified for the first time in battery research as effective co-diluents as they combine a broad electrochemical stability window with a low viscosity and intrinsic non-flammability. Incorporating these components, commonly used as heat transfer fluids, yields safer, less flammable electrolyte formulations with enhanced ion mobilities. In particular, the ternary co-diluent formulation shows improved ion mobility by reducing the electrolyte's viscosity while limiting excessive ion clustering. Based on the improved electrolyte transport kinetics, lower overvoltages and higher Coulombic efficiencies at current densities ≥ 1 mA cm −2 are achieved with the ternary co-diluent blend, resulting in markedly extended cycle life in an application-oriented zero-excess pouch cell compared with the baseline system. Complementary electrochemical and ex situ analysis of harvested electrodes at moderate current densities reveals no discernible differences in interphase morphology and composition, suggesting enhanced ion mobility as the primary cause of the improved high-rate performance.

electrolyte diluent

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