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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Three Dimensional Numerical Simulation and Characterization of Crack Growth in the Weld Region of a Friction Stir Welded Structure

Friction stir welding (FSW) fabrication technology is being adopted in aerospace applications. The use of this technology can reduce production cost, lead-times, reduce structural weight and need for fasteners and lap joints, which are typically the primary locations of crack initiation and multi-site fatigue damage in aerospace structures. FSW is a solid state welding process that is well-suited for joining aluminum alloy components; however, the process introduces residual stresses (both tensile and compressive) in joined components. The propagation of fatigue cracks in a residual stress field and the resulting redistribution of the residual stress field and its effect on crack closure have to be estimated. To insure the safe insertion of complex integral structures, an accurate understanding of the fatigue crack growth behavior and the complex crack path process must be understood. A life prediction methodology for fatigue crack growth through the weld under the influence of residual stresses in aluminum alloy structures fabricated using FSW will be detailed. The effects and significance of the magnitude of residual stress at a crack tip on the estimated crack tip driving force are highlighted. The location of the crack tip relative to the FSW and the effect of microstructure on fatigue crack growth are considered. A damage tolerant life prediction methodology accounting for microstructural variation in the weld zone and residual stress field will lead to the design of lighter and more reliable aerospace structures

Seshadri, Banavara R.↗
DOE OSTI2026

Evaluating Fracture Behavior of Bioinspired Alumina-YSZ Composites through Static and Dynamic Mechanical Testing

Ceramic materials are known for their high hardness and strength but are limited by low toughness and sudden failure. Inspired by the microstructure of dental enamel, which features undulating rods that promote crack deflection and energy absorption, architected specimens were designed and manufactured. A bespoke ball-on-ring (BoR) testing apparatus measured the static biaxial rupture strength of these bioinspired materials. Additionally, impact testing with spherical steel projectiles assessed their behavior under dynamic conditions. Monolithic alumina specimens were compared to architected alumina reinforced with yttria-stabilized zirconia (YSZ) rods produced via direct ink write 3D printing. BoR tests revealed that monolithic alumina had a Weibull modulus of 10.53 and a characteristic strength of 372.3 MPa, while the composite specimens showed a Weibull modulus of 5.46 and a characteristic strength of 213.3 MPa. Although the composite failed at lower stresses and projectile velocities, it exhibited notable crack deflection and fracture resistance, with cracks being effectively interrupted by the rod inclusions. The composite specimens also resulted in fewer and larger fragments upon impact. Finite element simulations confirmed the effectiveness of the rod structures in enhancing fracture resistance.

Bioinspired materials↗
DOE OSTI2020

Bridging Scales to Model Reactive Diffusive Transport in Porous Media

Published by ECS. Two novel scale-bridging algorithms to model reaction-diffusion transport in porous media are presented. The algorithms are based on direct numerical simulations and couple the information of a micro-scale model, which accounts for the large field of view provided by micro X-ray computed tomography (X-ray CT), and a nano-scale model, which locally resolves transport in the fine structure extracted from nano X-ray CT. The micro-scale model is discretized in the through-plane direction into a 1D grid, where effective properties and internal boundaries are determined based on the results from the nano-scale model. The validated algorithms are used to examine transport of oxygen in precious group metal-free electrodes considering both zero- and first-order kinetics. Unlike conventional methods, the results show that the effective diffusivity is not a passive property but increases in regions where the reaction-rate coefficient is large. The proposed algorithms account for the multiscale coupling of reaction-diffusion transport and material microstructure, thus improving the predictions compared to conventional methods.

42 ENGINEERING↗
DOE OSTI2021

How Austenitic Is a Martensitic Steel Produced by Laser Powder Bed Fusion? A Cautionary Tale

