Anisotropic power diagrams for polycrystal modelling: Efficient generation of curved grains via optimal transport
Not Available
Buze, M. (ORCID:000000018617974X)↗
SEARCH · Search NASA
Results for “Polycrystal Modeling”
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.
Data-Driven Constitutive Model for the Inelastic Response of Metals: Application to 316H Steel
Here, predictions of the mechanical response of structural elements are conditioned by the accuracy of constitutive models used at the engineering length-scale. In this regard, a prospect of mechanistic crystal-plasticity-based constitutive models is that they could be used for extrapolation beyond regimes in which they are calibrated. However, their use for assessing the performance of a component is computationally onerous. To address this limitation, a new approach is proposed whereby a surrogate constitutive model (SM) of the inelastic response of 316H steel is derived from a mechanistic crystal plasticity-based polycrystal model tracking the evolution of dislocation densities on all slip systems. The latter is used to generate a database of the expected plastic response and dislocation content evolution associated with several instances of creep loading. From the database, a SM is developed. It relies on the use of orthogonal polynomial regression to describe the evolution of the dislocation content. The SM is then validated against predictions of the dead load creep response given by the polycrystal model across a range of temperatures and stresses. When the SM is used to predict the response of 316H during complex non monotonic loading, extrapolating to new loading conditions, it is found that predictions compare particularly well against those from the physics-based polycrystal model.
36 MATERIALS SCIENCE↗
Grain Size-Dependent Crystal Plasticity Constitutive Model for Polycrystal Materials
A new method to introduce grain size-dependence in classical crystal plasticity constitutive model is developed by considering the resistance to dislocation motion in the grain boundary influence region as equivalent to that of a work hardening. A general framework for the size-dependent constitutive model is derived by implementing this method on a core and mantle model. The work hardening, equivalent to grain boundary effect, is realized by introducing a fictitious, pre-existing, plastic strain grain boundary influence region (mantle) following the principles of classical crystal plasticity. This fictitious plastic strain, in effect, increases the yield strength and decreases the initial hardening coefficient of the grain. With the thickness of grain boundary influence region and the distribution of introduced plastic strain remaining the same, the grain boundary effect increased as the grain size becomes smaller. A simplified model that considers the grain boundary effect on a grain average sense is also developed under this general framework. Implementation of this general framework to the specific cases of crystal plasticity constitutive models is demonstrated by considering the case of power law flow rule and hyperbolic-secant hardening rule. Finally, the grain size-dependent constitutive model is validated by comparing the predicted stress-strain behavior of polycrystal samples with different average grain sizes under uniaxial loading with the experimental results.
Grain Size↗
Multiscale Modeling of Structurally-Graded Materials Using Discrete Dislocation Plasticity Models and Continuum Crystal Plasticity Models
A multiscale modeling methodology that combines the predictive capability of discrete dislocation plasticity and the computational efficiency of continuum crystal plasticity is developed. Single crystal configurations of different grain sizes modeled with periodic boundary conditions are analyzed using discrete dislocation plasticity (DD) to obtain grain size-dependent stress-strain predictions. These relationships are mapped into crystal plasticity parameters to develop a multiscale DD/CP model for continuum level simulations. A polycrystal model of a structurally-graded microstructure is developed, analyzed and used as a benchmark for comparison between the multiscale DD/CP model and the DD predictions. The multiscale DD/CP model follows the DD predictions closely up to an initial peak stress and then follows a strain hardening path that is parallel but somewhat offset from the DD predictions. The difference is believed to be from a combination of the strain rate in the DD simulation and the inability of the DD/CP model to represent non-monotonic material response.
Saether, Erik↗
Sensitivity Study of Multiscale and Phenomenological Elasto-Viscoplastic Grade 91 Material Models for Component-Scale Response
Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.
42 ENGINEERING↗
Size effect in stainless steel thin wires under tension
Miniaturization of modern devices drives the exploration into how materials behave in small dimensions. However, the mechanical behaviors of micro-scale specimens under varied loading conditions are still challenging to characterize and predict. In this work, specimen size-dependent weakening effect was experimentally revealed in 316L stainless steel thin wires under tension with wire diameter decreasing from 200 μm down to 61 μm. The underlying mechanisms were primarily associated with surface grain softening and strain localization, which were largely dependent on the number of grains across wire diameter. Through an analog composite model with surface and interior grains that follows the Hall-Petch relation, the extent of surface grain softening was found to coincide with the flow stress reduction due to lack of boundary strengthening. Further using elasto-plastic self-consistent polycrystal modeling, the stress-strain behaviors of interior and surface grains were well described in analogy to those of bulk materials with similar-size grains and extremely large grains, respectively. Subsequently, the analog composite model well simulated the specimen size-dependent weakening effect under tension. This article provided more insights into the size effect exploration and proposed a practical model for describing size-dependent mechanical behaviors of small-scale materials to be applied to miniaturized devices.
316L stainless steel↗
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↗
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↗
Achieving High Efficiency in Reduced Order Modeling for Large Scale Polycrystal Plasticity Simulations
Reduced order models for the nonlinear response of heterogeneous microstructures typically require a construction (or training) stage to build the reduced order basis. In this manuscript, an efficient model construction strategy for the eigenstrain homogenization method (EHM) is presented. The proposed strategy relies on a parallel, element-by-element, conjugate gradient solver. Near linear scaling has been achieved with respect to the number of degrees of freedom used to resolve the microstructure. Linear scaling with respect to the number of pre-analyses required to construct the reduced order model (ROM) follows from the EHM formulation. Furthermore, a parallel implementation for fast evaluation of the constructed ROM has been developed using shared memory parallelization. It has been shown that for large microstructures with ≈ 10,000 grains, the total computational cost of evaluating the nonlinear response of a polycrystal could be reduced by approximately an order of magnitude using 32 cores with respect to serial ROM simulation. The present methodology has been verified using an additively manufactured polycrystalline microstructure of a nickel-based superalloy, Inconel 625. The capability of the developed framework to construct a ROM for such large microstructures, as well as the ability of the ROM to predict average and local quantities of interest has been demonstrated.
microscale↗
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↗
Mechanistic Modeling of Creep and Aging for Time Extrapolation in LPBF 316 SS - LANL
A polycrystal constitutive model was extended to quantify the effects heterogeneity in the microstructure, porosity, and internal stresses on the mechanical response (creep, tensile) of 316H steel processed by additive manufacturing. The starting point is a recently developed and fully calibrated advanced microstructuresensitive mechanistic constitutive model fit to literature data pertaining to conventionally processed steels. The model is embedded within the full field fast Fourier transform framework (LApx) and several hundred simulations were performed. Synthetic microstructures with different grain aspect ratios were generated using the Neper package to study the role of grain morphology on the mechanical responses. The effect of initial cell interior dislocation density, cell-wall density, precipitate content, porosity, and internal stresses on creep rates are quantified for different polycrystalline microstructures. The study allowed to prospectively quantify the role of distinct features of the microstructure on the mechanical response under tension and creep loading scenarios. Critically it is found that minor anisotropies in the tensile response of samples, emanating from the aspect ratio of grains correlated into significant differences in creep rate as a function of loading directions. The study also allows to identify initial microstructure porosities below which the materials response only moderately degraded by the initial microstructure.
36 MATERIALS SCIENCE↗
Using machine-learning to understand complex microstructural effects on the mechanical behavior of Ti-6Al-4V alloys
Structural materials properties are highly dependent on their microstructure. Their microstructure is in turn affected by multiple fabrication and thermo-mechanical treatment parameters, all of which conform a highly-dimensional parametric space with often hidden correlations that are difficult to extract by experimentation alone. This is particularly true for alloys of the dual-phase Ti-6Al-4V family, with their greatly complex and rich microstructures, which combine several intrinsic length scales associated with multiple grain and subgrain structures, grains with different crystal lattices (α and β phases), and complex chemistry. In this paper we use a comprehensive set of machine learning techniques to develop predictive tools relating the yield strength and hardening rate of these alloys to a set of input parameters covering extensive ranges. The data generator is a finite-element crystal plasticity model for polycrystal deformation that takes into account slip anisotropy and employs standard dislocation evolution models for the α and β phases of Ti-based alloys. Our dataset includes over two thousand independent simulations and is used to train the machine learning models, which are then used to establish correlations between microstructural parameters and the alloys’ mechanical response. Our results point to the most influential parameters affecting yield strength and hardening rate, information that can then be used to guide experimental synthesis and characterization efforts to save time and resources.
36 MATERIALS SCIENCE↗
Analytic model of dislocation density evolution in fcc polycrystals accounting for dislocation generation, storage, and dynamic recovery mechanisms
Here, an analytic model of the evolution of dislocation density in fcc polycrystals is described. The evolution equations approximately account for most known dislocation storage, dynamic recovery, and dislocation generation mechanisms in fcc polycrystals. Specifically, the model incorporates network (forest) and grain boundary storage, mobile-network and mobile–mobile annihilation, screw–screw annihilation via athermal and thermal single cross-slip, generation by double cross-slip (Koehler mechanism, including dipole formation), Frank-Read sources, grain boundary nucleation, and mobile–immobile dislocation nucleation due to shock loading. Single cross-slip is assumed to proceed through the Friedel–Escaig (FE) mechanism; the corresponding activation energy is calculated using a modified FE model. The activation energy for double cross-slip is calculated for the first time by extending the FE model. The exact evolution equations are integro-differential equations, and as such are difficult to implement in a code; hence, the evolution equations are simplified by making several approximations. Preliminary results on copper are presented, including comparisons to experimental data.
36 MATERIALS SCIENCE↗
Constitutive Models for Shape Memory Alloy Polycrystals
Shape memory alloys (SMA) exhibiting the superelastic or one-way effects can produce large recoverable strains upon application of a stress. In single crystals this stress and resulting strain are very orientation dependent. We show experimental stress/strain curves for a Ni-Al single crystal for various loading orientations. Also shown are model predictions; the open and closed circles indicate recoverable strains obtained at various stages in the transformation process. Because of the strong orientation dependence of shape memory properties, crystallographic texture can be expected to play an important role in the mechanical behavior of polycrystalline SMA. It is desirable to formulate a constitutive model to better understand and exploit the unique properties of SMA.
Comstock, R. J., Jr.↗
An FFT-based approach for Bloch wave analysis: application to polycrystals
A method based on the Fast Fourier Transform is proposed to obtain the dispersion relation of acoustic waves in heterogeneous periodic media with arbitrary microstructures. The microstructure is explicitly considered using a voxelized Representative Volume Element (RVE). The dispersion diagram is obtained solving an eigenvalue problem for Bloch waves in Fourier space. To this aim, two linear operators representing stiffness and mass are defined through the use of differential operators in Fourier space. The smallest eigenvalues are obtained using the implicitly restarted Lanczos and the subspace iteration methods, and the required inverse of the stiffness operator is done using the conjugate gradient with a preconditioner. The method is used to study the propagation of acoustic waves in elastic polycrystals, showing the strong effect of crystal anistropy and polycrystaline texture on the propagation. It is shown that the method combines the simplicity of classical Fourier series analysis with the versatility of Finite Elements to account for complex geometries proving an efficient and general approach which allows the use of large RVEs in 3D.
97 MATHEMATICS AND COMPUTING↗
Deformation of Hydrous Phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and Stishovite [Si 1‐n H 4n O 2 ] Relevant to Anisotropy of the Earth's Mantle
Seismic anisotropy of the Earth's mantle has been mostly attributed to crystallographic preferred orientation (CPO) generated during subduction and convection of an anhydrous mantle. But some hydrous phases are also stable at mantle conditions. Here we present results from diamond-anvil cell deformation experiments at high pressure and temperature on hydrous phases Egg [AlSiO 3 (OH)], δ [AlO(OH)] and hydrous stishovite [Si 1-n H 4n O 2 ], transformed from the clay mineral kaolinite. They develop strong CPO during axial compression, suggesting that they likely contribute to seismic anisotropy and heterogeneity in the mantle. Comparing experimental results with viscoplastic polycrystal plasticity models suggest that phase Egg deforms dominantly by (001) slip, δ by (010) slip and stishovite by {100} slip which could be incorporated in future models of mantle geodynamics.
58 GEOSCIENCES↗