Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reproducible science”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

C3MechLite: An integrated component library of compact kinetic mechanisms for low-carbon, carbon neutral and zero-carbon fuels

Based on our latest detailed chemical reaction mechanism, C3MechV4.0, we have developed two reduced reaction mechanisms—C3MechLite and C3MechCore—targeting C 0 –C 3 chemical species including NH 3 . C3MechLite (61 species), contains a number of species comparable to GRI-Mech (53 species), that can accurately predict the combustion characteristics of hydrogen, carbon monoxide, ammonia, methane, natural gas, nitrogen oxides, and their mixtures for a wide range of conditions. C3MechCore (118 species) targets a more comprehensive range of C 0 –C 3 fuels, including ammonia, methanol, ethanol, and dimethyl ether. Both mechanisms demonstrate predictive accuracy comparable to C3MechV4.0 for the combustion characteristics of the target fuels. C3MechLite is designed with a component library structure, enabling further reduction in mechanism size depending on the fuel(s) of interest for 2D/3D numerical simulations. Various combinations of component libraries were validated, and the average prediction error remains within 1 % compared to C3MechLite. Furthermore, the mechanism was applied to 3D LES simulations of H 2 lifted flames and was confirmed to reproduce flame characteristics with high accuracy. C3MechLite and its component library structure enable high-fidelity and computationally efficient chemical kinetic mechanisms, paving the way for application in more complex combustion simulations.

Ammonia↗
DOE OSTI2022

Characterizing vertical upper ocean temperature structures in the European Arctic through unsupervised machine learning

In-situ observations of subsurface ocean temperatures are, in many regions, inconsistently distributed in time and space. These spatio-temporal inconsistencies in the observational network lead to difficulties in utilizing those observations effectively for ocean model evaluation or understanding larger-scale ocean characteristics. Model accuracy of subsurface ocean characteristics is especially important within regions that contain complex ocean structures. One such region is the European Arctic which not only contains several types of water masses with unique characteristics, but also wintertime sea ice coverage and complex bathymetry. This study presents an unsupervised neural networking technique that can be used in combination with traditional ocean model evaluation techniques to provide additional information on the accuracy of modeled vertical ocean temperature profiles. Self-organizing maps is an unsupervised machine learning technique that we apply to approximately twenty thousand Argo and CTD temperature profiles from 2012 to 2020 in the European Arctic to categorize the observed vertical ocean temperature structures in the top 150 m. The observed ocean profile categories, or neurons, defined by the self-organizing map show strong spatial and temporal dependencies. We then use the neuron weights, or the learned temperature profile structure of each neuron, to validate the spatial and temporal variability of modeled vertical temperature structures. This analysis gives us new insights about the model’s capabilities to reproduce specific vertical structures of the top-most ocean layer within different regions and seasons. Mapping modeled ocean temperature profiles onto the neuron-space of the observationally-defined self organized map highlights the potential of this method to advance our understanding of model deficiencies in that region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An investigation of irreproducibility in maximum likelihood phylogenetic inference

Phylogenetic trees are essential for studying biology, but their reproducibility under identical parameter settings remains unexplored. Here, we find that 3515 (18.11%) IQ-TREE-inferred and 1813 (9.34%) RAxML-NG-inferred maximum likelihood (ML) gene trees are topologically irreproducible when executing two replicates (Run1 and Run2) for each of 19,414 gene alignments in 15 animal, plant, and fungal phylogenomic datasets. Notably, coalescent-based ASTRAL species phylogenies inferred from Run1 and Run2 sets of individual gene trees are topologically irreproducible for 9/15 phylogenomic datasets, whereas concatenation-based phylogenies inferred twice from the same supermatrix are reproducible. Our simulations further show that irreproducible phylogenies are more likely to be incorrect than reproducible phylogenies. These results suggest that a considerable fraction of single-gene ML trees may be irreproducible. Increasing reproducibility in ML inference will benefit from providing analyses’ log files, which contain typically reported parameters (e.g., program, substitution model, number of tree searches) but also typically unreported ones (e.g., random starting seed number, number of threads, processor type).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Existence of La-site antisite defects in LaMO 3 (M = Mn, Fe, and Co) predicted with many-body diffusion quantum Monte Carlo

