Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “DISLOCATIONS”

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 199 records · Page 11

A model of thermal creep and annealing in finite domains based on coupled dislocation climb and vacancy diffusion

Here, we develop a framework to investigate thermal creep and annealing in finite domains, where the climb motion of discrete dislocations is coupled to the diffusion of a continuum vacancy field. The model is first formulated in a continuum finite-deformation setting. All governing equations and boundary conditions are obtained from a unified irreversible thermodynamics principle. The resulting model couples a mechanical boundary value problem (BVP), a vacancy diffusion BVP, and the climb and glide motion of the discrete dislocation network within the crystal. The framework is then linearized for implementation in three-dimensional (3D) discrete dislocation dynamics (DDD) simulations for arbitrary anisotropic crystals. A solution scheme is developed based on the superposition principle, which is imposed weakly on the dislocation network to obtain a Galerkin solution for the nodal climb velocities. The framework includes diffusional (Nabarro–Herring) creep deformation as well as dislocation creep by climb-assisted-glide. The method is applied to simulate the annealing of vacancy loops in Al, with good agreement to experimental measurements by Silcox and Hirsch. We further consider the effects of annealing under stress, and of the proximity of the vacancy loops to loaded and free boundaries Simulations in polycrystalline materials are carried out to highlight the effects of the grain size on dislocation climb and vacancy loop annealing. The method is also applied to estimate the creep rate due to climb-assisted glide of jogged-screw dislocations in γ-TiAl, and results are compared to experiments by Viswanathan et al. Finally, we discuss the effects of uniaxial and hydrostatic stresses on the two diffusive deformation pathways of the material, namely Nabarro–Herring creep and dislocation climb.

36 MATERIALS SCIENCE↗

Significant disparity of non-basal dislocation activities in hot-rolled highly-textured Mg and Mg-3Al-1Zn alloy under tension

It is well-established that Mg-3Al-1Zn (AZ31) alloy exhibits much better tensile ductility than pure Mg. However, the underlying mechanisms for such difference still remain relatively unexplored at the dislocation level. In this work, we deformed hot-rolled highly-textured pure Mg and AZ31 samples with similar initial microstructures (i.e., grain size and texture) under tension along the rolling direction. Apparent shear banding was only observed in the pure Mg samples from the early stage of the deformation. Shear bands could act as easy paths for crack propagation, leading to low ductility in pure Mg. Apparent shear banding was not noted in deformed AZ31 samples. We then investigated the deformation mechanisms at the dislocation level using transmission electron microscopy. Systematic tilting experiments and statistical analyses of multiple grains at different strain levels revealed a significant disparity of non-basal dislocation activities between pure Mg and AZ31. For pure Mg, dislocations were activated since the early stage of plastic deformation. For AZ31, dislocations were mostly absent at all strain levels, even in the strain-to-failure samples. Non-basal dislocations, including prismatic and pyramidal dislocations, were observed. The promotion of the non-basal dislocation activities and the suppression of dislocations in AZ31 are expected to offer more sustainable hardening, which could elucidate the absence of apparent shear banding and much-improved ductility in AZ31 compared to pure Mg.

Zhao, Dexin↗

Screw dislocation impingement and slip transfer at fcc-bcc semicoherent interfaces

Understanding of dislocation-interface interactions at the atomistic level could provide useful insights into the mechanisms behind slip transfer, a key information for higher length-scale formulations. Here, a systematic study of screw dislocation impingement on both Fe/Al and Fe/Cu interfaces with Kurdjumov-Sachs orientation relationship has been carried out using atomistic simulations. For the “strong” Fe/Al interface, there is appreciable dislocation delocalization, in which the change in dislocation content provided by Nye tensor analysis. Interestingly, the delocalization is significant enough that large simulation cells must be employed in order to avoid boundary effects. When the impinging dislocations are far away from the boundaries, the large accumulation of shear strains nearby the impinging dislocations at the interface causes dislocation nucleation/emission at these places. While for the “weak” Fe/Cu interface, the dislocation delocalization at the interface is notably long-range, extending more than 20 nm away from where the initial dislocations entered the interface.

36 MATERIALS SCIENCE↗

Sequential slip transfer of mixed-character dislocations across Σ3 coherent twin boundary in FCC metals: a concurrent atomistic-continuum study

