Search

The first stellarator design was a simple tube of plasma twisted and closed on itself in the form of a figure-8. The line of such devices, however, was quickly ended over concerns related to plasma stability. We revisit the figure-8 concept, re-imagined as a modern optimized stellarator, and find the potential for a high degree of stability, as well as exceptionally simple construction. In particular, the design that we find admits planar coils, and is the first quasi-isodynamic stellarator design to have this property. Our work is made possible by recent theoretical progress in the near-axis theory of quasi-isodynamic stellarators, combined with fundamental progress in the numerical solution of three-dimensional magnetohydrodynamic equilibria that cannot be well represented using traditional cylindrical coordinates.

Tokamaks and stellarators are the leading two magnetic confinement devices for producing fusion energy, begging the question of whether the strengths of the two could be merged into a single concept. To meet this challenge, we propose a first-of-its kind optimized stellarator-tokamak hybrid. Compared to a typical tokamak coil set, only a single simple type of stellarator coil has to be added which leads to a compact, volume- and transport-preserving magnetic field, with an added rotational transform that reaches levels thought to enhance stability.

Upper bounds on the growth of instabilities in gyrokinetic systems have recently been derived by considering the optimal perturbations that maximise the growth of a chosen energy norm. This technique has previously been applied to two-species gyrokinetic systems with fully kinetic ions and electrons. However, in tokamaks and stellarators, the expectation from linear instability analyses is that the most important kinetic-electron contribution to ion-scale modes comes from the trapped electrons, which bounce faster than the timescale upon which instabilities evolve. As a result, a fully-kinetic electron response is not required to describe unstable modes in most cases. Here, we apply the optimal mode analysis to a reduced two-species system that consists of fully gyrokinetic ions and bounce-averaged electrons with the aim of finding a tighter bound on ion-scale instabilities in toroidal geometry. This analysis yields bounds that are greatly reduced in comparison to the earlier two-species result. Moreover, if the energy norm is properly chosen, wave-particle resonance effects can be captured, reproducing the stabilisation of density-gradient-driven instabilities in maximum-$J$ devices. The optimal mode analysis also reveals that the maximum-$J$ property has an additional stabilising effect on ion-temperature-gradient-driven instabilities, even in the absence of an electron-free energy source. This effect is explained in terms of the concept of mode inertia, making it distinct from other mechanisms.

We investigate the confinement properties of two recently devised quasi-helically symmetric stellarator configurations, HSK and QSTK. Both have been optimized for large critical gradients of the ion temperature gradient mode, which is an important driver of turbulent transport in magnetic confinement fusion devices. To predict the resulting core plasma profiles, we utilize an advanced theoretical framework based on the gyrokinetic codes GENE and GENE-3D, coupled to the transport code TANGO. Compared to the HSX stellarator, both HSK and QSTK achieve significantly higher core-to-edge temperature ratios, partly thanks to their smaller aspect ratios, with the other part due to more detailed shaping of the magnetic geometry achieved during optimization. The computed confinement time, however, is less sensitive to core temperature than edge temperature, simply due to the disproportionate influence the edge has on stored plasma energy. We therefore emphasize the possible benefits of further optimizing turbulence in the outer core region, and the need to include accurate modelling of confinement in the edge region in order to assess overall plasma performance of turbulence optimized stellarators.

The gyrokinetic theory of the residual flow, in the electrostatic limit, is revisited, with optimized stellarators in mind. We consider general initial conditions for the problem, and identify cases that lead to a non-zonal residual electrostatic potential, i.e. one having a significant component that varies within a flux surface. We investigate the behavior of the ``intermediate residual'' in stellarators, a measure of the flow that remains after geodesic acoustic modes have damped away, but before the action of the slower damping that is caused by unconfined particle orbits. The case of a quasi-isodynamic stellarator is identified as having a particularly large such residual, owing to the small orbit width achieved by optimization.

Integrating turbulence into stellarator optimization is shown by targeting the onset for the ion-temperature-gradient mode, highlighting effects of parallel connection length, local magnetic shear, and flux surface expansion. The result is a compact quasihelically symmetric stellarator configuration, admitting a set of uncomplicated coils, with significantly reduced turbulent heat fluxes compared to a known stellarator. The new configuration combines low values of neoclassical transport, good alpha particle confinement, and Mercier stability at a plasma beta of almost 2$\%$., Comment: 5 pages, 5 figures. Phys. Rev. Research 5, L032030 (2023)

