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

Spectral Mesh-Free Quadrature for Planar Regions Bounded by Rational Parametric Curves

This article presents spectral, mesh-free, Green’s theorem-based numerical quadrature schemes for integrating functions over planar regions bounded by rational parametric curves. Our algorithm proceeds in two steps: (1) We first find intermediate quadrature rules for line integrals along the region’s boundary curves corresponding to Green’s theorem. (2) We then use a high-order quadrature rule to compute the numerical antiderivative of the integrand along a coordinate axis, which is used to evaluate the Green’s theorem line integral. We present two methods to compute the intermediate quadrature rule. The first is spectrally accurate (it converges faster than any algebraic order with respect to number of quadrature points) and is relatively easy to implement, but has no guarantee of polynomial exactness. The second guarantees exactness for polynomial integrands up to a pre-specified degree with an a priori-known number of quadrature points and retains the convergence properties of the first, but is slightly more complicated. The quadrature schemes have applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order meshes, and initialization of multi-material simulations with rational geometry. We compare the quadrature schemes produced using our method to other methods in the literature and show that they are much more efficient both in terms of number of quadrature points and computational time. We provide an open-source implementation of the algorithm in MATLAB.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. Additionally, the scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys

We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials’ design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum–niobium–titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.

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

Fused x-ray and fast neutron CT reconstruction for imaging large and dense objects

Megavolt x-ray computed tomography (CT) is a powerful tool for three-dimensional characterization. However, its utility is limited for large objects composed of high-atomic number (Z) materials, where x rays fail to penetrate. Information from fast neutron CT (FNCT) can complement x-ray CT reconstructions since fast neutrons can more readily penetrate high-Z objects. In this work, we demonstrate a method for combining FNCT and x-ray CT data to create a single reconstruction, more accurate than could be achieved with either x rays or fast neutrons alone. The algorithm was tested on an exemplar comprising multiple concentric, nested cylinders of different materials. Simulated and empirical x-ray CT data were acquired for the exemplar using a 9 MV bremsstrahlung spectrum. Additional simulated and empirical FNCT data were acquired using an accelerator based fast neutron source. The FNCT data were used to synthesize x-ray CT data and augment the x-ray CT data missing due to lack of penetration. This approach mitigates artifacts that would otherwise negatively affect the accuracy and resolution of a single-modality reconstructed volume.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Probing the state of hydrogen in δ − AlOOH at mantle conditions with machine learning potential

Hydrous and nominally anhydrous minerals are a fundamental class of solids of enormous significance to geophysics. They are the water carriers in the deep geological water cycle and impact structural, elastic, plastic, and thermodynamic properties and phase relations in Earth's forming aggregates (rocks). They play a critical role in the geochemical and geophysical processes that shape the planet. Their complexity has prevented predictive calculations of their properties, but progress in materials simulations ushered by machine-learning potentials is transforming this state of affairs. Here, we adopt a hybrid approach that combines deep learning potentials (DPs) with the strongly constrained and appropriately normed meta-generalized gradient approximation functional to simulate a prototypical hydrous system. We illustrate the success of this approach to simulate δ − AlOOH ( δ ) , a phase capable of transporting water down to near the core-mantle boundary of the Earth ( ∼ 2900 km depth and ∼ 135 GPa ) in subducting slabs. A high-throughput sampling of phase space using molecular dynamics simulations with DPs sheds light on the hydrogen-bond behavior and proton diffusion at geophysical conditions. These simulations provide a pathway for a deeper understanding of these crucial components that shape Earth's internal state. Published by the American Physical Society 2024

Luo, Chenxing (ORCID:0000000341166851)↗
DOE OSTI2020

Cross-Cutting Software Solutions in Support of Experimental Analysis Challenges at National Scattering Facilities

Rapid improvements in instrumentation at neutron and light source combined with new capabilities in computing, data analytics, and machine learning offer unprecedented scientific capabilities to address fundamental and applied science questions in a wide area of subjects from quantum materials to biological systems. In order to make these capabilities accessible to the broader scientific community, significant developments are needed related to workflow spanning from experimental control to advanced materials simulations as well as the underpinning scientific software. In this chapter, we explore these needs and recent developments at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory. A central part of the workflow at the SNS is the data reduction framework Mantid which is developed as a partnership between European and US neutron-scattering facilities.

Proffen, Thomas↗
DOE OSTI2021

Axom

Axom is an open-source library of "building block" software components that provide core infrastructure capabilities for HPC applications. Such capabilities include: input file parsing and verification for various file formats, format specification for material shape input to multi-material simulations and tools for placing those shapes on meshes, parallel distributed coordination of diagnostic messages, geometric primitives, spatial queries and spatial search acceleration data structures, in memory key-value data store for managing simulation data and parallel file I/O. Axom components are used widely across LLNL Advanced Simulation and Computing program applications. It is also part of the LLNL institutionally-supported RADIUSS project, which promotes and funds its adoption by projects across LLNL.