Sequential slip transfer across grain boundaries (GB) has an important role in size-dependent propagation of plastic deformation in polycrystalline metals. For example, the Hall–Petch effect, which states that a smaller average grain size results in a higher yield stress, can be rationalised in terms of dislocation pile-ups against GBs. In spite of extensive studies in modelling individual phases and grains using atomistic simulations, well-accepted criteria of slip transfer across GBs are still lacking, as well as models of predicting irreversible GB structure evolution. Slip transfer is inherently multiscale since both the atomic structure of the boundary and the long-range fields of the dislocation pile-up come into play. In this work, concurrent atomistic-continuum simulations are performed to study sequential slip transfer of a series of curved dislocations from a given pile-up on Σ3 coherent twin boundary (CTB) in Cu and Al, with dominant leading screw character at the site of interaction. A Frank-Read source is employed to nucleate dislocations continuously. It is found that subject to a shear stress of 1.2 GPa, screw dislocations transfer into the twinned grain in Cu, but glide on the twin boundary plane in Al. Moreover, four dislocation/CTB interaction modes are identified in Al, which are affected by (1) applied shear stress, (2) dislocation line length, and (3) dislocation line curvature. Our results elucidate the discrepancies between atomistic simulations and experimental observations of dislocation-GB reactions and highlight the importance of directly modeling sequential dislocation slip transfer reactions using fully 3D models.

36 MATERIALS SCIENCE↗

The core structure and recombination energy of a copper screw dislocation: a Peierls study

The recombination process of dislocations is central to cross-slip, and transmission through Σ3 grain boundaries among other fundamental plastic deformation processes. Despite its importance, a detailed mechanistic understanding remains lacking. In this paper, we apply a continuous dislocation model, inspired by Peierls and Nabarro, complete with an ab-initio computed -surface and continuous units of infinitesimal dislocation slip, towards computing the stress-dependent recombination path of both an isotropic and anisotropic Cu screw dislocation. Under no applied stress, our model reproduces the stacking fault width between Shockley partial dislocations as predicted by discrete linear elasticity. Upon application of a compressive Escaig stress, the two partial dislocations coalesce to a separation of ~|b|. Upon increased loading the edge components of each partial dislocation recede, leaving behind a spread Peierls screw dislocation, indicating the recombined state. We demonstrate that the critical stress required to achieve the recombined state is independent of the shear modulus. Rather the critical recombination stress depends on an energy difference between an unstable fault energy (γτ) and the intrinsic stacking fault energy (γτ-γisf). We report recombination energies of ΔW = 0.168 eV/Å and ΔW = 0.084 eV/Å, respectively, for the Cu screw dislocation within isotropic and anisotropic media. Finally, we develop an analytic model which provides insight into our simulation results which compare favourably with other (similar) models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First-principles study of atomic and electronic structures of 60° perfect and 30°/90° partial glide dislocations in CdTe

The atomic and electronic structures of 60° glide perfect and 30°/90° glide partial dislocations in CdTe are studied using combined semi-empirical and density functional theory calculations. The calculations predict that the dislocation cores tend to undergo significant reconstructions along the dislocation lines from the singly-periodic (SP) structures, yielding either doubly-periodic (DP) ordering by forming a dimer or quadruply-periodic (QP) ordering by alternating a dimer and a missing dimer. Charge modulation along the dislocation line, accompanied by the QP reconstruction for the Cd-/Te-core 60° perfect and 30° partials or the DP reconstruction for the Cd-core 90° partial, results in semiconducting character, as opposed to the metallic character of the SP dislocation cores. Dislocation-induced defect states for the 60° Cd-/Te-core are located relatively close to the band edges, whereas the defect states lie in the middle of the band gap for the 30° Cd-/Te-core partial dislocations. In addition to the intracore charge modulation within each QP core, the possibility of intercore charge transfer between two different dislocation cores when they are paired together in the same system is discussed. As a result, the analysis of the electronic structures reveals the potential role of the dislocations on charge transport in CdTe, particularly in terms of charge trapping and recombination.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Characterization of the tip field of a discrete dislocation pileup for the development of physically based micromechanics

It is shown, on the basis of calculations by Eshelby et al. (1951), Armstrong et al. (1966), and Chou and Li (1969), that a single parameter, such as the force on the leading dislocation (F), the crack extension force, or the stress intensity factor, is capable of characterizing uniquely the entire tip field of a discrete dislocation pileup, including the positions of mobile dislocations behind the locked leading dislocation at the tip. Conversely, the position of the i-th mobile dislocation X(i) is related to the value of F and is capable of characterizing the entire stress, strain, and displacement fields at the tip of a discrete dislocation pileup. If the interactions between dislocations are linear elastic, the measured positions of the mobile dislocations can be used to determine the value of F, which can then be used as a quantitative measure of the strength of a dislocation barrier resisting the propagation of a microslip or the nucleation of a microfracture.

