Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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 271 records · Page 15

NEML2: An efficient and modular multiphysics constitutive modeling library for hybrid computing environments

This paper presents NEML2, an open-source, high-performance library developed for constitutive material modeling, designed to support the flexible and modular development of models for complex material behavior. Building on the foundational structure of its predecessor, NEML, the NEML2 library introduces significant improvements, including enhanced vectorization, automatic differentiation, and seamless integration with PyTorch, facilitating the application of machine learning techniques in material simulations. NEML2 provides a C++ backend with Python bindings, enabling users to create custom material models that can be executed efficiently on both CPU and GPU platforms. The library also supports coupling with Multiphysics simulation frameworks like MOOSE, making it suitable for realistic simulations involving coupled physical processes. Rigorous quality assurance through unit and regression testing ensures the reliability of results, while the extensible, user-friendly design encourages collaboration and reproducibility across the scientific community. This paper provides an overview of NEML2’s architecture, core features, and applications, highlighting its impact on accelerating material qualification and advancing computational methods in materials science.

GPU↗

Deep Generative Models for Materials Discovery and Machine Learning-Accelerated Innovation

Machine learning and artificial intelligence (AI/ML) methods are beginning to have significant impact in chemistry and condensed matter physics. For example, deep learning methods have demonstrated new capabilities for high-throughput virtual screening, and global optimization approaches for inverse design of materials. Recently, a relatively new branch of AI/ML, deep generative models (GMs), provide additional promise as they encode material structure and/or properties into a latent space, and through exploration and manipulation of the latent space can generate new materials. These approaches learn representations of a material structure and its corresponding chemistry or physics to accelerate materials discovery, which differs from traditional AI/ML methods that use statistical and combinatorial screening of existing materials via distinct structure-property relationships. However, application of GMs to inorganic materials has been notably harder than organic molecules because inorganic structure is often more complex to encode. In this work we review recent innovations that have enabled GMs to accelerate inorganic materials discovery. We focus on different representations of material structure, their impact on inverse design strategies using variational autoencoders or generative adversarial networks, and highlight the potential of these approaches for discovering materials with targeted properties needed for technological innovation.

36 MATERIALS SCIENCE↗

Multiscale plasticity of geomaterials predicted via constrained optimization‐based granular micromechanics

Abstract A general framework to derive nonlinear elastic and elastoplastic material models from granular micromechanics is proposed, where a constraint‐based variational structure is introduced to classical grain contact‐based homogenization methods of hyperelasticity. Like the classical hyperelastic methods, reference solutions for closed‐form hyperelastic material models are analytically derived from the grain‐scale contact mechanics. However, unlike prior methods, the proposed homogenization framework defines closed‐form hyperelastoplastic material models that extend multiscale variational methods to granular plasticity. The proposed framework is used to develop novel granular micromechanics‐based macroscopic models for a Mises type solid, Drucker–Prager type plasticity, and grain‐contact cohesive‐debonding with a deviatorically and volumetrically coupled nonlinearly elastic response. Macroscopic plastic parameters and yield criteria are explicitly related to their microscale counterparts, for example, the friction coefficient governing intergranular slip. Numerical examples and comparison to measurements from the literature, including triaxial compaction of concrete, are provided to investigate model predictions and demonstrate calibration to experimental data.

Bryant, E. C.↗

Nuclear Material Process Modeling at the Y-12 National Security Complex

Dynamic simulation modeling is used at Y-12 to evaluate and forecast nuclear material inventories and production capacities to ensure that future supply can meet mission demand. Model outputs are analyzed by numerous Y-12 organizations and programs and coordinated with NNSA’s Office of Secondary Stage Production Modernization. Data-driven decisions for both short-term and long-term strategic planning for Y-12 mission execution are informed by the model. Dynamic simulation modeling capabilities for Y-12 nuclear material production continue to be expanded and refined.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Schema Elements for Granta Annual Report: FY2024