The properties of LaMO 3 (M: 3d transition metal) perovskite crystals are significantly dependent on point defects, whether introduced accidentally or intentionally. The most studied defects in La-based perovskites are the oxygen vacancies and doping impurities on the La and M sites. Here, we identify that intrinsic antisite defects, the replacement of La by the transition metal, M, can be formed under M-rich and O-poor growth conditions, based on results of an accurate many-body ab initio approach. Our fixed-node diffusion Monte Carlo (FNDMC) calculations of LaMO 3 (M = Mn, Fe, and Co) find that such antisite defects can have low formation energies and are magnetized. Complementary density functional theory (DFT)-based calculations show that Mn antisite defects in LaMnO 3 may cause the p-type electronic conductivity. These features could affect spintronics, redox catalysis, and other broad applications. Our bulk validation studies establish that FNDMC reproduces the antiferromagnetic state of LaMnO 3 , whereas DFT with PBE (Perdew–Burke–Ernzerhof), SCAN (strongly constrained and appropriately normed), and the LDA+U (local density approximation with Coulomb U) functionals all favor ferromagnetic states, at variance with experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Inferring and evaluating satellite-based constraints on NO x emissions estimates in air quality simulations

Satellite observations of tropospheric NO 2 columns can provide top-down observational constraints on emissions estimates of nitrogen oxides (NO x ). Mass-balance-based methods are often applied for this purpose but do not isolate near-surface emissions from those aloft, such as lightning emissions. Here, we introduce an inverse modeling framework that couples satellite chemical data assimilation to a chemical transport model. In the framework, satellite-constrained emissions totals are inferred using model simulations with and without data assimilation in the iterative finite-difference mass-balance method. The approach improves the finite-difference mass-balance inversion by isolating the near-surface emissions increment. We apply the framework to separately estimate lightning and anthropogenic NO x emissions over the Northern Hemisphere for 2019. Using overlapping observations from the Ozone Monitoring Instrument (OMI) and the Tropospheric Monitoring Instrument (TROPOMI), we compare separate NO x emissions inferences from these satellite instruments, as well as the impacts of emissions changes on modeled NO 2 and O 3 . OMI inferences of anthropogenic emissions consistently lead to larger emissions than TROPOMI inferences, attributed to a low bias in TROPOMI NO 2 retrievals. Updated lightning NO x emissions from either satellite improve the chemical transport model's low tropospheric O 3 bias. The combined lighting and anthropogenic emissions updates improve the model's ability to reproduce measured ozone by adjusting natural, long-range, and local pollution contributions. Thus, the framework informs and supports the design of domestic and international control strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

32 examples of LLM applications in materials science and chemistry: towards automation, assistants, agents, and accelerated scientific discovery

Abstract Large language models (LLMs) are reshaping many aspects of materials science and chemistry research, enabling advances in molecular property prediction, materials design, scientific automation, knowledge extraction, and more. Recent developments demonstrate that the latest class of models are able to integrate structured and unstructured data, assist in hypothesis generation, and streamline research workflows. To explore the frontier of LLM capabilities across the research lifecycle, we review applications of LLMs through 32 total projects developed during the second annual LLM hackathon for applications in materials science and chemistry, a global hybrid event. These projects spanned seven key research areas: (1) molecular and material property prediction, (2) molecular and material design, (3) automation and novel interfaces, (4) scientific communication and education, (5) research data management and automation, (6) hypothesis generation and evaluation, and (7) knowledge extraction and reasoning from the scientific literature. Collectively, these applications illustrate how LLMs serve as versatile predictive models, platforms for rapid prototyping of domain-specific tools, and much more. In particular, improvements in both open source and proprietary LLM performance through the addition of reasoning, additional training data, and new techniques have expanded effectiveness, particularly in low-data environments and interdisciplinary research. As LLMs continue to improve, their integration into scientific workflows presents both new opportunities and new challenges, requiring ongoing exploration, continued refinement, and further research to address reliability, interpretability, and reproducibility.

Computer Science↗
DOE OSTI2023

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accelerating discrete dislocation dynamics simulations with graph neural networks

