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 631 records · Page 35

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗

USE OF TWO-PISTON SPLAT QUENCHING TO INVESTIGATE & CHARACTERIZE THE IMPACT OF COMPOSITIONAL VARIATIONS ON RAPID SOLIDIFICATION MICROSTRUCTURES & SUB-MICROSCALE FEATURES IN STAINLESS STEEL ALLOY.

The objective of this dissertation was to use two-piston splat quenching (SQ) to investigate the impact of compositional modifications on the solidification and microstructure of rapidly solidified austenitic stainless steels (SS) and to demonstrate the ability of SQ to quickly and effectively simulate rapid solidification conditions similar to those found in powder bed fusion (PBF) additive techniques. PBF techniques like laser powder bed fusion (LPBF) are being implemented across a breadth of research and industrial applications to create parts with complex geometries and performance capabilities while pushing the current limits of processing conditions and understanding of material systems. In this work, SQ was used to experimentally produce rapid solidification in 20+ unique austenitic SS compositions with systematic variations of the chrome and nickel equivalency ratio (Cr/Nieq) through targeted compositional modifications. From the targeted change of Cr, Ni, and Mo concentrations in rapidly solidified SS alloys, the ferrite solidification mode was found to be the primary solidification mode at significantly lower Cr/Nieq than previously predicted for RS. Also, decreasing concentrations of Fe at a constant Cr/Nieq ratio (i.e., different Fe isopleths), or increased Mo concentrations at a constant Cr/Nieq ratio were found to suppress the ferrite to austenite massive transformation when compared to alloys with lower concentrations at the same Cr/Nieq. Using an established empirical relationship between cell size and cooling rate, the SQ technique was estimated to produce cooling rates between 106 and 108 K/s. Thermal gradients were extracted from 2-D heat transfer simulations of the SQ solidification event and used with these cooling rates to produce solidification rate estimates for SQ which were between 0.4-1.6m/s. The primary solidification mode was observed to be the determining factor in which elements segregated to the cell boundaries during RS, for which Cr and Mo were the main elements to segregate during primary austenite solidification and Ni during primary ferrite solidification. Finally, the solidification rates and conditions produced by SQ experiments resulted in similar microstructures, features, and microsegregation to what was found in LPBF samples of the same feedstock.

Hasenbusch, Zachary↗

Using machine-learning to understand complex microstructural effects on the mechanical behavior of Ti-6Al-4V alloys

Structural materials properties are highly dependent on their microstructure. Their microstructure is in turn affected by multiple fabrication and thermo-mechanical treatment parameters, all of which conform a highly-dimensional parametric space with often hidden correlations that are difficult to extract by experimentation alone. This is particularly true for alloys of the dual-phase Ti-6Al-4V family, with their greatly complex and rich microstructures, which combine several intrinsic length scales associated with multiple grain and subgrain structures, grains with different crystal lattices (α and β phases), and complex chemistry. In this paper we use a comprehensive set of machine learning techniques to develop predictive tools relating the yield strength and hardening rate of these alloys to a set of input parameters covering extensive ranges. The data generator is a finite-element crystal plasticity model for polycrystal deformation that takes into account slip anisotropy and employs standard dislocation evolution models for the α and β phases of Ti-based alloys. Our dataset includes over two thousand independent simulations and is used to train the machine learning models, which are then used to establish correlations between microstructural parameters and the alloys’ mechanical response. Our results point to the most influential parameters affecting yield strength and hardening rate, information that can then be used to guide experimental synthesis and characterization efforts to save time and resources.

36 MATERIALS SCIENCE↗

Microstructural Model for Creep-Fatigue Interaction in Grade 91 Steel

This report describes the completion of a crystal plasticity finite element simulation framework for creep-fatigue deformation and damage in Grade 91 steel. In addition to a description of the framework and a set of verification calculations, the report details an initial set of simulations applying the framework to examine low cycle fatigue and creep-fatigue interaction in Grade 91 at 600° C. Grade 91 is a likely material for future advanced microreactors and creep-fatigue will be a limiting structural damage mechanism, particularly for mobile reactors or load following systems. The framework described here could provide a better understanding of creep-fatigue damage mechanisms and directly examine creep, fatigue, and creep-fatigue interaction for Grade 91 at realistic component service conditions. This in turn could lead to improved design methods and component life estimation approaches, leading to more efficient and economical component designs.