Free energy, widely used as a measure of turbulence intensity in weakly collisional plasmas, has been recently found to be a suitable basis to describe both linear and nonlinear growth in a wide class gyrokinetic systems. The simplicity afforded by this approach is accompanied by some drawbacks, notably the lack of any explicit treatment of wave-particle effects, which makes the theory unable to describe things like stability thresholds or dependence on the geometry of the background magnetic field. As a step toward overcoming these limitations, we propose an extension of the theory based on a generalization of free energy. With this it is demonstrated that resonance effects are recovered, and the bounds on growth are significantly reduced. The simplicity and efficient computation of the associated "optimal" growth rates makes the theory potentially applicable to stellarator optimization.

We present a stellarator configuration optimized for a large threshold (``critical gradient'') for the onset of the ion temperature gradient (ITG) driven mode, which achieves the largest critical gradient we have seen in any stellarator. Above this threshold, gyrokinetic simulations show that the configuration has low turbulence levels over an experimentally relevant range of the drive strength. The applied optimization seeks to maximize the drift curvature, leading to enhanced local-shear stabilization of toroidal ITG modes, and the associated turbulence. These benefits are combined with excellent quasisymmetry, yielding low neoclassical transport and vanishingly small alpha particle losses. Analysis of the resulting configuration suggests a trade-off between magnetohydrodynamic (MHD) and ITG stability., Comment: 11 pages, 7 figures. Submitted to the Journal of Plasma Physics

We present a modified gyrokinetic theory to predict the critical gradient that determines the linear onset of the ion temperature gradient (ITG) mode in stellarator plasmas. A coarse-graining technique is applied to the drift curvature, entering the standard gyrokinetic equations, around local minima. Thanks to its simplicity, this novel formalism yields an estimate for the critical gradient with a computational cost low enough for application to stellarator optimization. On comparing against a gyrokinetic solver, our results show good agreement for an assortment of stellarator designs. Insight gained here into the physics of the onset of the ITG driven instability enables us to devise a compact configuration, similar to the Wendelstein 7-X device, but with almost twice the ITG linear critical gradient, an improved nonlinear critical gradient, and reduced ITG mode transport above the nonlinear critical gradient., Comment: 5 pages, 4 figures. Accepted for publication in Physical Review Research

A single-field-period quasi-isodynamic stellarator configuration is presented. This configuration, which resembles a twisted strip, is obtained by the method of direct construction, that is, it is found via an expansion in the distance from the magnetic axis. Its discovery, however, relied on an additional step involving numerical optimization, performed within the space of near-axis configurations defined by a set of adjustable magnetic-field parameters. This optimization, completed in 30 seconds on a single cpu core using the SIMSOPT code, yields a solution with excellent confinement, as measured by the conventional figure of merit for neoclassical transport, effective ripple, at a modest aspect ratio of eight. The optimization parameters that led to this configuration are described, its confinement properties are assessed, and a set of magnetic-field coils is found. The resulting transport at low collisionality is much smaller than that of W7-X, and the device needs significantly fewer coils thanks to the reduced number of field periods., Comment: 13 pages, 9 figures

We introduce modes of instantaneous optimal growth of free energy for the fully electromagnetic gyrokinetic equations. We demonstrate how these "optimal modes" arise naturally from the free energy balance equation, allowing its convenient decomposition, and yielding a simple picture of energy flows. Optimal modes have a number of other favorable features, such as their low-dimensionality, efficiency of computation, and the fact that their growth rates provide a rigorous and "tight" upper bound on both the nonlinear growth rate of energy, and the linear growth rate of traditional gyrokinetic (normal mode) instabilities. We provide simple closed form solutions for the optimal growth rates in a number of asymptotic limits, and compare these with our previous bounds.

Upper bounds on the growth of free energy in gyrokinetics are derived. These bounds apply to all local gyrokinetic instabilities in the geometry of a flux tube, i.e. a slender volume of plasma aligned with the magnetic field, regardless of the geometry of field, the number of particle species, or collisions. The results apply both to linear instabilities and to the nonlinear growth of finite-amplitude fluctuations.