Accurate phase fraction analysis is an essential element of the microstructural characterization of alloys and often serves as a basis to quantify effects such as heat treatment or mechanical deformation. Additive manufacturing (AM) of metals, due to the intrinsic nonequilibrium solidification and spatial variability, creates additional challenges for the proper quantification of phase fraction. Such challenges are exacerbated when the alloy itself is prone to deformation-induced phase transformation. Using commonly available in-house X-ray diffraction (XRD) and electron backscatter diffraction (EBSD) and less commonly used synchrotron-based high-energy X-ray diffraction, we characterized nitrogen-atomized 17-4 precipitation-hardening martensitic stainless steel, a class of AM alloy that has received broad attention within the AM research community. On the same build, our measurements recovered the entire range of reported values on the austenite phase fractions of as-built AM 17-4 in literature, from ≈100% martensite to ≈100% austenite. Aided by Calphad simulation, our experimental findings established that our as-built AM 17-4 is almost fully austenitic and that in-house XRD and EBSD measurements are subject to significant uncertainties created by the specimen’s surface finish. Hence, measurements made using these techniques must be understood in their correct context. Our results carry significant implications, not only to AM 17-4 but also to AM alloys that are susceptible to deformation-induced structure transformation and suggest that characterizations with less accessible but bulk sensitive techniques such as synchrotron-based high energy X-ray diffraction or neutron diffraction may be required for proper understanding of these materials.

17-4 steel↗
DOE OSTI2021

A constitutive model for glass–ceramic materials

Glass–ceramics have received recent attention for use in glass–ceramic to metal hermetic seals. Due to their heterogeneous microstructure, these materials exhibit a number of advantageous responses over conventional glass based seals. Key amongst them is the possibility of a controllable thermal strain response and apparent coefficient of thermal expansion which may be used to minimize thermally induced residual stresses for aforementioned seals. Furthermore, these behaviors result from an inorganic glass matrix and variety of crystalline ceramic phases including silica polymorph(s) that may undergo reversible solid-to-solid transformations with associated inelastic strain. Correspondingly, these materials exhibit complex thermomechanical responses associated with multiple inelastic mechanisms (viscoelasticity and phase transformation). While modeling these behaviors is essential for developing and analyzing the corresponding applications, no such model exists. Therefore, in this work a three-dimensional continuum constitutive model for glass–ceramic materials combining these various inelastic mechanisms is developed via an internal state variable approach. A corresponding fully implicit three dimensional numerical formulation is also proposed and implemented. The model is used to simulate existing experiments and validate the proposed formalism. As an example, the simple seal problem of a glass–ceramic seal inside a concentric metal shell is explored. Finally, the impact of the cooling rates, viscoelastic shift factors, and inelastic strain on final residual stress state are all investigated and the differing contributions highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stability and Mobility of Disconnections in Solute Atmospheres: Insights from Interfacial Defect Diagrams

This Letter explores the stability of disconnections (step-dislocation defects) at grain boundaries in binary alloys. We introduce interfacial defect diagrams, derived from atomistic simulations and segregation theory, to predict the stability of disconnections in the temperature-solute concentration phase space and relate it to governing segregation mechanisms. These diagrams reveal multiple stability regimes influenced by solute-induced clustering and pinning effects impacting the thermal migration of disconnections and offering insights into their thermodynamics and kinetic properties.

Crystal structure↗
DOE OSTI2026

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗
DOE OSTI2022

Solute segregation and precipitation across damage rates in dual-ion–irradiated T91 steel

We report that dual-ion irradiations using 5.0 MeV defocused Fe 2+ ions and co-injected energy degraded 2.00 to 2.85 MeV He 2+ ions were conducted on a Fe9CrMo ferritic-martensitic steel T91 to 17 dpa at a damage rate range of 5 × 10 -5 dpa/s to 3 × 10 -3 dpa/s at 445°C, followed by characterization of the microstructure using conventional and scanning transmission electron microscopy. Radiation induced Ni/Si clusters and radiation induced segregation were quantified using energy dispersive X-ray spectroscopy at each condition and were compared with the same material irradiated in the BOR-60 reactor and in the as-received condition. No significant Cr segregation was found at lath boundaries after dual-ion irradiation, while Ni and Si enrichments both decreased with increasing damage rate leading to a sharp decrease in the density of Ni/Si clusters with damage rate. Increased point defect recombination at higher ion damage rates likely reduced the Ni/Si cluster density compared with BOR-60. Although the overall vacancy concentration and diffusion are enhanced by the irradiation damage rate, the lack of time for thermal diffusion and ballistic displacements of solutes are significant limiting factors for Ni/Si cluster formation. This work demonstrates the effect of irradiation damage rate on elemental segregation and clustering when using ion irradiation to simulate reactor irradiation.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Equiaxed Dendritic Solidification Experiment (EDSE)

