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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2006

Microstructural Study of Micron-Sized Craters Simulating Stardust Impacts in Aluminum 1100 Targets

Various microscopic techniques were used to characterize experimental micro- craters in aluminium foils to prepare for the comprehensive analysis of the cometary and interstellar particle impacts in aluminium foils to be returned by the Stardust mission. First, SEM (Scanning Electron Microscopy) and EDS (Energy Dispersive X-ray Spectroscopy) were used to study the morphology of the impact craters and the bulk composition of the residues left by soda-lime glass impactors. A more detailed structural and compositional study of impactor remnants was then performed using TEM (Transmission Electron Microscopy), EDS, and electron diffraction methods. The TEM samples were prepared by Focused Ion Beam (FIB) methods. This technique proved to be especially valuable in studying impact crater residues and impact crater morphology. Finally, we also showed that InfraRed microscopy (IR) can be a quick and reliable tool for such investigations. The combination of all of these tools enables a complete microscopic characterization of the craters.

Leroux, Hugues↗
NASA NTRS2024

Analytical Simulation of Effects of Local Mechanisms and Microstructure on Creep Response of Unidirectional Ceramic Matrix Composites

A micromechanics-based method has been developed to analyze the creep response of uncoated ceramic matrix minicomposites. Although the global stress level at which the creep response is analyzed is lower than the composite proportional limit, local stresses are assumed to be high enough that localized damage is present in the composite in the form of matrix microcracks. To model the composite creep response, including the effects of matrix cracking, a fiber shear lag-based methodology is employed. In this approach, stresses are assumed to vary in the fiber as a function of time and distance from the crack plane. The varying stresses are then used to compute the overall creep strain for the composite. Various assumptions regarding the level of matrix microcracking in the composite and the level of creep in the fiber and matrix in various portions of the composite unit cell are also examined. The creep response of the fiber is modeled using a linear Burgers model. The model is applied to a SiCf/SiC unidirectional minicomposite system. The computed creep results are compared to experimentally obtained values. The effects of the local fiber volume fraction on the overall creep response of the composite are also studied. This work will allow increased understanding of the key material damage mechanisms and load sharing that take place during creep conditions and can be expanded to provide improved analysis methods for full macrocomposites.

Micromechanics↗
DOE OSTI2022

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗
NASA NTRS2008

Multiscale Modeling of Damage Processes in Aluminum Alloys: Grain-Scale Mechanisms

This paper has two goals related to the development of a physically-grounded methodology for modeling the initial stages of fatigue crack growth in an aluminum alloy. The aluminum alloy, AA 7075-T651, is susceptible to fatigue cracking that nucleates from cracked second phase iron-bearing particles. Thus, the first goal of the paper is to validate an existing framework for the prediction of the conditions under which the particles crack. The observed statistics of particle cracking (defined as incubation for this alloy) must be accurately predicted to simulate the stochastic nature of microstructurally small fatigue crack (MSFC) formation. Also, only by simulating incubation of damage in a statistically accurate manner can subsequent stages of crack growth be accurately predicted. To maintain fidelity and computational efficiency, a filtering procedure was developed to eliminate particles that were unlikely to crack. The particle filter considers the distributions of particle sizes and shapes, grain texture, and the configuration of the surrounding grains. This filter helps substantially reduce the number of particles that need to be included in the microstructural models and forms the basis of the future work on the subsequent stages of MSFC, crack nucleation and microstructurally small crack propagation. A physics-based approach to simulating fracture should ultimately begin at nanometer length scale, in which atomistic simulation is used to predict the fundamental damage mechanisms of MSFC. These mechanisms include dislocation formation and interaction, interstitial void formation, and atomic diffusion. However, atomistic simulations quickly become computationally intractable as the system size increases, especially when directly linking to the already large microstructural models. Therefore, the second goal of this paper is to propose a method that will incorporate atomistic simulation and small-scale experimental characterization into the existing multiscale framework. At the microscale, the nanoscale mechanics are represented within cohesive zones where appropriate, i.e. where the mechanics observed at the nanoscale can be represented as occurring on a plane such as at grain boundaries or slip planes at a crack front. Important advancements that are yet to be made include: 1. an increased fidelity in cohesive zone modeling; 2. a means to understand how atomistic simulation scales with time; 3. a new experimental methodology for generating empirical models for CZMs and emerging materials; and 4. a validation of simulations of the damage processes at the nano-micro scale. With ever-increasing computer power, the long-term ability to employ atomistic simulation for the prognosis of structural components will not be limited by computation power, but by our lack of knowledge in incorporating atomistic models into simulations of MSFC into a multiscale framework.

