Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microstructure simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Model-informed, Adaptive Physical Vapor Deposition to Fabricate Hierarchical Binary-alloy Thin-films

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of microstructure on fatigue crack propagation in additive manufactured nickel-based superalloy Haynes 282: an experiment and crystal plasticity study

Haynes 282 is a γ' precipitation-strengthened nickel-based superalloy known for its exceptional high-temperature creep resistance and excellent fabricability. Recent advancements in powder-bed fusion-based additive manufacturing (PBF-AM) have enabled the fabrication of Haynes 282 with accurate, site-specific control of the grain orientation and morphology at the microscale. This ability opens up new avenues for microstructure design, to improve the material’s fatigue crack resistance and service life. This paper investigates the fatigue crack growth behavior of hybrid microstructure Haynes 282 fabricated via PBF-AM. Previous experiments revealed a higher crack propagation rate in the coarse columnar-grained microstructural regions when compared against fine-grained areas. Here, a strain gradient crystal plasticity model was adapted to study the fracture-related mechanical fields at the crack tip in the two microstructures. The simulation results showed a consistent influence of grain structure and texture on crack propagation, as was seen in the experiment. The model analysis revealed higher crack propagation driving force along crack direction in the coarse-grained sharply textured microstructure and higher driving force for crack kinking in fine-grained more diffusely textured microstructure, which is ascribed to the combined effect of yield stress, hardening rate, texture and grain morphology. The presented modeling approach will facilitate the development of the AM-based accurate microstructure design by deepening the fundamental study in AM-specific microstructure-properties relations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Age hardening response of Cu-30Ni alloys: The role of Si microalloying additions

Cu-30Ni alloys offer high strength and excellent corrosion resistance for maritime applications. While primarily a solid solution system, industrial alloys typically include microalloying additions of Fe, Mn, Nb and/or Si to enable age hardening. However, an understanding of the microstructural evolution during age hardening remains incomplete. Here, we conduct systematic ageing treatments and report the resulting hardness and microstructures. The Cu-30Ni alloy with Si aged at 650 °C for 6 h demonstrates significantly enhanced Vickers hardness, reaching ∼201 HV 1 compared to 103 HV 1 for the as-homogenised sample. Si-rich clusters and γ′ precipitates are identified, with their composition, size, and volume fraction quantified to determine their strengthening contributions. First-principles atomistic simulations elucidate the underlying formation mechanisms of these clusters and precipitates, highlighting the critical role of Si in driving their nucleation and growth. This study advances the design of high-strength, Cu alloys with the potential for corrosion resistance in demanding maritime environments.

Age hardening↗
DOE OSTI2022

A cycle-jump acceleration method for the crystal plasticity simulation of high cycle fatigue of the metallic microstructure

For metallic materials, high-cycle fatigue life is sensitive to underlying microstructure features including secondary phases, textures, grains morphology, etc. The traditional, data-based safe-life approaches for modeling fatigue don’t explicitly consider the microstructure and can’t guide study in microstructure modification for improved fatigue property. Crystal plasticity-based simulation provides increased model fidelity at the expense of immense computation time, making it inapplicable for high cycle fatigue. In this work, an acceleration method based on cycle-jump approach (Lesne and Savalle, 1989) was developed for microstructure-based high-cycle FE simulation using crystal plasticity constitutive-model. We report this method demonstrated high efficiency in benchmark tests of various conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Calibrating uncertain parameters in melt pool simulations of additive manufacturing

Melt pool scale numerical modeling of additive manufacturing (AM) processes can provide predictive capabilities and theoretical insight into the process-property-structure-performance relationships for AM parts. Despite capabilities of numerical models to solve complex multi-physics problems, it is often important to consider a tradeoff between detailed physics and computational cost. Therefore, sources of uncertainty in both experimental conditions and the parameters needed for modeling require models to be validated against empirical evidence. Here, a method is proposed to calibrate uncertain parameters used in continuum-scale melt pool models for powder bed fusion (PBF) AM. Both a simplified heat transfer model and a heat transfer and fluid flow model were investigated. A surrogate model and Markov chain-based optimization algorithm calibrated melt pool geometry for models within experimental variation of the target melt pool width and depth from the NIST AM-Bench 2018-02 dataset. The melt pool temperature distributions, solidification parameters, and simulated multi-layer solidification microstructures were compared between the two models. Similar results from both models indicate that calibrated, lower fidelity numerical models may be used in place of higher fidelity models to generate melt pool solidification data. Finally, these calibrated models therefore enable lower computational cost melt pool simulations without a noticeable decrease in simulation accuracy for grain-scale microstructure simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Under the microscope: Reduced activation ferritic martensitic steel Eurofer-97 following ion-Irradiation and high-temperature high-pressure water exposure

