Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userKinetic particle simulations in a global toroidal geometry
82
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The gyrokinetic toroidal code (GTC) has been upgraded for global simulations by coupling the core and scrape-off layer (SOL) regions across the separatrix with field-aligned particle-grid interpolations. A fully kinetic particle pusher for high frequency waves (ion cyclotron frequency and beyond) and a guiding center pusher for low frequency waves have been implemented using cylindrical coordinates in a global toroidal geometry. The two integrators correctly capture the particle orbits and agree well with each other, conserving energy and canonical angular momentum. As a verification and application of this new capability, ion guiding center simulations have been carried out to study ion orbit losses at the edge of the DIII-D tokamak for single null magnetic separatrix discharges. The ion loss conditions are examined as a function of the pitch angle for cases without and with a radial electric field. The simulations show good agreement with past theoretical results and with experimentally observed feature in which high energy ions flow out along the ion drift orbits and then hit the divertor plates. A measure of the ion direct orbit loss fraction shows that the loss fraction increases with the ion energy for DIII-D in the initial velocity space. Finally, as a further verification of the capability of the new code, self-consistent simulations of zonal flows in the core region of the DIII-D tokamak were carried out. All DIII-D simulations were perfomed in the absence of turbulence.
rate research
Read More
This paper presents the current state of the global gyrokinetic code ORB5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The ORB5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a ``hybrid model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.
We describe a systematic development of kinetic entropy as a diagnostic in fully kinetic particle-in-cell (PIC) simulations and use it to interpret plasma physics processes in heliospheric, planetary, and astrophysical systems. First, we calculate kinetic entropy in two forms -- the ``combinatorial form related to the logarithm of the number of microstates per macrostate and the ``continuous form related to $f ln f$, where $f$ is the particle distribution function. We discuss the advantages and disadvantages of each and discuss subtleties about implementing them in PIC codes. Using collisionless PIC simulations that are two-dimensional in position space and three-dimensional in velocity space, we verify the implementation of the kinetic entropy diagnostics and discuss how to optimize numerical parameters to ensure accurate results. We show the total kinetic entropy is conserved to three percent in an optimized simulation of anti-parallel magnetic reconnection. Kinetic entropy can be decomposed into a sum of a position space entropy and a velocity space entropy, and we use this to investigate the nature of kinetic entropy transport during collisionless reconnection. We find the velocity space entropy of both electrons and ions increases in time due to plasma heating during magnetic reconnection, while the position space entropy decreases due to plasma compression. This project uses collisionless simulations, so it cannot address physical dissipation mechanisms; nonetheless, the infrastructure developed here should be useful for studies of collisional or weakly collisional heliospheric, planetary, and astrophysical systems. Beyond reconnection, the diagnostic is expected to be applicable to plasma turbulence and collisionless shocks.
Toroidal drift waves with unconventional mode structures and non-ground eigenstates, which differ from typical ballooning structure mode, are found to be important recently by large scale global gyrokinetic simulations and especially become dominant at strong gradient edge plasmas [cf., Xie and Xiao, Phys. Plasmas, 22, 090703 (2015)]. The global stability and mode structures of drift wave in this steep edge density and temperature gradients are examined by both direct numerical solutions of a model two-dimensional eigen equation and analytical theory employing WKB-ballooning approach. Theory agrees with numerical solutions quite well. Our results indicate that (i) non-ground eigenstates and unconventional mode structures generally exist and can be roughly described by two parameters `quantum number $l$ and ballooning angle $vartheta_k$, (ii) local model can overestimate the growth rate largely, say, $>50%$, and (iii) the narrow steep equilibrium profile leads to twisting (triangle-like) radial mode structures. With velocity space integral, semi-local theory predicts that the critical jump gradient of the most unstable ion temperature gradient mode from ground state $l=0$ to non-ground state $l=1$ is $L_T^{-1}Rsim50$. These features can have important consequences to turbulent transport.
Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations.
Understanding the nature of the turbulent fluctuations below the ion gyroradius in solar-wind turbulence is a great challenge. Recent studies have been mostly in favor of kinetic Alfven wave (KAW) type of fluctuations, but other kinds of fluctuations with characteristics typical of magnetosonic, whistler and ion Bernstein modes, could also play a role depending on the plasma parameters. Here we investigate the properties of the sub-proton-scale cascade with high-resolution hybrid-kinetic simulations of freely-decaying turbulence in 3D3V phase space, including electron inertia effects. Two proton plasma beta are explored: the intermediate $beta_p=1$ and low $beta_p=0.2$ regimes, both typically observed in solar wind and corona. The magnetic energy spectum exhibits $k_perp^{-8/3}$ and $k_|^{-7/2}$ power laws at $beta_p=1$, while they are slightly steeper at $beta_p=0.2$. Nevertheless, both regimes develop a spectral anisotropy consistent with $k_|sim k_perp^{2/3}$ at $k_perprho_p>1$, and pronounced small-scale intermittency. In this context, we find that the kinetic-scale cascade is dominated by KAW-like fluctuations at $beta_p=1$, whereas the low-$beta$ case presents a more complex scenario suggesting the simultaneous presence of different types of fluctuations. In both regimes, however, a non-negligible role of ion Bernstein type of fluctuations at the smallest scales seems to emerge.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