The objective of the research is to quantitatively determine and understand the fundamental mechanisms that control the microstructural evolution during equiaxed dendritic solidification. A microgravity experiment will be conducted to obtain benchmark data on the transient growth and interaction of up to four equiaxed crystals of a pure and transparent metal analog (succinonitrile, SCN) under strictly diffusion-dominated conditions. Of interest in the experiment are the transient evolution of the primary and secondary dendrite tip speeds, the dendrite morphology and solid fraction, the tip selection criterion, and the temperature field in the melt for a range of interaction "strengths" between the crystals. The experiment extends the microgravity measurements of Glicksman and co-workers isothermal dendritic growth experiment (IDGE) for steady growth of a single dendrite to a case where growth transients are introduced due to thermal interactions between neighboring dendrites - a situation closer to actual casting conditions. Corresponding Earth-based experiments will be conducted to ascertain the influence of melt convection. The experiments are supported by a variety of analytical models and numerical simulations. The data will be used to develop and test theories of transient dendritic growth and the solidification of multiple interacting equiaxed crystals in a supercooled melt.

Beckermann, C.↗
DOE OSTI2026

Development of a coupled experimental–computational approach for engineering optimization of spout-fluidized bed particle coating systems

The design of spout-fluidized bed (SFB) coating systems for nuclear particle fuels typically relies on trial-and-error processes, comprising iterative and time-consuming coating deposition experiments and post-deposition characterization. At an engineering scale, this approach to guided SFB system design is inefficient, highlighting the need for streamlined experimental methodologies which can correlate fluidization conditions to downstream coating outcomes. In this study, we combine time-resolved particle image velocimetry (PIV) with CFD–DEM simulations to benchmark hydrodynamic behavior in a 3D spout-fluidized bed. By exploiting easily accessible optical measurements of particle motion at the bed wall and within the spouting region, we obtain quantitative velocity fields that can be directly compared with model predictions of the occluded bed region, without resorting to complex imaging and characterization techniques such as X-ray or magnetic resonance tomography. Experimental benchmarking reveals strong agreement between CFD–DEM and PIV in the spout and annulus regions, while discrepancies near the wall highlight areas for future model development. Here, the proposed integrated experimental–numerical framework will enable a direct connection between measured variables and numerically predicted fluidization performance of dense, surrogate nuclear particle fuel feedstock such that experimental SFB component design can be rapidly evaluated, informing design decisions for nozzle geometry and operating conditions. Future work will extend this framework by correlating quantified fluidization metrics across nozzle geometries and operating conditions with the resulting coating morphology, microstructure, and uniformity. Establishing these correlations will enable predictive links between hydrodynamic performance and coating quality, providing a rational, scalable basis for optimizing SFB design prior to coating deposition.

CFD/DEM↗
DOE OSTI2026

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗
DOE OSTI2023

Stochastic generation of electrolyzer anode catalyst layers

Here, we introduce a stochastic methodology to reproduce the complex pore structure observed in commercial iridium catalyst layers. This method preserves the α pore (pores smaller than 250 nm) and β pore (pores greater than or equal to 250 nm) regions of the catalyst layer. The morphology of the generated materials was validated by comparing the pore size distributions of generated materials against those obtained from commercial materials imaged using x-ray nano computed tomography. We further demonstrate that the pore size distributions of the generated materials are statistically indistinguishable from the imaged catalyst layers, indicating that the stochastic methodology is capable of accurately reproducing catalyst layer morphology. Pore network modelling was conducted on the generated catalyst materials to simulate single-phase permeability, electrical conductivity, and ionic conductivity, and these properties were found to be within experimentally measured ranges for electrolyzer catalyst layers. Additionally, simulations were performed on the generated materials with varying ionomer and iridium catalyst loadings. As the ionomer loading is added, proton conductivity increases exponentially, which demonstrates the importance of optimizing ionomer loading, considering that these effects will be exacerbated in the hydration and temperature conditions of operating electrolyzers. The stochastic material generation method presented in this work is a powerful tool for the development of novel low loading catalyst layers, where the effect of various structural parameters on electrolyzer performance characteristics can be explored.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accelerated fission rate irradiation design, pre-irradiation characterization, and adaptation of conventional PIE methods for U-10Mo and U-17Mo