36 MATERIALS SCIENCE↗

Meshfree simulation and prediction of recrystallized grain size in friction stir processed 316L stainless steel

Friction stir processing (FSP) is a promising solid-phase microstructural modification technique that can repair and enhance damaged stainless steel surfaces exposed to harsh environments. The quality of the repaired material is closely correlated to the recrystallized grain size in the stir zone (SZ), which is influenced by the thermomechanical conditions dictated by FSP process parameters. Thus, establishing a reliable relationship between these parameters and recrystallized grain size in the SZ is crucial for optimizing repair quality. However, existing experimental approaches often rely on indirect temperatures measured far from the SZ, along with rough strain rate estimations, which are imprecise and time-consuming. Meanwhile, existing mesh-based modeling methods usually face numerical challenges when dealing with the large material deformations inherent in FSP. Here, to address these issues, this study introduces a meshfree process model for FSP based on the smoothed particle hydrodynamics (SPH) method, aimed at predicting process conditions under different parameters. The model is validated using experimental data from 11 combinations of tool traverse and rotation speeds on 316 L stainless steel. Correlations between process parameters, material flow, temperature, strain, strain rate, and recrystallized grain size are revealed through SPH simulations and electron backscatter diffraction (EBSD) imaging. The results show that in situ SZ temperatures range from 1071 to 1322°C, which exceed the tool temperature by over 300°C. Furthermore, SZ temperature, strain rate, and grain size increase monotonically with higher tool temperature and faster traverse speed. A relationship is then established between the model-predicted Zener-Hollomon parameter and the recrystallized grain size based on EBSD data, expressed as ln(d) = -0.364 ln(Z) + 14.673. Finally, this relationship exhibits satisfactory accuracy with errors of less than 26.9% in predicting grain sizes at various SZ locations, which offers valuable insights for optimizing FSP repair processes for 316 L stainless steel.

316L stainless steel↗

Discerning the effect of various irradiation modes on the corrosion of Zircaloy-4

Using proton irradiation, this study investigates the individual influence of several factors on the corrosion kinetics of Zircaloy-4 in a hydrogenated water environment simulating a Pressurized Water Reactor (PWR). Using both simultaneous irradiation-corrosion and autoclave corrosion, we separately examine: (i) the effect of pre-irradiation on modifying the structure of the material, (ii) the impact of irradiation on creating defects in the growing oxide layer during corrosion, and (iii) the influence of irradiation on increasing the corrosion potential through radiolysis during corrosion. To replicate neutron-irradiated microstructure, two proton pre-irradiation schedules were employed: Schedule 1 (isothermal irradiation at 350 °C to 5 dpa) to simulate high-temperature PWR conditions, and Schedule 2 (two-step process: irradiation to 2.5 dpa at -10 °C followed by 2.5 dpa at 350 °C) to simulate lower temperature PWR and Boiling Water Reactor (BWR) conditions. Long-term autoclave corrosion testing over 360 days at 320 °C revealed no significant difference between unirradiated samples and those pre-irradiated according to either schedule, with all samples exhibiting sub-cubic kinetics within the pre-transition regime. Irradiated samples underwent Simultaneous Irradiation Corrosion (SIC) tests, corroding in 320 °C water while being irradiated with protons. Corrosion was found to accelerate in all SIC-tested samples relative to autoclave conditions, with the greatest increase observed in non-pre-irradiated samples. Pre-irradiation with either schedule resulted in a slower corrosion rate compared to non-pre-irradiated regions under SIC conditions. The degree of radiolysis observed in the SIC tests surpassed typical PWR conditions, approaching levels found in BWRs. Radiolysis products were identified as the primary contributors to accelerated corrosion, corroborated by radiolysis bar tests. Furthermore, these findings underscore the intricate interactions between irradiation, corrosion, and water chemistry in determining Zircaloy-4 corrosion kinetics within nuclear reactor environments.

Pressurized Water Reactor↗

Machine-learning-based, online estimation of ceramic’s microstructure upon the laser spot brightness during laser sintering

