Theoretical astrophysics and plasma physics at RPC

Abstract:

AbstractWe investigate two flow bursts in a series of Earthward bursty bulk flows (BBFs) observed by the Magnetospheric Multiscale spacecraft in Earth's magnetotail at (−24, 7, 4) RE in Geocentric Solar Magnetospheric coordinates. At the leading edges of the BBFs, we observe complex magnetic field structures. In particular, we focus on one BBF which contains large‐amplitude magnetic field fluctuations on the time scale of the proton gyroperiod, and another with a large scale dipolarization. For both events, the magnetic field structures are associated with flux increases of supra‐thermal ions with energies ≳100 keV. We observe that helium ions dominate the ion flux at energies ≳150 keV. We investigate the ion acceleration mechanism and its dependence on the mass and charge state of H+ and He2+ ions. We show that for both events, the ions with gyroradii smaller than the dawn‐dusk scale of the structure are accelerated by the ion bulk flow. For ions with larger gyroradii, the acceleration is likely due to a localized spatially limited electric field for the event with a large‐scale dipolarization. For the event with fluctuating magnetic field, the acceleration of ions with gyroradii comparable with the scale of the magnetic fluctuations can be explained by resonance acceleration.

Abstract:

The transition from strong (fluidlike) to nearly marginal (Floquet-type) regimes of ion-temperature-gradient (ITG) driven turbulence is studied in the stellarator Wendelstein 7-X by means of numerical simulations. Close to marginality, extended (along magnetic field lines) linearly unstable modes are dominant, even in the presence of kinetic electrons, and provide a drive that results in finite turbulent transport. A total suppression of turbulence above the linear stability threshold of the ITG modes, commonly present in tokamaks and known as the “Dimits shift,” is not observed. We show that this is mostly due to the peculiar radial structure of marginal turbulence, which is more localized than in the fluid case and therefore less likely to be stabilized by shearing flows.

Abstract:

A novel approach to global gyrokinetic simulation is implemented in the flux-tube code stella. This is done by using a subsidiary expansion of the gyrokinetic equation in the perpendicular scale length of the turbulence, originally derived by Parra and Barnes [Plasma Phys. Controlled Fusion, 57 054003, 2015], which allows the use of Fourier basis functions while enabling the effect of radial profile variation to be included in a perturbative way. Radial variation of the magnetic geometry is included by utilizing a global extension of the Grad-Shafranov equation and the Miller equilibrium equations which is obtained through Taylor expansion. Radial boundary conditions that employ multiple flux-tube simulations are also developed, serving as a more physically motivated replacement to the conventional Dirichlet radial boundary conditions that are used in global simulation. It is shown that these new boundary conditions eliminate much of the numerical artefacts generated near the radial boundary when expressing a non-periodic function using a spectral basis. We then benchmark the new approach both linearly and non-linearly using a number of standard test cases.