A fully-implicit Particle-In-Cell Monte Carlo Collision code for the simulation of inductively coupled plasmas

Stefano Mattei, Minh Quang Tran
2017
Journal paper

Abstract

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. (C) 2017 The Authors. Published by Elsevier Inc.

Official source

https://infoscience.epfl.ch/record/231744?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

Linac4 is the new negative hydrogen ion (H-) linear accelerator of the European Organization for Nuclear Research (CERN). Its ion source operates on the principle of Radio-Frequency Inductively Coupled Plasma (RF-ICP) and it is required to provide 50 mA of H- beam in pulses of 600 us with a repetition rate up to 2 Hz and within an RMS emittance of 0.25 pi mm mrad in order to fullfil the requirements of the accelerator. This thesis is dedicated to the characterization of the hydrogen plasma in the Linac4 H- ion source. We have developed a Particle-In-Cell Monte Carlo Collision (PIC-MCC) code to simulate the RF-ICP heating mechanism and performed measurements to benchmark the fraction of the simulation outputs that can be experimentally accessed. The code solves self-consistently the interaction between the electromagnetic field generated by the RF coil and the resulting plasma response, including a kinetic description of charged and neutral species. A fully-implicit implementation allowed to simulate the high density regime of the Linac4 H- ion source, ensuring the energy conservation while maintaining the computational resources tractable. We studied the capacitive to inductive transition characteristic of the initial phase of the pulsed discharge. The simulation results were confirmed by time-resolved photometry measurements and allowed quantifying the effect of the hydrogen pressure and of the external magnetic cusp field on the transition dynamics. This provided insights into possible modifications to the magnetic cusp field configuration to maximize the power deposited to the plasma. The optimal ion source configuration maximizes the density of volume produced H-, the flux of H0 atoms onto the cesiated molybdenum plasma electrode surface at the origin of H- emission, and minimizes the electron density and energy in the beam formation region. We simulated the high-density regime (10^19 m-3) representative of the nominal operation of the Linac4 ion source during beam extraction. We performed a parametric study of the RF current, hydrogen pressure and magnetic configuration (cusp and filter) to assess their impact on the plasma parameters. The simulation results allowed assessment of these parameters and provided guidelines for the optimization of the ion source operational and design parameters. The simulation results of electron density, electron energy and hydrogen dissociation degree showed excellent agreement with optical emission spectroscopy measurements both as a function of RF coil current and magnetic configuration. The outputs of these simulations provide crucial inputs to beam formation and extraction physics models. Dedicated PIC software packages are being developed that will eventually shed insight into essential beam parameters such as the intensity and emittance of the H° beam and the intensity of the co-extracted electrons.

EPFL2017

Numerical Simulation of the RF Plasma Discharge in the Linac4 H- Ion Source

Stefano Mattei, Minh Quang Tran

This paper presents a Particle-In-Cell Monte Carlo Collision simulation of the Radio-Frequency (RF) plasma heating in the Linac4 H- ion source at CERN. The model self-consistently takes into account the electromagnetic field generated by the RF coil, the external static magnetic fields and the resulting plasma response, including a kinetic description of the charged species ((e)-, H+, H-2(+), H-3(+), H-), as well as the atomic and molecular (vibrationally resolved) populations. The simulation is performed for the nominal operational condition of 40 kW RF power and 3 Pa H-2 pressure. Results show that the plasma spatial distribution is non-uniform in the plasma chamber, with a density peak of n(e) = 5 . 10(19) m(-3) in the RF coil region. In the filter field region the electron density drops by two orders of magnitude, with a substantial reduction of the electron energy as well. This results in a ratio e/H- approximate to 1 in the extraction region. The vibrational population is characterized by a two temperature distribution, with the high vibrational states showing a factor 2 higher termperature. A very good agreement is found between the simulation results and optical emission spectroscopy measurement performed on a dedicated test stand at CERN.

Amer Inst Physics2017

Model-based optimization of magnetic control in the TCV tokamak: design and experiments

Federico Pesamosca

Thermonuclear controlled fusion is a promising answer to the current energy and climate issues, providing a safe carbon-free source of energy which is virtually inexhaustible. In magnetic confinement thermonuclear fusion based on tokamak reactors, hydrogen fuel in the state of plasma is confined using a system of external and self-generated magnetic fields. This thesis contributes to the development of magnetic confinement fusion research by applying techniques derived from control engineering to designing magnetic controllers for tokamak plasmas. Specifically, a new design for feedback control of plasma shape and position in the TCV tokamak is provided, and its efficacy is studied in dedicated simulations and experiments.Elongated plasmas lead to improved plasma performance in tokamaks, which is key to sustaining fusion conditions in reactors. The required magnetic field for shaping the plasma column however results in an unstable equilibrium that makes feedback control of the vertical plasma position mandatory. Active stabilization of the axisymmetric plasma vertical instability is a standard feature of elongated tokamaks and will be a fundamental feature in the ITER magnetic control system, since a loss of vertical control and the subsequent plasma disruption can lead to unacceptable heat loads on the plasma facing components.The TCV tokamak is the ideal benchmark for investigating the effect of plasma shaping on tokamak physics and performance, with its system of 16 independently powered poloidal field coils and a vessel with an elongated cross section. Shape and position control are coupled problems in TCV as they share the same poloidal field coils as actuators, requiring a multivariable approach to designing magnetic controllers.In this thesis, controller design for TCV is based on a model for the coupled plasma-vessel-coils electromagnetic dynamics: the RZIp model. In this axisymmetric model, the plasma current distribution is fixed but is free to move radially and vertically in the poloidal plane of a toroidal reference frame. An extension to this model is suggested, relaxing the assumption of rigid displacement in the radial direction to include plasma shape deformation and leading to a semi-rigid RZIp model which better fits numerical equilibria.An improvement to the existing algorithm for shape and position control in TCV is then proposed. In this new approach, tuning of the plasma position controller, in charge of vertical stabilization, can be performed independently of the shape controller, which itself acts on a stable system. Static decoupling is achieved and the shape controller is designed on the basis of an improved model for the plasma deformation, which includes the plasma contribution to the static magnetic flux perturbation. Simulations in closed loop with the RZIp model are provided to evaluate several optimized schemes.Finally, the vertical controller is optimized including the plasma dynamics as part of the controller design. Structured H-infinity, extending classical H-infinity to fixed-structure control systems, is applied to obtain a controller using all available coils for position control, and in particular a coil combination optimized for vertical stabilization. Closed-loop performance improvement is demonstrated in dedicated TCV experiments, confirming the simulation results and paving the way for the routine integration of the optimized position and shape controller in TCV discharges.

EPFL2021

EPFL2017

Model-based optimization of magnetic control in the TCV tokamak: design and experiments

Federico Pesamosca

EPFL2021