Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dislocations in crystals”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Oscillating-Crucible Technique for Silicon Growth

Technique yields better mixing of impurities and superior qualiity crystals. Accellerated motion stirs melt which reduces temperature gradients and decreases boundary layer for diffusion of impurities near growing surface. Results better mixing of impurities into melt, decrease in tendency for dendritic growth or cellular growth and crystals with low dislocation density. Applied with success to solution growth and Czochralski growth, resulting in large crystals of superior quality.

Daud, T.↗
DOE OSTI2021

Inhomogeneous time-reversal symmetry breaking in Sr 2 RuO 4

We show that the observed time-reversal symmetry breaking (TRSB) of the superconducting state in Sr 2 RuO 4 can be understood as originating from inhomogeneous strain fields near edge dislocations of the crystal. Specifically, we argue that, without strain inhomogeneities, Sr 2 RuO 4 is a single-component, time-reversal symmetric superconductor, likely with d x 2 -y 2 symmetry. However, due to the strong strain inhomogeneities generated by dislocations, a slowly decaying subleading pairing state contributes to the condensate in significant portions of the sample. As it phase winds around the dislocation, time-reversal symmetry is locally broken. Global phase locking and TRSB occur at a sharp Ising transition that is not accompanied by a change of the single-particle gap and yields a very small heat capacity anomaly. Our model thus explains the puzzling absence of a measurable heat capacity anomaly at the TRSB transition in strained samples and the dilute nature of the time-reversal symmetry broken state probed by muon spin rotation experiments. We propose that plastic deformations of the material may be used to manipulate the onset of broken time-reversal symmetry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Crack tip microplasticity mediated by microstructure gradients

Abstract Traditional fracture theories infer damage at cracks (local field) by surveying loading conditions away from cracks (far field). This approach has been successful in predicting ductile fracture, but it normally assumes isotropic and homogeneous materials. However, myriads of manufacturing procedures induce heterogeneous microstructural gradients that can affect the accuracy of traditional fracture models. This work presents a microstructure‐sensitive finite element approach to explore the shielding effects of grain size and crystallographic orientation gradients on crack tip microplasticity and blunting. A dislocation density‐based crystal plasticity model conveys texture evolution, grain size effects, and directional hardening by computing the constraint from dislocation structures. The results demonstrate that the microstructure can act as a buffer between the local and far fields that affects the crack tip microplasticity variability. For nominal opening loading, grain size and texture affect the local ductility and induce a non‐negligible multiaxial plastic deformation. Furthermore, driving forces based on measuring displacements away from the crack tip are less affected by the microstructure, which suggests that traditional experimental methods smear out important crack tip variability.

Castelluccio, Gustavo M.↗
NASA NTRS2014

Data Processing and Experimental Design for Micrometeorite Impacts in Small Bodies

Comets and asteroids have been altered from their original "pristine" state by impacts occurring throughout their 4.5 billion year lives: [1]. Proof of shock deformation has been detected in the crystal structure of several Stardust samples from Comet Wild 2 [2, 3]. Analyses indicated that the planar dislocations in the crystal structure of the minerals had been imparted by impacts sustained during their lives, and not due to the aerogel capture process. Distortions to crystal structure also affect the ideal absorption spectra in the infrared, and [4], thus providing indirect evidence of its impact history and a means of remotely investigating the impact history of small bodies through comparing laboratory spectra with spectra observed by telescopes or spacecraft. -The effects of impacts propagating shock waves through minerals were investigated through laboratory impact experiments. Utilizing NASA Johnson Space Center's Experimental Impact Laboratory, projectiles were fired from the vertical gun at velocities ranging from 2.0 to 2.8 km/sec, projected impact velocities between Kuiper Belt Objects. Two types of projectiles were used, including spherical alumina ceramic, whose density mimics that of rock, and cylinders made from the same material that they impacted. The target materials chosen for testing included: OLIVINES forsterite (Mg2SiO4) and fayalite, Fe2SiO4); PYROXENES enstatite (Mg2Si2O6) and diopside (MgCaSi2O6); and CARBONATES magnesite (MgCO3) and siderite (FeCO3). Targets were impacted at either 25 C or cooled to -20 C to examine the effects of temperature, if any, on lattice distortions during the shock propagation. As comets and asteroids can undergo a wide range of temperatures in their orbital lifetimes, the effect of temperature on the equation of state of minerals being shocked needs to be examined for interpreting the results of these experiments. The porosity of the target mineral is varied by either grinding it into a powder/granular texture or as whole mineral rocks to investigate the differences in shock propagation when voids are present. By varying velocity, ambient temperature, and porosity, we can investigate different variables affecting impacts in the solar system. -Data indicates that there is a non-linear relationship between peak shock pressure and the variation in infrared spectral absorbances by the distorted crystal structure. The maximum variability occurs around 37 GPa in enstatite and forsterite. The particle size distribution of the impacted material similarly changes with velocity/peak shock pressure. -The experiments described above are designed to measure the near- to mid-IR effects from these changes to the mineral structure. See Lederer et al., this meeting for additional experimental results.

