Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dislocation mobility”

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 55 records · Page 3

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

Helium causing disappearance of a/2<111> dislocation loops in binary Fe-Cr ferritic alloys

In this work, single and dual-beam self-ion irradiations were performed at 500°C on ultra-high purity Fe14%Cr alloy to ~0.33 displacements-per-atom (dpa), and 0 or 3030 atomic-parts-per-million (appm) helium/dpa, respectively. Using transmission electron microscopy, we reveal that helium can drastically modify the dislocation loop Burgers vector in Fe-Cr alloys. Helium co-implantation caused complete disappearance of a/2<111> type dislocation loops, and the microstructure consisted of only a<100> loops. Conversely, a/2<111> type loops were predominant without He co-implantation. The total loop density remained largely unaffected. The results strikingly contrast literature asserting that helium stabilizes a/2<111> type loops in bcc Fe alloys, based on low temperature irradiations. Collectively analyzing the results with literature suggest that the small positive interaction between helium and self-interstitial atoms (SIA) in Fe predicted by atomistic simulations maybe insufficient to holistically explain the dislocation loop microstructure development in presence of helium. Helium-SIA positive binding inadvertently implies elevated a/2<111> loop fraction and higher loop densities that the present results contradict. Helium induced high cavity density causing a preferential loss of highly glissile <111> clusters, leaving the matrix saturated with <100> type clusters is proposed as a potential mechanism. Further, the in-situ irradiations combined with Burgers vector analysis strengthened the evidence of Cr-induced dislocation loop mobility reduction that appears to stabilize the a/2<111> type loops and causes higher loop densities in Fe-Cr alloys.

36 MATERIALS SCIENCE↗

Chromium enhances the mechanical performance of 3d transition metal high entropy alloys

A comprehensive examination of the compositional effects on the deformation behavior of high entropy alloys (HEAs) was conducted through nanoindentation, indentation creep, and stress relaxation experiments at ambient temperature. The alloys investigated included NiCoFe, NiCoCr, NiCoCrFe, and NiCoCrFeMn, all of which underwent identical thermomechanical processing. From the experimental results, NiCoCr exhibited the highest maximum shear stress for dislocation nucleation (13 GPa) and nanoindentation hardness (3.8 GPa) among the alloys tested. It also showed the lowest stress sensitivity (0.02) and activation volume (3–4 b 3 ) under creep and stress relaxation conditions, indicative of its superior plastic flow properties. Density functional theory (DFT) calculations further revealed an uneven charge density and larger bonding-length variation in NiCoCr due to the presence of Cr, which increased lattice distortion, impeding dislocation movement. Furthermore, in this study, post-deformation microscopy using transmission electron microscopy (TEM) revealed a high density of stacking faults and twins in NiCoCr that effectively enhanced its creep strength. The results highlight the significance of specific elemental effects, particularly from Cr in this study, over configurational entropy effect (i.e., compositional complexity) in governing the deformation microstructure and mechanical properties of HEAs. Finally, these insights will be instrumental for the design and optimization of advanced alloys for load-bearing applications.

36 - MATERIALS SCIENCE↗

Prismatic-to-basal plastic slip transition in zirconium

Here, we model three-dimensional dislocation ensembles in Zr pillars using newly developed mobility laws for dislocations on prismatic and basal planes. The effects of the loading orientation and temperature on the micro-scale mechanical response of single crystals are systematically explored. Easy dislocation glide is observed on prismatic planes, while hard glide occurs on basal planes. By selection of sample temperature and loading orientation, we show that the transition from prismatic to basal glide of dislocation ensembles can be controlled. We also show that the presence of dislocation cross-slip strongly promotes prismatic glide as a result of prismatic/basal cross-slip energy asymmetry. Cross-slip does not alter the occurrence of plastic slip transition but only induces a shift towards higher temperatures. The preponderance of plastic slip on basal or prismatic planes is found to be mediated by a transitional dislocation, which is composed of glissile segments on parallel prismatic planes connected by glissile super-jogs on basal planes. These findings provide a mechanistic understanding of experiments and highlight the significance of transitional prismatic-basal dislocations on the macroscopic characteristics of plasticity in HCP crystals.

