Computational Analysis Division

BARC, India

Computational Analysis Division

BARC, India
Time filter
Source Type

Ikkurthi V.R.,Computational Analysis Division | Hemani H.,Computational Analysis Division | Sugandhi R.,Indian Institute for Plasma Research | Rawat S.,Indian Institute of Technology Kharagpur | And 3 more authors.
Procedia Engineering | Year: 2017

A multi-scale model has been developed to generate and validate spallation (dynamic fracture) model parameters for metallic single crystals. The model uses Molecular Dynamic (MD) simulations at the atomic scales to obtain parameters for a dynamic fracture model called the Nucleation and Growth (NAG) model. The parameters obtained from MD are then used to describe the nucleation and growth of voids in macroscopic hydrodynamic simulations of material impact and spallation. The spall strength and thickness have earlier been computed and compared with the experimental results for single crystal copper (Cu), Molybdenum (Mo) and Niobium (Nb). A good match with the experiments was obtained for Cu and Mo. The mismatch in case of Nb has been dealt in present work. The potential used in MD simulations for Nb has been replaced by a more sophisticated stiff potential and a more accurate Equation of State (EOS) has been used in the hydrodynamic simulations. The spall model parameters obtained this way result in a good match between hydrodynamic simulations and experiments. It is seen that the growth threshold for voids is the important factor in determining the spall strength of single crystals rather than the nucleation threshold for voids. © 2017 The Authors.

Bukkuru S.,Andhra University | Bhardwaj U.,Computational Analysis Division | Rao A.D.P.,Andhra University | Warrier M.,Computational Analysis Division
Journal of Physics: Conference Series | Year: 2016

Molecular Dynamics (MD) simulations are suitable to study the details of atomistic processes. Identification of interstitials is one of the main difficulties in the analysis of defect dynamics. To identify the interstitials, different methods use different criteria. All the methods, except the Wigner-Seitz method, assume some input parameters which may lead to an erroneous number of interstitials. We describe an unsupervised clustering algorithm, called Max Space clustering method, to identify interstitials without assuming any input parameters. We show that the algorithm is suited for identifying interstitials at high temperatures when amplitude of atomic oscillations become larger, a fact ignored by other models. © Published under licence by IOP Publishing Ltd.

Warrier M.,Computational Analysis Division | Bhardwaj U.,Computational Analysis Division | Hemani H.,Computational Analysis Division | Schneider R.,University of Greifswald | And 2 more authors.
Journal of Nuclear Materials | Year: 2015

We report on molecular Dynamics (MD) simulations carried out in fcc Cu and bcc W using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code to study (i) the statistical variations in the number of interstitials and vacancies produced by energetic primary knock-on atoms (PKA) (0.1-5 keV) directed in random directions and (ii) the in-cascade cluster size distributions. It is seen that around 60-80 random directions have to be explored for the average number of displaced atoms to become steady in the case of fcc Cu, whereas for bcc W around 50-60 random directions need to be explored. The number of Frenkel pairs produced in the MD simulations are compared with that from the Binary Collision Approximation Monte Carlo (BCA-MC) code SDTRIM-SP and the results from the NRT model. It is seen that a proper choice of the damage energy, i.e. the energy required to create a stable interstitial, is essential for the BCA-MC results to match the MD results. On the computational front it is seen that in-situ processing saves the need to input/output (I/O) atomic position data of several tera-bytes when exploring a large number of random directions and there is no difference in run-time because the extra run-time in processing data is offset by the time saved in I/O. © 2015 Elsevier B.V.

Bhardwaj U.,Computational Analysis Division | Bukkuru S.,Andhra University | Warrier M.,Computational Analysis Division
Journal of Computational Physics | Year: 2016

An algorithm to rigorously trace the interstitialcy diffusion trajectory in crystals is developed. The algorithm incorporates unsupervised learning and graph optimization which obviate the need to input extra domain specific information depending on crystal or temperature of the simulation. The algorithm is implemented in a flexible framework as a post-processor to molecular dynamics (MD) simulations. We describe in detail the reduction of interstitialcy diffusion into known computational problems of unsupervised clustering and graph optimization. We also discuss the steps, computational efficiency and key components of the algorithm. Using the algorithm, thermal interstitialcy diffusion from low to near-melting point temperatures is studied. We encapsulate the algorithms in a modular framework with functionality to calculate diffusion coefficients, migration energies and other trajectory properties. The study validates the algorithm by establishing the conformity of output parameters with experimental values and provides detailed insights for the interstitialcy diffusion mechanism. The algorithm along with the help of supporting visualizations and analysis gives convincing details and a new approach to quantifying diffusion jumps, jump-lengths, time between jumps and to identify interstitials from lattice atoms. © 2015 Elsevier Inc.

Warrier M.,Computational Analysis Division | Pahari P.,Computational Analysis Division | Chaturvedi S.,Computational Analysis Division
Journal of Molecular Modeling | Year: 2015

