Showing 1–50 of 51 results for author: Plunk, G G

A compact stellarator-tokamak hybrid

Authors: S. A. Henneberg, G. G. Plunk

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2024; originally announced June 2024.

Journal ref: Physical Review Research (Vol. 6, No. 2), (2024)

Quasi-isodynamic stellarators with low turbulence as fusion reactor candidates

Authors: Alan G. Goodman, Pavlos Xanthopoulos, Gabriel G. Plunk, Håkan Smith, Carolin Nührenberg, Craig D. Beidler, Sophia A. Henneberg, Gareth Roberg-Clark, Michael Drevlak, Per Helander

Abstract: The stellarator is a type of fusion energy device that - if properly designed - could provide clean, safe, and abundant energy to the grid. To generate this energy, a stellarator must keep a hot mixture of charged particles (known as a plasma) sufficiently confined by using a fully shaped magnetic field. If this is achieved, the heat from fusion reactions within the plasma can be harvested as ener… ▽ More The stellarator is a type of fusion energy device that - if properly designed - could provide clean, safe, and abundant energy to the grid. To generate this energy, a stellarator must keep a hot mixture of charged particles (known as a plasma) sufficiently confined by using a fully shaped magnetic field. If this is achieved, the heat from fusion reactions within the plasma can be harvested as energy. We present a novel method for designing reactor-relevant stellarator magnetic fields, which combine several key physical properties. These include plasma stability, excellent confinement of the fast moving particles generated by fusion reactions, and reduction of the turbulence that is known to limit the performance of the most advanced stellarator experiment in the world, Wendelstein 7-X. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Comments: 11 pages, 15 figures

Energetic bounds on gyrokinetic instabilities. Part 4. Bounce-averaged electrons

Authors: P. J. Costello, G. G. Plunk

Abstract: 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 th… ▽ More 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. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Assessing global ion thermal confinement in critical-gradient-optimized stellarators

Authors: A. Bañón Navarro, G. T. Roberg-Clark, G. G. Plunk, D. Fernando, A. Di Siena, F. Wilms, F. Jenko

Abstract: 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… ▽ More 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. △ Less

Submitted 28 October, 2023; originally announced October 2023.

The residual flow in well-optimized stellarators

Authors: G. G. Plunk, P. Helander

Abstract: 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 ste… ▽ More 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. △ Less

Submitted 11 April, 2024; v1 submitted 22 October, 2023; originally announced October 2023.

Journal ref: Journal of Plasma Physics , Volume 90 , Issue 2 , April 2024 , 905900205

Helicity of the magnetic axes of quasi-isodynamic stellarators

Authors: Katia Camacho Mata, Gabriel G. Plunk

Abstract: In this study, we explore the influence of the helicity of the magnetic axis-defined as the self-linking number of the curve-on the quality of quasi-isodynamic stellarator-symmetric configurations constructed using the near-axis expansion method (Camacho Mata et al. 2022; Plunk et al. 2019). A class of magnetic axes previously unexplored within this formalism is identified when analyzing the axis… ▽ More In this study, we explore the influence of the helicity of the magnetic axis-defined as the self-linking number of the curve-on the quality of quasi-isodynamic stellarator-symmetric configurations constructed using the near-axis expansion method (Camacho Mata et al. 2022; Plunk et al. 2019). A class of magnetic axes previously unexplored within this formalism is identified when analyzing the axis shape of the QIPC configuration (Subbotin et al. 2006): the case of half-helicity (per field period). We show these shapes are compatible with the near-axis formalism and how they can be used to construct near-axis stellarators with up-to 5 field-periods, $ε_{eff} \approx$ 1.3%, and similar rotational transform as existing conventionally optimized designs, without the need of a plasma boundary optimization. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 26 pages, 21 figures

Higher order theory of quasi-isodynamicity near the magnetic axis of stellarators

Authors: Eduardo Rodriguez, Gabe G. Plunk

