Abstract

A self-consistent model is presented for the simulation of a multi-component plasma in the tokamak boundary. A deuterium plasma is considered, with the plasma species that include electrons, deuterium atomic ions and deuterium molecular ions, while the deuterium atoms and molecules constitute the neutral species. The plasma and neutral models are coupled via a number of collisional interactions, which include dissociation, ionization, charge-exchange and recombination processes. The derivation of the three-fluid drift-reduced Braginskii equations used to describe the turbulent plasma dynamics is presented, including its boundary conditions. The kinetic advection equations for the neutral species are also derived, and their numerical implementation discussed. The first results of multi-component plasma simulations carried out by using the global Braginskii solver (GBS) code are then presented and analyzed, being compared with results obtained with the single-component plasma model.

Highlights

  • The boundary of a tokamak plays a crucial role in determining the overall performance of the device, as it sets the confinement of particles and heat, determines the heat exhaust to the vessel walls and controls the impurity level in the core [1]

  • We describe the development and numerical implementation in the GBS code of a multi-component model that addresses the turbulent multi-ion species plasma dynamics through a set of fluid drift-reduced Braginskii equations, while each multiple neutral species are simulated by solving a kinetic equation

  • We present the first results from simulations of turbulence in the tokamak boundary carried out by using the multi-component plasma model described in Secs

Read more

Summary

Introduction

The boundary of a tokamak plays a crucial role in determining the overall performance of the device, as it sets the confinement of particles and heat, determines the heat exhaust to the vessel walls and controls the impurity level in the core [1]. The resulting nHESEL [33, 26, 34] code allows for the simulation of a single-ion plasma including the interactions with three neutral species: cold hydrogen molecules puffed into the system, warm atoms resulting from the dissociation of the hydrogen molecules and hot hydrogen atoms generated by charge-exchange processes Such a model was used to study the plasma fueling in the presence of gas puffs and the formation of a density shoulder in the tokamak boundary at a high gas puffing rate. We describe the development and numerical implementation in the GBS code of a multi-component model that addresses the turbulent multi-ion species plasma dynamics through a set of fluid drift-reduced Braginskii equations, while each multiple neutral species are simulated by solving a kinetic equation. D presents the neutral system of equations in the matrix form implemented in GBS

Collisional processes in multi-component deuterium plasmas
The three-fluid drift-reduced Braginskii equations
Boundary conditions
The kinetic model for the neutral species and its formal solution
First simulation of a multi-component plasma with the GBS code
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call