Search NASA⌕ Search

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.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Learning grain boundary segregation energy spectra in polycrystals

The segregation of solute atoms at grain boundaries (GBs) can profoundly impact the structural properties of metallic alloys, and induce effects that range from strengthening to embrittlement. And, though known to be anisotropic, there is a limited understanding of the variation of solute segregation tendencies across the full, multidimensional GB space, which is critically important in polycrystals where much of that space is represented. Here we develop a machine learning framework that can accurately predict the segregation tendency—quantified by the segregation enthalpy spectrum—of solute atoms at GB sites in polycrystals, based solely on the undecorated (pre-segregation) local atomic environment of such sites. We proceed to use the learning framework to scan across the alloy space, and build an extensive database of segregation energy spectra for more than 250 metal-based binary alloys. The resulting machine learning models and segregation database are key to unlocking the full potential of GB segregation as an alloy design tool, and enable the design of microstructures that maximize the useful impacts of segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Diffusion and Deformation Mechanism Maps in Li Metal from Atomistic Simulations

The rate of Li transport in the Li metal anode is important for the operation of Li metal-solid state batteries. Here, transport rates due to diffusion and creep are predicted in Li using atomistic simulations. First, molecular dynamics is used to estimate the rate of Li diffusion along dislocations and in grain boundary triple junctions. By combining this data with that from a prior study of grain boundary diffusion the dominant mechanisms and rates of self-diffusion in Li polycrystals are predicted as a function of grain size, grain shape, dislocation density, and temperature. Second, the dominant creep mechanisms are predicted and used to estimate critical current densities and void annihilation times. Grain boundary sliding and coble creep are the dominant mechanisms for micron-sized grains. Lastly, a continuum model for interfacial contact loss reveals that high dislocation densities of ∼10 12 /cm 2 enable achieving battery performance targets for Li grain sizes of ∼10 μm.

Batteries↗
DOE OSTI2022

Development of a concurrent coupled Atomistic - Continuum model to predict the defect and grain structure for Additive Manufacturing process [Slides]

Development of a coupled atomistic-continuum model for metal AM process. Two-way coupling between atomistic (heat flux) and continuum (temperature) domain. Identification of grain and dislocation structure at the atomistic domain. Conversion of dislocation data from discrete atomistic to density form at the continuum. Obtain the spatial distribution of grain structure and dislocation density in additively manufactured metal polycrystal material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Transformation-induced fracture toughening in CuAlBe shape memory alloys: A phase-field study

Stress-induced transformation ahead of a crack tip can relax the stress concentration and produce fracture toughness in shape memory alloys (SMAs). In this manuscript, we utilize a nonisothermal phase-field model (PFM) to study the martensitic transformation induced crack tip toughening in CuAlBe SMA. The force-displacement curve, transformation zone and high stress zone in single crystalline samples show high dependency on the grain orientation with respect to the crack alignment. Comparison between isothermal and non-isothermal simulations reveals that the transformation-induced self-heating decreases the toughening capability by increasing the critical transformation stress. Investigating the high stress zones in CuAlBe and the nontransforming CuAl alloy shows that the toughening is obtained by redistributing the stress concentration far from the crack tip, as the high stress zone follows the tip of the transformation zone. Increasing the acuity of the crack tip is found to generate more symmetric martensite wings on both side of the crack axis and a more localized high stress zone. The polycrystal specimen displays higher toughening due to the internal constraints related to the presence of various grains with difference orientations. Depending on the orientation of the grain inclosing the crack tip, the toughening can be lower or higher. Finally, coincidence of the crack tip with a triple junction is found to improve the toughening behavior.

42 ENGINEERING↗
DOE OSTI2021

Evaluating the grain-scale deformation behavior of a single-phase FCC high entropy alloy using synchrotron high energy diffraction microscopy

