Search NASA⌕ Search

SEARCH · Search NASA

Results for “polycrystal”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Optical-Microphysical Cirrus Model

A model is presented that permits the simulation of the optical properties of cirrus clouds as measured with depolarization Raman lidars. It comprises a one-dimensional cirrus model with explicit microphysics and an optical module that transforms the microphysical model output to cloud and particle optical properties. The optical model takes into account scattering by randomly oriented or horizontally aligned planar and columnar monocrystals and polycrystals. Key cloud properties such as the fraction of plate-like particles and the number of basic crystals per polycrystal are parameterized in terms of the ambient temperature, the nucleation temperature, or the mass of the particles. The optical-microphysical model is used to simulate the lidar measurement of a synoptically forced cirrostratus in a first case study. It turns out that a cirrus cloud consisting of only monocrystals in random orientation is too simple a model scenario to explain the observations. However, good agreement between simulation and observation is reached when the formation of polycrystals or the horizontal alignment of monocrystals is permitted. Moreover, the model results show that plate fraction and morphological complexity are best parameterized in terms of particle mass, or ambient temperature which indicates that the ambient conditions affect cirrus optical properties more than those during particle formation. Furthermore, the modeled profiles of particle shape and size are in excellent agreement with in situ and laboratory studies, i.e., (partly oriented) polycrystalline particles with mainly planar basic crystals in the cloud bottom layer, and monocrystals above, with the fraction of columns increasing and the shape and size of the particles changing from large thin plates and long columns to small, more isometric crystals from cloud center to top. The findings of this case study corroborate the microphysical interpretation of cirrus measurements with lidar as suggested previously.

Reichardt, J.↗
DOE OSTI2026

Comparative analysis of plasticity-based GND density estimation methods in crystal plasticity finite element models

In crystal plasticity finite element (CPFE) simulations, accurately quantifying geometrically necessary dislocations (GNDs) is critical for capturing strain gradients in polycrystals. We compare different methods for quantifying GNDs, all of which originate from the Nye tensor, which is computed as the curl of the plastic deformation gradient. The projection technique directly decomposes the Nye tensor onto individual screw and edge dislocation components to compute GNDs. This approach requires converting a nine-component Nye tensor into densities for a larger number of dislocation systems, a fundamentally underdetermined (non-unique) process, which is resolved using L2 minimization. In contrast, when employing CPFE analysis, one could directly compute dislocation densities on each slip system using shear gradients. Projection and slip gradient methods are compared with respect to their prediction of GNDs with changing grain size, strain, and grain neighborhoods, including multigrain junctions. Although these techniques match analytical GND densities for single slip, single crystal deformation, and are consistent with anticipated overall GND trends, we find that the GND densities from projection techniques are significantly lower than those predicted from CPFE-based slip gradients in polycrystals. A suggested improvement of only using the active dislocation systems in the projection technique almost entirely resolved this mismatch.

Crystal plasticity↗
DOE OSTI2025

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Micromechanics of slip bands on a free surface

A micromechanics analysis for the formation and propagation of slip bands on the free surface of a polycrystal under monotonic loading is presented. For the growth of slip bands, the analysis satisfies the conditions of both equilibrium and displacement continuity, as well as the relation between slip and the resolved shear stress throughout the polycrystal. Numerical calculations show how the microstress field causes the concentration of plastic deformation in discrete sliding bands and give results which are in good qualitative agreement with known slip band observations on aluminum single crystals.

Lin, S. R.↗
NASA NTRS1984

Bubble formation in oxide scales on SiC

The oxidation of alpha-SiC single crystals and sintered alphaand beta-SiC polycrystals has been investigated at elevated temperatures. Bubble formation is commonly observed in oxide scales on polycrystalline SiC, but is rarely found on single-crystal scales; bubbles result from the preferential oxidation of C inclusions, which are abundant in SiC polycrystals. The absence of bubbles on single crystals, in fact, implies that diffusion of the gaseous species formed on oxidation, CO (or possibly SiO), controls the rate of oxidation of SiC.

Mieskowski, D. M.↗
NASA NTRS1992

Superplastic forming of ceramic insulation

