Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing methods”

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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Enabling deeper learning on big data for materials informatics applications

Abstract The application of machine learning (ML) techniques in materials science has attracted significant attention in recent years, due to their impressive ability to efficiently extract data-driven linkages from various input materials representations to their output properties. While the application of traditional ML techniques has become quite ubiquitous, there have been limited applications of more advanced deep learning (DL) techniques, primarily because big materials datasets are relatively rare. Given the demonstrated potential and advantages of DL and the increasing availability of big materials datasets, it is attractive to go for deeper neural networks in a bid to boost model performance, but in reality, it leads to performance degradation due to the vanishing gradient problem. In this paper, we address the question of how to enable deeper learning for cases where big materials data is available. Here, we present a general deep learning framework based on Individual Residual learning (IRNet) composed of very deep neural networks that can work with any vector-based materials representation as input to build accurate property prediction models. We find that the proposed IRNet models can not only successfully alleviate the vanishing gradient problem and enable deeper learning, but also lead to significantly (up to 47%) better model accuracy as compared to plain deep neural networks and traditional ML techniques for a given input materials representation in the presence of big data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A machine learning aided interpretable model for rupture strength prediction in Fe-based martensitic and austenitic alloys

Abstract The class of 9–12% Cr ferritic-martensitic alloys (FMA) and austenitic stainless steels have received considerable attention due to their numerous applications in high temperature power generation industries. To design high strength steels with prolonged service life requires a thorough understanding of the long-term properties, e.g., creep rupture strength, rupture life, etc., as a function of the chemical composition and processing parameters that govern the microstructural characteristics. In this article, the creep rupture strength of both 9–12% Cr FMA and austenitic stainless steel has been parameterized using curated experimental datasets with a gradient boosting machine. The trained model has been cross validated against unseen test data and achieved high predictive performance in terms of correlation coefficient ( $$R^{2} > 0.98 $$ R 2 > 0.98 for 9–12% Cr FMA and $$R^{2} > 0.95 $$ R 2 > 0.95 for austenitic stainless steel) thus bypassing the need for additional comprehensive tensile test campaigns or physical theoretical calculations. Furthermore, the feature importance has been computed using the Shapley value analysis to understand the complex interplay of different features.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Machine learning with persistent homology and chemical word embeddings improves prediction accuracy and interpretability in metal-organic frameworks

Machine learning has emerged as a powerful approach in materials discovery. Its major challenge is selecting features that create interpretable representations of materials, useful across multiple prediction tasks. We introduce an end-to-end machine learning model that automatically generates descriptors that capture a complex representation of a material’s structure and chemistry. This approach builds on computational topology techniques (namely, persistent homology) and word embeddings from natural language processing. It automatically encapsulates geometric and chemical information directly from the material system. We demonstrate our approach on multiple nanoporous metal–organic framework datasets by predicting methane and carbon dioxide adsorption across different conditions. Our results show considerable improvement in both accuracy and transferability across targets compared to models constructed from the commonly-used, manually-curated features, consistently achieving an average 25–30% decrease in root-mean-squared-deviation and an average increase of 40–50% in R 2 scores. A key advantage of our approach is interpretability: Our model identifies the pores that correlate best to adsorption at different pressures, which contributes to understanding atomic-level structure–property relationships for materials design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

New proposed ITER divertor design using carbon insert on tungsten to mitigate ELMs and secondary radiation effects on nearby components