We report although the deformation behavior of high-entropy alloys (HEAs) has been extensively studied at the macroscale, many important properties have yet to be explored for these alloys at the microscale, thus hampering accurate prediction of damage and failure processes. Synchrotron high-energy diffraction microscopy (HEDM) and fast-Fourier transform-based crystal plasticity modeling was conducted to investigate the three-dimensional (3D) grain-resolved micromechanical response for approximately 1,900 constituent grains within a single-phase FCC HEA up to 1% applied strain. The evolution of grain-resolved elastic strains, lattice reorientations, and maximum resolved shear stresses (mRSS) were evaluated to quantify elastic, yield, and fully plastic behavior. Overall, the initial critical resolved shear stress (CRSS), determined via in situ HEDM and companion modeling, was found to be > 20% higher than estimated using the classical polycrystalline Taylor factor (M = 3.06). However, a descriptive parameter based on the average grain-resolved Taylor factor ($\overline{M}$) was found to show excellent agreement with plastic yielding behavior observed within HEDM datasets. Noticeable deviations in HEDM lattice reorientations compared to both EVP-FFT simulations and classical predictions for FCC polycrystals were discovered, highlighting the complexity in correlating local lattice reorientations, Taylor, and Schmid factors with plastic response for this material at the grain-scale. Therefore, it is anticipated that the overall trends and parameter identification of 3D grain-resolved properties in this study can serve as an important foundation for continued mesoscale investigation on both well-established and newly developed Cantor-like HEAs.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

A Convolutional Neural Network for Enhancement of Multi-Scale Localization in Granular Metallic Representative Unit Cells

A convolutional neural network was used to enhance the localization of strain and stress for a generalized method of cells model of a metallic microstructure. Enhanced shear strains, measured in terms of the linear regression coefficients as a function of ground truth strains, were improved from inaccurate and uncorrelated (slope=0.003, Rsq=0.000) to accurate and well correlated (slope=0.890, Rsq=0.882) relative to ground truth (slope=1.0, Rsq=1.0). In applying the convolutional neural network, a convolutional stride of 1.0 (padding=’same’) was only modestly effective while strides of 2 or 3 were more effective yet at higher cost. Additional convolutional layers were generally more expensive than additional dense layers, often with limited benefit. The accuracy of enhanced localized shear strains and stress is expected to yield benefits for damage progression models, especially in the context of hierarchical multi-scale methods where the generalized method of cells is applied at the intermediate scale.

Machine Learning↗
DOE OSTI2024

Enabling resonant ultrasound spectroscopy in high magnetic fields

Resonant ultrasound spectroscopy (RUS) is a powerful method to determine elastic constants with high accuracy and precision from a single measurement of the mechanical resonances of a sample. Conventionally, the quantitative extraction of elastic moduli with RUS assumes free boundary conditions which can often lead to the adoption of unstable sample positioning between ultrasonic transducers that is incompatible with extreme environments like high magnetic fields. We show that, under specific conditions, introducing a small amount of adhesive between a RUS sample and ultrasonic transducers introduces a perturbation to the free resonance condition which can be accounted for by a simple model. This means elastic constants can be determined to within the uncertainty of conventional RUS, but with significant improvements including sample stability and control of sample orientation. We demonstrate the efficacy of this approach with measurements on a range of materials including room temperature measurements on polycrystalline metals, temperature-dependent measurements of the structural phase transition in strontium titanate single crystals, and magnetic field-dependent measurements of magnetic phase transitions in gadolinium polycrystals up to 14 T.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain size effects on slip band development

Crystallographic slip localizations, known as slip bands, concentrate stress in polycrystals, often leading to the nucleation of damage. Slip band development has been experimentally shown to be sensitive to grain size, tending to develop more frequently and with a greater intensity in large grains. In this work, we investigate the influence of grain size on the propensity for crystallographic slip band development. To this end, we employ the slip band-fast Fourier transform method (SB-FFT). SB-FFT is a 3D, full-field crystal plasticity model that permits the incremental development of discrete crystallographic slip bands according to microstructure and material properties. We present a model Inconel 718 tricrystal to isolate the effect of grain size. Our findings show that slip bands in large grains develop at lower applied strain levels and at a faster rate than slip bands in small grains. The grain size effect is due to a backstress produced by the interaction of the slip band and its neighboring grain. The backstress is most intense at small grain sizes, impeding slip activity within a developing slip band and immediately surrounding the slip band.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electronic structure and two-band superconductivity in unconventional high- T c cuprates Ba 2 CuO 3+δ

