Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical 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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molar-Mass-Dependent Partitioning of Polyethylene in Nanopores of Model Catalyst Supports from Small-Angle Neutron Scattering

Heterogeneous catalysis offers opportunities to enhance valorization of plastic waste via chemical recycling through control of the upcycled product distributions. Minimizing low-value light hydrocarbons is desired; however, fundamental insights into how to control selectivity are lacking. Here we use contrast variation with small-angle neutron scattering (SANS), model perdeuterated polyethylenes (dPEs), and a model liquid hydrocracking product (tetradecane) to quantify polymer partitioning within mesoporous silica (SBA-15). Polyethylene concentration within the mesopores is increased relative to the bulk solution, and this partitioning increases as the temperature increases. However, this polyethylene partitioning is maximized when the radius of gyration of the polymer chains is comparable to the SBA-15 pore size (10 nm). An increased partitioning at higher temperatures is attributed to entropically driven adsorption of PE within the mesopores. There is no observed preferential partitioning of hexatriacontane (a model oligomer) within the mesopores at the temperatures examined. Furthermore, these results suggest that pore size could promote the selective partitioning of polymer species into the mesopores by size. For plastic upcycling, pore-size-dependent partitioning should increase the probability for the reaction of long polymers over oligomeric and small-molecule polyolefin depolymerization products.

Adsorption↗
DOE OSTI2025

Dual aging pathways of Cu-SSZ-13 SCR catalysts: Hydrothermal vs. sulfur-induced deactivation

Hydrothermal aging (HTA) and chemical poisoning are two primary factors contributing to the real-world degradation of Cu-SSZ-13 SCR catalysts. Investigating field-returned samples offers valuable insights into performance degradation caused by these mechanisms. However, the simultaneous presence of both deactivation pathways complicates the isolation of their individual effects in post-mortem analyses. In this study, we separately prepared model Cu-SSZ-13 SCR catalysts subjected to hydrothermal-aging and sulfur-induced chemical poisoning. Using various characterization techniques, we elucidated the specific role of each aging process in catalyst deactivation and compared the results to real-world field-aged catalysts. Our findings show that hydrothermal aging at 650 °C for 100 h caused dealumination of the zeolite framework but no significant CuO x cluster formation. In contrast, sulfur aging (via sulfur exposure, calcination at 550 °C, and desulfation up to 750 °C) led to CuO x formation without any observable dealumination. On model catalysts, sulfur poisoning was found to reduce Cu mobility and the amount of active Cu sites, thus degrading catalyst activity. Although some activity was recovered upon desulfation, a portion of the initial catalyst activity remained irreversibly lost due to CuO x formation. We demonstrate that this occurs because sulfated species impede the ability of multi-nuclear Cu species (e.g., Cu dimers) to split back into their isolated form, leading to CuSO 4 -clusters that oxidatively desulfate to CuO x species. This degradation pathway explains the significant reduction in activity of field-aged samples, where substantial CuSO 4 -cluster accumulation leads to reduced active Cu and subsequent conversion to CuO x . Furthermore, the conclusions from model catalysts were extended directly to field-aged commercial samples, elucidating the decline in activity and chemical properties during field deployment.

Catalyst Deactivation↗
DOE OSTI2025

Probing Condensed-Phase Structure and Dynamics in Hierarchical Zeolites and Nanosheets for Catalytic Upgradation of Biomass (Final Report)

Understanding complex reaction pathways in systems governed by multi-scale collective interactions across time and length scales remains a central scientific challenge. This project was guided by the hypothesis that the interplay among oligomers, solvents, and active sites can be tuned by a suitable choice of solvation environment and pore architecture in solid-acid catalysts to direct chemical transformations relevant to biomass conversion. Zeolites and zeolite nanosheets were used as model platforms, allowing for the interaction of macromolecules with the surface of the zeolite nanosheets and with smaller pores that host catalytically active sites. To investigate these coupled phenomena, we employ a multi-scale computational framework that integrates molecular-level descriptions with advanced sampling approaches to capture key physical and chemical interactions. Our work through this project improved fundamental understanding of how reactants and solid-acid catalysts interact in solvent-rich environments, thereby enabling the rational design of catalytic systems that upgrade biomass with enhanced selectivity and energy efficiency. In addition, the project developed advanced sampling methodologies critical for disentangling complex, reactive processes in multi-component catalytic environments.

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

Effect of Solvents on Lignin–Surface Interactions via Molecular Dynamics Simulations

