Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15

Prediction of Probabilistic Shock Initiation Thresholds of Energetic Materials Through Evolution of Thermal-Mechanical Dissipation and Reactive Heating

The ignition threshold of an energetic material (EM) quantifies the macroscopic conditions for the onset of self-sustaining chemical reactions. The threshold is an important theoretical and practical measure of material attributes that relate to safety and reliability. Historically, the thresholds are measured experimentally. In this work, we present a new Lagrangian computational framework for establishing the probabilistic ignition thresholds of heterogeneous EM out of the evolutions of coupled mechanical-thermal-chemical processes using mesoscale simulations. Furthermore, the simulations explicitly account for microstructural heterogeneities, constituent properties, and interfacial processes and capture processes responsible for the development of material damage and the formation of hotspots in which chemical reactions initiate. The specific mechanisms tracked include viscoelasticity, viscoplasticity, fracture, post-fracture contact, frictional heating, heat conduction, reactive chemical heating, gaseous product generation, and convective heat transfer. To determine the ignition threshold, the minimum macroscopic loading required to achieve self-sustaining chemical reactions with a rate of reactive heat generation exceeding the rate of heat loss due to conduction and other dissipative mechanisms is determined. Probabilistic quantification of the processes and the thresholds are obtained via the use of statistically equivalent microstructure sample sets (SEMSS). The predictions are in agreement with available experimental data.

36 MATERIALS SCIENCE↗

Micro–macro finite element modeling method for rub response in abradable coating materials

Gas turbine engines experience “rub” when the rotating blades come in contact with a static abradable coating. This results in extreme strain rates and dynamics inside a high-temperature/high-pressure environment. Current rub models are phenomenological and do not reflect the underlying microstructures, thus limiting their prediction accuracy. In this work, a microstructure-informed, reduced order modeling framework is introduced for simulating abradable coating “rub" behavior. This framework comprises a microscale model constructed based on digitized abradable microstructure and explicitly simulates the mechanical behavior of each constituent phases and their interactions. After calibration and validation with experiment data, the calibrated microscale model is used to generate data across a vast range of applied strain rates and temperature with various load paths. Then, the virtually generated data are used to fit the macroscopic-reduced order model, which enables fast component scale rub simulation without compromising the integrity of the complex material behavior. In conclusion, the proposed effort will address the technical challenge of predicting abradable material behavior during rub through the application of multiscale modeling from microstructure to engines behavior, effectively reducing the development costs and time of new abradable material for better “rub” properties.

36 MATERIALS SCIENCE↗

Microstructure and Fission Product Distribution Examination in the UCO kernel of TRISO Fuel Particles

The effect of irradiation and fabrication process on the UCO fuel kernel, kernel-buffer interlayer and fission products, has not been fully studied for UCO fuel kernels of TRISO coated particles. Specifically, the potential correlation(s) with fission product retention, composition and distribution within these structures for the tristructural isotropic (TRISO) fuel performance at normal and accident conditions, are not well reported in literature. Structural and chemical information in the irradiated UCO fuel kernel, e.g., the stoichiometry variations, may further contribute to fuel performance modeling under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The microstructural and micro-chemical evolution of selected irradiated UCO fuel kernels, kernel buffer interfaces and recoil zones were characterized using Scanning Transmission Electron Microscopy (STEM) and Energy Dispersive Spectroscopy (EDS) elemental mapping. The fuel kernel of particles from Advanced Gas Reactor (AGR)-2 irradiated to 10.8% fissions per initial metal atom (FIMA), were examined with the following main findings: (1) large fission gas bubbles were mainly identified in the zirconium and molybdenum enriched but oxygen depleted phases in the center of fuel kernel, (2) significant inter-diffusion between the UCO fuel kernel and pyro-carbon buffer layer was observed as a high density of sub-micron sized carbon inclusions was identified at near the interface in the kernel side and some isolated UCO particles were located in the buffer layer, (3) pores filled with free carbon were also observed in the recoil zone specimen. Fuel kernels from two AGR-1 compacts irradiated to 11.4 and 17.4% FIMA respectively, were examined with the following main findings: (1) two primary phases consisting of a “high-Z (atomic mass)” UC(O) phase (rocksalt; a=0.496 nm) that also contained Zr, Mo, and/or Ru and a “low-Z” UO2(C) phase (fluorite structure; a = 0.547nm) in which trace amounts of I, Nd, Pr, and/or Er could be detected. This latter phase contained many precipitates of the “high-Z” phase. Conversely, this latter phase did not contain precipitates of the UO2-type material, (2) Xe bubbles were found almost exclusively in the “high-Z” phase and in the particular example shown, the Zr-containing variant. Further work in this area should address whether this potential Xe preference for just one of the high-Z phases is indeed a real effect.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Imaging broadband soft X-ray transmission-grating spectrograph for a wavelength range λ > 111 Å

