Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials simulations”

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 OSTI2022

Efficient Hierarchical State Vector Simulation of Quantum Circuits via Acyclic Graph Partitioning

Early but promising results in quantum computing have been enabled by the concurrent development of quantum algorithms, devices, and materials. Classical simulation of quantum programs has enabled the design and analysis of algorithms and implementation strategies targeting current and anticipated quantum device architectures. In this paper, we present a graph-based approach to achieve efficient quantum circuit simulation. Our approach involves partitioning the graph representation of a given quantum circuit into sub-graphs/circuits that exhibit better data locality. Simulation of each sub-circuit is organized hierarchically, with the iterative construction and simulation of smaller state vectors, improving overall performance. Also, this partitioning reduces the number of passes through data, improving the total computation time. We present three partitioning strategies and observe that acyclic graph partitioning typically results in the best time-to-solution. In contrast, other strategies reduce the partitioning time at the expense of potentially increased simulation times. Experimental evaluation demonstrates the effectiveness of our approach.

Fang, Bo↗
DOE OSTI2025

Insulation Activation and Contamination in Low-Power Experiments

This is a poster intended to be presented alongside a paper under the same name (INL/CON-24-82555 Rev:001) which displays some discussion and results from the paper. In short, an investigation of various industrial-grade high-temperature insulation materials were simulated to have been irradiated under conditions representative of a high-temperature nuclear experiment in order to identify materials with low activation potential.

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

NEML2 - THE NEW ENGINEERING MATERIAL MODEL LIBRARY, VERSION 2

SF-22-157 NEML2 is an extension of an earlier open source code developed at ANL (SF-xxxx) called NEML. These libraries provide tools for constructing and using models of engineering materials, with a focus on high temperature materials for nuclear reactor applications. The software include tools for setting up material models, calibrating models, storing models for later use in other simulations, and running simple "material point" simulations of the material behavior. The main difference between NEML and NEML2 is the new software has been complete restructured to run on GPUs and other accelerators in order to improve the performance of methods used to train constitutive models against data and to improve the performance of the final models in production environments, like finite element software. While over time we expect NEML2 to become more-or-less feature complete with NEML we are filing a new software disclosure because (1) the library has been complete restructured and (2) we want to host NEML2 in a separate repository.

MESSNER, MARK↗
DOE OSTI2025

Large scale polymer toughening of two-dimensional materials revealed by in situ TEM fracture tests and multiscale simulations

Two-dimensional (2D) materials offer significant potential for applications in energy-harvesting devices, batteries, sensors, and transistors. However, their intrinsic brittleness makes them prone to mechanical failure, limiting their practical use. In this work, we perform in situ transmission electron microscopy (TEM) fracture tests on monolayer MoSe2 and uncover an extrinsic toughening effect induced by an ultrathin adsorbed polystyrene adlayer. This adlayer substantially enhances the fracture resistance of the 2D flakes. Through a combination of molecular dynamics simulations and finite element analysis, we elucidate the molecular mechanism behind this toughening effect. Further, it arises from the active crack-bridging behavior of entangled polymer chains and the formation of a fracture process zone that stabilizes crack propagation and increases the energy required for crack extension. The proposed toughening mechanism offers a pathway to improving the mechanical reliability of 2D material-based devices by mitigating the risk of sudden failure.

2D materials↗
DOE OSTI2024

Robust training of machine learning interatomic potentials with dimensionality reduction and stratified sampling

Abstract Machine learning interatomic potentials (MLIPs) enable accurate simulations of materials at scales beyond that accessible by ab initio methods and play an increasingly important role in the study and design of materials. However, MLIPs are only as accurate and robust as the data on which they are trained. Here, we present DImensionality-Reduced Encoded Clusters with sTratified (DIRECT) sampling as an approach to select a robust training set of structures from a large and complex configuration space. By applying DIRECT sampling on the Materials Project relaxation trajectories dataset with over one million structures and 89 elements, we develop an improved materials 3-body graph network (M3GNet) universal potential that extrapolates more reliably to unseen structures. We further show that molecular dynamics (MD) simulations with the M3GNet universal potential can be used instead of expensive ab initio MD to rapidly create a large configuration space for target systems. We combined this scheme with DIRECT sampling to develop a reliable moment tensor potential for titanium hydrides without the need for iterative augmentation of training structures. This work paves the way for robust high-throughput development of MLIPs across any compositional complexity.

Qi, Ji (ORCID:0000000158089931)↗
DOE OSTI2020

