Determining the representative volume element size for three-dimensional microstructural material characterization. Part 1: Predictive models
Explore the source record for details and available documents.
Harris, William M.↗
SEARCH · Search NASA
Results for “microstructure model”
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.
In situ measurement of damage evolution in shocked magnesium as a function of microstructure
Accurate modeling and prediction of damage induced by dynamic loading in materials have long proved to be a difficult task. Examination of postmortem recovered samples cannot capture the time-dependent evolution of void nucleation and growth, and attempts at analytical models are hindered by the necessity to make simplifying assumptions, because of the lack of high-resolution, in situ, time-resolved experimental data. We use absorption contrast imaging to directly image the time evolution of spall damage in metals at ~1.6-μm spatial resolution. We observe a dependence of void distribution and size on time and microstructure. The insights gained from these data can be used to validate and improve dynamic damage prediction models, which have the potential to lead to the design of superior damage-resistant materials.
36 MATERIALS SCIENCE↗
Constrained GAN-Generated X-Ray CT Data For Self-Supervised And Foundation-Model Segmentation Of Concrete Microstructures
Three-dimensional characterization of materials using X-ray computed tomography (XCT) is challenging due to the complexity of internal structures, noise, and variations in resolution. Traditional computer vision models often struggle to accurately segment these images, particularly in domain-specific applications like materials science. While supervised deep learning approaches have been developed to address the limitations of conventional algorithms, they typically require large amounts of labeled training data and often fail to generalize across different datasets. Self-supervised, few-and zero-shot learning methods have gained prominence in natural image processing and segmentation tasks, but their application to scientific imaging remains limited due to the unique structural complexity, noise, and textural artifacts present in materials science data. In this work, we investigate how domain adaptation, leveraging physics-based and GAN-generated synthetic data, impacts segmentation performance. We introduce a modified Contrastive Unpaired Translation (CUT) model designed to generate realistic labeled data, which can be used for training, pre-training, and fine-tuning segmentation models for real XCT microstructure data. We evaluate the performance of two segmentation approaches: a self-supervised network (SSL-ALPNet) and a foundation model (Segment Anything Model), assessing their improvements when pre-trained and/or fine-tuned on the synthesized data. Our results demonstrate that leveraging synthetic data significantly enhances segmentation performance, particularly in challenging materials science applications.
Ziabari, Amir [ORNL] (ORCID:000000034776457X)↗
Direct numerical simulations in solid mechanics for quantifying the macroscale effects of microstructure and material model-form error
Two fundamental approximations in macroscale solid-mechanics modeling are (1) the assumption of scale separation in homogenization theory and (2) the use of a macroscopic plasticity material model that represents, in a mean sense, the multitude of inelastic processes occurring at the microscale. With the goal of quantifying the errors induced by these approximations on engineering quantities of interest, we perform a set of direct numerical simulations (DNS) in which polycrystalline microstructures are embedded throughout a macroscale structure. The largest simulations model over 50,000 grains. The microstructure is idealized using a randomly close-packed Voronoi tessellation in which each polyhedral Voronoi cell represents a grain. An face centered cubic crystal-plasticity model is used to model the mechanical response of each grain. The overall grain structure is equiaxed, and each grain is randomly oriented with no overall texture. The detailed results from the DNS simulations are compared to results obtained from conventional macroscale simulations that use homogeneous isotropic plasticity models. The macroscale plasticity models are calibrated using a representative volume element of the idealized microstructure. Furthermore, we envision that DNS modeling will be used to gain new insights into the mechanics of material deformation and failure.
36 MATERIALS SCIENCE↗
Elasto-viscoplastic fast Fourier transform modeling framework for assessing microstructural effects on stress intensity factors characterizing fracture toughness
A large-strain elasto-viscoplastic fast Fourier transform (LS-EVPFFT) model with non-periodic (NP) velocity-based boundary conditions is adapted to simulate the sensitivity of stress intensity factors on microstructure for 304L stainless steel. The material was characterized via electron backscattered diffraction (EBSD) serial-sectioning to obtain a measured 3-D microstructural cell to perform simulations. The NP-LS-EVPFFT model, including the simulation setup and boundary conditions, was verified using a crystal plasticity finite element (CPFE) model. To this end, the generation of meshes of notched specimens was developed, which involved creating Python scripts for mesh “cutting” in Abaqus, and Sculpt scripts in Cubit for meshing of the measured microstructural cell processed with DREAM.3D. The complexity of the mesh preparation highlighted the advantages of the FFT-based model, which circumvents the mesh generation process. Given the efficiency of the FFT-based model, statistical distribution of stress intensity factors in function of crystal orientation at the crack tip, grain structure, and crystallographic texture surrounding the crack tip were predicted. Further, the distributions reveal about 10% variation of stress intensity factors with microstructure with the most significant sensitivity found to be the crystal orientation at the crack tip. The methodology developed in this work is discussed as a practical simulation tool for predicting the sensitivity of stress intensity factors on microstructural variability in metallic materials.
36 MATERIALS SCIENCE↗
Emulation of neutron-irradiated microstructure of austenitic 21Cr32Ni model alloy using dual-ion irradiation
Here, in this study, the capability of heavy-ion irradiation to emulate neutron irradiation was demonstrated on an austenitic 21Cr32Ni type ternary model alloy. The model alloy used in this study is chemically analogous but compositionally simpler than of alloy 800H, which is a candidate austenitic Ni alloy which has been proposed for use in Generation IV reactors. The microstructure of the 21Cr32Ni model alloy irradiated in the BOR-60 fast reactor to 17.1 dpa and 35 dpa at ~380°C was characterized using transmission electron microscopy (TEM). The 17.1 dpa BOR-60 irradiated microstructure was then compared with the microstructure of the same material developed under dual-ion (DI) irradiation using various He/dpa ratios between 0.1 and 16.6 appm He/dpa in the temperature range of 430°C-500°C. The results showed that both neutron and DI irradiation of 21Cr32Ni model alloy produced dislocations in the form of a dislocation network as well as {111}-type faulted dislocation loops, cavities, and radiation-induced Ni enrichment at radiation-induced sinks. When the dose and the He/dpa ratio were kept similar to those in neutron irradiation, DI irradiation of the 21Cr32Ni model alloy at 460°C resulted in over-nucleation of small cavities and in a high density of faulted dislocation loops compared to those observed in the fast-neutron irradiated alloy of the same heat irradiated at ~380°C. The optimal condition for reproducing the neutron-irradiated microstructure was DI irradiation at 460°C and 0.1 appm He/dpa. In that case, the faulted loop and cavity size distributions in the BOR-60 irradiated 21Cr32Ni model alloy samples closely matched with those measured in the DI irradiated 21Cr32Ni model alloy sample. The fact that the He/dpa is an order of magnitude smaller than the helium generation rate for fast neutron irradiation, stops over nucleation and allows for the development of a similar microstructure as for neutron irradiation.
21Cr32Ni model alloy↗
A new smoothed particle hydrodynamics non-Newtonian model for friction stir welding: Process modeling and simulation of microstructure evolution in a magnesium alloy
We present a new smoothed particle hydrodynamics (SPH) model for friction stir welding (FSW). FSW has broad commercial application in the marine, aerospace, rail, and automotive industries. However, development of the FSW process for each new application has remained largely empirical. Few established numerical modeling techniques have been developed that can explain and predict important features of the process physics involved in FSW. This is particularly true in the areas of material flow and mixing mechanisms. In this paper, we present a novel modeling approach to simulate FSW that may have significant advantages over current finite element or finite difference based methods. Unlike traditional grid-based methods, Lagrangian particle methods such as SPH can simulate the dynamics of interfaces, large material deformations, and the material’s strain and temperature history without employing complex tracking schemes. Three-dimensional simulations of FSW on AZ31 Mg alloy are performed. The temperature history and distribution, grain size, microhardness as well as the texture evolution are presented. Numerical results are found to be in good agreement with experimental observations.
Pan, Wenxiao↗
Microstructural evolution in NF616 (P92) and Fe–9Cr–0.1C-model alloy under heavy ion irradiation
In this comparative study, in situ investigations of the microstructure evolution in a Fee9Cr ferritic emartensitic steel, NF616, and a Fee9Cre0.1C-model alloy with a similar ferriticemartensitic microstructure have been performed. NF616 and Fee9Cre0.1C-model alloy were irradiated to high doses (up to ~10 dpa) with 1 MeV Kr ions between 50 and 673 K. Defect cluster density increased with dose and saturated in both alloys. The average size of defect clusters in NF616 was constant between 50 and 573 K, on the other hand average defect size increased with dose in Fee9Cre0.1C-model alloy around ~1 dpa. At low temperatures (50e298 K), alignment of small defect clusters resulted in the formation of extensive defects in Fee9Cre0.1C-model alloy around ~2e3 dpa, while similar large defects in NF616 started to form at a high temperature of 673 K around ~5 dpa. Interaction of defect clusters with the lath boundaries were found to be much more noticeable in Fee9Cre0.1C-model alloy. Differences in the microstructural evolution of NF616 and Fee9Cre0.1C-model alloy are explained by means of the defect cluster trapping by solute atoms which depends on the solute atom concentrations in the alloys.
Ferritic-martensitic alloys↗
Novel method to control explosive shock sensitivity: Formulations, experiments and modeling
The microstructure of a heterogeneous high explosive (HE) affects its shock initiation sensitivity and detonation performance. Alteration of the void content and/or void structure (i.e., bulk heating or mechanical damage) therefore changes the shock initiation behavior. Controlling and predicting the change in shock sensitivity after an HE has undergone microstructural changes addresses an important and challenging goal for the design and understanding of novel energetic material formulations. In this study, we aimed to develop HE systems to precisely tune shock sensitivity by a thermal treatment prior to use. Specifically, we incorporated a small fraction (1 wt% or less) of thermally expandable microspheres (TEMs) during the formulation process of various plastic-bonded explosives (PBX). TEMs typically consist of a thermoplastic acrylonitrile shell (10-50 µm diameter), which encapsulates an inert low boiling hydrocarbon. Upon heating, the TEMs expand as the shell softens while the hydrocarbon gasifies, increasing the internal pressure and expanding the particle by as much as 120 vol% (irreversibly). Here, we present our progress on shock sensitivity comparisons of HE formulations doped with TEMs after heating and expansion. Experiments, using Composition C-4 thermally cycled to 120°C, showed an increase in shock sensitivity. Additionally, mesoscale modeling revealed that the TEM itself does not act as a hotspot, but instead has a secondary effect on the response of nearby voids. We conclude that TEMs indeed do provide a method to tune shock sensitivity. Further, our results emphasize the non-uniformity of shock waves at the mesoscale, and that upstream defects influence the reactivity of downstream defects.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Microstructure Scale Lithium-Ion Battery Modeling: Part III. When and Where Lithium Plating Occurs and its Correlation with the Electrode Microstructure
Li-ion battery performance and degradation are closely related to the cell’s underlying electrode microstructure. Electrode microstructures are typically characterized with volume-averaged properties that neglect the impact of local heterogeneities. However, local heterogeneities create hot spots that can trigger degradation onset. Herein, a microstructure scale electrochemical model is used to investigate the impact of microstructure heterogeneity on lithium plating. The model predicts lithium plating is not uniform, even when considering a relatively small portion of the electrode (a cross-sectional area of 154×144 µm 2 ), preferring to plate on larger particles as compared to smaller particles. While local heterogeneities control where plating occurs, the model predicts that volume-averaged properties control when plating occurs. Additionally, the model predicts that the active material specific surface area has a linear relationship with the plating onset. However, the linear relationship between increased active material surface area and delayed plating response appears to be sensitive to the microstructure feature used to increase the active interface area. Here, a comparative case-study is explored where the specific surface area is increased by either reducing the active material particle diameter, adding open-porosity cracks, or increasing the active material surface roughness. The model predicts that increasing the specific surface area by reducing the active material particle diameter is the most effective strategy for delaying lithium plating. At 6C, reducing particle size is shown to be 3 and 20 times more effective than, respectively, adding open-porosity cracks and increasing surface roughness. A dual-layer electrode architecture combining gradations both for average properties and uniformities is eventually proposed to improve homogeneous material utilization and reduce degradation at high charge rates.
25 ENERGY STORAGE↗
Creating physically-based three-dimensional microstructures: Bridging phase-field and crystal plasticity models.
In order to better incorporate microstructures in continuum scale models, we use a novel finite element (FE) meshing technique to generate three-dimensional polycrystalline aggregates from a phase field grain growth model of grain microstructures. The proposed meshing technique creates hexahedral FE meshes that capture smooth interfaces between adjacent grains. Three dimensional realizations of grain microstructures from the phase field model are used in crystal plasticity-finite element (CP-FE) simulations of polycrystalline a -iron. We show that the interface conformal meshes significantly reduce artificial stress localizations in voxelated meshes that exhibit the so-called "wedding cake" interfaces. This framework provides a direct link between two mesoscale models - phase field and crystal plasticity - and for the first time allows mechanics simulations of polycrystalline materials using three-dimensional hexahedral finite element meshes with realistic topological features.
36 MATERIALS SCIENCE↗
Mesoscale Modeling Approach for Quantifying Microstructure-Aware Micromechanical Responses in Metal Hydrides
Metal hydrides can undergo significant volume changes upon hydrogen uptake and release, which induce a mechanical response that depends not only on the evolving hydrogen composition but also on the microstructure. We present a comprehensive mesoscale modeling framework based on microelasticity theory to quantify the micromechanical responses of metal hydrides, specifically focusing on a hydrogenating polycrystalline MgH 2x particle within a host material as a model micromechanical system. Utilizing digitally generated realistic microstructures and density-functional-theory-derived parameters, we analyzed highly nonuniform local stress profiles in the polycrystalline hydrides under the clamping force exerted by the host during hydrogenation. Our framework also allows us to predict the corresponding strain energy accumulation and mechanical hot spots formation in the hydrides, highlighting their roles in thermodynamic destabilization and mechanical failure, respectively. Through extensive parametric simulations, we further quantified the influence of interface type, crystallinity, grain size, loading ratio, and host stiffness, providing practical guidance for optimizing microstructural design and host material selection. This proposed approach is broadly applicable to micromechanical systems with complex microstructural features involving chemical reaction- and/or phase-transformation-induced deformation.
36 MATERIALS SCIENCE↗
Graphite Oxidation Behavior: NRC Graphite Behavior Model
Graphite microstructure and oxidation Graphite crystal structure and basal planes Reactive surface area sites Zig-zag and Arm-chair Pore microstructure effects Interconnectivity Reactivity and oxygen transport Oxygen-graphitereaction Not C-O reaction! Temperature effects Oxygen diffusion to graphite Diffusion into the microstructure High temperature vs low temperature diffusion Oxidation: It’s a gradient of factors
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Preliminary Phase Field Computational Model Development
This interim report presents progress towards the development of meso-scale models of magnetic behavior that incorporate microstructural information. Modeling magnetic signatures in irradiated materials with complex microstructures (such as structural steels) is a significant challenge. The complexity is addressed incrementally, using the monocrystalline Fe (i.e., ferrite) film as model systems to develop and validate initial models, followed by polycrystalline Fe films, and by more complicated and representative alloys. In addition, the modeling incrementally addresses inclusion of other major phases (e.g., martensite, austenite), minor magnetic phases (e.g., carbides, FeCr precipitates), and minor nonmagnetic phases (e.g., Cu precipitates, voids). The focus of the magnetic modeling is on phase-field models. The models are based on the numerical solution to the Landau-Lifshitz-Gilbert equation. From the computational standpoint, phase-field modeling allows the simulation of large enough systems that relevant defect structures and their effects on functional properties like magnetism can be simulated. To date, two phase-field models have been generated in support of this work. First, a bulk iron model with periodic boundary conditions was generated as a proof-of-concept to investigate major loop effects of single versus polycrystalline bulk iron and effects of single non-magnetic defects. More recently, to support the experimental program herein using iron thin films, a new model was generated that uses finite boundary conditions representing surfaces and edges. This model has provided key insights into the domain structures observed in magnetic force microscopy (MFM) measurements. Simulation results for single crystal thin-film iron indicate the feasibility of the model for determining magnetic domain wall thickness and mobility in an externally applied field. Because the phase-field model dimensions are limited relative to the size of most specimens used in experiments, special experimental methods were devised to create similar boundary conditions in the iron films. Preliminary MFM studies conducted on single and polycrystalline iron films with small sub-areas created with focused ion beam have correlated quite well qualitatively with phase-field simulations. However, phase-field model dimensions are still small relative to experiments thus far. We are in the process of increasing the size of the models and decreasing specimen size so both have identical dimensions. Ongoing research is focused on validation of the phase-field model. Validation is being accomplished through comparison with experimentally obtained MFM images (in progress), and planned measurements of major hysteresis loops and first order reversal curves. Extrapolation of simulation sizes to represent a more stochastic bulk-like system will require sampling of various simulations (i.e., with single non-magnetic defect, single magnetic defect, single grain boundary, single dislocation, etc.) with distributions of input parameters. These outputs can then be compared to laboratory magnetic measurements and ultimately to simulate magnetic Barkhausen noise signals.
36 MATERIALS SCIENCE↗
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↗
An integrated two-dimensional modeling method for predicting ductility of thin-walled die cast magnesium
In this paper, a two-dimensional (2D) microstructure-based modeling method is developed in order to predict the ductility of a thin-walled high pressure die cast magnesium (Mg) by considering the three-dimensional (3D) thru-thickness pore distributions. For this purpose, a series of 3D synthetic microstructure-based finite element models and the corresponding 2D models are first generated with various pore volume fractions, pore size distributions and pore shapes. The input material properties for the 2D models are determined based on the 3D cubic model with a spherical pore and the generalized Neuber’s rule. Based on the resulting ductility of the 3D and 2D models, a 3D/2D ductility correlation curve is developed as a function of the characteristics/shape of the input fracture strain curve used in the 3D model. The validity of the ductility correlation curve is examined with casting samples with actual microstructures. Here, the results show that the suggested 2D modeling methodology can be used together with the ductility correlation curve in predicting the ductility of thin-walled Mg castings.
36 MATERIALS SCIENCE↗
Comprehensive characterization of irradiation induced defects in ceria: Impact of point defects on vibrational and optical properties
Validation of multiscale microstructure evolution models can be improved when standard microstructure characterization tools are coupled with methods sensitive to individual point defects. We demonstrate how electronic and vibrational properties of defects revealed by optical absorption and Raman spectroscopies can be used to compliment transmission electron microscopy (TEM) and x-ray diffraction (XRD) in the characterization of microstructure evolution in ceria under non-equilibrium conditions. Experimental manifestation of non-equilibrium conditions was realized by exposing cerium dioxide (CeO2) to energetic protons at elevated temperature. Two sintered polycrystalline CeO2 samples were bombarded with protons accelerated to a few MeVs. These irradiation conditions produced a microstructure with resolvable extended defects and a significant concentration of point defects. A rate theory (RT) model was parametrized using the results of TEM, XRD, and thermal conductivity measurements to infer point defect concentrations. An abundance of cerium sublattice defects suggested by the RT model is supported by Raman spectroscopy measurements, which show peak shift and broadening of the intrinsic T2g peak and emergence of new defect peaks. Additionally, spectroscopic ellipsometry measurements performed in lieu of optical absorption reveals the presence of Ce3+ ions associated with oxygen vacancies. This work lays the foundation for a coupled approach that considers a multimodal characterization of microstructures to guide and validate complex defect evolution models.
Chauhan, Vinay S.↗