Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models, Chemical”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Enabling accurate chemical modeling of shocked energetic materials using a machine learning interatomic potential

Understanding the complex chemistry of organic materials under dynamic compression is important for many applications, but it is challenging due to the large number of reactions occurring at various time scales. Here, in this study, we develop a machine learning potential based on Chebyshev polynomials to study the insensitive energetic material 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) under detonation. We discuss a strategy for constructing diverse training data needed to capture the complex chemistry of TATB. Our potential demonstrates strong transferability across a wide range of thermodynamic conditions and other explosives, enabling accurate and reliable chemical modeling of organic materials under extreme conditions. The efficiency of our approach allows for simulations over several nanoseconds and for large system sizes, providing detailed insights into the chemistry of shocked TATB. The model accurately reproduces experimental Hugoniot equation of state data, and our simulations reveal the rapid formation of nitrogen-rich carbon clusters following shock. The methods and datasets developed here offer a robust framework for accurate chemical modeling of other shocked organic energetic materials.

Chemistry↗
DOE OSTI2026

Predictive Chemical Kinetic Modeling: Where We Succeed, Where We Struggle, and What Comes Next

Chemical kinetic modeling plays a foundational role in fields ranging from energy to environmental science, pharmaceuticals, and advanced materials. The past two decades have seen remarkable progress, particularly in modeling gas-phase reactions for thermochemical processes, leading to impactful industrial applications such as steam cracking and air quality management. However, new challenges are emerging. The successful development of systematic methodologies for the description of gas-phase kinetics opens the possibility to apply the same approach to the study of more challenging systems. Here, we review recent advances, including ab initio transition state theory-based master equation estimation of elementary rates, automated mechanism generation, machine-learning-assisted kinetics, and uncertainty quantification, and discuss the advances needed to apply the same methodological approach in areas such as heterogeneous catalysis, electrochemistry, liquid-phase and solid-state reactivity, and multiscale model integration. We advocate for the development of targeted tools, especially methods that go beyond empirical tuning toward first-principles-based predictions. We highlight the need for accessible software and AIaugmented workflows to democratize modeling for industry and academia alike. In this perspective, we call attention to not only what has worked but also what remains unsolved, advocating to avoid overemphasizing successes in scientific works at the expense of realism. The next decade should focus on predictive capability, physical accuracy, and community infrastructure (e.g., databases and services) to enable innovation across diverse fields. We argue that kinetic modeling, properly equipped, can accelerate discovery far beyond its traditional domains.

ab initio calculations↗
DOE OSTI2024

Improvement of mechanistic fuel-cladding chemical interaction modeling in BISON

This report describes work performed during FY2024 under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to inform and improve mechanistic models of fuel-cladding chemical interaction (FCCI) in metallic fuel. For fuel-side FCCI, atomistic simulations were performed to determine the diffusivity of iron (Fe) in the $\alpha$ and $\gamma$ phases of uranium (U). A model of liquid penetration of cladding due to melting of the fuel-side FCCI region was updated to account for the finite size of the FCCI region, and the model was validated through comparison with tests performed in the Fuel Behavior Test Apparatus (FBTA). For cladding wastage formation, a reduced-order model was improved by comparison with a multi-scale mechanistic model to better quantify the ROM parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

A Vertically Resolved Canopy Improves Chemical Transport Model Predictions of Ozone Deposition to North Temperate Forests

Abstract Dry deposition is the second largest tropospheric ozone (O 3 ) sink and occurs through stomatal and nonstomatal pathways. Current O 3 uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O 3 fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O 3 fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O 3 deposition and its drivers. We show that GEOS‐Chem overpredicts observed O 3 fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O 3 deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O 3 flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O 3 deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O 3 fluxes; and (c) the CTM big‐leaf treatment overestimates O 3 ‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O 3 predictions.