The Dimits shift, an upshift in the onset of turbulence from the linear instability threshold, caused by self-generated zonal flows, can greatly enhance the performance of magnetic confinement plasma devices. Except in simple cases, using fluid approximations and model magnetic geometries, this phenomenon has proved difficult to understand and quantitatively predict. To bridge the large gap in complexity between simple models and realistic treatment in toroidal magnetic geometries (e.g. tokamaks or stellarators), the present work uses fully gyrokinetic simulations in Z-pinch geometry to investigate the Dimits shift through the lens of tertiary instability analysis, which describes the emergence of drift waves from a zonally dominated state. Several features of the tertiary instability, previously observed in fluid models, are confirmed to remain. Most significantly, an efficient reduced-mode tertiary model, which previously proved successful in predicting the Dimits shift in a gyrofluid limit (Hallenbert & Plunk, J. Plasma Phys., vol.87, issue 05, 2021, 905870508), is found to be accurate here, with only slight modifications to account for kinetic effects., Comment: 36 pages, 15 figures

A family of rigorous upper bounds on the growth rate of local gyrokinetic instabilities in magnetized plasmas is derived from the evolution equation for the Helmholtz free energy. These bounds hold for both electrostatic and electromagnetic instabilities, regardless of the number of particle species, their collision frequency, and the geometry of the magnetic field. A large number of results that have earlier been derived in special cases and observed in numerical simulations are thus brought into a unifying framework. These bounds apply not only to linear instabilities but also imply an upper limit to the nonlinear growth of the free energy.

The tertiary instability is believed to be important for governing magnetised plasma turbulence under conditions of strong zonal flow generation, near marginal stability. In this work, we investigate its role for a collisionless strongly driven fluid model, self-consistently derived as a limit of gyrokinetics. It is found that a region of absolute stability above the linear threshold exists, beyond which significant nonlinear transport rapidly develops. While within this range a complex pattern of transient zonal evolution is observed before a stable profile is found, the Dimits transition itself is found to coincide with a tertiary instability threshold so long as linear effects are included. Through a simple and readily extendable procedure tracing its origin to St-Onge 2017 (arXiv:1704.05406) the stabilising effect of the typical zonal profile can be approximated and the accompanying reduced mode estimate is found to be in good agreement with nonlinear simulations., Comment: 30 pages, 9 figures; improved the derivation of Appendix A, slightly modified several statements for accuracy, updated figures and notation

A first-principles method to calculate the critical temperature gradient for the onset of the ion-temperature-gradient mode (ITG) in linear gyrokinetics is presented. We find that conventional notions of the connection length previously invoked in tokamak research should be revised and replaced by a generalized correlation length to explain this onset in stellarators. Simple numerical experiments and gyrokinetic theory show that localized "spikes" in shear, a hallmark of stellarator geometry, are generally insufficient to constrain the parallel correlation length of the mode. ITG modes that localize within bad drift curvature wells that have a critical gradient set by peak drift curvature are also observed. A case study of nearly helical stellarators of increasing field period demonstrates that the critical gradient can indeed be controlled by manipulating magnetic geometry, but underscores the need for a general framework to evaluate the critical gradient. We conclude that average curvature and global shear set the correlation length of resonant ITG modes near the absolute critical gradient, the physics of which is included through direct solution of the gyrokinetic equation. Our method, which handles general geometry and is more efficient than conventional gyrokinetic solvers, could be applied to future studies of stellarator ITG turbulence optimization., Comment: 22 pages, 6 figures. Accepted for publication in The Journal of Plasma Physics (JPP)

In stellarators, zonal flow activity depends sensitively on geometry of the three dimensional magnetic field, via an interplay of mechanisms that is not fully understood. In this work, we investigate this by studying three magnetic configurations of the Wendelstein 7-X stellarator. We find that variation in linear zonal flow damping is accompanied by variation in nonlinear drive, and identify key geometric features that control these effects. Understanding the resulting balance is important for the development of reduced models of turbulent transport., Comment: 11 pages, 9 figures

