Hydrogen is known to embrittle austenitic stainless steels, which are widely used in high-pressure hydrogen storage and delivery systems, but the mechanisms that lead to such material degradation are still being elucidated. The current work investigates the deformation behavior of single crystal austenitic stainless steel 316L through combined uniaxial tensile testing, characterization and atomistic simulations. Thermally precharged hydrogen is shown to increase the critical resolved shear stress (CRSS) without previously reported deviations from Schmid's law. Molecular dynamics simulations further expose the statistical nature of the hydrogen and vacancy contributions to the CRSS in the presence of alloying. Slip distribution quantification over large in-plane distances (>1 mm), achieved via atomic force microscopy (AFM), highlights the role of hydrogen increasing the degree of slip localization in both single and multiple slip configurations. The most active slip bands accumulate significantly more deformation in hydrogen precharged specimens, with potential implications for damage nucleation. For 〈110〉 tensile loading, slip localization further enhances the activity of secondary slip, increases the density of geometrically necessary dislocations and leads to a distinct lattice rotation behavior compared to hydrogen-free specimens, as evidenced by electron backscatter diffraction (EBSD) maps. The results of this study provide a more comprehensive picture of the deformation aspect of hydrogen embrittlement in austenitic stainless steels.
The critical stress for cutting of a void and He bubble (generically referred to as a cavity) by edge and screw dislocations has been determined for FCC Fe0.70Cr0.20Ni0.10—close to 300-series stainless steel—over a range of cavity spacings, diameters, pressures, and glide plane positions. The results exhibit anomalous trends with spacing, diameter, and pressure when compared with classical theories for obstacle hardening. These anomalies are attributed to elastic anisotropy and the wide extended dislocation core in low stacking fault energy metals, indicating that caution must be exercised when using perfect dislocations in isotropic solids to study void and bubble hardening. In many simulations with screw dislocations, cross-slip was observed at the void/bubble surface, leading to an additional contribution to strengthening. We refer to this phenomenon as cavity cross-slip locking, and argue that it may be an important contributor to void and bubble hardening.
Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.
The growth of helium bubbles impacts structural integrity of materials in nuclear applications. Understanding helium bubble nucleation and growth mechanisms is critical for improved material applications and aging predictions. Systematic molecular dynamics simulations have been performed to study helium bubble nucleation and growth mechanisms in Fe70Ni11Cr19 stainless steels. First, helium cluster diffusivities are calculated at a variety of helium cluster sizes and temperatures for systems with and without dislocations. Second, the process of diffusion of helium atoms to join existing helium bubbles is not deterministic and is hence studied using ensemble simulations for systems with and without vacancies, interstitials, and dislocations. We find that bubble nucleation depends on diffusion of not only single helium atoms, but also small helium clusters. Defects such as vacancies and dislocations can significantly impact the diffusion kinetics due to the trapping effects. Vacancies always increase the time for helium atoms to join existing bubbles due to the short-range trapping effect. This promotes bubble nucleation as opposed to bubble growth. Interestingly, dislocations can create a long-range trapping effect that reduces the time for helium atoms to join existing bubbles. This can promote bubble growth within a certain region near dislocations.
Tritium population thermodynamics and transport kinetics critically define the tritium storage performance of zirconium tritides that can be used for a variety of nuclear applications including tritium-producing burnable absorber rods. Both thermodynamic and kinetic properties can be sensitive to grain sizes of materials and can be significantly altered by irradiated defects during operation under the reactor environments. A thorough experimental characterization of how these properties evolve under different reactor conditions and different initial grain structures is extremely challenging. Here molecular dynamics simulations are used to investigate tritium population and diffusion in zirconium with and without different planar symmetric and asymmetric tilt grain boundaries and irradiated defects. Here, we found that in addition to trapping tritium, the most significant effect of planar grain boundaries is to increase tritium diffusivity on the boundary plane. Furthermore, fine grain structures are found to mitigate the change of tritium diffusivity due to irradiated point defects as these point defects are likely to migrate to and sink at grain boundaries.
The interplay between hydrogen and dislocations (e.g., core and elastic energies, and dislocation-dislocation interactions) has implications on hydrogen embrittlement but is poorly understood. Continuum models of hydrogen enhanced local plasticity have not considered the effect of hydrogen on dislocation core energies. Energy minimization atomistic simulations can only resolve dislocation core energies in hydrogen-free systems because hydrogen motion is omitted so hydrogen atmosphere formation can’t occur. Additionally, previous studies focused more on face-centered-cubic than body-centered-cubic metals. Discrete dislocation dynamics studies of hydrogen-dislocation interactions assume isotropic elasticity, but the validity of this assumption isn’t understood. We perform time-averaged molecular dynamics simulations to study the effect of hydrogen on dislocation energies in body-centered-cubic iron for several dislocation character angles. We see atmosphere formation and highly converged dislocation energies. We find that hydrogen reduces dislocation core energies but can increase or decrease elastic energies of isolated dislocations and dislocation-dislocation interaction energies depending on character angle. We also find that isotropic elasticity can be well fitted to dislocation energies obtained from simulations if the isotropic elastic constants are not constrained to their anisotropic counterparts. These results are relevant to ongoing efforts in understanding hydrogen embrittlement and provide a foundation for future work in this field.
Comprehensive molecular dynamics tensile test simulations have been performed to study the delamination processes of seven different grain boundaries / cleavage planes (Σ1{111}, Σ3{111}, Σ5{100}, Σ7{111}, Σ9{411}, Σ11{311}, and R{100}/{411}) containing a helium bubble. Combinations of a variety of conditions are explored including different strain rates, system dimensions, bubble density, bubble radius, bubble pressure, and temperature. We found that in general, grain boundaries absorb less energies with decreasing strain rate but increasing bubble areal density, bubble pressure, bubble radius, and temperature. The propensity of grain boundary delamination is sensitive to grain boundary type: The random grain boundary R{100}/{411} is one of the most brittle boundaries whereas the Σ1{111} cleavage plane and the Σ3{111} twin boundary are two of the toughest boundaries. The sorted list of grain boundary fracture vulnerability obtained from our dynamic tensile test simulations differs from the one obtained from our decohesion energy calculations, confirming the important role of plastic deformation during fracture. Detailed mechanistic analyses are performed to interpret the simulated results.
We explore the character angle dependence of dislocation-solute interactions in a face-centered cubic random Fe0.70Ni0.11Cr0.19 alloy through molecular dynamics (MD) simulations of dislocation mobility. Using the MD mobility data, we determine the phonon and thermally activated solute drag parameters which govern mobility for each dislocation character angle. The resulting parameter set indicates that, surprisingly, the solute energy barrier does not depend on character angle. Instead, only the zero-temperature flow stress—which is dictated by the activation area for thermal activation—is dependent on character angle. By analyzing the line roughness from MD simulations and the geometry of a bowing dislocation line undergoing thermal activation, we conclude that the character angle dependence of the activation area in this alloy is governed by the dislocation line tension, rather than the dislocation-solute interaction itself. Our findings motivate further investigation into the line geometry of dislocations in solid solutions.
Tritium diffusion in α-Zr containing point defects such as vacancies or self-interstitial atoms (SIAs) is simulated using molecular dynamics. Point defects rapidly aggregate to form extended defects, such as 3D nanoclusters and Frank loops. The geometry of extended defects is affected by the presence of tritium. At low temperature and in the absence of tritium, vacancies aggregate to form stacking fault pyramids. Addition of tritium at these temperatures promotes aggregation of vacancies to form 3D nanoclusters, within which the tritium concentration can be sufficiently high to suggest that these defects may serve as nucleation sites for hydride precipitation. Trapping of tritium in vacancy nanocluster reduces the calculated bulk diffusivity by an amount proportional to the vacancy concentration. At high temperature, vacancy clusters change shape to form planar basal dislocation loops, which bind tritium less strongly, leading to a sharp reduction in the fraction of trapped tritium and a corresponding increase in tritium diffusivity at high temperature. In contrast, SIAs increase tritium diffusion through α-Zr. Analysis of atomic trajectories shows that tritium does not interact directly with SIAs. In conclusion, diffusion enhancement is instead related to expansion of the lattice.
Our ability to shape and finish a component by combined methods of fabrication including (but not limited to) subtractive, additive, and/or no theoretical mass-loss/addition during the fabrication is now popularly known as solid freeform fabrication (SFF). Fabrication of a telescope mirror is a typical example where grinding and polishing processes are first applied to shape the mirror, and thereafter, an optical coating is usually applied to enhance its optical performance. The area of nanomanufacturing cannot grow without a deep knowledge of the fundamentals of materials and consequently, the use of computer simulations is now becoming ubiquitous. This article is intended to highlight the most recent advances in the computation benefit specific to the area of precision SFF as these systems are traversing through the journey of digitalization and Industry-4.0. Specifically, this article demonstrates that the application of the latest materials modelling approaches, based on techniques such as molecular dynamics, are enabling breakthroughs in applied precision manufacturing techniques.
The fundamental interactions between an edge dislocation and a random solid solution are studied by analyzing dislocation line roughness profiles obtained from molecular dynamics simulations of Fe0.70Ni0.11Cr0.19 over a range of stresses and temperatures. These roughness profiles reveal the hallmark features of a depinning transition. Namely, below a temperature-dependent critical stress, the dislocation line exhibits roughness in two different length scale regimes which are divided by a so-called correlation length. This correlation length increases with applied stress and at the critical stress (depinning transition or yield stress) formally goes to infinity. Above the critical stress, the line roughness profile converges to that of a random noise field. Motivated by these results, a physical model is developed based on the notion of coherent line bowing over all length scales below the correlation length. Above the correlation length, the solute field prohibits such coherent line bow outs. Using this model, we identify potential gaps in existing theories of solid solution strengthening and show that recent observations of length-dependent dislocation mobilities can be rationalized.