Literature Review on the Theory and Measurement of Dry Sliding Friction of Metals

Understanding the friction behavior between two sliding bodies can inform the design of machines, processing of materials, and simulation of dynamic processes. Kinetic friction is a complex phenomenon that depends on a multitude of factors such as sliding velocity, normal force, contact area, surface roughness, material properties, lubrication conditions, and thermal effects. This literature review covers the major known effects of sliding velocity, normal load, and surface roughness on the measured kinetic friction coefficient in the context of microscopic friction phenomena. Classic macroscale friction models are reviewed to illustrate approaches for simulating friction behavior. Prominent experimental friction setups within the literature are discussed with respect to achievable velocity and pressure regimes. The background information gathered here will be used to inform experimental procedures and modeling strategies of the exploratory research (ER) project titled “Measurement of Dynamic Friction via Kolsky Bar” (20200418ER).

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum embedding theories to simulate condensed systems on quantum computers.

Quantum computers hold promise to improve the efficiency of quantum simulations of materials and to enable the investigation of systems and properties that are more complex than tractable at present on classical architectures. Here, we discuss computational frameworks to carry out electronic structure calculations of solids on noisy intermediate-scale quantum computers using embedding theories, and we give examples for a specific class of materials, that is, solid materials hosting spin defects. These are promising systems to build future quantum technologies, such as quantum computers, quantum sensors and quantum communication devices. Although quantum simulations on quantum architectures are in their infancy, promising results for realistic systems appear to be within reach.

Vorwerk, Christian↗
DOE OSTI2021

Surface prediction and measurement for modulated tool path (MTP) turning

Here, this paper describes a time-domain simulation for predicting surface finish in modulated tool path (MTP) turning, which uses sinusoidal axis motions in the feed direction to produce discontinuous chips for ductile workpiece materials. The simulation includes: the low frequency and low amplitude tool oscillation in the feed direction; the time-varying chip thickness, cutting force, and tool displacement; and the plastic side flow effect used to calibrate the effective tool nose radius. Comparisons between predicted and measured surface profiles are presented for turning a 6061-T6 aluminum cylinder as a function of the MTP oscillation frequency and amplitude with discontinuous chip formation.

42 ENGINEERING↗
DOE OSTI2020

Simulation Analysis of Geometry and Material Effects for Dropkinson Bar

The reported research is motivated by the need to address a key issue affecting the Dropkinson bar apparatus. This unresolved issue is the interference of the stress wave reflected from the bar-beam boundary with the measurement of the stress-strain response of a material tested in the apparatus. The purpose of the wave beam that is currently connected to the bar is to dissipate the stress wave, but the portion of the wave reflected from the bar-beam boundary is still significant. First, we focused on understanding which parameters affect the reflected wave's arrival time at a strain gauge. Specifically, we used finite-element numerical simulations with the Sierra/SM module to study the effects of various bar-beam connection fixities, alternative wave beam materials, and alternative geometries of the Dropkinson bar system based on a monolithic design. The conclusion of this study is that a partial reflection always occurs at the bar-beam boundary (or, for a monolithic design, at a point where the bar geometry changes). Therefore, given a fixed total length of the bar, it is impossible to increase the reflected wave's arrival time by any significant amount. After reaching this conclusion, we focused instead on trying to minimize the energy of the reflected stress wave circulating up and down through the bar over a relatively long period of time (10 ms). Once again, we used numerical simulations with the Sierra/SM module to investigate the effects of various bar-beam connection fixities, alternative wave beam materials, and parameters of an asymmetric monolithic design of the bar-and-beam system. This study demonstrated that various parameters can significantly affect the energy of the wave reflections, with the difference between best and worst configurations being about one order of magnitude in terms of energy. Based on the obtained results, we conclude with concrete takeaways for Dropkinson bar users and propose potential directions for future research and optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding the Breakdown Behavior of Ultrawide-Bandgap Boron Nitride Power Diodes Using Device Modeling

