Theoretical astrophysics and plasma physics at RPC

Abstract:

We comprehensively model the X-ray luminosity emergent from time-dependent relativistic accretion discs, developing analytical models of the X-ray luminosity of thermal disc systems as a function of black hole mass M, disc mass Md, and disc α-parameter. The X-ray properties of these solutions will be directly relevant for understanding tidal disruption event (TDE) observations. We demonstrate an extremely strong suppression of thermal X-ray luminosity from large mass black holes, LX ∼ exp (− m7/6), where m is a dimensionless mass, roughly the black hole mass in unity of 106M⊙. This strong suppression results in upper observable black hole mass limits, which we demonstrate to be of order Mlim ≃ 3 × 107M⊙, above which thermal X-ray emission will not be observable. This upper observable black hole mass limit is a function of the remaining disc parameters, and the full dependence can be described analytically (equation 82). We demonstrate that the current population of observed X-ray TDEs is indeed consistent with an upper black hole mass limit of order M ∼ 107M⊙, consistent with our analysis.

Abstract:

A new pseudolocal tomography algorithm is developed for soft X-ray(SXR) imaging measurements of the turbulent electron temperature fluctuations (δ Te) in tokamaks and stellarators. The algorithm overcomes the constraints of limited viewing ports on the vessel wall (viewing angle) and limited number of lines of sight (LOS). This is accomplished by increasing the number of LOS locally in a region of interest. Numerical modeling demonstrates that the wavenumber spectrum of the turbulence can be reliably reconstructed, with an acceptable number of viewing angles and LOS and suitable low SNR detectors. We conclude that a SXR imaging diagnostic for measurements of turbulent δ Te using a pseudolocal reconstruction algorithm is feasible

Abstract:

A new algorithm for toroidal flow shear in a linearly implicit, local δf gyrokinetic code is described. Unlike the current approach followed by a number of codes, it treats flow shear continuously in time. In the linear gyrokinetic equation, time-dependences arising from the presence of flow shear are decomposed in such a way that they can be treated explicitly in time with no stringent constraint on the time step. Flow shear related time dependences in the nonlinear term are taken into account exactly, and time dependences in the quasineutrality equation are interpolated. Test cases validating the continuous-in-time implementation in the code GS2 are presented. Lastly, nonlinear gyrokinetic simulations of a JET discharge illustrate the differences observed in turbulent transport compared with the usual, discrete-in-time approach. The continuous-in-time approach is shown, in some cases, to produce fluxes that converge to a different value than with the discrete approach. The new approach can also lead to substantial computational savings by requiring radially narrower boxes. At fixed box size, the continuous implementation is only modestly slower than the previous, discrete approach.

Abstract:

We have used the local-δf gyrokinetic code GS2 to perform studies of the effect of flux-surface shaping on two highly-shaped, low- and high-β JT-60SA-relevant equilibria, including a successful benchmark with the GKV code. We find that for a high-performance plasma, i.e. one with high plasma beta and steep pressure gradients, the turbulent outwards radial fluxes may be reduced by minimizing the elongation. We explain the results as a competition between the local magnetic shear and finite-Larmor-radius (FLR) stabilization. Electromagnetic studies indicate that kinetic ballooning modes are stabilized by increased shaping due to an increased sensitivity to FLR effects, relative to the ion-temperature-gradient instability. Nevertheless, at high enough β, increased elongation degrades the local magnetic shear stabilization that enables access to the region of ballooning second-stability.