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

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Simulating Micrometeoroid Bombardment on Mercury: Identifying New Space Weathering Features in the Laboratory

Airless bodies are continually exposed to the harsh interplanetary space environment, causing the alteration of their surfaces through a process known as space weathering (SW). SW is dominated by solar wind irradiation and micrometeoroid bombardment, which together alter the spectral, microstructural, and chemical characteristics of grains on the surface of airless bodies. The effects of SW are variable and depend on the heliocentric distance and the initial composition of the planetary surface, among other factors. On the Moon and S- type asteroids, it is well established that SW darkens and reddens the surface of these bodies and attenuates absorption bands across the Vis-NIR wavelengths. These spectral effects are largely the result of the production of metallic Fe nanoparticles (npFe). However, Mercury’s SW environment is unique compared to the Moon and S-type asteroids. Mercury is a geochemical endmember, with a surface composition low in Fe (<2 wt.%) and enriched in volatile components, including regions of the surface hypothesized to contain up to 4 wt.% graphite in the low reflectance material (LRM). In fact, the presence of graphite on Mercury has been hypothesized to act as a reducing agent for silicates during SW to produce Si-bearing Fe-rich metal. In addition, its location in the solar system exposes Mercury to an extreme SW environment, with the surface of the planet experiencing an intense solar wind flux and higher flux and velocity of micrometeoroid impactors compared to the Moon and S-type asteroids. The effects that such harsh SW has on materials with the unique compositional characteristics of Mercury are not well constrained. To better understand SW at Mercury, we must investigate these processes in the laboratory. Here, we present the results of our analyses of the spectral, microstructural and chemical characteristics of Mercury analog samples irradiated by pulsed laser to simulate the short duration, high temperature events associated with micrometeoroid impacts.

Nick Bott↗
DOE OSTI2020

Spatial correlation functions of paracrystals with radial symmetry

Here, we develop a phenomenological model to describe the structure of radially symmetric paracrystals whose long-range order are destroyed by propagation of particle fluctuations. General expressions are derived for the spatial correlation functions in one-, two-, and three-dimensional spaces. And the spatial correlation in paracrystals in reciprocal space is further discussed and clarified. The developed method can be used to quantitatively analyze the microstructure of paracrystalline materials in both real and reciprocal spaces via scattering experiments and computer simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anisotropic physics-regularized interpretable machine learning of microstructure evolution

Anisotropic Physics-Regularized Interpretable Machine Learning Microstructure Evolution (APRIMME) is a general-purpose machine learning solution for grain growth simulations. In prior work, PRIMME employed a deep neural network to predict site-specific migration as a function of its neighboring sites to model normal, isotropic, grain growth behavior. This work aims to extend this method by incorporating grain boundary misorientation-based grain growth behavior. APRIMME is trained on anisotropic simulations created using the Monte Carlo-Potts (MCP) model. Furthermore, the results of this work are compared statistically using grain radius, number of sides per grain, mean neighborhood misorientations, and the standard deviation of triple junction dihedral angles, and are found to match in most cases. The exceptions are small and seem to be related to two causes: (1) the deterministic model of APRIMME is learning from the stochastic simulations of MCP, which seems to accentuate triple junction behaviors; and, (2) a bias against very small grains is made evident in a quicker decrease in grains than expected at the beginning of an APRIMME simulation. APRIMME is also evaluated for its general ability to capture anisotropic grain growth behavior by first investigating different test case initial conditions, including a circle grain, three grain, and hexagonal grain microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Large Area Detection of Microstructural Defects with Multi-Mode Ultrasonic Signals

Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Austenite formation and cementite dissolution during intercritical annealing of a medium-manganese steel from a martensitic condition

