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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Coarse-grained simulation of colloidal self-assembly, cation exchange, and rheology in Na/Ca smectite clay gels

Knowledge Gap: The aggregation of clay minerals—layered silicate nanoparticles—strongly impacts fluid flow, solute migration, and solid mechanics in soils, sediments, and sedimentary rocks. Experimental and computational characterization of clay aggregation is inhibited by the delicate water-mediated nature of clay colloidal interactions and by the range of spatial scales involved, from 1 nm thick platelets to flocs with dimensions up to micrometers or more. Simulations: Using a new coarse-grained molecular dynamics (CGMD) approach, we predicted the microstructure, dynamics, and rheology of hydrated smectite (more precisely, montmorillonite) clay gels containing up to 2,000 clay platelets on length scales up to 0.1 μm. Further, simulations investigated the impact of simulation time, platelet diameters (6 to 25nm), and the ratio of Na to Ca exchangeable cations on the assembly of tactoids (i.e., stacks of parallel clay platelets) and larger aggregates (i.e., assemblages of tactoids). We analyzed structural features including tactoid size and size distribution, basal spacing, counterion distribution in the electrical double layer, clay association modes, and the rheological properties of smectite gels. Findings: Our results demonstrate new potential to characterize and understand clay aggregation in dilute suspensions and gels on a scale of thousands of particles with explicit representation of counterion clouds and with accuracy approaching that of all-atom molecular dynamics (MD) simulations. For example, our simulations predict the strong impact of Na/Ca ratio on clay tactoid formation and the shear-thinning rheology of clay gels.

42 ENGINEERING↗

DOE OSTI2025

Unravelling Microstructure Selection in an Additively Manufactured Eutectic High‐Entropy Alloy

High-entropy alloys (HEAs) are promising candidates for advanced structural applications due to their excellent mechanical properties. Additive manufacturing (AM), with its rapid solidification conditions, enables the creation of unique nonequilibrium microstructures. To fully leverage the synergy between AM and HEAs, understanding how processing affects structure and properties is essential. Here, how solidification rate influences microstructure evolution and phase transformation pathway in laser additively manufactured AlCrFe2Ni2 eutectic HEAs is investigated. By increasing the laser scan speed and hence the solidification rate, distinct solidification modes evolving from coupled eutectic to anomalous eutectic and eventually to single-phase solidification are revealed. These transitions result in distinct microstructures and a wide range of mechanical properties. Thermodynamic modeling and molecular dynamics simulations reveal that low cooling rates allow for sufficient atomic diffusion and phase separation, facilitating coupled eutectic growth. In contrast, rapid cooling suppresses diffusion and destabilizes the solid–liquid interface, promoting anomalous or single-phase solidification. This integrated experimental and computational approach provides a multiscale understanding of solidification mechanisms in HEAs and underscores how kinetic effects can over-ride thermodynamic predictions under nonequilibrium conditions. Furthermore, these results demonstrate that AM can serve as a powerful tool to design HEAs with tailored microstructures and properties.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Preliminary Results on Process Modeling Tools for Determining Variability in Additively Manufactured Stainless Steel 316 Parts

The Advanced Materials and Manufacturing Technologies program aims to accelerate the development, qualification, demonstration, and deployment of advanced materials and manufacturing technologies to enable reliable and economical nuclear energy. However, the distinct characteristics of additive manufacturing (AM) materials, stemming from their unique processing history, microstructure, and properties, pose significant challenges for the qualification and certification of nuclear components. These challenges primarily arise from component-scale variations in microstructure and properties influenced by local process conditions and geometry, which affect thermal history, melt pool dynamics, and microstructure evolution. Computational modeling tools can play a crucial role in predicting and controlling this variability. This report presents preliminary results on process modeling tools designed to predict microstructure variability in additively manufactured stainless steel 316 parts. It details the software packages and physical modeling approaches employed to simulate an AM component within an automated process modeling workflow. Initial results are demonstrated through comparisons between predicted microstructures and experimental measurements across various representative processing conditions. The report concludes by discussing the challenges inherent in process modeling of AM components and outlines a plan for future development needs.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

DOE OSTI2025