Lignin, an essential building block of lignocellulosic biomass, is a potential abundant source of aromatic monomers for the polymer and chemical industry. Reductive catalytic fractionation (RCF) is one promising process that can produce high yields of phenolic monomers and oligomers from lignin under different catalytic conditions. An important choice in optimizing RCF is the selection of solvent; however, detailed insights into solvent effects on lignin behaviors and interactions remain limited. Here, in this work, we perform all-atom molecular dynamics simulations to study the solvation of lignin, solvent-mediated conformational changes, and the interaction of solvated lignin oligomers with model surfaces. We focus on the behavior of an oligomeric lignin model compound in methanol, ethanol, a binary mixture of ethanol and water, and water at both the RCF reaction temperature (473 K) and room temperature. Analysis of structural features of lignin suggests that these three organic solvent systems favorably solvate lignin, resulting in a more extended conformation suitable for catalytic conversion to valuable chemicals. We further introduce model palladium (Pd) and carbon (C) surfaces to understand how solvent choice impacts adsorption onto a representative catalytic surface and support, and to quantify the competition among the reactant and solvent molecules for the surface. Unbiased simulations suggest strong adsorption of lignin on both Pd and C surfaces at 473 K, with notable solvent-mediated differences in adsorption energies. Additionally, our findings indicate that lignin adsorption is promoted by the entropy change resulting from the displacement of solvent molecules from the surface. This study provides a molecular perspective of adsorption of lignin onto varying surfaces, which is a step towards understanding and optimizing the catalytic conversion of lignin into valuable chemicals.

adsorption↗
NASA NTRS1980

A model of comet comae. II - Effects of solar photodissociative ionization

Improvements to a computer model of coma photochemistry are described. These include an expansion of the chemical reactions network and new rate constants that have been measured only recently. Photolytic reactions of additional molecules are incorporated, and photolytic branching ratios are treated in far greater detail than previously. A total of 25 photodissociative ionization (PDI) reactions are now considered (as compared to only 3 PDI reactions previously). Solar PDI of the mother molecule CO2 is shown to compete effectively with photoionization of CO in the production of observed CO(+). The CO(+) density peak predicted by the improved model, for CO2 or CO mother molecules, is deep in the inner coma, in better agreement with observation than the old CO2 model. However, neither CO2 nor CO mother molecule calculations reproduce the CO(+)/H2O(+) ratio observed in comet Kohoutek. PDI products of CO2, CO, CH4, and NH3 mother molecules fuel a complex chemistry scheme, producing inner coma abundances of CN, C2, and C3 much greater than previously calculated.

Huebner, W. F.↗
DOE OSTI2025

Data and code from: Multivariate bayesian regression model for predicting disposed ash composition at U.S. coal fired power stations

This dataset contains the code and data files needed for implementation of a Multivariate Bayesian Regression model, described in Jin et al. (2025), for the historical prediction of the chemical composition of disposed coal ash at U.S. coal fired power plants as a function of annualized coal purchase data. The integrated coal supply data file (CoalSupplyDataset.csv) represents a compilation of monthly fuel purchase records for the period 1973-2022 at major U.S. power stations. These records were obtained from the U.S. Energy Information Administration. The CSV file also contains, for each coal purchase record, the coal region of the mine as defined by the U.S. Geological Survey. Data entry errors and data gaps in the EIA records were corrected as described in Jin et al. This CSV file represents the integrated coal supply data after corrections were made. The model structure and fitting parameters are encoded in pickle file format (Bayesian.pkl). The model was developed with the coal supply data and coal ash composition data, apportioned according to the Stratified Shuffle Split for training and testing subsets. The model was built using Python and the PyMC library. Reference Publication: Jin, Z.; Huang, J.; Hower, J.C.; Hsu-Kim, H.(2025). Predictive Assessment of the Chemical Composition of Coal Ash in Reserve at U.S. Disposal Sites. Environmental Science & Technology.

Coal ash composition↗
DOE OSTI2026

Cathodic Protection Modeling for Hanford Underground Double-Shell Tank Farms

