Search NASASearch

SEARCH · Search NASA

Results for “Microstructure simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 55 records · Page 3

Three-dimensional in Situ Observations of Polycrystalline Microstructure Evolution During Directional Solidification of Transparent Alloys Aboard the ISS and Quantitative Comparison with Numerical Modeling

Cellular/dendritic microstructures formed during solidification have a crucial influence on the mechanical properties of a wide range of structural alloys. By minimizing the amount of gravity-induced convection in the liquid, directional solidification experiments using transparent organic alloys conducted in the DECLIC-DSI onboard the International Space Station have provided unique 3D in situ observations of the spatiotemporal evolution of the solid-liquid interface during the formation of cellular and dendritic microstructures under purely diffusive growth conditions in polycrystalline samples containing several grains with a small misorientation with respect to the temperature gradient. Those observations have made it possible to perform benchmark quantitative comparisons with the predictions of state-of-the-art phase-field simulations of microstructure formation in 3D on experimentally relevant length and time scales. This talk will report quantitative comparisons between microgravity experiments and phase-field simulations in succinonitrile-camphor alloys of two different compositions that shed new light on the role of subgrain boundaries in the spatiotemporal evolution of the primary cellular spacing and the selection of dendritic array structures.

Kaihua Ji

Physics-Based Modeling and Simulation of Self-Reacting Friction Stir Welding Using Computational Fluid Dynamics

A physics-based model was developed to simulate the behavior of material in a self-reacting friction stir welding (SR-FSW) process for the joining of metals. This steady-state model builds upon fundamental computational fluid dynamic (CFD) principles within Ansys Fluent to solve the discretized equations. The effective viscosity is calculated using a viscoplastic model using a Sheppard-Wright formulation of flow stress. Numerous advancements have been made in the incorporated physics including (1) temperature-dependent material properties; (2) locally adaptable flow and thermal boundary conditions; and (3) adapting material properties in nugget in response to microstructural changes. Simulation strategies to accelerate computation and improve numerical stability include adapting the mesh refinement and solver relaxation factors during simulation. The result is a highly robust and computationally efficient model capable of providing the material flow and temperature history across the domain. As material history determines the local microstructure and ultimately weld strength, an accurate and detailed physics-based model has the potential to accelerate SR-FSW process development. The model is highly adaptable to changes in process parameters, tool design, or alloy.

Process Modeling

Age hardening response of Cu-30Ni alloys: The role of Si microalloying additions

Cu-30Ni alloys offer high strength and excellent corrosion resistance for maritime applications. While primarily a solid solution system, industrial alloys typically include microalloying additions of Fe, Mn, Nb and/or Si to enable age hardening. However, an understanding of the microstructural evolution during age hardening remains incomplete. Here, we conduct systematic ageing treatments and report the resulting hardness and microstructures. The Cu-30Ni alloy with Si aged at 650 °C for 6 h demonstrates significantly enhanced Vickers hardness, reaching ∼201 HV 1 compared to 103 HV 1 for the as-homogenised sample. Si-rich clusters and γ′ precipitates are identified, with their composition, size, and volume fraction quantified to determine their strengthening contributions. First-principles atomistic simulations elucidate the underlying formation mechanisms of these clusters and precipitates, highlighting the critical role of Si in driving their nucleation and growth. This study advances the design of high-strength, Cu alloys with the potential for corrosion resistance in demanding maritime environments.

Age hardening

Under the microscope: Reduced activation ferritic martensitic steel Eurofer-97 following ion-Irradiation and high-temperature high-pressure water exposure

Here, this study is designed to characterise the microstructural behaviour of Eurofer-97 steel under ion irradiation and subsequent exposure to high-temperature high-pressure (HTHP) water. Eurofer-97, a ferritic-martensitic steel, has been developed to withstand the conditions of fusion reactors in the locations in contact with coolant with an elevated level of neutron flux, such as the breeder-wall blanket. The material has been studied after self-ion irradiation (using Fe ions) simulating the microstructural effects of neutron irradiation limited to the subsurface layer. The corrosion properties of the Eurofer-97 steel were studied by exposure to HTHP water up to 331 °C. Advanced microstructural characterisation using scanning, transmission electron and focused ion beam microscopy was performed on the as-received microstructure and after ion irradiation. This was then characterised after exposure for 240 h in high-temperature water. Eurofer-97 had a dense, columnar Cr-rich inner oxide, followed by a Fe-rich outer oxide layer. In the irradiated condition the grain structure and oxide itself was less ordered. No appreciable difference in oxide thickness was identified between the irradiated and unirradiated specimens after this short exposure time.