A fluid system is derived to describe electrostatic magnetized plasma turbulence at scales somewhat larger than the Larmor radius of a given species. It is related to the Hasegawa- Mima equation, but does not conserve enstrophy, and, as a result, exhibits a forward cascade of energy, to small scales. The inertial-range energy spectrum is argued to be shallower than a -11/3 power law, as compared to the -5 law of the Hasegawa-Mima enstrophy cascade. This property, confirmed here by direct numerical simulations of the fluid system, may help explain the fluctuation spectrum observed in gyrokinetic simulations of streamer-dominated electron-temperature-gradient driven turbulence [Plunk et al., 2019], and also possibly some cases of ion-temperature-gradient driven turbulence where zonal flows are suppressed [Plunk et al., 2017].

It is demonstrated that finite-pressure, approximately quasi-axisymmetric stellarator equilibria can be directly constructed (without numerical optimization) via perturbations of given axisymmetric equilibria. The size of such perturbations is measured in two ways, via the fractional external rotation and, alternatively, via the relative magnetic field strength, i.e. the average size of the perturbed magnetic field, divided by the unperturbed field strength. It is found that significant fractional external rotational transform can be generated by quasi-axisymmetric perturbations, with a similar value of the relative field strength, despite the fact that the former scales more weakly with the perturbation size. High mode number perturbations are identified as a candidate for generating such transform with local current distributions. Implications for the development of a general non-perturbative solver for optimal stellarator equilibria is discussed.

We study the effect of turbulent transport in different magnetic configurations of the Weldenstein 7-X stellarator. In particular, we performed direct numerical simulations with the global gyrokinetic code GENE-3D, modeling the behavior of Ion Temperature Gradient turbulence in the Standard, High-Mirror, and Low-Mirror configurations of W7-X. We found that the Low-Mirror configuration produces more transport than both the High-Mirror and the Standard configurations. By comparison with radially local simulations, we have demonstrated the importance of performing global nonlinear simulations to predict the turbulent fluxes quantitatively.

In fusion devices, the geometry of the confining magnetic field has a significant impact on the instabilities that drive turbulent heat loss. This is especially true of stellarators, where the "trapped electron mode" (TEM) is stabilised if specific optimisation criteria are satisfied, as in the Wendelstein 7-X experiment (W7-X). Here we find, by numerical simulation, that W7-X indeed has low TEM-driven transport, and also benefits from stabilisation of the ion-temperature-gradient mode, giving theoretical support for the existence of enhanced confinement regimes at finite density gradients.

The condition of omnigenity is investigated, and applied to the near-axis expansion of Garren and Boozer (1991a). Due in part to the particular analyticity requirements of the near-axis expansion, we find that, excluding quasi-symmetric solutions, only one type of omnigenity, namely quasi-isodynamicity, can be satisfied at first order in the distance from the magnetic axis. Our construction provides a parameterization of the space of such solutions, and the cylindrical reformulation and numerical method of Landreman and Sengupta (2018); Landreman et al. (2019), enables their efficient numerical construction.

An asymptotic expansion is performed to obtain quasi-axisymmetric magnetic configurations that are weakly non-axisymmetric. A large space of solutions is identified, which satisfy the condition of quasi-axisymmetry on a single magnetic flux surface, while (non-axisymmetric) globally quasi-axisymmetric solutions are shown to not exist, agreeing with the conclusions of previous theoretical work. The solutions found are shown to be geometrically constrained at low aspect ratio or high toroidal period number. Solutions satisfying the more general condition of omnigeneity (generalized quasi-axisymmetry) are also shown to exist, and it is found that quasi-axisymmetric deformations can be superposed with an omnigenous solution, while preserving the property of omnigeneity, effectively extending the space of "good" configurations. A numerical solution of the first order quasi-axisymmetry problem is demonstrated and compared with solutions found with a widely used MHD equilibrium solver, independently verifying that quasi-axisymmetry is satisfied at the appropriate order. It is thereby demonstrated that approximately quasi-axisymmetric solutions can be directly constructed, i.e. without using numerical search algorithms.

Optimised stellarators and other magnetic-confinement devices having the property that the average magnetic curvature is favourable for all particle orbits are called maximum-$J$ devices, and have recently been shown to be immune to trapped-particle instabilities driven by the density gradient. Gyrokinetic simulations reveal, however, that another instability can arise, which is also associated with particle trapping but causes less transport than typical trapped-electron modes. The nature of this instability is clarified here. It is shown to be similar to the "ubiquitous mode" in tokamaks, and is driven by ion free energy but requires trapped electrons to exist., Comment: Accepted to Journal of Plasma Physics