The ceramic microstructure strongly influences its properties. During manufacturing, the online monitoring of microstructure is critical to ensure the desired material properties. So far, the microstructure on the relevant scale is usually characterized offline using scanning electron microscopy (SEM), which is time and cost-consuming. In this work, we demonstrate a cost-effective, machine learning (ML)-based approach to simulate the SEM micrographs in real-time from the laser spot brightness. We experimentally observed a strong correlation between the laser spot brightness and the corresponding microstructure at the exact locations. The brightness values obtained from thermal emission images and the corresponding SEM micrographs were used in the training datasets. The ML algorithm was a style-based conditional generative adversarial network (CGAN). After training, the ML model could generate high-fidelity microstructure images within 0.1 seconds based on in-situ captured brightness at the laser sintering spot. We used the average grain sizes as the metric to evaluate the accuracy of the ML-predicted micrographs. Here, the ML-predicted microstructures were in good agreement, with less than 5% in difference from the real SEM images. In conclusion, we demonstrate the cost-effective, online microstructure estimation during laser sintering with a simple setup (a camera, a regular computer, and the ML model).

08 HYDROGEN↗

Effects of composition and canister centerline cooling on microstructure, phase distribution, and chemical durability of dehalogenated iron phosphate waste forms

This paper discusses the effects of composition and canister centerline cooling (CCC) on the microstructures, phase distribution, and chemical durability of dehalogenated iron phosphate waste forms starting from chloride-based salt simulants including KCl (Simple-1), NdCl3 (Simple-2), SrCl2 (Simple-3), LaCl3 (Simple-4), LiCl (Simple-5), CsCl (Simple-6), as well as a complex simulant containing LiCl, KCl, NaCl, CsI, SrCl2, CeCl3, and NdCl3 called ERV2. The simulants represent salt constituents present in salt wastes from electrorefiners used for electrochemical reprocessing. These salts are reacted with ammonium dihydrogen phosphate (ADP) resulting in a NH4Cl (dechlorination) byproduct that is captured and phosphate product containing oxides of the salt cations. This product is further stabilized using Fe2O3 through vitrification at higher temperatures. The results presented herein describe the properties of CCC-treated materials following dechlorination and vitrification including phase identification and quantification, elemental distributions, and chemical durabilities. The chemical durability in the CCC-cooled samples are lower durability than quenched materials for the ERV2 samples. The phases formed upon CCC cooling are very complex compositionally and microstructurally.

iron phosphate glass, iron phosphate waste forms, ↗

Microstructural Pattern Formation during Far-from-Equilibrium Alloy Solidification

We introduce a new phase-field formulation of rapid alloy solidification that quantitatively in- corporates non-equilibrium effects at the solid-liquid interface over a very wide range of interface velocities. Simulations identify a new dynamical instability of dendrite tip growth driven by solute trapping at velocities approaching the absolute stability limit. They also reproduce the formation of the widely observed banded microstructures, revealing how this instability triggers transitions between dendritic and microsegregation-free solidification. Predicted band spacings agree quantita- tively with observations in rapidly solidified Al-Cu thin films.

36 MATERIALS SCIENCE↗

Modeling Induction Stirring and Particle Tracking in Molten Uranium

Two independent numerical models have been developed to simulate the behavior of carbon impurities in molten uranium metal. Informed by experimental parameters, one model was created using the commercial software Star-CCM+ and compared with another developed using open-source codes, including OpenFOAM, Finite Element Method Magnetics (FEMM) and a First Passage Kinetic Monte Carlo (FPKMC) approach. The target experimental system features a 404g uranium metal charge containing an average carbon concentration of 139 ppm which was melted in a vacuum induction furnace at 1400° C then resolidified. The microstructures of the uranium and its impurities before and after melting have been characterized and reported separately. Prior to simulating the uranium-carbon system described, the numerical models were validated using a previously published nonradioactive experimental system to ensure agreement with expected output values. Focus has been placed on modeling velocity fields under induction stirring and impurity particle trajectories.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Investigation of the effect of etch pits on the material properties of carbon fiber structures

