Source author record

G. G. Plunk

G. G. Plunk appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

Energetic bounds on gyrokinetic instabilities. Part I: Fundamentals

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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].

preprint2020arXiv

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

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.

preprint2020arXiv

Suppression of turbulence by trapped electrons in optimised stellarators

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.

preprint2019arXiv

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

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.

preprint2015arXiv

Irreversible energy flow in forced Vlasov dynamics

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Resilience of quasi-isodynamic stellarators against trapped-particle instabilities

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.

preprint2015arXiv

The universal instability in general 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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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).

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Considering Fluctuation Energy as a Measure of Gyrokinetic Turbulence

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.

preprint2012arXiv

Freely decaying turbulence in two-dimensional electrostatic gyrokinetics

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*.

preprint2012arXiv

Landau Damping in a Turbulent Setting

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.

preprint2011arXiv

Energy transfer and dual cascade in kinetic magnetized plasma turbulence

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.

preprint2011arXiv

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

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.

preprint2010arXiv

Gyrokinetic simulation of entropy cascade in two-dimensional electrostatic turbulence

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.

preprint2010arXiv

Two-dimensional gyrokinetic turbulence

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.

preprint2009arXiv

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

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.