Granta: Materials Intelligence (Granta: MI) is a commercial database software distributed by Ansys, Inc. that is utilized by the Nuclear Security Enterprise (NSE) to organize and store relevant materials data. Lack of standard and well-documented database schema is the primary obstacle to an NSE materials data management solution, so the objective of this project is to create and document such a schema. In FY21, an approach for designing, documenting, and managing a standard database schema was described based on the creation of schema elements (collections of attributes used to describe particular aspects of the data) to be used as building blocks for creating various database tables without duplication. In FY22, these methods were applied through a multi-site collaboration to create and document the schema elements necessary to build a thermogravimetric analysis (TGA) testing table. In FY23 the schema was expanded to include elements for a differential scanning calorimetry (DSC) table, along with schema for supporting metadata tables including Instruments, Projects, Documents, and Testing Series. In FY24 the following progress was made, again through multi-site collaboration: • The existing schema elements were modified to accommodate thermomechanical analysis (TMA) data, and a table, Test Data: TMA, was created for managing TMA data. • The elements necessary for the following additive manufacturing (AM) data tables (directed at data specific to selective laser sintering AM technology) were created: • AM Builds • AM Processes • AM Part Designs • Built AM Parts • AM Feedstock Materials • AM Feedstock Material Batches • The elements necessary for creating a Calibrated Material Models table were created, and the Calibrated Material Models table was created. In FY25 the existing schema will be deployed on the production enterprise Granta instance on the enterprise secure network. Schema elements will be appended, and new elements created as necessary, to allow the creation of tables specifically to support materials testing, AM process development, and design and analysis for modernization programs.

36 MATERIALS SCIENCE↗

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma-catalyst modeling for materials selection: challenges and opportunities in nitrogen oxidation

Reliable kinetic models are essential to rationalize observations and to guide the selection of appropriate plasma-catalyst combinations. Kinetic models coupling plasma- and heterogeneous-catalytic-chemistries, however, are primitive. Here, we combine reduced plasma chemistries and density functional theory parameterized surface reactions to model NO formation from N2 oxidation, an alternative nitrogen fixation process. We first show plasma excited species can enhance turnover frequencies of both Pt and Au. The enhancing potential is a function of plasma species and their densities. We then compare NO production between plasma/Pt and plasma/Au with an integral reactor model and explore the dependence of NO production as a function of plasma conditions, reactor configurations and chemical compositions. Plasma/Pt generates more NO than plasma/Au under most thermal and plasma parametric conditions. Plasma/Au outperforms plasma/Pt under high plasma excitations at low temperatures. Furthermore, these results highlight that selection of plasma-catalyst combinations should consider the properties of plasmas, materials and the coupling of the two.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fatigue strength reduction model: RANDOM3 and RANDOM4 user manual. Appendix 2: Development of advanced methodologies for probabilistic constitutive relationships of material strength models

FORTRAN programs RANDOM3 and RANDOM4 are documented in the form of a user's manual. Both programs are based on fatigue strength reduction, using a probabilistic constitutive model. The programs predict the random lifetime of an engine component to reach a given fatigue strength. The theoretical backgrounds, input data instructions, and sample problems illustrating the use of the programs are included.

Boyce, Lola↗

Probabilistic constitutive relationships for cyclic material strength models

A methodology is developed that provides a probabilistic treatment for the lifetime of structural components of aerospace propulsion systems subjected to fatigue. Material strength degradation models, based on primitive variables, include both a fatigue strength reduction model and a fatigue crack growth model. Probabilistic analysis is based on simulation, and both maximum entropy and maximum penalized likelihood methods are used for the generation of probability density functions. The resulting constitutive relationships are included in several computer programs.

Boyce, L.↗

Benchmarking large language models for materials synthesis: The case of atomic layer deposition

In this work, we introduce an open-ended question benchmark, ALDbench, to evaluate the performance of large language models (LLMs) in materials synthesis, and, in particular, in the field of atomic layer deposition, a thin film growth technique used in energy applications and microelectronics. Our benchmark comprises questions with a level of difficulty ranging from the graduate level to domain expert current with the state of the art in the field. Human experts reviewed the questions along the criteria of difficulty and specificity, and the model responses along four different criteria: overall quality, specificity, relevance, and accuracy. We ran this benchmark on an instance of OpenAI’s GPT-4o. The responses from the model received a composite quality score of 3.7 on a 1–5 scale, consistent with a passing grade. However, 36% of the questions received at least one below average score. An in-depth analysis of the responses identified at least five instances of suspected hallucination. Finally, we observed statistically significant correlations between the difficulty of the question and the quality of the response, the difficulty of the question and the relevance of the response, the specificity of the question, and the accuracy of the response as graded by the human experts. Furthermore, this emphasizes the need to evaluate LLMs across multiple criteria beyond difficulty or accuracy.