We demonstrate that the scaling properties of slab ion and electron temperature gradient driven turbulence may be derived by dimensional analysis of a drift kinetic system with one kinetic species. These properties have previously been observed in gyrokinetic simulations of turbulence in magnetic fusion devices., Comment: To be published, Phys. Plasmas (2017)

In the complex 3D magnetic fields of stellarators, ion-temperature-gradient turbulence is shown to have two distinct saturation regimes, as revealed by petascale numerical simulations, and explained by a simple turbulence theory. The first regime is marked by strong zonal flows, and matches previous observations in tokamaks. The newly observed second regime, in contrast, exhibits small- scale quasi-two-dimensional turbulence, negligible zonal flows, and, surprisingly, a weaker heat flux scaling. Our findings suggest that key details of the magnetic geometry control turbulence in stellarators., Comment: Erratum added to end

We present a theory of the nonlinear growth of zonal flows in magnetized plasma turbulence, by the mechanism of secondary instability. The theory is derived for general magnetic geometry, and is thus applicable to both tokamaks and stellarators. The predicted growth rate is shown to compare favorably with nonlinear gyrokinetic simulations, with the error scaling as expected with the small parameter of the theory., Comment: New J. Phys. 2017

It is shown that in perfectly quasi-isodynamic stellarators, trapped particles with a bounce frequency much higher than the frequency of the instability are stabilizing in the electrostatic and collisionless limit. The collisionless trapped-particle instability is therefore stable as well as the ordinary electron-density-gradient-driven trapped-electron mode. This result follows from the energy balance of electrostatic instabilities and is thus independent of all other details of the magnetic geometry.

The "universal" instability has recently been revived by Landreman, Antonsen and Dorland [1], who showed that it indeed exists in plasma geometries with straight (but sheared) magnetic field lines. Here it is demonstrated analytically that this instability can be present in more general sheared and toroidal geometries. In a torus, the universal instability is shown to be closely related to the trapped-electron mode, although the trapped-electron drive is usually dominant. However, this drive can be weakened or eliminated, as in the case in stellarators with the maximum-$J$ property, leaving the parallel Landau resonance to drive a residual mode, which is identified as the universal instability.

We study a quasi-two-dimensional electrostatic drift kinetic system as a model for near-marginal ion temperature gradient (ITG) driven turbulence. A proof is given of the nonlinear stability of this system under conditions of linear stability. This proof is achieved using a transformation that diagonalizes the linear dynamics and also commutes with nonlinear $E\times B$ advection. For the case when linear instability is present, a corollary is found that forbids nonlinear energy transfer between appropriately defined sets of stable and unstable modes. It is speculated that this may explain the preservation of linear eigenmodes in nonlinear gyrokinetic simulations. Based on this property, a dimensionally reduced ($\infty\times\infty \rightarrow 1$) system is derived that may be useful for understanding dynamics around the critical gradient of Dimits.

A recent paper [Phys. Plasmas 20, 032304 (2013)] considered the forced linear Vlasov equation as a model for the quasi-steady state of a single stable plasma wavenumber interacting with a bath of turbulent fluctuations. This approach gives some insight into possible energy flows without solving for nonlinear dynamics. The central result of the present work is that the forced linear Vlasov equation exhibits asymptotically zero (irreversible) dissipation to all orders under a detuning of the forcing frequency and the characteristic frequency associated with particle streaming. We first prove this by direct calculation, tracking energy flow in terms of certain exact conservation laws of the linear (collisionless) Vlasov equation. Then we analyze the steady-state solutions in detail using a weakly collisional Hermite-moment formulation, and compare with numerical solution. This leads to a detailed description of the Hermite energy spectrum, and a proof of no dissipation at all orders, complementing the collisionless Vlasov result., Comment: Small changes for clarity

