Related papers: Test of some numerical limiters for the conservati…
A control oriented, lumped parameter model for the tokamak transformer including the slow flux penetration in the plasma (skin effect transformer model) is presented. The model does not require detailed or explicit information about plasma…
Full-f global gyrokinetic simulations of the plasma boundary have until now required heroic computational efforts and case-by-case expert intervention, precluding systematic parameter scans. Here we demonstrate a paradigm shift: hundreds of…
The linear drift-tearing mode is analyzed for different regimes of the plasma-$\beta$, ion-skin-depth parameter space with an unreduced, extended-MHD model. New dispersion relations are found at moderate plasma $\beta$ and previous…
Performing a stable, long duration simulation of driven MHD turbulence with a high thermal Mach number and a strong initial magnetic field is a challenge to high-order Godunov ideal MHD schemes because of the difficulty in guaranteeing…
High-temperature superconducting (HTS) magnets and other advances have led to renewed interest in magnetic mirrors for fusion energy. The non-Maxwellian nature of mirror plasmas necessitates kinetic modeling to predict, optimize and design…
We set up a mapping procedure able to translate the evolution of the radial profile of fast ions, interacting with Toroidal Alfv\'en Eigenmodes, into the dynamics of an equivalent one dimensional bump-on-tail system. We apply this mapping…
Second-order Volume-preserving algorithms (VPAs) for simulating charged particle motion in electromagnetic fields have been generalized to a rotating angle formulation by using the matrix decomposition methods. Based on this method, the…
We analyze why the discretization of linear transport with asymmetric Hermite basis functions can be instable in quadratic norm. The main reason is that the finite truncation of the infinite moment linear system looses the skew-symmetry…
A class of variational schemes for the hydrodynamic-electrodynamic model of lossless free-electron gas in a quasineutral background is developed for high-quality simulations of surface plasmon polaritons. The Lagrangian density of lossless…
This work reports on the development and numerical implementation of the linear electromagnetic gyrokinetic (GK) model in a tokamak flux-tube geometry using a moment approach based on the expansion of the perturbed distribution function on…
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to…
Plasma-jet-driven magneto-inertial fusion (PJMIF) is an alternative approach to controlled nuclear fusion which aims to utilize a line-replaceable dense plasma liner as a repetitive spherical compression driver. In this experiment, first…
Simulation of many-particle system evolution by molecular dynamics takes to decrease integration step to provide numerical scheme stability on the sufficiently large time interval. It leads to a significant increase of the volume of…
Understanding turbulent transport physics in the tokamak edge and scrape-off layer (SOL) is critical to developing a successful fusion reactor. The dynamics in these regions plays a key role in achieving high fusion performance by…
Abstract Simulation domain in field-aligned coordinates of the electrostatic gyrokinetic nonlinear turbulence global code, NLT, is extended to include the magnetic axis. The artificial boundary near the magnetic axis is replaced by the…
In this work, we present first-of-their-kind nonlinear local gyrokinetic simulations of electromagnetic turbulence at mid-radius in the burning plasma phase of the conceptual high-$\beta$, reactor-scale, tight-aspect-ratio tokamak STEP…
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO…
In the present article we describe a few simple and efficient finite volume type schemes on moving grids in one spatial dimension combined with appropriate predictor-corrector method to achieve higher resolution. The underlying finite…
Simulating fully kinetic, two-species plasmas is computationally challenging due to the stiff multiscale dynamics of electrons and ions. While enforcing a quasi-neutral time evolution mitigates this stiffness, it requires an electric…
Designing economical magnetic confinement fusion power plants motivates computational tools that can estimate plasma behavior from engineering parameters without direct reliance on experimental measurement of the plasma profiles. In this…