Kinetic Simulations of Neoclassical and Anomalous Transport Processes in Helical Systems

Drift kinetic and gyrokinetic theories and simulations are powerful means for quantitative predictions of neoclassical and anomalous transport fluxes in helical systems such as the Large Helical Device (LHD). The δf Monte Carlo particle simulation code, FORTEC-3D, is used to predict radial profiles of the neoclassical particle and heat transport fluxes and the radial electric field in helical systems. The radial electric field profiles in the LHD plasmas are calculated from the ambipolarity condition for the neoclassical particle fluxes obtained by the global simulations using the FORTEC-3D code, in which effects of ion or electron finite orbit widths are included. Gyrokinetic Vlasov simulations using the GKV code verify the theoretical prediction that the neoclassical optimization of helical magnetic configuration enhances the zonal flow generation which leads to the reduction of the turbulent heat diffusivity χi due to the ion temperature gradient (ITG) turbulence. Comparisons between results for the high ion temperature LHD experiment and the gyrokinetic simulations using the GKV-X code show that the χi profile and the poloidal wave number spectrum of the density fluctuation obtained from the simulations are in reasonable agreements with the experimental results. It is predicted theoretically and confirmed by the linear GKV simulations that the E × B rotation due to the background radial electric field Er can enhance the zonal-flow response to a given source. Thus, in helical systems, the turbulent transport is linked to the neoclassical transport through Er which is determined from the ambipolar condition for neoclassical particle fluxes and influences the zonal flow generation leading to reduction of the turbulent transport. In order to investigate the Er effect on the regulation of the turbulent transport by the zonal flow generation, the flux-tube bundle model is proposed as a new method for multiscale gyrokinetic simulations

各種文章の字種比率

国立国語研究所17の書名 : 国立国語研究所研究報告

表記のゆれを測る

国立国語研究

Radially local approximation of the drift kinetic equation

A novel radially local approximation of the drift kinetic equation is presented. The new drift kinetic equation that includes both E×B and tangential magnetic drift terms is written in the conservative form and it has favorable properties for numerical simulation that any additional terms for particle and energy sources are unnecessary for obtaining stationary solutions under the radially local approximation. These solutions satisfy the intrinsic ambipolarity condition for neoclassical particle fluxes in the presence of quasisymmetry of the magnetic field strength. Also, another radially local drift kinetic equation is presented, from which the positive definiteness of entropy production due to neoclassical transport and Onsager symmetry of neoclassical transport coefficients are derived while it sacrifices the ambipolarity condition for neoclassical particle fluxes in axisymmetric and quasi-symmetric systems

The Eulerian variational formulation of the gyrokinetic system in general spatial coordinates

The Eulerian variational formulation of the gyrokinetic system with electrostatic turbulence is presented in general spatial coordinates by extending our previous work [H. Sugama et al., Phys. Plasmas 25, 102506 (2018)]. The invariance of the Lagrangian of the system under an arbitrary spatial coordinate transformation is used to derive the local momentum balance equation satisfied by the gyrocenter distribution functions and the turbulent potential, which are given as solutions of the governing equations. In the symmetric background magnetic field, the derived local momentum balance equation gives rise to the local momentum conservation law in the direction of symmetry. This derivation is in contrast to the conventional method using the spatial translation in which the asymmetric canonical pressure tensor generally enters the momentum balance equation. In the present study, the variation of the Lagrangian density with respect to the metric tensor is taken to directly obtain the symmetric pressure tensor, which includes the effect of turbulence on the momentum transport. In addition, it is shown in this work how the momentum balance is modified when the collision and/or external source terms are added to the gyrokinetic equation. The results obtained here are considered useful for global gyrokinetic simulations investigating both neoclassical and turbulent transport processes even in general non-axisymmetric toroidal systems