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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

LASSO for CALPHAD Model Selection Enables Data-Efficient Thermodynamic Modeling: An Application in Thermochemical Hydrogen Production Materials

Phenomenological CALPHAD (CALculation of PHAse Diagrams) models, widely used for multicomponent materials, often contain a considerable number of parameters and require fitting using data from a relatively small number of experimental measurements or theoretical calculations. Sometimes these parameters are introduced for the purpose of improving model fits but without clear physical justification, which leads to overparametrized models with poor generalization performance. Automated approaches for optimal model selection based on the available data therefore become critical. Here, in this work, a least absolute shrinkage and selection operator (LASSO)-based approach is developed for model selection by leveraging the linearity of the CALPHAD model with respect to its parameters to convert the model selection and fitting to a LASSO minimization problem. We demonstrate its utility for thermodynamic modeling of thermochemical hydrogen (TCH) production materials using lanthanum strontium manganite (LSM) as an example. Various TCH-relevant properties, including oxygen stoichiometry as a function of oxygen partial pressure, enthalpy of reduction, and entropy of reduction, are successfully predicted with reasonable accuracy using a minimal set of model parameters. Importantly, the model selection and fitting involve minimal human decision; it can therefore be applied to high-throughput DFT defect calculations and yield efficient workflows for TCH material modeling and optimization.

CALPHAD↗
DOE OSTI2022

Accelerating the discovery of battery electrode materials through data mining and deep learning models

The availability of crystalline materials databases allows for building accurate machine learning (ML) models that can accelerate the exploration of materials chemical space for energy storage applications. In this work, we screen all inorganic materials included in the Materials Project and AFLOW databases as potential metal-ion battery electrodes. We develop an efficient protocol to mine and screen raw data in current databases and provide a new database of electrode materials by considering pairs of charged and discharged electrodes. This effort leads to a new database with over 190,000 instances, in contrast to the original battery database which contains about 5000. The expanded battery data set is then used to build regression-based deep neural network models for predicting average voltages and percentage volume changes upon charging and discharging, which present improvements of at least 28% for target properties with respect to previous models, and are now able to predict anode electrodes (low voltage region) as well as electrodes that will not work in electrochemical cells (negative voltages), overcoming the challenges identified in previous ML models for battery electrodes. Additionally, a further screening of the expanded database itself allowed us to identify 35 novel electrode candidates with excellent battery performance metrics.

25 ENERGY STORAGE↗
DOE OSTI2024

Next-Generation Materials Design: Quantum Mechanics and Data-Driven Modeling

The future of materials design is rapidly advancing through the combination of quantum mechanics and data-driven modeling. These approaches integrate quantum principles with advanced data analysis, enabling precise insights into material behavior. This talk will highlight recent progress in using these methods for computational design, particularly in high-entropy alloy catalysts, emphasizing the role of hierarchical machine-learning architectures for accurate predictions. Additionally, I will discuss our work on developing machine learning interatomic potentials (MLPs) for single-element metals, metal oxides, and alloys under extreme conditions, focusing on melting behavior and phase properties at high temperatures and pressures. We have also refined our MLP models to capture dynamic surface interactions, such as CO2 and CO adsorption on MgO, using both static and molecular dynamics simulations. These models maintain high accuracy while significantly reducing computational costs compared to first-principles calculations. By enabling efficient and accurate simulations, this work supports broader community adoption, optimizes datasets for materials discovery, and extends the accessible time, size, and environmental conditions beyond the limits of experiments and traditional simulations.

machine learning↗
DOE OSTI2024

Leveraging a Neural Network-Enhanced Reproducing Kernel Particle Method for Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and resulting in reduced performance and service life. A coupled electro-chemo-mechanical reproducing kernel particle method (RKPM) formulation has been developed to analyze this system. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based model construction by RKPM is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. Here, a neural network-enhanced reproducing kernel particle method (NN-RKPM) [1, 2] is introduced to effectively model damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RKPM is additionally used to inform how crack opening and closure in turn affect the coupled chemical equations and material microstructure. Reference: [1] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, pp 4422-4454, https://doi.org/10.1002/nme.7040, 2022. [2] Baek, J., Chen, J. S., "A Neural Network-Based Enrichment of Reproducing Kernel Approximation for Modeling Brittle Fracture", Computer Methods in Applied Mechanics and Engineering Vol. 410, 116590, 2024.