Molecular dynamics (MD) simulations of high velocity impact (1–6 km/s) of RDX crystal with a nanometer-sized void, has been carried out to understand the mechanism of increase in temperature at void locations under shock loading. Similar simulations are then carried out on single-crystal copper for better interpretation of the results. A reactive potential that can simulate chemical reactions (ReaxFF) has been used for RDX, whereas an EAM potential has been used for Cu. Increased temperature at the void locations are observed under shock loading. The atomic motion, temperature, average potential energy per atom (PE), and average kinetic energy per atom (KE) in and around the voids are closely monitored in order to understand the reason for temperature increase. We compare our results with existing proposed mechanisms and show that some of the proposed mechanisms are not necessary for increased temperature at a void location. It is shown that the directed particle velocity is efficiently is converted into randomized velocity due to the presence of voids thereby increasing the local temperature transiently. In this initial stage (few picoseconds) of the shock, chemical reactions of energetic materials do not play a part in the temperature rise. © 2015, Springer-Verlag Berlin Heidelberg.

Abhishek A.,Indian Institute for Plasma Research | Warrier M.,Computational Analysis Division | Ganesh R.,Indian Institute for Plasma Research
Transactions of the Indian Institute of Metals | Year: 2015

Helium has very low solubility in metals leading to formation of helium clusters and complexes with vacancies/interstitials. Clusters and complexes then diffuse and coalesce to form helium bubbles of nanometric size or higher. These bubbles are one of the reasons for the brittle failure of materials. In order to study the diffusion process, molecular dynamics simulations of helium diffusion in FeCr grain boundaries are carried out for different orientations of the bicrystal. A total of eight different configurations of the bicrystal, have been studied in the temperature range of 700–1,000 K. Potential energy analysis of host and grain boundary (GB) atoms predicts average GB potential energy in the range of −3.7 to −3.6 eV, which is 0.2 eV higher than that in the host matrix. The migration energy of helium in grain boundary is found to be 0.3–0.7 eV, which is an order of magnitude larger than that in the bulk crystal. The grain boundary width calculated for all the bi-crystals lie within 11–14 Å and the cage distance of helium atom is of the order of the bond length (2.87 Å) of the host atoms. For ∑9〈110〉{1 1 4} orientation more than 100 helium trajectories have been analysed to measure the statistical variation of migration energy. © The Indian Institute of Metals - IIM 2015.

Abhishek A.,Indian Institute for Plasma Research | Warrier M.,Computational Analysis Division | Ganesh R.,Indian Institute for Plasma Research | Caro A.,Los Alamos National Laboratory
Journal of Nuclear Materials | Year: 2016

Helium(He) produced by transmutation process inside structural material due to neutron irradiation plays a vital role in the degradation of material properties. We have carried out Molecular dynamics(MD) simulations to study the growth of He bubble in Iron-Chromium alloy. Simulations are carried out at two different temperatures, viz. 0.1 K and 800 K, upto He bubble radius of 2.5 nm. An equation for variation of volume of He bubbles with the number of He atoms is obtained at both the temperatures. Bubble pressure and potential energy variation is obtained with increasing bubble radius. Dislocations are also found to be emitted after the bubble reaches a critical radius of 0.39 nm at 800 K. Separate MD simulations of He with pre-created voids are also carried out to study the binding energies of He and Vacancy (V) to Hem-Vn cluster. Binding energies are found to be in the range of 1-5.5 eV. © 2016 Elsevier B.V.

Rawat S.,Refueling Technology Division | Warrier M.,Computational Analysis Division | Chaturvedi S.,Computational Analysis Division | Chavan V.M.,Refueling Technology Division
AIP Conference Proceedings | Year: 2011

Molecular Dynamics simulations have been performed for uniform triaxial expansion of single-crystal copper at high strain-rate (5×10 9/s) with different initial temperatures. We have used the LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) code for this study. It is found that the macroscopic Nucleation and Growth (NAG) model is valid at the atomistic scale. There is a monotonic decrease in the void nucleation and growth thresholds with increase in temperature. When we approach the melting point of copper, we observe a double-dip in the pressure-time profile. The analysis of this double-dip shows that the first minimum corresponds to the loss of long-range order but there is no void nucleation at this juncture. The second minimum corresponds to the void nucleation and incipient growth of voids. © 2011 American Institute of Physics.

Rawat S.,Refueling Technology Division | Warrier M.,Computational Analysis Division | Chaturvedi S.,Computational Analysis Division | Chavan V.M.,Refueling Technology Division
AIP Conference Proceedings | Year: 2012

Impact of single crystal copper plates at 1100 m/s impact velocity has been simulated using the molecular dynamics code LAMMPS [1]. We validate the void nucleation threshold obtained under tri-axial strain conditions [2] by high velocity impact of copper plates. In this set up, the flyer remains in intimate contact with the target before the impact. We have also performed the impact of non-contact configuration of the flyer-target system, where the flyer starts a small distance away from the target before impact. Due to small misalignment of crystal planes in the flyer & target, the shock produced by impact creates defects in the target during its propagation. These defects become preferred void nucleation sites at the time of interaction of release waves and lead to spall at a lower value of tensile pressure. © 2012 American Institute of Physics.

Loading Computational Analysis Division collaborators
Loading Computational Analysis Division collaborators