Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation decomposition”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Characterizing the Effects of Thermal Profile and Gas Environment on the Heat Treatment of JSC-1A Lunar Simulant Using Combined Thermogravimetry and Mass Spectroscopy

As NASA strives towards a long duration presence on the Moon, it has become increasingly important to learn how to better utilize resources from the lunar surface for everything from habitats, vehicle infrastructure, and chemical extraction. To that end, a variety of lunar simulants have been sourced from terrestrially available volcanic minerals and glass as Apollo regolith is a finite resource unavailable for experimentation needing large masses. However, while mineralogy and chemical composition can approach that of lunar material in these simulants, there are still distinct non-lunar phases such as hydrates, carbonates, sulfates, and clays that can cause simulants to behave distinctly non-lunar in a variety of processing conditions that would ideally be applied to lunar material. Notably, severe glassy bubbling has been documented in a variety of vacuum sintering experiments on JSC-1A lunar Maria simulant heated via microwaves. The origins of this outgassing have not been well understood, but are normally attributed to the decomposition of non-lunar contaminates intrinsic to virtually all terrestrially sourced simulants. As such, a series of controlled environmental tests were performed to ascertain the origins of the high temperature outgassing and to develop heat treatments that can drive JSC-1A closer to lunar composition and behavior. It was found that in JSC-1A that distinct gas evolutions of water, carbon dioxide, and sulfur dioxide occur in both inert gas and vacuum. Additionally, the presence of hydrogen during heat treatments was shown to dramatically change gas evolutions, leading to distinctly more lunar-like composition and behavior from JSC-1A simulant.

sinter↗
NASA NTRSDate not supplied

Characterizing the Effects of Thermal Profile and Gas Environment on the Heat Treatment of JSC-1A Lunar Simulant Using Combined Thermogravimetry and Mass Spectroscopy

As NASA strives towards a long duration presence on the Moon, it has become increasingly important to learn how to better utilize resources from the lunar surface for everything from habitats, vehicle infrastructure, and chemical extraction. To that end, a variety of lunar simulants have been sourced from terrestrially available volcanic minerals and glass as Apollo regolith is a finite resource unavailable for experimentation needing large masses. However, while mineralogy and chemical composition can approach that of lunar material in these simulants, there are still distinct non-lunar phases such as hydrates, carbonates, sulfates, and clays that can cause simulants to behave distinctly non-lunar in a variety of processing conditions that would ideally be applied to lunar material. Notably, severe glassy bubbling has been documented in a variety of vacuum sintering experiments on JSC-1A lunar Maria simulant heated via microwaves. The origins of this outgassing have not been well understood, but are normally attributed to the decomposition of non-lunar contaminates intrinsic to virtually all terrestrially sourced simulants. As such, a series of controlled environmental tests were performed to ascertain the origins of the high temperature outgassing and to develop heat treatments that can drive JSC-1A closer to lunar composition and behavior. It was found that in JSC-1A that distinct gas evolutions of water, carbon dioxide, and sulfur dioxide occur in both inert gas and vacuum. Additionally, the presence of hydrogen during heat treatments was shown to dramatically change gas evolutions, leading to distinctly more lunar-like composition and behavior from JSC-1A simulant.

sinter↗
DOE OSTI2022

Cleavage of an aromatic ring and radical migration

The present study undertakes a theoretical evaluation of thermal decomposition of aromatic-ring radicals. Potential energy surfaces and associated reaction rate coefficients were calculated for 1- and 2-naphthalenyl, acetanaphthylenyl, and pyrenyl radicals. Kinetic Monte-Carlo simulations were performed to examine the rates of ring cleavage in two sooting laminar premixed flames of ethylene. The simulations showed that the thermal aromatic-ring cleavage is comparable in rate to oxyradical decomposition in a heavier-sooting flame. The simulation also revealed, unexpectedly, fast internal ring radical migration, comparable in frequency to reaction events of aromatic growth.

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

A COUPLED DETERMINISTIC TRANSPORT CALCULATION IN COMSOL USING PROPER ORTHOGONAL DECOMPOSITION

A coupled physics problem consisting of radiation transport and heat transfer was modeled using COMSOL Multiphysics to simulate a criticality accident. Reduced order models (ROM) were created to accelerate the radiation transport calculation through Proper Orthogonal Decomposition (POD). Additionally, on-the-fly neutron cross sections were generated from a non-linear function through a feed forward neural network. The Livelink for MATLAB module allowed for data transfer between an otherwise independent solver and in-house modules. A data driven approach to POD was developed to act as a surrogate model to couple with heat transfer, with the POD module written in MATLAB coupled to an isolated heat transfer solve in COMSOL. A coupled transient solution resulted in a good approximation of the flux and temperature, with the maximum error being 3E-2 and 2.5E-3 respectively.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1986