Abstract Building a successful device for the magnetic fusion energy production is a great challenge. ITER is an international project of the tokamak based magnetic fusion design being developed for the demonstration of the feasibility of thermonuclear technologies for future realization of successful commercial fusion energy. A key obstacle to a successful magnetic fusion energy production is however, the performance during abnormal events including plasma disruptions and edge-localized modes (ELMs). A credible reactor design must tolerate at least a few of these transient events without serious consequences such as melting of the structure. This paper investigates and compares the performance of the current ITER tokamak design during two types of transient events, i.e., ELMs occurring at normal operation and disruptions during abnormal operation. We simulated the divertor components response using our integrated 3D HEIGHTS package. The simulations include self-consistent modeling of the interaction of the released core plasma particles with the initial solid divertor material, energy deposition processes, vaporization of divertor material, secondary plasma formation and MHD evolution, incident core particles collisions and scattering from this dense secondary plasma, photon radiation of secondary plasma, and the resulting heat loads on nearby components. Our simulations showed that using a small carbon insert around the strike point can significantly reduce the overall expected damage on the tungsten dome structure, reflector plates, and prevent tungsten vaporization and its potential core plasma contamination.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Improving deep learning model performance under parametric constraints for materials informatics applications

Abstract Modern machine learning (ML) and deep learning (DL) techniques using high-dimensional data representations have helped accelerate the materials discovery process by efficiently detecting hidden patterns in existing datasets and linking input representations to output properties for a better understanding of the scientific phenomenon. While a deep neural network comprised of fully connected layers has been widely used for materials property prediction, simply creating a deeper model with a large number of layers often faces with vanishing gradient problem, causing a degradation in the performance, thereby limiting usage. In this paper, we study and propose architectural principles to address the question of improving the performance of model training and inference under fixed parametric constraints. Here, we present a general deep-learning framework based on branched residual learning (BRNet) with fully connected layers that can work with any numerical vector-based representation as input to build accurate models to predict materials properties. We perform model training for materials properties using numerical vectors representing different composition-based attributes of the respective materials and compare the performance of the proposed models against traditional ML and existing DL architectures. We find that the proposed models are significantly more accurate than the ML/DL models for all data sizes by using different composition-based attributes as input. Further, branched learning requires fewer parameters and results in faster model training due to better convergence during the training phase than existing neural networks, thereby efficiently building accurate models for predicting materials properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Roadmap for transforming heterogeneous catalysis with artificial intelligence

Artificial intelligence (AI) is poised to transform heterogeneous catalysis, opening avenues for catalytic materials discovery. By uncovering intricate patterns in high-dimensional data, AI has been reshaping our pursuit of sustainable catalytic processes across the energy, environmental and chemical sectors. This promise, however, hinges on overcoming fundamental barriers, including limitations in data availability and quality, challenges in the generalizability and interpretability of data-augmented decisions, and the persistent gap between in silico predictions and experiments. Furthermore, we outline a forward-looking roadmap for deeply integrating AI into heterogeneous catalysis with an AI-ready data ecosystem, multimodal foundation models, and ultimately autonomous laboratories to accelerate the development of next-generation catalytic technologies via AI-empowered human–machine collaboration.

Computational methods↗
DOE OSTI2024

A generative artificial intelligence framework based on a molecular diffusion model for the design of metal-organic frameworks for carbon capture

Metal-organic frameworks (MOFs) exhibit great promise for CO 2 capture. However, finding the best performing materials poses computational and experimental grand challenges in view of the vast chemical space of potential building blocks. Here, we introduce GHP-MOFassemble, a generative artificial intelligence (AI), high performance framework for the rational and accelerated design of MOFs with high CO 2 adsorption capacity and synthesizable linkers. GHP-MOFassemble generates novel linkers, assembled with one of three pre-selected metal nodes (Cu paddlewheel, Zn paddlewheel, Zn tetramer) into MOFs in a primitive cubic topology. GHP-MOFassemble screens and validates AI-generated MOFs for uniqueness, synthesizability, structural validity, uses molecular dynamics simulations to study their stability and chemical consistency, and crystal graph neural networks and Grand Canonical Monte Carlo simulations to quantify their CO 2 adsorption capacities. We present the top six AI-generated MOFs with CO 2 capacities greater than 2m mol g -1 , i.e., higher than 96.9% of structures in the hypothetical MOF dataset.

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

ChemGraph as an agentic framework for computational chemistry workflows

