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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Microstructural characterization of the oxidation of nuclear graphite under chronic and accident conditions via XCT and synchrotron XCT

Graphite is a candidate material to moderate fast neutrons and for structural components in the US next-generation graphite-moderated reactors. A graphite core is conceived as a large formation of interconnected bricks that primarily serves as the moderator of fast neutrons, holds instrumentation, fuel elements, control rods, and is a receptacle for molten salts for Molten Salt Reactors (MSRs) designs. During the operation of a graphite-moderated power plant, graphite components might be subjected to chronic oxidation during normal operating conditions or aggressive oxidation as a result of accidental ingress of air that reacts vigorously with the graphite core. Using synchrotron, x-ray computed tomography (XCT), this research systematically characterized microstructural changes that accompanied these two oxidation scenarios. Chronic oxidation was studied by characterizing IG-110, PCEA, and NBG-18 specimens that were gradually oxidized in air at a low temperature (520°C). The accidental ingress of air into the graphite reactor was simulated by subjecting the grades mentioned above to high-temperature oxidation at approximately 750°C. This research is the first, in situ, systematic characterization of nuclear graphite microstructural evolution that can be associated with the two possible oxidation scenarios and provide insight into related repercussions. The results show that the microstructure and pore connectivity influence the rate of oxidation and evolution of the microstructure under the two oxidation regimes. These results are crucial to understanding which grades are more resilient to each type of oxidation and documenting the damage created in the graphite.

Arregui Mena, Jose'↗
DOE OSTI2026

Additively manufactured refractory high-entropy alloys with superior radiation resistance

Refractory high-entropy alloys (RHEAs) are promising candidates for next-generation nuclear and high-temperature applications. Among many approaches to manufacture RHEAs, additive manufacturing (AM) represents the most recent and advanced metal manufacturing method which allows near-net-shape manufacturing to reduce material waste and post-processing time. However, performance of AM RHEAs under complex irradiation conditions remains largely unexplored. Here, in this study, we demonstrate for the first time the response of directed energy deposition (DED) AM quaternary RHEAs (HfTaVW, CrTaVW) subjected to sequential dual-beam ion irradiation, consisting of helium pre-implantation followed by high-dose heavy ion bombardment. Compositions of DED AM RHEAs were selected using Monte Carlo (MC) simulations based on a cluster expansion (CE) Hamiltonian parameterized by density functional theory (DFT). Post-irradiation microstructural characterization revealed that the AM RHEA maintained remarkable stability, with suppressed helium bubble growth and reduced defect accumulation compared to conventional alloys. Even at high doses (∼100 dpa), the alloy exhibited no void swelling, a low density of dislocation loops, and no evidence of severe degradation. These results highlight the intrinsic ability of AM-derived microstructures and multicomponent chemistry to synergistically mitigate irradiation effects. Our findings establish AM RHEAs as a class of materials with superior resistance to radiation damage under conditions relevant to advanced fusion and fission environments and demonstrate the importance of sequential ion beam studies in evaluating their long-term performance.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct current microplasma formation around microstructure arrays

We demonstrate the formation and transition behaviors of a microplasma around microstructure arrays at different gas pressures via two-dimensional particle-in-cell/Monte Carlo collision simulations. It is found that the microdischarge occurs outside the cathode microcavities at the lowest pressure and starts penetrating the microcavities with a curved sheath edge as the pressure increases. At higher pressure, coupled periodic microhollow cathode discharges (MHCDs) are formed inside the microcavities. Further increasing the gas pressure results in the disappearance of the MHCDs, and the dominant discharge shifts outside of the microcavity, locating above the protrusion tips. The effect of the space charge shielding on the discharge and the conditions for MHCD formation are discussed. The macroscopic discharge parameter scalings with the gas pressure and the electron kinetics are also examined. The results are helpful for deeply understanding the microplasma formation with nonplanar electrodes, which inform the scaling, design, and optimization of microplasma array devices across a wide range of pressure regimes in practical applications.

42 ENGINEERING↗
DOE OSTI2022

Process-microstructure relationship of laser processed thermoelectric material Bi2Te3

