Optimizing stellarators for large flows

Plasma flow is damped in stellarators because they are not intrinsically ambipolar, unlike tokamaks, in which the flux-surface averaged radial electric current vanishes for any value of the radial electric field. Only quasisymmetric stellarators are intrinsically ambipolar, but exact quasisymmetry is impossible to achieve in non-axisymmetric toroidal configurations. By calculating the violation of intrinsic ambipolarity due to deviations from quasisymmetry, one can derive criteria to assess when a stellarator can be considered quasisymmetric in practice, i.e. when the flow damping is weak enough. Let us denote by $\alpha$ a small parameter that controls the size of a perturbation to an exactly quasisymmetric magnetic field. Recently, it has been shown that if the gradient of the perturbation is sufficiently small, the flux-surface averaged radial electric current scales as $\alpha^2$ for any value of the collisionality. It was also argued that when the gradient of the perturbation is large, the quadratic scaling is replaced by a more unfavorable one. In this paper, perturbations with large gradients are rigorously treated. In particular, it is proven that for low collisionality a perturbation with large gradient yields, at best, an $O(|\alpha|)$ deviation from quasisymmetry. Heuristic estimations in the literature incorrectly predicted an $O(|\alpha|^{3/2})$ deviation.Comment: 24 pages, 2 figures. To appear in Plasma Physics and Controlled Fusio

The effect of tangential drifts on neoclassical transport in stellarators close to omnigeneity

In general, the orbit-averaged radial magnetic drift of trapped particles in stellarators is non-zero due to the three-dimensional nature of the magnetic field. Stellarators in which the orbit-averaged radial magnetic drift vanishes are called omnigeneous, and they exhibit neoclassical transport levels comparable to those of axisymmetric tokamaks. However, the effect of deviations from omnigeneity cannot be neglected in practice. For sufficiently low collision frequencies (below the values that define the $1/\nu$ regime), the components of the drifts tangential to the flux surface become relevant. This article focuses on the study of such collisionality regimes in stellarators close to omnigeneity when the gradient of the non-omnigeneous perturbation is small. First, it is proven that closeness to omnigeneity is required to preserve radial locality in the drift-kinetic equation for collisionalities below the $1/\nu$ regime. Then, it is shown that neoclassical transport is determined by two layers in phase space. One of the layers corresponds to the $\sqrt{\nu}$ regime and the other to the superbanana-plateau regime. The importance of the superbanana-plateau layer for the calculation of the tangential electric field is emphasized, as well as the relevance of the latter for neoclassical transport in the collisionality regimes considered in this paper. In particular, the tangential electric field is essential for the emergence of a new subregime of superbanana-plateau transport when the radial electric field is small. A formula for the ion energy flux that includes the $\sqrt{\nu}$ regime and the superbanana-plateau regime is given. The energy flux scales with the square of the size of the deviation from omnigeneity. Finally, it is explained why below a certain collisionality value the formulation presented in this article ceases to be valid.Comment: 36 pages. Version to be published in Plasma Physics and Controlled Fusio

Flow damping in stellarators close to quasisymmetry

Quasisymmetric stellarators are a type of optimized stellarators for which flows are undamped to lowest order in an expansion in the normalized Larmor radius. However, perfect quasisymmetry is impossible. Since large flows may be desirable as a means to reduce turbulent transport, it is important to know when a stellarator can be considered to be sufficiently close to quasisymmetry. The answer to this question depends strongly on the size of the spatial gradients of the deviation from quasisymmetry and on the collisionality regime. Recently, formal criteria for closeness to quasisymmetry have been derived in a variety of situations. In particular, the case of deviations with large gradients was solved in the $1/\nu$ regime. Denoting by $\alpha$ a parameter that gives the size of the deviation from quasisymmetry, it was proven that particle fluxes do not scale with $\alpha^{3/2}$, as typically claimed, but with $\alpha$. It was also shown that ripple wells are not necessarily the main cause of transport. This paper reviews those works and presents a new result in another collisionality regime, in which particles trapped in ripple wells are collisional and the rest are collisionless.Comment: 14 pages, 2 figures. To appear in Plasma Physics and Controlled Fusio

Long-wavelength limit of gyrokinetics in a turbulent tokamak and its intrinsic ambipolarity

