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 577 records · Page 32

Ultra-fast interpretable machine-learning potentials

Abstract All-atom dynamics simulations are an indispensable quantitative tool in physics, chemistry, and materials science, but large systems and long simulation times remain challenging due to the trade-off between computational efficiency and predictive accuracy. To address this challenge, we combine effective two- and three-body potentials in a cubic B-spline basis with regularized linear regression to obtain machine-learning potentials that are physically interpretable, sufficiently accurate for applications, as fast as the fastest traditional empirical potentials, and two to four orders of magnitude faster than state-of-the-art machine-learning potentials. For data from empirical potentials, we demonstrate the exact retrieval of the potential. For data from density functional theory, the predicted energies, forces, and derived properties, including phonon spectra, elastic constants, and melting points, closely match those of the reference method. The introduced potentials might contribute towards accurate all-atom dynamics simulations of large atomistic systems over long-time scales.

36 MATERIALS SCIENCE↗

Quantum Communication Networks for Energy Applications: Review and Perspective

Abstract The energy sector is expected to undergo significant changes in the coming decades with the advent of new technologies, including smart grid development, microgrid expansion, increasing electric vehicle and renewable energy usage, and enhanced measures to minimize greenhouse gas emission, among others. In tandem, these changes are expected to create new opportunities for the deployment of quantum technologies within the energy sector. Building on the authors' previous reviews on the current state of and future opportunities for quantum sensing, quantum computing and quantum simulations for energy sector applications, this work provides an overview of recent progress in quantum networking and communications for the energy industry, with a focus on platforms, devices, and protocols, including quantum teleportation and quantum key distribution. Specific areas of relevance to the energy sector are then analyzed, including the role of quantum networks for greenhouse gas monitoring, secure data collection and transmission in smart grids, nuclear power plants’ safety, facilitating oil and gas exploration, and other energy‐relevant applications. This review concludes with a brief overview of areas for future innovation, including the need for platforms for simulating quantum networks, quantum material and platform design, and computational approaches to accelerate quantum protocol discovery and development.

Paudel, Hari P.↗

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗

Evolution of highly multimodal Rayleigh–Taylor instabilities

Rayleigh–Taylor (RT) instabilities are important fluid instabilities that arise in inertial confinement fusion (ICF) capsule implosions, and many other contexts. Multi-mode coupling is observed in experiments and plays a substantial role in material mix from RT instabilities. In this work, we study the evolution of highly multimodal perturbations (power law distribution) that approximate those found at manufactured material interfaces. We use simulations of over 2000 different perturbations in the LANL code xRAGE to identify distinct phases in the processes of bubble growth and bubble merger which can be visualized in a 2D phase portrait with clear regimes of mode growth and decay. Our results show that the dynamic evolution of the instability strongly depends on the mode of the perturbations and mode interactions. The merger process accelerates bubble growth. A non-Markovian region and a transition of the instability from: (1) initial exponential growth to (2) linear growth and to (3) quadratic growth and asymptotic behavior, are clearly captured in the phase space. We have developed a quantitative model of bubble growth that reproduces the dynamic behavior of ensembles of perturbations. Implications for ICF capsules designed for robustness against instabilities are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimal local truncation error method for 3-D elasticity interface problems

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10 5 -10 9 times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

42 ENGINEERING↗

A positivity-preserving and conservative intersection-distribution-based remapping algorithm for staggered ALE hydrodynamics on arbitrary meshes

In this study, we introduce new intersection-distribution-based remapping tools for indirect staggered arbitrary Lagrangian-Eulerian (ALE) simulations of multi-material shock hydrodynamics on arbitrary meshes. In addition to conserving momentum and total energy, the three-stage remapper proposed in this work preserves non-negativity of the internal energy. At the first stage, we construct slope-limited piecewise-linear reconstructions of all conserved quantities on zones of the source mesh and perform intersection-based remap to obtain bound-preserving zonal quantities on the target mesh. At the second stage, we define bound-preserving nodal quantities of the staggered ALE discretization as convex combinations of corner quantities. The nodal internal energy is corrected in a way which keeps it non-negative, while providing exact conservation of total energy. At the final stage, we distribute the non-negative nodal internal energy to corners, zones and materials using non-negative weights. Proofs of positivity preservation are provided for each stage. This work is a natural extension of our paper [14] in which a similar intersection-distribution-based remapping procedure was employed. The original version used a nodal kinetic energy fix which did not provably ensure positivity preservation for the zonal internal energy after the final distribution stage. The new algorithm cures this potential drawback by using ‘coordinated’ limiters for piecewise-linear reconstructions, remapping the internal energy to nodes and correcting it before redistribution. The effectiveness of the new nodal fix is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗

CALPHAD-based ICME design of single-step aging to enhance mechanical strength of WAAM Haynes 282

To match the strength of wire-arc additive manufactured Haynes 282 to its wrought counterpart via a single-step aging heat treatment, the CALPHAD (Calculation of Phase Diagrams) method is integrated with physics-based process-structure-property models and experimental validation. The integrated computational materials engineering (ICME) framework simulates the effects of aging on γ′ and M 23 C 6 precipitation and the resulting yield strength. To improve simulation reliability, the interfacial energies between γ/γ′ and γ/M 23 C 6 carbides were estimated by comparison with precipitation kinetic modeling and measured precipitate sizes. γ′ and M23C6 were found to precipitate simultaneously between 640 and 860 °C, producing microstructures similar to those produced by two-step aging. The optimal γ′ size for peak yield stress was calculated to be 20–23 nm. WAAM Haynes 282 aged at 780 °C for 50 h exceeded the mechanical performance of its wrought counterpart subjected to two-step aging, though desired properties can also be achieved at 800 °C for 16 h or less. The error in yield strength is less than 20 MPa, demonstrating good agreement between the modeling framework and experiments. Creep studies showed that WAAM Haynes 282 exceeded the calculated rupture time, reaching 481 h. This proposed methodology can accelerate the design of aging heat treatments for any γ′-strengthened nickel-base alloy, minimizing the resources required for trial-and-error experiments.

CALPHAD↗

Design and Application of Materials for Sequestration and Immobilization of 99 Tc

99 Technetium ( 99 Tc) is a hazardous radionuclide generated by the nuclear industry that poses a serious environmental threat. The wide variation and complex chemistries of nuclear waste streams containing 99 Tc often create unique, site specific challenges when sequestering and immobilizing the waste in a matrix suitable for long-term storage and disposal. Therefore, an effective management plan for 99 Tc containing radioactive waste will likely require a variety of suitable materials/matrices capable of adapting to and addressing these challenges. In this review, we discuss and highlight the key developments for effective removal and immobilization of 99 Tc in inorganic waste forms. Specifically, we review the synthesis, characterization, and application of materials for targeted removal of 99 Tc from (simulated) waste solutions under various experimental conditions. These materials include: i) layered double hydroxides (LDHs), ii) metal-organic frameworks (MOFs), iii) ion-exchange resins (IERs) as well as cationic organic polymers (COPs), and iv) surface modified natural clay materials (SMCMs). Secondly, we discuss some of the major and recent developments towards 99 Tc immobilization in i) glass, ii) cement, and iii) iron mineral waste forms. Finally, we present future challenges that need to be addressed for the design, synthesis, and selection of suitable matrices for the efficient sequestration and immobilization of 99 Tc from targeted wastes. Finally, the purpose of this review is to inspire the researches on the design and application of various suitable materials/matrices for selective removal of 99 Tc present globally in different radioactive wastes and its immobilization in stable/durable waste forms.

36 MATERIALS SCIENCE↗

Multiple Sources of Riparian Wetland Suspended Solids during Episodic Rain Events: Influence on Uranium Transport