We report that the recently discovered cuprate superconductor Ba 2 CuO 3+δ exhibits a high T c ≃73 K at δ ≃ 0.2. The polycrystal grown under high pressure has a structure similar to La 2 CuO 4 but with dramatically different lattice parameters due to the CuO 6 octahedron compression. The crystal field in the compressed Ba 2 CuO 4 leads to an inverted Cu 3d e g complex with the d x 2 -y 2 orbital sitting below the d 3z 2 - r 2 and an electronic structure highly unusual compared to the conventional cuprates. We construct a two-orbital Hubbard model for the Cu d 9 state at hole doping x = 2δ and study the orbital-dependent strong correlation and superconductivity. For the undoped case at x = 0 , we found that strong correlation drives an orbital-polarized Mott-insulating state with the spin-1/2 moment of the localized d 3z 2 - r 2 orbital. In contrast to the single-band cuprates where superconductivity is suppressed in the overdoped regime, hole doping the two-orbital Mott insulator leads to orbital-dependent correlations and the robust spin and orbital exchange interactions produce a high-T c antiphase d-wave superconductor even in the heavily doped regime at x = 0.4 . We conjecture that Ba 2 CuO 3+δ realizes mixtures of such heavily hole-doped superconducting Ba 2 CuO 4 and disordered Ba 2 CuO 3 chains in a single-layer or predominately separated bilayer structure. Our findings suggest that unconventional cuprates with liberated orbitals as doped two-band Mott insulators can be a direction for realizing high-T c superconductivity with enhanced transition temperature T c .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Computed entropy spectra for grain boundary segregation in polycrystals

Abstract Grain boundary solute segregation influences most bulk material properties, and understanding solute thermodynamics at grain boundaries is critical for engineering them. However, the vast grain boundary space in polycrystals is challenging to evaluate due to its size, especially for the intrinsically hard-to-compute segregation excess entropy. Here data science methods are used to generate a database of site-wise grain boundary segregation entropy spectra for 155 dilute binary alloys within the harmonic approximation. The spectral framework allows scale bridging between the calculated atomistic site-wise energy-entropy spectra and macroscopic segregation entropy estimates. The results affirm that macroscopic averaging is not sufficient: a spectral treatment of grain boundary segregation is needed to accurately model bulk temperature dependence of grain boundary solute segregation. The calculated spectral entropy database and thermodynamic framework can be applied for both understanding segregation experiments and alloy design exercises, paving the way to a finite-temperature grain boundary genome.

Chemistry↗
NASA NTRS1995

Investigation of Strain Aging in the Ordered Intermetallic Compound beta-NiAl

The phenomenon of strain aging has been investigated in polycrystalline and single crystal NiAl alloys at temperatures between 300 and 1200 K. Static strain aging studies revealed that after annealing at 1100 K for 7200 s (i.e., 2h) followed by furnace cooling, high purity, nitrogen-doped and titanium-doped polycrystalline alloys exhibited continuous yielding, while conventional-purity and carbon-doped alloys exhibited distinct yield points and Luders strains. Prestraining by hydrostatic pressurization removed the yield points, but they could be reintroduced by further annealing treatments. Yield points could be reintroduced more rapidly if the specimens were prestrained uniaxially rather than hydrostatically, owing to the arrangement of dislocations into cell structures during uniaxial deformation. The time dependence of the strain aging events followed at t(exp 2/3) relationship suggesting that the yield points observed in polycrystalline NiAl were the result of the pinning of mobile dislocations by interstitials, specifically carbon. Between 700 and 800 K, yield stress plateaus, yield stress transients upon a ten-fold increase in strain rate, work hardening peaks, and dips in the strain rate sensitivity (SRS) have been observed in conventional-purity and carbon-doped polycrystals. In single crystals, similar behavior was observed; in conventional-purity single crystals, however, the strain rate sensitivity became negative resulting in serrated yielding, whereas, the strain rate sensitivity stayed positive in high purity and in molybdenum-doped NiAl. These observations are indicative of dynamic strain aging (DSA) and are discussed in terms of conventional strain aging theories. The impact of these phenomena on the composition-structure-property relations are discerned. Finally, a good correlation has been demonstrated between the properties of NiAl alloys and a recently developed model for strain aging in metals and alloys developed by Reed-Hill et al.