Discrete dislocation dynamics (DDD) is a widely employed computational method to study plasticity at the mesoscale that connects the motion of dislocation lines to the macroscopic response of crystalline materials. However, the computational cost of DDD simulations remains a bottleneck that limits its range of applicability. Here, we introduce a new DDD-GNN framework in which the expensive time-integration of dislocation motion is entirely substituted by a graph neural network (GNN) model trained on DDD trajectories. As a first application, we demonstrate the feasibility and potential of our method on a simple yet relevant model of a dislocation line gliding through an array of obstacles. We show that the DDD-GNN model is stable and reproduces very well unseen ground-truth DDD simulation responses for a range of straining rates and obstacle densities, without the need to explicitly compute nodal forces or dislocation mobilities during time-integration. Our approach opens new promising avenues to accelerate DDD simulations and to incorporate more complex dislocation motion behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Additively manufactured strain sensors for in-pile applications

Accurate, real-time monitoring of strain in fuel, cladding, and structural components of nuclear reactors is critical to better understand radiation induced phenomena during reactor tests and operations. The data provided is crucial to verify physics-based, multiscale modeling and simulation efforts, which aim to shorten the timeline for the development of new nuclear materials. Resistive strain gauges have limited performance during in-pile experiments due to the harsh operating conditions and limited physical space between fuel and cladding components. In this work, aerosol jet printing using silver nanoparticle inks was used to fabricate interdigitated electrode capacitive strain gauges on aluminum alloy 6061 tensile specimens. Here to simulate the temperatures of a traditional light water reactor, the capacitive strain gauges were tested with a mechanical test frame up to 300 °C and compared to commercially available bondable resistive strain gauges. The printed capacitive strain gauges exhibited a gauge factor of 1.0 and showed higher reproducibility and predictability of strain sensing performance than the resistive strain gauge. The results demonstrate the potential of aerosol jet printing to fabricate strain sensors with predictable performance and reduced invasiveness for high-temperature applications with confined spacing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Impact of confined geometries on hopping and trapping of motile bacteria in porous media

We use a random walk particle-tracking (RWPT) approach to elucidate the impact of porous media confinement and cell-cell interactions on bacterial transport. The model employs stochastic alternating motility states consisting of hopping movement and trapping reorientation. The stochastic motility patterns are defined based on direct visualization of individual trajectory data. We validate our model against experimental data, at single-cell resolution, of bacterial E. coli motion in three-dimensional confined porous media. Results show that the model is able to efficiently simulate the spreading dynamics of motile bacteria as it captures the impact of cell-cell interaction and pore confinement, which marks the transition to a late-time subdiffusive regime. Furthermore, the model is able to qualitatively reproduce the observed directional persistence. Lastly, our RWPT model constitutes a meshless simple method which is easy to implement and does not invoke ad hoc assumptions but represents the basis for a multiscale approach to the study of bacterial dispersal in porous systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Understanding Recent Global Hydroclimate Change Using Multivariate Detection and Attribution Techniques and GCM Experiments

The observed slowdown in the rate of global average surface temperature (T) warming during the post-1998 decade (referred to at the time as the “global warming hiatus”) posed a challenge for climate researchers because it was not captured in global climate models projection. Observational and modeling evidence hinted that that leveling of the T trend was linked with the simultaneous appearance of cold surface temperature in the eastern equatorial Pacific that formed a protracted La Niña-like state, a phenomenon known as the Inter-decadal Pacific Oscillation (IPO). It was however not clear whether the concurrent equatorial Pacific cooling was a regional response to the rise in atmospheric greenhouse gas levels or anthropogenic aerosols that is not captured by the climate models or whether it was an unrelated, serendipitous manifestation of a natural phenomenon, internal to the climate system. To resolve this issue, we proposed in this project to develop and apply fingerprinting methods to differentiate between the observed climate response to various external forcing agents and the manifestations of internal climate variability. The project also sought to find out if climate models are able to reproduce the change signals detected in key observational phenomena and if they can be used to investigate the underlying processes. As such we examined changes in global hydroclimate variability and in an outstanding example of regional hydroclimate change, namely the multi-decadal 20 th century change in Sahel rainfall.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Cubic to hexagonal tuning in Fe 2 Mn(Si 1– x Ge x ) Heusler alloys