In the present study, in situ high-energy X-ray diffraction (HEXRD), dilatometry, and metallography were conducted for the purpose of understanding microstructure evolution during intercritical annealing of a Fe-0.19C-4.39Mn steel from a martensitic initial condition. Three different simulations for austenite growth and solute partitioning during isothermal holding were conducted using the DICTRA™ module of Thermo-Calc®. The microstructures after intercritical annealing exhibit film-like retained austenite; however, in situ HEXRD indicates that retained austenite in the initial martensitic microstructure had decomposed upon heating to the isothermal holding temperature, suggesting that film-like austenite may be generated during intercritical annealing via nucleation and growth of new austenite and without necessarily preserving initial retained austenite films in the martensitic microstructure. Metallography and in situ HEXRD also indicate that cementite had formed upon heating but did not readily dissolve during an 1800 s isothermal hold although considerable growth of austenite had occurred during the isothermal hold, suggesting that austenite nucleation and initial growth during intercritical annealing does not solely occur due to carbon partitioning. A simulation incorporating cementite suggests that the increase in austenite fraction during intercritical annealing initially occurs in association with Mn partitioning followed by austenite growth induced by gradual cementite dissolution.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Effective conductivity of hard-sphere dispersions

A simulation technique is used to compute zeta(2), the key three-point microstructural parameter which arises in three-point conductivity bounds on the effective conductivity sigma(e) of isotropic two-phase composites. Using the new algorithm, the zeta(2) parameter is computed for dispersions of hard spheres by Monte Carlo simulation, and the results are compared with those of several analytical calculations of zeta(2). The simulation data on zeta(2) are then used to compute three-point bounds on the conductivity sigma(e) and Torquato's (1985) approximation for sigma(e).

Miller, C. A.↗
DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Elucidating the effect of preheating temperature on melt pool morphology variation in Inconel 718 laser powder bed fusion via simulation and experiment

In laser powder bed fusion (L-PBF) additive manufacturing, the mechanical performance, microstructure and defects of fabricated parts are closely associated with the melt pool morphology, e.g., its dimension and shape through the building process. Past studies have largely focused on how the process parameters such as laser power and scan speed affect melt pool characteristics. In this study, the melt pool morphology variation as a function of preheating temperature in the conduction, transition, and keyhole regimes and the underlying mechanisms in each regime are investigated through ex-situ sample characterization and computation thermal fluid dynamics (CtFD) simulation. Single tracks with different combinations of laser power and scan speed are deposited on an Inconel 718 bare plate preheated to a temperature range of 100–500 °C in the experiment. Significant changes are observed in melt pool morphology as a function of preheating temperature from optical measurements of melt track cross sections. The depth of melt pool in the three regimes increases monotonically with preheating temperature, e.g., at 500 °C, the experimental melt pool depth is increased by 49% in conduction regime, 34% in transition regime and 33% in keyhole regime, respectively, while the variation of melt pool width in each regime does not all follow an increasing trend but depends on the melt pool regimes. Melt pool width variation in the conduction and transition regimes is found to depend on the enhanced heat conduction directly related to temperature dependent thermal properties. Through validated CtFD simulations, it is found that in the keyhole regime the evaporation mass, recoil pressure, and laser drilling effect is enhanced with higher preheating temperature, which gives rise to a deeper melt pool. The simulations also reveal that preheating temperature significantly elongates the melt track length due to the increased flow rate and strong recoil pressure that accelerates the backward flow.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Long term oxidation of NiCoCrAlY coated Ni-based superalloys: A comparison of observed and simulated interdiffusion

Lifetimes of MCrAlY-type coatings can easily surpass 25 kh when the criteria of ß-depletion is used. To reduce experimental effort, it’s necessary to develop models capable of simulating interdiffusion and oxidation. In the present study, NiCoCrAlY coatings were high velocity oxyfuel deposited on PWA 1483, MARM247, and CMSX-4 substrates, and samples were exposed at 900 °C for 5–20 kh in air+10 vol%H 2 O to study interdiffusion rates. Here, Thermo-Calc/DICTRA was used to reproduce the observed differences, and microstructures were compared. Based on the agreement between calculated and observed rates, simulations were projected past the experimental exposures, utilizing ß-depletion as lifetime criteria.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupling of Spark Plasma Sintering with Advanced Modeling to Enable Process Scale-Up: Presentation to DOE-NE [Slides]

