Related papers: BOUT++: a framework for parallel plasma fluid simu…
BOUT++ is a 3D nonlinear finite-difference plasma simulation code, capable of solving quite general systems of PDEs, but targeted particularly on studies of the edge region of tokamak plasmas. BOUT++ is publicly available, and has been…
The BOUT++ code is used to simulate ELMs in a shifted circle equilibrium. Reduced ideal MHD simulations are first benchmarked against the linear ideal MHD code ELITE, showing good agreement. Diamagnetic drift effects are included finding…
BOUT++ is a software package designed for solving plasma fluid models. It has been used to simulate a wide range of plasma phenomena ranging from linear stability analysis to 3D plasma turbulence, and is capable of simulating a wide range…
A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. This design enables the same code to be…
Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid…
Numerical simulation of plasma turbulence in the Large Plasma Device (LAPD) [Gekelman et al, Rev. Sci. Inst., 62, 2875, 1991] is presented. The model, implemented in the BOUndary Turbulence (BOUT) code [M. Umansky et al, Contrib. Plasma…
Tokamak plasmas are confined by a magnetic field that limits the particle and heat transport perpendicular to the field. Parallel to the field the ionised particles can move freely, so to obtain confinement the field lines are "closed" (ie.…
Challenging aspects of large-scale turbulent edge simulations in plasma physics include robust nonlinear solvers and efficient preconditioners. This paper presents recent advances in the scalable solution of nonlinear partial differential…
The transport of heat and particles in the relatively collisional edge regions of magnetically confined plasmas is a scientifically challenging and technologically important problem. Understanding and predicting this transport requires the…
Quantum computing holds great promise to accelerate scientific computations in fluid dynamics and other classical physical systems. While various quantum algorithms have been proposed for linear flows, developing quantum algorithms for…
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to…
By using a non-local model, fluid simulations can capture kinetic effects in the parallel electron heat-flux better than is possible using flux limiters in the usual diffusive models. Non-local and diffusive models are compared using a test…
The properties of linear instabilities in the Large Plasma Device [W. Gekelman et al., Rev. Sci. Inst., 62, 2875 (1991)] are studied both through analytic calculations and solving numerically a system of linearized collisional plasma fluid…
Self-consistent full-size turbulent-transport simulations of the divertor and SOL of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this…
Three dimensional electromagnetic gyrofluid simulations of the ideal ballooning mode blowout scenario for tokamak edge localized modes (ELMs) are presented. Special emphasis is placed on energetic diagnosis, examining changes in the growth…
A challenging and fundamental research problem is the better understanding and control of the turbulent transport of heat in present-day tokamak fusion experiments. Recent developments in numerical methods along with enormous gains in…
Non-local closures allow kinetic effects on parallel transport to be included in fluid simulations. This is especially important in the scrape-off layer, but to be useful there the non-local model requires consistent kinetic boundary…
Two ocean models are considered for geophysical flow simulations: the multi-layer shallow water equations and the multi-layer primitive equations. For the former, we investigate the parallel performance of exponential time differencing…
We conducted WMLES to examine the performance of a simple and widely used ODE-based equilibrium wall model in a spatially-developing 3D TBL inside a bent square duct (Schwarz and Bradshaw 1994) and 3D separated flows behind a skewed bump…
We review a scalable two- and three-dimensional computer code for low-temperature plasma simulations in multi-material complex geometries. Our approach is based on embedded boundary (EB) finite volume discretizations of the minimal…