Atomistic simulations are essential in chemistry and materials science but remain challenging to run due to the expert knowledge required for the setup, execution, and validation stages of these calculations. We present ChemGraph, an agentic framework powered by artificial intelligence and state-of-the-art simulation tools to streamline and automate computational chemistry and materials science workflows. ChemGraph leverages graph neural network-based foundation models for accurate yet computationally efficient calculations and large language models (LLMs) for natural language understanding, task planning, and scientific reasoning to provide an intuitive and interactive interface. We evaluate ChemGraph across 13 benchmark tasks and demonstrate that smaller LLMs (GPT-4o-mini, Claude-3.5-haiku, Qwen-2.5-14B) perform well on simple workflows, while more complex tasks benefit from using larger models. Importantly, we show that decomposing complex tasks into smaller subtasks through a multi-agent framework enables GPT-4o to reach perfect accuracy and smaller LLMs to match or exceed single-agent GPT-4o's performance in these benchmarks.

Computational chemistry↗
DOE OSTI2023

Ab-initio simulations of coherent phonon-induced pumping of carriers in zirconium pentatelluride

Abstract Laser-driven coherent phonons can act as modulated strain fields and modify the adiabatic ground state topology of quantum materials. Here we use time-dependent first-principles and effective model calculations to simulate the effect of the coherent phonon induced by strong terahertz electric field on electronic carriers in the topological insulator ZrTe 5 . We show that a coherent A 1g Raman mode modulation can effectively pump carriers across the band gap, even though the phonon energy is about an order of magnitude smaller than the equilibrium band gap. We reveal the microscopic mechanism of this effect which occurs via Landau–Zener-Stückelberg tunneling of Bloch electrons in a narrow region in the Brillouin zone center where the transient energy gap closes when the system switches from strong to weak topological insulator. The quantum dynamics simulation results are in excellent agreement with recent pump-probe experiments in ZrTe 5 at low temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Bioplastic design using multitask deep neural networks

Abstract Non-degradable plastic waste jeopardizes our environment, yet our modern lifestyle and current technologies are impossible to sustain without plastics. Bio-synthesized and biodegradable alternatives such as polyhydroxyalkanoates (PHAs) have the potential to replace large portions of the world’s plastic supply with cradle-to-cradle materials, but their chemical complexity and diversity limit traditional resource-intensive experimentation. Here, we develop multitask deep neural network property predictors using available experimental data for a diverse set of nearly 23,000 homo- and copolymer chemistries. Using the predictors, we identify 14 PHA-based bioplastics from a search space of almost 1.4 million candidates which could serve as potential replacements for seven petroleum-based commodity plastics that account for 75% of the world’s yearly plastic production. We also discuss possible synthesis routes for the identified promising materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modeling the electrical double layer at solid-state electrochemical interfaces

Models of the electrical double layer (EDL) at electrode/liquid-electrolyte interfaces no longer hold for all-solid-state electrochemistry. Here we show a more general model for the EDL at a solid-state electrochemical interface based on the Poisson–Fermi–Dirac equation. By combining this model with density functional theory predictions, the interconnected electronic and ionic degrees of freedom in all-solid-state batteries, including the electronic band bending and defect concentration variation in the space-charge layer, are captured self-consistently. Along with a general mathematical solution, the EDL structure is presented in various materials that are thermodynamically stable in contact with a lithium metal anode: the solid electrolyte Li 7 La 3 Zr 2 O 12 (LLZO) and the solid interlayer materials LiF, Li 2 O and Li 2 CO 3 . Here, the model further allows design of the optimum interlayer thicknesses to minimize the electrostatic barrier for lithium ion transport at relevant solid-state battery interfaces.

25 ENERGY STORAGE↗
DOE OSTI2025

Modeling the impact of structure and coverage on the reactivity of realistic heterogeneous catalysts