It is well known that etch pits are formed on carbon surfaces during their consumption in reactive gas environments such as during ablation of carbon-based thermal protection system (TPS) materials. These pits start out from atomic defects on the carbon surface, grow in size and then ultimately lead to the complete removal of the solid carbon material [1, 2]. However, the effect of such etch pits on the material properties of carbon ablators such as FiberForm are poorly understood, as the maximum size of pits are on the order of microns in size, which makes them hard to observe experimentally. Hence, we have developed a module within the Porous Microstructure Analysis code PuMA [3] that can generate etch pits on an arbitrary material surface, with a prescribed distribution of pit size and density. Figure 1 shows a FiberForm sample with artificially generated pits of varying radii and pit density. We use this newly developed module to calculate material properties of FiberForm, which is the base material for one of the most commonly used TPS materials - Phenolic Impregnated Carbon Ablator (PICA). The detailed microstructure of FiberForm in PuMA will either be generated synthetically or obtained from X-ray microtomography [4]. Our simulations show that material properties such as the thermal conductivity and tortuosity decrease with increasing degree of pitting. This will be helpful to more accurately predict the degradation of carbon-based TPS during ablation, and furthermore may allow us to identify previously unknown TPS failure modes due to pitting.

PuMA↗

Gradient microstructure induced shear band constraint, delocalization, and delayed failure in CuZr nanoglasses

The mechanical properties of metallic glasses (MGs) can be significantly tuned by introducing a nanoscale glass-glass interface microstructure. Here, we design three types of Cu 64 Zr 36 MG seamless gradient microstructures by combining glassy grain sizes from 3 to 15 nm. We perform tensile loading molecular dynamics simulations on the named gradient nanoglasses to characterize their mechanical behavior. In all cases, results indicate the generation of gradient plasticity following the gradient in grain size, i.e. local plasticity increases from large to small grain size regions. Loading parallel to the gradient direction results in initial broad necking followed by catastrophic failure. In contrast, loading perpendicular to the gradient direction results in diffused shear band propagation from large to small grain size regions, delaying the generation of a dominant shear band and further delocalizing plastic deformation. In conclusion, the results demonstrate the synergistic effects generated by heterogeneous gradient designs in the mechanical behavior of metallic glasses.

36 MATERIALS SCIENCE↗

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗

A study of subsurface crack initiation produced by rolling contact fatigue

Results of subsurface crack initiation studies produced by pure rolling contact fatigue in 7075-T6 aluminum alloy are presented in this article. Microstructural changes and subsequent crack initiation below the contacting surface in cylindrical test specimens subjected to repeated rolling contact are illustrated. The rolling conditions are simulated in a three-dimensional elastic-plastic finite element model in order to estimate the plastic strains and residual stresses in the test material. The numerically estimated distribution of plastic strains in the model correlate well with the extent of microstructural changes observed in the test specimen. Results also indicate that a combination of plastic strains and low values of residual stresses is conducive to subsurface crack initiation and growth.

Kumar, Arun M.↗

Multi-scale simulation of high burnup UO2 nuclear fuel during loss-of-coolant accident conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized that pulverization is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted, and the results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Micromechanical Surrogate Machine Learning Model for Creep Deformation Modeling

Process variability during the manufacture of gas turbine engine hot section components can significantly affect the material’s resulting microstructure. In casting, for instance, geometric variation within a component (thin sections versus thick sections, radial location) influences cooling rates and the resulting grain size. The high temperature creep response is known to be sensitive to grain size owing to a diffusional creep mechanism which occurs more readily along grain boundaries. Microstructural variation correspondingly drives mechanical behavior which propagates into component scale performance uncertainty. These factors are essential when planning inspection, maintenance, and repair strategies within a reliability framework. These benefits provide opportunities to increase overall energy efficiency through refined margins. Critically, there is an opportunity to bolster existing data-driven reliability models using physics-driven process-structure-property relations. Here we present recent work establishing a framework for evaluating the probabilistic creep performance of high-temperature materials. A novel microstructure-sensitive crystal plasticity finite element model is established that captures both grain boundary and crystallographic deformation effects. The computationally expensive physics model is calibrated using a statistical approach and this high-fidelity model is subsequently used to train a computationally efficient machine learning surrogate model. The surrogate model is essential for sampling a large ensemble of simulated structure-property pair results. The ensemble data are then mined to extract salient trends to be incorporated into a microstructure-sensitive reliability model. The proposed approach represents a novel way to capture microstructure-sensitive trends from physics-based models within a modern reliability framework.

Fernandez-Zelaia, Patxi [ORNL]↗