Gao, Q.↗

Modeling of dislocation channel width evolution in irradiated metals

Defect-free dislocation channel formation has been reported to promote plastic instability during tensile testing via localized plastic flow, leading to a distinct loss of ductility and strain hardening in many low-temperature irradiated materials. In order to study the underlying mechanisms governing dislocation channel width and formation, the channel formation process is modeled via a simple stochastic dislocation-jog process dependent upon grain size, defect cluster density, and defect size. Dislocations traverse a field of defect clusters and jog stochastically upon defect interaction, forming channels of low defect-density. And based upon prior molecular dynamics (MD) simulations and in-situ experimental transmission electron microscopy (TEM) observations, each dislocation encounter with a dislocation loop or stacking fault tetrahedron (SFT) is assumed to cause complete absorption of the defect cluster, prompting the dislocation to jog up or down by a distance equal to half the defect cluster diameter. Channels are predicted to form rapidly and are comparable to reported TEM measurements for many materials. Predicted channel widths are found to be most strongly dependent on mean defect size and correlated well with a power law dependence on defect diameter and density, and distance from the dislocation source. Due to the dependence of modeled channel width on defect diameter and density, maximum channel width is predicted to slowly increase as accumulated dose increases. The relatively weak predicted dependence of channel formation width with distance, in accordance with a diffusion analogy, implies that after only a few microns from the source, most channels observed via TEM analyses may not appear to vary with distance because of limitations in the field-of-view to a few microns. Furthermore, examinations of the effect of the so-called “source-broadening” mechanism of channel formation showed that its effect is simply to add a minimum thickness to the channel without affecting channel dependence on the given parameters.

36 MATERIALS SCIENCE↗

Statistics of dislocation pinning at localized obstacles

Pinning of dislocations at nanosized obstacles like precipitates, voids, and bubbles is a crucial mechanism in the context of phenomena like hardening and creep. The interaction between such an obstacle and a dislocation is often studied at fundamental level by means of analytical tools, atomistic simulations, and finite element methods. Nevertheless, the information extracted from such studies cannot be utilized to its maximum extent on account of insufficient information about the underlying statistics of this process comprising a large number of dislocations and obstacles in a system. Here, we propose a new statistical approach, where the statistics of pinning of dislocations by idealized spherical obstacles is explored by taking into account the generalized size-distribution of the obstacles along with the dislocation density within a three-dimensional framework. Starting with a minimal set of material parameters, the framework employs the method of geometrical statistics with a few simple assumptions compatible with the real physical scenario. The application of this approach, in combination with the knowledge of fundamental dislocation-obstacle interactions, has successfully been demonstrated for dislocation pinning at nanovoids in neutron irradiated type 316-stainless steel in regard to the non-conservative motion of dislocations. An interesting phenomenon of transition from rare pinning to multiple pinning regimes with increasing irradiation temperature is revealed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Strongly correlated breeding of high-speed dislocations

Under very high stresses, dislocations can be accelerated to approach the speed of shear wave over a distance as short as 10 1 nm. Our atomistic simulations demonstrate that dislocations with such high speeds often react in counter-intuitive manners that are beyond textbook descriptions of conventional dislocation behavior. A high-speed dislocation can “rebound” when hitting a free surface rather than simply annihilate. When two high-speed dislocations collide, they can “penetrate” through each other. An individual dislocation can even spontaneously generate multiple dislocations via self-dissociation. These anomalous mechanisms lead to rapid proliferation of dislocations that are strongly correlated both spatially and temporally, and as such may play a role in high-stress and high-strain-rate plastic deformation; a potentially related case is nanoscale pristine single crystals, which often yield via a large strain burst at ultrahigh stresses.

36 MATERIALS SCIENCE↗

Edge dislocations in dicalcium silicates: Experimental observations and atomistic analysis

Understanding defects and influence of dislocations on dicalcium silicates (Ca{sub 2}SiO{sub 4}) is a challenge in cement science. We report a high-resolution transmission electron microscopy image of edge dislocations in Ca{sub 2}SiO{sub 4}, followed by developing a deep atomic understanding of the edge dislocation-mediated properties of five Ca{sub 2}SiO{sub 4} polymorphs. By decoding the interplay between core dislocation energies, core structures, and nucleation rate of reactivity, we find that γ-C2S and α-C2S polymorphs are the most favorable polymorphs for dislocations in Ca{sub 2}SiO{sub 4}, mainly due to their large pore channels which take away majority of the distortions imposed by edge dislocations. Furthermore, in the context of edge dislocation, while α-C2S represents the most active polymorph for reactivity and crystal growth, β-C2S represents the most brittle polymorph suitable for grinding. This work is the first report on the atomistic-scale analysis of edge dislocation-mediated properties of Ca{sub 2}SiO{sub 4} and may open up new opportunities for tuning fracture and reactivity processes of Ca{sub 2}SiO{sub 4} and other cement components.