Here, the competition between the stability of the cubic and hexagonal full Heusler alloys and the implications concerning their magnetic properties were systematically studied through the detailed structural and magnetic characterization of the Fe 2 Mn(Si 1– x Ge x ) system. This system was specifically chosen as the parent compositions are cubic ( x = 0) and hexagonal ( x = 1). It is found that the formation of hexagonal phases occurs for the x ≥ 0.6 samples, whereas its phase fraction monotonically increases with x until the pure hexagonal Fe 2 MnGe is formed. The change in structure results in high sensitiveness of both the saturation of magnetization ($M_S$) and Curie temperature ($T_C$) with x values, related to a strong magnetocrystalline anisotropy of the hexagonal phase. Both cubic and hexagonal magnetic features were qualitatively reproduced by Density Functional Theory (DFT) calculations. This work provides an experimental and theoretical foundation for further design of Heusler systems with controlled structures and magnetic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Controlling Structure and Properties of Vapor-Deposited Glasses of Organic Semiconductors: Recent Advances and Challenges

The past decade has seen great progress in manipulating the structure of vapor-deposited glasses of organic semiconductors. Upon varying the substrate temperature during deposition, glasses with a wide range of density and molecular orientation can be prepared from a given molecule. Here, we review recent studies that show the structure of vapor-deposited glasses can be tuned to significantly improve the external quantum efficiency and lifetime of organic light-emitting diodes (OLEDs). We highlight the ability of molecular simulations to reproduce experimentally observed structures, setting the stage for in silico design of vapor-deposited glasses in the coming decade. Finally, we identify research opportunities for improving the properties of organic semiconductors by controlling the structure of vapor-deposited glasses.

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

Performance and Durability of Pure-Water-Fed Anion Exchange Membrane Electrolyzers Using Baseline Materials and Operation

Water electrolysis powered by renewable electricity produces green hydrogen and oxygen gas, which can be used for energy, fertilizer, and industrial applications and thus displace fossil fuels. Pure-water anion-exchange-membrane (AEM) electrolyzers in principle offer the advantages of commercialized proton-exchange-membrane systems (high current density, low cross over, output gas compression, etc.) while enabling the use of less-expensive steel components and nonprecious metal catalysts. AEM electrolyzer research and development, however, has been limited by the lack of broadly accessible materials that provide consistent cell performance, making it difficult to compare results across studies. Further, even when the same materials are used, different pretreatments and electrochemical analysis techniques can produce different results. Here, we report an AEM electrolyzer comprising commercially available catalysts, membrane, ionomer, and gas-diffusion layers operating near 1.9 V at 1 A cm –2 in pure water. After the initial break in, the performance degraded by 0.67 mV h –1 at 0.5 A cm –2 at 55 °C. We detail the key preparation, assembly, and operation techniques employed and show further performance improvements using advanced materials as a proof-of-concept for future AEM-electrolyzer development. Here, the data thus provide an easily reproducible and comparatively high-performance baseline that can be used by other laboratories to calibrate the performance of improved cell components, nonprecious metal oxygen evolution, and hydrogen evolution catalysts and learn how to mitigate degradation pathways.

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

Bandgap tuning of two-dimensional materials by sphere diameter engineering

Developing a precise and reproducible bandgap tuning method that enables tailored design of materials is of crucial importance for optoelectronic devices. Towards this end, we report a sphere diameter engineering (SDE) technique to manipulate the bandgap of two-dimensional (2D) materials. A one-to-one correspondence with an ideal linear working curve is established between the bandgap of MoS 2 and the sphere diameter in a continuous range as large as 360 meV. Fully uniform bandgap tuning of all the as-grown MoS 2 crystals is realized due to the isotropic characteristic of the sphere. More intriguingly, both a decrease and an increase of the bandgap can be achieved by constructing a positive or negative curvature. By fusing individual spheres in the melted state, post-synthesis bandgap adjustment of the supported 2D materials can be realized. In conclusion, this SDE technique, showing good precision, uniformity and reproducibility with high efficiency, may further accelerate the potential applications of 2D materials.

36 MATERIALS SCIENCE↗