Zagaris, George↗
DOE OSTI2022

cyqc

A collection of codes which implement quantum algorithms for correlated materials simulations.

Yao, Yongxin↗
DOE OSTI2023

mkite

mkite is a distributed computing platform for materials simulation. mkite is built with the server-client pattern, decoupling production databases from client runners. When used in combination with message brokers, mkite enables any available client to perform calculations without prior hardware specification on the server side. Furthermore, the software enables the creation of complex workflows with multiple inputs and branches, facilitating the exploration of combinatorial chemical spaces. The mkite suite provides recipes and tools to interact with package such as VASP, but is extensible to any other simulation package. Finally, mkite helps keeping the provenance of calculations in a SQL database. A complete description of the software is available at https://arxiv.org/abs/2301.08841.

Schwalbe Koda, Daniel↗
DOE OSTI2023

tih_vqe [SWR-23-32]

This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

Graf, Peter↗
DOE OSTI2022

Compressive strength of the plastic-bonded explosive surrogate idoxuridine under quasistatic and dynamic compression

Inert simulant materials, or “mocks”, are often used as surrogates for plastic-bonded explosives (PBX) in non-detonative tests in order to mitigate hazards. Mocks should reproduce as many non-detonative properties of the explosive as possible, including structural behavior in a variety of thermal and mechanical environments. Recently, the molecular crystal idoxuridine (IDOX) has been identified as an ideal mock for the main component in the explosive polymer-matrix composite PBX 9501, and has performed favorably under quasistatic loading conditions. Here, in order to assess robustness over a range of mechanical environments, plastic-bonded IDOX was compression tested from 0.001/s to 2000/s strain rates and compared to PBX 9501 historical data. Plastic-bonded IDOX showed good agreement to PBX 9501 across these strain rates, justifying continued development and production as a mock.

Burch, Alexandra C.↗
DOE OSTI2021

Embedding Sensors in 3D Printed Metal Structures

The Transformational Challenge Reactor (TCR) program is leveraging recent advances in modeling and simulation, materials, and additive manufacturing (AM) technologies to design a modern nuclear reactor. Some of the main TCR technologies include in situ monitoring and the integration of sensors during the manufacturing of quality-significant nuclear reactor components. This report describes the general procedure and process optimization for embedding sensors within generic stainless steel 316 (SS316) components using laser powder bed fusion (LPBF). A more detailed, quality-significant test plan and supporting procedures are available upon request (ORNL/TM-2021/2127). LPBF involves the use of a scanning laser to selectively melt regions of a powder bed, additively building a part layer by layer. This report describes the LPBF processing technique and discusses the effects of LPBF processing parameters on the success of the sensor embedding process. Experiments used machined cavities in the form of channels in an SS316 base for the sensors to lay in while material is additively built over the top, thereby embedding them in an SS316 matrix. A preliminary investigation involved using empty SS316 sheaths as surrogates to explore the effects of various LPBF processing parameters and the dimensional requirements of the machined channels. Microstructural investigations showed that a smaller channel width/depth combination closer to the sensor’s diameter was best for the embedding process. After the desired parameters were selected, Type-K thermocouples were embedded and evaluated post-embedding using nondestructive thermal testing, as well as destructive sectioning and microscopy. Post-build characterization showed that the thermocouples were well-bonded to the SS316 matrix and were fully functional after embedding. During thermal testing to temperatures up to 500 °C, the embedded thermocouples read consistently with one another and deviated only slightly from the readings of a nonembedded thermocouple located within the furnace. This slight discrepancy was most likely due to differences in the thermal time constants for a nonembedded thermocouple vs. a thermocouple embedded in a solid SS316 block. The results presented in this report will serve as the foundation for future work that will focus on embedding sensors in relevant TCR reactor components and eventually testing those components under neutron irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unpacking Modeling Analysis for the Circular Economy

This session will cover: (1) How models simulate material recovery, recycling, product life cycles and waste management. (2) Using life cycle analysis (LCA) to measure environmental impacts across product life cycles and guide decision-making. (3) Optimizing resource use, reducing energy consumption and minimizing waste in recycling systems. (4) Evaluating EPR and other interventions through scenario analysis to make informed decisions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Hydrodynamic Coupling to a Homogenized Radiation Transport Method based on Young Measures

Resolving radiation transport fields subject to opacity profiles with strong, oscillatory line structure while potentially falling under intermediate optical depth conditions presents a numerical challenge in radiation transport modeling. The Young measure-based homogenization technique formulated by Haut et al. (2017) was investigated as a candidate method for resolving radiation fields under such conditions more accurately. The method was compared against frequently-utilized mean opacity methods as the Rosseland and Planck formulations. In this work, all methods were tested through radiation slab calculations separately comprised of aluminum, copper, and krypton, each for different thermodynamic conditions. Following these offline radiation slab calculations, demonstrations shifted towards the SCEPTRE radiation transport code and, subsequently, the multiphysics ALEGRA code for approximately-coupled radiation-material simulations. Throughout all the simulations shown in this study, for a fixed computational cost, the homogenized method was observed to be more accurate than any of the solutions determined through traditional mean opacity approaches.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