Here, this study is designed to characterise the microstructural behaviour of Eurofer-97 steel under ion irradiation and subsequent exposure to high-temperature high-pressure (HTHP) water. Eurofer-97, a ferritic-martensitic steel, has been developed to withstand the conditions of fusion reactors in the locations in contact with coolant with an elevated level of neutron flux, such as the breeder-wall blanket. The material has been studied after self-ion irradiation (using Fe ions) simulating the microstructural effects of neutron irradiation limited to the subsurface layer. The corrosion properties of the Eurofer-97 steel were studied by exposure to HTHP water up to 331 °C. Advanced microstructural characterisation using scanning, transmission electron and focused ion beam microscopy was performed on the as-received microstructure and after ion irradiation. This was then characterised after exposure for 240 h in high-temperature water. Eurofer-97 had a dense, columnar Cr-rich inner oxide, followed by a Fe-rich outer oxide layer. In the irradiated condition the grain structure and oxide itself was less ordered. No appreciable difference in oxide thickness was identified between the irradiated and unirradiated specimens after this short exposure time.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructure and creep strength of simulated intercritical heat-affected zone of grade 91 steel

The intercritical heat affected zone (ICHAZ) has been reported as one of the most Type IV cracking susceptible regions in 9Cr creep-resistant steel weldments. However, creep degradation mechanisms within the ICHAZ itself need further clarifications. In this work, two ICHAZ specimens of Grade 91 steel, low-temperature ICHAZ (LT-ICHAZ) and high-temperature ICHAZ (HT-ICHAZ), were simulated using the Gleeble thermomechanical system by exposing to two peak temperatures (860 °C and 900 °C) between $A_{C1}$ and $A_{C3}$. Dramatically different creep strengths of two simulated ICHAZs were observed and studied. We report dilation curve analysis indicates a high fraction of newly transformed martensite formed in the HT-ICHAZ, which results in a much higher hardness (360 HV0.5) of the HT-ICHAZ than 266 HV0.5 of LT-ICHAZ. Precipitates, especially $M_{23}C_6$ carbides, were not fully dissolved in both ICHAZ specimens. After a typical postweld heat treatment (760 °C-2 hours), the faster recovery of low-carbon martensite and reduced precipitation strengthening due to $M_{23}C_6$ carbides coarsening in the HT-ICHAZ led to a significant hardness reduction. These microstructural degradations in the HT-ICHAZ made its creep lifetime about 36 times shorter than that of the LT-ICHAZ tested at 650 °C. The remaining tempered martensite from base metal in the LT-ICHAZ was the primary contributor to maintain its high creep resistance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗
DOE OSTI2024

GrainGNN: A dynamic graph neural network for predicting 3D grain microstructure

We propose GrainGNN, a surrogate model for the evolution of polycrystalline grain structure under rapid solidification conditions in metal additive manufacturing. High fidelity simulations of solidification microstructures are typically performed using multicomponent partial differential equations (PDEs) with moving interfaces. The inherent randomness of the PDE initial conditions (grain seeds) necessitates ensemble simulations to predict microstructure statistics, e.g., grain size, aspect ratio, and crystallographic orientation. Here, currently such ensemble simulations are prohibitively expensive and surrogates are necessary.In GrainGNN, we use a dynamic graph to represent interface motion and topological changes due to grain coarsening. We use a reduced representation of the microstructure using hand-crafted features; we combine pattern finding and altering graph algorithms with two neural networks, a classifier (for topological changes) and a regressor (for interface motion). Both networks have an encoder-decoder architecture; the encoder has a multi-layer transformer long-short-term-memory architecture; the decoder is a single layer perceptron.We evaluate GrainGNN by comparing it to high-fidelity phase field simulations for in-distribution and out-of-distribution grain configurations for solidification under laser power bed fusion conditions. GrainGNN results in 80%–90% pointwise accuracy; and nearly identical distributions of scalar quantities of interest (QoI) between phase field and GrainGNN simulations compared using Kolmogorov-Smirnov test. GrainGNN's inference speedup (PyTorch on single x86 CPU) over a high-fidelity phase field simulation (CUDA on a single NVIDIA A100 GPU) is 150×–2000× for 100-initial grain problem. Further, using GrainGNN, we model the formation of 11,600 grains in 220 seconds on a single CPU core.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Three-dimensional microstructure-explicit and void-explicit mesoscale simulations of detonation of HMX at millimeter sample size scale