Hanford stores millions of gallons of radioactive and chemically hazardous waste from the production of weapon materials in tank farms consisting of underground carbon-steel storage tanks surrounded by reinforced concrete. Six of these Hanford tank farms use double-shell storage tanks (DSTs). The DST farms were constructed from 1968 to 1986 with a planned 40–50 year design life, so some are already operating beyond their initial life expectancy. Ultrasonic testing (UT) has indicated significant thinning on the bottom of the secondary (outer) liner of these tanks, believed to arise from groundwater intrusion driving concrete side corrosion. There is no direct access to the steel/concrete interface between the tank and the concrete pad, making it difficult to apply a chemical-based mitigation strategy or to conduct repairs, but cathodic protection (CP) is a possible method to inhibit further concrete-side corrosion. Hanford already uses CP to protect below grade steel piping within the tank farms and connected to the tanks, but this system was not designed to protect the tank bottoms. CP design must account for the structures surrounding the DSTs, including the steel reinforcing bars (rebar) within the concrete pad and vault, various process lines, and the existing CP system. In this study, finite element analysis (FEA) modeling was carried out to simulate CP protection of 1) a single tank and CP anode to develop options for modeling the rebar and to compare to a simpler circuit model and 2) the entire Hanford AN tank farm as a representative example consisting of seven tanks, associated piping, and both existing and new CP anodes. Both circuit and FEA models predict that significant protective current could be delivered to the bottoms of the tanks with the addition of tank-protection anodes below the depth of the tanks. Simulations with only the existing pipe-protection anodes active confirmed that only a very small current to the tank bottoms is predicted under present conditions. Multiple simplified representations of the dome and wall rebar were tested to reduce the computational complexity of the tank-farm simulations, resulting in modeling the rebar as edge elements with a prescribed effective circumference that matches the real rebar surface area. The geometry of the rebar is also simplified into horizontal hoops around the tank walls and radial rebar over the dome with increased effective circumference to retain the target surface area. This simplification was found to greatly reduce the complexity and solution time of the models without large changes in current distributions, especially to the tank bottom. A range of values were tested for model parameters such as soil and concrete resistivities and polarization resistance to investigate their impact on the current and electric potential distributions. Depending on the parameters used, FEA simulations predict some risk of overprotection, particularly on the piping system; since overprotection can also lead to surface damage associated with hydrogen gas generation at the interface (e.g. hydrogen embrittlement or damage to coatings), this needs to be considered when refining the design of the new CP system. Comparison between the FEA models and the circuit model representation demonstrated that the circuit model could not match the predicted FEA current distribution, even when using the exact same surface areas. This discrepancy appeared to be at least partly attributable to the impact of the relative positions of the tank components and anodes to each other and to the ground surface. The FEA model accounts for the relative positions since it solves the governing equations in three dimensions, but the circuit model cannot account for the positioning. In particular, the circuit model underpredicts the current to the tank bottom and overpredicts the current to the dome compared to FEA for the baseline geometry. The FEA models omitted the electrically isolated rebar in the bottom concrete slab. However, a circuit based stray current model estimated that only 2.1% of the total current through the slab would stray into the rebar, corresponding to ~0.21 A for a target current density of 2 mA/ft2 to the tank bottom. The estimated corrosion driven by this amount of stray current is predicted to yield a lifetime of >400 years for the minimum rebar diameter, assuming an acceptable cross-section area loss of 10%.

d'Entremont, Anna [Savannah River National Laborat↗
DOE OSTI2026

Hierarchical Multi-agent Large Language Model Reasoning for Autonomous Heterogeneous Catalyst Discovery

Artificial intelligence is reshaping scientific exploration, but most methods automate procedural tasks without engaging in scientific reasoning, limiting autonomy in discovery. We demonstrate that hierarchical agentic large language model reasoning can efficiently drive simulation and scientific exploration. Across two chemical applications, CO adsorption on Cu surface transition metal adatoms and on M–N–C catalysts, reasoning-guided exploration reduces required atomistic simulations by up to 90% relative to heuristic or random selection. Comparisons across single-agent, multi-agent, and stochastic baselines show that hierarchical strategies yield more coherent and information-efficient search trajectories. Reasoning traces reveal chemically grounded decisions that cannot be explained by semantic bias or stochastic sampling. We realize these agentic reasoning strategies in Materials Agents for Simulation and Theory in Electronic-structure Reasoning (MASTER), a multimodal system that translates natural language into density functional theory workflows. Altogether, multi-agent collaboration accelerates heterogeneous catalyst discovery and marks a step toward more autonomous, reasoning-guided scientific exploration.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

TChem-atm (v2.0.0): scalable performance-portable multiphase atmospheric chemistry

