Theoretical astrophysics and plasma physics at RPC

Abstract:

In electrostatic simulations of MAST plasma at electron-gyroradius scales, using the local flux-tube gyrokinetic code GS2 with adiabatic ions, we find that the long-time saturated electron heat flux (the level most relevant to energy transport) decreases as the electron collisionality decreases. At early simulation times, the heat flux "quasi-saturates" without any strong dependence on collisionality, and with the turbulence dominated by streamer-like radially elongated structures. However, the zonal fluctuation component continues to grow slowly until much later times, eventually leading to a new saturated state dominated by zonal modes and with the heat flux proportional to the collision rate, in approximate agreement with the experimentally observed collisionality scaling of the energy confinement in MAST. We outline an explanation of this effect based on a model of ETG turbulence dominated by zonal-nonzonal interactions and on an analytically derived scaling of the zonal-mode damping rate with the electron-ion collisionality. Improved energy confinement with decreasing collisionality is favourable towards the performance of future, hotter devices.

Abstract:

It is shown that in low-beta plasmas, such as the solar corona, some instances of the solar wind, the aurora, inner regions of accretion discs, their coronae, and some laboratory plasmas, Alfvenic fluctuations produce no ion heating within the gyrokinetic approximation, i.e., as long as their amplitudes (at the Larmor scale) are small and their frequencies stay below the ion Larmor frequency. Thus, all low-frequency ion heating in such plasmas is due to compressive fluctuations: density perturbations and non-Maxwellian perturbations of the ion distribution function. Because these fluctuations energetically decouple from the Alfvenic ones already in the inertial range, the above conclusion means that the energy partition between ions and electrons in low-beta plasmas is decided at the outer scale, where turbulence is launched, and can in principle be determined from MHD models of the relevant astrophysical systems. Any additional ion heating must come from non-gyrokinetic mechanisms such as cyclotron heating or the stochastic heating owing to distortions of ions' Larmor orbits. An exception to these conclusions occurs in the Hall limit, i.e., when the ratio of the ion to electron temperatures is as low as the ion beta (equivalently, the electron beta is order unity). In this regime, compressive fluctuations (slow waves) couple to Alfvenic ones above the Larmor scale (viz., at the ion inertial or ion sound scale), the Alfvenic and compressive cascades join and then separate again into cascades of fluctuations that linearly resemble kinetic Alfven and (oblique) ion cyclotron waves, with the former heating electrons and the latter ions. The two cascades are shown to decouple, scalings for them are derived, and it is argued physically that the two species will be heated by them at approximately equal rates.

Abstract:

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle $\alpha$ between the wall and the magnetic field $\vec{B}$ is oblique. Here, we consider the fusion-relevant case of a shallow-angle, $\alpha \ll 1$, electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for $\alpha / \sqrt{\tau+1} \gg \sqrt{m_{\text{e}}/m_{\text{i}}}$, where $m_{\text{e}}$ is the electron mass, $m_{\text{i}}$ is the ion mass, $\tau = T_{\text{i}}/ZT_{\text{e}}$, $T_{\text{e}}$ is the electron temperature, $T_{\text{i}}$ is the ion temperature, and $Z$ is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii $\rho_{\text{s}} = \sqrt{m_{\text{i}} \left(ZT_{\text{e}} + T_{\text{i}} \right) } / ZeB$, where $e$ is the proton charge and $B = |\vec{B}|$ is the magnitude of the magnetic field. We study the dependence on $\tau $ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by $\tau$. The kinetic model is shown to be asymptotically equivalent to Chodura's fluid model at small ion temperature, $\tau \ll 1$, for $|\ln \alpha| > 3|\ln \tau | \gg 1$. In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyro-orbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, $\tau \gg 1$, relevant because $T_{\text{i}}$ is measured to be a few times larger than $T_{\text{e}}$ near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by $\sqrt{\alpha}$ or $1/\sqrt{\tau}$, depending on which is largest.

Abstract:

Supermassive black hole binaries (SMBHBs) produced in galaxy mergers are thought to complete their coalescence, below separations of r_GW=10^{-3} (M_BH/10^8 M_sun)^{3/4} pc, as their orbit decays due to the emission of gravitational waves (GWs). It may be possible to identify such GW-driven inspirals statistically in an electromagnetic (EM) survey for variable sources. A GW-driven binary spends a characteristic time T_GW at each orbital separation r_orb < r_GW that scales with the corresponding orbital time t_orb as T_GW = (const) t_orb^{8/3}. If the coalescing binary produces variations in the EM emission on this timescale, then it could be identified as a variable source with a characteristic period t_var = t_orb. The incidence rate of sources with similar inferred BH masses, showing near-periodic variability on the time-scale t_var, would then be proportional to t_var^{8/3}. Luminosity variations corresponding to a fraction f_Edd<0.01 of the Eddington luminosity would have been missed in current surveys. However, if the binary inspirals are associated with quasars, we show that a dedicated survey could detect the population of SMBHBs with a range of periods around tens of weeks. The discovery of a population of periodic sources whose abundance obeys N_var = (const) t_var^{8/3} would confirm (i) that the orbital decay is indeed driven by GWs, and (ii) that circumbinary gas is present at small orbital radii and is being perturbed by the BHs. Deviations from the t_var^{8/3} power-law could constrain the structure of the circumbinary gas disk and viscosity-driven orbital decay. We discuss constraints from existing data, and quantify the sensitivity and sky coverage that could yield a detection in future surveys.

Abstract:

KNOSOS (KiNetic Orbit-averaging SOlver for Stellarators) is a freely available, open-source code that calculates neoclassical transport in low-collisionality plasmas of three-dimensional magnetic confinement devices by solving the radially local drift-kinetic and quasineutrality equations. The main feature of KNOSOS is that it relies on orbit-averaging, which removes the dependence on the coordinate along the magnetic field line, and allows to solve the drift-kinetic equation very fast. KNOSOS treats rigorously the effect of the component of the magnetic drift that is tangent to magnetic surfaces, and of the component of the electrostatic potential that varies on the flux-surface, {\varphi}_1. Furthermore, the equation solved is linear in {\varphi}_1, which permits an efficient solution of the quasineutrality equation. As long as the radially local approach is valid, KNOSOS can be applied to the calculation of neoclassical transport in stellarators (helias, heliotrons, heliacs, etc.) and tokamaks with broken axisymmetry. In this paper, we show several calculations for the stellarators W7-X, LHD, NCSX and TJ-II that provide benchmark with standard local codes and demonstrate the advantages of this approach.