Metallic U alloys have high U density and thermal conductivity and thus have been explored since the beginning of nuclear power research. Alloys of U with modest amounts of Mo, such as U-10 wt % Mo (U-10Mo), are of particular interest because the γ-U crystal structure in this alloying addition shows prolonged stability in reactor service. Historically, radiation data on U-10Mo fuels were collected in Na fast reactors or lower temperature research reactor conditions, but little is known about irradiation behavior, particularly swelling and creep, at irradiation temperatures between 250 and 500°C. This work discusses the methodology and pre-irradiation characterization results from a U-Mo irradiation campaign performed in the High Flux Isotope Reactor at Oak Ridge National Laboratory. U-10Mo and U-17Mo samples irradiations are being completed at temperatures ranging from 250 to 500°C to three targeted fission densities between 2 × 10 20 and 1.5 × 10 21 fissions per cubic centimeter. Swelling measurement of the specimen sizes studied here required development and assessment of new methods for volume determination before and after irradiation. Laser profilometry and X-ray computation tomography (XCT) were used to provide preirradiation characterization of samples to determine the error and applicability of each to determine swelling following irradiation. These outcomes are contextualized through use of BISON simulations performed to assess the predicted expansion of U-Mo fuels subjected to the irradiation conditions of this work. Use of existing BISON fuel performance models predicted a maximum of 7% swelling under the irradiation conditions of this study. Pre-irradiation characterization revealed the as-cast U-Mo fuel samples were uniformly large-grained fully cubic U crystals with small U-C/N bearing precipitates and pores distributed throughout. Samples were found to contain a bulk porosity between .4 and 3% because of the casting process. Local porosity in areas far from large, interconnected pores was found by Slice-and-View to be under .2%. Nanometer-sized precipitates rich in C and N were identified in all samples, likely because of impurities during the fabrication process. Dendritic bands were also observed throughout the samples. These bands were characterized by variable Mo content that deviated from the overall Mo content by 2–3 wt %. No other microstructural features were correlated to these bands. Mechanical properties were found to be slightly strengthened compared to literature reports of bulk U-Mo fuels due to the nano-scale precipitates throughout the sample.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

SA508 low alloy steel to 316L stainless steel dissimilar metal joint made by powder metallurgy hot isostatic pressing

Joining ferritic SA508 low alloy steel (LAS) and austenitic 316 L stainless steel (SS) via powder metallurgy hot isostatic pressing (PM-HIP) was evaluated as an alternative method to welding. This study investigated the mechanical and microstructural evolutions of the bimetallic interface under different joint designs and heat treatments. Here, the direct joining of dissimilar metal alloys by PM-HIP method resulted in two designs: 1) powder SA508 to solid bar 316 L (P508–B316L) and 2) powder 316 L solid bar SA508 (P316L-B508). In both cases, P508–B316L and P316L-demonstrated satisfactory tensile strength, however, high hardness and severe depreciation in toughness were located on the bimetallic interface. The mechanisms responsible for the detrimental mechanical properties were verified. Large oxides were observed in P508–B316L due to the prior powder boundary (PPB) oxides present in SA508 powder. The intense sensitization occurred from the formation of M 23 C 6 carbides, consequently from the slow cooling after PM-HIP in P316L-B508. Post-HIP heat treatments were explored to reduce the distance of carbide formation; however, the heat treatment could not eliminate the carbides. The experimental results were compared to the diffusion couple simulation as a function of carbide formation with distance. The analysis also showed the high hardness at the bimetallic interface was primarily contributed by solid solution strengthening. In conclusion, the direct joining of P316L-B508 and P508–B316L via PM-HIP was deemed to be unfeasible, and a transitional material is necessary to impede the diffusion of carbon.

316L stainless Steel↗
DOE OSTI2025

A Robust Data-Driven Approach for Mechanical Serial Sectioning

Mechanical serial sectioning (MSS) provides detailed microstructural information across large length scales. By repeatedly removing thin layers of material and imaging the exposed surface, a 3D representation of a specimen’s internal structure can be constructed, enabling failure analysis and feature identification that are otherwise inaccessible via conventional 2D or nondestructive evaluation techniques. Achieving consistent and accurate material removal can be challenging due to system variability, requiring an experienced operator to manually adjust parameters, prolonging data collection times and necessitating post-processing routines to standardize the data. Here, to address these challenges, this paper presents the employment of a one-step model predictive control (MPC) framework tailored to a run-to-run (R2R) controller. The R2R-MPC controller automates the parameter selection process, improving the consistency of material removal through iterative feedback for disturbance rejection and accurate tracking of the target removal rate. Using a data-driven approach, the controller robustly adapts to changing material characteristics. The effectiveness of the R2R-MPC controller is demonstrated through simulation and experimental results and compared to previous data collection procedures.

