A non-twisting flux tube for local gyrokinetic simulations
arXiv:2012.04785 · doi:10.1088/1361-6587/abf8f4
Abstract
Local gyrokinetic simulations use a field-aligned domain that twists due to the magnetic shear of the background magnetic equilibrium. However, if the magnetic shear is strong and/or the domain is long, the twist can become so extreme that it fails to properly resolve the turbulence. In this work, we derive and implement the "non-twisting flux tube," a local simulation domain that remains rectangular at all parallel locations. Convergence and runtime tests indicate that it can calculate the heat flux more efficiently than the conventional flux tube. For one test case, it was 30 times less computationally expensive and we found no case for which it was more expensive. It is most advantageous when the magnetic shear is high and the domain includes at least two regions of turbulent drive (e.g. stellarator simulations, pedestal simulations, tokamak simulations with several poloidal turns). Additionally, it more accurately models the inboard midplane when the magnetic shear is large. Lastly, we show how the non-twisting flux tube can be generalized to allow further optimization and control of the simulation domain.
arXiv comment49 pages, 13 figures
References in corpus (6)
- Multiscale Gyrokinetics for Rotating Tokamak Plasmas: Fluctuations, Transport and Energy Flows
- Critically balanced ion temperature gradient turbulence in fusion plasmas
- Up-down symmetry of the turbulent transport of toroidal angular momentum in tokamaks
- How eigenmode self-interaction affects zonal flows and convergence of tokamak core turbulence with toroidal system size
- The Parallel Boundary Condition for Turbulence Simulations in Low Magnetic Shear Devices
- Eliminating turbulent self-interaction through the parallel boundary condition in local gyrokinetic simulations