Exascale granular microstructure reconstruction in 3D volumes of arbitrary geometries with generative learning

Reconstructing 3D granular microstructures within volumes of arbitrary geometries from limited 2D image data is crucial for predicting the material properties, as well as performances of structural components accounting for material microstructural effects. We present a novel generative learning framework that enables exascale reconstruction of granular microstructures within complex 3D geometric volumes. Building upon existing transfer learning techniques using pre-trained convolutional neural networks (CNN), we introduce several key innovations to overcome the difficulties inherent in arbitrary geometries. Our framework incorporates periodic boundary conditions using circular padding techniques, ensuring continuity and representativeness of the reconstructed microstructures. We also introduce a novel seamless transition reconstruction (STR) method that creates statistically equivalent transition zones to integrate multiple pre-existing 3D microstructure volumes. Based on STR, we propose a cost-effective strategy for reconstructing microstructures within complex geometric volumes, minimizing computational waste. Validation through numerical experiments using kinetic Monte Carlo simulations demonstrates accurate reproduction of grain statistics, including grain size distributions and morphology. A case study involving the reconstruction of a 4-blade propeller microstructure illustrates the method’s capability to efficiently handle complex geometries. In conclusion, the proposed framework significantly reduces computational demands while maintaining high reconstruction quality, paving the way for scalable microstructure reconstruction in materials design and analysis.

36 MATERIALS SCIENCE↗

NASA NTRSDate not supplied

Simulating Micrometeoroid Bombardment on Mercury: Identifying New Space Weathering Features in the Laboratory

Airless bodies are continually exposed to the harsh interplanetary space environment, causing the alteration of their surfaces through a process known as space weathering (SW). SW is dominated by solar wind irradiation and micrometeoroid bombardment, which together alter the spectral, microstructural, and chemical characteristics of grains on the surface of airless bodies. The effects of SW are variable and depend on the heliocentric distance and the initial composition of the planetary surface, among other factors. On the Moon and S- type asteroids, it is well established that SW darkens and reddens the surface of these bodies and attenuates absorption bands across the Vis-NIR wavelengths. These spectral effects are largely the result of the production of metallic Fe nanoparticles (npFe). However, Mercury’s SW environment is unique compared to the Moon and S-type asteroids. Mercury is a geochemical endmember, with a surface composition low in Fe (<2 wt.%) and enriched in volatile components, including regions of the surface hypothesized to contain up to 4 wt.% graphite in the low reflectance material (LRM). In fact, the presence of graphite on Mercury has been hypothesized to act as a reducing agent for silicates during SW to produce Si-bearing Fe-rich metal. In addition, its location in the solar system exposes Mercury to an extreme SW environment, with the surface of the planet experiencing an intense solar wind flux and higher flux and velocity of micrometeoroid impactors compared to the Moon and S-type asteroids. The effects that such harsh SW has on materials with the unique compositional characteristics of Mercury are not well constrained. To better understand SW at Mercury, we must investigate these processes in the laboratory. Here, we present the results of our analyses of the spectral, microstructural and chemical characteristics of Mercury analog samples irradiated by pulsed laser to simulate the short duration, high temperature events associated with micrometeoroid impacts.

Nick Bott↗

NASA NTRS1990

Effective conductivity of hard-sphere dispersions

A simulation technique is used to compute zeta(2), the key three-point microstructural parameter which arises in three-point conductivity bounds on the effective conductivity sigma(e) of isotropic two-phase composites. Using the new algorithm, the zeta(2) parameter is computed for dispersions of hard spheres by Monte Carlo simulation, and the results are compared with those of several analytical calculations of zeta(2). The simulation data on zeta(2) are then used to compute three-point bounds on the conductivity sigma(e) and Torquato's (1985) approximation for sigma(e).

Miller, C. A.↗

NASA NTRS1986

A study of spectrum fatigue crack propagation in two aluminum alloys. I - Spectrum simplification. II - Influence of microstructures

