Engineering topics

Rogers, B. N.

Publications and source records attributed to Rogers, B. N..

4 records

DOE OSTI2020

Fluid and gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Corrigendum: Thermal dynamics in the flux-coordinate independent turbulence code GRILLIX

Here, we have discovered that one of the original articles findings, that the actual choice of the dimensionless parameter χ ∥i0 for the ion heat conductivity seemed to play no role for the simulation results, was not physical. Instead, due to a mistake in the handling of input parameters, χ ∥i0 was actually not used in the code – the dimensionless electron heat conductivity χ ∥e0 was used instead.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Multigrid treatment of implicit continuum diffusion

Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Global 3D Braginskii simulations of the tokamak edge region of IWL discharges

A study of plasma turbulence and profile evolution in conditions of low (L-mode) and high (H-mode) confinement at the edge of an axisymmetric, nested circular flux-surface approximation to an inner wall limited (IWL) Alcator C-Mod discharge is presented, using numerical simulations with the global drift-ballooning (GDB) code. Further, GDB solves drift-reduced Braginskii two-fluid equations for electromagnetic low-frequency turbulence in a 3D annulus centered on the last closed flux-surface (LCFS). Three simulations that investigate the conditions of a reference L-mode, a high density, and a high temperature (or H-mode-like) shot were performed using realistic parameters. L-mode transport appears to be largely driven by drift resistive ballooning structures. Its pressure profile exhibits a near-SOL breakpoint that Mirror Langmuir Probes (MLP) detect in C-Mod. The high density simulation sees an increase in the size of convective cells and enhanced turbulent transport, while H-mode conditions develop improved confinement, balanced $E$ x $B$ and ion diamagnetic drifts in the closed-flux region, and spontaneous generation of temperature pedestal with a density pedestal remaining absent. A statistical characterization of the turbulence both in the SOL and the closed-flux region is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