How to evaluate neoclassical transport coefficients by a single 
δf simulation

The [Formula: see text] method is extended to calculate all possible components of neoclassical transport matrices, including non-diagonal components and interspecies components of any species combinations in an arbitrary magnetic field configuration, by a single simulation. The method is implemented in a global neoclassical code FORTEC-3D and applied to cases of multi-ion species plasmas in a tokamak and a helical configuration, respectively. The results are benchmarked against analytic solutions and local simulation results. It is found that the results obtained with the newly developed method provide reliable solutions that satisfy the Onsager symmetry relation as well as the intrinsic ambipolarity in a tokamak configuration.