Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mechanistic Modeling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Documentation, assessment, and improvements of fission gas modeling capabilities in BISON

Fission gas behavior critically affects nuclear fuel performance, influencing thermal conductivity, microstructure, swelling, and fuel-cladding interactions. The Simple Integrated Fission Gas Release and Swelling (Sifgrs) model, implemented within the BISON fuel performance code, provides advanced capabilities for modeling these phenomena. This manuscript reviews and assesses Sifgrs’ modular structure and submodels describing intragranular, intergranular, and release mechanisms. The modular implementation enables rapid improvements as our understanding advances and new fuel forms emerge. After an overview of capabilities, we present a comprehensive validation suite, comparing predictions to experimental data from tests under normal and transient conditions, including reproduction of the Vitanza curve. Results demonstrate Sifgrs’ accuracy, robustness, and applicability. Current limitations and ongoing efforts are discussed, including mechanistic, multiscale modeling, uncertainty quantification, and validation within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, supporting the safe and efficient use of nuclear fuels.

BISON↗
DOE OSTI2026

Toward the performance assessment of advanced nuclear waste forms: temperature dependence of lanthanide borosilicate glass dissolution

Lanthanide borosilicate (LaBS) glasses are among the most promising waste forms for the immobilization of high-level radioactive waste generated from advanced nuclear fuel cycles. However, the temperature dependence of their dissolution kinetics remains poorly understood and constrained, limiting the integration of these materials into established performance assessment models. Here, we investigate the dissolution behavior of the legacy AmCm2-19 LaBS glass and the benchmark alkali aluminoborosilicate ISG-1 in deionized water between 50 °C and 250 °C using ASTM C1285 (Product Consistency Test-B) protocols. For AmCm2-19 LaBS glass, normalized elemental release rates for boron and silicon increase with temperature before plateauing near 150 °C, consistent with solubility-limited behavior. From data obtained at 50 °C and 100 °C, Arrhenius analysis yields activation energies of E a (B) = 24.8 ± 0.3 kJ mol⁻¹ and E a (Si) = 14.4 ± 0.2 kJ mol⁻¹, similar or slightly lower than those previously reported for two other compositions of LaBS glasses. No secondary phases or alteration layers were detected by SEM-EDX or pXRD. These results establish one of the first temperature-dependent kinetic datasets for LaBS glass dissolution, providing quantitative parameters to inform mechanistic corrosion models and predictive simulations of glass degradation in geological disposal environments.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Advanced Fuels Campaign Execution Plan