An investigation of the fatigue crack propagation FCP behavior of two aluminum alloys is performed to simulate spectrum loading conditions found at critical locations in high performance fighter aircraft. Negative loads are shown to be eliminated for the tension-compression spectrum for low to intermediate maximum stress intensities, and load interactions are found to be more significant at higher stress intensities and with more plasticity at the crack tip. In the second part, the influence of microstructural features including grain size, inclusions, and dispersoids on constant amplitude and spectrum crack growth behavior in aluminum alloys is studied. At low stress intensities the I/M alloy demonstrated better FCP resistance than the P/M 7091 alloy for both constant amplitude and spectrum testing, and the inhomogeneous planar slip and large grain size of 7050 limit dislocation interactions, thereby improving FCP performance.

Telesman, J.↗

DOE OSTI2026

Dynamic strength of iron under pressure-temperature conditions of Earth’s inner core

Iron (Fe) is a primary constituent of terrestrial planetary cores, yet its rheological properties under extreme conditions remain uncertain. Here we present direct measurements of Fe strength at 310-430 GPa pressures and 3700-5800 K temperatures, obtained using Rayleigh-Taylor (RT) instability experiments at the National Ignition Facility. Single-crystal α-Fe samples with [001] and [111] orientations are shock-ramp compressed past the α-ε transition along paths approaching Earth’s inner core conditions. We find that ε-Fe derived from [001] α-Fe is consistently stronger (11-20 GPa) than that from [111] α-Fe (8-18 GPa), contrary to the trend at ambient conditions. Large-scale molecular dynamics simulations reproduce this atypical strength anisotropy and attribute it to microstructural evolution during the phase transition and subsequent ε-phase plasticity. Ripple growth analysis further constrains viscosities of 100-170 Pa·s under the driven conditions. Our results provide experimental benchmarks for Fe rheology at inner-core conditions, with implications for seismic anisotropy and the geodynamo.

Condensed-matter physics↗

NASA NTRS1985

Liquid Phase Miscibility Gap Materials

The manner in which the microstructural features of liquid-phase miscibility gap alloys develop was determined. This will allow control of the microstructures and the resultant properties of these alloys. The long-duration low gravity afforded by the shuttle will allow experiments supporting this research to be conducted with minimal interference from buoyancy effects and gravitationally driven convection currents. Ground base studies were conducted on Al-In, Cu-Pb, and Te-Tl alloys to determine the effect of cooling rate, composition, and interfacial energies on the phase separation and solidification processes that influence the development of microstructure in these alloys. Isothermal and directional cooling experiments and simulations are conducted. The ground based activities are used as a technological base from which flight experiments formulated and to which these flight experiments are compared.

Gelles, S. H.↗

NASA NTRSDate not supplied

Overview of Ablative TPS Modeling at NASA Ames

Over the past decade, NASA has invested in efforts to build predictive thermal protection system (TPS) material models from the micro-scale to the macro-scale. To complement the mission design cycle process and reduce the need for extensive testing, NASA is developing modeling and simulation tools that enable characterizing material properties and response to hot plasma experienced during atmospheric entry. Traditional material response and ablation modeling tools, such as the heritage code FIAT, and its multidimensional siblings, TITAN and 3dFIAT, are being complemented with newly developed software such as Icarus and PATO. Both of these programs are three-dimensional, finite-volume solvers that use unstructured meshes and 21st century programming paradigms to allow for efficient parallel simulations. FIAT and Icarus are also used for TPS sizing purposes. Today, these traditional tools are being supplemented with computational materials models at the atomistic level. The scales of interest range from computational chemistry (Density Functional Theory [DFT]), to atomistic simulations (Molecular Dynamics [MD]), to the microscale with the Porous Microstructure Analysis (PuMA) software that was recently awarded the 2022 NASA Software of the Year award. Finally, thermo-structural modeling is also of interest to the TPS Materials branch and done using commercial tools such as MSC MARC, MENTAT, NASTRAN and PATRAN. The present talk will also link the use of these computational tools to current NASA missions and projects associated with challenging and complex vehicles entries/reentries.

materials modeling↗

DOE OSTI2024

Adaptively remeshed multiphysical modeling of resistance forge welding with experimental validation of residual stress fields and measurement processes

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Wrought FeCrAl alloy (C26M) cladding behavior and burst under simulated loss-of-coolant accident conditions

