The impact of water on slip system activity in olivine and the formation of bimodal crystallographic preferred orientations
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Dislocations are line defects in crystals that multiply and self-organize into a complex network during strain hardening. The length of dislocation links, connecting neighboring nodes within this network, contains crucial information about the evolving dislocation microstructure. By analyzing data from Discrete Dislocation Dynamics (DDD) simulations in face-centered cubic (fcc) Cu, we characterize the statistical distribution of link lengths of dislocation networks during strain hardening on individual slip systems. Here, our analysis reveals that link lengths on active slip systems follow a double-exponential distribution, while those on inactive slip systems conform to a single-exponential distribution. The distinctive long tail observed in the double-exponential distribution is attributed to the stress-induced bowing out of long links on active slip systems, a feature that disappears upon removal of the applied stress. We further demonstrate that both observed link length distributions can be explained by extending a one-dimensional Poisson process to include different growth functions. Specifically, the double-exponential distribution emerges when the growth rate for links exceeding a critical length becomes super-linear, which aligns with the physical phenomenon of long links bowing out under stress. This work advances our understanding of dislocation microstructure evolution during strain hardening and elucidates the underlying physical mechanisms governing its formation.
Single crystal superalloys have become increasingly popular for turbine blade and vane applications due to their high strength, and creep and fatigue resistance at elevated temperatures. The crystallographic orientation of a single crystal material greatly affects its material properties, including elastic modulus, shear modulus, and ductility. These directional properties, along with the type of loading and temperature, dictate an anisotropic response in the yield strength, creep resistance, creep rupture ductility, fatigue resistance, etc. A significant amount of research has been conducted to determine the material properties in the <001> orientation, yet the material properties deviating from the <001> orientation have not been assessed for all cases. Based on the desired application and design criteria, a crystal orientation is selected to yield the maximum properties. Currently, single crystal manufacturing is able to control the primary crystallographic orientation within 15 of the target orientation, which is an acceptable deviation to meet both performance and cost guidelines; the secondary orientation is rarely specified. A common experiment is the standard load-controlled tensile test, in which specimens with different orientations can be loaded to observe the material response. The deformation behavior of single-crystal materials under tension and compression is known to be a function of not only material orientation, but also of varying microdeformation (i.e. dislocation) mechanisms. The underlying dislocation motion causes deformation via slip, and affects the activation of specific slip systems based on load and orientation. The slip can be analyzed by observing the visible traces left on the surface of the specimen from the slip activity within the single crystal material. The goal of this thesis was to predict the slip systems activated in three-dimensional stress fields of a notched tensile specimen, as a function of crystal orientation, using finite element analysis without addressing microstructural deformation mechanisms that govern their activation. Out of three orientations tested, the specimen with a [110] load orientation and a [001] growth direction had the lowest maximum resolved shear stress; this specimen orientation appears to be the best design candidate for a tensile application.
FeCrAl cladding is one of the candidate materials for the near‐term accident‐tolerant fuel technologies under development. Research on high‐temperature mechanical behaviors of single‐crystal FeCrAl alloy is rather limited. Previous studies have reported the mechanical property of low‐index orientation in single‐crystal FeCrAl alloy at room temperature. However, the critical resolved shear stress to activate slip systems can be orientation and temperature dependent. Here, single‐crystal grains in a coarse‐grained FeCrAl alloy with different crystallographic orientations are selected to preferentially activate {110}<111> slip systems or {112}<111> slip systems. Micropillars are fabricated in the selected single‐crystal grains and tested at elevated temperatures in situ in a scanning electron microscope. The critical resolved shear stresses of {110}<111> slip systems and {112}<111> slip systems are determined at various temperatures. The critical resolved shear stress shows a temperature dependence and orientation independence. This study provides important insight for understanding the deformation mechanisms of FeCrAl alloys at elevated temperatures.
Fatigue crack growth tests were carried out on large-grain Al 7029 aluminum alloy and the finite element method was used to calculate the stress field near the tip of a zigzag crack. The resolved shear stresses on all 12 slip systems were computed, and the resolved shear stress intensity coefficient (RSSIC) was defined. The RSSIC was used to analyze the irregular crack path and was correlated with the rate of single-slip-plane shear crack growth. Fatigue crack growth was found to be caused primarily by shear decohesion at a crack tip. When the RSSIC on a single-slip system was much larger than all the others, the crack followed a single-slip plane. When the RSSICs on two conjugate slip systems were comparable, a crack grew in a zigzag manner on these planes and the macrocrack-plane bisected the two active slip planes. The maximum RSSIC on the most active slip system is proposed as a parameter to correlate with the shear fatigue crack growth rate in large 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.
Nonequilibrium molecular dynamics simulations of shock loaded single-crystal Al in the $\langle$100$\rangle$, $\langle$110$\rangle$, $\langle$111$\rangle$, and $\langle$123$\rangle$ orientations are conducted to study elastic and plastic shockwave formation and details associated with dislocation activity. A computer vision-based approach is implemented to capture the presence of dislocations and describe their spatial characteristics in the zone of nucleation behind the propagating shockwave. The methodology developed relies on the sequences of images extracted during shock loading that show dislocation activity within a cross section of the sample. Results reveal that the spacing between activated slip systems is orientation dependent and exhibits a modest reduction for the $\langle$100$\rangle$ and $\langle$111$\rangle$ orientations as shock pressure increases. Comparisons are made to existing theoretical models. Such relationships between shock pressure and dislocation activity, extracted from molecular dynamics simulations, can be used to inform higher length scale simulations or modeling of dislocation-based plasticity during shock.
Twinning is a key deformation mechanism in hexagonal close-packed (hcp) metals, which are typified by a lack of easily activated slip systems that can accommodate a general state of loading. Pragmatically, the nucleation and evolution of twin domains occurs concomitantly with slip, such that the eventual twin network is conditioned by both external and internal stresses resulting, among others, from the evolution of dislocations. However our understanding of the interplay between dislocations, and twin nucleation and stability remains limited. This work focuses on elucidating the influence of dislocation-mediated plasticity on the formation, growth, and stability of twin embryos. First, a new mesoscale spectral crystal plasticity-twinning framework is developed and used to quantify the change in the free energy landscape following the nucleation {$10\bar{1}2$} of twins in Mg. Representative twin morphologies are modeled under conditions of limited and profuse slip activity, which are emulative of small- and bulk-scale samples, respectively. Then, the driving traction profiles around twin embryos are investigated via a sharp interface approach to obtain insights into how concurrent slip-mediated plasticity can influence the growth/stabilization of nanometric twins. The driving traction profiles are further utilized to determine the stability of twin embryos post loading. The initial dislocation density in the samples, and within the different domains (i.e., twin vs. parent), is seen to have a significant effect on the twin nucleation stress. Namely, the activation of high levels of concomitant plasticity in the parent grain is seen to significantly drive the formation of nanometric twin nuclei at stresses as low as 250 MPa. Further, the sharp interface analysis reveals that profuse plasticity in the parent grain simultaneously alters the forward and back stresses, such that the magnitude/polarity of the driving tractions become increasingly favorable for nanometric twin growth when slip is active. Finally, prior plasticity in the parent grain is seen to result in favorable driving tractions for nanometric twin growth, even at applied stresses as low as ~ 100 MPa. In conclusion, these results are in stark contrast to a case without any dislocations, wherein applied stresses as high as ~ 500 MPa are necessary to grow the twin domains.
The work focuses on the determination of the critical resolved shear stress (CRSS) in titanium (Ti) and titanium-aluminum (Ti-Al) alloys, influenced by an array of factors such as non-symmetric fault energies and minimum energy paths, dislocation core-widths, short-range order (SRO) effects which alter the local atomic environment, and tension-compression (T-C) asymmetry affected by intermittent slip motion. To address these multifaceted complexities, an advanced theory has been developed, offering an in-depth understanding of the mechanisms underlying slip behavior. The active slip systems in these materials are basal, prismatic, and pyramidal planes, with the latter involving both $\langle$a$\rangle$ and $\langle$c + a$\rangle$ dislocations. Each slip system is characterized by distinct Wigner-Seitz cell configurations for misfit energy calculations, varying partial dislocation separation distances, and unique dislocation trajectories—all critical to precise CRSS calculations. The theoretical CRSS results were validated against a comprehensive range of experimental data, demonstrating a strong agreement and underscoring the model's efficacy.
Damage evolution in engineering metal alloys at the grain scale exhibits significant microstructural heterogeneity and anisotropy. These heterogeneities create local hotspots for stress and strain localization, leading to void nucleation. Crystal orientation influences the active slip systems around voids, affecting lattice rotation and potentially forming discontinuities. At low triaxiality, voids may change shape due to lower stress, rotation, elongation, and coalescence. At high triaxiality, the correlation between crystal orientation and void growth rate becomes stronger, resembling the behavior observed in isolated single crystals. Therefore, understanding the effects of crystal orientation, heterogeneous strain, and defect evolution is crucial for single crystal fracture characterization. Here, in this work, a coupled phase-field damage (PFD) and crystal plasticity (CP) model is implemented within a finite element framework to analyze crystal deformation and failure. The CP method employs a dislocation density-based constitutive model, while intragranular failure is modeled using an anisotropic PFD method. The PFD model considers both the stored energy due to elastic stretching and the energy release due to defect formation and crack formation. A single crystal Al2219 with an intracrystalline spherical void is chosen to analyze fracture. The study finds that fracture propagation is strongly correlated with crystal orientations. This coupled CP-PFD model provides accurate failure prediction in crystalline materials by incorporating the effects of crystal orientations and existing voids. This study demonstrates how the local microstructure and defects influence plastic deformation and failure mechanisms in metal alloys.
Additive manufacturing (AM) potentially enables fabrication and repair of plasma facing components (PFCs) of tungsten. However, deformation mechanisms in AM-fabricated tungsten under service-relevant thermomechanical conditions remain unknown majorly due to challenges in achieving crack-free W. This study reveals the deformation mechanisms operative under tension at 800 °C and 1200 °C in the electron beam melting powder bed fusion (EBM-PBF) fabricated W. Textured columnar grains with mixed <001 >/<111 >|| build-direction were observed. Strong anisotropy in tensile properties persisted across test temperatures, and all specimens exhibited extensive deformation-induced banding phenomena. Grain boundary character analysis of the deformed specimens indicated that these deformation bands were kink bands with tilt grain boundaries prominently present at the band-matrix interface. Interestingly, the material within the kink bands rotated toward higher resolved shear stress, providing insights into the origins of kink band formation in body-centered cubic refractory metals. Intragranular misorientation axis analysis revealed that plasticity was dominated by {110} <111 > slip at 800 °C, whereas additional (213)[11̄1] and (112)[11̄1] slip systems activated at 1200 °C. A dense low-angle boundary networks observed near fractured surface indicated the onset of dynamic recrystallization. Results reveal plasticity governing mechanisms in AM-fabricated W under power plant-relevant conditions, and provide insights into kink band attributes in refractory metals.
A study of deformation mechanisms in the single crystal superalloy PWA 1480 was conducted. Monotonic and cyclic tests were conducted from 20 to 1093 C. Both (001) and near-(123) crystals were tested, at strain rates of 0.5 and 50 percent/minute. The deformation behavior could be grouped into two temperature regimes: low temperatures, below 760 C; and high temperatures, above 820 to 950 C depending on the strain rate. At low temperatures, the mechanical behavior was very anisotropic. An orientation dependent CRSS, a tension-compression asymmetry, and anisotropic strain hardening were all observed. The material was deformed by planar octahedral slip. The anisotropic properties were correlated with the ease of cube cross-slip, as well as the number of active slip systems. At high temperatures, the material was isotropic, and deformed by homogeneous gamma by-pass. It was found that the temperature dependence of the formation of superlattice-intrinsic stacking faults was responsible for the local minimum in the CRSS of this alloy at 400 C. It was proposed that the cube cross-slip process must be reversible. This was used to explain the reversible tension-compression asymmetry, and was used to study models of cross-slip. As a result, the cross-slip model proposed by Paidar, Pope and Vitek was found to be consistent with the proposed slip reversibility. The results were related to anisotropic viscoplastic constitutive models. The model proposed by Walter and Jordan was found to be capable of modeling all aspects of the material anisotropy. Temperature and strain rate boundaries for the model were proposed, and guidelines for numerical experiments were proposed.
Not Available
The microstructural characteristics and morphology of interfaces in metals can be crucial in governing the initiation of plasticity and early deformation mechanisms under extreme conditions. During high-strain-rate deformation, these interfaces significantly affect dislocation nucleation, twinning, and other mechanisms that directly impact material strength and failure. Despite their importance, a substantial knowledge gap remains between the observed macroscopic material behavior and the underlying role of bimetal interfaces in plasticity initiation. Here, to address this gap, large-scale molecular dynamics simulations are performed on Al/Ti bimetal structures to examine the effect of interface characteristics on the onset of plasticity under uniaxial compression. Specifically, this study investigates how interface roughness (flat vs waveform interfaces) modifies the initiation of plastic events in the microstructure. Atomistic simulations indicate that interface roughness (a microscopic behavior) reduces the stress required for dislocation nucleation, thereby reducing the peak stress relative to a flat interface. For the square interface, the step height strongly influences plasticity initiation by setting the separation of locally flat regions. Varying the interface rotation relative to the loading direction (macroscopic behavior) reveals that the peak stress for both flat and waveform interfaces initially decreases and then increases with rotation, accompanied by a shift in dominant mechanisms—from Al twinning (flat) to interface sliding and Ti-dominated phase transformation and twinning at lower angles. The change in the stress-strain slope during initial compression reflects how rotation alters the resolved shear stress and activates different slip systems. Overall, the study provides valuable insights into the role of bimetal interfaces in controlling plasticity initiation and early deformation pathways for Al/Ti under extreme loading conditions.
Metals and their alloys, except for a few intermetallics, are inherently ductile, i.e. plastic deformation precedes fracture in these materials. Therefore, resistance to fracture is directly related to the development of the plastic zone at the crack tip. Recent studies indicate that the fracture toughness of single crystals depends on the crystallographic orientation of the notch as well as the loading direction. In general, the dependence of crack propagation resistance on crystallographic orientation arises from the anisotropy of (i) elastic constants, (ii) plastic deformation (or slip), and (iii) the weakest fracture planes (e.g. cleavage planes). Because of the triaxial stress state at the notch tips, many slip systems that otherwise would not be activated during uniaxial testing, become operational. The plastic zone formation in single crystals has been tackled theoretically by Rice and his co-workers and only limited experimental work has been conducted in this area. The study of the stresses and strains in the vicinity of a FCC single crystal notch tip is of relatively recent origin. We present experimental and numerical investigation of 3D stress fields and evolution of slip sector boundaries near notches in FCC single crystal tension test specimens, and demonstrate that a 3D linear elastic finite element model that includes the effect of material anisotropy is shown to predict active slip planes and sectors accurately. The slip sector boundaries are shown to have complex curved shapes with several slip systems active simultaneously near the notch. Results are presented for surface and mid-plane of the specimens. The results demonstrate that accounting for 3D elastic anisotropy is very important for accurate prediction of slip activation near FCC single crystal notches loaded in tension. Results from the study will help establish guidelines for fatigue damage near single crystal notches.
Parallel experiments have been performed in order to develop a comprehensive model for stress cracking (SCC) in structural materials. The central objective is to determine the relationship between the activity and selectivity of the microstructure of structural materials to their dissolution kinetics and experimentally measured SCC kinetics. Zinc was chosen as a prototype metal system. The SCC behavior of two oriented single-crystal disks of zinc in a chromic oxide/sodium sulfate solution (Palmerton solution) were determined. It was found that: (1) the dissolution rate is strongly (hkil)-dependent and proportional to the exposure time in the aggressive environment; and (2) a specific slip system is selectively active to dissolution under applied stress and this slip line controls crack initiation and propagation. As a precursor to potential microgrvity experiments, electrophoretic mobility measurements of zinc particles were obtained in solutions of sodium sulfate (0.0033 M) with concentrations of dissolved oxygen from 2 to 8 ppm. The equilibrium distribution of exposed oriented planes as well as their correlation will determine the particle mobility.
Nickel-based single-crystal (Ni-SX) superalloys under cyclic stress are susceptible to cracking at stress-concentration sites, eventually leading to low-cycle fatigue (LCF) failure. LCF cracks typically originate from intrinsic defects (e.g., voids and carbides) within solidified dendrites. However, systematic quantitative experimental analyses of defect-mediated local damage remain limited. To thoroughly understand the microscopic origins and evolution of LCF damage, correlated 3D mapping of dendrites across various regions is essential. Here, in this study, macroscale micro-computed tomography (μ-CT) was initially used to capture internal interdendritic secondary cracks within bulk DD413 superalloy after LCF testing at 760 °C. Subsequently, a multimodal methodology combining synchrotron 3D microdiffraction (3D-μXRD), high-resolution μ-CT, and electron microscopy was established. This approach allowed precise localization of internal damage zones near interdendritic secondary cracks and detailed mapping of the 3D correlated distributions of dendrites, defects, and residual stress/strain fields within these zones at submicron spatial resolution. Finally, the same approach was applied to specimens subjected to interrupted loading at approximately 40 % of the fatigue life to uncover the early damage states of dendrites. The dendrite cores (DCs) and interdendritic regions (IDs) exhibit microscale heterogeneous mechanical responses: nearly defect-free DCs accumulate local irreversible slip along specific slip systems to generate slip bands, while the IDs containing various defects accommodate local microplasticity through the activation of multiple slip systems around these defects. The local tensile stress near defects in the IDs exceeds that in the DC slip band regions by more than threefold, leading to the generation of local damage zones within the IDs. Chain-like defect distributions facilitate the interconnection of these local zones into a continuous damage region, further elevating the overall tensile stress in the IDs. Additionally, geometrically necessary dislocations alone are insufficient as indicators of LCF damage; both the internal stress state and its magnitude must be considered. These experimental results provide critical data and insights for the development of multi-physics fatigue models.
Various methods for calibrating crystal plasticity finite element (CPFE) models lead to non-unique input parameter values, which subsequently introduce uncertainty in the predicted mechanical response. Sensitivity analysis (SA) conducted on crystal plasticity models is used to identify how variability in these parameters contribute to output uncertainty. Traditional SA on CPFE parameters uses simplified surrogate models to save computational time. However, the accuracy of the surrogate models depends on the quantity of training data used, and any modeling error can propagate into the SA results, potentially affecting their reliability. In this work, the elementary effects test (EET) method, a global SA technique using direct CPFE simulations was employed, and the results obtained were compared with the First Order Second Moment (FOSM) method. ExaConstit, an open-source GPU-enabled CPFE code, was used to perform the simulations and direct SA. The EET method was accurately able to capture the non-linear effects of all the input parameters on the output and is a valuable approach for reliably attributing parameter sensitivities in CPFE models. Based on the results, efficient strategies to perform future parameter calibration and SA are discussed. Additionally, the SA trends observed in different single crystal orientations closely mirrored the activity of the slip systems.