Abstract: The condition of quasi-isodynamicity is derived to second order in the distance from the magnetic axis. We do so using a formulation of omnigenity that explicitly requires the balance between the radial particle drifts at opposite bounce points of a magnetic well. This is a physically intuitive alternative to the integrated condition involving distances between bounce points, used in previous work… ▽ More The condition of quasi-isodynamicity is derived to second order in the distance from the magnetic axis. We do so using a formulation of omnigenity that explicitly requires the balance between the radial particle drifts at opposite bounce points of a magnetic well. This is a physically intuitive alternative to the integrated condition involving distances between bounce points, used in previous works. We investigate the appearance of topological defects in the magnetic field strength (``puddles''). A hallmark of quasi-isodynamic fields, the curved contour of minimum field strength, is found to be inextricably linked to these defects. Our results pave the way to constructing solutions that satisfy omnigenity to a higher degree of precision, and also to simultaneously consider other physical properties, like shaping and stability. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Critical gradient turbulence optimization toward a compact stellarator reactor concept

Authors: G. T. Roberg-Clark, G. G. Plunk, P. Xanthopoulos, C. Nührenberg, S. A. Henneberg, H. M. Smith

Abstract: 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… ▽ More 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$\%$. △ Less

Submitted 6 October, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 5 pages, 5 figures. Phys. Rev. Research 5, L032030 (2023)

Energetic bounds on gyrokinetic instabilities. Part III. Generalized free energy

Authors: G. G. Plunk, P. Helander

Abstract: 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 d… ▽ More 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. △ Less

Submitted 15 September, 2023; v1 submitted 3 January, 2023; originally announced January 2023.

Journal ref: Journal of Plasma Physics, 89(4), 905890419 (2023)

Reduction of electrostatic turbulence in a quasi-helically symmetric stellarator via critical gradient optimization

Authors: G. T. Roberg-Clark, P. Xanthopoulos, G. G. Plunk

Abstract: 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. T… ▽ More 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. △ Less

Submitted 10 November, 2022; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: 11 pages, 7 figures. Submitted to the Journal of Plasma Physics

Coarse-grained gyrokinetics for the critical ion temperature gradient in stellarators

Authors: G. T. Roberg-Clark, G. G. Plunk, P. Xanthopoulos

Abstract: 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… ▽ More 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. △ Less

Submitted 11 August, 2022; originally announced August 2022.

Comments: 5 pages, 4 figures. Accepted for publication in Physical Review Research

Direct construction of stellarator-symmetric quasi-isodynamic magnetic configurations

Authors: Katia Camacho Mata, Gabriel G. Plunk, Rogerio Jorge

Abstract: We develop the formalism of the first order near-axis expansion of the MHD equilibrium equations described in Garren & Boozer (1991), Plunk et al. (2019) and Plunk et al. (2021), for the case of a quasi-isodynamic, N-field period, stellarator symmetric, single-well magnetic field equilibrium. The importance of the magnetic axis shape is investigated, and we conclude that control of the curvature a… ▽ More We develop the formalism of the first order near-axis expansion of the MHD equilibrium equations described in Garren & Boozer (1991), Plunk et al. (2019) and Plunk et al. (2021), for the case of a quasi-isodynamic, N-field period, stellarator symmetric, single-well magnetic field equilibrium. The importance of the magnetic axis shape is investigated, and we conclude that control of the curvature and torsion is crucial to obtain omnigenous configurations with finite aspect ratio and low effective ripple, especially for a higher number of field periods. For this reason a method is derived to construct classes of axis shapes with favourable curvature and torsion. Solutions are presented, including a three-field-period configuration constructed at an aspect ratio of A=20, with a maximum elongation of e=3.2 and an effective ripple under 1%, which demonstrates that high elongation is not a necessary feature of QI stellarators. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 28 pages, 19 figures

A single-field-period quasi-isodynamic stellarator

Authors: R. Jorge, G. G. Plunk, M. Drevlak, M. Landreman, J. -F. Lobsien, K. Camacho Mata, P. Helander

Abstract: 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 configuration… ▽ More 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. △ Less

Submitted 1 September, 2022; v1 submitted 11 May, 2022; originally announced May 2022.

Comments: 13 pages, 9 figures

Energetic bounds on gyrokinetic instabilities. Part II. Modes of optimal growth

Authors: G. G. Plunk, Per Helander

