Related papers: GX: a GPU-native gyrokinetic turbulence code for t…
To faithfully simulate ITER and other modern fusion devices, one must resolve electron and ion fluctuation scales in a five-dimensional phase space and time. Simultaneously, one must account for the interaction of this turbulence with the…
The gyrokinetic turbulence code GS2 was used to investigate the effects of plasma beta on linear, collisionless ion temperature gradient (ITG) modes and trapped electron modes (TEM) in National Compact Stellarator Experiment (NCSX)…
Turbulent flows and fluid-structure interactions (FSI) are ubiquitous in scientific and engineering applications, but their accurate and efficient simulation remains a major challenge due to strong nonlinearities, multiscale interactions,…
Waves are all around us--be it in the form of sound, electromagnetic radiation, water waves, or earthquakes. Their study is an important basic tool across engineering and science disciplines. Every wave solver serving the computational…
We present a GPU-accelerated version of a high-order discontinuous Galerkin discretization of the unsteady incompressible Navier-Stokes equations. The equations are discretized in time using a semi-implicit scheme with explicit treatment of…
Rarefied gas effects are of critical importance for the aerodynamic performance of hypersonic vehicles operating at high altitudes. In these scenarios, conventional computational fluid dynamics (CFD) solvers break down as the linear…
In this work, a comparison of the global gyrokinetic codes EUTERPE and GENE-3D in stellarator configurations of LHD and W7-X is carried out. In linear simulations with adiabatic electrons, excellent agreement is found in the mode numbers,…
I present HPRMAT, a high-performance solver library for the linear systems arising in R-matrix coupled-channel scattering calculations in nuclear physics. Designed as a drop-in replacement for the linear algebra routines in existing…
We propose an efficient, accurate and robust IMEX solver for the compressible Navier-Stokes equation describing non-ideal gases with general equations of state. The method, which is based on an $h-$adaptive Discontinuos Galerkin spatial…
High order accurate and explicit time-stable solvers are well suited for hyperbolic wave propagation problems. As a result of the complexities of real geometries, internal interfaces and nonlinear boundary and interface conditions,…
Numerical radiation-hydrodynamics (RHD) for non-relativistic flows is a challenging problem because it encompasses processes acting over a very broad range of timescales, and where the relative importance of these processes often varies by…
Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become…
We present KinetiX, a software toolkit to generate computationally efficient fuel-specific routines for the chemical source term, thermodynamic and mixture-averaged transport properties for use in combustion simulation codes. The C++…
Nuclear fusion plays a pivotal role in the quest for reliable and sustainable energy production. A major roadblock to viable fusion power is understanding plasma turbulence, which significantly impairs plasma confinement, and is vital for…
We present a set of rules to guide the design of GPU algorithms. These rules are grounded on the principle of reducing waste in GPU utility to achieve good speed up. In accordance to these rules, we propose GPU algorithms for 2D…
Over the past decades, multiple gyrokinetic codes have shown to be able to simulate turbulence and associated transport in the core of Tokamak devices. However, their application to the edge and scrape-off layer (SOL) region presents…
For magnetic confinement fusion in tokamak plasmas, some of the limitations to the particle and energy confinement times are caused by turbulence and collisions between particles in toroidal geometry, which determine the "anomalous" and the…
We present herein the results of a linear gyrokinetic analysis of electromagnetic microinstabilites in the conceptual high-$\beta$, reactor-scale, tight-aspect-ratio tokamak STEP (Spherical Tokamak for Energy Production,…
One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are…
The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, kinetic Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy…