Related papers: AstroGK: Astrophysical Gyrokinetics Code
The triangular mesh-based gyrokinetic scheme enables comprehensive axis-to-edge studies across the entire plasma volume. Our approach employs triangular finite elements with first-derivative continuity (C1), building on previous work to…
A comprehensive gyrokinetic simulation model has been implemented in the global toroidal gyrokinetic code (GTC) and verified for studying low-frequency waves and turbulence in magnetic fusion plasmas by treating all kinetic-MHD processes on…
In this work, we have formulated and implemented a mixed unstructured mesh-based finite element (FE)-Fourier decomposition scheme for gyrokinetic simulations in realistic tokamak geometry. An efficient particle positioning…
We present a new numerical code designed to solve the Einstein field equations for axisymmetric spacetimes. The long term goal of this project is to construct a code that will be capable of studying many problems of interest in axisymmetry,…
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…
We present a new model for simulating the electromagnetic fluctuations with frequencies much lower than the ion cyclotron frequency in plasmas confined in general magnetic configurations. This novel model (termed as GK-E&B) employs…
A unified gas kinetic scheme (UGKS) for multi-scale and multi-component plasma transport is constructed. The current scheme is a direct modeling method, where the time evolution solutions from the Vlasov-BGK equations of electron and ion…
In many astrophysical plasmas, the Coulomb collision is insufficient to maintain an isotropic temperature, and the system is driven to the anisotropic regime. In this case, magnetohydrodynamic (MHD) models with anisotropic pressure 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 discontinuous Galerkin (DG) method suitable for large-scale astrophysical simulations on Cartesian meshes as well as arbitrary static and moving Voronoi meshes is presented. Most major astrophysical fluid dynamics codes use a finite…
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…
GX is a code for solving the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for…
This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of…
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…
Solving the Euler equations of ideal hydrodynamics as accurately and efficiently as possible is a key requirement in many astrophysical simulations. It is therefore important to continuously advance the numerical methods implemented in…
We introduce a new code, ECOSMOG, to run N-body simulations for a wide class of modified gravity and dynamical dark energy theories. These theories generally have one or more new dynamical degrees of freedom, the dynamics of which are…
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…
The explicit quasi-monotonic conservative TVD scheme and numerical method for the solution of the gravitational MHD equations are developed. The 2D numerical code for the simulation of multidimensional selfgravitating MHD flows on the…
The properties of the boundary plasma in a tokamak are now recognized to play a key role in determining the achievable fusion power and the lifetimes of plasma-facing components. Accurate quantitative modeling and improved qualitative…
In this paper the global mode structures of linear ion-temperature-gradient (ITG) modes in tokamak plasmas are obtained by combining results from the local gyrokinetic code GS2 with analytical theory. Local gyrokinetic calculations, using…