degradation↗
DOE OSTI2024

Equivariant graph convolutional neural networks for the representation of homogenized anisotropic microstructural mechanical response

Composite materials with different microstructural material symmetries are common in engineering applications where grain structure, alloying and particle/fiber packing are optimized via controlled manufacturing. In fact these microstructural tunings can be done throughout a part to achieve functional gradation and optimization at a structural level. To predict the performance of particular microstructural configuration and thereby overall performance, constitutive models of materials with microstructure are needed. In this work we provide neural network architectures that provide effective homogenization models of materials with anisotropic components. These models satisfy equivariance and material symmetry principles inherently through a combination of equivariant and tensor basis operations. We demonstrate them on datasets of stochastic volume elements with different textures and phases where the material undergoes elastic and plastic deformation, and show that the these network architectures provide significant performance improvements.

anisotropy↗
DOE OSTI2021

Structure Development of the Interphase between Drying Cellulose Materials Revealed by In Situ Grazing-Incidence Small-Angle X-ray Scattering

The nano- to microscale structures at the interface between materials can define the macroscopic material properties. These structures are extremely difficult to investigate for complex material systems, such as cellulose-rich materials. The development of new model cellulose materials and measuring techniques has opened new possibilities to resolve this problem. We present a straightforward approach combining micro-focusing grazing-incidence small-angle X-ray scattering and atomic force microscopy (AFM) to investigate the structural rearrangements of cellulose/cellulose interfaces in situ during drying. Based on the results, we propose that molecular interdiffusion and structural rearrangement play a major role in the development of the properties of the cellulose/cellulose interphase; this model is representative of the development of the properties of joint/contact points between macroscopic cellulose fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Heat pump water heater enhanced with phase change materials thermal energy storage: Modeling study

A promising solution to improve the first hour rating (FHR) of a heat pump water heater (HPWH) involves employing a secondary tank which contains phase change material (PCM) capsules. To better understand the influence of PCM thermal storage on the HPWH operational performance, a dynamic model was developed to simulate and analyze the behavior of a newly developed HPWH technology that incorporates PCM storage into a standard HPWH to optimize key parameters such as the uniform energy factor and FHR. Mathematical models of several key components of the proposed HPWH-PCM integrated thermal energy storage (TES) system, e.g., water heater tank, PCM TES tank, evaporator, compressor, and expansion valve, have been elaborated. Also, a model-based control co-simulation platform was developed to integrate a embed PCM storage HPWH dynamic model with a control model for better supporting control design, analysis, verification, and validation. The model accuracy has been validated through comparing simulation results with lab test results, with a mean average percentage error of <5.5% for most of the selected performance variables. In addition, using the developed co-simulation platform, the demand response control strategy was studied to evaluate the load flexibility of the combined HPWH-PCM storage system by shifting system power usage to outside of the 3.5 h peak load period.

25 ENERGY STORAGE↗
DOE OSTI2020

Emerging materials intelligence ecosystems propelled by machine learning