Fully three-dimensional (3D) microstructure-explicit and void-explicit mesoscale simulations of the shock-to-detonation (SDT) process of pressed granular HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine) are performed. The overall size scale of the models is up to 3 × 3 × 15 mm3, with ∼30 000 grains and 206 265 voids. The models account for the heterogeneous material microstructure, constituent distribution, constituent morphology, and voids. Loading conditions considered involve piston velocities in the range of 600–1200 m/s or pressures in the range of 4–8 GPa. The focus is on analyzing the SDT process and the effects of microstructure and voids on the run-to-detonation distance (RDD). Companion two-dimensional (2D) simulations are also carried out to assess the differences between 2D and 3D. Statistically equivalent microstructure sample sets (SEMSSs) are generated and used for both 2D and 3D, allowing the prediction of the statistical and probabilistic Pop plots (PPs). The predictions are in general agreement with trends in available experimental data in the literature. It is found that both the microstructure (heterogeneous grain size, morphology, and size distribution) and voids significantly affect the RDD and the PPs. These effects are systematically delineated and quantified via the use of SEMSSs with different combinations of attributes. A recently developed probabilistic formulation for the PPs is used to characterize the results, allowing uncertainties in the relations between the shock pressure and RDD arising from material heterogeneities to be quantified. The probabilistic formulation is further used to quantify the confidence levels in the ranked order of influences of different combinations of microstructure and voids on the PPs.

Miller, Christopher (ORCID:0000000206844381)↗
DOE OSTI2022

Mechanical Behavior of Additively Manufactured Molybdenum and Fabrication of Microtextured Composites

Refractory metals are a class of high-melting-temperature materials suitable for use in extreme environment applications. Interestingly, during additive manufacturing many pure refractory metals exhibit a switch from (001) to (111) build direction fiber preference with increasing surface energy density. Here we exploit this solidification physics to fabricate material with “mesoscale composite” engineered structures consisting of features with contrasting (001) and (111) build direction microtextures. Separately, elevated temperature tensile testing of EBM fabricated material with a randomized distribution of mixed (001)/(111)-fiber grains is shown to exhibit excellent properties. These results are utilized to build a crystal plasticity model for evaluating the local inelastic response of the composite mesoscale structures. Analysis of printed microstructures and microstructure-scale simulations indicate that both macro-scale and localized material behavior may be tailored. This strategy can be potentially used to synthesize materials with optimized performance for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Text Mining for Process–Structure–Properties Relationships in Metals

With the advent of large language models (LLMs), the vast unstructured text within millions of academic papers is increasingly accessible for materials discovery—although significant challenges remain. While LLMs offer promising few- and zero-shot learning capabilities, particularly valuable in the materials domain where expert annotations are scarce, general-purpose LLMs often fail to address key materials-specific queries without further adaptation. To bridge this gap, fine-tuning LLMs on human-labeled data is essential for effective structured knowledge extraction (Liu in The Importance of Human-Labeled Data in the Era of LLMs, 2023). Here, in this study, we introduce a novel annotation schema designed to extract generic process–structure–properties relationships from scientific literature. We demonstrate the utility of this approach using a dataset of 128 abstracts, with annotations drawn from two distinct domains: high-temperature materials (Domain I) and uncertainty quantification in simulating materials microstructure (Domain II). Initially, we developed a conditional random field (CRF) model based on MatBERT—a domain-specific BERT variant—and evaluated its performance on Domain I. Subsequently, we compared this model with a fine-tuned LLM (GPT-4o from OpenAI) under identical conditions. Our results indicate that fine-tuning LLMs can significantly improve entity extraction performance over the BERT-CRF baseline on Domain I. However, when additional examples from Domain II were incorporated, the performance of the BERT-CRF model became comparable to that of the GPT-4o model. These findings underscore the potential of our schema for structured knowledge extraction and highlight the complementary strengths of both modeling approaches.

Materials science↗
DOE OSTI2023

Flow Strength Measurements of Wrought and AM SS304L via Pressure Shear Plate Impact Experiments

Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Part-scale microstructure prediction for laser powder bed fusion Ti-6Al-4V using a hybrid mechanistic and machine learning model

Laser powder bed fusion (LPBF) Ti-6Al-4V is widely studied for use in structural applications in aerospace and medical industries, but mechanical anisotropy and microstructural inhomogeneity prohibits its wider adoption. Although successful microstructure prediction models have been developed, a remaining challenge is their limited integration across length/time scales and validation by experimental studies. Here, this work proposes a physics-augmented machine learning surrogate model to unite predictions of LPBF temperature, β phase morphology and texture, and α/α’ formation into a single framework that is calibrated and validated with experiments. First, a phase field (PF) model of the martensitic β→α’ transformation is developed and calibrated using data from in-situ synchrotron cyclic heating/cooling studies quantifying the variation of α phase fraction with time. In parallel, an established finite difference-Monte Carlo (FDMC) model predicts the part-scale temperature profile and β grain formation during solidification. A dataset is developed using LPBF cyclic temperature descriptors from the FDMC model as inputs and corresponding α/α’ phase fraction and width from the PF model as outputs. Five machine learning (ML) regression models are tested and optimized, having mean absolute error in testing ≤ 4 %, and the k-nearest neighbors (KNN) model is selected as the best performing. The KNN model is called at the nodal level during post-processing of the FDMC model to replace and downscale the response of the PF model. The combined agility and accuracy of the hybrid FDMC-ML model enables part-scale microstructure predictions that can be further used for property predictions to accelerate AM process optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparison of three state-of-the-art crystal plasticity based deformation twinning models for magnesium alloys

