Source author record

W. Dorland

W. Dorland 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

20works
8topics
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

20 published item(s)

preprint2022arXiv

Energy partition between Alfvénic and compressive fluctuations in magnetorotational turbulence with near-azimuthal mean magnetic field

The theory of magnetohydrodynamic (MHD) turbulence predicts that Alfvénic and slow-mode-like compressive fluctuations are energetically decoupled at small scales in the inertial range. The partition of energy between these fluctuations determines the nature of dissipation, which, in many astrophysical systems, happens on scales where plasma is collisionless. However, when the magnetorotational instability (MRI) drives the turbulence, it is difficult to resolve numerically the scale at which both types of fluctuations start to be decoupled because the MRI energy injection occurs in a broad range of wavenumbers, and both types of fluctuations are usually expected to be coupled even at relatively small scales. In this study, we focus on collisional MRI turbulence threaded by a near-azimuthal mean magnetic field, which is naturally produced by the differential rotation of a disc. We show that, in such a case, the decoupling scales are reachable using a reduced MHD model that includes differential-rotation effects. In our reduced MHD model, the Alfvénic and compressive fluctuations are coupled only through the linear terms that are proportional to the angular velocity of the accretion disc. We numerically solve for the turbulence in this model and show that the Alfvénic and compressive fluctuations are decoupled at the small scales of our simulations as the nonlinear energy transfer dominates the linear coupling below the MRI-injection scale. In the decoupling scales, the energy flux of compressive fluctuations contained in the small scales is almost double that of Alfvénic fluctuations. Finally, we discuss the application of this result to prescriptions of ion-to-electron heating ratio in hot accretion flows.

preprint2022arXiv

Three-Dimensional Inhomogeneity of Electron-Temperature-Gradient Turbulence in the Edge of Tokamak Plasmas

Nonlinear multiscale gyrokinetic simulations of a Joint European Torus edge pedestal are used to show that electron-temperature-gradient (ETG) turbulence has a rich three-dimensional structure, varying strongly according to the local magnetic-field configuration. In the plane normal to the magnetic field, the steep pedestal electron temperature gradient gives rise to anisotropic turbulence with a radial (normal) wavelength much shorter than in the binormal direction. In the parallel direction, the location and parallel extent of the turbulence are determined by the variation in the magnetic drifts and finite-Larmor-radius (FLR) effects. The magnetic drift and FLR topographies have a perpendicular-wavelength dependence, which permits turbulence intensity maxima near the flux-surface top and bottom at longer binormal scales, but constrains turbulence to the outboard midplane at shorter electron-gyroradius binormal scales. Our simulations show that long-wavelength ETG turbulence does not transport heat efficiently, and significantly decreases overall ETG transport -- in our case by $\sim$40 \% -- through multiscale interactions.

preprint2020arXiv

Ion versus electron heating in compressively driven astrophysical gyrokinetic turbulence

The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear hybrid gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_\rmi/Q_\rme$ as a function of the compressive-to-Alfvénic driving power ratio $P_\compr/P_\AW$, of the ratio of ion thermal pressure to magnetic pressure $β_\rmi$, and of the ratio of ion-to-electron background temperatures $T_\rmi/T_\rme$. It is shown that $Q_\rmi/Q_\rme$ is an increasing function of $P_\compr/P_\AW$. When the compressive driving is sufficiently large, $Q_\rmi/Q_\rme$ approaches $\simeq P_\compr/P_\AW$. This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfvén waves is absent for both low and high $β_\rmi$, meaning that the compressive driving is directly connected to the ion entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low $β_\rmi$ and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

preprint2020arXiv

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$-pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal-temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations. A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear. The break up of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system -- the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol{E}\times\boldsymbol{B}$ flow. By analysing the linear ITG modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

preprint2016arXiv

Phase mixing vs. nonlinear advection in drift-kinetic plasma turbulence