The Advanced Fuels Campaign (AFC) Execution Plan outlines the strategy, mission, scope, near-term and long-term goals, structure, and organization associated with nuclear fuels and materials research, development, and demonstration activities within the Department of Energy’s (DOE) Nuclear Fuel Cycle and Supply Chain (NFCSC) program. NFCSC has been given responsibility to identify and mature advanced fuel technologies for the DOE using a science-based approach, focused on developing a fundamental understanding of nuclear fuels and materials to drive development of integrated nuclear fuel and materials technology. This science-based approach combines theory, experiments, and multiscale modeling and simulation to achieve a predictive understanding of relevant behaviors ranging from fuel fabrication processes (and their resulting fuel microstructures) through fuel/cladding performance under irradiation (in contrast to more empirical, observation-based approaches frequently used in fuel performance modeling and fuel qualification). The traditional scope of AFC includes the evaluation and development of multiple fuel forms to support two fuel cycle options: once-through and full recycle. The word “fuel” is used generically to include conventional fuels, transmutation targets, and any associated cladding or duct materials. The once-through fuel cycle addresses advanced light water reactor fuels with enhanced performance, extended burnup, and reduced waste generation. In fiscal year (FY) 2012, AFC’s scope expanded to include research, development, and demonstration (RD&D) for light water reactor (LWR) fuels with enhanced accident tolerance. Fuel fabrication activities include the development of innovative methods to enhance process efficiencies, reduce waste, and improve control over as-fabricated fuel microstructural properties to achieve desired in-reactor performance. Using modern modeling and simulation approaches, the objective is to predict fresh fuel properties given the feedstock characteristics and fabrication process parameters. The performance-related activities include small-scale, in-reactor, and out-of-reactor phenomenological testing (distinct from, but synergistic with, integral prototypic testing) and extensive, quantitative characterization (focusing on characterization of fuel and cladding materials at the scale of microstructure) both before and after testing. Larger-scale, prototypic experiments are conducted in concert with phenomenological testing to drive a Fuel Development and Qualification program, incorporating a fundamental understanding of fuel behavior performance characteristics. Then, using the tools developed under the productive science-based approach, fuels will be optimized to meet specific performance requirements, thereby minimizing the need to repeatedly perform large-scale, integral experiments over a wide parametric range as a means of experimental exploration. Two significant initiatives are underway within AFC. First, a gap analysis completed in early FY 2019 identified critical irradiation testing needs that are lacking within the national light water reactor (LWR) fuels testbed since the shutdown of the Halden Reactor in 2018. The identified gaps are for instrumented, prototypic testing of LWR fuels, especially under boiling water reactor conditions, ramp conditions, and conditions leading to fuel failure; these needs exist for supporting current LWR fuels and their possible extension to higher burnups, but are especially urgent relative to near-term development and qualification of accident-tolerant fuels. Recommendations that resulted from the Halden Gap Analysis focused on enhancements at Advanced Test Reactor (ATR) and Transient Reactor Test Facility (TREAT) to fill gaps in testing capabilities relative to these needs. Second, a concerted effort to develop and demonstrate a systematic approach to accelerating the development, testing, and qualification of new fuel systems has been initiated. This is highlighted by a test strategy that combines the considerable advances in multiscale, mechanistic fuel modeling of recent years with a MiniFuel separate effects test program in the High Flux Isotope Reactor (HFIR) and a Fission Accelerated Steady-state Testing (FAST) semi-integral accelerated test program in ATR. This approach is being tested/demonstrated using the metallic fuel system, but if successful it is expected to be extensible to multiple fuel types and diverse applications. This document includes an overview of the NFCSC program, a definition of science-based development of nuclear fuels, near-term goals for Advanced LWR fuels (ALFs), and longer-term goals for Advanced Reactor Fuels (ARFs) RD&D. This includes the activities that will be conducted to achieve success toward the grand challenge, as well as the goals and milestones to be achieved over the next few decades of research and development. Long-term goals are based on the DOE Office of Nuclear Energ

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Pyrite (001) Interface Chemistry is Controlled by a Sulfoxy Termination

Pyrite (FeS 2 ) is the most common sulfide mineral on Earth, forming through inorganic reactions in the crust and oceanic hydrothermal systems and via microbially driven processes in anaerobic sediments. The pyrite–water interface is the site of a wide range of adsorption and reaction processes in Earth systems including oxidation that dramatically affects the geochemistry of surface waters and influences global carbon and oxygen cycles. Mechanistic geochemical models of pyrite interfacial reactivity, however, are limited by the lack of experimentally derived atomistic structures of the reduced and reacting surfaces. Here, in this work, we reveal the atomic-scale structure of the pyrite (001)-water interface that forms at very low oxygen partial pressures, relevant to suboxic environments in Earth. The interface structure and surface speciation were obtained using the crystal truncation rod method supported by ambient-pressure photoelectron spectroscopy and density functional theory calculations. The surface is dominantly composed of disulfide groups bound to a single oxygen atom, forming a sulfoxy group that has no known molecular or bulk mineral analog. This surface is interpreted as the first step in the oxidative dissolution of pyrite. The sulfoxy group is readily protonated through surface acid–base reactions that alter the structure of interfacial water and the free energy of interfacial reactions. Surface iron sites are not oxidized. Surprisingly, this interface can likely develop in equilibrium with bulk pyrite in some reducing and acidic solutions. This termination is therefore likely representative of pyrite surfaces under a vast range of experimental, industrial and Earth conditions.

oxidation↗
DOE OSTI2026

From microbial diversity to functional potential using dimensionality reduction

The high dimensionality of microbial diversity data from ‘omics observations can be reduced using Machine Learning, with many recent studies showcasing ML utility for exploratory ecological feature finding and process prediction. Here, we compare the Self Organizing Map (SOM) dimensionality reduction method to the well-documented sample-based Principal Coordinate Analysis (PCoA) and taxa-based Weighted Gene Correlation Network Analysis (WGCNA) using near daily 16S rRNA gene amplicon sequencing data from the 2019 to 2020 MOSAiC International Arctic Drift Expedition. We then map k-means clustering outputs from each method to available metagenomes, extracting functionally distinct seasonal microbial ecotypes in the surface Arctic Ocean. Our results indicate the SOM method better represented expected seasonal transitions and identified a greater number of metabolically distinct functional groups than the more traditional PCoA ordination. Ultimately, we identified four community ecotypes with distinct taxonomic and functional cut-offs driven by seasonality, water mass, and substrate turnover, highlighting the importance of succession in functional diversity for the central Arctic Ocean. These results reinforce ML dimensionality reduction as a meaningful translator in the mining of historical amplicon datasets to address modern mechanistic questions and potentially provide ’omics informed ecotype diversity to leverage in mechanistic biogeochemical models.