Abstract: 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, efficien… ▽ More 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. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Energetic bounds on gyrokinetic instabilities. Part I: Fundamentals

Authors: P. Helander, G. G. Plunk

Abstract: 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-ampli… ▽ More 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. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Predicting the Z-pinch Dimits shift through gyrokinetic tertiary instability analysis of the entropy mode

Authors: A. Hallenbert, G. G. Plunk

Abstract: 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 com… ▽ More 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. △ Less

Submitted 15 November, 2021; originally announced November 2021.

Comments: 36 pages, 15 figures

Upper bounds on gyrokinetic instabilities

Authors: P. Helander, G. G. Plunk

Abstract: 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… ▽ More 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. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Predicting the Dimits shift through reduced mode tertiary instability analysis in a strongly driven gyrokinetic fluid limit

Authors: A. Hallenbert, G. G. Plunk

Abstract: 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, beyon… ▽ More 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. △ Less

Submitted 9 August, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 30 pages, 9 figures; improved the derivation of Appendix A, slightly modified several statements for accuracy, updated figures and notation

Calculating the linear critical gradient for the ion-temperature-gradient mode in magnetically confined plasmas

Authors: G. T. Roberg-Clark, G. G. Plunk, P. Xanthopoulos

Abstract: 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… ▽ More 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. △ Less

Submitted 27 April, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Comments: 22 pages, 6 figures. Accepted for publication in The Journal of Plasma Physics (JPP)

Journal ref: Journal of Plasma Physics, Volume 87, Issue 3, June 2021, 905870306

Exploring zonal flow mediated saturation on stellarators

Authors: C. D. Mora Moreno, J. H. E. Proll, G. G. Plunk, P. Xanthopoulos

Abstract: 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… ▽ More 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. △ Less

Submitted 30 September, 2020; originally announced September 2020.

Comments: 11 pages, 9 figures

Enstrophy non-conservation and the forward cascade of energy in two-dimensional electrostatic magnetized plasma turbulence

Authors: G. G. Plunk

Abstract: 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 t… ▽ More 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]. △ Less

Submitted 7 August, 2020; originally announced August 2020.

Perturbing an axisymmetric magnetic equilibrium to obtain a quasi-axisymmetric stellarator

Authors: G. G. Plunk

Abstract: 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 pertur… ▽ More 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. △ Less

Submitted 20 August, 2020; v1 submitted 6 May, 2020; originally announced May 2020.

Global gyrokinetic simulations of ITG turbulence in the configuration space of the Wendelstein 7-X stellarator

Authors: A. Bañón Navarro, G. Merlo, G. G. Plunk, P. Xanthopoulos, A. von Stechow, A. Di Siena, M. Maurer, F. Hindenlang, F. Wilms, F. Jenko

Abstract: 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 mor… ▽ More 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. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Suppression of turbulence by trapped electrons in optimised stellarators

Authors: J. H. E. Proll, P. Xanthopoulos, P. Helander, G. G. Plunk, B. J. Faber, T. Görler, H. M. Smith, M. J. Pueschel

Abstract: 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-d… ▽ More 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. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Direct construction of optimized stellarator shapes. III. Omnigenity near the magnetic axis

Authors: G. G. Plunk, M. Landreman, P. Helander

Abstract: 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 provid… ▽ More 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. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Direct construction of optimized stellarator shapes. II. Numerical quasisymmetric solutions

Authors: Matt Landreman, Wrick Sengupta, Gabriel G Plunk

Abstract: Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to minimize symmetry-breaking Fourier modes of the field stre… ▽ More Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding center particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to minimize symmetry-breaking Fourier modes of the field strength $B$. Here instead we directly construct magnetic fields in cylindrical coordinates that are quasisymmetric to leading order in distance from the magnetic axis, without using optimization. The method involves solution of a 1-dimensional nonlinear ordinary differential equation, originally derived by Garren and Boozer [Phys. Fluids B 3, 2805 (1991)]. We demonstrate the usefulness and accuracy of this optimization-free approach by providing the results of this construction as input to the codes VMEC and BOOZ_XFORM, confirming the purity and scaling of the magnetic spectrum. The space of magnetic fields that are quasisymmetric to this order is parameterized by the magnetic axis shape along with three other real numbers, one of which reflects the on-axis toroidal current density, and another one of which is zero for stellarator symmetry. The method here could be used to generate good initial conditions for conventional optimization, and its speed enables exhaustive searches of parameter space. △ Less

