Related papers: Development of a gyrokinetic-MHD energetic particl…
We present a high-order spatial discretization of a continuum gyrokinetic Vlasov model in axisymmetric tokamak edge plasma geometries. Such models describe the phase space advection of plasma species distribution functions in the absence of…
We describe a Godunov-type magnetohydrodynamic (MHD) code based on the Miyoshi and Kusano (2005) solver which can be used to solve various astrophysical hydrodynamic and MHD problems. The energy equation is in the form of entropy…
In the merging-compression method of plasma start-up, two flux-ropes with parallel toroidal current are formed around in-vessel poloidal field coils, before merging to form a spherical tokamak plasma. This start-up method, used in the…
We describe a novel splitting approach to numerical relativistic magnetohydrodynamics (RMHD) designed to expand its applicability to the domain of ultra-high magnetisation (high-$\sigma$). In this approach, the electromagnetic field is…
We present a new magnetohydrodynamic-particle-in-cell (MHD-PIC) code integrated into the Athena++ framework. It treats energetic particles as in conventional PIC codes while the rest of thermal plasmas are treated as background fluid…
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…
Various MHD (magnetohydrodynamic) equilibrium tools, some of which being recently developed or considerably updated, are used on the COMPASS tokamak at IPP Prague. MHD equilibrium is a fundamental property of the tokamak plasma, whose…
In preparation for extending the JOREK nonlinear MHD code to stellarators, a hierarchy of stellarator-capable reduced and full MHD models has been derived and tested. The derivation was presented at the EFTC 2019 conference. Continuing this…
It is shown that a considerable improvement in the global gyrokinetic electromagnetic simulations can be achieved by a slight modification of the simulation scheme. The new scheme is verified, simulating a Toroidal Alfv\'en Eigenmode in…
A new class of static magnetohydrodynamic (MHD) magnetic island bifurcations is identified in rotating spherical tokamak plasmas during single- and two-fluid resistive MHD simulations. As the magnitude of an externally applied…
The design of commercially feasible magnetic confinement fusion reactors strongly relies on the reduced turbulent transport in the plasma edge during operation in the high confinement mode (H-mode). We present first global turbulence…
We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which…
The equilibrium code, VMEC, is used to study external kinks in low $\beta$ tokamaks and $l=2$ stellarators. The applicability of the code when modelling nonlinear MHD effects is explored in an attempt to understand and predict how the…
Understanding plasma instabilities is essential for achieving sustainable fusion energy, with large-scale plasma simulations playing a crucial role in both the design and development of next-generation fusion energy devices and the…
The induced electric field in a tokamak drives a parallel electron current flow. In an inhomogeneous, finite beta plasma, when this electron flow is comparable to the ion thermal speed, the Alfven mode wave solutions of the electromagnetic…
Simulations of decaying magnetohydrodynamic (MHD) turbulence are performed with a fluid and a kinetic code. The initial condition is an ensemble of long-wavelength, counter-propagating, shear-Alfv\'{e}n waves, which interact and rapidly…
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…
A general formulation of the problem of calculating the spectrum of stable and unstable eigenmodes of linearized perturbations about a magnetically confined toroidal plasma is presented. The analysis is based on a new hydromagnetic…
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…
Runaway electrons (REs) generated during disruption events in tokamaks can carry mega-Ampere level currents, potentially causing damage to plasma-facing components. Understanding RE evolution during disruption events is important for…