A family of aperiodic Mo/Be multilayer mirrors was designed for maximum uniform reflectivity in the ranges 111 – 138, 111 – 150, 111 – 200 and 111 – 222 Å at normal radiation incidence. The simulations took into account the existence of Mo-on-Be and Be-on-Mo transition layers. It was shown that the reflectivity may be as high as 10% in an octave wavelength range of 111 – 222 Å. We demonstrate the operation of an imaging transmission-grating spectrograph with a broadband Mo/Be multilayer mirror with a uniform reflectivity in the wavelength range 111 – 138 Å synthesised at the Institute for the Physics of Microstructures, RAS. The spectrograph was used to record the line spectra of multiply charged fluorine and magnesium ions from a laser-produced plasma. The role of the grating support structure in the formation of spectral source images is discussed. The spatial resolution of 50 μm is demonstrated. (laser applications and other topics in quantum electronics)

36 MATERIALS SCIENCE↗

Damage initiation and evolution in Al-Si layered microstructures under shock loading conditions at atomic scales

Designing materials with microstructural features like multiphase interfaces have shown significant promise for usage in the next generation defense and nuclear applications. The dynamic response of these interfaces in extreme environments is observed to vary with the deformability of individual phases, which could subsequently alter the favored damage nucleation sites related to spall failure. Here, we investigate the role of spacing of interfaces in a nanocrystalline layered Aluminum-Silicon system on the shock wave propagation behavior, and microstructural and defect evolution using classical molecular dynamics simulations. The simulations are carried out with different hypothetical Al/Si microstructures including variations in the distribution of Si grains as layers in a nanocrystalline Al matrix. The molecular dynamics simulations suggest that spall failure is preferentially initiated at Al/Si interfaces. In addition, for the same system volume and same concentration of Si, a higher number of layers is marked by an increase in shock wave velocity and reduced resistance to spall failure.

36 MATERIALS SCIENCE↗

Parametric Reconstruction of Glass Fiber-reinforced Polymer Composites from X-ray Projection Data—A Simulation Study

We present a new approach to estimate geometry parameters of glass fibers in glass fiber-reinforced polymers from simulated X-ray micro-computed tomography scans. Traditionally, these parameters are estimated using a multi-step procedure including image reconstruction, pre-processing, segmentation and analysis of features of interest. Each step in this chain introduces errors that propagate through the pipeline and impair the accuracy of the estimated parameters. In the approach presented in this paper, we reconstruct volumes from a low number of projection angles using an iterative reconstruction technique and then estimate position, direction and length of the contained fibers incorporating a priori knowledge about their shape, modeled as a geometric representation, which is then optimized. Using simulation experiments, we show that our method can estimate those representations even in presence of noisy data and only very few projection angles available.

42 ENGINEERING↗

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

Introduction The Porous Microstructure Analysis (PuMA) software is an open source framework for image-based simulation, primarily used to determine effective properties based on material microstructure. PuMA was originally developed for the study of NASA thermal protection materials; however, many of the solvers in PuMA have applicability to a broad range of materials science applications. PuMA version 3.2 computes material surface area, pore diameters, effective thermal conductivity, continuum and rarefied tortuosity, and permeability. For anisotropic materials, PuMA can estimate material orientation and compute anisotropic thermal conductivity and elasticity. In this talk, a brief overview of the PuMA software and underlying methods will be presented, as well as some recent and ongoing developments, including the use of immersed boundary methods for image-based simulation and the development of a new weave segmentation tool, called TomoSAM. Cut-Cell method for heat and mass transfer For simulations on complex microstructures, traditional unstructured meshing techniques often prove to be difficult and time-intensive. Voxel-based solvers, which represent the surface as a staircase structure, are relatively simple to implement but can lose accuracy when feature resolution is poor. In this work, we present a novel 3D cut-cell method for solving the variable coefficient Poisson equation on complex microstructures, suitable for the determination of effective thermal conductivity or tortuosity of a material. The method uses a Marching Cubes/Marching Squares surface reconstruction to create cut-cells and determine geometric quantities. A flux-correction method is extended to 3D, with least squares gradient reconstruction, to solve for the boundary fluxes in the cut-cells. Verification cases show the solver achieves globally 2nd order accuracy on complex microstructures. TomoSAM TomoSAM, a module of the PuMA software, has been developed as a plugin for 3D Slicer, a software platform used for 3D image processing and visualization. It utilizes the Segment Anything Model (SAM), a deep learning model capable of identifying objects and generating image masks based on minimal user input. This feature enables efficient segmentation of complex 3D datasets, particularly of woven materials, from tomography or similar imaging methods, reducing the need for manual segmentation.

