Predictive Crystal Plasticity Modeling of Single Crystal Nickel Based on First-Principles Calculations
Not provided.
Materials Science↗
SEARCH · Search NASA
Results for “Single crystal model”
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.
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↗
Constitutive modeling for single crystal superalloys
A crystallographic approach to constitutive modeling of single crystal superalloys is discussed. The approach is based on identifying the active slip planes and slip directions. The shear stresses are computed on each of the slip planes from applied stress components. The slip rate is then computed on each slip system and the microscopic inelastic strain rates are the sum of the slip in the individual slip systems. The constitutive model was implemented in a finite element code using twenty noted isoparametric solid elements. Constants were determined for octahedral and cube slip systems. These constants were then used to predict tension-compression asymmetry and fatigue loops. Other data was used to model the tensile and creep response.
Stouffer, Donald C.↗
A viscoplastic model for single crystals
A viscoplastic constitutive model is described in which deformation behavior is postulated on representative slip systems and the behavior of the entire crystal is determined by summing the slip on the active slip systems. By building in the slip geometry known from the metallurgical literature, it is possible to predict the anisotropic deformation behavior and to model in a straightforward manner other phenomena which have been described by metallurgists in crystallographic terms. Elevated temperature tension-torsion tests were run and used to verify the model's predictive abilities. Ratchetting behavior under thermomechanical loading conditions is specifically addressed.
Jordan, E. H.↗
Unified constitutive model for single crystal deformation behavior with applications
Single crystal materials are being used in gas turbine airfoils and are candidates for other hot section components because of their increased temperature capabilities and resistance to thermal fatigue. Development of a constitutive model which assesses the inelastic behavior of these materials has been studied in 2 NASA programs: Life Prediction and Constitutive Models for Engine Hot Section Anisotropic Materials and Biaxial Constitutive Equation Development for Single Crystals. The model has been fit to a large body of constitutive data for single crystal PWA 1480 material. The model uses a unified approach for computing total inelastic strains (creep plus plasticity) on crystallographic slip systems reproducing observed directional and strain rate effects as a natural consequence of the summed slip system quantities. The model includes several of the effects that have been reported to influence deformation in single crystal materials, such as shear stress, latent hardening, and cross slip. The model is operational in a commercial Finite Element code and is being installed in a Boundary Element Method code.
Walker, K. P.↗
Constitutive and life modeling of single crystal blade alloys for root attachment analysis
Work to develop fatigue life prediction and constitutive models for uncoated attachment regions of single crystal gas turbine blades is described. At temperatures relevant to attachment regions, deformation is dominated by slip on crystallographic planes. However, fatigue crack initiation and early crack growth are not always observed to be crystallographic. The influence of natural occurring microporosity will be investigated by testing both hot isostatically pressed and conventionally cast PWA 1480 single crystal specimens. Several differnt specimen configurations and orientations relative to the natural crystal axes are being tested to investigate the influence of notch acuity and the material's anisotropy. Global and slip system stresses in the notched regions were determined from three dimensional stress analyses and will be used to develop fatigue life prediction models consistent with the observed lives and crack characteristics.
Meyer, T. G.↗
Constitutive modelling of single crystal and directionally solidified superalloys
Successful attempts were made to model the deformation behavior of nickel base superalloys to be used in gas turbine engines based on both a macroscopic constitutive model and a micromechanical formulation based on crystallographic slip theory. These models were programmed as FORTRAN subroutines, are currently being used to simulate thermomechanical loading predictions expected at the fatigue critical locations on a single crystal turbine blade. Such analyses form a natural precursor to the application of life prediction methods to gas turbine airfoils.
Walker, Kevin P.↗
Constitutive modelling of single crystal and directionally solidified superalloys
The trend towards improved engine efficiency and durability places increasing demands on materials that operate in the hot section of the gas turbine engine. These demands are being met by new coatings and materials such as single crystal and directionally solidified nickel-base superalloys which have greater creep/fatigue resistance at elevated temperatures and reduced susceptibility to grain boundary creep, corrosion and oxidation than conventionally cast alloys. Work carried out as part of a research program aimed at the development of constitutive equations to describe the elevated temperature stress-strain-time behavior of single crystal and directionally solidified turbine blade superalloys is discussed. The program involves both development of suitable constitutive models and their verification through elevated temperature tension-torsion testing of single crystals of PWA 1480.
Jordan, E. H.↗
Electro-chemo-mechanical finite-element model of single-crystal and polycrystalline NMC cathode particles embedded in an argyrodite solid electrolyte
Not Available
Taghikhani, Kasra (ORCID:0000000277638315)↗
Modeling Single-Crystal Battery Materials: From Fundamental Understanding to Performance Evaluation
The performance of rechargeable batteries is fundamentally influenced by the physicochemical properties and microstructural features of their key material components. Recent experimental advancements have highlighted the potential of single-crystal (SC) morphologies to address inherent limitations of polycrystalline (PC) electrodes and solid-state electrolytes, offering tunable charge transport kinetics and improved cell cycling performance. Here, this review examines how state-of-the-art computational modeling, from atomistic and mesoscale to continuum-level approaches, including machine learning methodologies, has been utilized to investigate the critical factors governing the electrochemical behavior of SC battery materials. We explore how predictive modeling can elucidate the processing–structure–property–performance relationships of SC cathodes, anodes, and solid-state electrolytes, with a focus on unique SC characteristics such as crystallographic anisotropy, size effects, and facet-dependent properties. Additionally, we identify limitations in commonly used modeling techniques and discuss strategies to address these challenges. By integrating high-fidelity simulations with experimental insights, this review aims to outline a clear path for the rational design and optimization of SC battery components, paving the way for accelerated advancements in energy storage technologies.
Materials science↗
Constitutive modeling for single crystal superalloys
The inelastic response of single crystal gamma/gamma prime superalloys is quite different from the behavior of polycrystalline nickel base superalloys. Upto a critical temperature the yield stress of single crystal alloys is a function of the material orientation relative to the direction of the applied stress and the material exhibits significant tension/compression asymmetry. This behavior is primarily due to slip on the octahedral slip system. Above the critical temperature there is a sharp drop in the yield stress, cube slip becomes more predominant and the tension/compression asymmetry is reduced. Similar orientation and tension/compression asymmetry is observed in creep and secondary creep above the critical temperature is inferred to occur by octahedral slip. There are two exceptions to this behavior. First, loading near the (111) orientation exhibits cube slip at all temperatures, and; second, loading near the (001) orientation produces only octahedral slip at all temperatures. The constitutive model is based on separating the total global strain into elastic and inelastic components. This model is developed and briefly discussed.
Stouffer, D. C.↗
Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX
Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Fatigue damage modeling for coated single crystal superalloys
A high temperature, low-cycle fatigue life prediction method for coated single crystal nickel-base superalloys is being developed. The method is being developed for use in predicting crack initiation life of coated single crystal turbine airfoils. Although the models are being developed using coated single crystal PWA 1480, they should be readily adaptable to other coated nickel-base single crystal materials. The coatings choosen for this effort were of two generic types: a low pressure plasma sprayed NiCoCrAlY overlay, designated PWA 286, and an aluminide diffusion, designated PWA 273. In order to predict the useful crack initiation life of airfoils, the constitutive and failure behavior of the coating/substrate combination must be taken into account. Coatings alter the airfoil surface microstructure and are a primary source from which cracks originate. The adopted life prediction approach addresses this complexity by separating the coating and single crystal crack initiation regimes. This provides a flexible means for using different life model formulations for the coating and single crystal materials. At the completion of this program, all constitutive and life model formulations will be available in equation form and as software. The software will use the MARC general purpose finite element code to drive the constitutive models and calculate life parameters.
Nissley, David M.↗
Fully implicit crystal plasticity models representing orientations with modified Rodrigues parameters
Here, this work describes a crystal plasticity formulation combining several mathematical, numerical, and implementation choices to produce a highly efficient model. Specifically, the key choices in the implementation are (1) representing orientations with modified Rodrigues parameters, (2) implementing a fully coupled implicit time integration for the elastic stretch, the crystal orientations, and the model internal variables, (3) implementing the model in the NEML2 constitutive modeling framework, based on PyTorch, to vectorize the calculations and port the computation to GPUs and other hardware accelerators, and (4) an exact implementation of the consistent tangent matrix, even for arbitrary coupling to other field variables beyond the displacements, like temperature, neutron fluence, etc. The first two features of the model are, to our knowledge, novel. The paper considers each of these choices individually as well as the final model as a whole. This includes a full description of modified Rodrigues parameters, their advantages over other representations of orientations, the mathematical formulae and tools required to implement a model with modified Rodrigues parameters, and a detailed description of the geometry of the space of modified Rodrigues parameters (in an appendix). It also includes a description of a fully implicit time integration scheme for the orientations and the advantages in representing orientations with modified Rodrigues parameters in implementing such a model. The work then assess, via numerical examples, the advantages of fully coupled implicit time integration versus more common decoupled and explicit time integration schemes. These studies demonstrate the computational advantages of fully coupled integration versus other time integration algorithms, though the performance of the competing models depends on the complexity of the underlying single crystal model. The study concludes by demonstrating that the choice of time integration method affects the sharpness of the predicted texture, with explicit methods for integrating the orientations overestimating texture sharpness and implicit methods underestimating texture sharpness.
Crystal plasticity↗
Finite element (MARC) solution technologies for viscoplastic analyses
A need for development of realistic constitutive models for structural components operating at high temperatures, accompanied by appropriate solution technologies for stress/life analyses of these components is studied. Viscoplastic models provide a better description of inelastic behavior of materials, but their mathematical structure is very complex. The highly nonlinear and stiff nature of the constitutive equations makes analytical solutions difficult. Therefore, suitable solution, finite element or other numerical, technologies must be developed to make these models adaptable for better and rational designs of components. NASA-Lewis has developed several solution technologies and successfully applied them to the solution of a number of uniaxial and multiaxial problems. Some of these solution technologies are described along with the models and representative results. The solution technologies developed and presented encompass a wide range of models, such as, isotropic, anisotropic, metal matrix composites, and single crystal models.
Arya, V. K.↗
Anisotropic constitutive modeling for nickel base single crystal superalloy Rene N4 at 982 C
A back stress/drag stress constitutive model based on a crystallographic approach to model single crystal anisotropy is presented. Experimental results demonstrated the need for the back stress variable in the inelastic flow equations. Experimental findings suggested that back stress is orientation dependent and controls both strain hardening and recovery characteristics. Due to the observed stable fatigue loops at 1800 F, drag stress is considered constant for this temperature. The constitutive model operated with constraints determined only from tensile data was extensively tested from simple tensile and fatigue to complicated strain hold tests. The model predicted very well under those conditions.
Stouffer, D. C.↗
High dopant activation in arsenic doped single-crystal CdTe thin films: Insights from MBE growth and rapid thermal processing
Single-crystal model systems are valuable tools to investigate fundamental material properties. In this work, we use molecular beam epitaxy to deposit in situ arsenic (As) doped single-crystal CdTe films on large area Si substrates to better understand As doping for photovoltaic applications. We found that As incorporation is highly temperature dependent: a substrate temperature difference of 50 °C can lead to several orders of magnitude difference in As concentration. Cd overpressure during in situ doping may limit out-diffusion of As but decrease As incorporation, especially at lower growth temperatures. Carrier concentrations greater than 10 16 cm −3 can be achieved with or without Cd overpressure when annealed at temperatures above 500 °C. However, unlike the low (∼1% to 5%) dopant activation commonly observed in polycrystalline CdTe, our films achieve significantly higher activation ratios—exceeding 50%, and in some cases approaching 80%. These values are consistent with or exceed prior reports in single-crystal CdTe systems. In addition to as-deposited arsenic concentrations, we also consider arsenic distribution after different rapid thermal processing temperatures. We propose a detailed definition and description of how arsenic incorporation is considered and calculated. Due to carrier concentration saturation, As incorporation also needs to be controlled to average levels of 10 17 cm −3 to achieve high activation. These findings suggest that higher annealing temperature regimes may be beneficial to polycrystalline CdTe based PV devices.
36 MATERIALS SCIENCE↗