We report that the age of cognitive computing and artificial intelligence (AI) is just dawning. Inspired by its successes and promises, several AI ecosystems are blossoming, many of them within the domain of materials science and engineering. These materials intelligence ecosystems are being shaped by several independent developments. Machine learning (ML) algorithms and extant materials data are utilized to create surrogate models of materials properties and performance predictions. Materials data repositories, which fuel such surrogate model development, are mushrooming. Automated data and knowledge capture from the literature (to populate data repositories) using natural language processing approaches is being explored. The design of materials that meet target property requirements and of synthesis steps to create target materials appear to be within reach, either by closed-loop active-learning strategies or by inverting the prediction pipeline using advanced generative algorithms. AI and ML concepts are also transforming the computational and physical laboratory infrastructural landscapes used to create materials data in the first place. Surrogate models that can outstrip physics-based simulations (on which they are trained) by several orders of magnitude in speed while preserving accuracy are being actively developed. Automation, autonomy and guided high-throughput techniques are imparting enormous efficiencies and eliminating redundancies in materials synthesis and characterization. The integration of the various parts of the burgeoning ML landscape may lead to materials-savvy digital assistants and to a human-machine partnership that could enable dramatic efficiencies, accelerated discoveries and increased productivity. Here, we review these emergent materials intelligence ecosystems and discuss the imminent challenges and opportunities. The materials research landscape is being transformed by the infusion of approaches based on machine learning. This Review discusses the emerging materials intelligence ecosystems and the potential of human-machine partnerships for fast and efficient virtual materials screening, development and discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Kernel Enriched Meshfree Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Laboratory of the Rockies (NLR), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. References: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 Mark-I Station Blackout: Part II – Uncertainty Analysis

Single case comparisons between severe accident simulations can provide detailed insights into severe accident model behavior, however, they cannot offer insights into model uncertainty, sensitivity to uncertain parameters, or underlying model biases.Here in this analysis, the single case benchmark comparison of the MELCOR material interaction models for a station blackout (SBO) scenario of a boiling water reactor (BWR) using representative Fukushima Daiichi Unit 1 boundary conditions is expanded to include an uncertainty analysis. As part of this uncertainty analysis, 1200 simulations are performed for each material interaction model (2400 total), with random sampling of 14 uncertain MELCOR input parameters. Input parameters are selected for their impact on models representing core degradation processes. These include candling, fuel rod failure, debris quenching and dryout. The analysis performed here is not a traditional “best-estimate” uncertainty analysis that uses best-estimate parameters or identifies best-estimate figure of merit distributions. Instead, it is an exploratory uncertainty analysis that identifies and interrogates underlying model form biases of the two material interaction models (eutectics and interactive materials models). Uniform distributions are applied to all uncertain parameters to ensure coverage of the model parameter uncertainty space. Key findings from this study include underlying model form biases exhibited by material interaction models, and notable differences in accident progression outcomes between the material interaction models. This uncertainty study extends and confirms the conclusions from the first part of this study, which compared the impact of material interaction modeling on simulation of a short-term station blackout scenario with representative Fukushima Daiichi Unit I boundary conditions. In particular, this study confirms that the eutectics model generally exhibits accelerated degradation and failure of fuel components, the core plate, and the lower head. The eutectics model also has a tendency to exhibit a greater degree of core degradation, greater debris mass formation, and larger debris mass ejection. Finally, the eutectics model exhibits higher maximum temperatures for fuel, cladding, particulate debris, oxidic molten pool, and metallic molten pool components than the interactive materials model; interactive materials model simulations exhibit a soft “limitation” on maximum temperatures that is related to the temperature at which material relocation occurs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Design and Analysis of a Mutual Inductance Level Sensor for Liquid Metals