36 MATERIALS SCIENCE↗

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The effects of stress, temperature and facet structure on growth of { 10 1 &#xAF; 2 } twins in Mg: a molecular dynamics and phase field study

The present study focuses on the roles of structure, stress and temperature on the mobility of facets bounding twin domains and on the overall impact on the growth kinetics of { 10 1 &#xAF; 2 } twins in Mg. We first use molecular dynamics simulations (MD) to study the relationship between interface structures, stress, temperature and mobility. As facets can exhibit distinct atomic scale structures, we also quantify the effect of stress relaxation mediated by misfit dislocations on the mobility of the basal-prismatic interface. The information gathered for ten different facets is used to calibrate an anisotropic phase field model predicting the effects of temperature and stress on the kinetics and shape of a growing { 10 1 &#xAF; 2 } 3D twin, by comparison with MD data. The comparative analysis suggests that, as compared to simulations of individual facets, the effective mobility of facets can be significantly different during growth of a 3D twin. This can be attributed to the complex three-dimensional internal stress state of the twin and to coupling in facet motion triggered by facet junctions. As such, by directly fitting both the twin growth kinetics and morphology as predicted by MD for different stresses and temperatures, we can identify a limited set of rate-limiting facets that control twin growth and morphology, and measure their effective mobility.

36 MATERIALS SCIENCE↗

Incorporating point defect generation due to jog formation into the vector density-based continuum dislocation dynamics approach

During plastic deformation of crystalline materials, point defects such as vacancies and interstitials are generated by jogs on moving dislocations. A detailed model for jog formation and transport during plastic deformation was developed within the vector density-based continuum dislocation dynamics framework. Here as a part of this model, point defect generation associated with jog transport was formulated in terms of the volume change due to the non-conservative motion of jogs. Balance equations for the vacancies and interstitials including their rate of generation due to jog transport were also formulated. A two-way coupling between point defects and dislocation dynamics was then completed by including the stress contributed by the eigen-strain of point defects. A jog drag stress was further introduced into the mobility law of dislocations to account for the energy dissipation during point defects generation. A number of test problems and a fully coupled simulation of dislocation dynamics and point defects generation and diffusion were performed. The results show that there is an asymmetry of vacancy and interstitial generation due to the different formation energies of the two types of defects. The results also show that a higher hardening rate and a higher dislocation density are obtained when the point defect generation mechanism is coupled to dislocation dynamics.

36 MATERIALS SCIENCE↗

Simulation of dislocation evolution in microparticle impacts over a wide range of impact velocities

We report the dynamic loading environment created by high velocity microparticle impacts is very difficult to model. The extremely high strain rates, large deformations, and high temperatures affect the flow stress of the material in ways that standard flow stress models cannot capture. Plastic deformation and strain hardening in metals is controlled by the motion of dislocations. Dislocations can be nucleated, stored as forest dislocations, or be annihilated as loading progresses. A comprehensive accounting of dislocation density change is needed to accurately describe rate and temperature dependent dislocation glide and evolution across the wide rage of loading conditions present in microparticle impact problems. Therefore, we implement and apply a newly proposed flow stress model (Hunter and Preston, 2015, 2022) that is then coupled with mobile and immobile dislocation density evolution equations. This model is implemented in Los Alamos National Laboratory’s hydrodynamics code, FLAG, to model copper-on-copper microparticle impacts. This new strength model allows for accurate simulation of particle rebound and flattening across a wide range of impact velocities.

36 MATERIALS SCIENCE↗

Mechanical properties and deformation mechanisms of single crystal Mg micropillars subjected to high-strain-rate C-axis compression

Here, the mechanical properties and deformation mechanisms of single crystal magnesium under c-axis quasi-static and high-strain rate compressions are investigated through in situ scanning electron microscope (SEM) experiments and post-mortem transmission electron microscope (TEM) characterization. The findings revealed that ductility and high rates of hardening are preserved for pillars as large as 15 μm. Furthermore, rate effects result in a mild increase in flow stress with plastic deformations controlled primarily by the slip of type dislocations. Importantly and in contrast to other literature reports, plastic deformation occurs in the absence of twining. As the strain increases and plastic deformation exceeds about 4%, crystal rotation activates basal slip, <$\mathrm{a}$> type dislocations, resulting in a more rate independent flow stress. TEM observation on micropillars compressed at a strain rate of 250/s, revealed the activation of {${11}$$\bar{2}$$\bar{2}$} < $\bar{1}$$\bar{1}23$ > slip systems and high mobility of screw dislocations as major contributors to plastic strains in excess of 10% without fracture. These findings are relevant to the design of lightweight materials used in transportation systems, e.g., selection of material grain size. Moreover, the experimental data here reported provides the materials science community with a unique opportunity to validate discrete dislocation dynamics (DDD) formulations employed in multiscale design of materials.

36 MATERIALS SCIENCE↗

In situ TEM annealing of neutron-irradiated Ti reveals a two-stage mechanism for elevated temperature radiation damage recovery

Understanding how irradiation-induced defects evolve at elevated temperatures is of critical importance to predicting materials' behavior under steady-state and accident scenarios. However, such mechanistic insight into microstructural evolution is limited by the nature of ex situ annealing and subsequent imaging. Here we show direct observation and quantification of defect recovery in neutron-irradiated Ti using in situ transmission electron microscopy (TEM) annealing experiments. In agreement with our prior work, and at temperatures below the irradiation temperature (T irr = 300 °C), dislocation loops are observed to glide. At elevated temperatures (>500 °C), dislocation lines become mobile and promote significant recovery of the microstructure. These mechanisms challenge the established electron irradiation-based model for radiation damage recovery, which originally suggests dissolution of static defect clusters, and demonstrates the importance of in situ characterization in understanding defect evolution in irradiated materials.

Hirst, Charles A. [Univ. of Michigan, Ann Arbor, M↗

Exploring the relation between transonic dislocation glide and stacking fault width in FCC metals

Theory predicts limiting gliding velocities that dislocations cannot overcome. Computational and recent experiments have shown that these limiting velocities are soft barriers and dislocations can reach transonic speeds in high rate plastic deformation scenarios. In this paper we systematically examine the mobility of edge and screw dislocations in several face centered cubic (FCC) metals (Al, Au, Pt, and Ni) in the extreme large-applied-stress regime using molecular dynamics simulations. Our results show that edge dislocations are more likely to move at transonic velocities due to their high mobility and lower limiting velocity than screw dislocations. Importantly, among the considered FCC metals, the dislocation core structure determines the dislocation’s ability to reach transonic velocities. This is likely due to the variation in stacking fault width due to relativistic effects near the limiting velocities.

36 MATERIALS SCIENCE↗

Validating continuum theory for Cottrell atmosphere solute drag by molecular dynamics simulations

When a dislocation moves through a field of mobile solute atoms, solutes segregate to the dislocation and form Cottrell atmospheres which exert solute drag forces. Over the last 70 years, continuum theory has been used extensively to estimate these drag forces and their dependence on the dislocation velocity, however few prior works have validated the accuracy of continuum theories. Here, in this work, molecular dynamics (MD) simulations of dislocation motion in face-centered cubic Ni containing interstitial H atoms were performed in order to test the accuracy of continuum theory predictions. Our results demonstrate that continuum theory provides an accurate estimate of the solute drag force for velocities below the critical velocity at which the solute drag force is maximized. Above the critical velocity, continuum theory systematically deviates from MD. Additional analysis reveals that this deviation results from the multi-valued and unstable nature of solute drag under load control, and also from the transition to random solute drag which occurs at high velocities.

Dislocations↗

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗

Quantitative insights into the dislocation source behavior of twin boundaries suggest a new dislocation source mechanism

Pop-in statistics from nanoindentation with spherical indenters are used to determine the stress required to activate dislocation sources in twin boundaries (TBs) in copper and its alloys. The TB source activation stress is smaller than that needed for bulk single crystals, irrespective of the indenter size, dislocation density and stacking fault energy. Because an array of pre-existing Frank partial dislocations is present at a TB, we propose that dislocation emission from the TB occurs by the Frank partials splitting into Shockley partials moving along the TB plane and perfect lattice dislocations, both of which are mobile. The proposed mechanism is supported by recent high resolution transmission electron microscopy images in deformed nanotwinned (NT) metals and may help to explain some of the superior properties of nanotwinned metals (e.g. high strength and good ductility), as well as the process of detwinning by the collective formation and motion of Shockley partial dislocations along TBs.

36 MATERIALS SCIENCE↗

Statistical mechanics of kinks on a gliding screw dislocation

The ability of a body-centered-cubic metal to deform plastically is limited by the thermally activated glide motion of screw dislocations, which are line defects with a mobility exhibiting complex dependence on temperature, stress, and dislocation segment length. We derive an analytical expression for the velocity of dislocation glide, based on a statistical mechanics argument, and identify an apparent phase transition marked by a crossover temperature above which the activation energy for glide effectively halves, changing from the formation energy of a double kink to that of a single kink. The analysis is in quantitative agreement with direct kinetic Monte Carlo simulations.

36 MATERIALS SCIENCE↗

Unfaulting mechanisms of interstitial Frank loops in fluorite-structured ThO 2

Unfaulting of Frank loops in irradiated fluorite-structured oxides is of significance for understanding microstructural evolution. However, the mechanisms have not yet been directly observed. To this end, we utilize molecular dynamics to reveal the atomistic details related to the unfaulting process of interstitial Frank loops in ThO 2 , which involve a single pair or multiple pairs of co-propagating Shockley partials to annihilate the stacking faults. Here, we find that the unfaulting is achieved via a synchronous shear of the partial pairs to remove the extrinsic stacking fault in the cation sublattice and the intrinsic stacking fault in the anion sublattice. The high mobility of oxygen at the dislocation core may reduce the activation barriers of dislocation nucleation and migration. These findings provide a fundamental understanding of the transformation of faulted loops in irradiated ThO 2 , and could be transferable to other fluorite-structured oxides.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Atomic-Scale Tracking of Topological Defect Motion and Incommensurate Charge Order Melting

Charge order pervades the phase diagrams of quantum materials where it competes with superconducting and magnetic phases, hosts electronic phase transitions and topological defects, and couples to the lattice generating intricate structural distortions. Incommensurate charge order is readily stabilized in manganese oxides, where it is associated with anomalous electronic and magnetic properties, but its nanoscale structural inhomogeneity complicates precise characterization and understanding of its relationship with competing phases. Leveraging atomic-resolution variable-temperature cryogenic scanning transmission electron microscopy, we characterize the thermal evolution of charge order as it transforms from its ground state in a model manganite system. We find that mobile networks of discommensurations and dislocations generate phase inhomogeneity and induce global incommensurability in an otherwise lattice-locked modulation. Driving the order to melt at high temperatures, the discommensuration density grows, and regions of order locally decouple from the lattice periodicity. Published by the American Physical Society 2025

Schnitzer, Noah (ORCID:0000000210421146)↗

Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb

Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.

36 MATERIALS SCIENCE↗