36 MATERIALS SCIENCE

In Situ Observation of Growth Dynamics in DECLIC Directional Solidification Insert Onboard ISS: DSI-R Flight Campaign

The study of solidification microstructure formation is of utmost importance for materials design and processing, as solid-liquid interface patterns largely govern mechanical and physical properties. Pattern selection occurs under dynamic conditions of growth in which the initial morphological instability evolves nonlinearly and undergoes a reorganization process. The dynamic and nonlinear nature of this instability renders in situ observation of the interface an invaluable tool to gain knowledge on the time-evolution of the interface pattern. Transparent organic analogs, which solidify like metallic alloys, allow direct visualization of interface dynamics. Extensive ground-based studies of both metallic and organic bulk samples have established the presence of significant convection during solidification processes that alters the formation of interfacial microstructures. A reduced-gravity environment is therefore mandatory for fluid flow elimination in bulk samples. In the framework of the CNES project MISOL3D (MIcrostrutures de SOLidification 3D) and the NASA projects DSIP (Dynamical Selection of 3D Interface Patterns), SPADES (SPAtiotemporal Evolution of three-dimensional DEndritic array Structures) and CAMUS (ComputAtional Studies of MicrostrUcture Formation During Alloy Solidification in Microgravity), we participated in the development of the Directional Solidification Insert (DSI) of the DEvice for the study of Critical Liquids and Crystallization (DECLIC). The DECLIC-DSI is dedicated to in situ and real-time characterization of solid-liquid interface patterns during directional solidification of transparent alloys in diffusive transport regime. Between April 2010 and March 2011, the first ISS campaign (DSI) explored the entire range of microstructures resulting in unprecedented observations. A second campaign (DSI-R), performed between October 2017 and December 2018, in which the insert contained an alloy of higher solute concentration, allowed to complete the benchmark database. The increase of solute concentration resulted in well-developed dendritic patterns at lower velocities (lower interface curvature and larger tip radius). The microstructure resulting from dendritic growth is dominant in metallurgy so that it is fundamental to understand the mechanisms of its formation. The main aims of this experimental campaign are to understand: the mechanisms of the cell to dendrite transition, the fundamental mechanisms of sidebranching formation, the dependence of dendrite tip shapes on growth conditions, the interaction of primary array and secondary sidebranches, and the influence of subgrain boundaries on the spatiotemporal organization of the array structure. Preparation, analysis and interpretation of the experiments performed onboard ISS are considerably enhanced by experiments performed on ground using thin-samples (Pr. Trivedi’s group) and phase-field simulations of microstructure formation in a diffuse growth regime (Pr. Karma’s group). In this summary, we will present an initial assessment of the results obtained during the DSI-R campaign. Nomenclature Vp: pulling velocity G: thermal gradient L: solidified length rtip: tip curvature radius Acronyms/Abbreviations CADMOS: Centre d'Aide au Développement des activités en Micro-pesanteur et des Opérations Spatiales CNES: Centre National d’Études Spatiales DECLIC: Device for the study of Critical LIquids and Crystallization DSI: Directional Solidification Insert DSI-R: Directional Solidification Insert - Refurbish ISS: International Space Station NASA: National Aeronautics and Space Administration PF: phase-field SCN: Succinonitrile 3D: three-dimensional

Mota, F.L.

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Text Mining for Process–Structure–Properties Relationships in Metals