36 MATERIALS SCIENCE↗

Field Dislocation Mechanics for heterogeneous elastic materials: A numerical spectral approach

Spectral methods using Fast Fourier Transform (FFT) algorithms have recently seen a surge in interest in the mechanics of materials community. The present work addresses the critical question of determining accurate local mechanical fields using FFT methods without artificial fluctuations arising from materials and defects induced discontinuities. Precisely, this work introduces a numerical approach based on intrinsic discrete Fourier transforms for the simultaneous treatment of material discontinuities arising from the presence of dislocations and from elastic stiffness heterogeneities. To this end, the elasto-static equations of the field dislocation mechanics theory for periodic heterogeneous materials are numerically solved with FFT in the case of dislocations in proximity of inclusions of varying stiffness. An optimal intrinsic discrete Fourier transform method is sought based on two distinct schemes. A centered finite difference scheme for differential rules are used for numerically solving the Poisson-type equation in the Fourier space, while centered finite differences on a rotated grid is chosen for the computation of the modified Fourier–Green’s operator associated with the Lippmann–Schwinger-type equation. By comparing different methods with analytical solutions for an edge dislocation in a composite material, it is found that the present spectral method is accurate, devoid of any numerical oscillation, and efficient even for an infinite phase elastic contrast like a hole embedded in a matrix containing a dislocation. The present FFT method is then used to simulate physical cases such as the elastic fields of dislocation dipoles located near the matrix/inclusion interface in a 2D composite material and the ones due to dislocation loop distributions surrounding cubic inclusions in 3D composite material. In these configurations, the spectral method allows investigating accurately the elastic interactions and image stresses due to dislocation fields in the presence of elastic inhomogeneities.

36 MATERIALS SCIENCE↗

Effect of a micro-scale dislocation pileup on the atomic-scale multi-variant phase transformation and twinning

In this paper, we perform concurrent atomistic–continuum (CAC) simulations to assess the contribution of the internal stress induced by the microscale dislocation pileup at an atomically structured interface 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) the pileup-induced internal stress concentration profile follows neither the classical Eshelby model nor the super-dislocation model alone, but a combination of them; (c) 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. The two variants of the hexagonal phases form a twin with respect to each other; (d) 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 (e) the microscale dislocation pileup-induced internal shear stress and the structural change at the atomic-scale interface 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. In conclusion, the gained knowledge will advance our understanding of how the multi-phase material behaves in many complex physical processes, such as the synthesis of multi-phase high-entropy alloys or superhard ceramics under high-pressure torsion, deep mantle earthquakes in geophysics, and so on, which all involve dislocation slip, PTs, twinning, and their interactions across from the atomistic to the microscale and beyond.

36 MATERIALS SCIENCE↗

Coupling continuum dislocation transport with crystal plasticity for application to shock loading conditions

Here, we have developed a multi-physics modeling approach that couples continuum dislocation transport, nonlinear thermoelasticity, crystal plasticity, and consistent internal stress and deformation fields to simulate the single-crystal response of materials under extreme dynamic conditions. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. Nonlinear thermoelasticity provides a thermodynamically consistent equation of state to relate stress (including pressure), temperature, energy densities, and dissipation. Crystal plasticity is coupled to dislocation transport via Orowan's expression where the constitutive description makes use of recent advances in dislocation velocity theories applicable under extreme loading conditions. The configuration of geometrically necessary dislocation density gives rise to an internal stress field that can either inhibit or accentuate the flow of dislocations. An internal strain field associated with the internal stress field contributes to the kinematic decomposition of the overall deformation. The paper describes each theoretical component of the framework, key aspects of the constitutive theory, and some details of a one-dimensional implementation. Results from single-crystal copper plate impact simulations are discussed in order to highlight the role of dislocation transport and pile-up in shock loading regimes. The main conclusions of the paper reinforce the utility of the modeling approach to shock problems.

36 MATERIALS SCIENCE↗

Comparative modeling of the disregistry and Peierls stress for dissociated edge and screw dislocations in Al