Tomography↗

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is an open source framework for image-based simulation, primarily used to determine effective properties based on material microstructure. PuMA was originally developed for the study of NASA thermal protection materials; however, many of the solvers in PuMA have applicability to a broad range of materials science applications. PuMA version 3.2 computes material surface area, pore diameters, effective thermal conductivity, continuum and rarefied tortuosity, and permeability. For anisotropic materials, PuMA can estimate material orientation and compute anisotropic thermal conductivity and elasticity. In this talk, a brief overview of the PuMA software and underlying methods will be presented, as well as some recent and ongoing developments, including the use of immersed boundary methods for image-based simulation and the development of a new weave segmentation tool, called TomoSAM. Cut-Cell method for heat and mass transfer For simulations on complex microstructures, traditional unstructured meshing techniques often prove to be difficult and time-intensive. Voxel-based solvers, which represent the surface as a staircase structure, are relatively simple to implement but can lose accuracy when feature resolution is poor. In this work, we present a novel 3D cut-cell method for solving the variable coefficient Poisson equation on complex microstructures, suitable for the determination of effective thermal conductivity or tortuosity of a material. The method uses a Marching Cubes/Marching Squares surface reconstruction to create cut-cells and determine geometric quantities. A flux-correction method is extended to 3D, with least squares gradient reconstruction, to solve for the boundary fluxes in the cut-cells. Verification cases show the solver achieves globally 2nd order accuracy on complex microstructures. TomoSAM TomoSAM, a module of the PuMA software, has been developed as a plugin for 3D Slicer, a software platform used for 3D image processing and visualization. It utilizes the Segment Anything Model (SAM), a deep learning model capable of identifying objects and generating image masks based on minimal user input. This feature enables efficient segmentation of complex 3D datasets, particularly of woven materials, from tomography or similar imaging methods, reducing the need for manual segmentation.

Tomography↗

Li-ion battery design through microstructural optimization using generative AI

Lithium-ion batteries are used across various applications, necessitating tailored cell designs to enhance performance. Optimizing electrode manufacturing parameters is a key route to achieving this, as these parameters directly influence the microstructure and performance of the cells. However, linking process parameters to performance is complex, and experimental or modeling campaigns are often slow and expensive. This study introduces a fast computational optimization framework for electrode manufacturing parameters. A generative model, trained on a small dataset of microstructural images associated with different manufacturing parameters, efficiently generates representative microstructures for new parameters. This model is integrated into a Bayesian optimization loop that includes microstructure generation, characterization, and simulation, aiming to find optimal manufacturing parameters for a particular application. Significant improvement in the energy density of a 4680 cell is achieved through bespoke cell design, highlighting the importance of cell-scale normalization. The framework’s modularity allows its application to various advanced materials manufacturing scenarios.

batteries↗

Multiscale and Machine Learning Modeling for Process-informed Microstructure Prediction in Additively Manufactured Materials using MALAMUTE

The Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy aims to develop and qualify additively manufactured materials for nuclear applications. One key challenge to this is the microstructural variability observed in the additively manufactured products and their impact on the properties and performance of the material in extreme environments. AMMT is using a combination of high-throughput experimental and modeling techniques to accelerate qualification. Conventionally, in-situ and ex-situ characterizations and testing are performed to correlate different aspects of the additive manufacturing process to the final product and its performance. However, adopting a trial-and-error approach to experimentally evaluate the vast range of process parameters required to capture microstructural variability is cost-prohibitive. Modeling and simulation provide a comparatively inexpensive way to understand and correlate the microstructural evolution to the processing conditions. The modeling and simulation work-packages within the AMMT program aims to use physics-based and machine learning models to develop a digital twin for additive manufacturing that can correlate the process conditions to the final product and establish a process-structure-property-performance (PSPP) correlation. The melting and subsequent solidification that occurs during the additive process is a complex phenomenon that requires multiscale multiphysics analysis. This work package focuses on understanding the role of process variabilities on the unique microstructural characteristics of additively manufactured materials. Microstructural features at the subgrain level, such as compositional micro-heterogeneity and dislocation cells, are of particular interest here since they can influence the creep properties and radiation performance. Idaho National Laboratory’s Multiphysics Object-Oriented Simulation Environment (MOOSE), specifically the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) software, provides an ideal platform for developing the multiphysics multiscale model to explore the intricacies of the microstructural evolution during the AM processes within a single framework. Furthermore, given that such full-fidelity simulations can be computationally intensive, reduced order models are necessary to explore the PSPP space for additively manufactured materials in an efficient, reliable, and cost-effective way. This work focuses on capturing the microstructural variabilities at the subgrain level that are often missing in the part-scale models. In fiscal year 2025, we significantly advanced upon our work in the last fiscal year, in terms of the predictive capabilities of the physics-based and ML models, by adding the capabilities to capture subgrain-level micro-segregation during solidification using phase-field model and to predict the time-dependent dynamics of the AM process through the MOGPAR model. The alloy solidification model in MOOSE incorporates the thermodynamic properties and free energy relevant to 316 stainless steel. The model demonstrates the Cr and Ni segregation that occurs during solidification, including that the rate of solidification. The microstructural evolution model is connected to the process conditions via the surrogate model developed in this work. This enables predictions of the final microstructure in conjunctions with the manufacturing process. This work supports AMMT's rapid qualification goals by laying the foundation for an efficient and cost-effective model establishing the PSPP correlation for AM. The generated microstructures and predicted micro-segregation can be used by other work packages under AMMT to evaluate the properties and environmental response of the material at the mesoscale. Thus, this work helps to identify the key microstructural features at the subgrain level that are significant in property and performance predictions of additively manufactured components. This work will also provide inputs to the large-scale process variability models to reevaluate and validate assumptions and simplifications made in the part-scale models. Furthermore, through active learning this work can help identify the data need from both modeling and experimental sides for development of a robust digital twin for additive manufacturing and accelerate the AMMT's qualification efforts.

36 - MATERIALS SCIENCE↗

Quantifying uncertainty in the process-structure relationship for Al–Cu solidification

Phase field method (PFM) is a simulation tool to predict microstructural evolution during solidification and is helpful to establish the process-structure relationship for alloys. The robustness of the relationship however is affected by model-form and parameter uncertainties in PFM. In this paper, the uncertainty associated with the thermodynamic and process parameters of PFM is studied and quantified. Surrogate modeling is used to interpolate four quantities of interests (QoIs), including dendritic perimeter, area, primary arm length, and solute segregation, as functions of thermodynamic and process parameters. A sparse grid approach is applied to mitigate the curse-of-dimensionality computational burden in uncertainty quantification. Polynomial chaos expansion is employed to obtain the probability density functions of the QoIs. The effect of parameter uncertainty on the Al–Cu dendritic growth during solidification simulation are investigated. Here, the results show that the dendritic morphology varies significantly with respect to the interface mobility and the initial temperature.

36 MATERIALS SCIENCE↗

Phase field simulation of martensitic transformation in pre-strained nanocomposite shape memory alloys

We show in this paper how strain engineering alters the fundamental characteristic of a martensitic transformation (MT) and gives it a new set of properties including large quasi-linear elastic strain response with nearly vanishing hysteresis and low elastic modulus. The work is motivated and inspired by a recent experimental study on elastic and inelastic (transformation) strain matching in a pre-strained nano-composite with Nb nanowires embedded in a NiTi shape memory alloy matrix. In particular, we demonstrate by computer simulation that dislocations at Nb/NiTi interfaces produced by the pre-straining are responsible for the unprecedented properties. Microstructural evolution captured in the simulations reveals that local stress fields associated with the dislocations regulate the nucleation and growth of martensite, turning the otherwise sharp, strong first-order transition into a continuous, high-order like transition. The simulations predict that the stress-strain hysteresis and modulus of the composite decrease with increasing amount of pre-strain, which agrees well with the experimental measurement. This study suggests a design strategy by introducing non-uniform stress fields for enhanced properties of shape memory alloys.

36 MATERIALS SCIENCE↗

Continuum stress intensity factors from atomistic fracture simulations

