Elsevier

Fusion Engineering and Design

Volume 100, November 2015, Pages 220-225
Fusion Engineering and Design

Molecular dynamics simulations of interactions between energetic dust and plasma-facing materials

https://doi.org/10.1016/j.fusengdes.2015.05.079Get rights and content

Abstract

The interactions between dust and plasma-facing material (PFM) relate to the lifetime of PFM and impurity production. Series results have been obtained theoretically and experimentally but more detailed studies are needed. In present research, we investigate the evolution of kinetic, potential and total energy of plasma-facing material (PFM) in order to understand the dust/PFM interaction process. Three typical impacting energy are selected, i.e., 1, 10 and 100 keV/dust for low-, high- and hyper-energy impacting cases. For low impacting energy, dust particles stick on PFM surface without damaging it. Two typical time points exist and the temperature of PFM grows all the time but PFM structure experience a modifying process. Under high energy case, three typical points appear. The temperature curve fluctuates in the whole interaction process which indicates there are dust/PFM and kinetic/potential energy exchanges. In the hyper-energy case in present simulation, the violence dust/PFM interactions cause sputtering and crater investigating on energy evolution curves. We further propose the statistics of energy distribution. Results show that about half of impacting energy consumes on heating plasma-facing material meanwhile the other half on PFM structure deformation. Only a small proportion becomes kinetic energy of interstitial or sputtering atoms.

Introduction

Over the past few years progress towards the understanding of dust behaviors in fusion devices has been significant through experiments, modelling and model validation in it [1], [2], [3]. Many plasma processes such as flaking of blisters [4], volume polymerization [5], cracking development [6], arcing [7], brittle destruction [3] and surface melting [9] even ELM and disruption events [10], [11] contribute to the production of dust particles. Dust particles increases the operation difficulties of plasma since the dusts can contaminate the core plasma [12]. The contamination can increase the amount of impurity concentration in plasma or even lead to disruption [13], [14]. Dust particles can also worsen the operation safety of fusion devices, including the risk of exploding and tritium leaking. As a result dust generation and its transport have been attached great importance to.

Because of the dust-PFM (plasma-facing materials) interaction associated with the life time of PFM and impurity generations, it is also a crucial issue in fusion science. Impacting velocities, angle, materials of dust, number of impacting dust, etc. are the factors influencing dust-PFM interaction. Although some studies are being performed for the interaction, it remains a major unsolved problem. One of the reasons is lacking of diagnose methods on dust-PFM interaction process. To the author's knowledge, few effective diagnose techniques can investigate dust-PFM interaction at present. Our knowledge on the interaction is mainly from numeric simulation results, such as finite element method or molecular dynamics [15], [16], [17]. Moreover, previous dust studies basically focus on carbon dust owing to the widely use of carbon plasma-facing materials (such as CFC or graphite) [18], [19]. However, the characteristics of carbon and metal dust are quite different. Carbon materials are active in hydrogen's isotopes which is used as fusion fuel [14], [21], [22], [23], [24], [25]. The chemical erosion is the primary erosion mechanism for carbon materials. But metal materials are more stable in hydrogen environment (especially tungsten) thus the physical sputtering caused by heavy ions or dust particles becomes a more important erosion process than chemical erosion comparing to carbon PFM. As metal materials will be used in ITER, the research on metal dust-PFM interactions has become necessary and urgent.

The rest of this paper is organized as follows. Section 2, we briefly overview the molecular dynamics (MD) simulation methods. In Section 3, we present dust-PFM interaction process and investigate the effect of tungsten (W) dust particles impacting into W surfaces. Finally, we summarize our results and provide concluding remarks in Section 4.

Section snippets

Methods

The interactions between dust particle and PFM surface were simulated with classical molecular dynamics (MD) employing many-body interatomic potentials implemented in the LAMMPS code (Large-scale Atomic/Molecular Massively Parallel Simulator) [26]. Both of the materials of dust particles and surface are W. The embedded atom method (EAM) potential, namely the Finnis–Sinclair (FS) potential [27] is used to describe the W–W interactions.

The total energy of FS potential is

UFS=UN+UP

M. W. Finnis and

Energy evolutions

We investigate the temporal evolution of three average energies of PFM atoms: kinetic energy 〈E〉, potential energy 〈U〉 and total energy 〈T〉. These energies are defined as follows:

E(t)=1NPFMiNPFMEi(t)U(t)=1NPFMiNPFMUi(t)T(t)=E(t)+U(t)here N is the total number of PFM atoms, Ei, Ui and Ti are the kinetic, potential and total energy of atom i.

