Molecular dynamics simulations of interactions between energetic dust and plasma-facing materials
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
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:
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)
- et al.
A review of dust in fusion devices: implications for safety and operational performance
Fusion Eng. Des.
(2002) - et al.
Iter, safety and licensing
Fusion Eng. Des.
(2007) - et al.
Dust generation mechanisms under powerful plasma impacts to the tungsten surfaces in iter elm simulation experiments
J. Nucl. Mater.
(2013) - et al.
Correlation of submicron dust production in DIII-D to impulsive wall heating from elms
J. Nucl. Mater.
(2009) - et al.
Vertical target and FW erosion during off-normal events and impurity production and transport during elms typical for ITER-FEAT
J. Nucl. Mater.
(2001) - et al.
Dust in magnetic confinement fusion devices and its impact on plasma operation
J. Nucl. Mater.
(1999) - et al.
Modeling of velocity distributions of dust in Tokamak edge plasmas and dust-wall collisions
J. Nucl. Mater.
(2009) - et al.
Carbon dust formation from re-deposited layers in high-density hydrogen/helium plasmas in the NAGDIS-II device
J. Nucl. Mater.
(2005) - et al.
Deuterium retention in carbon dust and carbon-tungsten mixed dust prepared by deuterium arc discharge
J. Nucl. Mater.
(2004) - et al.
Sputtering and reflection under cluster bombardment of solids
Nucl. Instrum. Methods Phys. Res. Sect. B-Beam Interact. Mater. Atoms
(2013)
Dust in magnetic fusion devices
Plasma Phys. Control. Fusion
Suppression of blister formation and deuterium retention on tungsten surface due to mechanical polishing and helium pre-exposure
Nucl. Fusion
Growth precursors and dynamics of dust particle formation in the ar/ch4 and ar/c2h2 plasmas
Plasma Sources Sci. Technol.
Collection strategy, inner morphology, and size distribution of dust particles in ASDEX upgrade
Nucl. Fusion
Divertor tungsten tile melting and its effect on core plasma performance
Nucl. Fusion
Dust studies in DIII-D and textor
Nucl. Fusion
Cited by (4)
Molecular dynamics simulation and modelling on the dust-wall interactions in tokamak
2022, Journal of Nuclear MaterialsCitation 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 AtomsCitation 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 AtomsCitation 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.
Analysis of hypervelocity impacts: The tungsten case
2022, Nuclear Fusion