The research goal of this project is to develop at Idaho National Laboratory (INL) a first-of-its-kind Multiphysics Object-Oriented Simulation Environment (MOOSE)-based, multiscale, multiphysics spark plasma sintering (SPS) modeling and simulation code application, termed “Freya.” Freya will simulate the thermo-mechanical-electrical aspects of the SPS fabrication process and will be paired with lower length scale sub-models, such as phase-field, to predict the resulting microstructure. SPS is an advanced manufacturing process that can be used to solve a variety of material manufacturing challenges; however, this process is an extremely challenging problem for modeling and simulation. The SPS process is inherently multiphysics and multi-scale, with the macroscale electro-thermo-mechanical behavior linked intricately to the microstructure evolution of the part being sintered. Accurate modeling and simulation tools, specifically geared towards the SPS process, are needed to predict the influence of the multiple variables involved in the manufacturing process. Modeling and simulation accuracy is achieved and demonstrated through comparison to multiple validation experiments. The validation efforts for Freya include both separate effects and complete multiphysics SPS process experiments. One of the key benefits this Laboratory Directed Research & Development (LDRD) project offers stems from the emphasis placed on experimental validation of the Freya models, both on the individual length scales and of the final coupled multiscale multiphysics simulations. Experimental validation of Freya’s multiscale coupling capability provides the technical credibility necessary for potential future industry and research partners to accept the simulation predictions.

36 MATERIALS SCIENCE↗
NASA NTRS1986

A study of spectrum fatigue crack propagation in two aluminum alloys. I - Spectrum simplification. II - Influence of microstructures

An investigation of the fatigue crack propagation FCP behavior of two aluminum alloys is performed to simulate spectrum loading conditions found at critical locations in high performance fighter aircraft. Negative loads are shown to be eliminated for the tension-compression spectrum for low to intermediate maximum stress intensities, and load interactions are found to be more significant at higher stress intensities and with more plasticity at the crack tip. In the second part, the influence of microstructural features including grain size, inclusions, and dispersoids on constant amplitude and spectrum crack growth behavior in aluminum alloys is studied. At low stress intensities the I/M alloy demonstrated better FCP resistance than the P/M 7091 alloy for both constant amplitude and spectrum testing, and the inhomogeneous planar slip and large grain size of 7050 limit dislocation interactions, thereby improving FCP performance.

Telesman, J.↗
DOE OSTI2023

An asynchronous parallel high-throughput model calibration framework for crystal plasticity finite element constitutive models

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. Here in this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

304L stainless steel↗
DOE OSTI2021

Role of scan strategies and heat treatment on grain structure evolution in Fe-Si soft magnetic alloys made by laser-powder bed fusion

A major goal in printing soft magnetic Fe-Si steels using additive manufacturing is to take advantage of the potential for complex geometric designs and site-specific grain control. One major step in the processing of these alloys is understanding how processing parameters might impact how the as-built microstructure responds to annealing (i.e. the annealing response). The impact of scan strategy on the annealing response for thin wall geometries is specifically explored. Two scan strategies were explored for a thin wall geometry that produced a strongly columnar grain structure and equiaxed grain structure. Additionally, samples from both scan strategies annealed at 1200 °C showed a marked difference in annealing response with the more equiaxed sample seeing full recrystallization and grain growth, while the more columnar grain structure saw little change in microstructure. After analysis through characterization techniques and thermal-mechanical simulations Differences in internal energy within the grains were ruled out because calculated GND density values were similar for both samples. The formation of secondary particles was ruled out as a contributing factor due to the type of oxide formations and their size. It was concluded that the contributing factor to the difference in the annealing response were a difference in the resulting grain size and the density of high angle grain boundaries. These two differences were largely attributed to differences in the thermal gradient conditions due to grains preferentially growing in the direction of the steepest thermal gradient.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗
DOE OSTI2020

Characterization of PyC/SiC interfaces with FIB-SEM tomography

Tristructural-isotropic (TRISO) fuel is being considered for multiple reactor designs due to its exceptional performance. This work studies pyrocarbon (PyC) and silicon carbide (SiC) substrates that are near-representative of the inner pyrocarbon (IPyC) and SiC in TRISO fuel. Focused ion beam (FIB)-scanning electron microscopy (SEM) tomography was used to understand the influence of different processing routes on the PyC/SiC interface structure. Additionally, this work outlines a methodology for characterizing the IPyC/SiC interface when considering fission product accommodation and retention in TRISO fuel design. The high-fidelity images collected with FIB-SEM tomography can be used to achieve a more in-depth interpretation of microstructural parameters and pore-shape analysis that can inform simulations and processing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