A scaling theory of long-wavelength electrostatic turbulence in a magnetised, weakly collisional plasma (e.g., ITG turbulence) is proposed, with account taken both of the nonlinear advection of the perturbed particle distribution by fluctuating ExB flows and of its phase mixing, which is caused by the streaming of the particles along the mean magnetic field and, in a linear problem, would lead to Landau damping. It is found that it is possible to construct a consistent theory in which very little free energy leaks into high velocity moments of the distribution function, rendering the turbulent cascade in the energetically relevant part of the wave-number space essentially fluid-like. The velocity-space spectra of free energy expressed in terms of Hermite-moment orders are steep power laws and so the free-energy content of the phase space does not diverge at infinitesimal collisionality (while it does for a linear problem); collisional heating due to long-wavelength perturbations vanishes in this limit (also in contrast with the linear problem, in which it occurs at the finite rate equal to the Landau-damping rate). The ability of the free energy to stay in the low velocity moments of the distribution function is facilitated by the "anti-phase-mixing" effect, whose presence in the nonlinear system is due to the stochastic version of the plasma echo (the advecting velocity couples the phase-mixing and anti-phase-mixing perturbations). The partitioning of the wave-number space between the (energetically dominant) region where this is the case and the region where linear phase mixing wins its competition with nonlinear advection is governed by the "critical balance" between linear and nonlinear timescales (which for high Hermite moments splits into two thresholds, one demarcating the wave-number region where phase mixing predominates, the other where plasma echo does).

preprint2014arXiv

Collisionless Reconnection in the Large Guide Field Regime: Gyrokinetic Versus Particle-in-Cell Simulations

Results of the first validation of large guide field, $B_g / δB_0 \gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $β_i = 0.01$ and solar wind relevant $β_i = 1$ are presented, where $δB_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $β_i B_g / δB_0 \gtrsim 1$ and $B_g / δB_0 \gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.

preprint2013arXiv

An Oscillating Langevin Antenna for Driving Plasma Turbulence Simulations

A unique method of driving Alfvenic turbulence via an oscillating Langevin antenna is presented. This method of driving is motivated by a desire to inject energy into a finite domain numerical simulation in a manner that models the nonlinear transfer of energy from fluctuations in the turbulent cascade at scales larger than the simulation domain.. The oscillating Langevin antenna is shown to capture the essential features of the larger scale turbulence and efficiently couple to the plasma, generating steady-state turbulence within one characteristic turnaround time. The antenna is also sufficiently flexible to explore both strong and weak regimes of Alfvenic plasma turbulence.

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

Turbulent transport and heating of trace heavy ions in hot, magnetized plasmas

Scaling laws for the transport and heating of trace heavy ions in low-frequency, magnetized plasma turbulence are derived and compared with direct numerical simulations. The predicted dependences of turbulent fluxes and heating on ion charge and mass number are found to agree with numerical results for both stationary and differentially rotating plasmas. Heavy ion momentum transport is found to increase with mass, and heavy ions are found to be preferentially heated, implying a mass-dependent ion temperature for very weakly collisional plasmas and for partially-ionized heavy ions in strongly rotating plasmas.

preprint2012arXiv

Zero-Turbulence Manifold in a Toroidal Plasma

Sheared toroidal flows can cause bifurcations to zero-turbulent-transport states in tokamak plasmas. The maximum temperature gradients that can be reached are limited by subcritical turbulence driven by the parallel velocity gradient. Here it is shown that q/ε(magnetic field pitch/inverse aspect ratio) is a critical control parameter for sheared tokamak turbulence. By reducing q/ε, far higher temperature gradients can be achieved without triggering turbulence, in some instances comparable to those found experimentally in transport barriers. The zero-turbulence manifold is mapped out, in the zero-magnetic-shear limit, over the parameter space (γ_E, q/ε, R/L_T), where γ_E is the perpendicular flow shear and R/L_T is the normalised inverse temperature gradient scale. The extent to which it can be constructed from linear theory is discussed.

preprint2011arXiv

A Weakened Cascade Model for Turbulence in Astrophysical Plasmas

A refined cascade model for kinetic turbulence in weakly collisional astrophysical plasmas is presented that includes both the transition between weak and strong turbulence and the effect of nonlocal interactions on the nonlinear transfer of energy. The model describes the transition between weak and strong MHD turbulence and the complementary transition from strong kinetic Alfven wave (KAW) turbulence to weak dissipating KAW turbulence, a new regime of weak turbulence in which the effects of shearing by large scale motions and kinetic dissipation play an important role. The inclusion of the effect of nonlocal motions on the nonlinear energy cascade rate in the dissipation range, specifically the shearing by large-scale motions, is proposed to explain the nearly power-law energy spectra observed in the dissipation range of both kinetic numerical simulations and solar wind observations.