Cladding burst experiments for FeCrAl cladding were performed in the Severe Accident Test Station facility at Oak Ridge National Laboratory. These experiments were simulated using the BISON fuel performance code to better understand the cladding plastic behavior and failure under simulated loss-of-coolant accident conditions. 3D cladding surface boundary conditions were generated using composite axial and azimuthal profiles from experiment thermocouple data. To improve the simulation analysis capabilities in BISON for cladding burst behavior, new thermal creep, plasticity, and failure stress models specific to C26M, a wrought FeCrAl alloy, were developed and implemented. Initial cladding burst results indicated a general underprediction in the failure temperature of the six cladding burst simulations versus the observed failure temperatures. Close investigation of the experiment timing versus the underlying tensile test data revealed that, compared with the tensile specimens, the cladding tubes did not experience the same long holding time at high temperatures. New tensile tests were performed at high temperatures using a temperature ramp similar to the simulated loss-of-coolant accident experiments. These new tensile curves showed an approximately 80% increase in the ultimate tensile strength of the C26M alloy, indicating that a holding time of 10 min at 700 °C and 800 °C allows annealing to change the material microstructure. Using the updated tensile properties, the burst temperatures and stresses from the simulations showed remarkable agreement with the experimental results. This study was then extended by varying the initial pressure to highlight the burst temperature difference between standard Zircaloy-4 and C26M cladding under equivalent conditions. The results show that C26M has a burst temperature that is approximately 70–130 K greater than that of Zircaloy-4. In conclusion, these modeling predictions can be further improved by collecting high-temperature tensile data for C26M beyond the temperature ranges used in this work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Uncovering grain and subgrain microstructure at the scale of additive manufacturing melt tracks with a scalable cellular automaton solidification model

Metal additive manufacturing, characterized by rapid solidification, yields refined grains with a distinctive cellular subgrain microstructure that plays a pivotal role in determining material properties. Due to the significant computational expense demanded to simulate the required physics with submicron spatial resolution, their numerical simulations have been limited to proof-of-concept studies to either 2D or small subregions of a melt pool. In this study, an open-source, scalable, solidification code, muMatScale, based on the cellular automaton method, has been developed to predict the grain and the underlying subgrain microstructure over an entire melt pool. The model incorporates flexible parallelization schemes, utilizing MPI and OpenMP GPU Offloading, in addition to appropriate multi-physics specific to non-equilibrium rapid solidification in AM. The impact of nucleation parameters on grain microstructures was investigated with a focus on grain size variations and morphology transitions. With selected nucleation parameters, the simulation predicted the grain size, subgrain morphology, crystallographic orientation, and microsegregation aligned with experimental measurements. The model demonstrates that epitaxial grain growth is a dominant factor at the melt pool boundary, influencing grain size variation under different grain sizes in the build plate while maintaining consistent primary dendrite arm spacing under identical thermal conditions. Here, the highly efficient numerical model enables large-scale simulations with a spatial resolution of 100 nm or less, unveiling unprecedented insights into thermal and solutal diffusion driven grain growth, and the subgrains with microsegregation within grains in 3D across scales. muMatScale will enable the linking of submicron length-scale microstructure to part-level material behavior by investigating fundamental solidification problems at the intercellular scale in many-track and many-layer builds.

36 MATERIALS SCIENCE↗

NASA NTRS2004

High-Pressure Minerals in Meteorites: Constraints on Shock Conditions and Duration