Superplasticity has been demonstrated in many fine-grained structural ceramics and ceramic composites, including yttria-stabilized tetragonal zirconia polycrystal (YTZP), alumina, and Al2O3-reinforced zirconia (Al2O3/YTZ) duplex composites and SiC-reinforced Si3N4. These superplastic ceramics obviously offer the potential benefit of forming net shape or near net shape parts. This could be particularly useful for forming complicated shapes that are difficult to achieve using conventional forming techniques, or require elaborate, subsequent machining. In the present study, we successfully demonstrated the following: (1) superplastic 3Y-TXP and 20 percent Al2O3/YTZ composite have for the first time been successfully deformed into hemispherical caps via a biaxial gas-pressure forming technique; (2) no experimental difficulty was encountered in applying the required gas pressures and temperatures to achieve the results, thus, it is certain that higher rates of deformation than those presented in this study will be possible by using the current test apparatus at higher temperatures and pressures; and (3) an analytical model incorporating material parameters, such as variations during forming in the strain rate sensitivity exponent and grain growth-induced strain hardening, is needed to model accurately and therefore precisely control the biaxial gas-pressure forming of superplastic ceramics. Based on the results of this study, we propose to fabricate zirconia insulation tubes by superplastic extrusion of zirconia polycrystal. This would not only reduce the cost, but also improve the reliability of the tube products.

Nieh, T. G.↗
DOE OSTI2026

Mapping the Residual Stress Distribution in Polycrystalline Quartz Tiger's Eye Using Raman Spectroscopy

Stress distribution maps in polycrystalline materials are needed to reveal stress pathways caused by short‐and long‐range interactions of crystallites. Stress induced changes to wavenumbers in Raman spectroscopy is a well‐known phenomenon that occurs when a crystalline material undergoes elastic strain. Using spatially resolved Raman spectroscopy, we have developed a technique to experimentally measure grain‐scale residual stresses across a sample of polycrystal quartz Tiger's Eye. The results of this technique are validated using multiple criteria, including evaluation of the residual stress map itself and the accuracy of the measurements taken. The residual map of Tiger's Eye shows striking similarity to observed characteristics of other polycrystalline residual stress distributions, revealing heterogenous areas of high and low magnitude stresses. The stress magnitudes are consistent with residual stress magnitudes previously measured from polycrystalline quartz. Characterization of Tiger's Eye quartz also revealed the geometric nature of the iron oxide inclusions, not previously observed in Tiger's Eye quartz. With this technique, we present a residual stress map of Tiger's Eye quartz that shows with accuracy a heterogeneous stress state with inter‐ and intra‐granular detail.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complete development of LAROMance-type surrogate models for Zircaloy, as well as initial integration and testing in MOOSE, which addresses industry needs for licensing advanced LWR fuels

This milestone report details the progress in the constitutive surrogate modeling framework for Zircaloy’s mechanical behavior under extreme conditions, within the NEAMS program throughout the U.S. fiscal year 2024 (FY24). The presented work expands on the spectral based mechanical solver and the mechanistic material model for Zircaloy, which were introduced in the NEAMS program during FY22 and detailed in the previous fiscal year (FY23). As part of the NEAMS program, the current effort described in this report focused on the development of a novel viscplastic (creep) surrogate for the mechanistic, polycrystal-based model. To this end, the work was subdivided into two interdependent components related to the necessary improvements to the LApx constitutive model and algorithm and, to the development and use of a surrogate model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

The Influence of Mesoscale Particle Structure on the Electrode Degradation and Resultant Electrochemistry of Lithium Ion Cells With Nickel Rich Lithium Nickel Manganese Cobalt Oxide Positive Electrodes

Nickel-rich lithium nickel manganese cobalt oxide, NMC (LiNi x Mn y Co z O 2 ), materials are desirable positive electrodes in lithium ion batteries, providing high capacity and energy density. The mesoscale structure of NMC materials is commonly a polycrystalline aggregate, providing opportunity for short lithium ion transport distance of the small primary particles, yet facile material handling due to the larger secondary particles. On (de)lithiation the NMC unit cell changes volume, where the anisotropic strain can result in secondary particle fracture. This secondary particle fracture process during cycling has been associated with several degradation modes of NMC materials in LIBs. In this work, a milling process was determined whereby the secondary particles could be pre-fractured with retention of the parent primary particle crystallographic structure, crystallite size, and morphology, providing the ability to unambiguously determine the influence of NMC secondary particle size and mesostructure on resultant functional behavior. The bulk and interfacial properties and electrochemistry of as-received commercially obtained polycrystal NMC811 (LiNi 0.8 Mn 0.1 Co 0.1 O 2 , PC NMC) and its milled MPC NMC counterpart were compared. The smaller secondary particle size, higher surface area MPC NMC was found to result in 1) greater cathode tortuosity, 2) reduced surface Ni on the cycled MPC cathodes consistent with surface reconstruction, and 3) increased Ni deposition on the anodes from cathode-anode crosstalk. These factors manifested in unfavorable electrochemical behavior of MPC with lower functional capacity at both 1C and C/10 rates and higher impedance over extended moderate voltage (dis)charge cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructure Clones