Suspended solids can be the primary vector for transporting contaminants in streams. The objective of this study was to determine whether changes in the properties of suspended solids during rain events impacted contaminant transport. Stream water was collected during five episodic events downstream from a U-contaminated wetland located in South Carolina, USA. The suspended particles were initially composed of Fe-flocs (particles formed in situ prior to the rain event) that had significantly greater Fe, Mn, organic-C, and U content than particles collected later during a sampling rain event. XANES and EXAFS revealed that U in the Fe-flocs was U(VI) and that it was not incorporated in a mineral structure but existed as inner- or outer-sphere adsorbed uranyl species associated with organic matter and Fe-oxides. The uranyl had an extraordinarily high affinity for the suspended solids, with solid to liquid U ratios of >72,000 (μg/kg)/(μg/L). After the initial flush of Fe-flocs, a greater fraction of the suspended solids had lower organic-C, Fe, Mn, and amorphous phases and were composed of more quartz, kaolinite, and gibbsite, resulting in lower U concentrations than those in the solids collected earlier in the rain event. This study highlights the importance of understanding suspended solids as transport vectors and their potential dynamic nature during rain events.

computer simulations↗

Rapid Adaptation of Chemical Named Entity Recognition Using Few-Shot Learning and LLM Distillation

Named entity recognition (NER) has been widely used in chemical text mining for the automatic identification and extraction of chemical entities. However, existing chemical NER systems primarily focus on scenarios with abundant training data, requiring significant human effort on annotations. This poses challenges for applications in the chemical field, such as catalysis, where many advancements have traditionally relied on trial-and-error investigations and incremental adjustment of variables. This hinders catalysis science and technology progress in addressing emerging energy and environmental crises. In this work, we propose a few-shot NER model that can quickly adapt to extract new types of chemical entities by using only a limited number of annotated examples. Our model employs a metric-learning approach to transfer entity similarity knowledge from high-resource chemical domains (with abundant annotations) to enable effective entity recognition in low-resource specialized domains (limited annotation). We validate the effectiveness of our model on a few-shot chemical NER benchmark built based on six existing chemical NER data sets. Experiments show that the proposed few-shot NER model can achieve reasonable performance with only 5 examples per entity type and shows consistent improvement as the number of examples increases. Furthermore, we demonstrate how the proposed model can be trained with large language model (LLM) annotated data, opening a new pathway for rapid adaptation of NER systems. Furthermore, our approach leverages the knowledge broadness of large language models for chemistry while distilling this knowledge into a lightweight model suitable for efficient and in-house use.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent-Treated Zirconium-Based Nanoporous UiO-66 Metal–Organic Frameworks for Enhanced CO 2 Capture

Metal–organic framework (MOF)-derived nanoscale porous materials are widely deployed in carbon capture, but the CO 2 capacity is varied by different post-treatments, and its mechanism is still unclear. In this work, we prepared UiO-66 and treated it with methanol solvent and thermal activation approaches, which showed ~3 times enhanced CO 2 capacities from 15.1 to 45.0 mg/g and excellent recyclability of 10 cycles. The methanol treatment efficiently removes the residual guest molecules, including N,N-dimethylformamide, dangling organic linkers, and their derivatives, in the micropores (~0.8 nm) of UiO-66 and improves the surface area, pore volume, and void fraction to enhance the CO 2 capacity close to the ideal UiO-66 materials. The molecular dynamics simulation also proved a good linear relationship between the surface areas, void fraction, and CO 2 capacity. This work provides a deep understanding of the MOF’s activation mechanism and its applications in CO 2 capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Generating multi-scale Li-ion battery cathode particles with radial grain architectures using stereological generative adversarial networks

Abstract Understanding structure-property relationships of Li-ion battery cathodes is crucial for optimizing rate-performance and cycle-life resilience. However, correlating the morphology of cathode particles, such as in LiNi0.8Mn0.1Co0.1O2 (NMC811), and their inner grain architecture with electrode performance is challenging, particularly, due to the significant length-scale difference between grain and particle sizes. Experimentally, it is not feasible to image such a high number of particles with full granular detail. A second challenge is that sufficiently high-resolution 3D imaging techniques remain expensive and are sparsely available at research institutions. Here, we present a stereological generative adversarial network-based model fitting approach to tackle this, that generates representative 3D information from 2D data, enabling characterization of materials in 3D using cost-effective 2D data. Once calibrated, this multi-scale model can rapidly generate virtual cathode particles that are statistically similar to experimental data, and thus is suitable for virtual characterization and materials testing through numerical simulations. A large dataset of simulated particles with inner grain architecture has been made publicly available.

25 ENERGY STORAGE↗