Hochhalter, J. D.↗
NASA NTRS1993

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗
DOE OSTI2023

Optimizing the Composite Cathode Microstructure in All‐Solid‐State Batteries by Structure‐Resolved Simulations

Abstract All‐solid‐state batteries are considered as an enabler for applications requiring high energy and power density. However, they still fall short of their theoretical potential due to various limitations. One issue is poor charge transport kinetics resulting from both material inherit limitations and non‐optimized design. Therefore, a better understanding of the relevant properties of the cathode microstructure is necessary to improve cell performance. In this article, we identify optimization potentials of the composite cathode by structure‐resolved electrochemical 3D‐simulations. In our simulation study, we investigate the influence of cathode active material fraction, density, particle size, and active material properties on cell performance. Special focus is set on the impact of grain boundaries on the cathode design. Based on our simulation results, we can predict target values for cell manufacturing and reveal promising optimization strategies for an improved cathode design.

25 ENERGY STORAGE↗
DOE OSTI2020

A Potts Model parameter study of particle size, Monte Carlo temperature, and “Particle-Assisted Abnormal Grain Growth”

A Potts Model was proposed to account for temperature and particle heterogeneity-dependent Zener Pinning. This was accomplished by relating the random fluctuations of grain boundary position allowed by the Potts Model switching probability at finite simulation temperature to experimentally observed grain growth stagnation behavior. We assume that these fluctuations arise from random fluctuations in the thermodynamic and kinetic properties of the grain boundaries and/or interfaces as they change with temperature. As an application of this model, the grain growth kinetics of U-10 wt. % Mo nuclear fuels were simulated with the input of microstructural images during different heat treatment processes. Simulated average grain growth behavior is in good agreement with experiments.

Frazier, William E.↗
DOE OSTI2022

A micromorphic filter for determining stress and deformation measures from direct numerical simulations of lower length scale behavior

A micromorphic filter is presented for the extraction of stress and deformation measures from underlying direct numerical simulation (DNS) of microstructural mechanical response. The filter is consistent with the micromorphic continuum theory of Eringen and Suhubi (1964), and allows the interrogation of a DNS to generate higher order constitutive models. The filter adopts aspects of the overlap coupling techniques for atomistic-continuum scale-bridging in order to update the nodal degrees of freedom of the finite-element-based filter. Here, two methods for the extraction of quantities of interest from the underlying DNS are presented. The first is a least squares fitting technique which is appropriate for finite-element-style macro-scale domains, and the second is a variationally based approach which generalizes the method to any macro-scale domain which has a variational form. No restrictions on the form of the micro-scale simulation are inherent to the approach, provided that the required quantities (stress, position, etc.) can be provided. The filter is applied to a simple homogeneous medium undergoing homogeneous deformation to demonstrate the recovery of a classical continuum response, what we call the null case. We then apply the framework to the same homogeneous medium with a central and edge penny-shaped crack to show micromorphic effects. We demonstrate that varying the size of the filter causes the micromorphic behaviors to arise and then dissipate as the filter becomes smaller than the crack's length scale. Additionally we show the effect of filter size on a chiral element network along with a plate perforated with evenly spaced and perturbed holes through the thickness. The micromorphic effects are shown to be present for the chiral elements under any loading but become most pronounced when a single element is located within each averaging domain. Micromorphic effects become pronounced for the plate upon the onset of localization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

On the quantitative characterization of weld microstructures

