The Parallel Boundary Condition for Turbulence Simulations in Low Magnetic Shear Devices

Abstract

Flux tube simulations of plasma turbulence in stellarators and tokamaks typically employ coordinates which are aligned with the magnetic field lines. Anisotropic turbulent fluctuations can be represented in such field-aligned coordinates very efficiently, but the resulting non-trivial boundary conditions involve all three spatial directions, and must be handled with care. The standard "twist-and-shift" formulation of the boundary conditions [Beer, Cowley, Hammett \textit{Phys. Plasmas} \textbf{2}, 2687 (1995)] was derived assuming axisymmetry and is widely used because it is efficient, as long as the global magnetic shear is not too small. A generalization of this formulation is presented, appropriate for studies of non-axisymmetric, stellarator-symmetric configurations, as well as for axisymmetric configurations with small global shear. The key idea is to replace the "twist" of the standard approach (which accounts only for global shear) with the integrated local shear. This generalization allows one significantly more freedom when choosing the extent of the simulation domain in each direction, without losing the natural efficiency of field-line-following coordinates. It also corrects errors associated with naive application of axisymmetric boundary conditions to non-axisymmetric configurations. Simulations of stellarator turbulence that employ the generalized boundary conditions require much less resolution than simulations that use the (incorrect, axisymmetric) boundary conditions. We also demonstrate the surprising result that (at least in some cases) an easily implemented but manifestly incorrect formulation of the boundary conditions does {\it not} change important predicted quantities, such as the turbulent heat flux.

References in corpus (2)

• Residual zonal flows in tokamaks and stellarators at arbitrary wavelengths
• Optimisation of confinement in a fusion reactor using a nonlinear turbulence model

Cited by in corpus (13)

• $\text{stella}$: a mixed implicit-explicit, delta-f gyrokinetic code for general magnetic field configurations
• Eliminating turbulent self-interaction through the parallel boundary condition in local gyrokinetic simulations
• Prevention of core particle depletion in stellarators by turbulence
• A non-twisting flux tube for local gyrokinetic simulations
• Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators
• Dimits shift, avalanche-like bursts, and Solitary propagating structures in the two-field Flux-Balanced Hasegawa-Wakatani model for plasma edge turbulence
• A Flux-Balanced Fluid Model for Collisional Plasma Edge Turbulence: Numerical Simulations with Different Aspect Ratios
• Ion-temperature-gradient stability near the magnetic axis of quasisymmetric stellarators
• How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning
• A phase-shift-periodic parallel boundary condition for low-magnetic-shear scenarios
• Gyrokinetic simulations in stellarators using different computational domains
• Global linear drift-wave eigenmode structures on flux surfaces in stellarators: ion temperature gradient mode
• Microtearing Turbulence and Its Role in High-Density-Gradient Plasmas in Wendelstein 7-X