Related papers: Test of some numerical limiters for the conservati…

Stable, accurate, divergence-free simulation of magnetized supersonic turbulence is a severe test of numerical MHD schemes and has been surprisingly difficult to achieve due to the range of flow conditions present. Here we present a new,…

High-resolution simulations of particle-based kinetic plasma models typically require a high number of particles and thus often become computationally intractable. This is exacerbated in multi-query simulations, where the problem depends on…

We present a new Vlasov code for collisionless plasmas in the nonrelativistic regime. A Darwin approximation is used for suppressing electromagnetic vacuum modes. The spatial integration is based on an extension of the flux-conservative…

The hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell…

A new magnetohydrodynamics (MHD) code based on initial value approach, GMEC_I, has been developed for simulating various MHD physics in tokamak plasmas, as the MHD foundation of the gyrokinetic-MHD energetic particle simulation code (GMEC)…

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…

A new version of GBS (Ricci et al. Plasma Phys. Control. Fusion 54, 124047, 2012; Halpern et al. J. Comput. Phys. 315, 388-408, 2016; Paruta et al. Phys. Plasmas 25, 112301, 2018) is described. GBS is a three-dimensional, flux-driven,…

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of…

The transport of heat out of tokamak plasmas by turbulence is the dominant mechanism limiting the performance of fusion reactors. Turbulence can be driven by the ion temperature gradient (ITG) and suppressed by toroidal sheared flows.…

A challenging and fundamental research problem is the better understanding and control of the turbulent transport of heat in present-day tokamak fusion experiments. Recent developments in numerical methods along with enormous gains in…

The simulation of turbulence in the boundary region of a tokamak is crucial for understanding and optimizing the performance of fusion reactors. In this work, the use of low-rank linear algebra techniques is shown to enhance the efficiency…

The goal of this paper is to understand from a mathematical point of view the magnetic confinement of plasmas for fusion. Following Fr\'enod and Sonnendr\"ucker \cite{FS2}, we first use two-scale convergence tools to derive a gyrokinetic…

Tokamak fusion reactors are actively studied as a means of realizing energy production from plasma fusion. However, due to the substantial cost and time required to construct fusion reactors and run physical experiments, numerical…

Kinetic simulations of collisionless plasmas are computationally challenging due to phase space mixing and filamentation, resulting in fine-scale velocity structures. This study compares three methods developed to reduce artifacts related…

Astrophysical plasmas in relativistic spacetimes, such as black hole accretion flows, are often weakly collisional and require kinetic modeling to capture non-local transport and particle acceleration. However, the extreme scale separation…

Strong $E\times B$ plasma flow shear is beneficial for reducing turbulent transport. However, traditional methods of driving flow shear do not scale well to large devices such as future fusion power plants. In this paper, we use a large…

We perform a numerical analysis of a class of randomly perturbed {H}amiltonian systems and {P}oisson systems. For the considered additive noise perturbation of such systems, we show the long time behavior of the energy and quadratic…

In this work, we propose a new numerical method for the Vlasov-Poisson system that is both asymptotically consistent and stable in the quasineutral regime, i.e. when the Debye length is small compared to the characteristic spatial scale of…

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling…

Spherical tokamaks (STs) have many desirable features that make them a suitable choice for fusion power plants. To understand their confinement properties, accurate calculation of turbulent micro-instabilities is necessary for tokamak…