We present a fully-implicit electromagnetic Particle-In-Cell Monte Carlo collision code, called NINJA, written for the simulation of inductively coupled plasmas. NINJA employs a kinetic enslaved ...Jacobian-Free Newton Krylov method to solve self-consistently the interaction between the electromagnetic field generated by the radio-frequency coil and the plasma response. The simulated plasma includes a kinetic description of charged and neutral species as well as the collision processes between them. The algorithm allows simulations with cell sizes much larger than the Debye length and time steps in excess of the Courant–Friedrichs–Lewy condition whilst preserving the conservation of the total energy. The code is applied to the simulation of the plasma discharge of the Linac4 H− ion source at CERN. Simulation results of plasma density, temperature and EEDF are discussed and compared with optical emission spectroscopy measurements. A systematic study of the energy conservation as a function of the numerical parameters is presented.
The present status of kinetic modeling of particle dynamics in hydrogen negative ion (H−) source plasmas and their comparisons with experiments are reviewed and discussed with some new results. The ...main focus is placed on the following topics, which are important for the research and development of H− sources for intense and high-quality H− ion beams: (i) effects of non-equilibrium features of electron energy distribution function on volume and surface H− production, (ii) the origin of the spatial non-uniformity in giant multi-cusp arc-discharge H− sources, (iii) capacitive to inductive (E to H) mode transition in radio frequency-inductively coupled plasma H− sources and (iv) extraction physics of H− ions and beam optics, especially the present understanding of the meniscus formation in strongly electronegative plasmas (so-called ion-ion plasmas) and its effect on beam optics. For these topics, mainly Japanese modeling activities, and their domestic and international collaborations with experimental studies, are introduced with some examples showing how models have been improved and to what extent the modeling studies can presently contribute to improving the source performance. Close collaboration between experimental and modeling activities is indispensable for the validation/improvement of the modeling and its contribution to the source design/development.
Abstract
It is reported from the measurement of a negative ion beamlet emittance that the characteristic three-Gaussian components are observed in the emittance diagram. The negative ion trajectories ...are calculated using 3D PIC-MCC simulation in order to clarify the origin of this complicated phase space structure and underlying physical mechanism. The characteristic three-Gaussian components can be reproduced in the emittance diagram from the simulation result. It is verified that the Gaussian components are caused by the negative ions extracted from the different plasma meniscus region, that is, the central region or the region near the edges of the meniscus.
Abstract
Sequences of hydrogen (H) and deuterium (D) experiments have been done by NIFS research and development negative ion source (RNIS) for the deuterium NBI development. In the experiments, the ...co-extracted electron current with the negative ions and the electron density in the plasma generation region in the D experiment have been around three times higher than that in the H experiment. To explain the difference of the electron density in the RNIS driver region, a zero-dimensional numerical model is developed in the present study. The model only focuses on the isotope effect for vibrationally excited level of electronically grounded state molecules and its relevant cross-sections. The calculation results show that difference of the ionization channel numbers via molecular vibrationally excited states could be a reason to enhance ionization rate in D plasma.
Our previous study by two dimension in real space and three dimension in velocity space-particle in cell model shows that the curvature of the plasma meniscus causes the beam halo in the negative ion ...sources. The negative ions extracted from the periphery of the meniscus are over-focused in the extractor due to the electrostatic lens effect, and consequently become the beam halo. The purpose of this study is to verify this mechanism with the full 3D model. It is shown that the above mechanism is essentially unchanged even in the 3D model, while the fraction of the beam halo is significantly reduced to 6%. This value reasonably agrees with the experimental result.
The difference of tungsten impurity transport characteristics between a high recycling regime and a partially detached regime has been studied with the IMPGYRO code. Background plasma profiles from a ...JT-60U model geometry, computed from SOLPS, have been used. To obtain such characteristic regimes, we have changed electron (Qe) and ion (Qi) input powers at the core boundary. In the high-recycling regime, the tungsten impurities are transported toward the upstream of the SOL. On the other hand, in the partially detached regime, most tungsten impurities are localized near the inner and outer divertors. These features are mainly related to the background plasma temperature and ion flow.