Recently, the electrostatic gyrokinetic Hamiltonian and change of coordinates have been computed to order $\epsilon^2$ in general magnetic geometry. Here $\epsilon$ is the gyrokinetic expansion parameter, the gyroradius over the macroscopic scale length. Starting from these results, the long-wavelength limit of the gyrokinetic Fokker-Planck and quasineutrality equations is taken for tokamak geometry. Employing the set of equations derived in the present article, it is possible to calculate the long-wavelength components of the distribution functions and of the poloidal electric field to order $\epsilon^2$. These higher-order pieces contain both neoclassical and turbulent contributions, and constitute one of the necessary ingredients (the other is given by the short-wavelength components up to second order) that will eventually enter a complete model for the radial transport of toroidal angular momentum in a tokamak in the low flow ordering. Finally, we provide an explicit and detailed proof that the system consisting of second-order gyrokinetic Fokker-Planck and quasineutrality equations leaves the long-wavelength radial electric field undetermined; that is, the turbulent tokamak is intrinsically ambipolar.Comment: 70 pages. Typos in equations (63), (90), (91), (92) and (129) correcte

Electrostatic potential variations on stellarator magnetic surfaces in low collisionality regimes

The component of the neoclassical electrostatic potential that is non-constant on the magnetic surface, that we denote by $\tilde\varphi$, can affect radial transport of highly charged impurities, and this has motivated its inclusion in some modern neoclassical codes. The number of neoclassical simulations in which $\tilde\varphi$ is calculated is still scarce, partly because they are usually demanding in terms of computational resources, especially at low collisionality. In this paper the size, the scaling with collisionality and with aspect ratio, and the structure of $\tilde\varphi$ on the magnetic surface are analytically derived in the $1/\nu$, $\sqrt{\nu}$ and superbanana-plateau regimes of stellarators close to omnigeneity; i. e. stellarators that have been optimized for neoclassical transport. It is found that the largest $\tilde\varphi$ that the neoclassical equations admit scales linearly with the inverse aspect ratio and with the size of the deviation from omnigeneity. Using a model for a perturbed omnigeneous configuration, the analytical results are verified and illustrated with calculations by the code KNOSOS. The techniques, results and numerical tools employed in this paper can be applied to neoclassical transport problems in tokamaks with broken axisymmetry.Comment: 30 pages, 12 figures, 1 table. Published versio

COLA. III. Radio Detection of Active Galactic Nucleus in Compact Moderate Luminosity Infrared Galaxies

We present results from 4.8 GHz Very Large Array (VLA) and global very long baseline interferometry (VLBI) observations of the northern half of the moderate FIR luminosity (median L_(IR) = 10^(11.01) L_☉) COLA sample of star-forming galaxies. VLBI sources are detected in a high fraction (20/90) of the galaxies observed. The radio luminosities of these cores (~10^(21) W Hz^(–1)) are too large to be explained by radio supernovae or supernova remnants and we argue that they are instead powered by active galactic nuclei (AGNs). These sub-parsec scale radio cores are preferentially detected toward galaxies whose VLA maps show bright 100-500 parsec scale nuclear radio components. Since these latter structures tightly follow the FIR to radio-continuum correlation for star formation, we conclude that the AGN-powered VLBI sources are associated with compact nuclear starburst environments. The implications for possible starburst-AGN connections are discussed. The detected VLBI sources have a relatively narrow range of radio luminosity consistent with models in which intense compact Eddington-limited starbursts regulate the gas supply onto a central supermassive black hole. The high incidence of AGN radio cores in compact starbursts suggests little or no delay between the starburst phase and the onset of AGN activity

Wave attenuation and dispersion due to floating ice covers

Experiments investigating the attenuation and dispersion of surface waves in a variety of ice covers are performed using a refrigerated wave flume. The ice conditions tested in the experiments cover naturally occurring combinations of continuous, fragmented, pancake and grease ice. Attenuation rates are shown to be a function of ice thickness, wave frequency, and the general rigidity of the ice cover. Dispersion changes were minor except for large wavelength increases when continuous covers were tested. Results are verified and compared with existing literature to show the extended range of investigation in terms of incident wave frequency and ice conditions