The purpose of this paper is to understand how low energy plasmaspheric electrons respond to ULF waves excited by interplanetary shocks impinging on magnetosphere. It is found that both energy and pitch angle disperse...The purpose of this paper is to understand how low energy plasmaspheric electrons respond to ULF waves excited by interplanetary shocks impinging on magnetosphere. It is found that both energy and pitch angle dispersed plasmaspheric electrons with energy of a few eV to tens of eV can be generated simultaneously by the interplanetary shock. The subsequent period of successive dispersion signatures is around 40 s and is consistent with the ULF wave period(third harmonic). By tracing back the energy and pitch angle dispersion signatures, the position of the electron injection region is found to be off-equator at around -32° in the southern hemisphere. This can be explained as the result of injected electrons being accelerated by higher harmonic ULF waves(e.g. third harmonic) which carry a larger amplitude electric field off-equator. The dispersion signatures are due to the flux modulations(or accelerations) of " local" plasmaspheric electrons rather than electrons from the ionosphere. With the observed wave-borne large electric field excited by the interplanetary shock impact, the kinetic energy can increase to a maximum of 23 percent in one bouncing cycle for plasmaspheric electrons satisfying the drift-bounce resonance condition by taking account of both the corotating drift and bounce motion of the local plasmaspheric electron.展开更多
Nonlinear interaction between tearing modes(TM) and slab ion-temperature-gradient(ITG) modes is numerically investigated by using a Landau fluid model. It is observed that the energy spectra with respect to wavenumber...Nonlinear interaction between tearing modes(TM) and slab ion-temperature-gradient(ITG) modes is numerically investigated by using a Landau fluid model. It is observed that the energy spectra with respect to wavenumbers become broader during the transition phase from the ITG-dominated stage to TM-dominated stage. Accompanied with the fast growth of the magnetic island, the frequency of TM/ITG with long/short wavelength fluctuations in the electron/ion diamagnetic direction decreases/increases respectively. The decrease of TM frequency is identified to result from the effect of the profile flattening in the vicinity of the magnetic island, while the increase of the frequencies of ITG fluctuations is due to the eigenmode transition of ITG induced by the large scale zonal flow and zonal current related to TM. Roles of zonal current induced by the ITG fluctuations in the instability of TM are also analyzed. Finally, the electromagnetic transport features in the vicinity of the magnetic island are discussed.展开更多
Cyrokinetics is widely applied in plasma physics. However, this framework is limited to weak turbulence levels and low drift-wave frequencies because high-frequency gyro-motion is reduced by the gyro-phase averaging. ...Cyrokinetics is widely applied in plasma physics. However, this framework is limited to weak turbulence levels and low drift-wave frequencies because high-frequency gyro-motion is reduced by the gyro-phase averaging. In order to test where gyrokinetics breaks down, Waltz and Zhao developed a new theory, called cyclokinetics [R. E. Waltz and Zhao Deng, Phys. Plasmas 20, 012507 (2013)]. Cyclokinetics dynamically follows the high-frequency ion gyro-motion which is nonlineaxly coupled to the low-frequency drift-waves interrupting and suppressing gyro-averaging. Cyclokinetics is valid in the high-frequency (ion cyclotron frequency) regime or for high turbulence levels. The ratio of the cyclokinetic perturbed distribution function over equilibrium distribution function δf/F can approach 1. This work presents, for the first time, a numerical simulation of nonlinear cyclokinetic theory for ions, and describes the first attempt to completely solve the ion gyro-phase motion in a nonlinear turbulence system. Simulations are performed [Zhao Deng and R. E. Waltz, Phys. Plasmas 22(5), 056101 (2015)] in a local flux-tube geometry with the parallel motion and variation suppressed by using a newly developed code named rCYCLO, which is executed in parallel by using an implicit time-advanced Eulerian (or continuum) scheme [Zhao Deng and R. E. Waltz, Comp. Phys. Comm. 195, 23 (2015)]. A novel numerical treatment of the magnetic moment velocity space derivative operator guarantee saccurate conservation of incremental entropy. By comparing the more fundamental cyclokinetic simulations with the corresponding gyrokinetic simulations, the gyrokinetics breakdown condition is quantitatively tested. Gyrokinetic transport and turbulence level recover those of cyclokinetics at high relative ion cyclotron frequencies and low turbulence levels, as required. Cyclokinetic transport and turbulence level are found to be lower than those of gyrokinetics at high turbulence levels and low-Ω* values with stable ion cyclotron modes. The gyrokinetic approximation is found to break down when the density perturbation exceeds 20%, or when the ratio of nonlinear E x B frequency over ion cyclotron frequency exceeds 20%. This result indicates that the density perturbation of the Tokamak L-mode near-edge is not sufficiently large for breaking the gyro-phase averaging. For cyclokinetic simulations with sufficiently unstable ion cyclotron (IC) modes and sufficiently low Ω^* ~10, the high-frequency component of the cyclokinetic transport can exceed that of the gyrokinetic transport. However, the low-frequency component of the cyclokinetic transport does not exceed that of the gyrokinetic transport. For higher and more physically relevant Ω^* ≥50 values and physically realistic IC driving rates, the low-frequency component of the cyclokinetic transport remains smaller than that of the gyrokinetic transport. In conclusion, the "L-mode near-edge short-fall" phenomenon, observed in some low-frequency gyrokinetic turbulence transport simulations, does not arise owing to the nonlinear coupling of high-frequency ion cyclotron motion to low-frequency drift motion.展开更多
基金supported by National Natural Science Foundation of China National Natural Science Foundation of China (41421003 and 41627805)
文摘The purpose of this paper is to understand how low energy plasmaspheric electrons respond to ULF waves excited by interplanetary shocks impinging on magnetosphere. It is found that both energy and pitch angle dispersed plasmaspheric electrons with energy of a few eV to tens of eV can be generated simultaneously by the interplanetary shock. The subsequent period of successive dispersion signatures is around 40 s and is consistent with the ULF wave period(third harmonic). By tracing back the energy and pitch angle dispersion signatures, the position of the electron injection region is found to be off-equator at around -32° in the southern hemisphere. This can be explained as the result of injected electrons being accelerated by higher harmonic ULF waves(e.g. third harmonic) which carry a larger amplitude electric field off-equator. The dispersion signatures are due to the flux modulations(or accelerations) of " local" plasmaspheric electrons rather than electrons from the ionosphere. With the observed wave-borne large electric field excited by the interplanetary shock impact, the kinetic energy can increase to a maximum of 23 percent in one bouncing cycle for plasmaspheric electrons satisfying the drift-bounce resonance condition by taking account of both the corotating drift and bounce motion of the local plasmaspheric electron.
基金Project supported by the National Key R&D Program of China(Grant Nos.2017YFE0301100 and 2017YFE0300500)the National Natural Science Foundation of China(Grant Nos.11675038,11775069,and 11305027)the Fundamental Research Funds for the Central Universities of China(Grant No.DUT17RC(4)54)
文摘Nonlinear interaction between tearing modes(TM) and slab ion-temperature-gradient(ITG) modes is numerically investigated by using a Landau fluid model. It is observed that the energy spectra with respect to wavenumbers become broader during the transition phase from the ITG-dominated stage to TM-dominated stage. Accompanied with the fast growth of the magnetic island, the frequency of TM/ITG with long/short wavelength fluctuations in the electron/ion diamagnetic direction decreases/increases respectively. The decrease of TM frequency is identified to result from the effect of the profile flattening in the vicinity of the magnetic island, while the increase of the frequencies of ITG fluctuations is due to the eigenmode transition of ITG induced by the large scale zonal flow and zonal current related to TM. Roles of zonal current induced by the ITG fluctuations in the instability of TM are also analyzed. Finally, the electromagnetic transport features in the vicinity of the magnetic island are discussed.
文摘Cyrokinetics is widely applied in plasma physics. However, this framework is limited to weak turbulence levels and low drift-wave frequencies because high-frequency gyro-motion is reduced by the gyro-phase averaging. In order to test where gyrokinetics breaks down, Waltz and Zhao developed a new theory, called cyclokinetics [R. E. Waltz and Zhao Deng, Phys. Plasmas 20, 012507 (2013)]. Cyclokinetics dynamically follows the high-frequency ion gyro-motion which is nonlineaxly coupled to the low-frequency drift-waves interrupting and suppressing gyro-averaging. Cyclokinetics is valid in the high-frequency (ion cyclotron frequency) regime or for high turbulence levels. The ratio of the cyclokinetic perturbed distribution function over equilibrium distribution function δf/F can approach 1. This work presents, for the first time, a numerical simulation of nonlinear cyclokinetic theory for ions, and describes the first attempt to completely solve the ion gyro-phase motion in a nonlinear turbulence system. Simulations are performed [Zhao Deng and R. E. Waltz, Phys. Plasmas 22(5), 056101 (2015)] in a local flux-tube geometry with the parallel motion and variation suppressed by using a newly developed code named rCYCLO, which is executed in parallel by using an implicit time-advanced Eulerian (or continuum) scheme [Zhao Deng and R. E. Waltz, Comp. Phys. Comm. 195, 23 (2015)]. A novel numerical treatment of the magnetic moment velocity space derivative operator guarantee saccurate conservation of incremental entropy. By comparing the more fundamental cyclokinetic simulations with the corresponding gyrokinetic simulations, the gyrokinetics breakdown condition is quantitatively tested. Gyrokinetic transport and turbulence level recover those of cyclokinetics at high relative ion cyclotron frequencies and low turbulence levels, as required. Cyclokinetic transport and turbulence level are found to be lower than those of gyrokinetics at high turbulence levels and low-Ω* values with stable ion cyclotron modes. The gyrokinetic approximation is found to break down when the density perturbation exceeds 20%, or when the ratio of nonlinear E x B frequency over ion cyclotron frequency exceeds 20%. This result indicates that the density perturbation of the Tokamak L-mode near-edge is not sufficiently large for breaking the gyro-phase averaging. For cyclokinetic simulations with sufficiently unstable ion cyclotron (IC) modes and sufficiently low Ω^* ~10, the high-frequency component of the cyclokinetic transport can exceed that of the gyrokinetic transport. However, the low-frequency component of the cyclokinetic transport does not exceed that of the gyrokinetic transport. For higher and more physically relevant Ω^* ≥50 values and physically realistic IC driving rates, the low-frequency component of the cyclokinetic transport remains smaller than that of the gyrokinetic transport. In conclusion, the "L-mode near-edge short-fall" phenomenon, observed in some low-frequency gyrokinetic turbulence transport simulations, does not arise owing to the nonlinear coupling of high-frequency ion cyclotron motion to low-frequency drift motion.