Related papers: Some numerical aspects of the conservative PSM sch…
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…
The simulation of problems in kinetic plasma physics are often challenging due to strongly coupled phenomena across multiple scales. In this work, we propose a wavelet-based coarse-grained numerical scheme, based on the framework of…
Verification and validation of the internal kink instability in tokamak have been performed for both gyrokinetic (GTC) and kinetic-MHD codes (GAM-solver, M3D-C1-K, NOVA, XTOR-K). Using realistic magnetic geometry and plasma profiles from…
We present the TRIDK code, which is written in Python, that solves Triad Reduced Interactions for a Drift Kinetic system of equations. The four-dimensional drift kinetic system of equations captures nonlinear electromagnetic (or…
In this work, we focus on the development of the particle-in-cell scheme and the application to the studies of Alfv\'en waves and energetic particle physics in tokamak plasmas. The $\delta f$ and full $f$ schemes are formulated on the same…
We have developed a deterministic conservative solver for the inhomogeneous Fokker-Planck-Landau equation coupled with the Poisson equation, which is a {classical mean-field} primary model for collisional plasmas. Two subproblems, i.e. the…
It is shown that a considerable improvement in the global gyrokinetic electromagnetic simulations can be achieved by a slight modification of the simulation scheme. The new scheme is verified, simulating a Toroidal Alfv\'en Eigenmode in…
Extended magnetohydrodynamic (MHD) simulations of tokamak plasmas regularly produce outputs in multi-dimensional, multiple-field formats; these code-specific formats make it difficult to do cross-code validation/coupling and analyze at a…
Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the turbulent $\boldsymbol{E}\times\boldsymbol{B}$…
A largely unsolved theoretical issue in controlled fusion research is the consistent \textit{kinetic} treatment of slowly-time varying plasma states occurring in collisionless and magnetized axisymmetric plasmas. The phenomenology may…
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 rampdown phase of a tokamak pulse is difficult to simulate and often exacerbates multiple plasma instabilities. To reduce the risk of disrupting operations, we leverage advances in Scientific Machine Learning (SciML) to combine physics…
Simulation of plasmas in electromagnetic fields requires numerical solution of a kinetic equation that describes the time evolution of the particle distribution function. In this paper we propose a finite volume scheme based on integral…
Classically simulating circuits with bosonic codes is challenging due to the prohibitive cost of simulating quantum systems with many, possibly infinite, energy levels. We propose an algorithm to simulate circuits with encoded…
A conservative formulation of the drift-reduced fluid plasma model is constructed by analytically inverting the implicit relation defining the polarisation velocity as a function of the time-derivative of the electric field. The obtained…
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 scheme that tightly couples kinetic turbulence codes across a spatial interface is introduced. This scheme evolves from considerations of competing strategies and down-selection. It is found that the use of a composite kinetic…
Continuum computational kinetic plasma models evolve the distribution function of a plasma species $f_s$ on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence,…
We develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is…
Two-dimensional gyrokinetics is a simple paradigm for the study of kinetic magnetised plasma turbulence. In this paper, we present a comprehensive theoretical framework for this turbulence. We study both the inverse and direct cascades (the…