Search

We present the first-of-its-kind coupling of a continuum full-f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principles understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and NSTX parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

The effect of varying magnetic shear on scrape-off layer turbulence and profiles is studied via electromagnetic gyrokinetic simulations of a helical scrape-off layer model. We develop a model helical geometry with magnetic shear and a corresponding field-aligned coordinate system, which is used for simulations with the Gkeyll code. We find that perpendicular transport is reduced in cases with stronger shear, resulting in higher peak particle and heat fluxes to the endplates. Electromagnetic effects slightly increase transport in strong shear cases., Comment: 11 pages, 9 figures. arXiv admin note: text overlap with arXiv:2112.06880

We demonstrate that cross-field transport in the scrape-off layer (SOL) can be moderately increased by electromagnetic effects in high-beta regimes, resulting in a broader electron heat-flux width on the endplates. This conclusion is taken from full-$f$ electromagnetic gyrokinetic simulations of a helical SOL model that roughly approximates the SOL of the National Spherical Torus Experiment (NSTX). The simulations have been performed with the Gkeyll code, which recently became the first code to demonstrate the capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines with sheath boundary conditions. We scan the source rate and thus $\beta$ so that the normalized pressure gradient (the MHD ballooning parameter $\alpha \propto \partial \beta / \partial r \propto \beta / L_p$) is scanned over an experimentally-relevant range, $\alpha = 0.3-1.5$. While there is little change in the pressure gradient scale length $L_p$ near the midplane as beta is increased, a 10% increase in cross-field transport near the midplane results in an increase in the electron heat-flux width $\lambda_q$ and a 25% reduction of the peak electron heat flux to the endplates., Comment: 13 pages, 10 figures

We present an energy-conserving discontinuous Galerkin scheme for the full-$f$ electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for $\partial A_\parallel/\partial t$, the inductive component of the parallel electric field, using a generalized Ohm's law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Amp\`ere cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation., Comment: 38 pages, 11 figures

Turbulent dynamics in the scrape-off layer (SOL) of magnetic fusion devices is intermittent with large fluctuations in density and pressure. Therefore, a model is required that allows perturbations of similar or even larger magnitude to the time-averaged background value. The fluid-turbulence code GRILLIX is extended to such a global model, which consistently accounts for large variation in plasma parameters. Derived from the drift reduced Braginskii equations, the new GRILLIX model includes electromagnetic and electron-thermal dynamics, retains global parametric dependencies and the Boussinesq approximation is not applied. The penalisation technique is combined with the flux-coordinate independent (FCI) approach [F. Hariri and M. Ottaviani, Comput.Phys.Commun. 184:2419, (2013); A. Stegmeir et al., Comput.Phys.Commun. 198:139, (2016)], which allows to study realistic diverted geometries with X-point(s) and general boundary contours. We characterise results from turbulence simulations and investigate the effect of geometry by comparing simulations in circular geometry with toroidal limiter against realistic diverted geometry at otherwise comparable parameters. Turbulence is found to be intermittent with relative fluctuation levels of up to 40% showing that a global description is indeed important. At the same time via direct comparison, we find that the Boussinesq approximation has only a small quantitative impact in a turbulent environment. In comparison to circular geometry the fluctuations are reduced in diverted geometry, which is related to a different zonal flow structure. Moreover, the fluctuation level has a more complex spatial distribution in diverted geometry. Due to local magnetic shear, which differs fundamentally in circular and diverted geometry, turbulent structures become strongly distorted in the perpendicular direction and are eventually damped away towards the X-point.

We present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker-Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the `Lenard-Bernstein,' or `Dougherty,' operator. Several novel algorithmic innovations are reported. The concept of weak-equality is introduced and used to define weak-operators to compute primitive moments needed in the updates. Weak-equality is also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integration by parts are used to construct the discrete weak-form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the $v^2$ term is not included in the basis-set used in the discretization. A series of benchmark problems are presented and show that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping. We conclude with an outline of future work, in particular with some indications of how the algorithms presented here can be extended to use the Rosenbluth potentials to compute the drag and diffusion coefficients., Comment: Source code and input files available on bitbucket. See appendix