Artificial intelligence↗

Parallel simulation via SPPARKS of on-lattice kinetic and Metropolis Monte Carlo models for materials processing

Abstract SPPARKS is an open-source parallel simulation code for developing and running various kinds of on-lattice Monte Carlo models at the atomic or meso scales. It can be used to study the properties of solid-state materials as well as model their dynamic evolution during processing. The modular nature of the code allows new models and diagnostic computations to be added without modification to its core functionality, including its parallel algorithms. A variety of models for microstructural evolution (grain growth), solid-state diffusion, thin film deposition, and additive manufacturing (AM) processes are included in the code. SPPARKS can also be used to implement grid-based algorithms such as phase field or cellular automata models, to run either in tandem with a Monte Carlo method or independently. For very large systems such as AM applications, the Stitch I/O library is included, which enables only a small portion of a huge system to be resident in memory. In this paper we describe SPPARKS and its parallel algorithms and performance, explain how new Monte Carlo models can be added, and highlight a variety of applications which have been developed within the code.

36 MATERIALS SCIENCE↗

Material Response Modeling of Ablative Thermal Protection Systems using PATO

Developing thermal protection systems (TPS) for future space vehicles involves an extensive design and test cycle. A key component of the cycle is determining margins for a safe design that rely on confidence in the performance of the TPS. Predicting the complicated multiphysics phenomena that occur during atmospheric entry requires high-fidelity modeling tools. To this end, the Porous-material Analysis Toolbox based on OpenFOAM (PATO) has been developed. PATO is an open-source software for computing material response of reactive porous materials submitted to high-temperature environments. Recent applications of PATO include computation of the full heatshield 3D material response from the Mars Science Laboratory atmospheric entry and ablation during arc jet testing. Current efforts are underway to loosely couple PATO with other discipline specialized codes including hypersonic computational fluid dynamics (CFD) to assess the effects of pyrolysis-gas blowing into the boundary layer and computational solid mechanics to address modeling of mechanical erosion. Surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation are also being added.

Thermal Protection Systems↗

Material Response Modeling of Ablative Thermal Protection Systems using PATO

Developing thermal protection systems (TPS) for future space vehicles involves an extensive design and test cycle. A key component of the cycle is determining margins for a safe design that rely on confidence in the performance of the TPS. Predicting the complicated multiphysics phenomena that occur during atmospheric entry requires high-fidelity modeling tools. To this end, the Porous-material Analysis Toolbox based on OpenFOAM (PATO) has been developed. PATO is an open-source software for computing material response of reactive porous materials submitted to high-temperature environments. Recent applications of PATO include computation of the full heatshield 3D material response from the Mars Science Laboratory atmospheric entry and ablation during arc jet testing. Current efforts are underway to loosely couple PATO with other discipline specialized codes including hypersonic computational fluid dynamics (CFD) to assess the effects of pyrolysis-gas blowing into the boundary layer and computational solid mechanics to address modeling of mechanical erosion. Surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation are also being added.

Thermal Protection Systems↗

Study of materials performance model for aircraft interiors

A demonstration version of an aircraft interior materials computer data library was developed and contains information on selected materials applicable to aircraft seats and wall panels, including materials for the following: panel face sheets, bond plies, honeycomb, foam, decorative film systems, seat cushions, adhesives, cushion reinforcements, fire blocking layers, slipcovers, decorative fabrics and thermoplastic parts. The information obtained for each material pertains to the material's performance in a fire scenario, selected material properties and several measures of processability.

Leary, K.↗

Generalizing the Gurson model using symbolic regression and transfer learning to relax inherent assumptions

Abstract To generate material models with fewer limiting assumptions while maintaining closed-form, interpretable solutions, we propose using genetic programming based symbolic regression (GPSR), a machine learning (ML) approach that describes data using free-form symbolic expressions. To maximize interpretability, we start from an analytical, derived material model, the Gurson model for porous ductile metals, and systematically relax inherent assumptions made in its derivation to understand each assumption’s contribution to the GPSR model forms. We incorporate transfer learning methods into the GPSR training process to increase GPSR efficiency and generate models that abide by known mechanics of the system. The results show that regularizing the GPSR fitness function is critical for generating physically valid models and illustrate how GPSR allows a high level of interpretability compared with other ML approaches. The method of systematic assumption relaxation allows the generation of models that address limiting assumptions found in the Gurson model, and the symbolic forms allow conjecture of decreased material strength due to void interaction and non-symmetric void shapes.

