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.

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↗

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↗

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↗

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.↗

Geometric Requirements for Modeling Polycrystalline Representative Volume Elements Using the Generalized Method of Cells

To accurately predict the mechanical response of advanced metallic parts, the grain structure of the metal must be considered and coupled to the macroscale in a multiscale model. The size of the microscale representative volume element must be sufficiently large, and the grid used to discretize the microscale must be adequately refined, in order to ensure an accurate and converged solution. The generalized method of cells is used to predict the mechanical properties of polycrystalline representative volume elements of copper containing randomly oriented grains with the properties of a single face-centered-cubic crystal. Moreover, the results of a parametric study are used to conclude the geometric requirements for an objective microscale model, with respect to the elastic stiffness.

Multiscale Modeling↗

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.↗

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗

Polycrystal Snowfakes

Building on the previous success of a mono-crystal snowflake growth model, a team of NASA researchers has invented and implemented a poly-crystal variant of the model to simulate simultaneous crystalline growth along lattices of different orientations. Inferring snowfall intensity and amount using remote sensing relies on knowing the scattering properties of snowing particles. However, these scattering properties depend strongly on the geometric shapes of the snowing particles. Poly-crystal snow particles, such as bullet rosettes, appear frequently in nature. Simulating the growth of both mono-crystal and poly-crystal particles and obtaining accurate scattering properties are paramount for furthering accurate quantitative estimates of snowfalls.

Clune, Tom↗

Critical experiments of the self-consistent model for polycrystalline Hastelloy-X

A viscoplastic constitutive model is presented for the estimation of the overall mechanical response of Hastelloy-X polycrystalline metals from a knowledge of single crystal behavior. The behavior of polycrystal is derived from that of single crystals using a self-consistent formulation. The single crystal behavior which has been used was developed by summing postulated slip on crystallographic slip systems. The plasticity and creep are treated coupledly using unified viscoplastic model which includes the interaction effects between rapid and slow deformation at elevated temperature. The validity of the model is directly tested by experiments on Hastelloy-X in both single crystal and polycrystalline versions.

Shi, Shixiang↗

Exposure damage mechanisms for KCl windows in high power laser systems

An experimental study of the 10.6 micrometer and 0.6328 micrometer optical properties of single crystal and europium doped polycrystal is described. Significant variations in the optical properties are observed over periods of exposure up to 100 hours. Models are proposed to predict the 10.6 micrometer absorptivity for long exposure periods. Mechanical creep has been detected in both materials at high temperature.

Blaszuk, P. R.↗

Highly Supercooled Riming and Unusual Triple-Frequency Radar Signatures Over McMurdo Station, Antarctica

Riming of ice crystals by supercooled water droplets is an efficient ice growth process, but its basic properties are still poorly known. While it has been shown to contribute significantly to surface precipitation at mid-latitudes, little is known about its occurrence at high latitudes. In Antarctica, two competing effects can influence the occurrence of riming: (i) the scarcity of supercooled liquid water clouds due to the extremely low tropospheric temperatures and (ii) the low aerosol concentration, which may lead to the formation of fewer and larger supercooled drops potentially resulting in an enhanced riming efficiency. In this work, by exploiting the deployment of an unprecedented number of multiwavelength remote sensing systems (including triple-frequency radar measurements) in West Antarctica, during the Atmospheric Radiation Measurements West Antarctic Radiation Experiment (AWARE) field campaign, we evaluate the riming incidence at McMurdo Station and find that riming occurs at much lower temperatures when compared to previous results in the mid-latitudes. This suggests the possible occurrence of a common atmospheric state over Antarctica that includes a rather stable atmosphere inhibiting turbulent mixing, and a high riming efficiency driven by large cloud droplets. We then focus on a peculiar case study featuring a persistent layer with a particularly pronounced riming signature in triple-frequency radar data but only a relatively modest amount of supercooled liquid water. In-depth analysis of the radar observations suggests that such signatures can only be explained by the combined effects of moderately rimed aggregates or similarly shaped florid polycrystals and a narrow particle size distribution (PSD). Simulations of this case study performed with a one-dimensional bin model indicate that similar triple-frequency radar observations can be reproduced when narrow PSDs are simulated. Such narrow PSDs can in turn be explained by two key factors: (i) the presence of a shallow homogeneous droplet or humidified aerosol freezing layer aloft seeding an underlying supercooled liquid layer, and (ii) the absence of turbulent mixing throughout a stable polar atmosphere that sustains narrow PSDs, as hydrometeors grow from the nucleation region aloft to ice particles of several millimeters in size, by vapor deposition and then riming.

Riming↗

Impact Testing on Reinforced Carbon-Carbon Flat Panels with Ice Projectiles for the Space Shuttle Return to Flight Program

Following the tragedy of the Orbiter Columbia (STS-107) on February 1, 2003, a major effort commenced to develop a better understanding of debris impacts and their effect on the space shuttle subsystems. An initiative to develop and validate physics-based computer models to predict damage from such impacts was a fundamental component of this effort. To develop the models it was necessary to physically characterize reinforced carbon-carbon (RCC) along with ice and foam debris materials, which could shed on ascent and impact the orbiter RCC leading edges. The validated models enabled the launch system community to use the impact analysis software LS-DYNA (Livermore Software Technology Corp.) to predict damage by potential and actual impact events on the orbiter leading edge and nose cap thermal protection systems. Validation of the material models was done through a three-level approach: Level 1--fundamental tests to obtain independent static and dynamic constitutive model properties of materials of interest, Level 2--subcomponent impact tests to provide highly controlled impact test data for the correlation and validation of the models, and Level 3--full-scale orbiter leading-edge impact tests to establish the final level of confidence for the analysis methodology. This report discusses the Level 2 test program conducted in the NASA Glenn Research Center (GRC) Ballistic Impact Laboratory with ice projectile impact tests on flat RCC panels, and presents the data observed. The Level 2 testing consisted of 54 impact tests in the NASA GRC Ballistic Impact Laboratory on 6- by 6-in. and 6- by 12-in. flat plates of RCC and evaluated three types of debris projectiles: Single-crystal, polycrystal, and "soft" ice. These impact tests helped determine the level of damage generated in the RCC flat plates by each projectile and validated the use of the ice and RCC models for use in LS-DYNA.

Melis, Matthew E.↗

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.↗

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗

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↗

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↗

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.↗