Simulator for multilevel optimization research

A computer program designed to simulate and improve multilevel optimization techniques is described. By using simple analytic functions to represent complex engineering analyses, the simulator can generate and test a large variety of multilevel decomposition strategies in a relatively short time. This type of research is an essential step toward routine optimization of large aerospace systems. The paper discusses the types of optimization problems handled by the simulator and gives input and output listings and plots for a sample problem. It also describes multilevel implementation techniques which have value beyond the present computer program. Thus, this document serves as a user's manual for the simulator and as a guide for building future multilevel optimization applications.

Padula, S. L.↗
NASA NTRS2021

A Phenomenon-Based Decomposition of Model-Based Estimates of Boreal Winter ENSO Variability

Climate models are now routinely being used to simulate and predict climate variability on time scales ranging from sub-seasonal to seasonal and longer. As such, there are now long histories of such simulations and predictions spanning multiple decades and multiple ensemble members, both of which are crucial for separating climate signal from climate noise. A key focus of such runs has been the El Niño-Southern Oscillation (ENSO), spurred by recent improvements in our ability to predict such events, though questions remain as to how well climate models do beyond simply always predicting the “canonical” atmospheric response to an ENSO event—something simple statistical models already do reasonably well. This is a critical issue that needs addressing, given the importance of event-to-event differences for predicting regional impacts of ENSO teleconnections, and the need to justify the expense of running sophisticated climate models. Unfortunately, current diagnostic tools are not well suited for quantifying the different sources of variability associated with specific phenomena such as ENSO. More generally, while much effort has focused on addressing model bias, less has been done to address errors in second moment statistics—an issue whose importance is gaining increased attention particularly as we build climate prediction systems capable of taking advantage of forecasts of opportunity—a capability that requires reliable estimates of forecast uncertainty. In this report, we outline a phenomenon-based statistical decomposition of climate variance(in essence a detailed variance budget)that is specifically tailored to address the above questions by separating the variability (both the signal and noise) into that tied to the long-term average impact of a particular phenomenon(the composite mean) and the event-to-event(E2E) variability about the composite mean. In addition, we provide related decompositions of the correlations that allow us to quantify how much of the agreement with observations (the skill) comes from the composite mean as opposed to from the E2Evariability. As an example, we present the results of such a decomposition for ENSO based on simulations with the GEOS atmospheric general circulation model (AGCM), with a focus on the monthly mean impacts over North America during boreal winter(December –March). Here we take advantage of existing GEOS AGCM simulations that were produced as companion simulations to MERRA-2for the period 1980-2016. Comparisons are made throughout with MERRA-2.

ENSO↗
DOE OSTI2021

Stabilizer Rank Strong Simulation v.0.1

This software improves upon the Clifford gate+T gate magic state strong simulation algorithm by Brayvi and Gosset. It provides a better stabilizer state decomposition that scales as 2^{0.463t}, an improvement over 2^{0.468t}, for t T gate magic states. This code enables simulations of universal quantum computing with that asymptotic scaling cost. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525 SAND2020-14293 O

SciDAC↗
DOE OSTI2020

Decomposition Algorithms for Solving NP-hard Problems on a Quantum Annealer

NP-hard problems such as the maximum clique or minimum vertex cover problems, two of Karp’s 21 NP-hard problems, have several applications in computational chemistry, biochemistry and computer network security. Adiabatic quantum annealers can search for the optimum value of such NP-hard optimization problems, given the problem can be embedded on their hardware. However, this is often not possible due to certain limitations of the hardware connectivity structure of the annealer. This paper studies a general framework for a decomposition algorithm for NP-hard graph problems aiming to identify an optimal set of vertices. Our generic algorithm allows us to recursively divide an instance until the generated subproblems can be embedded on the quantum annealer hardware and subsequently solved. Furthermore, the framework is applied to the maximum clique and minimum vertex cover problems, and we propose several pruning and reduction techniques to speed up the recursive decomposition. The performance of both algorithms is assessed in a detailed simulation study.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimized Lie–Trotter–Suzuki decompositions for two and three non-commuting terms