preprint2011arXiv

Gyrokinetic simulations of the tearing instability

Linear gyrokinetic simulations covering the collisional -- collisionless transitional regime of the tearing instability are performed. It is shown that the growth rate scaling with collisionality agrees well with that predicted by a two-fluid theory for a low plasma beta case in which ion kinetic dynamics are negligible. Electron wave-particle interactions (Landau damping), finite Larmor radius, and other kinetic effects invalidate the fluid theory in the collisionless regime, in which a general non-polytropic equation of state for pressure (temperature) perturbations should be considered. We also vary the ratio of the background ion to electron temperatures, and show that the scalings expected from existing calculations can be recovered, but only in the limit of very low beta.

preprint2011arXiv

Simulating Gyrokinetic Microinstabilities in Stellarator Geometry with GS2

The nonlinear gyrokinetic code GS2 has been extended to treat non-axisymmetric stellarator geometry. Electromagnetic perturbations and multiple trapped particle regions are allowed. Here, linear, collisionless, electrostatic simulations of the quasi-axisymmetric, three-field period National Compact Stellarator Experiment (NCSX) design QAS3-C82 have been successfully benchmarked against the eigenvalue code FULL. Quantitatively, the linear stability calculations of GS2 and FULL agree to within ~10%.

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.

preprint2009arXiv

Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas

We present a theoretical framework for plasma turbulence in astrophysical plasmas (solar wind, interstellar medium, galaxy clusters, accretion disks). The key assumptions are that the turbulence is anisotropic with respect to the mean magnetic field and frequencies are low compared to the ion cyclotron frequency. The energy injected at the outer scale scale has to be converted into heat, which ultimately cannot be done without collisions. A KINETIC CASCADE develops that brings the energy to collisional scales both in space and velocity. Its nature depends on the physics of plasma fluctuations. In each of the physically distinct scale ranges, the kinetic problem is systematically reduced to a more tractable set of equations. In the "inertial range" above the ion gyroscale, the kinetic cascade splits into a cascade of Alfvenic fluctuations, which are governed by the RMHD equations at both the collisional and collisionless scales, and a passive cascade of compressive fluctuations, which obey a linear kinetic equation along the moving field lines associated with the Alfvenic component. In the "dissipation range" between the ion and electron gyroscales, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like Electron RMHD equations and a passive phase-space cascade of ion entropy fluctuations. The latter cascade brings the energy of the inertial-range fluctuations that was damped by collisionless wave-particle interaction at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for these cascades. Astrophysical and space-physical applications are discussed in detail.

preprint2009arXiv

Direct multiscale coupling of a transport code to gyrokinetic turbulence codes

Direct coupling between a transport solver and local, nonlinear gyrokinetic calculations using the multiscale gyrokinetic code TRINITY [M. Barnes, Ph.D. thesis, arxiv:0901.2868] is described. The coupling of the microscopic and macroscopic physics is done within the framework of multiscale gyrokinetic theory, of which we present the assumptions and key results. An assumption of scale separation in space and time allows for the simulation of turbulence in small regions of the space-time grid, which are embedded in a coarse grid on which the transport equations are implicitly evolved. This leads to a reduction in computational expense of several orders of magnitude, making first-principles simulations of the full fusion device volume over the confinement time feasible on current computing resources. Numerical results from TRINITY simulations are presented and compared with experimental data from JET and ASDEX Upgrade plasmas.

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.

preprint2008arXiv

Linearized model Fokker-Planck collision operators for gyrokinetic simulations. II. Numerical implementation and tests

A set of key properties for an ideal dissipation scheme in gyrokinetic simulations is proposed, and implementation of a model collision operator satisfying these properties is described. This operator is based on the exact linearized test-particle collision operator, with approximations to the field-particle terms that preserve conservation laws and an H-Theorem. It includes energy diffusion, pitch-angle scattering, and finite Larmor radius effects corresponding to classical (real-space) diffusion. The numerical implementation in the continuum gyrokinetic code GS2 is fully implicit and guarantees exact satisfaction of conservation properties. Numerical results are presented showing that the correct physics is captured over the entire range of collisionalities, from the collisionless to the strongly collisional regimes, without recourse to artificial dissipation.