As we have confirmed that the velocity of mass-center of PFM is zero, 〈E〉 shows the temperature of PFM atoms. 〈U〉 is closely related to the

Conclusions

In present research, we investigate the temporal energy evolution of plasma-facing material in order to learn about the details of dust/plasma-facing material interactions. Three typical impacting energies are selected for low-, high- and hyper-velocity dust. Under different impacting energies, the temporal evolutions are not the same. For low impacting energy, two typical time points exist and the temperature of PFM grows all the time but PFM structure experiences one modification. Under high-

Acknowledgement

This work was supported by National Magnetic Confinement Fusion Science Program under Contracts Nos. 2013GB105001, 2013GB105002 and 2013GB109001, National Natural Science Foundation of China (NSFC) with Grant Nos. 11205198, 11305213 and 11405201, as well as Technological Development Grant of Hefei Science Center of CAS under contract No. 2014TDG-HSC003. The authors gratefully acknowledge the operation team of HIRFL for their help. Authors would also like to thank Yukihiro Tomita and Gakushi

References (34)

  • S.I. Krasheninnikov et al.

    Dust in magnetic fusion devices

    Plasma Phys. Control. Fusion

    (1999)
  • D. Nishijima et al.

    Suppression of blister formation and deuterium retention on tungsten surface due to mechanical polishing and helium pre-exposure

    Nucl. Fusion

    (2005)
  • S. Hong et al.

    Growth precursors and dynamics of dust particle formation in the ar/ch4 and ar/c2h2 plasmas

    Plasma Sources Sci. Technol.

    (2003)
  • M. Balden et al.

    Collection strategy, inner morphology, and size distribution of dust particles in ASDEX upgrade

    Nucl. Fusion

    (2014)
  • B. Lipschultz et al.

    Divertor tungsten tile melting and its effect on core plasma performance

    Nucl. Fusion

    (2012)
  • A.A. Tuccillo, A. Alekseyev, B. Angelini, S.V. Annibaldi, M.L. Apicella, G. Apruzzese, J. Berrino, E. Barbato, A....
  • D.L. Rudakov et al.

    Dust studies in DIII-D and textor

    Nucl. Fusion

    (2009)
  • Cited by (4)

    • Molecular dynamics simulation and modelling on the dust-wall interactions in tokamak

      2022, Journal of Nuclear Materials
      Citation Excerpt :

      Smoothed particle hydrodynamics method has been employed to study DWI [17], which indicated that beryllium and tungsten dust grains can produce significant damage on beryllium wall. Molecular dynamics (MD) simulation has been applied for tungsten dust with the velocity from several hundreds of m/s to 3 km/s [18,19]. The modification of surfaces and shallow layers by energetic dust is stated [20].

    • Study of crater formation and its characteristics due to impact of a cluster projectile on a metal surface by molecular dynamics approach

      2017, Nuclear Instruments and Methods in Physics Research, Section B: Beam Interactions with Materials and Atoms
      Citation Excerpt :

      The authors studied the carbon trapping rate, implantation depth and sputtering yield of tungsten due to the bombardment of deuterium and carbon ions at varying energies. To understand the dust and plasma facing material interaction Guo-jian Niu et al. [19] investigated the evolution of kinetic energy, potential energy and total energy of the tungsten based plasma facing material bombarded with tungsten dust particle consists of 541 atoms with energies 1, 10 and 100 keV/dust particle. In the present work, we have investigated crater formation on a tungsten plasma facing material and its sputtering characteristics by bombarding it with tungsten dust/cluster particles of three different sizes (formed of 14, 32, 64 atoms) at varying energies ranging from 0.01–10 keV/dust particle.

    • The influence of landing points on the sputtering of mono-crystal solids due to cluster impacting

      2017, Nuclear Instruments and Methods in Physics Research, Section B: Beam Interactions with Materials and Atoms
      Citation Excerpt :

      The hyper-velocity impact of interstellar nano-grains is also a potential field of sputtering theory for the erosion features of the surface of spacecraft [9], thrusters [10] and other astronomical objects [11]. The plasma contamination owning to the sputtering of plasma-facing materials is a focused problem for fusion science[12]. Although both experiments [13–15] and simulations [16–19] have been performed for the interactions between various clusters and solids [20–22], the understanding of cluster sputtering process is still not clear enough comparing to the well-developed mono-atomic sputtering theory.

    View full text