We propose a quantitative model of ion temperature gradient driven turbulence in toroidal magnetized plasmas. In this model, the turbulence is regulated by zonal flows, i.e. mode saturation occurs by a zonal-flow-mediated energy cascade ("shearing"), and zonal flow amplitude is controlled by nonlinear decay. Our model is tested in detail against numerical simulations to confirm that both its assumptions and predictions are satisfied. Key results include (1) a sensitivity of the nonlinear zonal flow response to the energy content of the linear instability, (2) a persistence of zonal-flow-regulated saturation at high temperature gradients, (3) a physical explanation of the nonlinear saturation process in terms of secondary and tertiary instabilities, and (4) dependence of heat flux in terms of dimensionless parameters., Comment: Final journal version. Some clarifications and a new Fig. 4

We investigate the linear theory of the ion-temperature-gradient (ITG) mode, with the goal of developing a general understanding that may be applied to stellarators. We highlight the Wendelstein 7X (W7-X) device. Simple fluid and kinetic models that follow closely from existing literature are reviewed and two new first-principle models are presented and compared with results from direct numerical simulation. One model investigates the effect of regions of strong localized shear, which are generic to stellarator equilibria. These "shear spikes" are found to have a potentially significant stabilizing affect on the mode; however, the effect is strongest at short wavelengths perpendicular to the magnetic field, and it is found to be significant only for the fastest growing modes in W7-X. A second model investigates the long-wavelength limit for the case of negligible global magnetic shear. The analytic calculation reveals that the effect of the curvature drive enters at second order in the drift frequency, confirming conventional wisdom that the ITG mode is slab-like at long wavelengths. Using flux tube simulations of a zero-shear W7-X configuration, we observe a close relationship to an axisymmetric configuration at a similar parameter point. It is concluded that scale lengths of the equilibrium gradients constitute a good parameter space to characterize the ITG mode. Thus, to optimize the magnetic geometry for ITG mode stability, it may be fruitful to focus on local parameters, such as the magnitude of bad curvature, connection length, and local shear at locations of bad curvature (where the ITG mode amplitude peaks)., Comment: v3: final version for journal

This is the first of two papers about collisionless, electrostatic micro-instabilities in stellarators, with an emphasis on trapped-particle modes. It is found that, in so-called maximum-$J$ configurations, trapped-particle instabilities are absent in large regions of parameter space. Quasi-isodynamic stellarators have this property (approximately), and the theory predicts that trapped electrons are stabilizing to all eigenmodes with frequencies below the electron bounce frequency. The physical reason is that the bounce-averaged curvature is favorable for all orbits, and that trapped electrons precess in the direction opposite to that in which drift waves propagate, thus precluding wave-particle resonance. These considerations only depend on the electrostatic energy balance, and are independent of all geometric properties of the magnetic field other than the maximum-$J$ condition. However, if the aspect ratio is large and the instability phase velocity differs greatly from the electron and ion thermal speeds, it is possible to derive a variational form for the frequency showing that stability prevails in a yet larger part of parameter space than what follows from the energy argument. Collisionless trapped-electron modes should therefore be more stable in quasi-isodynamic stellarators than in tokamaks., Comment: 9 pages, 1 figure

This paper presents a complete theoretical framework for plasma turbulence and transport in tokamak plasmas. The fundamental scale separations present in plasma turbulence are codified as an asymptotic expansion in the ratio of the gyroradius to the equilibrium scale length. Proceeding order-by-order in this expansion, a framework for plasma turbulence is developed. It comprises an instantaneous equilibrium, the fluctuations driven by gradients in the equilibrium quantities, and the transport-timescale evolution of mean profiles of these quantities driven by the fluctuations. The equilibrium distribution functions are local Maxwellians with each flux surface rotating toroidally as a rigid body. The magnetic equillibrium is obtained from the Grad-Shafranov equation for a rotating plasma and the slow (resistive) evolution of the magnetic field is given by an evolution equation for the safety factor q. Large-scale deviations of the distribution function from a Maxwellian are given by neoclassical theory. The fluctuations are determined by the high-flow gyrokinetic equation, from which we derive the governing principle for gyrokinetic turbulence in tokamaks: the conservation and local cascade of free energy. Transport equations for the evolution of the mean density, temperature and flow velocity profiles are derived. These transport equations show how the neoclassical corrections and the fluctuations act back upon the mean profiles through fluxes and heating. The energy and entropy conservation laws for the mean profiles are derived. Total energy is conserved and there is no net turbulent heating. Entropy is produced by the action of fluxes flattening gradients, Ohmic heating, and the equilibration of mean temperatures. Finally, this framework is condensed, in the low-Mach-number limit, to a concise set of equations suitable for numerical implementation., Comment: 113 pages, 3 figures

