skip to main content


This content will become publicly available on August 1, 2024

Title: Effect of a Long-Range Dislocation Pileup on the Atomic-Scale Hydrogen Diffusion near a Grain Boundary in Plastically Deformed bcc Iron

In this paper, we present concurrent atomistic-continuum (CAC) simulations of the hydrogen (H) diffusion along a grain boundary (GB), nearby which a large population of dislocations are piled up, in a plastically deformed bi-crystalline bcc iron sample. With the microscale dislocation slip and the atomic structure evolution at the GB being simultaneously retained, our main findings are: (i) the accumulation of tens of dislocations near the H-charged GB can induce a local internal stress as high as 3 GPa; (ii) the more dislocations piled up at the GB, the slower the H diffusion ahead of the slip–GB intersection; and (iii) H atoms diffuse fast behind the pileup tip, get trapped within the GB, and diffuse slowly ahead of the pileup tip. The CAC simulation-predicted local H diffusivity, Dpileup−tip, and local stresses, σ, are correlated with each other. We then consolidate such correlations into a mechanics model by considering the dislocation pileup as an Eshelby inclusion. These findings will provide researchers with opportunities to: (a) characterize the interplay between plasticity, H diffusion, and crack initiation underlying H-induced cracking (HIC); (b) develop mechanism-based constitutive rules to be used in diffusion–plasticity coupling models for understanding the interplay between mechanical and mass transport in materials at the continuum level; and (c) connect the atomistic deformation physics of polycrystalline materials with their performance in aqueous environments, which is currently difficult to achieve in experiments.

 
more » « less
Award ID(s):
2322675
NSF-PAR ID:
10486510
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Crystals
Volume:
13
Issue:
8
ISSN:
2073-4352
Page Range / eLocation ID:
1270
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. In this paper, in situ high-resolution electron backscattered diffraction (EBSD) is combined with concurrent atomistic-continuum (CAC) simulations to study the interactions between dislocation-mediated slip and grain boundaries (GBs) in Ni. It is found that the local stress associated with slip-GB intersections first increases upon the pileup of dislocations, then remains high even after the nucleation of dislocations in the neighboring grain, only relaxing after the nucleated dislocations propagate away from the GB due to more incoming dislocations participating in the pileup. The local stress relaxation is accompanied by an atomic-scale GB structure reconfiguration, which affects not only the subsequent dislocation transmission, but also the configuration of those dislocations away from the GB. These findings demonstrate the importance of incorporating local stress history at higher length scale models, such as crystal plasticity finite element. 
    more » « less
  2. In this paper, we perform concurrent atomistic-continuum (CAC) simulations to (i) characterize the internal stress induced by the microscale dislocation pileup at an atomically structured interface; (ii) decompose this stress into two parts, one of which is from the dislocations behind the pileup tip according to the Eshelby model and the other is from the dislocations at the pileup tip according to a super-dislocation model; and (iii) assess how such internal stresses contribute to the atomic-scale phase transformations (PTs), reverse PTs, and twinning. The main novelty of this work is to unify the atomistic description of the interface and the coarse-grained (CG) description of the lagging dislocations away from the interface within one single framework. Our major findings are: (a) the interface dynamically responds to a pileup by forming steps/ledges, the height of which is proportional to the number of dislocations arriving at the interface; (b) when the pre-sheared sample is compressed, a direct square-to-hexagonal PT occurs ahead of the pileup tip and eventually grows into a wedge shape; (c) upon a further increase of the loading, part of the newly formed hexagonal phase transforms back to the square phase. The square product phase resulting from this reverse PT forms a twin with respect to the initial square phase. All phase boundaries (PBs) and twin boundaries (TBs) are stationary and correspond to zero thermodynamic Eshelby driving forces; and (d) the stress intensity induced by a pileup consisting of 16 dislocations reduces the stress required for initiating a PT by a factor of 5.5, comparing with that in the sample containing no dislocations. This work is the first characterization of the behavior of PTs/twinning resulting from the reaction between a microscale dislocation slip and an atomically structured interface. 
    more » « less
  3. Recent developments in generalized continuum modeling methods ranging from coarse-grained atomistics to micromorphic theory offer potential to make more intimate physical contact with dislocation field problems framed at length scales on the order of microns. We explore a range of discrete dynamical and continuum mechanics approaches to crystal plasticity that are relevant to modeling behavior of populations of dislocations. Predictive atomistic and coarse-grained atomistic models are limited in terms of length and time scales that can be accessed; examples of the latter are discussed in terms of interactions of multiple dislocations in heterogeneous systems. Generalized continuum models alleviate restrictions to a significant extent in modeling larger scales of dislocation configurations and reactions, and are useful to consider effects of dislocation configuration on strength at characteristic length scales of sub-micron and above; these models require a combination of bottomup models and top-down experimental information to inform parameters and model form. The concurrent atomistic-continuum (CAC) method is extended to model complex multicomponent alloy systems using an average atom approach. Examples of CAC are presented, along with potential to assist in informing parameters of a recently developed micropolar crystal plasticity model based on a set of sub-micron dislocation field problems. Prospects for further developments are discussed. 
    more » « less
  4. The paper presents a multiscale study of the kinetic processes of the heteroepitaxial growth of the PbSe/PbTe (111) and PbTe/PbSe(001) systems, using the Concurrent Atomistic-Continuum (CAC) method as the simulation tool. The CAC simulations have reproduced the Stranski–Krastanov growth mode and the layer-by-layer growth mode of the two systems, respectively; the pyramid-shaped island morphology of the PbSe epilayer on PbTe (111), the square-like misfit dislocation networks within the PbTe/PbSe(001) interface, and the critical thickness for the PbTe/PbSe(001) system at which coherent interfaces transit to semi-coherent interfaces with the formation of misfit dislocations, all in good agreement with experimental observations. Four types of misfit dislocations are found to form during the growth of the two PbTe/PbSe heterosystems, and hexagonal-like misfit dislocation networks are observed within the PbSe/PbTe(111) interfaces. The growth processes, including the formation of misfit dislocations, have been visualized. Dislocation half-loops have been observed to nucleate from the epilayer surfaces. These half-loops extend towards the interface by climb or glide motions, interact with other half-loops, and form misfit dislocation networks at the interfaces and threading dislocations extending from interfaces to epilayer surfaces. The dominant types of misfit dislocations in both systems are found to be those with Burgers vectors parallel to the interfaces, whereas the misfit dislocations with Burgers vectors inclined to the interface have a low likelihood of generation and tend to annihilate. The size of the substrate is demonstrated to have a significant effect on the formation, evolution, and distribution of dislocations on the growth of PbSe on PbTe(111). 
    more » « less
  5. In this work, we investigate misfit dislocations in PbTe/PbSe heteroepitaxial systems using the concurrent atomistic–continuum (CAC) method. A potential model containing the long-range Coulombic interaction and short-range Buckingham potential is developed for the system. By considering the minimum potential energy of relaxed interface structures for various initial conditions and PbTe layer thicknesses, the equilibrium structure of misfit dislocations and the dislocation spacings in PbTe/PbSe(001) heteroepitaxial thin films are obtained as a function of the PbTe layer thicknesses grown on a PbSe substrate. The critical layer thickness above which misfit dislocations inevitably form, the structure of the misfit dislocations at the interfaces, and the dependence of average dislocation spacing on PbTe layer thickness are obtained and discussed. The simulation results provide an explanation for the narrowing of the spread of the distribution of misfit dislocation spacing as layer thickness increases in PbTe/PbSe(001) heteroepitaxy. 
    more » « less