Search NASA⌕ Search

Engineering topics

Mathew, Nithin

Publications and source records attributed to Mathew, Nithin.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Structure and migration of heavily irradiated grain boundaries and dislocations in Ni in the athermal limit

The microstructural evolution at and near preexisting grain boundaries (GBs) and dislocations in materials under high radiation doses is still poorly understood. In this work, we use the creation relaxation algorithm (CRA) developed for atomistic modeling of high-dose irradiation in bulk materials to probe the athermal limit of saturation of GB and dislocation core regions under irradiation in fcc Ni. We find that, upon continuously subjecting a single dislocation or GB to Frenkel pair creation in the athermal limit, a local steady-state disordered defect structure is reached with excess properties that fluctuate around constant values. Case studies are given for a straight screw dislocation which elongates into a helix under irradiation and several types of low- and high-angle GBs, which exhibit coupled responses such as absorption of extrinsic dislocations, roughening, and migration. A positive correlation is found between the initial GB energy and the local steady-state GB energy under irradiation across a wide variety of GB types. Metastable GB structures with similar density in the defect core region but different initial configurations are found to converge to the same limiting structure under CRA. The mechanical responses of pristine and irradiated dislocations and GB structures are compared under an applied shear stress. Irradiated screw and edge dislocations are found to exhibit a hardening response, migrating at larger flow stresses than their pristine counterparts. Mobile GBs are found to exhibit softening or hardening responses depending on GB character. Although some GBs recover their initial pristine structures upon migration outside of the radiation zone, many GBs sustain different flow stresses corresponding to altered mobile core structures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-driven modeling of dislocation mobility from atomistics using physics-informed machine learning

Dislocation mobility, which dictates the response of dislocations to an applied stress, is a fundamental property of crystalline materials that governs the evolution of plastic deformation. Traditional approaches for deriving mobility laws rely on phenomenological models of the underlying physics, whose free parameters are in turn fitted to a small number of intuition-driven atomic scale simulations under varying conditions of temperature and stress. This tedious and time-consuming approach becomes particularly cumbersome for materials with complex dependencies on stress, temperature, and local environment, such as body-centered cubic crystals (BCC) metals and alloys. In this paper, we present a novel, uncertainty quantification-driven active learning paradigm for learning dislocation mobility laws from automated high-throughput large-scale molecular dynamics simulations, using Graph Neural Networks (GNN) with a physics-informed architecture. We demonstrate that this Physics-informed Graph Neural Network (PI-GNN) framework captures the underlying physics more accurately compared to existing phenomenological mobility laws in BCC metals.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A combined kinetic Monte Carlo and phase field approach to model thermally activated dislocation motion

In this work, a kinetic Monte Carlo (KMC) scheme is integrated into a phase field dislocation dynamics (PFDD) model to account for thermal activation. While the implementation is general and could be applied to many different thermally activated processes, we apply the model to simulate a dislocation transmitting through an interface. Two cases are analyzed: a dislocation transmitting through a strong interface, and a dislocation transmitting through a more realistic grain boundary in tungsten. In conclusion, results and trends for both cases are discussed in detail, and ultimately show that the combined KMC-PFDD algorithm is effective in capturing thermally activated dislocation processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transition in helium bubble strengthening of copper from quasi-static to dynamic deformation

Damage from low-temperature irradiation and the subsequent degradation of materials performance pose significant challenges for the storage of radioactive materials and for peripheral components in some nuclear reactor designs. Fully understanding the mechanical behavior of such materials requires test data for strain rates in both the quasi-static (< 10/s) and dynamic (>> 10/s) regimes. While dynamic testing has generally been avoided in the past for neutron irradiated (contamination concerns) and ion irradiated (insufficient volume) materials, surface-sensitive Richtmyer-Meshkov instability (RMI) tests were used in the present work to overcome these limitations. Here, nanopillar compression, nanoindentation, and RMI testing data from a helium implanted surface layer (~10 µm thick) were compiled to explore the effects of helium bubbles on the materials strength of high-purity copper at strain rates of 0.001/s – 10 8 /s. While nano-mechanical testing revealed increases in yield strength and hardness with increasing helium dose from 1000 to 4000 appm He, RMI indicated no significant changes in strength as compared to unimplanted copper. Here, this discrepancy in behavior was rationalized through a combination of recent literature and follow-on molecular dynamics (MD) simulations, leading to the conclusion that the nanoscale helium bubbles acting as dispersed barriers to dislocation motion at quasi-static strain rates collapse under shock loading and cease to be effective barriers at high strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interstitial hydrogen enhances the mobility of some grain boundaries in tungsten

Segregation of interstitials at a grain boundary (GB) is known to generally lower its mobility. This phenomenon, called 'solute-drag', has important ramifications on the process of recrystallization and microstructural evolution. Here, we present predictions from molecular dynamics (MD) simulations which demonstrate that interstitial hydrogen in tungsten can in fact increase the mobility of some GBs which exhibit shear coupling. Assuming a disconnection-based mechanism, activation energies and pre-factors for disconnection nucleation are predicted from simulations of shear-coupled motion. In GBs where enhanced mobility is predicted, interstitial H reduces both the activation energy and the pre-factor for disconnection nucleation, thus effectively increasing the mobility. For GBs with diminished mobility, MD predicts that presence of interstitial H reduces the pre-factor and, in some cases, increases the activation energy. The reduction in the activation energy inferred from MD simulations are confirmed by nudged elastic band calculations. Temperature-dependent structural transitions are observed for some GBs, and the effect of interstitial H is found to change with the changes in structure. The effect of interstitial H is predicted to be complex and highly variable, providing some plausible explanations for experimental observations on the recrystallization of tungsten in presence of H-loaded plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Thermal gradient effect on helium and self-interstitial transport in tungsten

First-wall materials in a fusion reactor are expected to withstand harsh conditions, with high heat and particle fluxes that modify the materials microstructure. These fluxes will create strong gradients of temperature and concentration of diverse species. Besides the He ash and the hydrogenic species, neutron particles generated in the fusion reaction will collide with the material creating intrinsic defects, such as vacancies, self-interstitials atoms (SIAs), and clusters of such point defects. These defects and the He atoms will then migrate in the presence of the aforementioned gradients. In this study, we use nonequilibrium molecular dynamics to analyze the transport of He and SIAs in the presence of a thermal gradient in tungsten. We observe that, in all cases, the defects and impurity atoms tend to migrate toward the hot regions of the tungsten sample. The resulting species concentration profiles are exponential distributions, rising toward the hot regions of the sample, in agreement with irreversible thermodynamics analysis. For both He atoms and SIAs, we find that the resulting species flux is directed opposite to the heat flux, indicating that species transport is governed by a Soret effect (thermal-gradient-driven diffusion) characterized by a negative heat of transport that drives species diffusion uphill (from the cooler to the hot regions of the sample). Here, we demonstrate that the steady-state species profiles obtained accounting for the Soret effect vary significantly from those where temperature-gradient-driven transport is not considered and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene

All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

A 3D phase field dislocation dynamics model for body-centered cubic crystals

In this work, we present a 3D Phase Field Dislocation Dynamics (PFDD) model for body-centered cubic (BCC) metals. The model formulation is extended to account for the dependence of the Peierls barrier on the line-character of the dislocation. Simulations of the expansion of a dislocation loop belonging to {110} < 111 > the slip system are presented with direct comparison to Molecular Statics (MS) simulations. The extended PFDD model is able to capture the salient features of dislocation loop expansion predicted by MS simulations. In conclusion, the model is also applied to simulate the motion of a straight screw dislocation through kink-pair motion.

36 MATERIALS SCIENCE↗