3D Materials Science↗
DOE OSTI2023

A new efficient grain growth model using a random Gaussian-sampled mode filter

This paper presents the use of a Gaussian neighborhood mode filter for predicting grain growth in a manner similar to the solutions obtained by a Monte Carlo Potts model. This flexible grain growth model can quickly utilize modern, computationally optimized data science strategies on graphics processing units to simulate grain growth up to 100 times faster than the state-of-the-art, publicly available Monte Carlo Potts model. We show that, given the correct neighborhood, the mode filter can replicate normal grain growth in two or three dimensions. In addition, the paper briefly demonstrates the ability to model limited anisotropic in grain boundary energy and mobility. Anisotropic grain boundary energy is modeled by defining a weighted mode filter operation. Anisotropic grain boundary mobility is modeled by scaling and orienting the Gaussian neighborhood in a particular direction.

Anisotropy↗
DOE OSTI2021

The effects of free surfaces on deformation twinning in HCP metals

Deformation twinning is a predominant mode of plastic deformation for hexagonal close packed metals, like Mg and Ti. The heterogenous microstructure and the local stresses associated with twinning play a key role in their mechanical response and fracture. Surface analyses, like electron microscopy, are frequently employed to spatially map microstructure and micromechanical fields in order to study twinning behavior. However, these measurements are inherently influenced by the vicinity of the free surface. In this work, an elasto-visco-plastic fast-Fourier-transform (EVP-FFT) polycrystal modeling approach is employed to investigate the effects of free surfaces on twin development before and after loading. We compare calculated micromechanical fields on free surfaces with those calculated inside the bulk and, in some cases, experimental surface measurements. The results indicate that the creation of free surfaces can promote twin propagation and growth and can influence twin morphology by causing a twin lamella to become larger, more blunted and irregular. The structure along the twin boundaries are also affected, due to the higher driving stresses that extend prismatic-basal and basal-prismatic facets. Furthermore, free surfaces invoke different slip activities in the twin and the surrounding parent crystal by enhancing basal, prismatic and pyramidal slip in some localized regions, while reducing slip in others. We demonstrate that the simulated free-surface effects lead to better qualitative and quantitative agreement with experimental measurements from scanning electron microscopy and digital image correlation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Interpreting experimental measurements of helium bubbles using stochastic cluster dynamics models of heterogeneous nucleation and growth in irradiated ferritic alloys

Among a number of other advantageous features, ferritic/martensitic steels are being considered as fusion reactor structural materials due to their low intrinsic swelling under irradiation. However, under high-energy neutron irradiation, He produced through (n, α) reactions stabilizes vacancy clusters, which then act as seeds for bubble formation and growth, negating the intrinsic swelling resistance of these alloys. Standard models of irradiation damage accumulation and microstructural evolution consider homogeneous nucleation as the basis for bubble formation and growth. However, this generally does not explain the large bubble densities and sizes observed experimentally under a number of different conditions. Here, we propose a more realistic physical model of bubble nucleation, formation, and growth designed to capture recent experimental measurements of He-bubble formation and evolution during co-implantation of He and Fe ions in model ferritic alloys. We find that experimental results are explained only when the following three features are considered simultaneously: (i) heterogeneous nucleation of He-vacancy bubbles at defect sinks (e.g., dislocations, grain boundaries, and second-phase precipitates), (ii) nucleation and growth of bubbles via the ‘trap mutation’ mechanism (i.e., spontaneous production of Frenkel pairs due to absorption of He atoms), and (iii) transition from a growth-limited, He-stabilized bubble regime to a ‘runaway’ void-type growth scenario. The model is implemented into a stochastic cluster dynamics framework that takes advantage of cluster size grouping methods to accelerate the simulations, allowing us to reach 10 dpa of simulated irradiated dose, and to capture cluster sizes in excess of 20 nm. Further, a careful extrapolation exercise conducted assuming classical nucleation theory leads to excellent agreement with the experimental measurements at 50 dpa of irradiation.

36 MATERIALS SCIENCE↗