Submitted 3 December, 2018; v1 submitted 26 September, 2018; originally announced September 2018.

Quasi-axisymmetric magnetic fields: weakly non-axisymmetric case in a vacuum

Authors: G. G. Plunk, P. Helander

Abstract: 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… ▽ More 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. △ Less

Submitted 8 March, 2018; v1 submitted 9 January, 2018; originally announced January 2018.

Collisionless microinstabilities in stellarators. IV. The ion-driven trapped-electron mode

Authors: G. G. Plunk, J. W. Connor, P. Helander

Abstract: 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 wi… ▽ More 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. △ Less

Submitted 14 August, 2017; originally announced August 2017.

Comments: Accepted to Journal of Plasma Physics

On the scaling of ion and electron temperature gradient driven turbulence in slab geometry

Authors: G. G. Plunk

Abstract: 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. 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. △ Less

Submitted 9 March, 2017; originally announced March 2017.

Comments: To be published, Phys. Plasmas (2017)

Distinct turbulence saturation regimes in stellarators

Authors: G. G. Plunk, P. Xanthopoulos, P. Helander

Abstract: 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 quas… ▽ More 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. △ Less

Submitted 29 September, 2017; v1 submitted 9 March, 2017; originally announced March 2017.

Comments: Erratum added to end

Journal ref: Phys. Rev. Lett. 118, 105002 (2017)

Nonlinear growth of zonal flows by secondary instability in general magnetic geometry

Authors: G. G. Plunk, A. Bañón Navarro

Abstract: 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 paramete… ▽ More 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. △ Less

Submitted 18 February, 2017; originally announced February 2017.

Comments: New J. Phys. 2017

Resilience of quasi-isodynamic stellarators against trapped-particle instabilities

Authors: J. H. E. Proll, P. Helander, J. W. Connor, G. G. Plunk

Abstract: 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 ba… ▽ More 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. △ Less

Submitted 14 September, 2015; originally announced September 2015.

Journal ref: Phys. Rev. Lett. 108, 245002 (2012)

The universal instability in general geometry

Authors: P. Helander, G. G. Plunk

Abstract: 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 tra… ▽ More 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. △ Less

Submitted 30 June, 2015; originally announced June 2015.

On the nonlinear stability of a quasi-two-dimensional drift kinetic model for ion temperature gradient turbulence

Authors: G. G. Plunk

Abstract: 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 whe… ▽ More 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. △ Less

Submitted 15 April, 2015; v1 submitted 21 January, 2015; originally announced January 2015.

Generalized universal instability: Transient linear amplification and subcritical turbulence

Authors: Matt Landreman, Gabriel G. Plunk, William Dorland

Abstract: In this work we numerically demonstrate both significant transient (i.e. non-modal) linear amplification and sustained nonlinear turbulence in a kinetic plasma system with no unstable eigenmodes. The particular system considered is an electrostatic plasma slab with magnetic shear, kinetic electrons and ions, weak collisions, and a density gradient, but with no temperature gradient. In contrast to… ▽ More In this work we numerically demonstrate both significant transient (i.e. non-modal) linear amplification and sustained nonlinear turbulence in a kinetic plasma system with no unstable eigenmodes. The particular system considered is an electrostatic plasma slab with magnetic shear, kinetic electrons and ions, weak collisions, and a density gradient, but with no temperature gradient. In contrast to hydrodynamic examples of non-modal growth and subcritical turbulence, here there is no sheared flow in the equilibrium. Significant transient linear amplification is found when the magnetic shear and collisionality are weak. It is also demonstrated that nonlinear turbulence can be sustained if initialized at sufficient amplitude. We prove these two phenomena are related: when sustained turbulence occurs without unstable eigenmodes, states that are typical of the turbulence must yield transient linear amplification of the gyrokinetic free energy. △ Less

Submitted 13 January, 2015; originally announced January 2015.

