Using the gyrocenter-gauge kinetic theory,an electromagnetic version of the high frequency gyrokinetic numerical algorithm for particle-in-cell simulation has been developed.The new algorithm,being an alternative to a...Using the gyrocenter-gauge kinetic theory,an electromagnetic version of the high frequency gyrokinetic numerical algorithm for particle-in-cell simulation has been developed.The new algorithm,being an alternative to a direct Lorentz-force simulation,offers an efficient way to simulate the dynamics of plasma heating and current drive with radio frequency waves.Gyrokinetic formalism enables separation of gyrocenter and gyrophase motions of a particle in a strong magnetic field.From this point of view,a particlemay be viewed as a combination of a slow gyrocenter and a quickly changing Kruskal ring.In this approach,the nonlinear dynamics of high frequency waves is described by the evolution of Kruskal rings based on first principles physics.The efficiency of the algorithm is due to the fact that the simulation particles are advanced along the slow gyrocenter orbits,while the Kruskal rings capture fast gyrophase physics.Moreover,the gyrokinetic formalism allows separation of the cold response from kinetic effects in the current,which allows one to use much smaller number of particles than what is required by a direct Lorentz-force simulation.Also,the new algorithm offers the possibility to have particle refinement together with mesh refinement,when necessary.To illustrate the new algorithm,a simulation of the electromagnetic low-hybrid wave propagating in inhomogeneous magnetic field is presented.展开更多
Plasmas in modern tokamak experiments contain a significant fraction of impurity ion species in addition to main deuterium background.A new unlike-particle collision operator forδf particle simulation has been develo...Plasmas in modern tokamak experiments contain a significant fraction of impurity ion species in addition to main deuterium background.A new unlike-particle collision operator forδf particle simulation has been developed to study the nonlocal effects of impurities due to finite ion orbits on neoclassical transport in toroidal plasmas.A new algorithm for simulation of cross-collisions between different ion species includes test-particle and conserving field-particle operators.An improved field-particle operator is designed to exactly enforce conservation of number,momentum and energy.展开更多
基金This work was supported by the Multi-Scale Gyrokinetics project as a part of the U.S.DoE ASCR Multiscale Mathematics Research and Education Program.
文摘Using the gyrocenter-gauge kinetic theory,an electromagnetic version of the high frequency gyrokinetic numerical algorithm for particle-in-cell simulation has been developed.The new algorithm,being an alternative to a direct Lorentz-force simulation,offers an efficient way to simulate the dynamics of plasma heating and current drive with radio frequency waves.Gyrokinetic formalism enables separation of gyrocenter and gyrophase motions of a particle in a strong magnetic field.From this point of view,a particlemay be viewed as a combination of a slow gyrocenter and a quickly changing Kruskal ring.In this approach,the nonlinear dynamics of high frequency waves is described by the evolution of Kruskal rings based on first principles physics.The efficiency of the algorithm is due to the fact that the simulation particles are advanced along the slow gyrocenter orbits,while the Kruskal rings capture fast gyrophase physics.Moreover,the gyrokinetic formalism allows separation of the cold response from kinetic effects in the current,which allows one to use much smaller number of particles than what is required by a direct Lorentz-force simulation.Also,the new algorithm offers the possibility to have particle refinement together with mesh refinement,when necessary.To illustrate the new algorithm,a simulation of the electromagnetic low-hybrid wave propagating in inhomogeneous magnetic field is presented.
基金supported by U.S.DOE Contract No.DE-AC02-09CH11466.
文摘Plasmas in modern tokamak experiments contain a significant fraction of impurity ion species in addition to main deuterium background.A new unlike-particle collision operator forδf particle simulation has been developed to study the nonlocal effects of impurities due to finite ion orbits on neoclassical transport in toroidal plasmas.A new algorithm for simulation of cross-collisions between different ion species includes test-particle and conserving field-particle operators.An improved field-particle operator is designed to exactly enforce conservation of number,momentum and energy.