Lie–Trotter–Suzuki decompositions are an efficient way to approximate operator exponentials exp ( t H ) when H is a sum of n (non-commuting) terms which, individually, can be exponentiated easily. They are employed in time-evolution algorithms for tensor network states, digital quantum simulation protocols, path integral methods like quantum Monte Carlo, and splitting methods for symplectic integrators in classical Hamiltonian systems. Here, we provide optimized decompositions up to order t 6 . The leading error term is expanded in nested commutators (Hall bases) and we minimize the 1-norm of the coefficients. For n = 2 terms, several of the optima we find are close to those in McLachlan (1995). Generally, our results substantially improve over unoptimized decompositions by Forest, Ruth, Yoshida, and Suzuki. We explain why these decompositions are sufficient to efficiently simulate any one- or two-dimensional lattice model with finite-range interactions. This follows by solving a partitioning problem for the interaction graph.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Towards robust surrogate models: Benchmarking machine learning approaches to expediting phase field simulations of brittle fracture

Data-driven approaches have the potential to make modeling complex, nonlinear physical phenomena significantly more computationally tractable. For example, computational modeling of fracture is a core challenge where machine learning techniques have the potential to provide a much needed speedup that would enable progress in areas such as multi-scale modeling and uncertainty quantification. Currently, phase field modeling (PFM) of fracture is one such approach that offers a convenient variational formulation to model crack nucleation, branching and propagation. To date, machine learning techniques have shown promise in approximating PFM simulations. While standard fracture benchmarks represent realistic scenarios frequently observed in practice, they typically do not provide sufficiently challenging tests for data-driven methods. Here, to address this gap, we introduce a challenging dataset based on PFM simulations designed to benchmark and advance ML methods for fracture modeling. This dataset includes three energy decomposition methods, two boundary conditions, and 1000 random initial crack configurations for a total of 6000 simulations. Each sample contains 100 time steps capturing the temporal evolution of the crack field. Alongside this dataset, we also implement and evaluate Physics Informed Neural Networks (PINN), Fourier Neural Operators (FNO), and UNet models as baselines, and explore the impact of ensembling strategies on prediction accuracy. With this combination of our dataset and baseline models drawn from the literature we aim to provide a standardized and challenging benchmark for evaluating machine learning approaches to solid mechanics. Our results highlight both the promise and limitations of popular current models, and demonstrate the utility of this dataset as a testbed for advancing machine learning in fracture mechanics research.

Benchmark dataset↗
DOE OSTI2020

Effects of charged interfaces on electrolyte decomposition at the lithium metal anode

Lithium–Sulfur batteries are promising candidates to substitute conventional Li-ion batteries due to their higher energy density and reduced cost. However, several challenges related to the reactivity of the lithium metal anode have prevented this technology from becoming broadly commercialized. Lithium's high reactive nature leads to the continuous decomposition of the electrolyte and the formation of the solid-electrolyte interphase (SEI) layer. Thus, a comprehensive understanding of how the SEI film is formed is crucial to help elucidate improvements for this battery technology. In this work, we use density functional theory (DFT) based computational methods to investigate the effect of charged interfaces on the electrolyte degradation at the lithium anode. Here, several electrolyte mixtures were considered including 1,2-dimethoxyethane (DME) and 1,3-dioxolane (DOL) as solvents, and lithium bis(fluorosulfonyl)imide (LiFSI) and lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as salts. It is found that the extent of salt decomposition is higher when the interface is charged. In addition, two types of solvent reduction mechanisms are identified: C–O bond cleavage and radical attacks. Charge transfer and evolution analysis are characterized in detail. Finally, solvent decomposition reactions are found to become more thermodynamically favorable and their activation barriers to be diminished under the effect of constant potentials and electric fields.

25 ENERGY STORAGE↗
DOE OSTI2025

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2025

Multi-Scale Modeling for Plasma-Enhanced Ammonia Decomposition over Carbides and Nitrides

Ammonia is a carbon-free hydrogen carrier, but its decomposition typically requires high temperatures over costly Ru-based catalysts due to the large barrier for N≡N bond formation. We develop a multiscale framework combining density functional theory, zero-dimensional plasma kinetics, and microkinetic modeling to elucidate how non-thermal plasma (NTP) enables low-temperature NH 3 decomposition over Co-based carbides and nitrides, benchmarked against Ru and Co. Under thermal conditions, all catalysts are limited by N≡N bond formation, with Co 3 C(001) most active owing to its negatively charged surface, strong N* binding, and low activation barriers of N≡N bond formation. Plasma-induced vibrational excitation of NH 3 and its reactive radicals promotes a radical-driven •NH 2 –N* coupling pathway that dominates on Co 3 C(001) and Co 3 N(001), shifting the rate-limiting step to NH 3 (v1) dissociation, increasing turnover frequencies by up to 6 orders of magnitude, and reducing the temperature needed to reach a turnover frequency of 5 s –1 from >680 °C (Ru and Co under thermal condition) to 267 °C (Co 3 C) and 415 °C (Co 3 N). These results identify Co-based carbides and nitrides as promising plasma-active catalysts for energy-efficient hydrogen production from ammonia.