Stress intensity factors (SIFs) are utilized in continuum fracture mechanics to quantify the stress fields surrounding a crack in a homogeneous material in the linear elastic regime. Critical values of the SIFs define an intrinsic measure of the resistance of a material to propagate a crack. At atomic scales, yet, fracture occurs as a series of atomic bonds breaking, differing from the continuum description. As a consequence, a formal analog of the continuum SIFs calculated from atomistic simulations can have spatially localized, microstructural contributions that originate from varying bond configurations. The ability to characterize fracture at the atomic scale in terms of the SIFs offers both an opportunity to probe the effects of chemistry, as well as how the addition of a microstructural component affects the accuracy. We present a novel numerical method to determine SIFs from molecular dynamics (MD) simulations. The accuracy of this approach is first examined for a simple model, and then applied to atomistic simulations of fracture in amorphous silica. MD simulations provide time and spatially dependent SIFs, with results that are shown to be in good agreement with experimental values for fracture toughness in silica glass.

42 ENGINEERING↗

Digital polycrystalline microstructure generation using diffusion probabilistic models

Accurate micromechanical simulation of polycrystalline materials requires a realistic digital representation of the grain scale microstructure. Here, this work demonstrates the use of a generative diffusion probabilistic model for synthesizing single phase polycrystalline realizations. The model performs well and is capable of producing realistic microstructures consisting of not just simple equiaxed structures but also structures exhibiting more complex spatial arrangements. Masked microstructure generation reveals that the model is context aware of morphological descriptors which may be encoded in the latent space. Training on more diverse data sets, with scaled up architectures, may enable development of future models capable of synthesizing even more complex microstructural features.

36 MATERIALS SCIENCE↗

Hydrogen charging and desorption from microstructural viewpoint: A method for deconvoluting hydrogen desorption spectra and unveiling the hydrogen-microstructure interaction

Understanding the interaction of hydrogen with microstructural features in metallic materials is crucial for designing hydrogen-resistant alloys. Although thermal desorption spectroscopy (TDS) is widely used for investigating the hydrogen binding behavior of various microstructural features, its application to face-centered cubic (fcc) metals and alloys that exhibit low hydrogen diffusivity is limited due to the lumped TDS desorption signals. This paper shows that, by coupling a Sofronis–McMeeking type hydrogen transport model with a microstructure-informed finite-element model, TDS data can be deconvoluted to reveal the underlying adsorption–diffusion–desorption processes, hydrogen diffusivity, and trap-binding energies. In conclusion, the austenitic steel SS316L in solution-annealed condition is used as a demonstration material, and we focused on investigating the interaction of deuterium (hydrogen isotope) with grain boundaries, which is difficult to investigate from experiments alone but critical for design of alloys for hydrogen infrastructure.

Finite element simulation, Polycrystalline Microst↗

Understanding grain refinement and intergranular gas bubble evolution in U-10Mo fuel using phase-field modeling

Monolithic U-10Mo fuel undergoes significant microstructural changes in the form of grain refinement and gas bubble formation during burnup, which degrades its mechanical properties. In this talk, I present a phase-field model for microstructure evolution in U-10Mo developed using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Simulations demonstrate that grain refinement initiates at pre-existing grain boundaries (GBs) to eliminate the lattice distortion energy caused by the accumulation of self-interstitial loops. By employing an equation of state for xenon gas, we simulate the evolution of gas bubbles in the polycrystal microstructure. Large, interconnected bubbles are found to form along the triple junctions. The effects of defect production rate, diffusivities and GB mobility on the microstructure evolution are systematically studied. Homogenization is employed on the microstructures to obtain effective elastic constants and diffusivity as a function of fission density. The simulations provide critical insights on microstructure and property degradation in U-10Mo fuel.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Sputter deposition of ultrathick Bi coatings onto rotating substrates for inertial confinement fusion

Sputter-deposited Bi is a potential candidate material for the production of inertial confinement fusion (ICF) hohlraums, which are sphero-cylindrical canisters that serve as the housing for hydrogen-fuel-filled capsules in ICF implosions. The fabrication of hohlraums involves deposition onto rotating sphero-cylindrical substrates and, hence, requires an understanding of oblique angle deposition phenomena. Here, we systematically study the effect of dynamic substrate tilt on properties of Bi films deposited onto rotating planar substrates. Results show that film morphology, porosity, and electrical resistivity depend on substrate tilt. In comparison with films deposited onto substrates adhered to a stationary sample holder, films deposited on rotating substrates have more uniform microstructure with fewer void defects. Based on Monte Carlo simulations, these changes in film morphology and microstructure are attributed to averaging and geometric effects related to sample position and rotation. In conclusion, the fabrication of an ∼ 50-µm-thick Bi film on a rotating sphero-cylindrical mandrel demonstrates a significant step toward the production of Bi hohlraums.

Bismuth↗

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