FY25 Task 5: Small-Scale Mixing

The U.S. Department of Energy (DOE) Hanford Site has 177 underground storage tanks that contain a complex and diverse mix of chemical and radioactive wastes from past nuclear fuel reprocessing and waste management operations. The strategy of the DOE Hanford Field Office is to retrieve this waste, ~20 vol% of which is in the form of insoluble undissolved solids (UDS) or sludge, and treat it via immobilization at the Hanford Waste Treatment and Immobilization Plant (WTP). The diverse properties and characteristics of Hanford tank waste lead to major challenges related to its transport from the underground tanks to the WTP. These challenges, however, can be addressed by investigating the behavior of tank waste samples and simulant materials and evaluating their behavior against the capabilities of mixing and transport system designs that may be incorporated to retrieve and treat the waste.

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

An Alternative Conceptual Model for the Spent Nuclear Fuel–Water Interaction in Deep Geologic Disposal Conditions

For the long-term safety assessment of direct disposal of spent nuclear fuel in deep geologic repositories, knowledge on the radionuclide release rate from the UO 2 matrix is essential. This work provides a conceptual model to explain the results of leaching experiments involving used nuclear fuel or simulant materials in confirmed reducing conditions. Key elements of this model are: direct effect of radiation from radiolytic species (including defects and excited states) in the solid and in the first water layers in contact with its surface; and excess H 2 may be produced due to processes occurring at the surface of the spent fuel and in confined water volumes, which may also play a role in keeping the spent fuel surface in a reduced state. The implication is that the fractional radionuclide release rate used in most long-term safety assessments (10 -7 year -1 ) is over estimated because it assumes that there is net UO 2 oxidation caused by radiolysis, in contrast with the alternative conceptual model presented here. Furthermore, conventional water radiolysis models and radiation chemical yields published in the literature are not directly applicable to a heterogeneous system such as the spent fuel–water interface. Suggestions are provided for future work to develop more reliable models for the long-term safety assessment of spent nuclear fuel disposal.

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

Sorption and Desorption Study of a Nerve-Agent Simulant from Office Materials for Forensic Applications

Dimethyl methylphosphonate (DMMP) was used as a simulant to investigate the sorption and desorption of G-series nerve agents such as sarin (GB) released into an office space. DMMP was disseminated systematically as an aerosol and vapor on three common office materials, cardboard, polyurethane foam, and cured acrylic-based paint. The amount of DMMP initially captured on the office materials and the amount retained over a period of 10 hours was tracked gravimetrically and by solvent extraction and gas chromatography/mass spectrometry (GC/MS). Physicochemical properties of the office media including the specific surface area per mass, polarity, and Henry’s constants were determined to help explain the capture and release of DMMP that was observed. Polyurethane foam was found to have the largest sorption capacity for DMMP in a low-humidity setting, however, cardboard was found to have the largest sorption capacity at higher humidity. A fraction of the DMMP collected desorbed with first-order kinetics from all of the office materials in an open-air atmosphere over 8 hours. However, importantly for forensic purposes relevant to an indoor attack, a significant quantity (most notably on cardboard) of the nerve agent simulant was retained after 6 weeks and could be analyzed efficiently using solvent extraction followed by GC/MS.

chemical forensics↗
DOE OSTI2025

Programmable simulations of molecules and materials with reconfigurable quantum processors

Simulations of quantum chemistry and quantum materials are believed to be among the most important applications of quantum information processors. However, realizing practical quantum advantage for such problems is challenging because of the prohibitive computational cost of programming typical problems into quantum hardware. Here we introduce a simulation framework for strongly correlated quantum systems represented by model spin Hamiltonians that uses reconfigurable qubit architectures to simulate real-time dynamics in a programmable way. Our approach also introduces an algorithm for extracting chemically relevant spectral properties via classical co-processing of quantum measurement results. We develop a digital–analogue simulation toolbox for efficient Hamiltonian time evolution using digital Floquet engineering and hardware-optimized multi-qubit operations to accurately realize complex spin–spin interactions. As an example, we propose an implementation based on Rydberg atom arrays. In addition, we show how detailed spectral information can be extracted from the dynamics through snapshot measurements and single-ancilla control, enabling the evaluation of excitation energies and finite-temperature susceptibilities from a single dataset. To illustrate the approach, we show how to use the method to compute key properties of a polynuclear transition-metal catalyst and two-dimensional magnetic materials.

74 ATOMIC AND MOLECULAR PHYSICS↗