With the advent of large language models (LLMs), the vast unstructured text within millions of academic papers is increasingly accessible for materials discovery—although significant challenges remain. While LLMs offer promising few- and zero-shot learning capabilities, particularly valuable in the materials domain where expert annotations are scarce, general-purpose LLMs often fail to address key materials-specific queries without further adaptation. To bridge this gap, fine-tuning LLMs on human-labeled data is essential for effective structured knowledge extraction (Liu in The Importance of Human-Labeled Data in the Era of LLMs, 2023). Here, in this study, we introduce a novel annotation schema designed to extract generic process–structure–properties relationships from scientific literature. We demonstrate the utility of this approach using a dataset of 128 abstracts, with annotations drawn from two distinct domains: high-temperature materials (Domain I) and uncertainty quantification in simulating materials microstructure (Domain II). Initially, we developed a conditional random field (CRF) model based on MatBERT—a domain-specific BERT variant—and evaluated its performance on Domain I. Subsequently, we compared this model with a fine-tuned LLM (GPT-4o from OpenAI) under identical conditions. Our results indicate that fine-tuning LLMs can significantly improve entity extraction performance over the BERT-CRF baseline on Domain I. However, when additional examples from Domain II were incorporated, the performance of the BERT-CRF model became comparable to that of the GPT-4o model. These findings underscore the potential of our schema for structured knowledge extraction and highlight the complementary strengths of both modeling approaches.

Materials science

Computer Simulation of Spatial Arrangement and Connectivity of Particles in Three-Dimensional Microstructure: Application to Model Electrical Conductivity of Polymer Matrix Composite

Computer simulation is a powerful tool for analyzing the geometry of three-dimensional microstructure. A computer simulation model is developed to represent the three-dimensional microstructure of a two-phase particulate composite where particles may be in contact with one another but do not overlap significantly. The model is used to quantify the "connectedness" of the particulate phase of a polymer matrix composite containing hollow carbon particles in a dielectric polymer resin matrix. The simulations are utilized to estimate the morphological percolation volume fraction for electrical conduction, and the effective volume fraction of the particles that actually take part in the electrical conduction. The calculated values of the effective volume fraction are used as an input for a self-consistent physical model for electrical conductivity. The predicted values of electrical conductivity are in very good agreement with the corresponding experimental data on a series of specimens having different particulate volume fraction.

Louis, P.

GrainPaint: A multi-scale diffusion-based generative model for microstructure reconstruction of large-scale objects

Simulation-based approaches to microstructure generation can suffer from a variety of limitations, such as high memory usage, long computational times, and difficulties in generating complex geometries. Generative machine learning models present a way around these issues, but they have previously been limited by the fixed size of their generation area. Here, we present a new microstructure generation methodology leveraging advances in inpainting using denoising diffusion models to overcome this generation area limitation. We show that microstructures generated with the presented methodology are statistically similar to grain structures generated with a kinetic Monte Carlo simulator, SPPARKS.

36 MATERIALS SCIENCE

Multiscale Analysis of Structurally-Graded Microstructures Using Molecular Dynamics, Discrete Dislocation Dynamics and Continuum Crystal Plasticity

A multiscale modeling methodology is developed for structurally-graded material microstructures. Molecular dynamic (MD) simulations are performed at the nanoscale to determine fundamental failure mechanisms and quantify material constitutive parameters. These parameters are used to calibrate material processes at the mesoscale using discrete dislocation dynamics (DD). Different grain boundary interactions with dislocations are analyzed using DD to predict grain-size dependent stress-strain behavior. These relationships are mapped into crystal plasticity (CP) parameters to develop a computationally efficient finite element-based DD/CP model for continuum-level simulations and complete the multiscale analysis by predicting the behavior of macroscopic physical specimens. The present analysis is focused on simulating the behavior of a graded microstructure in which grain sizes are on the order of nanometers in the exterior region and transition to larger, multi-micron size in the interior domain. This microstructural configuration has been shown to offer improved mechanical properties over homogeneous coarse-grained materials by increasing yield stress while maintaining ductility. Various mesoscopic polycrystal models of structurally-graded microstructures are generated, analyzed and used as a benchmark for comparison between multiscale DD/CP model and DD predictions. A final series of simulations utilize the DD/CP analysis method exclusively to study macroscopic models that cannot be analyzed by MD or DD methods alone due to the model size.

Saether, Erik

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM

Additive Manufactured Composite Phase-Change Material for Thermal Energy Storage Applications