Many elementary deformation processes in metals involve the motion of dislocations. The planes of glide and specific processes dislocations prefer depend heavily on their atomic core structures. Atomistic simulations are desirable for dislocation modeling but their application to even sub-micron scale problems is in general computationally costly. Accordingly, continuum-based approaches, such as the phase-field microelasticity, phase-field dislocation dynamics (PFDD), generalized Peierls–Nabarro (GPN) models, and the concurrent atomistic–continuum (CAC) method, have attracted increasing attention in the field of dislocation modeling because they well represent both short-range cores interactions and long-range stress fields of dislocations. To better understand their similarities and differences, it is useful to compare these methods in the context of benchmark simulations and predictions. In this paper, we apply the CAC method and different PFDD variants – one of them is equivalent to a GPN model – to simulate an extended (i.e., dissociated) dislocation in Al with initially pure edge or pure screw character in terms of the disregistry. CAC and discrete forms of PFDD are also employed to calculate the Peierls stress. By conducting comprehensive convergence studies, we quantify the dependence of these measures on time/grid resolution and simulation cell size. Several important but often overlooked differences between PFDD/GPN variants are clarified. In conclusion, our work sheds light on the advantages and limitations of each method, as well as the path towards enabling them to effectively model complex dislocation processes at larger length scales.

36 MATERIALS SCIENCE↗

Dislocation drag and its influence on elastic precursor decay

Plastic deformation is mediated by the creation and movement of dislocations, and at high stress the latter is dominated by dislocation drag from phonon wind. Here, by simulating a 1-D shock impact problem we analyze the importance of accurately modeling dislocation drag and dislocation density evolution in the high stress regime. Dislocation drag is modeled according to a first-principles derivation as a function of stress and dislocation character, and its temperature and density dependence are approximated to the extent currently known. Much less is known about dislocation density evolution, leading to far greater uncertainty in these model parameters. In studying anisotropic fcc metals with character dependent dislocations, the present work generalizes similar earlier studies by other authors.

36 MATERIALS SCIENCE↗

A dislocation density-based continuum model of the anisotropic shock response of single crystal α-cyclotrimethylene trinitramine

Here, we have developed a model for the finite deformation thermomechanical response of α-cyclotrimethylene trinitramine (RDX). Our model accounts for nonlinear thermoelastic lattice deformation through a free energy-based equation of state developed by Cawkwell et al. (2016) in combination with temperature and pressure dependent elastic constants, as well as dislocation-mediated plastic slip on a set of slip systems motivated by experimental observation. The kinetics of crystal plasticity are modeled using the Orowan equation relating slip rate to dislocation density and the dislocation velocity developed by Austin and McDowell (2011), which naturally accounts for transition from thermally activated to dislocation drag limited regimes. Evolution of dislocation density is specified in terms of local ordinary differential equations reflecting dislocation–dislocation interactions. This paper presents details of the theory and parameterization of the model, followed by discussion of simulations of flyer plate impact experiments. Impact conditions explored within this combined simulation and experimental effort span shock pressures ranging from 1 to 3 GPa for four crystallographic orientations and multiple specimen thicknesses. Simulation results generated using this model are shown to be in strong agreement with velocimetry measurements from the corresponding plate impact experiments. Finally, simulation results are used to motivate conclusions about the nature of dislocation-mediated plasticity in RDX.

36 MATERIALS SCIENCE↗

A dislocation-based crystal plasticity framework for dynamic ductile failure of single crystals

We developed a framework for dislocation-based viscoplasticity and dynamic ductile failure to model high strain rate deformation and damage in single crystals. The rate-dependence of the crystal plasticity formulation is based on the physics of relativistic dislocation kinetics suited for extremely high strain rates. The damage evolution is based on the dynamics of void growth, which are governed by both micro-inertia as well as dislocation kinetics and dislocation substructure evolution. Furthermore, an averaging scheme is proposed in order to approximate the evolution of the dislocation substructure in both the macroscale as well as its spatial distribution at the microscale. In addition, a concept of a single equivalent dislocation density that effectively captures the collective influence of dislocation density on all active slip systems is proposed here. Together, these concepts and approximations enable the use of semi-analytic solutions for void growth dynamics developed in [J. Wilkerson and K. Ramesh. A dynamic void growth model governed by dislocation kinetics. J. Mech. Phys. Solids, 70:262–280, 2014.], which greatly reduce the computational overhead that would otherwise be required. The resulting homogenized framework has been implemented into a commercially available finite element package, and a validation study against a suite of direct numerical simulations was carried out.

36 MATERIALS SCIENCE↗