Here, a new technique to quantify grain shape and grain boundary curvature in the fusion zone (FZ) and the heat-effected zone (HAZ) of refractory metals and alloys has been developed. The application of this technique to quantify the effect of welding speed on these microstructural characteristics in simulated molybdenum welds is highlighted in this study. A 50% increase in weld speed more than doubles the fraction of grain segments in the fusion zone with orientations between 00 to 300 to the normal to the welding direction, with a corresponding 50% decrease in the fraction of grain segments oriented along the welding direction. This increase in weld speed also results in a reduction in grain size from approximately 120 µm to 84 µm along the WD at the FZ centerline. Grain boundary curvature is negative near the weld centerline and the edge of the fusion zone, indicating that with increasing distance from the weld centerline, the grain segments progressively rotate towards a direction normal to the welding direction at both locations. Presence of regions with relatively straight grain segments indicated by curvatures close to zero were identified in between the weld centerline and the edge of the FZ at all weld speeds.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Comparison of Commonly-Used Microwave Radiative Transfer Models for Snow Remote Sensing

This paper reviews four commonly-used microwave radiative transfer models that take different electromagnetic approaches to simulate snow brightness temperature (T(sub B)): the Dense Media Radiative Transfer - Multi-Layer model (DMRT-ML), the Dense Media Radiative Transfer - Quasi-Crystalline Approximation Mie scattering of Sticky spheres (DMRT-QMS), the Helsinki University of Technology n-Layers model (HUT-nlayers) and the Microwave Emission Model of Layered Snowpacks (MEMLS). Using the same extensively measured physical snowpack properties, we compared the simulated T(sub B) at 11, 19 and 37 GHz from these four models. The analysis focuses on the impact of using different types of measured snow microstructure metrics in the simulations. In addition to density, snow microstructure is defined for each snow layer by grain optical diameter (Do) and stickiness for DMRT-ML and DMRT-QMS, mean grain geometrical maximum extent (D(sub max)) for HUT n-layers and the exponential correlation length for MEMLS. These metrics were derived from either in-situ measurements of snow specific surface area (SSA) or macrophotos of grain sizes (D(sub max)), assuming non-sticky spheres for the DMRT models. Simulated T(sub B) sensitivity analysis using the same inputs shows relatively consistent T(sub B) behavior as a function of Do and density variations for the vertical polarization (maximum deviation of 18 K and 27 K, respectively), while some divergences appear in simulated variations for the polarization ratio (PR). Comparisons with ground based radiometric measurements show that the simulations based on snow SSA measurements have to be scaled with a model-specific factor of Do in order to minimize the root mean square error (RMSE) between measured and simulated T(sub B). Results using in-situ grain size measurements (SSA or D(sub max), depending on the model) give a mean T(sub B) RMSE (19 and 37 GHz) of the order of 16-26 K, which is similar for all models when the snow microstructure metrics are scaled. However, the MEMLS model converges to better results when driven by the correlation length estimated from in-situ SSA measurements rather than D(sub max) measurements. On a practical level, this paper shows that the SSA parameter, a snow property that is easy to retrieve in-situ, appears to be the most relevant parameter for characterizing snow microstructure, despite the need for a scaling factor.

Snow↗
DOE OSTI2025

Simultaneously Improving Process Efficiency and Mechanical Properties in Aluminum Alloys with Applied Magnetic Fields

In the aluminum heat treatment sector, large resistive batch furnaces are used for many heat treatment processes and are generally always running, operating at approximately 20–30% efficiency. The primary avenue for energy savings lies in shortening heat treatment times, which can be 24+ hours for solutionizing (450–625 °C) and aging treatments (150–250 ℃), depending on the alloy and specification requirements. The application of high magnetic fields (<9 T) has shown to significantly reduce aluminum alloy heat treatment times, while achieving enhanced mechanical properties. Here, we explore the effects of thermomagnetic processing on the heat treatment response of a commercial aluminum alloy A390.1. Heat treatment times were reduced by >80%, while achieving ~10% improvement in microhardness, and these improvements persisted after longer thermal exposures simulating use. The microstructural evolution and resulting properties will be discussed.

Kesler, Michael [ORNL] (ORCID:0000000161423128)↗
DOE OSTI2024

Integration of CeO 2 -Based Memristor with Vertically Aligned Nanocomposite Thin Film: Enabling Selective Conductive Filament Formation for High-Performance Electronic Synapses