Weaver, Mark Lovell↗
DOE OSTI2021

Understanding magnetic phase coexistence in Ru 2 Mn 1-x Fe x Sn Heusler alloys: A neutron scattering, thermodynamic, and phenomenological analysis

The random substitutional solid solution between the antiferromagnetic (AFM) full-Heusler alloy Ru 2 MnSn and the ferromagnetic (FM) full-Heusler alloy Ru 2 FeSn provides a rare opportunity to study FM-AFM phase competition in a near-lattice-matched, cubic system, with full solubility. At intermediate x in Ru 2 Mn 1-x Fe x Sn this system displays suppressed magnetic ordering temperatures, spatially coexisting FM and AFM order, and strong coercivity enhancement, despite rigorous chemical homogeneity. Here, we construct the most detailed temperature- and x-dependent understanding of the magnetic phase competition and coexistence in this system to date, combining wide-temperature-range neutron diffraction and small-angle neutron scattering with magnetometry and specific heat measurements on thoroughly characterized polycrystals. A complete magnetic phase diagram is generated, showing FM-AFM coexistence between x ≈ 0.30 and x ≈ 0.70. Furthermore, important insight is gained from the extracted length scales for magnetic phase coexistence (25–100 nm), the relative magnetic volume fractions and ordering temperatures, and remarkable x-dependent trends in magnetic and electronic contributions to specific heat. An unusual feature in the magnetic phase diagram (an intermediate FM phase) is also shown to arise from an extrinsic effect related to a minor Ru-rich secondary phase. The established magnetic phase diagram is then discussed with the aid of phenomenological modeling, clarifying the nature of the mesoscale phase coexistence with respect to the understanding of disordered Heisenberg models.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Subcritical Fracture Propagation in Rocks: an Examination Using the Methods of Fracture Mechanics and Non-Destructive Testing

An experimental investigation of tensile rock fracture is presented with an emphasis on characterizing time dependent crack growth using the methods of fracture mechanics. Subcritical fracture experiments were performed in moist air on glass and five different rock types at crack velocities using the double torsion technique. The experimental results suggest that subcritical fracture resistance in polycrystals is dominated by microstructural effects. Evidence for gross violations of the assumptions of linear elastic fracture mechanics and double torsion theory was found in the tests on rocks. In an effort to obtain a better understanding of the physical breakdown processes associated with rock fracture, a series of nondestructive evaluation tests were performed during subcritical fracture experiments on glass and granite. Comparison of the observed process zone shape with that expected on the basis of a critical normal principal tensile stress criterion shows that the zone is much more elongated in the crack propagation direction than predicted by the continuum based microcracking model alone.

Swanson, P. L.↗
DOE OSTI2023

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Triple junction solute segregation in Al-based polycrystals

Solute segregation is a crucial means of stabilizing nanostructured alloys, and at very small grain sizes, this requires consideration of triple junctions (TJs), which attain a meaningful volume fraction and thus become relevant for bulk material behavior. Here, in this study, the solute segregation spectra for grain boundary (GB) and TJ sites are calculated for a large number of dilute Al-based binary systems with available interatomic potentials. A defect-identifying algorithm is applied to quantify the average GB thickness and classify the intergranular site spectra into GB and junction subspectra. The algorithm is also applied to a hybrid electronic-level database for GB segregation of various solutes in Al, yielding polycrystalline TJ solute segregation spectra from first principles. The results suggest that TJ segregation is alloy or interatomic potential dependent and can exhibit either boundary or junction preference. With these spectra as inputs, the spectral GB segregation model gives quantitative predictions of segregation as a function of grain size, temperature, and total solute concentration, suitable for alloy screening and design.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