The objective of this research was to better understand the conditions and duration of shock metamorphism in meteorites through microstructural and microanalytical characterization of high-pressure minerals. A) Continue to investigate the mineralogy and microstructures of melt-veins in a suite of chondritic samples ranging from shock grades S3 through S6 to determine how the mineral assemblages that crystallize at high-pressure and are related to shock grade. B) Investigate the chemical, mineralogical, and microstructural heterogeneities that occur across melt veins to interpret crystallization histories. C) Use static high-pressure experiments to simulate crystallization of melt veins for mineralogical and textural comparisons with the melt veins of naturally shocked samples. D) Characterize the compositions and defect microstructures of polycrystalline ringwoodite, wadsleyite, majorite, (Mg,Fe)Si03-ilmenite and (Mg,Fe)SiO3-perovskite in S6 samples to understand the mechanisms of phase transformations that occur during shock. These results will combined with kinetic data to constrain the time scales of kinetic processes. E) Investigate the transformations of metastable high-pressure minerals back to low- pressure forms to constrain post-shock temperatures and estimates of the peak shock pressure. Of these objectives, we have obtained publishable data on A, B and D. I am currently doing difficult high-pressure melting and quench experiments on an L chondrite known as Mbale. These experiments will provide additional constraints on the mineral assemblages that are produced during rapid quench of an L chondrite at pressures of 16 to 25 GPa. Results from published or nearly published research is presented below. Lists of theses, dissertations and publications are given below.

Sharp, Thomas G.↗

NASA NTRSDate not supplied

Uncertainty Quantification and Sensitivity Analysis in Process-Structure-Property Simulations for Laser Powder Bed Fusion Additive Manufacturing

Process variations and process-induced defects like porosity cause significant uncertainty in the microstructure and mechanical behavior of additively manufactured metals. Establishing process-structure-property (PSP) relationships and quantifying uncertainty using experiments alone is costly, especially for structural applications where mechanical allowables must be established for qualification and certification. This work presents a PSP simulation framework for laser powder bed fusion with a focus on uncertainty quantification through probabilistic calibration and multi-fidelity uncertainty propagation. Motivated by phenomenological input parameters related to grain nucleation and growth that are difficult to characterize, a global sensitivity analysis (GSA) is completed. Through GSA, the most important input parameters are identified based on their influence on the statistical distributions of microstructural metrics that influence mechanical behavior, including grain size, morphology, and crystallographic texture. The results provide insight on what experiments are necessary to quantify and control PSP uncertainties, particularly those associated with the more challenging input parameters.

additive manufacturing↗

DOE OSTI2024

Location-Specific Microstructure Characterization Within AM Bench 2022 Nickel Alloy 718 3D Builds

Abstract The Additive Manufacturing Benchmark Test Series (AM Bench) is a broad effort to produce rigorous measurement datasets for validating AM computer simulations across the range of processing, structure, and properties, for many additive manufacturing (AM) build methods and material classes. Here, the microstructures of nickel alloy 718 AM Bench 2022 test artifacts produced using laser-based powder bed fusion (PBF-LB), in both as-built and fully heat-treated conditions, are examined. Cross sections are primarily characterized using large area scanning electron microscopy (SEM) electron backscatter diffraction (EBSD) and example analyses of the crystallographic textures are described. These data are part of a large set of in situ and ex situ measurements from both three-dimensional builds and laser tracks on bare plates. All the measurement data are available online with download links at www.nist.gov/ambench .

Levine, L. E. (ORCID:0000000334484229)↗

DOE OSTI2025

Optimizing Bendability and Hardness of Age-Hardenable Aluminum Sheets through Local Thermo-Mechanical Processing

This paper introduces a novel thermo-mechanical process to modify the local mechanical properties of 6xxx aluminum (Al-Mg-Si-Cu) alloy sheets. In this process, two pairs of rollers travel along the length of a sheet while locally bending and unbending it so that the final shape and thickness of the sheet remains unchanged. Room temperature bending/unbending (B/U) produces a deformation gradient through thickness and local hardening in both T4 and T6 sheets (with 41% and 18% greater Vickers hardness, respectively). High temperature bending/unbending performed at ~ 500°C via induction heating, produces significant improvements in local formability as measured by bend testing. Formability levels equivalent to the as-received T4 temper are achieved within the processed zones of T6 sheets with bend angles of ~ 150°, without disturbing the T6 temper in the remainder of the sheet. Finite element simulations and characterization explained the mechanical property improvements and the heterogenous microstructure with weakened Cube texture. Due to its similarity with the roller hemming process, and its compact dimensions, the apparatus developed for B/U has the potential for seamless integration with roller hemming robots to selectively process high strength aluminum sheet materials within mass-production settings.

Efe, Mert [BATTELLE (PACIFIC NW LAB)]↗