Here, this article describes the design and analysis of an electromagnetic level sensor for use in high-temperature liquid metal systems. The mutual inductance level sensor (MILS) described in this work was fabricated using two single-conductor mineral insulated cables wrapped in a bifilar fashion around a stainless steel tube core and was housed in an isolating thimble that preserved the pressure boundary of the test vessel. Two sensor variations were fabricated that differ only in active length, 1016 and 1778 mm. Experimental data were collected using the 1016-mm sensor (MILS-MKII-040) that demonstrated a sensitivity of 9.2 μ V/mm in a room temperature testing stand that used solid aluminum as a surrogate for liquid metal. Experimental data were collected using the 1778-mm sensor (MILS-MKII-070) that demonstrated a sensitivity of 6.9 μ V/mm in the high-temperature (300 ° C) liquid sodium environment at the mechanisms engineering test loop (METL) of Argonne National Laboratory. The sensor performance was found to be repeatable over the course of several months, with roughly ±1% deviation from nominal output. Finite element models were developed in COMSOL Multiphysics that fully describe each test setup, and the models were validated using experimental data. The validated COMSOL models were used to perform an array of analyses that examined the performance of the sensor in differing environments. Maximizing the coil diameter inside the isolating thimble was found to maximize the signal and sensitivity of the sensor. An optimal operating frequency was found near 1000 Hz using both experimental data and COMSOL. The influence of a metallic thimble surrounding the sensor and a metallic sensor core was quantified and found to be negligible at the optimal operating frequency. The sensitivity of the sensor was quantified when monitoring the level of additional liquid metals. These include lead, lead-bismuth eutectic (LBE), sodium-potassium alloy (NaK), and lithium (in addition to sodium). The sensitivities were quantified using liquid metal properties at 350 ° C and 650 ° C. The geometry of the test stand model, all material properties used in the model, and the results are presented in a manner that allows the reader can replicate the model and perform additional analyses.

COMSOL↗
DOE OSTI2022

Advanced Physics Models for Particle-Particle Interactions (UCD final report)

The LDRD "Advanced Physics Models for Particle-Particle Interactions” collaborative research with Professor Jennifer Sinclair Curtis and PhD student Chanh Nguyen at the University of California, Davis involved meso-scale, discrete element method (DEM) modeling of micron-sized particles impacting a substrate. Using DEM, physical insight into the details of particle response upon impact can be obtained. Such details provide the type of information necessary to enhance existing continuum or macro-scale simulation tools so that modeling for material strength and particle breakage models is based on particle physics versus empirical relationships. While molecular dynamics simulations of atoms can also provide physical insight into the particle response to impact, the scale of particles that are treated is on the order of angstroms, not microns. This study gives evidence to the ability of the DEM simulation technique to provide such fundamental, detailed insight for micron-sized particles as a function of particle properties and impact velocity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Additive Manufacturing Evaporative Casting

Traditional lost foam casting has been around for decades. The process uses foam forms blown by an injection molding like process or CNC milled into the desired shape, then pouring molten metal over the foam pattern to create a metallic object. Additive Manufacturing Evaporative Casting (AMEC) is a new process that eliminates foam forms by using additive manufacturing to 3D print the desired cast geometry. This not only saves time and money but allows for more advanced and complex designs that can’t be achieved by carving foam. Additionally, AMEC doesn’t require molds and tooling like other casting and foundry options. Because the AMEC process is new, extensive testing is needed to develop a better understanding of the process to minimize defects, quantify material properties, and start computer modeling for the process. This CRADA (collaborative research and development agreement) between ORNL and Skuld seeks to improve the process, develop a computational model, characterize material properties, and explore new applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty

Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Superconductor to exciton condensate transition in a model copper-oxide material

Abstract Superconductivity and exciton condensation are fundamental phenomena in condensed matter physics, associated with the condensation of electron–electron and electron–hole pairs, respectively, into coherent quantum states. In this study, we present evidence of a superconductor to exciton condensate transition within the context of the three-band Hubbard model of copper-oxide-like materials. As the electron–electron repulsion increases, the superconducting phase is superseded by exciton condensation. In support of theoretical predictions—not yet realized experimentally—we observe the coexistence of the two condensates in the vicinity of the transition where the quantum states become a superposition of electron–electron and electron–hole condensates. Coexistence is rigorously computed from large eigenvalues and their eigenvectors in both the two-electron reduced density matrix (2-RDM) and the particle-hole RDM, which we obtain from a direct variational ground-state energy minimization with respect to the 2-RDM by semidefinite programming. We further discern that adjacent d orbitals and intervening p orbitals facilitate electron–electron pairing between copper orbitals, thereby supporting the superexchange mechanism for superconductivity. These observations suggest the feasibility of witnessing a superconductor to exciton condensate transition in copper-oxide analogs, bearing significant implications for identifying materials conducive to efficient transport processes.

Schouten, Anna O. (ORCID:0000000264255219)↗