Ion temperature anisotropy in the tokamak scrape-off layer

Understanding tokamak exhaust-power heat loads on divertor plates depends critically on having a realistic model of the scrape-off layer (SOL) plasma. The Braginskii fluid model is often solved to understand the SOL plasma behavior. This model is based on the collisional limit for transport along the magnetic field B⃗. The ions and electron gyrofrequencies are assumed to be much larger than the Coulomb collision frequencies, which are nonetheless, sufficiently large to yield common parallel and perpendicular temperatures for each species, i.e., the temperatures are assumed to be isotropic. In certain circumstances such as encountered for the tokamak H-mode, the ion temperature can be quite anisotropic. In this work, the anisotropy effects are implemented in the 2D transport code UEDGE. Various geometries (1D slab, 2D slab and a toroidal tokamak geometry) are used to study the 2D structure of ion temperature anisotropy and its effects on plasma transport in detail. Results show that the effects of ion temperature anisotropy on the plasma parallel transport are substantial near the magnetic X-point, which leads to different steady state density profiles in the divertor regions. The extra mirror force introduced by ion temperature anisotropy can be one of the main forces contributing to the plasma flow in the SOL.

Development of advanced linearized gyrokinetic collision operators using a moment approach

The derivation and numerical implementation of a linearized version of the gyrokinetic (GK) Coulomb collision operator (Jorge et al., J. Plasma Phys., vol. 85, 2019, 905850604) and of the widely used linearized GK Sugama collision operator (Sugama et al., Phys. Plasmas, vol. 16, 2009, 112503) is reported. An approach based on a Hermite–Laguerre moment expansion of the perturbed gyrocentre distribution function is used, referred to as gyromoment expansion. This approach allows the considering of arbitrary perpendicular wavenumber and expressing the two linearized GK operators as a linear combination of gyromoments where the expansion coefficients are given by closed analytical expressions that depend on the perpendicular wavenumber and on the temperature and mass ratios of the colliding species. The drift-kinetic (DK) limits of the GK linearized Coulomb and Sugama operators are also obtained. Comparisons between the gyromoment approach and the DK Coulomb and GK Sugama operators in the continuum GK code GENE are reported, focusing on the ion-temperature-gradient instability and zonal flow damping, finding an excellent agreement. It is confirmed that stronger collisional damping of the zonal flow residual by the Sugama GK model compared with the GK linearized Coulomb (Pan et al., Phys. Plasmas, vol. 27, 2020, 042307) persists at higher collisionality. Finally, we show that the numerical efficiency of the gyromoment approach increases with collisionality, a desired property for boundary plasma applications.

Magnetic field generation from composition gradients in inertial confinement fusion fuel

Experimental asymmetries in fusion implosions can lead to magnetic field generation in the hot plasma core. For typical parameters, previous studies found that the magnetization Hall parameter, given by the product of the electron gyro-frequency and Coulomb collision time, can exceed one. This will affect the hydrodynamics through inhibition and deflection of the electron heat flux. The magnetic field source is the collisionless Biermann term, which arises from the Debye shielding potential in electron pressure gradients. We show that there is an additional source term due to the Z dependence of the Coulomb collision operator. If there are ion composition gradients, such as jets of carbon ablator mix entering the hot-spot, this source term can rapidly exceed the Biermann fields. In addition, the Biermann fields are enhanced due to the increased temperature gradients from carbon radiative cooling. With even stronger self-generated fields, heat loss to the carbon regions will be reduced, potentially reducing the negative effect of carbon mix. This article is part of a discussion meeting issue ‘Prospects for high gain inertial fusion energy (part 1)’.

Femtosecond pumping of nuclear isomeric states by the Coulomb collision of ions with quivering electrons

Efficient production of metastable quantum states of nuclei (isomers) is critical for exotic applica- tions, like nuclear clocks, nuclear batteries, clean nuclear energy, and nuclear gamma-ray lasers[1–6]. However, due to low reaction cross sections and quick decay, it is extremely difficult to acquire sig- nificant amount of isomers with short lifetimes via traditional accelerators or reactors. Here, we present femtosecond pumping of nuclear isomeric states by the Coulomb excitation of ions with the quivering electrons induced by laser fields for the first time. Nuclear isomers populated on the second excited state of 83Kr, are generated with a rate of 3.84 × 10^17 per second from a table-top hundreds-TW laser system. This high efficiency of isomer production can be explained by Coulomb collision[7] of ions with the quivering electrons during the laser-cluster interactions at nearly solid densities.

Nonlinear gyrokinetic Coulomb collision operator

A gyrokinetic Coulomb collision operator is derived, which is particularly useful to describe the plasma dynamics at the periphery region of magnetic confinement fusion devices. The derived operator is able to describe collisions occurring in distribution functions arbitrarily far from equilibrium with variations on spatial scales at and below the particle Larmor radius. A multipole expansion of the Rosenbluth potentials is used in order to derive the dependence of the full Coulomb collision operator on the particle gyroangle. The full Coulomb collision operator is then expressed in gyrocentre phase-space coordinates, and a closed formula for its gyroaverage in terms of the moments of the gyrocentre distribution function in a form ready to be numerically implemented is provided. Furthermore, the collision operator is projected onto a Hermite–Laguerre velocity space polynomial basis and expansions in the small electron-to-ion mass ratio are provided.