In magnesium alloys, deformation twinning and its interactions with dislocation slip are responsible for a sigmoidal shape stress–strain behavior and an asymmetrical tension–compression yield strength in magnesium alloys. The sensitivity of twinning to the underlying microstructure renders the crystal plasticity method the most commonly adopted modeling approach for magnesium-twinning. This paper compares three state-of-the-art crystal plasticity-based twinning models from the literature, namely the elastic-viscoplastic self-consistent twinning-detwinning (EVPSC-TDT) model, crystal plasticity finite element model based on enhanced predominate twin reorientation approach (CPFE-ePTR), and the crystal plasticity finite element model based on “discrete twinning” approach (CPFE-DT). A polycrystalline microstructure is simulated with all three methods to compare the resulting stress–strain curves and lattice strains to those from the experimentally measured data. All three methods showed the capability of capturing the experimental results with different levels of accuracy. Additionally, the EVPSC-TDT method avoids solving the finite element matrices and showed the highest computational efficiency. The CPFE-ePTR model shows a higher accuracy in capturing the lattice strain. The CPFE-DT relies on high-resolution finite element mesh and is much slower than the other two methods, but it captured the local deformation concentration and stress reversal phenomena near the twin band, which was not possible with the other two methods. Based on the comparison, guidance for the selection of the appropriate model based on the specific modeling target is provided in this paper.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Data-driven Whitney forms for structure-preserving control volume analysis

Control volume analysis models physics via the exchange of generalized fluxes between subdomains. Here, we introduce a scientific machine learning framework adopting a partition of unity architecture to identify physically-relevant control volumes, with generalized fluxes between subdomains encoded via Whitney forms. The approach provides a differentiable parameterization of geometry which may be trained in an end-to-end fashion to extract reduced models from full field data while exactly preserving physics. The architecture admits a data-driven finite element exterior calculus allowing discovery of mixed finite element spaces with closed form quadrature rules. An equivalence between Whitney forms and graph networks reveals that the geometric problem of control volume learning is equivalent to an unsupervised graph discovery problem. The framework is developed for manifolds in arbitrary dimension, with examples provided for H(div) problems in $\mathbb{R}$ establishing convergence and structure preservation properties. Finally, we consider a lithium-ion battery problem where we discover a reduced finite element space encoding transport pathways from high-fidelity microstructure resolved simulations. The approach reduces the 5.89M finite element simulation to 136 elements while reproducing pressure to under 0.1% error and preserving conservation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Identifying multiple synergistic factors on the susceptibility to stress relaxation cracking in variously heat-treated weldments

The 347H austenitic stainless steel has been widely used for pressure vessels and pipeline (PVP) applications due to its excellent creep and corrosion resistance, which fit ideally to the harsh conditions in petrochemical industries, fossil fuel or nuclear power plants, and modern energy storages. However, a failure mode has been commonly observed with cracks emerging at the heat affected zone (HAZ) of weldments during post-weld heat treatment (PWHT) or under intermediate to high temperature service conditions. This phenomenon is termed as Stress Relaxation Cracking (SRC) since the purpose of PWHT is to relieve the welding-induced residual stress fields, or as Stress Age Cracking (SAC) if failure happens during service. A leading literature explanation of this failure suggests that the residual stress relaxation and the precipitation dissolution and/or re-precipitation occur in the same temperature range, which can lead to locally high strains and thus to crack at the grain boundaries. Since in situ spatial measurements of residual stress fields, microstructural evolution, and failure processes are nearly infeasible, this work recourses to a micromechanical finite element framework that models the high temperature failure as the nucleation and growth of grain boundary cavities, whereas various parameters such as thermomechanical loading history and its evolution, the competition of grain-interior dislocation creep and grain-boundary diffusion in failure lifetime, and microstructural heterogeneities (such as the precipitate free zone near grain boundaries) can be quantitatively incorporated. It can be concluded from these microstructure-explicit simulations that an accurate knowledge of residual stress evolution and a carefully calibrated set of material constitutive parameters are the essential prerequisites for lifetime predictions. The understanding of individual governing factors also leads to a mechanistic interpretation of the observed SRC susceptibility C-curves. In conclusion, these results suggest that the criticality of residual stress evolution, but not the precipitation-induced local strains, be the leading factor for SRC.

347H stainless steel weldments↗