Additive manufacturing allows fabrication of custom-shaped thermoelectric materials while minimizing waste, reducing processing steps, and maximizing integration compared to conventional methods. Establishing the process-structure-property relationship of laser additive manufactured thermoelectric materials facilitates enhanced process control and thermoelectric performance. This research focuses on laser processing of bismuth telluride (Bi 2 Te 3 ), a well-established thermoelectric material for low temperature applications. Single melt tracks under various parameters (laser power, scan speed and number of scans) were processed on Bi 2 Te 3 powder compacts. A detailed analysis of the transition in the melting mode, grain growth, balling formation, and elemental composition is provided. Rapid melting and solidification of Bi 2 Te 3 resulted in fine-grained microstructure with preferential grain growth along the direction of the temperature gradient. Experimental results were corroborated with simulations for melt pool dimensions as well as grain morphology transitions resulting from the relationship between temperature gradient and solidification rate. Samples processed at 25 W, 350 mm/s with 5 scans resulted in minimized balling and porosity, along with columnar grains having a high density of dislocations.

Oztan, Cagri↗
NASA NTRSDate not supplied

Cathode Modeling of Solid-State Batteries

The search for safe, reliable, and compact high-capacity energy storage devices has led to increased interest in all-solid-state battery research. The use of solid electrolytes provides enhanced safety and durability due to their reduced flammability and increased mechanical strength compared to organic liquid electrolytes. Still, the use of solid electrolytes remains challenging. Computational modeling plays a substantial role in addressing these challenges. A particle dynamics electromechanical model for simulating electrochemical processes in a solid-state battery cathode will be presented. The model presents cathode microstructure at the particle level as a mixture of ionically conductive solid electrolyte particles, electrically conductive carbon additives, and cathodic reactant particles. After densification, the particle connectivity is analyzed to reconstruct the complex electric network connecting reactant particles with an anodic and cathodic current collectors through the electrolyte and carbon particles. The Kirchhoff’s matrix equation describing this electric network, is solved to obtain values of various critical parameters, such as the overall conductivity of the cathode for lithium ions and electrons, cathodic reactant material utilization, and the distribution of the electric current and voltages within the cathode. In addition, by representing the reactant particles as electrolyte or galvanic microcells governed by the Butler-Volmer electrochemical equation, the overall performance of battery cells during charge or discharge processes, respectively, can be predicted for a given cathodic powder composition. The presented model, executed on a high-performance computing architecture, essentially provides a valuable guidance in designing and developing future solid-state batteries.

solid-state battery↗
DOE OSTI2024

Image-Based Failure Assessment of Li-Ion Batteries

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. In this work, a digital twin is created to investigate the performance of a heterogeneous Li-ion battery cathode and simulate degradation accumulation. Pixel-based model construction is used to represent the complex material geometries from microstructural images supplied by the National Renewable Energy Laboratory (NREL). Because of the expected large deformation and crack opening, the reproducing kernel particle method (RKPM), a meshfree method with discretization at the image pixels, is used to approximate the field variables: electrostatic potential, concentration, and displacement. An interface modified reproducing kernel (IM-RK) is 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. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. IM-RK is additionally used to inform how crack evolution in turn affects the coupled electro-chemo-mechanical behavior of the Li-ion battery cathode.

image-based modeling↗
DOE OSTI2023

Investigation of the Effect of Chrome and Nickel Concentrations During Two-Piston Splat Quenching of Austenitic Stainless Steels

For solidification rates at or near equilibrium solidification conditions the effects of chrome (Cr) and nickel (Ni) on stainless steel solidification modes and microstructures are well detailed. However, fusion-based additive manufacturing (FBAM) processes that rely on faster, more rapid solidification rates call for a more in-depth understanding of the effects of rapid solidification on the solidification behavior and how these change with variations in alloy composition. Eleven custom stainless steel (SS) alloys with unique compositions based on 316L were made using targeted alloying element additions to generate feedstock with Cr/Ni eq ratios ranging from 0.9 to 2.1. Two-piston splat quenching (SQ) was used to produce rapid solidification conditions similar to those achieved in powder bed fusion processes in a fraction of the time and cost of traditional methods. Employing heat transfer simulations, SEM, TEM, and STEM characterization techniques, the SQ process was found to consistently produce solidification rates estimated to be between ~ 0.4 and 1.6 m/s. Five unique solidification microstructures were identified within the rapidly solidified SQ samples (primary austenite, massively transformed ferrite to austenite, primary austenite and massively transformed ferrite to austenite, primary ferrite, and primary ferrite and massively transformed ferrite to austenite). Both the primary ferrite solidification mode and primary ferrite phase were found to form at lower Cr/Ni eq values than previously predicted. SQ samples with compositions that were within the compositional specifications for 316L SS, an alloy which is commonly used in FBAM processes, demonstrated a wider range of potential solidification modes and microstructures that could form at rapid solidification rates than expected. Finally, by using SQ as a means to simulate rapid solidification conditions similar to those observed in FBAM processes, potential new alloy compositions can be screened faster and more cost effectively than purchasing and running test batches of the metallic powder.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Transmural Remodeling of Cardiac Microstructure in Aged Spontaneously Hypertensive Rats by Diffusion Tensor MRI