In magnetized plasmas, a turbulent cascade occurs in phase space at scales smaller than the thermal Larmor radius ("sub-Larmor scales") [Phys. Rev. Lett. 103, 015003 (2009)]. When the turbulence is restricted to two spatial dimensions perpendicular to the background magnetic field, two independent cascades may take place simultaneously because of the presence of two collisionless invariants. In the present work, freely decaying turbulence of two-dimensional electrostatic gyrokinetics is investigated by means of phenomenological theory and direct numerical simulations. A dual cascade (forward and inverse cascades) is observed in velocity space as well as in position space, which we diagnose by means of nonlinear transfer functions for the collisionless invariants. We find that the turbulence tends to a time-asymptotic state, dominated by a single scale that grows in time. A theory of this asymptotic state is derived in the form of decay laws. Each case that we study falls into one of three regimes (weakly collisional, marginal, and strongly collisional), determined by a dimensionless number D*, a quantity analogous to the Reynolds number. The marginal state is marked by a critical number D* = D0 that is preserved in time. Turbulence initialized above this value become increasingly inertial in time, evolving toward larger and larger D*; turbulence initialized below D0 become more and more collisional, decaying to progressively smaller D*., Comment: 12 pages, 12 figures; replaced to match published version

To address the problem of Landau damping in kinetic turbulence, the forcing of the linearized Vlasov equation by a stationary random source is considered. It is found that the time-asymptotic density response is dominated by resonant particle interactions that are synchronized with the source. The energy consumption of this response is calculated, implying an effective damping rate, which is the main result of this paper. Evaluating several cases, it is found that the effective damping rate can differ from the Landau damping rate in magnitude and also, remarkably, in sign. A limit is demonstrated in which the density and current become phase-locked, which causes the effective damping to be negligible; this potentially resolves an energy paradox that arises in the application of critical balance to a kinetic turbulence cascade., Comment: Introduction significantly expanded to help contextualize results. Calculations unchanged

In gyrokinetic theory there are two quadratic measures of fluctuation energy, left invariant under nonlinear interactions, that constrain the turbulence. The recent work of Plunk and Tatsuno [Phys. Rev. Lett. 106, 165003 (2011)] reported on the novel consequences that this constraint has on the direction and locality of spectral energy transfer. This paper builds on that work. We provide detailed analysis in support of the results of Plunk and Tatsuno but also significantly broaden the scope and use additional methods to address the problem of energy transfer. The perspective taken here is that the fluctuation energies are not merely formal invariants of an idealized model (two-dimensional gyrokinetics) but are general measures of gyrokinetic turbulence, i.e. quantities that can be used to predict the behavior of the turbulence. Though many open questions remain, this paper collects evidence in favor of this perspective by demonstrating in several contexts that constrained spectral energy transfer governs the dynamics., Comment: Final version as published. Some cosmetic changes and update of references

We respond to the recent comment [arXiv:1105.1593] on our Letter [G. G. Plunk and T. Tatsuno, Phys. Rev. Lett. 106, 165003 (2011)]. The comment claims that our argument for spectral transfer direction is incomplete. The comment gives an incomplete account of our argument. We explain how the argument of our Letter, in full, is self-contained consistent and well-supported by numerical evidence. Other technical statements made in the comment are explained to be incorrect. We direct interested readers to our forthcoming article, which should clarify any remaining technical ambiguities and also illuminate deeper physical issues.

The question of how nonlinear interactions redistribute the energy of fluctuations across available degrees of freedom is of fundamental importance in the study of turbulence and transport in magnetized weakly collisional plasmas, ranging from space settings to fusion devices. In this letter, we present a theory for the dual cascade found in such plasmas, which predicts a range of new behavior that distinguishes this cascade from that of neutral fluid turbulence. These phenomena are explained in terms of the constrained nature of spectral transfer in nonlinear gyrokinetics. Accompanying this theory are the first observations of these phenomena, obtained via direct numerical simulations using the gyrokinetic code {\tt AstroGK}. The basic mechanisms that are found provide a framework for understanding the turbulent energy transfer that couples scales both locally and non-locally., Comment: Reduced size for journal format