We present TChem-atm, a performance-portable approach that enables efficient simulation of chemically detailed and multiphase atmospheric chemistry on modern heterogeneous computing architectures. Unlike previous efforts that rely on architecture-specific code or focus exclusively on gas-phase chemistry, TChem-atm supports fully coupled gas–aerosol systems with execution across CPUs, NVIDIA GPUs, and AMD GPUs through the Kokkos programming model. It integrates the flexible multiphase capabilities of the Community Atmospheric Model Chemistry Package (CAMP) with the high-performance kinetic routines of TChem, and includes automatic Jacobian construction with support for a range of stiff ODE solvers. In a proof-of-concept integration with the particle-resolved model PartMC, TChem-atm reproduces the existing PartMC–CAMP implementation within solver tolerances and delivers substantial GPU speedups, especially for large particle populations. Performance benchmarks reveal substantial speedups on GPU platforms, particularly for large particle populations, with consistent results across hardware backends. TChem-atm enables performance-portable execution across CPUs and GPUs, though optimal efficiency may require modest architecture-specific tuning (e.g., team and vector sizes), with up to a twofold improvement on the NVIDIA H100. It directly supports sectional and particle-resolved host models, while modal aerosol schemes require minor adaptation to provide particle-scale quantities such as representative diameters. By enabling chemically detailed, multiphase simulations with performance portability and host-model flexibility, TChem-atm facilitates the incorporation of advanced chemistry into atmospheric models.

