Publication

Global Flux-Driven Simulations of Ion Temperature-Gradient and Trapped-Electron Modes Driven Turbulence with an Improved Multithreaded Gyrokinetic PIC Code

Emmanuel Lanti
2019
EPFL thesis
Abstract

Performance of tokamak fusion plasmas is heavily linked to the radial heat and particle transport, which is known to be mainly produced by turbulence driven by micro-instabilities. Understanding such processes is thus of key importance for the design and operation of future fusion reactors such as ITER. One of the most commonly used tool to study turbulent transport in the very hot core of tokamaks is gyrokinetic theory, which aims at reducing the physical complexity by averaging out the fast gyromotion, while keeping its averaged effect. The \orb{} program is a numerical code solving the gyrokinetic Maxwell-Boltzmann system of equations to study various phenomena such as turbulent transport. It is a global particle-in-cell code using a finite element representation with B-spline basis function for the fields and various noise control schemes.

The main objective of this work is to improve the performance and physical model of \orb{} in order to be able to carry out simulations relevant for the Tokamak à Configuration Variable (TCV) at the Swiss Plasma Center. This is done in three major steps.

First, the code undergoes a complete refactoring and optimization process with the aim of increasing its maintability and performance. Using a simple test bed being a simplification of the PIC method while retaining the key elements, various algorithm and parallelization schemes are developed. Furthermore, it is shown that one of the most critical kernels of the PIC algorithm, the charge deposition, can be accelerated by a factor 2.4 in \orb{}.

Second, the hybrid electron model, which is used to simulate certain electron modes at a lower numerical cost than the fully kinetic electron model, needs to be corrected. Indeed, in this model, trapped electrons are kinetic while passing electrons have an adiabatic response. This is problematic because it does not satisfy the ambipolarity condition and produces spurious sources of particles and momentum. A corrected model in which the flux-surface averaged contribution of the passing electrons is treated as kinetic is developed and implemented in \orb{}. As a first step, an approximation of the flux-surface average consisting of a Fourier filter retaining only the n=m=0n=m=0 mode is tested. Linearly, both the proper flux-surface average and its n=m=0n=m=0 approximation lead to an accurate description of the geodesic acoustic mode dynamics. Nonlinearly, the n=m=0n=m=0 model exhibit interesting physical features such as sharp radial structures at the mode rational surface q=1q=1 and a transport regime transition.

Finally, a moment-conserving heat source is implemented in the \orb{} code to allow for flux-driven simulations. The heat source conservation properties are verified and it is shown that conservation of density, parallel momentum, and zonal flows is satisfied up to machine precision. A mode-switching'' approach consisting in using an effective heat source computed from a gradient-driven simulation as an input for the flux-driven run is also tested and consistent results are obtained between gradient-driven, flux-driven, and mode-switched'' simulations. Finally, the heat source is characterized, along with variations, using a 2D phase-space diagnostics. It is shown that in some cases, the use of this fixed heat source can lead to a sampling problem causing an artificial source/sink of moments.

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 (36)
Tokamak
A tokamak (ˈtoʊkəmæk; токамáк) is a device which uses a powerful magnetic field to confine plasma in the shape of a torus. The tokamak is one of several types of magnetic confinement devices being developed to produce controlled thermonuclear fusion power. , it was the leading candidate for a practical fusion reactor. Tokamaks were initially conceptualized in the 1950s by Soviet physicists Igor Tamm and Andrei Sakharov, inspired by a letter by Oleg Lavrentiev. The first working tokamak was attributed to the work of Natan Yavlinsky on the T-1 in 1958.
Magnetic confinement fusion
Magnetic confinement fusion is an approach to generate thermonuclear fusion power that uses magnetic fields to confine fusion fuel in the form of a plasma. Magnetic confinement is one of two major branches of fusion energy research, along with inertial confinement fusion. The magnetic approach began in the 1940s and absorbed the majority of subsequent development. Fusion reactions combine light atomic nuclei such as hydrogen to form heavier ones such as helium, producing energy.
Simulation
A simulation is the imitation of the operation of a real-world process or system over time. Simulations require the use of models; the model represents the key characteristics or behaviors of the selected system or process, whereas the simulation represents the evolution of the model over time. Often, computers are used to execute the simulation. Simulation is used in many contexts, such as simulation of technology for performance tuning or optimizing, safety engineering, testing, training, education, and video games.
Show more
Related publications (233)

Inter-discharge optimization for fast, reliable access to ASDEX Upgrade advanced tokamak scenario

Olivier Sauter, Federico Alberto Alfredo Felici, Cassandre Ekta Contré, Simon Van Mulders, Hartmut Zohm

Rapid inter-discharge simulation and optimization using the RAPTOR code have allowed the development of a reliable and reproducible early heating strategy for an advanced tokamak (AT) scenario on ASDEX Upgrade. Solving for electron heat and current diffusi ...
Bristol2024

Multi-machine benchmark of the self-consistent 1D scrape-off layer model DIV1D from stagnation point to target with SOLPS-ITER

Holger Reimerdes

This paper extends a 1D dynamic physics-based model of the scrape-off layer (SOL) plasma, DIV1D, to include the core SOL and possibly a second target. The extended model is benchmarked on 1D mapped SOLPS-ITER simulations to find input settings for DIV1D th ...
Bristol2024

First-principle based predictions of the effects of negative triangularity on DTT scenarios

Olivier Sauter, Gabriele Merlo, Alberto Mariani

Plasmas with negative triangularity (NT) shape have been recently shown to be able to achieve H-mode levels of confinement in L-mode, avoiding detrimental edge localised modes. Therefore, this plasma geometry is now studied as a possible viable option for ...
Bristol2024
Show more
Related MOOCs (32)
Plasma Physics: Introduction
Learn the basics of plasma, one of the fundamental states of matter, and the different types of models used to describe it, including fluid and kinetic.
Plasma Physics: Introduction
Learn the basics of plasma, one of the fundamental states of matter, and the different types of models used to describe it, including fluid and kinetic.
Plasma Physics: Applications
Learn about plasma applications from nuclear fusion powering the sun, to making integrated circuits, to generating electricity.
Show more