Arctic Ocean↗
DOE OSTI2025

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

AFC Metal Fuel Workshop Presentations

Compilation of PowerPoint slide decks from presentations given at the 3rd Annual AFC Metal Fuels Workshop. The purpose of this workshop is to bring the national metallic fuel stakeholders together to discuss status and needs of metallic fuel and review the DOE AFC Metal Fuel R&D. Desired Outcomes - Advance Metallic Fuel Technology: 1. Stakeholder Community Building 2. National Program R&D Accountability, Planning Support, and Increased Stakeholder Advocation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

TES: Modelling Microbes to Predict Post-Fire Carbon Cycling in the Boreal Forest across Burn Severities (Exploratory Proposal)

Wildfires in boreal ecosystems represent a globally significant ecological process that is expected to be sensitive to changing climate as well as forest management strategies, but is currently insufficiently understood and represented in models. We sought to determine whether linking belowground microbial community composition, size, and activity to aboveground properties of burn severity and plant community composition (upland jack pine and upland spruce) would allow us to better model post-fire soil CO 2 fluxes using the Carbon, Organisms, Rhizosphere and Protection in the Soil Environment (CORPSE) model. We used an integrated modelling experimental approach, with mechanistic laboratory experiments designed to inform models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Assessing shellfish water exposure to fecal bacteria pollution in Salish Sea: three-dimensional modeling and implications for monitoring

Fecal bacteria (FB) contamination poses significant risks to shellfish safety and management in coastal and estuarine waters. Despite extensive pollution identification and correction efforts, FB contamination in shellfish-growing areas persists in the Salish Sea, highlighting the need to identify overlooked sources and better understand FB transport from riverine and shoreline inputs to shellfish beds. To address this, a high-resolution three-dimensional hydrodynamic model coupled with FB kinetics was developed and applied to a case study site in Salish Sea—Portage Bay—to simulate freshwater plume circulation, flushing dynamics, and bacterial transport. Daily FB loading from the major freshwater inflow—Nooksack River was generated by both linear interpolation and integrating a machine learning approach (XGBoost), trained on historical hydrological and meteorological data. The model successfully reproduced both the magnitude and seasonal variation of FB concentrations in Portage Bay for the year of 2021, demonstrating that simplified FB kinetics with first-order decay due to mortality was effective in this dynamic coastal environment with short flushing time. Model results identified the Nooksack River as the dominant far-field FB source, while scenario simulations showed that near-field coastal stormwater outfalls elevated local FB levels following rainfall, particularly under low-flow conditions. The XGBoost prediction provided comparable or superior accuracy to linear interpolation, particularly during periods of missing observational data, by capturing short-term variability and event-driven loading more effectively. Integrating data-driven riverine FB inputs with mechanistic coastal numerical modeling provides a robust framework for operational forecasting of shellfish bed exposure risk and supports adaptive monitoring and management of shellfish growing areas in the Salish Sea and similar coastal systems.

Salish Sea↗
DOE OSTI2026

From Cell to System: Accelerated hpc Simulations of BESS Aging under Frequency Regulation and Arbitrage use cases

Lithium-ion battery energy storage systems (BESS) packs have emerged as a leading solution for grid-scale energy storage, enhancing resiliency and balancing load fluctuations. Yet, experimental characterization of large-format LIB packs-particularly to assess performance and degradation over hundreds of cycles - demands substantial hardware investment and multi-year testing campaigns. In this work, we couple a hierarchical, physics-based modeling framework agnostic to electrode chemistries with high-performance computing to accelerate systems level evaluation by upto two orders of magnitude. Building on the open-source liionpack platform, we implement cell, module, and pack-scale electrochemical models enriched with mechanistic aging mechanisms and deploy them on an HPC cluster to simulate 150−200kWh systems over 500 - 1,000 cycles with in days. We subject these virtual B ESS to both constant-current cycling and realistic grid service profiles spanning frequency regulation, ramp-rate support, and energy arbitrage-and quantify the resulting degradation patterns. Our results reveal that localized cell aging can induce substantial nonuniformity at module and pack levels, with service-specific cycling protocols driving distinct aging modes. This rapid, multiscale modeling approach provides a powerful design-space exploration tool for optimizing electrical architecture, control strategies, and operational schedules to prolong pack lifetime and lower total cost of ownership.