Two-dimensional electrostatic turbulence in magnetized weakly-collisional plasmas exhibits a cascade of entropy in phase space [Phys. Rev. Lett. 103, 015003 (2009)]. At scales smaller than the gyroradius, this cascade is characterized by the dimensionless ratio D of the collision time to the eddy turnover time measured at the scale of the thermal Larmor radius. When D >> 1, a broad spectrum of fluctuations at sub-Larmor scales is found in both position and velocity space. The distribution function develops structure as a function of v_{perp}, the velocity coordinate perpendicular to the local magnetic field. The cascade shows a local-scale nonlinear interaction in both position and velocity spaces, and Kolmogorov's scaling theory can be extended into phase space., Comment: 8 pages, 10 figures, Conference paper presented at 2009 Asia-Pacific Plasma Theory Conference. Ver.2 includes corrected typos & updated references

Two-dimensional gyrokinetics is a simple paradigm for the study of kinetic magnetised plasma turbulence. In this paper, we present a comprehensive theoretical framework for this turbulence. We study both the inverse and direct cascades (the `dual cascade'), driven by a homogeneous and isotropic random forcing. The key characteristic length of gyrokinetics, the Larmor radius, divides scales into two physically distinct ranges. For scales larger than the Larmor radius, we derive the familiar Charney--Hasegawa--Mima (CHM) equation from the gyrokinetic system, and explain its relationship to gyrokinetics. At scales smaller than the Larmor radius, a dual cascade occurs in phase space (two dimensions in position space plus one dimension in velocity space) via a nonlinear phase-mixing process. We show that at these sub-Larmor scales, the turbulence is self-similar and exhibits power law spectra in position and velocity space. We propose a Hankel-transform formalism to characterise velocity-space spectra. We derive the exact relations for third-order structure functions, analogous to Kolmogorov's four-fifths and Yaglom's four-thirds laws and valid at both long and short wavelengths. We show how the general gyrokinetic invariants are related to the particular invariants that control the dual cascade in the long- and short-wavelength limits. We describe the full range of cascades from the fluid to the fully kinetic range., Comment: Being considered for publication in the Journal of Fluid Mechanics. Version 3 contains more cross referencing including forward referencing, to make the paper easier to read. Final pre-production version

Electrostatic turbulence in weakly collisional, magnetized plasma can be interpreted as a cascade of entropy in phase space, which is proposed as a universal mechanism for dissipation of energy in magnetized plasma turbulence. When the nonlinear decorrelation time at the scale of the thermal Larmor radius is shorter than the collision time, a broad spectrum of fluctuations at sub-Larmor scales is numerically found in velocity and position space, with theoretically predicted scalings. The results are important because they identify what is probably a universal Kolmogorov-like regime for kinetic turbulence; and because any physical process that produces fluctuations of the gyrophase-independent part of the distribution function may, via the entropy cascade, result in turbulent heating at a rate that increases with the fluctuation amplitude, but is independent of the collision frequency., Comment: Revtex, 4 pages, 3 figures; replaced to match published version

This paper describes a conceptual framework for understanding kinetic plasma turbulence as a generalized form of energy cascade in phase space. It is emphasized that conversion of turbulent energy into thermodynamic heat is only achievable in the presence of some (however small) degree of collisionality. The smallness of the collision rate is compensated by the emergence of small-scale structure in the velocity space. For gyrokinetic turbulence, a nonlinear perpendicular phase mixing mechanism is identified and described as a turbulent cascade of entropy fluctuations simultaneously occurring at spatial scales smaller than the ion gyroscale and in velocity space. Scaling relations for the resulting fluctuation spectra are derived. An estimate for the collisional cutoff is provided. The importance of adequately modeling and resolving collisions in gyrokinetic simulations is biefly discussed, as well as the relevance of these results to understanding the dissipation-range turbulence in the solar wind and the electrostatic microturbulence in fusion plasmas., Comment: iop revtex style, 14 pages, 1 figure; submitted to PPCF; invited talk for EPS Conference on Plasma Physics, Crete, June 2008; Replaced to match published version