The long-standing high blood pressure (also known as hypertension) overworks the heart. Microstructural remodeling is a key factor of hypertensive heart disease progression. Diffusion tensor magnetic resonance imaging (DT-MRI) is a powerful tool for the rapid noninvasive nondestructive delineation of the cardiomyocyte organization. The spontaneously hypertensive rat (SHR) is a well-established model of genetic hypertension. The goal of this study was to employ high-resolution DT-MRI and the SHR animal model to assess the transmural layer-specific remodeling of myocardial microstructure associated with hypertension. Ex vivo experiments were performed on excised formalin-fixed hearts of aged SHRs (n = 4) and age-matched controls (n = 4). The DT-MRI-derived fractional anisotropy (FA), longitudinal diffusivity (λ L ), transversal diffusivity (λ T ), and mean diffusivity (MD) served as the readout parameters investigated at three transmural zones (i.e., endocardium, mesocardium, and epicardium). The helix angles (HAs) of the aggregated cardiomyocytes and the orientation of laminar sheetlets were also studied. Compared with controls, the SHRs exhibited decreased epicardial FA, while FA changes in the other two transmural regions were insignificant. No substantial differences were observed in the diffusivity parameters and the transmural course of HAs between the two groups. A consistent distribution pattern of laminar sheetlet orientation was not identified for either group. Our findings are in line with the known cellular microstructure from early painstaking histological studies. Biophysical explanations of the study outcomes are provided. In conclusion, our experimental findings indicate that the epicardial microstructure is more vulnerable to high blood pressure leading to more pronounced changes in this region during remodeling. DT-MRI is well-suited for elucidating these alterations. The revealed transmural nonuniformity of myocardial reorganization may shed light on the mechanisms of the microstructure-function relationship in hypertension progression. Our results provide insights into the management of patients with systemic arterial hypertension, thus prevent the progression toward heart failure. The findings of this study should be acknowledged by electromechanical models of the heart that simulate the specific cardiac pathology.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

In Situ X-ray Radiography and Computational Modeling to Predict Grain Morphology in $\beta$-Titanium during Simulated Additive Manufacturing

The continued development of metal additive manufacturing (AM) has expanded the engineering metallic alloys for which these processes may be applied, including beta-titanium alloys with desirable strength-to-density ratios. To understand the response of beta-titanium alloys to AM processing, solidification and microstructure evolution needs to be investigated. In particular, thermal gradients (Gs) and solidification velocities (Vs) experienced during AM are needed to link processing to microstructure development, including the columnar-to-equiaxed transition (CET). In this work, in situ synchrotron X-ray radiography of the beta-titanium alloy Ti-10V-2Fe-3Al (wt.%) (Ti-1023) during simulated laser-powder bed fusion (L-PBF) was performed at the Advanced Photon Source at Argonne National Laboratory, allowing for direct determination of Vs. Two different computational modeling tools, SYSWELD and FLOW-3D, were utilized to investigate the solidification conditions of spot and raster melt scenarios. The predicted Vs obtained from both pieces of computational software exhibited good agreement with those obtained from in situ synchrotron X-ray radiography measurements. The model that accounted for fluid flow also showed the ability to predict trends unobservable in the in situ synchrotron X-ray radiography, but are known to occur during rapid solidification. A CET model for Ti-1023 was also developed using the Kurz–Giovanola–Trivedi model, which allowed modeled Gs and Vs to be compared in the context of predicted grain morphologies. Both pieces of software were in agreement for morphology predictions of spot-melts, but drastically differed for raster predictions. The discrepancy is attributable to the difference in accounting for fluid flow, resulting in magnitude-different values of Gs for similar Vs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Neural network–based pore flow field prediction in porous media using super resolution