The CeO 2 -based memristor has attracted significant attention due to its intrinsic resistive switching (RS) properties, large on/off ratio, and great plasticity, making it a promising candidate for artificial synapses. However, significant challenges such as high power consumption and poor device reliability hinder its broad application in neuromorphic microchips. To tackle these issues, in this work, we design a novel bilayer (BL) memristor by integrating a CeO 2 -based memristor with a Co-CeO 2 vertically aligned nanocomposite (VAN) layer and compare it with the single layer (SL) memristor. Preliminary electrical testing reveals that the BL memristor offers a reduced set/reset voltage (~67% lower), a higher on/off ratio (~5 × 10 2 ), enhanced device reliability, and improved device-to-device variation compared to the SL memristor. Insight from COMSOL simulation, coupled with microstructural analysis, provides a comprehensive elucidation on how the VAN layer facilitates the selective conductive filament (CF) formation. Subsequently, the plasticity of the BL memristor is evaluated through long-term potentiation/depression (LTP/LTD), paired-pulse facilitation (PPF), and spike-time-dependent plasticity (STDP). The spiking neural network (SNN) built upon the BL memristor achieves remarkable accuracy (~94%) after only 12 iterations, underscoring its potential for high-performance neural networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models

GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.

42 ENGINEERING↗
NASA NTRS2012

Characterization and Computational Modeling of Minor Phases in Alloy LSHR

The minor phases of powder metallurgy disk superalloy LSHR were studied. Samples were consistently heat treated at three different temperatures for long times to approach equilibrium. Additional heat treatments were also performed for shorter times, to assess minor phase kinetics in non-equilibrium conditions. Minor phases including MC carbides, M23C6 carbides, M3B2 borides, and sigma were identified. Their average sizes and total area fractions were determined. CALPHAD thermodynamics databases and PrecipiCalc(TradeMark), a computational precipitation modeling tool, were employed with Ni-base thermodynamics and diffusion databases to model and simulate the phase microstructural evolution observed in the experiments with an objective to identify the model limitations and the directions of model enhancement.

Jou, Herng-Jeng↗
DOE OSTI2020

A FEASIBILITY STUDY ON CRACK REPAIR IN AUSTENITIC STAINLESS STEEL DRY STORAGE CANISTERS USING ISOTHERMAL FRICTION STIR WELDING

In recent time, solid state crack repair techniques in spent nuclear fuel dry storage canisters (SNF-DSC) have garnered significant interest. Plates of austenitic stainless steels are usually arc welded to fabricate the cylindrical canisters which are prone to chloride-induced stress corrosion cracking. Friction stir welding (FSW) is considered to be a potential candidate for damage repair. In this work, electrical discharge machining was used to simulate cracks in 304L stainless steel plates of 12.7 mm thickness. Subsequently, isothermal FSW was carried out at two different temperatures (725 and 825°C) in order to repair the simulated crack. Microstructure-mechanical property correlations in FSWed plates were studied in detail. Significant grain refinement observed in the stir zone (SZ) was attributed to dynamic recrystallization occurring during FSW. A comparison of Vickers microhardness showed higher hardness and homogeneity in the 725°C SZ than its 825°C counterpart. The fraction of S3 boundary is found to be low in the SZ as compared to base metal (BM). Yield strength of the joints was found to be approximately 100 MPa higher than that of the BM accompanied by a drop in ductility by a factor of 2. Residual stress profile across 825°C weld was measured using x-ray diffraction.

Bhattacharyya, Madhumanti↗
NASA NTRS2009

A Method for Combining Experimentation and Molecular Dynamics Simulation to Improve Cohesive Zone Models for Metallic Microstructures

Fracture processes within a material begin at the nanometer length scale at which the formation, propagation, and interaction of fundamental damage mechanisms occur. Physics-based modeling of these atomic processes quickly becomes computationally intractable as the system size increases. Thus, a multiscale modeling method, based on the aggregation of fundamental damage processes occurring at the nanoscale within a cohesive zone model, is under development and will enable computationally feasible and physically meaningful microscale fracture simulation in polycrystalline metals. This method employs atomistic simulation to provide an optimization loop with an initial prediction of a cohesive zone model (CZM). This initial CZM is then applied at the crack front region within a finite element model. The optimization procedure iterates upon the CZM until the finite element model acceptably reproduces the near-crack-front displacement fields obtained from experimental observation. With this approach, a comparison can be made between the original CZM predicted by atomistic simulation and the converged CZM that is based on experimental observation. Comparison of the two CZMs gives insight into how atomistic simulation scales.

Hochhalter, J. D.↗