Journal ref: J. Plasma Phys. 81 (2015) 905810501

Irreversible energy flow in forced Vlasov dynamics

Authors: G. G. Plunk, J. T. Parker

Abstract: 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… ▽ More 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. △ Less

Submitted 7 January, 2015; v1 submitted 28 February, 2014; originally announced February 2014.

Comments: Small changes for clarity

Journal ref: Eur. Phys. J. D (2014) 68: 296

Understanding nonlinear saturation in zonal-flow-dominated ion temperature gradient turbulence

Authors: G. G. Plunk, A. Bañón Navarro, F. Jenko

Abstract: 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 assum… ▽ More 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. △ Less

Submitted 6 March, 2015; v1 submitted 31 January, 2014; originally announced January 2014.

Comments: Final journal version. Some clarifications and a new Fig. 4

Journal ref: Plasma Phys. Control. Fusion 57 045005 (2015)

Collisionless microinstabilities in stellarators III - the ion-temperature-gradient mode

Authors: G. G. Plunk, P. Helander, P. Xanthopoulos, J. W. Connor

Abstract: 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… ▽ More 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). △ Less

Submitted 25 March, 2014; v1 submitted 9 December, 2013; originally announced December 2013.

Comments: v3: final version for journal

Journal ref: Phys. Plasmas 21, 032112 (2014)

Collisionless microinstabilities in stellarators I - analytical theory of trapped-particle modes

Authors: P. Helander, J. H. E. Proll, G. G. Plunk

Abstract: 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 a… ▽ More 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. △ Less

Submitted 14 November, 2013; v1 submitted 13 November, 2013; originally announced November 2013.

Comments: 9 pages, 1 figure

Multiscale Gyrokinetics for Rotating Tokamak Plasmas: Fluctuations, Transport and Energy Flows

Authors: I. G. Abel, G. G. Plunk, E. Wang, M. Barnes, S. C. Cowley, W. Dorland, A. A. Schekochihin

Abstract: 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 e… ▽ More 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. △ Less

Submitted 2 December, 2013; v1 submitted 21 September, 2012; originally announced September 2012.

Comments: 113 pages, 3 figures

Journal ref: Rep. Prog. Phys. 76 116201, (2013)

Freely decaying turbulence in two-dimensional electrostatic gyrokinetics

Authors: T. Tatsuno, G. G. Plunk, M. Barnes, W. Dorland, G. G. Howes, R. Numata

Abstract: 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 th… ▽ More 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*. △ Less

Submitted 11 December, 2012; v1 submitted 7 August, 2012; originally announced August 2012.

Comments: 12 pages, 12 figures; replaced to match published version

Journal ref: Phys. Plasmas 19, 122305 (2012)

Landau Damping in a Turbulent Setting

Authors: G. G. Plunk

Abstract: 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 m… ▽ More 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. △ Less

Submitted 3 December, 2012; v1 submitted 15 June, 2012; originally announced June 2012.

Comments: Introduction significantly expanded to help contextualize results. Calculations unchanged

Considering Fluctuation Energy as a Measure of Gyrokinetic Turbulence

Authors: G. G. Plunk, T. Tatsuno, W. Dorland

Abstract: 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… ▽ More 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. △ Less

Submitted 18 October, 2012; v1 submitted 4 June, 2012; originally announced June 2012.

Comments: Final version as published. Some cosmetic changes and update of references

Journal ref: New J. Phys. 14 (2012) 103030

Kuramoto model with coupling through an external medium

Authors: David J. Schwab, Gabriel G. Plunk, Pankaj Mehta

Abstract: Synchronization of coupled oscillators is often described using the Kuramoto model. Here we study a generalization of the Kuramoto model where oscillators communicate with each other through an external medium. This generalized model exhibits interesting new phenomena such as bistability between synchronization and incoherence and a qualitatively new form of synchronization where the external medi… ▽ More Synchronization of coupled oscillators is often described using the Kuramoto model. Here we study a generalization of the Kuramoto model where oscillators communicate with each other through an external medium. This generalized model exhibits interesting new phenomena such as bistability between synchronization and incoherence and a qualitatively new form of synchronization where the external medium exhibits small-amplitude oscillations. We conclude by discussing the relationship of the model to other variations of the Kuramoto model including the Kuramoto model with a bimodal frequency distribution and the Millennium Bridge problem. △ Less

