Related papers: Electromagnetic full-$f$ gyrokinetics in the tokam…
A conservative discontinuous Galerkin scheme for a nonlinear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it…
We are developing a new continuum gyrokinetic code, Gkeyll, for use in edge plasma simulations, and here present initial simulations of turbulence on open field lines with model sheath boundary conditions. The code implements an energy…
Boundary plasma physics plays an important role in tokamak confinement, but is difficult to simulate in a gyrokinetic code due to the scale-inseparable nonlocal multi-physics in magnetic separatrix and open magnetic field geometry. Neutral…
A conservative scheme has been formulated and verified for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from drift kinetic equation is used to…
The kinetic effects of electrons are important to long wavelength magnetohydrodynamic(MHD)instabilities and short wavelength drift-Alfvenic instabilities responsible for turbulence transport in magnetized plasmas, since the non-adiabatic…
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…
Gyrokinetic simulations have greatly improved our theoretical understanding of turbulent transport in fusion devices. Most gyrokinetic models in use are delta-f simulations in which the slowly varying radial profiles of density and…
The gyrokinetic particle simulation serves as a powerful tool for the studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector…
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…
A high-order piecewise field-aligned triangular finite element method is developed and implemented for global electromagnetic gyrokinetic particle-in-cell simulations of tokamak plasmas with open field lines. The approach combines locally…
A framework for performing dynamic mesh adaptation with the discontinuous Galerkin method (DGM) is presented. Adaptations include modifications of the local mesh step size (h-adaptation) and the local degree of the approximating polynomials…
In this work, the non-ideal Ohm's law is integrated in the mixed variable-pullback scheme for the gyrokinetic particle simulations. This scheme captures the evolution of the symplectic solution of the gyrokinetic model accurately not only…
This paper generalizes the earlier work on the energy-based discontinuous Galerkin method for second-order wave equations to fourth-order semilinear wave equations. We first rewrite the problem into a system with a second-order spatial…
Geodesic acoustic modes (GAMs) are studied by means of the gyrokinetic global particle-in-cell code ORB5. Linear electromagnetic simulations in the low electron beta limit have been performed, in order to separate acoustic and Alfv\'enic…
Relativistic plasmas around compact objects can sometimes be approximated as being force-free. In this limit, the plasma inertia is negligible and the overall dynamics is governed by global electric currents. We present a novel numerical…
In this paper we use the GeneralizedMultiscale Finite ElementMethod (GMsFEM) framework, introduced in [20], in order to solve nonlinear elliptic equations with high-contrast coefficients. The proposed solution method involves linearizing…
We generalize the energy-based discontinuous Galerkin method proposed in [SIAM J. Num. Anal., 53(6):2705-2726, 2015.] to second-order semilinear wave equations. A stability and convergence analysis is presented along with numerical…
A newly developed weak Galerkin method is proposed to solve parabolic equations. This method allows the usage of totally discontinuous functions in approximation space and preserves the energy conservation law. Both continuous and…
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…
Recently, the electrostatic gyrokinetic Hamiltonian and change of coordinates have been computed to order $\epsilon^2$ in general magnetic geometry. Here $\epsilon$ is the gyrokinetic expansion parameter, the gyroradius over the macroscopic…