Related papers: A conservative scheme for electromagnetic simulati…
A gyrokinetic model is presented that can properly describe strong flows, large and small amplitude electromagnetic fluctuations occurring on scale lengths ranging from the electron Larmor radius to the equilibrium perpendicular pressure…
The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy…
Particle-In-Cell (kinetic) simulations of shear Alfv\'en wave interaction with one dimensional, across the uniform magnetic field, density inhomogeneity (phase mixing) in collisionless plasma were performed for the first time. As a result a…
For more than half a century, most of the plasma scientists have encountered a violation of the conservation laws of charge, momentum, and energy whenever they have numerically solve the first-principle equations of kinetic plasmas, such as…
Gyrokinetic theory is a basis for treating magnetised plasma dynamics slower than particle gyrofrequencies where the scale of the background is larger than relevant gyroradii. The energy of field perturbations can be comparable to the…
Fluid discontinuities, such as shock fronts and vortex sheets, can reflect waves and become unstable to corrugation. Analytical calculations of these phenomena are tractable in the simplest cases only, while their numerical simulations are…
3D2V continuum gyrokinetic simulations of electrostatic plasma turbulence in a straight, open-field-line geometry have been performed using the full-$f$ discontinuous-Galerkin code Gkeyll. These simulations include the basic elements of a…
A surrogate model for particle-in-cell plasma simulations based on a graph neural network is presented. The graph is constructed in such a way as to enable the representation of electromagnetic fields on a fixed spatial grid. The model is…
The robustness and accuracy of marginally resolved discontinuous Galerkin spectral element computations are evaluated for the standard formulation and a kinetic energy conserving split form on complex flow problems of physical and…
This paper describes a new code for simulating astrophysical plasmas that solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal electron fluid (ITEF) [A. Schekochihin et al., Astrophys. J. Suppl. \textbf{182}, 310…
Vortices, turbulence, and unsteady non-laminar flows are likely both prominent and dynamically important features of astrophysical disks. Such strongly nonlinear phenomena are often difficult, however, to simulate accurately, and are…
A variational method is used to derive a self-consistent macro-particle model for relativistic electromagnetic kinetic plasma simulations. Extending earlier work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp. 376-398,…
Nonlinear wave-wave coupling constitutes an important route for the turbulence spectrum evolution in both space and laboratory plasmas. For example, in a reactor relevant fusion plasma, a rich spectrum of symmetry breaking shear Alfv\'en…
We have developed a hybrid code GMEC: Gyro-kinetic Magnetohydrodynamics (MHD) Energetic-particle Code that can numerically simulate energetic particle-driven Alfv\'en eigenmodes and energetic particle transport in tokamak plasmas. In order…
In this paper, a physics-oriented stochastic kinetic scheme will be developed that includes random inputs from both flow and electromagnetic fields via a hybridization of stochastic Galerkin and collocation methods. Based on the BGK-type…
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…
We present a novel method for guaranteeing linear momentum in learned physics simulations. Unlike existing methods, we enforce conservation of momentum with a hard constraint, which we realize via antisymmetrical continuous convolutional…
We present a numerical model and a set of conservative algorithms for Non-Maxwellian plasma kinetics with inelastic collisions. These algorithms self-consistently solve for the time evolution of an isotropic electron energy distribution…
Low density plasmas are characterized by a large scale separation between the gyromotion of particles around local magnetic fields and the macroscopic scales of the system, often making global kinetic simulations computationally…
Split form schemes for Euler and Navier-Stokes equations are useful for computation of turbulent flows due to their better robustness. This is because they satisfy additional conservation properties of the governing equations like kinetic…