Jensen, E.↗
DOE OSTI2024

Relative molecular orientation can impact the onset of plasticity in molecular crystals

Abstract Creating or moving dislocations is the first step to dissipating mechanical energy via plastic deformation under contact loading. In molecular crystals there is both a lattice that defines crystal orientation and a relative orientation of the basis of the molecules. We define a normalization parameter which relates strain at yield, the hardness of the bulk crystal, and a distance parameter analogous to a Burgers vector that nominally predicts the relative ease of initiating plasticity in this broad class of materials. Analyzing the yield behavior of 10 different molecular crystals of varying space groups shows the inter-molecular orientation predicts the experimentally observed applied stress needed to nucleate dislocations. When molecules are oriented ‘parallel’ relative to one another the normalized maximum shear stress at the onset of plasticity is on the order of 3–5 times lower than when molecules within the crystal are ‘anti-parallel’, and molecules with a more equiaxed shape fall in between these bounds. This provides an initial indication of a structural feature which predicts the relative ease of initiating plasticity during contact loading in molecular crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

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↗
DOE OSTI2025

Coupled cluster and dislocation dynamics modeling of microstructure evolution in irradiated materials

We develop here a coupled cluster and dislocation dynamics framework to study the microstructure evolution of irradiated materials. The framework not only accounts for the three dimensional diffusion of radiation-generated clusters, but also their interaction with dislocation networks and the resultant climb motion of discrete dislocations within finite crystals. The framework is solved with a superposition solution scheme, and is applied to investigate the evolution of the irradiation-induced dislocation loops in zirconium (Zr), considering the effects of various bias factors including the diffusion anisotropy difference (DAD) of interstitials and interstitial clusters, the dislocation bias of defects to discrete dislocation segments, and the production bias of defects from the radiation cascade. We find that the DAD is the most critical factor influencing the kinetics of the loop evolution in Zr, while the recombination/interaction of mobile defects can induce a strong spatial dependence of the loop evolution together with the DAD. Here, the method is also adopted to study the evolution of interstitial $\langle$a$\rangle$ and vacancy $\langle$c$\rangle$ dislocation loop ensembles consistent with the microstructure observed during irradiation-induced growth of Zr. Our findings not only reveal the spatial dependence of the size and ellipticity of the dislocation loops, but also suggest a limit on the anisotropy factor of interstitials to reproduce the co-growth of $\langle$a$\rangle$ and $\langle$c$\rangle$ loops in zirconium, in good agreement with experimental observations and other simulation results.

Bias factors↗
DOE OSTI2022

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↗
DOE OSTI2020

Operando Bragg Coherent Diffraction Imaging of LiNi 0.8 Mn 0.1 Co 0.1 O 2 Primary Particles within Commercially Printed NMC811 Electrode Sheets

Due to complex degradation mechanisms, disparities between the theoretical and practical capacities of lithium-ion battery cathode materials persist. Specifically, Ni-rich chemistries such as LiNi 0.8 Mn 0.1 Co 0.1 O 2 (or NMC811) are one of the most promising choices for automotive applications; however, they continue to suffer severe degradation during operation that is poorly understood, thus challenging to mitigate. Here in this paper we use operando Bragg coherent diffraction imaging for 4D analysis of these mechanisms by inspecting the individual crystals within primary particles at various states of charge (SoC). Although some crystals were relatively homogeneous, we consistently observed non-uniform distributions of inter- and intracrystal strain at all measured SoC. Pristine structures may already possess heterogeneities capable of triggering crystal splitting and subsequently particle cracking. During low-voltage charging (2.7–3.5 V), crystal splitting may still occur even during minimal bulk deintercalation activity; and during discharging, rotational effects within parallel domains appear to be the precursor for the nucleation of screw dislocations at the crystal core. Ultimately, this discovery of the central role of crystal grain splitting in the charge/discharge dynamics may have ramifications across length scales that affect macroscopic performance loss during real-world battery operation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Bridgman-type apparatus for the study of growth-property relationships - Arsenic vapor pressure-GaAs property relationship