Herein, a device study using technology computer-aided design simulation to theoretically analyze the electrical performance of ultrawide-bandgap boron nitride (BN)-based vertical junction devices is performed, including h-BN Schottky diode, h-BN pn diode, and h-BN/AlN pn diode; this is also the first demonstration of the BN power devices in simulation. The material properties of BN are defined with recently reported data, and the physical mechanisms of the device performance are systematically investigated. Additionally, the h-BN junctions in this simulation shows excellent performance, especially for breakdown behaviors. Schottky diode shows a turn-on voltage of 0.6 V for Pt Schottky contact and breakdown voltages over 450 V for 5 μm, 6 × 10 15 cm –3 p-type-doped drift layer; The h-BN pn diode shows a turn-on voltage of 6 V and breakdown voltages over 3 kV with a critical electric field of 13.6 MV cm –1 for 2.5 μm, 2 × 10 16 cm –3 p-type-doped drift layer. The h-BN/AlN heterojunction pn diode shows a turn-on voltage of 5.8 V and breakdown voltage over 2 kV for 2.5 μm, 2 × 10 16 cm –3 n-type-doped AlN drift layer. Herein, an understanding of the device principles of vertical BN junctions is provided, which can serve as a reference for the future development of robust BN power electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗
DOE OSTI2021

Automated Parameterized CFD Simulations of Phase-Change Material Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.

25 ENERGY STORAGE↗
DOE OSTI2020

Benchmark Solutions for Radiation Transport in Stochastic Media with Inhomogeneous Material Statistics

Accurately solving implicit Monte Carlo (IMC) thermal photon transport problems with mixed material cells is important in realistic applications. The production IMC package at LLNL treats mixed material cells arising from ALE remap and hydrodynamics using the same approximate model. The new Imp IMC thermal photon transport package currently under development has both a material interface reconstruction (MIR) algorithm and a Levermore-Pomraning (LP) stochastic medium algorithm for treating mixed material cells. Existing stochastic medium algorithms for treating mixed material cells in IMC lack a complete theoretical basis. The IMC LP algorithm implementation has been demonstrated to reproduce published deterministic LP solutions for the particular case of spatially homogeneous material statistics. Realistic simulations will include spatially inhomogeneous material statistics (material mean chord lengths). In a previous investigation, the LP-model for transport in binary stochastic media in rod geometry was generalized to accommodate spatially varying material chord lengths, i.e., the mixing statistics were allowed to be nonhomogeneous. Analytical solutions were obtained and used to produce a verifi cation suite for the Imp IMC Levermore-Pomraning implementation for different spatial variations of the chord lengths. However, the accuracy of the LP model when the mixing statistics are nonhomogeneous has not been assessed and leaves open the question of whether local accuracy is improved or further degraded when chord lengths are not uniform. This shortcoming is rectifi ed here by developing benchmark analytic solutions for transport in binary Markovian stochastic mixtures in rod geometry with nonhomogeneous mixing statistics, using spatially varying chord lengths considered in the previous investigation based on the LP model. Methods for sampling a nonhomogeneous Poisson process (NHPP) are first described and used to construct individual realizations of the binary mixtures in rod geometry. Analytic solutions are then obtained for the forward and backward directed fluxes on a given realization, now viewed as a deterministic medium with alternating layers of the two materials with known interface locations. Finally, material averaged scalar fluxes are obtained using these sampling schemes with spatially linear and quadratic chord lengths and used to assess the accuracy of the previously obtained LP-model results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Quantum-classical embedding via ghost Gutzwiller approximation for enhanced simulations of correlated electron systems

Simulating correlated materials on present-day quantum hardware remains challenging due to limited quantum resources. Quantum embedding methods offer a promising route by reducing computational complexity through the mapping of bulk systems onto effective impurity models, allowing more feasible simulations on pre- and early-fault-tolerant quantum devices. Here, this work develops a quantum-classical embedding framework based on the ghost Gutzwiller approximation to enable quantum-enhanced simulations of ground-state properties and spectral functions of correlated electron systems. Circuit complexity is analyzed using an adaptive variational quantum algorithm on a statevector simulator, applied to the infinite-dimensional Hubbard model with increasing ghost mode numbers from 3 to 5, resulting in circuit depths growing from 16 to 104. Noise effects are examined using a realistic error model, revealing significant impact on the spectral weight of the Hubbard bands. To mitigate these effects, the Iceberg quantum error detection code is employed, achieving up to 40% error reduction in simulations. Finally, the accuracy of the density matrix estimation and the derived spectral function is benchmarked on IBM and Quantinuum quantum hardware, featuring distinct qubit-connectivity and employing multiple levels of error mitigation techniques.

