Related papers: GX: a GPU-native gyrokinetic turbulence code for t…
Finite element schemes based on discontinuous Galerkin methods possess features amenable to massively parallel computing accelerated with general purpose graphics processing units (GPUs). However, the computational performance of such…
Differential rotation is known to suppress linear instabilities in fusion plasmas. However, even in the absence of growing eigenmodes, subcritical fluctuations that grow transiently can lead to sustained turbulence. Here transient growth of…
We study parallel particle-in-cell (PIC) methods for low-temperature plasmas (LTPs), which discretize kinetic formulations that capture the time evolution of the probability density function of particles as a function of position and…
We propose an optimization approach for determining both hardware and software parameters for the efficient implementation of a (family of) applications called dense stencil computations on programmable GPGPUs. We first introduce a simple,…
With the recent proliferation of heterogeneous, GPU-accelerated supercomputers, high-order computational fluid dynamics (CFD) simulations of complex, turbulent flows are more accessible than ever. To leverage the computing power of these…
Discontinuous Galerkin (DG) methods for the numerical solution of partial differential equations have enjoyed considerable success because they are both flexible and robust: They allow arbitrary unstructured geometries and easy control of…
Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks in stellarators. In this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric…
The Gyrokinetic Toroidal Code at Princeton (GTC-P) is a highly scalable and portable particle-in-cell (PIC) code. It solves the 5D Vlasov-Poisson equation featuring efficient utilization of modern parallel computer architectures at the…
The effect of toroidal rotation on both turbulent and neoclassical transport of tungsten (W) in tokamaks is investigated using the flux-driven, global, nonlinear 5D gyrokinetic code GYSELA. Nonlinear simulations are carried out with…
For large-scale graph analytics on the GPU, the irregularity of data access and control flow, and the complexity of programming GPUs, have presented two significant challenges to developing a programmable high-performance graph library.…
This paper presents an overview of the functionalities and applications of Exasim, an open-source code for generating high-order discontinuous Galerkin codes to numerically solve parametrized partial differential equations (PDEs). The…
A new global two-fluid electromagnetic turbulence code, CENTORI, has been developed for the purpose of studying magnetically-confined fusion plasmas on energy confinement timescales. This code is used to evolve the combined system of…
Implicit methods and GPU parallelization are two distinct yet powerful strategies for accelerating high-order CFD algorithms. However, few studies have successfully integrated both approaches within high-speed flow solvers. The core…
Predictive modeling of stellarator plasmas is crucial for advancing nuclear fusion energy, yet it faces unique computational difficulties. One of the main challenges is accurately simulating the dynamics of specific particle species that…
In the numerical solution of partial differential equations using a method-of-lines approach, the availability of high order spatial discretization schemes motivates the development of sophisticated high order time integration methods. For…
A full-F, isothermal, electromagnetic, gyro-fluid model is used to simulate plasma turbulence in a COMPASS-sized, diverted tokamak. A parameter scan covering three orders of magnitude of plasma resistivity and two values for the ion to…
This paper describes a conceptual framework for understanding kinetic plasma turbulence as a generalized form of energy cascade in phase space. It is emphasized that conversion of turbulent energy into thermodynamic heat is only achievable…
Gaussian processes (GPs) are a widely used regression tool, but the cubic complexity of exact solvers limits their scalability. To address this challenge, we extend the GPRat library by incorporating a fully GPU-resident GP prediction…
The influence of collisions on the growth rate of trapped-electron modes (TEMs) in core plasmas is assessed through both analytical linear gyrokinetics and linear gyrokinetic simulations. Both methods are applied to the magnetic geometry of…
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…