Background: A material’s microstructure drives its material performance. Contemporary crystal plasticity experiments compare full-field strain measurements of polycrystal specimens to models. Because each specimen is unique, it is impossible to know which features of the observed deformation are deterministic vs statistical; thus, differences between model and experiment may or may not be significant. Objective: This paper introduces the invention of microstructure clones. Microstructure clones are 2D oligocrystal specimens that have nearly identical microstructures to remedy the aforementioned experimental limitations. Having specimens with nearly identical microstructures will allow for multiple destructive tests of a microstructure (either as repeats or intentionally different experiments), an ability to “see the future” by providing insight into how a specimen will deform, variability quantification, and experimental investigations of response to small microstructural changes. Methods: This work introduces microstructure clones. Repeatability of these clones is demonstrated in tensile bars of pure nickel. Local strain measurements from digital image correlation are compared between clone specimens and compared to results from a crystal plasticity finite element model. Results: Two sets of microstructure clones were tested in this study and displayed very consistent deformation responses within each clone set. Small observed differences in deformation invite investigation into microstructure stochasticity and the effect of small microstructural and loading differences. Conclusions: Microstructure clones represent a significant shift in understanding structure–property relationships. This work reshapes experimental crystal plasticity to allow for experiments that control for specific variables, quantification of microstructural stochasticity (and other sources of stochasticity), and opportunities for replicating experiments.

Crystal Plasticity↗
DOE OSTI2024

Data-Driven Surrogate Modeling with Microstructure-Sensitivity of Viscoplastic Creep in Grade 91 Steel

Abstract To support the development of advanced steel alloys tailored to withstand extreme conditions, it is imperative to account for the mechanical performance of components, while considering the influence of local microstructure on the macroscopic response. To this end, this study focuses on the development of microstructure-sensitive constitutive models for the mechanical response of Grade 91 steel exposed to extreme thermo-mechanical environments. Polynomial chaos expansion (PCE) surrogates are used to emulate high-fidelity polycrystal simulations of the viscoplastic response of Grade 91 steel as a function of the microstructure fingerprint (e.g., dislocations and precipitates). To cover a wide temperature–stress domain, two separate PCE surrogates—one that captures softening and the other that captures hardening behavior—are combined using another (sparse) Gaussian process regression model. The resulting constitutive creep surrogate model is integrated within the MOOSE finite element framework to simulate the intricate effects of microstructure, in particular MX-phase precipitates, on a component with a graded microstructure. Surrogate sensitivity analysis is applied to quantify the relevant impact of spatially varying microstructure on the creep response in a test-case involving a Grade 91 alloy with a prototypical weld.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Grain boundary solute segregation across the 5D space of crystallographic character

Solute segregation in materials with grain boundaries (GBs) has emerged as a popular method to thermodynamically stabilize nanocrystalline structures. However, the impact of varied GB crystallographic character on solute segregation has never been thoroughly examined. Here, this work examines Co solute segregation in a dataset of 7272 Al bicrystal GBs that span the 5D space of GB crystallographic character. Considerable attention is paid to verification of the calculations in the diverse and large set of GBs. In addition, the results of this work are favorably validated against similar bicrystal and polycrystal simulations. As with other work, we show that Co atoms exhibit strong segregation to sites in Al GBs and that segregation correlates strongly with GB energy and GB excess volume. Segregation varies smoothly in the 5D crystallographic space but has a complex landscape without an obvious functional form.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond curvature-driven grain growth: Insights from fully anisotropic Monte Carlo Potts simulations

Grain boundary (GB) motion away from the center of curvature, termed anti-curvature behavior, has recently been observed in 3D experiments but is not predicted by classical grain growth theory. In this study, we investigate this behavior using a novel, fully anisotropic Monte Carlo Potts (MCP) model that incorporates both misorientation and inclination dependencies of GB energy. We perform 3D grain growth simulations with isotropic and anisotropic GB energies to explore the relationship between GB velocity and curvature. Contrary to the classical relation that velocity is a product of reduced mobility and mean curvature, we observe no consistent correlation between velocity and mean curvature for individual GBs, even under isotropic conditions, though there is correlation between the average velocity and curvature for some cases. The 3D simulations exhibit frequent anti-curvature motion regardless of GB energy anisotropy including with isotropic GBs, though larger curvatures occur with anisotropic GB energy functions that promote low-energy GBs. In 2D simulations, anti-curvature behavior only occurs with the anisotropic functions that favor low-energy GBs. This difference between 2D and 3D results suggests that anti-curvature behavior results in part from the increased freedom of motion intrinsic in 3D GB networks. Furthermore, our results support recent experimental observations that demonstrate that simple curvature-driven models are insufficient for describing GB migration in polycrystals.