Chen, I-Chi [Ames Laboratory (AMES), Ames, IA (Uni↗
DOE OSTI2020

The TRANSCEND University Consortium: Theme 2 - Site Decommissioning, Deconstruction and Remediation - 20329

This paper will discuss the aims, objectives and progress to-date under Theme 2 of the TRANSCEND consortium project: Site decommissioning, deconstruction and remediation. Decommissioning nuclear sites involves waste retrieval, decontamination, deconstruction and, where necessary, containment and/or remediation of the remaining structure and surrounding land. Critical to management of these processes is limiting radiation exposure for the workforce, restricting the spread of radionuclides in groundwater, surface water and airborne particulates, and minimising the volume of contaminated waste for disposal. The aim of Theme 2 research is to develop new technologies for monitoring, remediation and containment that serve to minimise the volume of radioactively contaminated waste for disposal, for application prior to, during and after retrieval, deconstruction and decontamination operations. Prior research, conducted under the previous DISTINCTIVE project, demonstrated that colloidal silica grout can penetrate low permeability materials (including cement) for hydraulic barrier formation, and improved sorption capacity. The silica grout can be injected at surface using extremely small (potentially gravity-driven) fluid pressures, without the need for borehole drilling. Current work is investigating the erodibility of silica-grouted soils for inhibition of airborne and water-borne particulates as well as enhancing the grout's sorption capacity by addition of other materials to provide a chemical, as well as hydraulic, barrier to subsurface migration. EK remediation uses low voltage DC current to control migration of contaminants in porous media as well as to remove or degrade them. Researchers in the consortium have already demonstrated that low-energy ex-situ EK techniques can be used to provide remediation and volume minimisation for AWE legacy wastes in the UK. Under the current project researchers are building EK test cells containing simulated site materials at laboratory and intermediate(m)-scales to: remove, focus or degrade contaminants (remediation or waste minimisation); and direct subsurface water, chemical and colloid flow (fencing/containment or forced migration). This lab-based research is being informed by numerical models of EK processes that can subsequently be used to design full-scale on-site applications by nuclear site holders. Research under theme 2 of TRANSCEND will ultimately combine novel non-invasive detection technologies with EK techniques and colloidal silica grout barriers, to optimise the containment of radionuclide contamination in soils. Thus, allowing us to detect in-situ contamination, mobilise it to a selected location and grout it in-situ, prior to the initiation of decommissioning and deconstruction operations. (authors)

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

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2025

A foundation model for atomistic materials chemistry

Atomistic simulations of matter, especially those that leverage first-principles (ab initio) electronic structure theory, provide a microscopic view of the world, underpinning much of our understanding of chemistry and materials science. Over the last decade or so, machine-learned force fields have transformed atomistic modeling by enabling simulations of ab initio quality over unprecedented time and length scales. However, early machine-learning (ML) force fields have largely been limited by (i) the substantial computational and human effort required to develop and validate potentials for each particular system of interest and (ii) a general lack of transferability from one chemical system to the next. Here, we show that it is possible to create a general-purpose atomistic ML model, trained on a public dataset of moderate size, that is capable of running stable molecular dynamics for a wide range of molecules and materials. We demonstrate the power of the MACE-MP-0 model-and its qualitative and at times quantitative accuracy-on a diverse set of problems in the physical sciences, including properties of solids, liquids, gases, chemical reactions, interfaces, and even the dynamics of a small protein. The model can be applied out of the box as a starting or "foundation" model for any atomistic system of interest and, when desired, can be fine-tuned on just a handful of application-specific data points to reach ab initio accuracy. Establishing that a stable force-field model can cover almost all materials changes atomistic modeling in a fundamental way: experienced users obtain reliable results much faster, and beginners face a lower barrier to entry. Foundation models thus represent a step toward democratizing the revolution in atomic-scale modeling that has been brought about by ML force fields.

Batatia, Ilyes↗
DOE OSTI2024

Energy-dependent implementation of secondary electron emission models in continuum kinetic sheath simulations

The plasma-material interactions present in multiple fusion and propulsion concepts between the flow of plasma through a channel and a material wall drive the emission of secondary electrons. This emission is capable of altering the fundamental structure of the sheath region, significantly changing the expected particle fluxes to the wall. The emission spectrum is separated into two major energy regimes, a peak of elastically backscattered primary electrons at the incoming energy, and cold secondary electrons inelastically emitted directly from the material. The ability of continuum kinetic simulations to accurately represent the secondary electron emission is limited by relevant models being formulated in terms of monoenergetic particle interactions which cannot be applied directly to the discrete distribution function. As a result, rigorous implementation of energy-dependent physics is often neglected in favor of simplified, constant models. We present here a novel implementation of semi-empirical models in the boundary of continuum kinetic simulations which allows the full range of this emission to be accurately captured in physically-relevant regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