Phase-change materials play a critical role in industrial energy storage applications to drive efficiency improvements, thermal energy management, and carbon emissions reductions. Recently, it has been shown that rapid solidification of alloys with metastable immiscibility in the liquid phase has the potential to form unique microstructures in which a low-melting phase is uniformly distributed in a high-melting matrix. This feature can be exploited using additive manufacturing to produce components with complex geometries containing such unique phase-change microstructures. Phase-field simulations utilizing high-performance computing were used to provide a detailed description of the evolution of the active phase during service in terms of their morphology and composition in different polycrystalline matrix grain morphologies that are typically produced during additive manufacturing. Phase field simulations were performed using, MEUMAPPS-SL (Microstructure Evolution Using Massively Parallel Phase-field Simulations – Solid Liquid) code that was developed in-house by the Oak Ridge National Laboratory. The simulations utilized the capabilities of the Kestrel supercomputer at the National Renewable Energy Laboratory. The simulation results were compared with experimental results generated at Siemens Energy, Inc. The results indicate that the kinetics of liquid spreading along grain boundaries is largely determined by the mobility of the triple line along the intersection of the grain boundary liquid and the grain boundary plane.

25 ENERGY STORAGE

Direct Ab Initio Simulation of the Synthesis of BaZrO 3 and the Microstructure Impacts on Proton Transport

Controlling and predicting the processing-structure-performance relationship in functional materials is a grand challenge in materials science, with important implications for a wide range of emerging applications; a high fidelity understanding of the performance impact of microstructures formed under synthesis conditions is required to develop advanced materials, such as solid-state fuel cells and electrolyzers. Using the ceramic BaZrO 3 as a case study, we directly simulate the synthesis and investigate how proton transport is dictated by microstructures. We develop a framework that couples density functional theory (DFT), machine-learning interatomic potential (MLIP) driven molecular dynamics, and grand canonical Monte Carlo to perform large-scale, microstructure-resolved, atomistic simulations of proton transport in experimentally representative polycrystalline structures. Our fully ab initio approach, using a MLIP as a proxy for DFT, allows us to quantify the competition between two distinct diffusion mechanisms: one associated with grain-boundary regions and another within grains. When the impacts of grain boundaries are taken into account, proton transport exhibits substantial deviation from the bulk oxide limit. This addresses long-standing discrepancies between theory and experiments. Our integrated approach provides atomistic insight into microstructure-dependent proton pathways in BaZrO 3 and establishes a general protocol for predicting processing-structure-performance relationships.

organic

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material

Particle-Based Fiber Models of Woven Materials for Earth Entry Thermal Protection

Applications requiring materials with layer-to-layer strength, from basketry to thermal protection systems, use interlaced, three-dimensional woven materials. NASA is developing and deploying woven material heat shields for missions, including Artemis I (3D-MAT for compression pads) and potentially Mars Sample Return - Earth Entry System. These materials are complex, hierarchical and must protect from extreme environments and phenomena, such as deformation, impact and high-enthalpy heating. Woven material performance depends on microstructure, damage and weave geometry. Therefore, fiber-specific models are needed to simulate fiber contacts within the weave hierarchical geometry (fiber to tow to yarn to weave) and the inherent directionality of fibers. Explicit fiber models can simulate how weave microstructure evolution affects thermal and mechanical properties. We parameterize a discrete element bonded particle models (DEM-BPM) of fibers to capture thermal and mechanical behavior within and between fibers, with bonded and contact forces, respectively. We study the proportion of heat transfer and stress via the contact network, fiber bonds and the weave geometry, for example, with respect to yarn warp-weft identity (whether it interlaces weave layers). Our results demonstrate the importance of explicit fiber modeling for connecting microstructure with thermal and mechanical properties.

woven material

Simulation and analysis of small angle scattering (SAS) patterns of Ni-based superalloy microstructures generated by a phase-field model

This paper investigates the relationship between microstructural features and small-angle scattering (SAS) patterns in Ni-based superalloys using a combined phase-field and SAS simulation approach coupled with microstructure analyses. The simulated SAS patterns accurately capture key experimental observations previously reported in the literature, including the time-dependent transition from circular to square-shaped precipitates and the development of anisotropic SAS patterns. Importantly, our analysis reveals the correlations between characteristic length scales extracted from SAS profiles and microstructural descriptors, such as precipitate size and inter-precipitate distance. These findings provide a comprehensive understanding of the link between SAS profiles and microstructure evolution in Ni-based superalloys, offering valuable insights for materials characterization and design.

Microstructure