Oscillation-like diffusion of two-dimensional liquid dusty plasmas on one-dimensional periodic substrates with varied widths

The long-time diffusion of two-dimensional dusty plasmas on a one-dimensional periodic substrate with varying widths is investigated using Langevin dynamical simulations. When the substrate is narrow and the dust particles form a single row, the diffusion is the smallest in both directions. In this study, we find that as the substrate width gradually increases to twice its initial value, the long-time diffusion of the two-dimensional dusty plasmas first increases, then decreases, and finally increases again, giving an oscillation-like diffusion with varying substrate widths. When the width increases to a specific value, the dust particles within each potential well arrange themselves in a stable zigzag pattern, greatly reducing the diffusion and leading to the observed oscillation in the diffusion with the increase in the width. In addition, the long-time oscillation-like diffusion is consistent with the number of dust particles that are hopping across the potential wells of the substrate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

X-ray continuum spectroscopy of inertial confinement fusion implosions at the National Ignition Facility

A methodology for measuring x-ray continuum spectra of inertial confinement fusion (ICF) implosions is described. The method relies on the use of ConSpec, a high-throughput spectrometer using a highly annealed pyrolytic graphite crystal , which measures the spectra in the ≃20–30 keV range. Due to its conical shape, the crystal is sagittally focusing a Bragg-reflected x-ray spectrum into a line, which enhances the recorded x-ray emission signal above the high neutron-induced background accompanying ICF implosions at the National Ignition Facility. To improve the overall measurement accuracy, the sensitivity of the spectrometer measured in an off-line x-ray laboratory setting was revised. The error analysis was expanded to include the accuracy of the off-line measurements, the effect of the neutron-induced background, as well as the influence of possible errors in alignment of the instrument to the ICF target. Here we demonstrate how the improved methodology is applied in the analysis of ConSpec data with examples of a relatively low-neutron-yield implosion using a tritium–hydrogen–deuterium mix as a fuel and a high-yield deuterium–tritium (DT) implosion producing high level of the background. In both cases, the shape of the measured spectrum agrees with the exponentially decaying spectral shape of bremsstrahlung emission to within ±10%. In the case of the high-yield DT experiment, non-monotonic deviations slightly exceeding the measurement uncertainties are observed and discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Paradigms of frustration in superionic solid electrolytes

Superionic solid electrolytes have widespread use in energy devices, but the fundamental motivations for fast ion conduction are often elusive. In this Perspective, we draw upon atomistic simulations of a wide range of superionic conductors to illustrate some ways frustration can lower diffusion cation barriers in solids. Based on our studies of halides, oxides, sulfides and hydroborates and a survey of published reports, we classify three types of frustration that create competition between different local atomic preferences, thereby flattening the diffusive energy landscape. These include chemical frustration, which derives from competing factors in the anion–cation interaction; structural frustration, which arises from lattice arrangements that induce site distortion or prevent cation ordering; and dynamical frustration, which is associated with temporary fluctuations in the energy landscape due to anion reorientation or cation reconfiguration. For each class of frustration, we provide detailed simulation analyses of various materials to show how ion mobility is facilitated, resulting in stabilizing factors that are both entropic and enthalpic in origin. We propose the use of these categories as a general construct for classifying frustration in superionic conductors and discuss implications for future development of suitable descriptors and improvement strategies.

25 ENERGY STORAGE↗

Implementing a microphysics model in hydrodynamic simulations to study the initial plasma formation in dielectric ablator materials for direct-drive implosions

A microphysics model to describe the photoionization and impact ionization processes in dielectric ablator materials like plastic has been implemented into the 1-D hydrodynamic code LILAC for planar and spherical targets. Currently, the initial plasma formation during the early stages of a laser drive are modeled in an ad hoc manner, until the formation of a critical surface. Implementation of the physics-based models predict higher values of electron density, electron temperature, and pressure than the ad hoc model. Moreover, the numerical predictions are consistent with previous experimental observations of the shinethrough mechanism in plastic ablators. For planar targets, a decompression of the rear end of the target was observed, that is similar to recent experiments. An application of this model is to understand the laser-imprint mechanism that is caused by nonuniform laser irradiation due to single beam speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