Direct pore-scale simulations of fluid flow through porous media are computationally expensive to perform for realistic systems. Previous works have demonstrated using the geometry of the microstructure of porous media to predict the velocity fields therein based on neural networks. However, such trained neural networks do not perform well for unseen porous media with a large degree of heterogeneity. In this study we propose that incorporating a coarse velocity field in the input of neural networks is an effective way to improve the prediction performance. The coarse velocity field can be simulated with a low computational cost and provides global information to regularize the ill-posedness of the learning problem, which is usually caused by the use of local geometries due to the computational resource constraints. We show that incorporating the coarse-mesh velocity field significantly improves the prediction accuracy of the fine-mesh velocity field by comparison to the prediction that relies on geometric information alone, especially for the porous medium with a large interior vuggy pore space. We also show the flexibility of training the network in using coarse velocity fields with various resolutions. The results suggest that even using coarse velocity field with a very low resolution, the predictions are still enhanced and close to the ground truths. The feasibility of the method is further demonstrated by testing the trained network on real rocks. This study highlights the merits of incorporating a coarse-mesh velocity field into the input for neural networks, which provides global, physics-based information for the model, thereby improving the model's generalization capability.

42 ENGINEERING↗
DOE OSTI2023

A combined experimental and crystal plasticity study of grain size effects in magnesium alloys

This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.

Crystal plasticity↗
DOE OSTI2024

Microstructure and Chemical Composition of Al2O3-, Cr2O3- and MgO-rich Refractories Exposed to Plastic Ashes

Plastics are widely integrated with today’s lifestyle due to their low cost, durability, lightweight, water resistance, and flexibility. However, tremendous plastic waste cause severe threats to the environment. Thus, sustainable recycling of plastic waste is critical, where plastic waste gasification technology with carbon management becomes one of the most effective recycling methods. We simulate the gasification process by exposing the Al2O3-, Cr2O3- and MgO-rich refractories to synthetic plastic ashes at 1500 ℃ for 50 h. A comprehensive study on the microstructure and chemical composition change before/after exposure is conducted using scanning electron microscopy. Meanwhile, the differences among these three refractories after exposure are revealed for quantifying the performance in plastic ashes. Discussion of the results will establish an understanding of the influence of plastic ashes on the microstructure, composition, and phase of different types of refractories, providing insights into designing new refractory materials for future recycling of plastic waste.

Fang, Xiaotian↗
DOE OSTI2020

Frontiers in the Simulation of Dislocations

Dislocations play a vital role in the mechanical behavior of crystalline materials during deformation. To capture dislocation phenomena across all relevant scales, a multiscale modeling framework of plasticity has emerged, with the goal of reaching a quantitative understanding of microstructure–property relations, for instance, to predict the strength and toughness of metals and alloys for engineering applications. This review describes the state of the art of the major dislocation modeling techniques, and then discusses how recent progress can be leveraged to advance the frontiers in simulations of dislocations. Furthermore, the frontiers of dislocation modeling include opportunities to establish quantitative connections between the scales, validate models against experiments, and use data science methods (e.g., machine learning) to gain an understanding of and enhance the current predictive capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Volumetric analysis and mesh generation of real and artificial microstructural geometries

Producing a viable finite element mesh of realistic microstructural structural geometry is a critical step in analyzing the thermo-mechanical behavior of complex multi-material composites. Advancements in imaging technology such as micro computed tomography have allowed modelers to access high resolution mesoscale geometries for direct numerical simulation. However, converting from voxel based 3D images to usable finite element meshes has been challenging. A robust method including algorithms and software scripts for generating finite element meshes from 3D imaged microstructures is presented in this paper. It includes a routine for inserting cohesive elements around material interfaces to enable modeling of interface properties including delamination and damage. The algorithms and procedures presented in this method leverage currently available software packages for processing surface based geometry into volume based meshes. In addition to converting real geometry from physical imaging systems, algorithms for producing numerically generated and statistically equivalent microstructural geometry are also included. These artificial microstructures can be a valuable resource for modelers when physical specimens do not exist or are limited in quantity: Method establishes a workflow from voxel data to viable finite element mesh including interface information, and Includes a method for synthetic geometry generation based on metrics from real microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructure-Sensitive Uncertainty Quantification for Crystal Plasticity Finite Element Constitutive Models Using Stochastic Collocation Methods