36 MATERIALS SCIENCE↗

Higher-order factorization machine for accurate surrogate modeling in material design

Efficient and robust optimization is important in material science for identifying optimal structural parameters and enhancing material performance. Surrogate-based active learning algorithms have recently gained great attention for their ability to efficiently navigate large, high-dimensional design spaces. Among surrogate models, 2 nd -order factorization machine (FM) models are widely employed as the surrogate model in active learning algorithms due to their balance between simplicity and effectiveness. However, their quadratic nature limits their capacity to capture complex, higher-order interactions among variables, often leading to suboptimal solutions. To overcome this limitation, we propose an active learning scheme integrating a 3 rd -order FM model, capable of modeling three-variable interactions and more intricate relationships in material systems. We comprehensively evaluate the surrogate modeling performance of the 3 rd -order FM case using various objective functions. Furthermore, we examine the optimization reliability and efficiency of the 3 rd -order FM-based active learning in a real-world material design task (e.g., nanophotonic structures for transparent radiative cooling). Our study shows that the 3 rd -order FM outperforms the 2 nd -order model in both surrogate accuracy and optimization performance, highlighting higher-order models’ promises for material design and optimization problems.

Factorization machine↗

HydraGNN_GFM_FineTuning4Materials v1.0

This repository enables fine-tuning of the HydraGNN Predictive GFM 2026 — an open-source ensemble of pre-trained graph foundation models for atomistic materials modeling, developed at Oak Ridge National Laboratory. The GFM 2026 is freely available and downloadable via Globus from the OLCF Data Constellation (DOI: 10.13139/OLCF/2562660). Starting from these pre-trained weights, this repository provides a complete transfer learning pipeline for adapting the GFM ensemble to domain-specific molecular and materials property prediction tasks. It includes: 1) Utilities for ensemble fine-tuning with task-specific output heads 2) Example pipelines for eight widely-used materials and molecular datasets 3) Tools for model adaptation and head configuration 4) Data preprocessing utilities for each supported dataset 5) Benchmarking and evaluation scripts

Ungerboeck, Linda↗

Hydrostatic Stress Effects Incorporated Into the Analysis of the High-Strain-Rate Deformation of Polymer Matrix Composites

Procedures for modeling the effect of high strain rate on composite materials are needed for designing reliable composite engine cases that are lighter than the metal cases in current use. The types of polymer matrix composites that are likely to be used in such an application have a deformation response that is nonlinear and that varies with strain rate. The nonlinearity and strain rate dependence of the composite response is primarily due to the matrix constituent. Therefore, in developing material models to be used in the design of impact-resistant composite engine cases, the deformation of the polymer matrix must be correctly analyzed. However, unlike in metals, the nonlinear response of polymers depends on the hydrostatic stresses, which must be accounted for within an analytical model. An experimental program has been carried out through a university grant with the Ohio State University to obtain tensile and shear deformation data for a representative polymer for strain rates ranging from quasi-static to high rates of several hundred per second. This information has been used at the NASA Glenn Research Center to develop, characterize, and correlate a material model in which the strain rate dependence and nonlinearity (including hydrostatic stress effects) of the polymer are correctly analyzed. To obtain the material data, Glenn s researchers designed and fabricated test specimens of a representative toughened epoxy resin. Quasi-static tests at low strain rates and split Hopkinson bar tests at high strain rates were then conducted at the Ohio State University. The experimental data confirmed the strong effects of strain rate on both the tensile and shear deformation of the polymer. For the analytical model, Glenn researchers modified state variable constitutive equations previously used for the viscoplastic analysis of metals to allow for the analysis of the nonlinear, strain-rate-dependent polymer deformation. Specifically, we accounted for the effects of hydrostatic stresses. An important discovery in the course of this work was that the hydrostatic stress effects varied during the loading process, which needed to be accounted for within the constitutive equations. The model is characterized primarily by shear data, with tensile data used to characterize the hydrostatic stress effects.

Goldberg, Robert K.↗