Díaz-Ibarra, Oscar Homero [Sandia National Laborat↗
DOE OSTI2026

Opportunities and challenges in process modeling and simulation of electrochemical systems

Electrochemical technologies have garnered intense interest in both academic research and industrial applications due to their potential to increase energy efficiency and reduce carbon footprint. However, electrochemical process fundamentals have been absent from the chemical process simulators available to millions of chemical engineers worldwide. To expedite process research and development of electrochemical technologies in the chemical industry, it is imperative that essential electrochemical process fundamentals be incorporated into process simulators to support modeling and simulation of electrochemical processes. Here, this study examines three process fundamentals key to the research and development of electrochemical processes: electrochemical reaction kinetics, electrolyte thermodynamics, and heat and mass transfer. It further illustrates application of these process fundamentals with a case study modeling an electrochemical process for the conversion of acrylonitrile to adiponitrile. The modeling example highlights the roles of applied voltage on reaction rates, and electron flow rates on the performance of electrochemical processes. It suggests that the applied voltage and electron flow rates are two unique concepts that should be included in simulators for electrochemical processes.

09 BIOMASS FUELS↗
NASA NTRS1981

Atmospheric structure of the outer planets from thermal emission data

Methods for determining atmospheric structure exploit the opacities provided by the collision induced H2 dipole and the nu4 fundamental of CH4. In addition to earth-based observations, useful measurements of thermal emission from Jupiter and Saturn have been or soon will be made by several spacecraft, with results cross-checked with independent radio occultaion results. For Uranus and Neptune, only a limited set of whole-disk earth-based data exists. All the outer planets show evidence for stratospheric temperature inversions; temperature minima range from about 105 K for Jupiter and 87 K for Saturn, to roughly 55 K for Uranus and Neptune. Remaining problems may be resolved by better quantitative understanding of gas and aerosol absorption and scattering properties, chemical composition, and non-LTE source functions. Ultimately, temperature structure results must be supplemented by quantitative energy equilibrium models which will allow some meaning to be given to the relationships between such characteristics as temperature, clouds, incident solar and planetary radiation, and chemical composition.

Orton, G. S.↗
NASA NTRS1972

The fluid mechanics of thrombus formation

Experimental data are presented for the growth of thrombi (blood clots) in a stagnation point flow of fresh blood. Thrombus shape, size and structure are shown to depend on local flow conditions. The evolution of a thrombus is described in terms of a physical model that includes platelet diffusion, a platelet aggregation mechanism, and diffusion and convection of the chemical species responsible for aggregation. Diffusion-controlled and convection-controlled regimes are defined by flow parameters and thrombus location, and the characteristic growth pattern in each regime is explained. Quantitative comparisons with an approximate theoretical model are presented, and a more general model is formulated.

Source record↗
NASA NTRS1984

Nitrogen on Triton

The tentative identification of liquid nitrogen as the source of the 2.16-micron absorption band in the spectrum of Triton has suggested a surface model in which nitrogen, together with the previously identified methane, occur in some kind of chemical association. Water frost may be another surface chemistry component, and dark, photochemically-produced particles occur as a minor contaminant of the ices or liquid. Temperature is noted to be a critical factor in the surface model, because the calculated subsolar surface temperature of Triton lies near nitrogen's melting-freezing temperature. The present modeling of Triton suggests nitrogen's occurrence in the form of fine grains.

Clark, R. N.↗
DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗
DOE OSTI2025

Interface Nucleus Templating of Modular Intermetallic Morphologies: Chemical Pressure Complementarity, Columnar Domains, and Complex Disorder in Y 13 Ag 42.7 Zn 29.7

One element of the diversity of intermetallic phases is the formation of complex structures from the assembly of fragments of simpler structures. Recently, we devised the Interface Nuclear Approach as a model for understanding such modular arrangements, in which the intergrowth of different structures is driven by chemical pressure (CP) relief at shared motifs at the domain interfaces, referred to as interface nuclei. In this Article, we present the synthesis, crystal structure, and CP analysis of a new compound that expands on this theme, Y 13 Ag 42.7 Zn 29.7 . Its hexagonal structure contains interpenetrating domains based on the CaPd 5+x and EuMg 5 types. The CaPd 5+x -based regions are reminiscent of the lamellar intergrowth structures previously observed in the Y−Ag−Zn system. In Y 13 Ag 42.7 Zn 29.7 , however, the domains have a different morphology, forming columns that adopt a hexagonal rod-packing. The geometrical features of the remaining spaces are assigned, using the program GrowDomain, to the cores of trigonal units of the EuMg 5 type, while layers of disordered atoms occur at heights along z where the parent structures are mismatched. At the CaPd 5+x -type/EuMg 5 -type interfaces, simple interface nucleus motifs with strong CP-complementarity can be identified, while their distribution within the parent structures supports the notion of templated architectures in modular intermetallics.

Chemical structure↗
NASA NTRS1976

Chemical mapping of planetary surfaces

Two instruments, the gamma-ray spectrometer and the X-ray fluorescence spectrometer, are uniquely suited to the chemical mapping of planetary surfaces from orbit. Through their detection of characteristic line spectra they measure the concentrations of a suite of elements in each area overflown. Multielement chemical maps derived from these remote measurements are used in the construction of evolutionary models of planetary bodies and of the solar system as a whole. The NaI(T1) gamma-ray spectrometer and a gas proportional X-ray spectrometer were flown over 20 percent of the lunar surface during the Apollo 15 and 16 missions. These instruments measured chemical differences across the boundaries of known lunar provinces and revealed several new features of lunar-surface composition. Advanced spectrometers which are under development for future missions are able to educe much more information in a given time span than the Apollo instruments. They may be used in possible future missions such as Lunar Polar Orbiter, a Mars orbiter, a Mercury orbiter, outer planet satellite missions, rendezvous with asteroids and cometary nuclei, and surface-penetrating planetary probes.

Haines, E. L.↗
DOE OSTI2025

Galactic Chemical Evolution of Short-lived Radioactive Isotopes Produced by Explosive Nucleosynthesis: 60 Fe and 53 Mn

Several short-lived radionuclides (SLRs) are know to have existed in the early solar system (ESS). These species, which typically decay with half-lives of the order of a few million years, can be used to probe the timescales of events preceding the birth of the Sun. We investigate the ESS origin of 53 Mn, produced by core-collapse supernovae (CCSNe) and Type Ia supernovae (SNe Ia), and 60 Fe, produced exclusively by CCSNe. We model the evolution of the radioactive-to-stable abundance ratios of these SLRs with a Galactic chemical evolution (GCE) framework accounting for different SN yields, SN Ia delay times, and other Galactic features (K). A further set of models is calculated assuming that SNe Ia did not contribute any 53 Mn to the ESS. The predicted ratios are compared to meteoritic ratios to derive a distribution of solar isolation times that includes uncertainties due to stochastic chemical enrichment and measurements of the ESS values. The isolation times are then compared to those of 107 Pd and 182 Hf calculated in previous work. A self-consistent solution can be found within the current uncertainties, especially when using the GCE setups with K = 1.6 and 2.3, although the maximum likelihood for the 60 Fe distribution is typically ∼4–5 Myr shorter than for 53 Mn. The predicted 60 Fe/ 53 Mn ratio, instead, is completely inconsistent with the ESS value; this could be resolved using a larger fraction of faint CCSNe than usually considered in GCE models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

High-Fidelity Multiphysics Assessment of the Molten Chloride Reactor Experiment (MCRE) Integrating Neutronics, Thermal-hydraulics and Thermochemistry

This work presents a MOOSE-based Multiphysics model of the LOTUS MSR reactor system. The model includes the coarse-mesh turbulent thermal-hydraulics including passive advection of delayed neutron precursors of the primary system through Pronghorn, the power density and neutronics calculation employing Griffin, and the chemical interactions through Gibbs energy minimization and redox potential with Thermochimica. Incorporating a plate-out model, the full multi-physics model is used to evaluate the effects of deposition and removal of solid compounds formed in the molten salt in the reactor surfaces, evaluating the concentrations and regions in which nickel, chromium and iron would be deposited. Results indicate an almost isothermal state of the fuel during operation at 25 kW that do not increase the plate-out effects due to low temperature gradients.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