Uncertainty quantification (UQ) plays a major role in verification and validation for computational engineering models and simulations, and establishes trust in the predictive capability of computational models. In the materials science and engineering context, where the process-structure-property-performance linkage is well known to be the only road mapping from manufacturing to engineering performance, numerous integrated computational materials engineering (ICME) models have been developed across a wide spectrum of length-scales and time-scales to relieve the burden of resource-intensive experiments. Within the structure-property linkage, crystal plasticity finite element method (CPFEM) models have been widely used since they are one of a few ICME toolboxes that allows numerical predictions, providing the bridge from microstructure to materials properties and performances. Several constitutive models have been proposed in the last few decades to capture the mechanics and plasticity behavior of materials. While some UQ studies have been performed, the robustness and uncertainty of these constitutive models have not been rigorously established. In this work, we apply a stochastic collocation (SC) method, which is mathematically rigorous and has been widely used in the field of UQ, to quantify the uncertainty of three most commonly used constitutive models in CPFEM, namely phenomenological models (with and without twinning), and dislocation-density-based constitutive models, for three different types of crystal structures, namely face-centered cubic (fcc) copper (Cu), body-centered cubic (bcc) tungsten (W), and hexagonal close packing (hcp) magnesium (Mg). Our numerical results not only quantify the uncertainty of these constitutive models in stress-strain curve, but also analyze the global sensitivity of the underlying constitutive parameters with respect to the initial yield behavior, which may be helpful for robust constitutive model calibration works in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Super-resolution model for overlapping peak detection and improved spatial resolution in high-energy diffraction microscopy

Reconstruction quality in Far-field High-Energy Diffraction Microscopy (FF-HEDM) is limited by the spatial resolution of area detectors and the frequent occurrence of overlapping diffraction spots. To address these challenges, we developed a super-resolution (SR) framework using convolutional neural networks (CNNs) to recreate 2D diffraction peaks at up to ×8 resolution from raw detector data. A specialized simulation tool was created to generate synthetic training datasets with varying degrees of peak overlap. Integrated into the Microstructural Imaging using Diffraction Analysis Software (MIDAS), the SR model improves the spatial accuracy and precision of 3D grain reconstruction by an order of magnitude. This approach provides a robust solution for investigation of complex micromechanical states and material classes where the analysis is limited by the presence of overlapping peaks. Furthermore, the methodology developed here can potentially be extended to other techniques that require sub-pixel accuracy for high-fidelity data analysis.

High-energy diffraction microscopy↗
DOE OSTI2023

A mechanistic model for creep and thermal aging in Alloy 709

This report describes a physics-based model for creep and thermal aging in Alloy 709. Alloy 709 is an advanced austenitic alloy, targeted for use in future Sodium Fast Reactors (SFRs) and other advanced reactors. The material has superior high temperature properties compared to currently qualified 316 and 304 stainless steels. However, the available creep and thermal aging test database for Alloy 709 is significantly more limited compared to the historical materials. The physics-based model developed here is one way to accelerate the qualification of the material by providing more accurate long-term predictions for creep properties and thermal aging, compared to current empirical time-extrapolate techniques. The crystal plasticity finite element model is used to predict the deformation and failure of alloy 709. The same setup for the CPFE model is used in both the baseline model calibration process and the simulation campaigns for parameter inference. Specific constitutive choices are made for Alloy 709 to capture the primary deformation mechanisms. The dislocation creep formulation developed by Hu and Cocks is extended to account for coupled precipitation formation and the grain boundary cavitation model developed by Sham, Needleman, et al. is used to model grain boundary cavitation-induced failure. A novel update algorithm is proposed to render the semi-discrete constitutive update for the Sham-Needleman model unconditionally stable. A progressive calibration approach is adopted based on the observations that several types of material responses can be effectively decoupled. A surrogate model is trained based on full-fledged CPFE simulations to accelerate the forward model evaluations, and stochastic variational inference (SVI) is used to calibrate the unknown microstructural model parameters. The calibrated mechanistic model is used to predict the long-term creep life of Alloy 709, and the predictions are compared against classical empirical approaches.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Research study on materials processing in space, M566 experiment

Specimens of the aluminum-33 wt% copper eutectic partially melted and resolidified in the low effective gravity of the orbiting Skylab were examined and characterized with respect to microstructural defects and thermal conductivity values. The results obtained were compared with similar evaluations of ground-based simulation melt-resolidification experiments and as-prepared unidirectionally solidified specimens. Thermal conductivity data and electrical resistivity data at temperatures from 25 C to 400 C did not show significant differences between ground and space processed specimens. A methology of evaluating the defects in the Al-Al2Cu structure was implemented. A specimen from Skylab 3 showed signs of instability in growth and several grains were found in the ingot. The specimen from Skylab 4 did not show such marked instability in growth and was found to contain fewer defects than the ground-processed specimens. This agrees with data from Georgia Institute of Technology which showed that there were fewer defects in both their Skylab 3 and 4 specimens than in ground processed specimens.

Douglas, F. C.↗