Submitted 13 December, 2011; originally announced December 2011.

Comments: 9 pages, 3 figures

Response: Comment on "Energy Transfer and Dual Cascade in Kinetic Magnetized Plasma Turbulence"

Authors: G. G. Plunk, T. Tatsuno

Abstract: 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. Othe… ▽ More 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. △ Less

Submitted 5 June, 2011; originally announced June 2011.

Energy transfer and dual cascade in kinetic magnetized plasma turbulence

Authors: G. G. Plunk, T. Tatsuno

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 15 April, 2011; v1 submitted 27 July, 2010; originally announced July 2010.

Comments: Reduced size for journal format

Gyrokinetic simulation of entropy cascade in two-dimensional electrostatic turbulence

Authors: T. Tatsuno, M. Barnes, S. C. Cowley, W. Dorland, G. G. Howes, R. Numata, G. G. Plunk, A. A. Schekochihin

Abstract: 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 flu… ▽ More 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. △ Less

Submitted 16 August, 2010; v1 submitted 20 March, 2010; originally announced March 2010.

Comments: 8 pages, 10 figures, Conference paper presented at 2009 Asia-Pacific Plasma Theory Conference. Ver.2 includes corrected typos & updated references

Journal ref: J. Plasma Fusion Res. SERIES 9, 509 (2010)

Two-dimensional gyrokinetic turbulence

Authors: G. G. Plunk, S. C. Cowley, A. A. Schekochihin, T. Tatsuno

Abstract: 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 tw… ▽ More 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. △ Less

Submitted 10 July, 2010; v1 submitted 1 April, 2009; originally announced April 2009.

Comments: 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

The theory of gyrokinetic turbulence: A multiple-scales approach

Authors: Gabriel G. Plunk

Abstract: Gyrokinetics is a rich and rewarding playground to study some of the mysteries of modern physics. In this thesis I present work, motivated by the quest for fusion energy, which seeks to uncover some of the inner workings of turbulence in magnetized plasmas. I begin with the fundamental theory of gyrokinetics, and a novel formulation of its extension to the equations for mean-scale transport --… ▽ More Gyrokinetics is a rich and rewarding playground to study some of the mysteries of modern physics. In this thesis I present work, motivated by the quest for fusion energy, which seeks to uncover some of the inner workings of turbulence in magnetized plasmas. I begin with the fundamental theory of gyrokinetics, and a novel formulation of its extension to the equations for mean-scale transport -- the equations which must be solved to determine the performance of magnetically confined fusion devices. The second project presents gyrokinetic secondary instability theory as a mechanism to bring about saturation of the basic instabilities that drive gyrokinetic turbulence. Emphasis is put on the ability for this analytic theory to predict basic properties of the nonlinear state, which can be applied to a mixing length phenomenology of transport. The final project is an application of the methods from inertial range understanding of fluid turbulence, to describe the stationary state of fully developed two-dimensional gyrokinetic turbulence. This work explores the relatively new idea of a phase-space cascade, whereby fine scales are nonlinearly generated in both position space and velocity space, and ultimately smoothed by collisional entropy production. This process constitutes the thermodynamic balance which occurs in the true steady state of a turbulent plasma, including those found in fusion devices. △ Less

Submitted 23 April, 2009; v1 submitted 5 March, 2009; originally announced March 2009.

Comments: The PhD dissertation of Gabriel Galad Plunk, UCLA (2009)

Nonlinear phase mixing and phase-space cascade of entropy in gyrokinetic plasma turbulence

Authors: T. Tatsuno, W. Dorland, A. A. Schekochihin, G. G. Plunk, M. Barnes, S. C. Cowley, G. G. Howes

Abstract: 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 i… ▽ More 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. △ Less

Submitted 30 June, 2009; v1 submitted 16 November, 2008; originally announced November 2008.

Comments: Revtex, 4 pages, 3 figures; replaced to match published version

Journal ref: Phys Rev Lett 103, 015003 (2009)