Atomic-Scale Investigation of Dislocations and Interfacial Line Defects in Layered Chalcogenide Materials.
Abstract not provided.
SEARCH · Engineering Papers
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.
Abstract not provided.
Journal article based on previous SDRD work is to be submitted to Springer Nature
Abstract not provided.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The research under this project focused on a theoretical and computational modeling of dislocation dynamics of mesoscale deformation of metal single crystals. Specifically, the work aimed to implement a continuum statistical theory of dislocations to understand strain hardening and cell structure formation under monotonic loading. These aspects of crystal deformation are manifestations of the evolution of the underlying dislocation system under mechanical loading. The project had three research tasks: 1) Investigating the statistical characteristics of dislocation systems in deformed crystals. 2) Formulating kinetic equations of dislocations and coupling these kinetics equations and crystal mechanics. 3) Computational solution of coupled crystal mechanics and dislocation kinetics. Comparison of dislocation dynamics predictions with experimental results in the area of statistical properties of dislocations and their field was also a part of the proposed effort. In the first research task, the dislocation dynamics simulation method was used to investigate the spatial, orientation, velocity, and temporal statistics of dynamical dislocation systems, and on the use of the results from this investigation to complete the kinetic description of dislocations. The second task focused on completing the formulation of a kinetic theory of dislocations that respects the discrete nature of crystallographic slip and the physics of dislocation motion and dislocation interaction in the crystal. Part of this effort also targeted the theoretical basis for establishing the connection between discrete and continuum representation of dislocations and the analysis of discrete dislocation simulation results within the continuum framework. This part of the research enables the enrichment of the kinetic description with information representing the discrete dislocation systems behavior. The third task focused on the development of physics-inspired numerical methods of solution of the coupled dislocation kinetics and crystal mechanics framework. To a large extent, this task has also been successfully started. We have developed a custom finite-element approach with mesh points being a subset of the underlying crystal structure. When used to predict the evolution of the dislocation system, the planar motion of dislocations is naturally captured for all slip systems, thus minimizing numerical errors and providing simple ways to investigate cross slip and dislocation reactions. Preliminary results in this direction show that we are closer than ever in building a predictive framework for dislocation dynamics and mesoscale plasticity based on the first principles of dislocation dynamics. The rest of the report gives and overview of the research performed under this project and highlights the key results and open questions left for future investigations.
The ability to predict the behavior of structural components in a nuclear power plant is critical to the nuclear industry. Structural metals in the primary loop of nuclear power plants must endure challenges such as irradiation and mechanical and thermal loading, and these structural metal components must continue to function in potential transient and accident conditions throughout the operational lifetime of the power plant. This extreme operational environment changes the metal microstructure by creating additional defects. The physical interactions of dislocations with these defects govern how the metal will respond to future conditions. Therefore predicting the mechanical response of these metals requires a set of physically based and reliable models of dislocation and defect interactions. These microstructure elements include glide mobile and immobile dislocations, geometrically necessary dislocations, twinning dislocations, irradiation defects, and thermal aging defects. We present here a continuum dislocation dynamics crystal plasticity framework to capture the interaction mechanisms of these dislocations and defects, verified with a combination of benchmark problems and comparisons with experimental data for two different types of structural metals: alpha iron and nickel-based alloys. In our simulations of $\alpha$ iron we highlight the advantages of applying a Monte Carlo stochastic model of cross slip dislocation motion and show the importance of capturing the 3D nature of glide dislocation and self-interstitial atom loop radiation defect interactions. We demonstrate coupling of glide dislocations with geometrically necessary dislocations to capture the influence of lattice bending, including the sensitivity of the geometrically necessary dislocations to changes in the grain boundary angle. We further examine the interaction of glide dislocations with the twin dislocations and thermally aged defects which have been observed in a nickel-based alloy with additional models. Finally we assess the reliability of this crystal plasticity framework by comparing two dislocation glide velocity models across the range of normal operation temperatures. In successfully applying our crystal plasticity framework to multiple metals, we provide further evidence of the reliability of our approach. The results of this mechanism-based continuum dislocation dynamics crystal plasticity framework can be used to inform engineering scale models throughout the nuclear industry.
This report describes the work in NEAMS (Nuclear Energy Advanced Modeling and Simulation) to develop a mechanistic fission gas model for uranium nitride fuels in BISON. The existing Sifgrs (Simple integrated fission gas release and swelling) model tracks the average properties of two bubble populations in the bulk and at the grain boundaries. It was recognized that dislocations play a crucial role in the fission gas swelling of UN, and an irradiation-induced dislocation density model was needed, as well as a model describing how fission gas interacts with dislocations creating a third population of bubbles along dislocations. A mechanistic model that tracks dislocation bubbles was implemented in Sifgrs. This model was used to simulate fission gas swelling and release in UN. Lower-length-scale calculations and experimental observations from carbide fuel were leveraged to help populate the model with essential parameters. As a placeholder, an empirical function was formulated for the evolution of the dislocation network. Because of the difference in evolution of this network at different temperatures, the dislocation bubbles are able to capture behavior that the bulk intragranular bubbles cannot. It was found that the bulk bubbles dominate microscopic swelling at low temperatures, and the dislocation bubbles dominate at higher temperatures. Based on this model, the transition between the two bubble types is the main factor behind the breakaway swelling phenomenon in UN. In order to test the model, a couple of assessments were run. For the lower temperature JOYO pins, the model produces reasonable fission gas release and swelling values. For the higher temperature SP1 pin, the model dramatically underestimates the fission gas release. This issue is attributed to the gas being trapped inside dislocation bubbles, unable to escape to grain boundaries to cause fission gas release, and could be remedied by a mechanistic dislocation model that allows the dislocation density to decrease at very high temperatures. A preliminary mechanistic dislocation model was developed supported by first-principle calculations. These calculations provided valuable insight into how interstitial defects cluster in UN in the {110} orientation, and may eventually form dislocation loops, leading to the conclusion that dislocation loops may nucleate from these clusters. An estimate of the dislocation line energy in UN was obtained and will be improved in future work. The free-energy cluster dynamics code Centipede was used to track interstitial and vacancy absorption at dislocation loops and calculate their growth. In addition, improvements to Centipede were made including new convergence criteria for transient simulations and sink driving force updates.
This report describes the current efforts to implement a model to predict the engineering scale behavior of irradiated and aged Ni 2 Cr model alloy and Alloy 690 into the Grizzly software. Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components. Engineering scale behavior of materials are strongly impacted by the underlying microstructure: the distribution and interaction of dislocations and defects within the individual grains dictates how the material will respond to loading. The evolution of these dislocations and defects determines how the engineering scale material will perform under different loading and environmental conditions. Mesoscale models can function as a bridge to connect microstructure evolution models to engineering models of the ductile-brittle transition temperature (DBTT) curve and fracture. The toughness of a material, which governs the engineering scale fracture behavior, is a function of both elastic deformation and plastic deformation. While the yield stress measurements capture the elastic contribution to a material’s toughness, information about the dislocation movement is required to understand the plastic component of toughness. Crystal plasticity models predict dislocation density evolution and thereby connect microstructure evolution to engineering scale behavior predictions. These models calculate the plastic strain within a metal as a function of dislocations and other crystal defects, including precipitates and interstitial loops, among others, by tracking dislocation movement and the interaction of these dislocations with the crystal defects. Within NiCr alloys, long-range ordered (LRO) precipitates form under aging conditions. These precipitates act as barriers to dislocation motion and are either sheared by dislocations or pin dislocations. Additional defects, such as stacking fault tetrahedra (SFTs) from irradiation damage, also impede dislocation motion in these alloys. A crystal plasticity mode for Ni 2 Cr model alloy and Alloy 690 is being currently implemented in the Grizzly code at INL. This report includes a description of the theoretical models incorporated into the Grizzly crystal plasticity code. These theoretical models include mechanisms for dislocation glade and twinning as well as dislocation-defect interaction relationship for both LRO precipitates and SFTs. We have selected these dislocation models based on the deformation mechanisms observed in experimental micropillar compression tests conducted on the Ni 2 Cr model alloy. The models are then extended to Alloy 690 with ongoing parameter fitting.
The objectives of this work are the following: (1) The application of micro and mesoscale modeling techniques to study dislocation properties in ferritic and W-based materials; and (2) The development of computational models and tools to study damage accumulation in >1 dpa (fusion-like) conditions, both for Fe and W-based alloys. The high-temperature strength of structural ferritic alloys (ferritic/martensitic steels, ODS steels, bcc refractory alloys) hinges on the thermal stability of second phase particles and their interactions with dislocations. Irradiation damage can modify the structure and stability of both the particles and dislocations, particularly by the introduction of gas atoms, point defects and point defect clusters. The three aspects of materials strength that we are studying are: (a) Computation of dislocation mobility functions (stress-velocity relations) as a function of temperature and dislocation character. This will be done via molecular dynamics (MD) simulations of single dislocation motion under applied shear stress. This is a fundamental input to dislocation dynamics (DD) simulations and also provides fundamental insights into the high-temperature plastic behavior of ferritic materials. (b) Simulations of dislocation-obstacle interactions using MD and DD. This subtask includes simulating the effect on dislocation glide of precipitates (e.g., α' Cr precipitates), ODS particles, and irradiation induced defect clusters (e.g. voids, dislocation loops, etc.). (c) Implementation of this information (dislocation mobilities and dislocation-defect interaction rules) into DD codes that will allow us to study plasticity of single crystals Fe alloys under relevant irradiation conditions.
The multiaxial ratcheting behaviors of polycrystalline 20 ordinary carbon steel were investigated at room temperature. The macroscopic experimental results showed that the studied multiaxial ratcheting depends greatly on the mean stress, stress amplitude and loading path. The axial ratcheting strain increased with the increase of applied mean stress and stress amplitude. Apparent additional hardening was observed in the non-proportionally multiaxial cyclic loading. The multiaxial ratcheting of 20 carbon steel was lower than the corresponding uniaxial one and varies with different loading paths. Dislocation patterns and their evolutions of the multiaxial ratcheting of different loading paths were then investigated using transmission electron microscopy. The obtained images showed that, with the increasing number of loading cycles, the dislocation patterns evolved from dislocation lines and networks to dislocation tangles, walls and cells. After certain cycles, sub-grains were formed because of the re-arrangement of dislocations in the walls of cells and inside the cells since the cross slip of dislocations can be easily activated for the 20 carbon steel, a kind of body-centered cubic metal. The dislocation evolution of the multiaxial ratcheting is much quicker than that of the uniaxial one. With the reference to the uniaxial one of 20 carbon steel, the macroscopic multiaxial ratcheting behaviors can be qualitatively correlated with the microscopic observation of the dislocation patterns and their evolution. - Highlights: • Multiaxial loading hardly changes the cyclic stable feature of 20 carbon steel. • Multiaxial ratcheting of 20 carbon steel depends greatly on the load path. • Dislocation patterns evolve quicker in the multiaxial case. • The stabilized dislocation pattern is sub-grain, rather than the dislocation cell. • Sub-grains formed after certain cycles make the stable ratcheting strain rate large.
In this work we study the microstruture evolution of a newly developed 12% Cr martensitic/ferritic steel in as-received condition and after creep at 650 °C under 130 MPa and 80 MPa. The microstructure is described as consisting of mobile dislocations, dipole dislocations, boundary dislocations, precipitates, lath boundaries, block boundaries, packet boundaries and prior austenitic grain boundaries. The material is characterized employing light optical microscopy (LOM), scanning electron microscopy (SEM), transmission electron microscopy (TEM), X-ray diffraction (XRD) and electron backscatter diffraction (EBSD). TEM is used to characterize the dislocations (mobile + dipole) inside the subgrains and XRD measurements are used to the characterize mobile dislocations. Based on the subgrain boundary misorientations obtained from EBSD measurements, the boundary dislocation density is estimated. The total dislocation density is estimated for the as-received and crept conditions adding the mobile, boundary and dipole dislocation densities. Additionally, the subgrain size is estimated from the EBSD measurements. In this publication we propose the use of three characterization techniques TEM, XRD and EBSD as necessary to characterize all type of dislocations and quantify the total dislocation densty in martensitic/ferritic steels. - Highlights: • Creep properties of a novel 12% Cr steel alloyed with Ta • Experimental characterization of different types of dislocations: mobile, dipole and boundary • Characterization and interpretation of the substructure evolution using unique combination of TEM, XRD and EBSD.
Abstract Direct Molecular Dynamics (MD) simulations are being increasingly employed to model dislocation-mediated crystal plasticity with atomic resolution. Thanks to the dislocation extraction algorithm (DXA), dislocation lines can be now accurately detected and positioned in space and their Burgers vector unambiguously identified in silico, while the simulation is being performed. However, DXA extracts static snapshots of dislocation configurations that by themselves present no information on dislocation motion. Referred to as a sweep-tracing algorithm (STA), here we introduce a practical computational method to observe dislocation motion and to accurately quantify its important characteristics such as preferential slip planes (slip crystallography). STA reconnects pairs of successive snapshots extracted by DXA and computes elementary slip facets thus precisely tracing the motion of dislocation segments from one snapshot to the next. As a testbed for our new method, we apply STA to the analysis of dislocation motion in large-scale MD simulations of single crystal plasticity in BCC metals. We observe that, when the crystal is subjected to uniaxial deformation along its [001] axis, dislocation slip predominantly occurs on the {112} maximum resolved shear stress plane under tension, while in compression slip is non-crystallographic (pencil) resulting in asymmetric mechanical response. The marked contrast in the observed slip crystallography is attributed to the twinning/anti-twinning asymmetry of shears in the {112} planes relatively favoring dislocation motion in the twinning sense while hindering dislocations from moving in the anti-twinning directions.