Anti-curvature↗
DOE OSTI2025

A comparison of surrogate constitutive models for viscoplastic creep simulation of HT-9 steel

Mechanistic microstructure-informed constitutive models for the mechanical response of polycrystals are a cornerstone of computational materials science. However, as these models become increasingly more complex – often involving coupled differential equations describing the effect of specific deformation modes – their associated computational costs can become prohibitive, particularly in optimization or uncertainty quantification tasks that require numerous model evaluations. To address this challenge, surrogate constitutive models that balance accuracy and computational efficiency are highly desirable. Data-driven surrogate models, that learn the constitutive relation directly from data, have emerged as a promising solution. In this work, we develop two local surrogate models for the viscoplastic response of a steel: a piecewise response surface method and a mixture of experts model. These surrogates are designed to adapt to complex material behavior, which may vary with material parameters or operating conditions. The surrogate constitutive models are applied to creep simulations of HT-9 steel, an alloy of considerable interest to the nuclear energy sector due to its high tolerance to radiation damage, using training data generated from viscoplastic self-consistent (VPSC) simulations. In conclusion, we define a set of test metrics to numerically assess the accuracy of our surrogate models for predicting viscoplastic material behavior, and show that the mixture of experts model outperforms the piecewise response surface method in terms of accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling inter- and intra-granular dislocation transport using crystal plasticity

Here, this work presents the development of a crystal plasticity material model that incorporates both dislocation transport within grains and dislocation transfer across grain boundaries. This model has been implemented in the open-source finite element code MOOSE. In addition, a novel geometry-based criterion is developed to determine the direction of dislocation transfer across grain boundaries. The transfer criterion incorporates the geometric features of the grain boundary, such as the grain boundary plane normal, and its misorientation, which is accounted for through the orientation of the incoming and outgoing slip systems. The model is tested with several cases, including a copper single crystal, bi-crystal, and polycrystal. The development of the transfer criterion, implementation of the model, and its application to these test cases are discussed in detail.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Statistical evaluation of microscale stress conditions leading to void nucleation in the weak shock regime

Here, we investigate the heterogeneity of the stress state driven by anisotropic deformation response at the single crystal level through five statistical volume element (SVE) calculations of polycrystalline BCC tantalum. This work focuses on grain boundaries as a prominent material defect type prone to void nucleation based upon experimental observations of predominantly intergranular void nucleation in this material. The SVEs are constructed to be statistically representative of larger volumes of material and are meshed such that mean and standard deviation of grain size and orientation information is reconstructed. The computational meshes feature hexahedral (brick) elements and smooth conformal grain boundaries where significant stress concentration is known to occur, a tail effect of interest in the extreme events process of dynamic ductile damage. An existing micromechanical crystallographic plasticity model shown to capture the single crystal behavior of BCC tantalum well is used to perform the polycrystal calculations. The model includes representation of the non-Schmid effect of non-planar screw dislocation kinetics in tantalum. A three-dimensional stress state time profile predicted by damage modeling of a flyer plate impact experiment is applied as boundary conditions to each SVE. Resulting grain boundary stress state statistics are strongly non-Gaussian. Significant structural evolution is observed within the compressive hold before unloading into tension in the stress profile. Strong angular dependence of grain boundary traction magnitude with shock direction is observed. Non-Schmid effects continue to suggest their influence on propensity of microstructural defect types to nucleate voids. A general void nucleation criterion is proposed using probability theory. The general framework is specified to polycrystalline BCC tantalum in the weak shock regime to include the SVE calculations and literature molecular dynamics calculations of grain boundary void nucleation strength. Probability density functions (PDFs) are used to describe the interaction between the local stress state heterogeneity and the distributed grain boundary void nucleation strength state. A causation entropy maximization procedure removes the requirement for ad hoc selection of a PDF functional form and provides a rigorous procedure for data-based PDF determination. The resulting physically informed PDF describes the spatial appearance frequency of nucleated voids as a function of applied macroscale pressure. Lower length scale physics are thus packaged in a precise and computationally efficient way to provide computational plasticity insight to macroscale dynamic ductile damage models.

36 MATERIALS SCIENCE↗