Search NASA⌕ Search

Engineering topics

El-Azab, Anter

Publications and source records attributed to El-Azab, Anter.

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗

Evolution of dislocation loops in irradiated α-Uranium: An atomistically-informed cluster dynamics investigation

An atomistically-informed mean field cluster dynamics model has been developed to investigate the nucleation and growth of prismatic loops in irradiated a-Uranium. TEM analysis of neutron irradiated a-Uranium shows the evolution of self-interstitial atom and vacancy loops on (010) and (100) crystallographic planes, respectively, resulting in an anisotropic lattice swelling of its face-centered orthorhombic crystal. To provide model parameters, the crystallography of loops and the binding energy of point defects to these loops were studied using an angular dependent EAM potential and classical molecular dynamics (MD) simulations. Furthermore, using the bond-boost hyperdynamics method, the anisotropic diffusion of interstitials and vacancies in a-Uranium was investigated. Here, the mechanisms of point defect diffusion and the associated migration energies were reported and compared with previous DFT studies. The energetics and kinetic quantities mentioned above were used in the cluster dynamics model to predict the cluster density as a function of dose rate, dose and temperature and the results were compared to the reported neutron irradiation experiments. The model predictions reveal an accumulation of small sized vacancy loops along with a population of large and growing self-interstitial loops, which closely corresponds to the TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evolution of dislocation loops and voids in post-irradiation annealed ThO 2 : A combined in-situ TEM and cluster dynamics investigation

The effect of isochronal annealing on the evolution of dislocation loop and void population in proton irradiated ThO 2 has been investigated. Post-irradiation annealing in other actinide oxides like UO 2 shows significant loop coarsening. ThO 2 samples were irradiated with 2 MeV protons up to a dose of 0.1 dpa at 600°C. Post-irradiation isochronal annealing was performed at 600, 800, 1000 and 1100°C for 1 h at each temperature using in-situ TEM. Only faulted 1/3<111> type dislocation loops were observed, and their sizes and distribution were characterized. Further, the population of self-interstitial atom (SIA) dislocation loops did not show any significant growth and coarsening. Additionally, nanometric voids were observed at annealing temperatures of 1000 and 1100°C. Using cluster dynamics (CD), we have studied the nucleation and growth of point defects and defect clusters, i.e., SIA prismatic dislocation loops and nanometric and sub-nanometric voids in proton irradiated ThO 2 . The CD model was further utilized to predict the growth and coarsening of loops and voids during isochronal annealing at the experimental and higher temperatures. The model did not predict significant SIA loop growth which closely corresponds to the TEM observations. CD predicted SIA loop coarsening is insignificant even at high annealing temperature of 1500°C because the model only considers the growth of defect clusters by absorption of like point defects, i.e., SIA loops absorb interstitials and voids absorb vacancies, and cannot account for their migration and coalescence due to elastic interaction. The CD model also predicts the evolution of nanometric voids having mean size within the error bounds of TEM observations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Plasticity in irradiated FeCrAl nanopillars investigated using discrete dislocation dynamics

Here, in this paper, we investigate plasticity in irradiated FeCrAl nanopillars using discrete dislocation dynamics simulations (DDD), with comparisons to transmission electron microscopic (TEM) in situ tensile tests of ion- and neutron-irradiated commercial C35M FeCrAl alloy. The effects of irradiation-induced defects, such as a/2<111> and a<100> type loops and composition fluctuations representative of phase separation in irradiated FeCrAl alloys, are investigated separately as well as superposed together in simulations. We explore the effects of defects on the stress-strain behavior, specifically yield strength and hardening response, of FeCrAl nanopillars. Our simulations confirm the widely accepted fact that irradiated alloys exhibit a stress-strain response with higher yield strength as compared to unirradiated alloys. However, our DDD calculations reveal an atypical superposition of the hardening contributions due to composition inhomogeneity and irradiation loops wherein composition inhomogeneity annihilates the hardening due to irradiation loops at small scales. As a result, we observe that the yield strength in irradiated alloys, after taking into consideration the effects of both composition inhomogeneity and irradiation loops, is smaller than the yield strength of the alloys with only irradiation loops and is approximately same for the alloy with composition inhomogeneity alone. This is referred to as “destructive interference” between the hardening contributions due to composition fluctuations and irradiation loops in the paper. We also identify this destructive interference in the superposition in our parallel TEM in situ tensile tests on unirradiated, ion-irradiated, and neutron-irradiated C35M FeCrAl alloy. This destructive interference in the hardening contributions contrasts with the dispersed barrier hardening (DBH) models widely utilized by the experimental community to model the hardening contributions due to different irradiation induced defects. The effects of the loading orientations on the yield strength and hardening are investigated and the mechanisms for the hardening in irradiated FeCrAl alloys are also reported.

36 MATERIALS SCIENCE↗

A data driven approach for cross-slip modelling in continuum dislocation dynamics

Cross-slip is a thermally activated process by which screw dislocation changes its glide plane to another slip plane sharing the same Burgers vector. The rate at which this process happens is determined by a Boltzmann type expression that is a function of the screw segment length and the stress acting on the dislocation. In continuum dislocation dynamics (CDD), the information regarding the length of the screw dislocation segment and local stress state on dislocations are lost due to the coarse-grained representation of the density. Here, in this work, a data driven approach to characterize the lost information by analyzing the discrete dislocation configurations is proposed to enable cross-slip modeling in the CDD framework in terms of the coarse-grained dislocation density and stress fields. The analysis showed that the screw segment length follows an exponential distribution, and the stress fluctuations, defined as the difference between the stress on the dislocations and the mean field stress in CDD, follows a Lorentzian distribution. A novel approach for cross slip implementation in CDD employing the screw segment length and stress fluctuation statistics was proposed and rigorously tested by comparing the CDD cross-slip rates with discrete dislocation dynamics (DDD) rates. This approach has been applied in conjunction with three cross-slip models used in DDD simulations differing mainly in the functional form of cross slip activation energy. It was found that different cross-slip activation energy formulations yielded different cross-slip rates, yet the effect on mechanical stress-strain response and dislocation density evolution was minimal for the [001] type loading.

42 ENGINEERING↗

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↗