A precision Bridgman-type apparatus is described which was designed and constructed for the investigation of relationships between crystal growth parameters and the properties of GaAs crystals. Several key features of the system are highlighted, such as the use of a heat pipe for precise arsenic vapor pressure control and seeding without the presence of a viewing window. Pertinent growth parameters, such as arsenic source temperature, thermal gradients in the growing crystal and in the melt, and the macroscopic growth velocity can be independently controlled. During operation, thermal stability better than + or - 0.02 C is realized; thermal gradients can be varied up to 30 C/cm in the crystal region, and up to 20 C/cm in the melt region; the macroscopic growth velocity can be varied from 50 microns/hr to 6.0 cm/hr. It was found that the density of dislocations depends critically on As partial pressure; and essentially dislocation-free, undoped, crystals were grown under As pressure precisely controlled by an As source maintained at 617 C. The free carrier concentration varied with As pressure variations. This variation in free carrier concentration was found to be associated with variations in the compensation ratio rather than with standard segregation phenomena.

Parsey, J. M.↗
DOE OSTI2024

Exascale Dislocation Simulator

ExaDiS, or Exascale Dislocation Simulator, is a set of software modules written to enable numerical simulations of large groups of moving and interacting dislocations, line defects in crystals responsible for crystal plasticity. By tracking time evolution of sufficiently large dislocation ensembles, ExaDiS predicts plasticity response and plastic strength of crystalline materials. ExaDiS is built around a library of core functions for Discrete Dislocation Dynamics method specifically written to perform efficient computations on GPUs.

Bertin, Nicolas↗
NASA NTRS1998

New AFM Techniques for Investigating Molecular Growth Mechanisms of Protein Crystals

Atomic Force Microscopy (AFM) has emerged as a powerful technique for investigating protein crystal growth. Earlier AFM studies were among the first to demonstrate that these crystals grew by dislocation and 2D nucleation growth mechanisms [1]. These investigations were restricted to the micron range where only surface features, such as dislocation hillocks and 2D islands are visible. Most AFM instruments can scan at higher resolutions and have the potential to resolve individual protein molecules at nanometer ranges. Such scans are essential for determining the molecular packing arrangements on crystal faces and for probing the growth process at the molecular level. However, at this resolution the AFM tip influences the image produced, with the resulting image being a convolution of the tip shape and the surface morphology [2]. In most studies this problem is resolved by deconvoluting the image to obtain the true surface morphology. Although deconvolution routines work reasonably well for simple one- dimensional shapes, for complex surfaces this approach does not produce accurate results. In this study we devised a new approach which takes advantage of the precise molecular order of crystal surfaces, combined with the knowledge of individual molecular shapes from the crystallographic data of the protein and the AFM tip shape. This information is used to construct expected theoretical AFM images by convoluting the tip shape with the constructed crystal surface shape for a given surface packing arrangement. By comparing the images from actual AFM scans with the constructed ones for different possible surface packing arrangements, the correct packing arrangement can be conclusively determined. This approach was used in this study to determine the correct one from two possible packing arrangements on (I 10) faces of tetragonal lysozyme crystals. Another novel AFM technique was also devised to measure the dimension of individual growth units of the crystal faces. Measuring these units was not attempted before and most studies have assumed that the growth unit consisted of individual protein molecules. The linescan mode of AFM instruments allows the crystal surface to be scanned along a single line. By scanning across a growth step an image showing the motion of the step is obtained. Normally such an image shows a straight line for continuous and constant step velocity. In this study by increasing the scan rate and by decreasing the step velocity (by decreasing the supersaturation), we were able to capture images of individual growth events, shown by jump discontinuities in the step line. By suitable integration of the image the growth unit dimension in the scanned direction can be obtained. Since multiple units can be involved in the growth process it is necessary to collect a statistically relevant sample before drawing conclusions about the growth mechanism. This technique was successfully employed to obtain the dimensions of growth units for the (110) face, showing that they consisted of various aggregates corresponding to the 43 helices in the crystal structure.

Li, Huayu↗
DOE OSTI2023

Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Microstructural evolution on crystallizing the glassy phase in a 6 weight percent Y2O3-Si3N4 ceramic

X-ray diffraction and analytical electron microscopy have been used to study the crystallization of the grain-boundary glass in a 6 wt pct Y2O3-Si3N4 ceramic. Upon crystallization, high densities of dislocations formed in the Si3N4 grains and remained after 5 h at temperature. However, prolonged holds at the crystallization temperature effectively annealed out the dislocations. Other features present in the microstructure are characterized.

Lee, W. E.↗
DOE OSTI2022

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