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.
Equivalence of two different approaches to global δf gyrokinetic simulations