Vermeuel, Michael P. [Department of Soil, Water, a↗
DOE OSTI2024

Quantifying the impacts of marine aerosols over the southeast Atlantic Ocean using a chemical transport model: implications for aerosol–cloud interactions

The southeast Atlantic region, characterized by persistent stratocumulus clouds, has one of the highest uncertainties in aerosol radiative forcing and significant variability across climate models. In this study, we analyze the seasonally varying role of marine aerosol sources and identify key uncertainties in aerosol composition at cloud-relevant altitudes over the southeast Atlantic using the GEOS-Chem chemical transport model. We evaluate simulated aerosol optical depth (AOD) and speciated aerosol concentrations against those collected from ground observations and aircraft campaigns such as LASIC, ORACLES, and CLARIFY, conducted during 2017. The model consistently underestimates AOD relative to AERONET, particularly at remote locations like Ascension Island. However, when compared with aerosol mass concentrations from aircraft campaigns during the biomass burning period, it performs adequately at cloud-relevant altitudes, with a normalized mean bias (NMB) between -3.5 % (CLARIFY) and -7.5 % (ORACLES). At these altitudes, in the model, organic aerosols (63 %) dominate during the biomass burning period, while sulfate (41 %) prevails during austral summer, when dimethylsulfide (DMS) emissions peak in the model. Our findings indicate that marine sulfate can account for up to 69 % of total sulfate during the high-DMS period. Sensitivity analyses indicate that refining DMS emissions and oxidation chemistry may increase sulfate aerosol produced from marine sources, highlighting that there remains large uncertainty as to the role of DMS emissions in the marine boundary layer. Additionally, we find marine primary organic aerosol emissions may substantially increase total organic aerosol concentrations, particularly during austral summer. This study underscores the imperative need to refine marine emissions and their chemical transformations, as aerosols from marine sources are a major component of total aerosols at cloud-relevant altitudes and may impact uncertainties in aerosol radiative forcing over the southeast Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling Approach for the Aluminum-clad Dry Storage Pilot using HFIR Fuel

To confirm that the dry storage of aluminum-clad research reactor spent nuclear fuel (ASNF) will remain within the safety envelope after applied drying schemes and that the resulting evolution of the gas space composition, temperature, and pressure conditions are understood, a dry storage pilot project is being established. The pilot will incorporate an instrumented lid for discrete interval or for on-demand gas composition and temperature monitoring of two DOE Standard Canisters (DSCs) loaded with three High Flux Isotope Reactor (HFIR) inner cores per DSC. Each DSC would be subjected to a separate alternative candidate drying scheme. Canisters will undergo 1 to 5 years of monitoring, including internal temperature and gas sampling to track pressure and composition changes. This report outlines the approach for modeling the ASNF-in-canister behavior in terms of evolving gas space conditions for the ASNF dry storage pilot using HFIR fuel. The ASNF has an adherent surface oxyhydroxide layer comprised of boehmite/bayerite that generates hydrogen when subjected to irradiation. Three-dimensional multi-physics computational fluid dynamics simulations will be executed to compute the thermal field within the DSC and provide inputs to a chemical model employed to compute pressure buildup as hydrogen is generated in the system. Implemented in Cantera, the chemical model solves gas phase and aluminum oxyhydroxide surface-mediated radiolysis reactions. Gas phase reactions are sourced from Wittman and Hanson (2015), whereas surface-mediated reactions are incorporated by fitting experimental data using an optimization algorithm (Abboud, 2023). Water radiolysis reactions from Wren and Ball (2001) are adopted with modifications as described in Abboud (2023c). Understanding the effect of the hydrogen buildup over time is important for long-term storage safety considerations. Modeling results will include the canister pressure, temperature, and composition evolution from the initial helium backfill with the addition of radiolytically-evolved chemical species (e.g., hydrogen and oxygen). The specific HFIR cores for the pilot program have not yet been selected, and the overall design is still in development. The CFD-chemical model used for this work will be based on prior models with necessary updates to allow for improved accuracy and efficiency. The experimental data obtained from the HFIR demonstration will be used to improve and validate the computational models to predict the ASNF-in-canister behavior.

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

Flame Chemistry Workshop: a perspective on challenges and strategic actions in combustion experiments and chemical kinetics modeling

Continued progress in the development of predictive models for combustion chemistry—including ignition and flame behavior, species evolution, and combustion system performance—relies on overcoming enduring and emerging challenges in experimental measurements, theoretical formulations, and chemical kinetics mechanism construction. As combustion science continues to coincide with advances in sustainable fuels development, plasma technologies, and automated modeling capabilities, the need for coordinated, community-driven strategies is essential. The Flame Chemistry Workshop (FCWS), held biennially before the International Symposium on Combustion, serves as a dedicated platform to identify, consolidate, and address these challenges in a structured and collaborative manner. This perspective arises from discussions at the 7th FCWS in Milan, Italy (2024), and presents a collective view of the critical barriers currently limiting progress. Across the five technical domains discussed during the 7th FCWS – sustainable fuels combustion, advanced diagnostics for combustion measurements, experiments and modeling in plasma combustion, artificial intelligence and automated methods for theory and mechanisms generation, and chemical kinetic models—a series of persistent and emerging scientific challenges were identified, highlighting the need for deeper integration between three areas: theory, experiments, and modeling. In conclusion, the present article concisely describes present challenges that were identified in each of the technical domains in an effort to streamline and coordinate solutions to accelerate progress in combustion science.

Chemical kinetics↗
DOE OSTI2025

Hyperplane decision trees as piecewise linear surrogate models for chemical process design

Recent trends in chemical engineering research point towards an increasing reliance on data-driven modeling approaches. Neural networks, for instance, have proven to be accurate when data is plentiful and high-dimensional, but in many cases, they require computationally-intensive training procedures. Here, in this work, we describe hyperplane decision trees (HT) as a highly expressive and low-compute machine learning model architecture. These models are locally linear and have linear decision boundaries, resulting in a piecewise linear model of the data. This property allows them to be converted into mixed-integer linear constraints which can be globally optimized. Our open-source PyTorch implementation of this method is a fast, flexible, and accessible way to build accurate piecewise linear models of data.

Decision trees↗
DOE OSTI2024

Toward first principles-based simulations of dense hydrogen

Accurate knowledge of the properties of hydrogen at high compression is crucial for astrophysics (e.g., planetary and stellar interiors, brown dwarfs, atmosphere of compact stars) and laboratory experiments, including inertial confinement fusion. There exists experimental data for the equation of state, conductivity, and Thomson scattering spectra. However, the analysis of the measurements at extreme pressures and temperatures typically involves additional model assumptions, which makes it difficult to assess the accuracy of the experimental data rigorously. On the other hand, theory and modeling have produced extensive collections of data. They originate from a very large variety of models and simulations including path integral Monte Carlo (PIMC) simulations, density functional theory (DFT), chemical models, machine-learned models, and combinations thereof. At the same time, each of these methods has fundamental limitations (fermion sign problem in PIMC, approximate exchange–correlation functionals of DFT, inconsistent interaction energy contributions in chemical models, etc.), so for some parameter ranges accurate predictions are difficult. Recently, a number of breakthroughs in first principles PIMC as well as in DFT simulations were achieved which are discussed in this review. Here we use these results to benchmark different simulation methods. We present an update of the hydrogen phase diagram at high pressures, the expected phase transitions, and thermodynamic properties including the equation of state and momentum distribution. Furthermore, we discuss available dynamic results for warm dense hydrogen, including the conductivity, dynamic structure factor, plasmon dispersion, imaginary-time structure, and density response functions. We conclude by outlining strategies to combine different simulations to achieve accurate theoretical predictions that are based on first principles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Phase Field Modeling of Chemical Reaction Related Damage Evolution in Environmental Barrier Coatings

The advent of next-generation engines necessitates materials capable of withstanding temperatures beyond the reach of current superalloys. SiC-based ceramic matrix composites, augmented with environmental barrier coatings (EBCs), present a promising materials solution. Given the active search for effective and durable EBCs, there is a pressing need for modeling tools to understand and predict damage evolution in these materials to help accelerate their development. This study introduces a phase-field model (PFM) designed to simulate the thermally grown oxides (TGO) and phase transformation in the degradation and failure of EBCs. The model accounts for the severe volume expansion due to oxidation, alongside phase transformations and porosity evolution during thermal cycling, offering a comprehensive view of the damage processes. Simulation results are validated against experimental findings reported in the literature, establishing the model's potential as a significant tool for understanding and improving the resilience of EBCs in cyclic oxidative environments.

fast-diffusion path↗
DOE OSTI2025

Formation of late-generation atmospheric compounds inhibited by rapid deposition

Reactive organic carbon species are important fuel for atmospheric chemical reactions, including the formation of secondary organic aerosol. However, in parallel to atmospheric oxidation processes, deposition can remove compounds from the atmosphere and impact downstream environments. To understand the impact of deposition on atmospheric oxidation, we present a framework for predicting and visualizing the fate of a molecule on the basis of the physicochemical properties of compounds (Henry’s law constant, vapour pressure and reaction rate constants), which are used to estimate timescales for oxidation and deposition. Further, by implementing our deposition rates in chemical models, we show that deposition substantially suppresses atmospheric reactivity and aerosol formation by removing early-generation products and preventing the formation of large fractions (up to 90%) of downstream, late-generation compounds. Deposition is frequently missing in the laboratory experiments and detailed chemical modelling, which probably biases our understanding of atmospheric composition.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluating the Use of Foundational Chemical Language Models in Multimodal Graph Fusion

Rapid and accurate prediction of the physicochemical properties of molecules given their structures remains a key challenge in cheminformatics. Machine learning approaches offer high-throughput options, but the optimality of inductive biases and data representations are up for debate. For example, BERT-based masked language models (MLMs) can be trained in a self-supervised way on hundreds of millions to billions of readily available SMILES strings. Another option is graph neural networks (GNNs), which can operate directly on molecular structures. Yet, generating accurate molecular geometry is computationally expensive, leading to a relative scarcity in data compared to SMILES strings. It is attractive to combine these two paradigms by pre-training an LM on a large corpus of SMILES strings and embedding these representation into a geometric graph neural network. Despite the promise of such an approach, and contrary to previous studies, we find mixed results with the combination of the LMs and GNNs on several molecule datasets. In particular, we found evidence for improvement on the FreeSolv and QM7 benchmarks, but degraded performance on the ESOL, LIPO and QM9 datasets compared to a GNN baseline.

Francel, Collin [University of Alabama]↗
DOE OSTI2025

bmdrc: Python package for quantifying phenotypes from chemical exposures with benchmark dose modeling

Though chemical exposures are known to potentially have negative impacts on health, including contributing to chronic diseases such as cancer, the quantitative contribution of risk is not fully understood for every chemical. A commonly used approach to quantify levels of risk is to measure the proportion of organisms (such as a total number of zebrafish on a plate or mice in a cage) with abnormal behavioral responses or morphology at increasing concentrations of chemical exposure. A particular challenge with processing the proportional data from these assays is the appropriate estimation of chemical concentration levels that result in malformations or acute toxicity, as these values typically vary between experimental measurements. The recommended approach by the Environmental Protection Agency (EPA) is to fit benchmark dose curves with specific filters and model fitting steps, which are crucial to properly processing the proportional data. Several tools exist for the fitting of benchmark dose response curves, but none are standalone Python libraries built to process both morphological and behavioral data as proportions with all the EPA recommended filters, filter parameters, models, and model parameters. Thus, here we present the benchmark dose response curve (bmdrc) Python library, which was built to closely follow these EPA guidelines with helpful visualizations of filters and fitted model curves, and reports for reproducibility purposes. bmdrc is open-source and has demonstrated utility as a support package to an existing web portal for information on chemicals (https://srp.pnnl.gov). Our package will support any toxicology analysis where the response is a proportional value at increasing levels of a concentration of a chemical or chemical mixture.

Superfund↗
DOE OSTI2026

SmileyLlama: modifying large language models for directed chemical space exploration

Here we show that large language models (LLMs) can be transformed via supervised fine-tuning of engineered prompts into SmileyLlama for exploring the chemical space of drug molecules. We benchmark SmileyLlama against pretrained LLMs and chemical language models trained from scratch for generating valid and novel drug-like molecules, and use direct preference optimization to both improve SmileyLlama’s adherence to a prompt and as part of the iMiner reinforcement learning framework to predict molecules with optimized three-dimensional conformations and high binding affinity to drug targets. By training an LLM to speak directly as a chemical language model, while retaining most of its natural language capabilities, we show that SmileyLlama can reliably generate molecules with user-specified properties rather than acting only as a chatbot with knowledge of chemistry or as a virtual assistant. While SmileyLlama is geared toward drug discovery, the supervised fine-tuning/direct preference optimization/LLM framework can be extended to other chemical, biological and materials applications.

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

Photochemically Induced Acousto-optics Fluid Simulations

PIAFS is a finite-difference code to solve the compressible Navier-Stokes equations with chemical heating on Cartesian grids. It models chemical reactions of air (oxygen and carbon dioxide) with ozone subject to radiation. It uses a high-order WENO spatial discretization and explicit Runge-Kutta time integration. It is capable of parallel simulations using MPI. The code is written in C/C++.

Oudin, AlbertineN [Lawrence Livermore National Lab↗
DOE OSTI2026

Modeling Radiolysis and Chemical Reactions during Dry Storage of Aluminum-clad Spent Nuclear Fuel

After aluminum-clad spent nuclear fuel (ASNF) is removed from the reactor, it is initially stored in spent fuel pools, which are specially designed water-filled basins that provide temporary cooling to reduce the temperature of the fuel assemblies and provide radiation shielding. ASNF continues to generate heat due to the radioactive decay of elements within the fuel, which persists for many years post-shutdown as the residual radioactive products decay into more stable elements. During the wet storage period, an oxyhydroxide layer composed of boehmite/bayerite forms on the surfaces of the aluminum cladding from exposure to water in the pools. Road-ready packaging for long-term disposition of the ASNF involves dry storage in helium backfilled DOE standard canisters (DSCs). When the ASNF is removed from water storage and dried, most of the water is removed, but some physisorbed and chemisorbed water remains in the oxyhydroxide layers. This residual water can produce hydrogen when exposed to radiation from the ASNF during dry storage. Predicting hydrogen accumulation over time in the DSCs is critical for long-term storage considerations. Previous modeling efforts have developed coupled computational fluid dynamics (CFD)-chemical models to simulate temperature, pressure, and gas phase concentrations within the DSCs. These models use the thermal field predicted by CFD as input to a radiolysis model for the gas phase and the surface oxyhydroxide layer chemistry. Given the long storage period of the DSCs and the impracticality of long-term experiments, a simulation-based approach is necessary to assess chemical evolution within the canisters. This study advances the development of a modeling framework designed to simulate the chemical evolution of spent fuel canisters. Both thermal and radiation-driven reactions are considered, with radiation kinetics quantified using G-values. Sensitivity analysis identifies key parameters influencing species composition. Reaction pathway diagrams offer insight into dominant species formation routes, enabling more effective comparisons between model predictions and experimental observations, particularly regarding the production of hydrogen. Results show that the model predicts significant hydrogen gas production with minimal oxygen generation, primarily due to hydrogen formation via boehmite pathways. These findings underscore the importance of accurately characterizing surface-bound species and radiolysis kinetics. A deeper understanding of these mechanisms is critical for evaluating the long-term safety of nuclear waste storage.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2025

Enhancing Drinking Water Quality Modeling: Leveraging Physics Informed Neural Networks for Learning with Imperfect Reaction Models and Partial Data

Chemical kinetics models, typically formulated as systems of ordinary or partial differential equations, are valuable tools for simulating drinking water quality. However, these models often face inaccuracies due to discrepancies between the laboratory and the real-world conditions, as well as limitations in experimental analytical methods, hindering the accurate representation of the true underlying chemical mechanisms. In this study, we propose a Physics Informed Neural Network (PINN), using the eXtreme Theory of Functional Connections, to improve the prediction of chemical concentrations over time. The PINN method accounts for imperfect chemical models and incorporates partial data to improve predictions. Focusing on reactions describing water disinfection residual and disinfectant byproduct formation, which are crucial for public health and regulatory compliance, we demonstrate that the PINN model is able to accurately predict the concentrations of chemical species across various pH values. Notably, the model extends its accuracy to predict concentrations of chemical species not originally included in its training data. The developed method can be extended to a variety of chemical systems, offering a wide array of potential applications.

13 HYDRO ENERGY↗