ammonia decomposition↗
DOE OSTI2023

The Energetic Origins of Pi–Pi Contacts in Proteins

Accurate potential energy models of proteins must describe the many different types of noncovalent interactions that contribute to a protein’s stability and structure. Pi–pi contacts are ubiquitous structural motifs in all proteins, occurring between aromatic and nonaromatic residues and play a nontrivial role in protein folding and in the formation of biomolecular condensates. Guided by a geometric criterion for isolating pi–pi contacts from classical molecular dynamics simulations of proteins, we use quantum mechanical energy decomposition analysis to determine the molecular interactions that stabilize different pi–pi contact motifs. We find that neutral pi–pi interactions in proteins are dominated by Pauli repulsion and London dispersion rather than repulsive quadrupole electrostatics, which is central to the textbook Hunter–Sanders model. This results in a notable lack of variability in the interaction profiles of neutral pi–pi contacts even with extreme changes in the dielectric medium, explaining the prevalence of pi-stacked arrangements in and between proteins. We also find interactions involving pi-containing anions and cations to be extremely malleable, interacting like neutral pi–pi contacts in polar media and like typical ion–pi interactions in nonpolar environments. Like-charged pairs such as arginine–arginine contacts are particularly sensitive to the polarity of their immediate surroundings and exhibit canonical pi–pi stacking behavior only if the interaction is mediated by environmental effects, such as aqueous solvation.

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

Electrochemical formation of bis(fluorosulfonyl)imide-derived solid-electrolyte interphase at Li-metal potential

Lithium bis(fluorosulfonyl)imide-based liquid electrolytes are promising for realizing high Coulombic efficiency and long cycle life for next-generation Li-metal batteries. However, the role of anions in the formation of solid-electrolyte interphase remains unclear. Herein, we combine electrochemical analyses and X-ray photoelectron spectroscopy measurements both with and without sample washing, together with computational simulations, to propose the reaction pathways of electrolyte decomposition and correlate the interphase component solubility with the efficacy of passivation. Additionally, we discovered that not all the products derived from interphase-forming reactions were incorporated into the resulting passivation layer, with a notable portion present in the liquid electrolyte. Further, we found that the high-performance electrolytes can afford a sufficiently passivating interphase with minimized electrolyte decomposition, by incorporating more anion-decomposition products. Overall, this work presents a systematic approach of coupling electrochemical and surface analyses to paint a comprehensive picture of solid-electrolyte interphase formation, while identifying key attributes of high-performance electrolytes for guiding future designs.

Yu, Weilai [Stanford University, CA (United States↗
DOE OSTI2025

Status of GPU capabilities within the Shift Monte Carlo radiation transport code

Shift is a general-purpose Monte Carlo (MC) radiation transport code for fission, fusion, and national security applications. Shift has been adapted to efficiently run on GPUs in order to leverage leadership-class supercomputers. This work presents Shift’s current GPU capabilities. These include core radiation transport capabilities for eigenvalue and fixed-source simulations, and support for non-uniform domain decomposition, Doppler broadening, free-gas elastic scattering, general-purpose geometry, hybrid MC/deterministic transport, and depletion. Transport results demonstrate a 2–5× GPU-to-CPU speedup on a per-node basis for an eigenvalue problem on the Frontier supercomputer and a 28× speedup for a fixed-source problem on the Summit supercomputer.

Biondo, Elliott [ORNL] (ORCID:0000000290881360)↗
DOE OSTI2022

Nonlinear spectral analysis of ion acoustic solitons arising from a streaming charged object using the numerical inverse scattering transform

Observations of nonlinear coherent plasma phenomena in spacecraft and terrestrial experiments often rely on visual identification of solitary modes because nonlinear coherent modes are broadband in a Fourier transform based analysis. We implement an alternative spectral decomposition known as the inverse scattering transform and demonstrate its ability to successfully isolate nonlinear modes on the homogeneous and forced Korteweg–De Vries equation which models these nonlinear modes. We also demonstrate for the first time that this decomposition is useful when forcing is applied. This is because the stable modes generated by localized forcing are similar to the homogeneous solutions where the inverse scattering transform is a rigorous decomposition. This spectral technique is then applied to simulations of ion acoustic waves generated by a 10 mm spherical debris object interacting with the ionospheric plasma orbiting at 2000 km altitude. The algorithm is found to successfully detect the resultant solitons. This demonstrates the feasibility of using this spectral technique as a real time analysis tool for screening spacecraft data for nonlinear solitary modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