Adsorbates often cover the surfaces of catalysts densely as they carry out reactions, dynamically altering their structure and reactivity. Understanding adsorbate-induced phenomena and harnessing them in our broader quest for improved catalysts is a substantial challenge that is only beginning to be addressed. Here, in this work, we chart a path toward a deeper understanding of such phenomena by focusing on emerging in silico modeling methodologies, which will increasingly incorporate machine learning techniques. We first examine how adsorption on catalyst surfaces can lead to local and even global structural changes spanning entire nanoparticles, and how this affects their reactivity. We then evaluate current efforts and the remaining challenges in developing robust and predictive simulations for modeling such behavior. Last, we provide our perspectives in four critical areas—integration of artificial intelligence, building robust catalysis informatics infrastructure, synergism with experimental characterization, and adaptive modeling frameworks—that we believe can help surmount the remaining challenges in rationally designing catalysts in light of these complex phenomena.

catalytic mechanisms↗
DOE OSTI2023

Gaussian approximation potentials for accurate thermal properties of two-dimensional materials

Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. Here, in this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

2D Materials↗
DOE OSTI2022

Fossil group origins

Aims: We analyse the large-scale structure out to 100 Mpc around a sample of 16 confirmed fossil systems using spectroscopic information from the Sloan Digital Sky Survey Data Release 16. Methods: We computed the distance between our fossil groups (FGs) and the centres of filaments and nodes from the literature. We also studied the density of bright galaxies, since this parameter is thought to be a good mass tracers, as well as the projected over-densities of galaxies. Finally, we applied a friends-of-friends (FoF) algorithm to detect virialised structures around our FGs and obtain an estimate of the mass available in their surroundings. Results: We find that FGs are mainly located close to filaments, with a mean distance of 3.7 ± 1.1 R 200 and a minimum distance of 0.05 R 200 . On the other hand, none of our FGs were found close to intersections, with a mean and minimum distance of 19.3 ± 3.6 and 6.1 R 200 , respectively. There is a correlation that indicates FGs at higher redshifts are found in denser regions, when we use bright galaxies as tracers of the mass. At the same time, FGs with the largest magnitude gaps (Δm 12 > 2.5) are found in less dense environments and tend to host (on average) smaller central galaxies. Conclusions: Our results suggest that FGs formed in a peculiar position within the cosmic web, close to filaments and far from nodes, whereby their interaction with the cosmic web itself may be limited. We deduce that FGs with brightest central galaxies (BCGs) that are relatively faint, high values of Δm 12 , and low redshifts could, in fact, be systems that are at the very last stage of their evolution. Moreover, we confirm theoretical predictions that systems with the largest magnitude gap are not massive.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

MAPS: the MFEM Anisotropic Plasma Solver

Simulating magnetically confined fusion plasmas presents a uniquely challenging problem due to the nonlinear anisotropic heat conduction. We introduce the MAPS (MFEM Anisotropic Plasma Solver) tool, which uses a high-order finite element method to compute transport solutions on unstructured meshes. We show results for a set of three 2-D verification tests, two of which demonstrate the expected convergence properties for various mesh resolutions and polynomial degrees. We then discuss the convergence rate for the third test.

Barnett, Rhea [ORNL] (ORCID:0000000317527979)↗
DOE OSTI2020

Modeling DC electrical breakdown using a truncated emission spectrum for trapped radiation

Spontaneously emitted radiation from excited atoms can be of principle importance in certain modes of electrical breakdown, especially positive streamers and some regimes of Townsend breakdown. The electrostatic particle-in-cell code Aleph utilizes the direct simulation Monte Carlo method to compute radiation transport. When there is strong radiation trapping, this approach is limited in that it must resolve the timescale associated with self-absorption. This renders many cases computationally intractable as sub-femtosecond time steps can be required to compute solutions for phenomena that occur over nanoseconds or microseconds. For two specific cases which exhibit strong radiation trapping, we find that spontaneous emissions having a frequency near the line center are inactive in the breakdown process and can be neglected. This enables larger time steps and a computational speedup of up to two orders of magnitude is observed. Some considerations for determining the validity of making such an approximation for Townsend breakdown problems and positive ionization wave problems are presented.

Roberds, N. A.↗