A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and nonzonal components of all...A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and nonzonal components of all perturbed quantities are solved together. Using this new conservative scheme, linear simulation of kinetic ballooning mode has been successfully benchmarked with the electromagnetic hybrid model. Simulations of nonlinear interactions between magnetic islands and the ion temperature gradient(ITG) mode in a tokamak show that the islands rotate at the electron diamagnetic drift velocity. The linear ITG structure shifts from the island O-point toward the X-point due to the pressure flattening effect inside the islands, and the nonlinear ITG structure peaks along the magnetic island separatrix because of the increased pressure gradient there.展开更多
基金supported by the China National Magnetic Confinement Fusion Science Program (Grant No. 2018YFE0304100)the US Department of Energy, Office of Science, Office of Advanced Scientific Computing Research and Office of Fusion Energy Sciences, Scientific Discovery through Advanced Computing (Sci DAC) program under Award Number DE-SC0018270 (Sci DAC ISEP Center)+2 种基金the China Scholarship Council (Grant No. 201306010032)resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory (DOE Contract No. DE-AC05-00OR22725)the National Energy Research Scientific Computing Center (DOE Contract No. DE-AC02-05CH11231)
文摘A conservative scheme of kinetic electrons for gyrokinetic simulations in the presence of magnetic islands has been implemented and verified in the gyrokinetic toroidal code, where zonal and nonzonal components of all perturbed quantities are solved together. Using this new conservative scheme, linear simulation of kinetic ballooning mode has been successfully benchmarked with the electromagnetic hybrid model. Simulations of nonlinear interactions between magnetic islands and the ion temperature gradient(ITG) mode in a tokamak show that the islands rotate at the electron diamagnetic drift velocity. The linear ITG structure shifts from the island O-point toward the X-point due to the pressure flattening effect inside the islands, and the nonlinear ITG structure peaks along the magnetic island separatrix because of the increased pressure gradient there.