Ayalasomayajula, Surya [ORNL] (ORCID:0009000860788↗
DOE OSTI2024

EF-Hand Battle Royale: Hetero-ion Complexation in Lanmodulin

The lanmodulin (LanM) protein has emerged as an effective means for rare earth element (REE) extraction and separation from complex feedstocks without the use of organic solvents. Whereas the binding of LanM to individual REEs has been well characterized, little is known about the thermodynamics of mixed metal binding complexes (i.e., heterogeneous ion complexes), which limits the ability to accurately predict separation performance for a given metal ion mixture. In this paper, we employ the law of mass action to establish a theory of perfect cooperativity for LanM-REE complexation at the two highest-affinity binding sites. The theory is then used to derive an equation that explains the nonintuitive REE binding behavior of LanM, where separation factors for binary pairs of ions vary widely based on the ratio of ions in the aqueous phase, a phenomenon that is distinct from single-ion-binding chemical chelators. We then experimentally validate this theory and perform the first quantitative characterization of LanM complexation with heterogeneous ion pairs using resin-immobilized LanM. Importantly, the resulting homogeneous and heterogeneous constants enable accurate prediction of the equilibrium state of LanM in the presence of mixtures of up to 10 REEs, confirming that the perfect cooperativity model is an accurate mechanistic description of REE complexation by LanM. We further employ the model to simulate separation performance over a range of homogeneous and heterogeneous binding constants, revealing important insights into how mixed binding differentially impacts REE separations based on the relative positioning of the ion pairs within the lanthanide series. In addition to informing REE separation process optimization, these results provide mathematical and experimental insight into competition dynamics in other ubiquitous and medically relevant, cooperative binding proteins, such as calmodulin.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

When can we detect lianas from space? Toward a mechanistic understanding of liana‐infested forest optics

Abstract Lianas, woody vines acting as structural parasites of trees, have profound effects on the composition and structure of tropical forests, impacting tree growth, mortality, and forest succession. Remote sensing could offer a powerful tool for quantifying the scale of liana infestation, provided the availability of robust detection methods. We analyze the consistency and global geographic specificity of spectral signals—reflectance across wavelengths—from liana‐infested tree crowns and forest stands, examining the underlying mechanisms of these signals. We compiled a uniquely comprehensive database, including leaf reflectance spectra from 5424 leaves, fine‐scale airborne reflectance data from 999 liana‐infested canopies, and coarse‐scale satellite reflectance data covering 775 ha of liana‐infested forest stands. To unravel the mechanisms of the liana spectral signal, we applied mechanistic radiative transfer models across scales, establishing a synthesis of the relative importance of different mechanisms, which we corroborate with field data on liana leaf chemistry and canopy structure. We find a consistent liana spectral signal at canopy and stand scales across globally distributed sites. This signature mainly arises at the canopy level due to direct effects of more horizontal leaf angles, resulting in a larger projected leaf area, and indirect effects from increased light scattering in the near and short‐wave infrared regions, linked to lianas' less costly leaf construction compared with trees on average. The existence of a consistent global spectral signal for lianas suggests that large‐scale quantification of liana infestation is feasible. However, because the traits responsible for the liana canopy‐reflectance signal are not exclusive to lianas, accurate large‐scale detection requires rigorously validated remote sensing methods. Our models highlight challenges in automated detection, such as potential misidentification due to leaf phenology, tree life history, topography, and climate, especially where the scale of liana infestation is less than a single remote sensing pixel. The observed cross‐site patterns also prompt ecological questions about lianas' adaptive similarities in optical traits across environments, indicating possible convergent evolution due to shared constraints on leaf biochemical and structural traits.

Environmental Sciences & Ecology↗
DOE OSTI2025

Isolation and characterization of 24 phages infecting the plant growth-promoting rhizobacterium Klebsiella sp. M5al

Bacteriophages largely impact bacterial communities via lysis, gene transfer, and metabolic reprogramming and thus are increasingly thought to alter nutrient and energy cycling across many of Earth’s ecosystems. However, there are few model systems to mechanistically and quantitatively study phage-bacteria interactions, especially in soil systems. Here, we isolated, sequenced, and genomically characterized 24 novel phages infectingKlebsiellasp. M5al, a plant growth-promoting, nonencapsulated rhizosphere-associated bacterium, and compared many of their features against all 565 sequenced, dsDNAKlebsiellaphage genomes. Taxonomic analyses revealed that theseKlebsiellaphages belong to three known phage families (Autographiviridae,Drexlerviridae, andStraboviridae) and two newly proposed phage families (CandidatusMavericviridaeand Ca.Rivulusviridae). At the phage family level, we found that core genes were often phage-centric proteins, such as structural proteins for the phage head and tail and DNA packaging proteins. In contrast, genes involved in transcription, translation, or hypothetical proteins were commonly not shared or flexible genes. Ecologically, we assessed the phages’ ubiquity in recent large-scale metagenomic datasets, which revealed they were not widespread, as well as a possible direct role in reprogramming specific metabolisms during infection by screening their genomes for phage-encoded auxiliary metabolic genes (AMGs). Even though AMGs are common in the environmental literature, only one of our phage families,Straboviridae, contained AMGs, and the types of AMGs were correlated at the genus level. Host range phenotyping revealed the phages had a wide range of infectivity, infecting between 1–14 of our 22 bacterial strain panel that included pathogenicKlebsiellaandRaoultellastrains. This indicates that not all capsule-independent Klebsiella phages have broad host ranges. Together, these isolates, with corresponding genome, AMG, and host range analyses, help build theKlebsiellamodel system for studying phage-host interactions of rhizosphere-associated bacteria.

Science & Technology - Other Topics↗
DOE OSTI2025

The effect of heat treatment on the defect evolution in LPBF 316H stainless steel

This work investigates the effect of processing and heat treatment on defect evolution and irradiation response of laser powder bed fusion (LPBF) 316H stainless steel, with comparisons to LPBF 316L and wrought 316L/316H. Using in-situ and ex-situ ion irradiations across a wide parameter space—temperature (300–675 °C), dose (0.2–25 dpa), dose rate (10 -3 –10 -5 dpa/s), and helium co-implantation (20–2500 appm)—we correlated void swelling, dislocation loop evolution, and segregation behavior with pre-irradiation microstructures (as-built, stress-relieved, solution-annealed, and cold-worked). Results show that swelling is strongly controlled by dislocation density: intermediate densities maximize swelling, while solution annealing or cold working reduce susceptibility to levels comparable to wrought alloys. Low dose rates and helium both promote cavity nucleation and lower the incubation barrier, with helium suppressing the role of dislocation density and driving swelling behavior toward wrought-like response. Loop evolution in SA LPBF 316H resembles wrought 316L but shows localized denuded zones near low-angle grain boundaries. STEM-EDS mapping further revealed Ni segregation at voids and sparse Al-rich oxides, without evidence of Ni–Si precipitates. Collectively, these findings identify dislocation density, helium content, and dose rate as the factors governing swelling in LPBF 316H and provide mechanistic datasets for model validation, supporting LAIN and the qualification of AM austenitic steels for nuclear service.

316 stainless steel↗
DOE OSTI2025

Long-term measurements of ice nucleating particles at Atmospheric Radiation Measurement (ARM) sites worldwide

Ice nucleating particles (INPs) play a critical role in cloud microphysics and precipitation formation, yet long-term, spatially extensive observational datasets remain limited. Here, we present one of the most comprehensive publicly available datasets of immersion-mode INP concentrations using a single analytical method, generated through the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility. INP filter samples have been collected across a broad range of environments – including agricultural plains, Arctic coastlines, high-elevation mountain sites, marine regions, and urban areas – via fixed observatories, mobile facility deployments, and vertically-resolved tethered balloon system operations. We describe the standardized processing and quality assurance pipeline, from filter collection and processing using the Ice Nucleation Spectrometer to final data products archived on the ARM Data Discovery portal. The dataset includes both total INP concentrations and selectively treated samples, allowing for classification of biological, organic, and inorganic INP types. It features a continuous 5-year record of INP measurements from a central U.S. site, with data collection still ongoing. Seasonal and site-specific differences in INP concentrations are illustrated through intercomparisons at −10 and −20 °C, revealing distinct regional sources and atmospheric drivers. We also outline mechanisms for researchers to access existing data, request additional sample analyses, and propose future field campaigns involving ARM INP measurements. This dataset supports a wide range of scientific applications, from observational and mechanistic studies to model development, and provides critical constraints on aerosol-cloud interactions across diverse atmospheric regimes (Creamean et al., 2024, 2020b; https://doi.org/10.5439/1770816).

Creamean, Jessie M. [Colorado State Univ., Fort Co↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