Source author record

William Dorland

William 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

12works
7topics
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

12 published item(s)

preprint2026arXiv

Hidden free energy released by explicit parity-time-symmetry breaking

It is shown that the familiar two-stream instability is the result of spontaneous parity-time (PT)-symmetry breaking in a conservative system, and more importantly, explicit PT-symmetry breaking by viscosity can destabilize the system in certain parameter regimes that are stable when viscosity vanishes. This reveals that complex systems may possess hidden free energies protected by PT-symmetry and viscosity, albeit dissipative, can expose the systems to these freed energies by breaking PT-symmetry explicitly. Such a process is accompanied by instability and total variation growth.

preprint2016arXiv

Fundamental form of the electrostatic $δf$-PIC algorithm and discovery of a converged numerical instability

The $δf$ particle-in-cell algorithm has been a useful tool in studying the physics of plasmas, particularly turbulent magnetized plasmas in the context of gyrokinetics. The reduction in noise due to not having to resolve the full distribution function indicates an efficiency advantage over standard ("full-$f$") particle-in-cell. Despite its successes, the algorithm behaves strangely in some circumstances. In this work, we document a fully-resolved numerical instability that occurs in the simplest of multiple-species test cases: the electrostatic $Ω_H$ mode. There is also a poorly-understood numerical instability that occurs when one is under-resolved in particle number, which may require a prohibitively large number of particles to stabilize. Both of these are independent of the time-stepping scheme, and we conclude that they exist if the time advancement were exact. The exact analytic form of the algorithm is presented, and several schemes for mitigating these instabilities are presented.

preprint2016arXiv

Transport and deceleration of fusion products in microturbulence

The velocity-space distribution of alpha particles born in fusion devices is subject to modification at moderate energies due to turbulent transport. Therefore, one must calculate the evolution of an equilibrium distribution whose functional form is not known \emph{a priori}. Using a novel technique, applicable to any trace impurity, we have made this calculation not only possible, but particularly efficient. We demonstrate a microturbulence-induced departure from the local slowing-down distribution, an inversion of the energy distribution, and associated modifications to the alpha heating and pressure profiles in an ITER-like scenario.

preprint2015arXiv

Generalized universal instability: Transient linear amplification and subcritical turbulence

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.

preprint2015arXiv

Universal instability for wavelengths below the ion Larmor scale

We demonstrate that the universal mode driven by the density gradient in a plasma slab can be absolutely unstable even in the presence of reasonable magnetic shear. Previous studies from the 1970s that reached the opposite conclusion used an eigenmode equation limited to $L_x \gg ρ_i$, where $L_x$ is the scale length of the mode in the radial direction, and $ρ_i$ is the ion Larmor radius. Here we instead use a gyrokinetic approach which does not have this same limitation. Instability is found for perpendicular wavenumbers $k_y$ in the range $0.7 \lesssim k_y ρ_i \lesssim 100$, and for sufficiently weak magnetic shear: $L_s / L_n \ge 17$, where $L_s$ and $L_n$ are the scale lengths of magnetic shear and density. Thus, the plasma drift wave in a sheared magnetic field may be unstable even with no temperature gradients, no trapped particles, and no magnetic curvature.

preprint2014arXiv

Intrinsic momentum transport in up-down asymmetric tokamaks

Recent work demonstrated that breaking the up-down symmetry of tokamak flux surfaces removes a constraint that limits intrinsic momentum transport, and hence toroidal rotation, to be small. We show, through MHD analysis, that ellipticity is most effective at introducing up-down asymmetry throughout the plasma. We detail an extension to GS2, a local $δf$ gyrokinetic code that self-consistently calculates momentum transport, to permit up-down asymmetric configurations. Tokamaks with tilted elliptical poloidal cross-sections were simulated to determine nonlinear momentum transport. The results, which are consistent with experiment in magnitude, suggest that a toroidal velocity gradient, $(\partial u_{ζi} / \partial ρ) / v_{th i}$, of 5% of the temperature gradient, $(\partial T_{i} / \partial ρ) / T_{i}$, is sustainable. Here $v_{th i}$ is the ion thermal speed, $u_{ζi}$ is the ion toroidal mean flow, $ρ$ is the minor radial coordinate normalized to the tokamak minor radius, and $T_{i}$ is the ion temperature. Since other intrinsic momentum transport mechanisms scale poorly to larger machines, these results indicate that up-down asymmetry is the most feasible method to generate the current experimentally-measured rotation levels in reactor-sized devices.

preprint2014arXiv

Validating modelling assumptions of alpha particles in electrostatic turbulence

To rigorously model fast ions in fusion plasmas, a non-Maxwellian equilibrium distribution must be used. In the work, the response of high-energy alpha particles to electrostatic turbulence has been analyzed for several different tokamak parameters. Our results are consistent with known scalings and experimental evidence that alpha particles are generally well-confined: on the order of several seconds. It is also confirmed that the effect of alphas on the turbulence is negligible at realistically low concentrations, consistent with linear theory. It is demonstrated that the usual practice of using a high-temperature Maxwellian gives incorrect estimates for the radial alpha particle flux, and a method of correcting it is provided. Furthermore, we see that the timescales associated with collisions and transport compete at moderate energies, calling into question the assumption that alpha particles remain confined to a flux surface that is used in the derivation of the slowing-down distribution.

preprint2013arXiv

Alfven Wave Collisions, The Fundamental Building Block of Plasma Turbulence II: Numerical Solution

This paper presents the numerical verification of an asymptotic analytical solution for the nonlinear interaction between counterpropagating Alfven waves, the fundamental building block of astrophysical plasma turbulence. The analytical solution, derived in the weak turbulence limit using the equations of incompressible MHD, is compared to a nonlinear gyrokinetic simulation of an Alfven wave collision. The agreement between these methods signifies that the incompressible solution satisfactorily describes the essential dynamics of the nonlinear energy transfer, even under the weakly collisional plasma conditions relevant to many astrophysical environments.

preprint2011arXiv

Gyrokinetic Simulations of Solar Wind Turbulence from Ion to Electron Scales

The first three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of $k^{-2.8}$ as observed in \emph{in situ} spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic \Alfven wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at sub-ion-Larmor radius scales, which provides the first evidence of the ion entropy cascade in an electromagnetic turbulence simulation.

preprint2010arXiv

AstroGK: Astrophysical Gyrokinetics Code

The gyrokinetic simulation code AstroGK is developed to study fundamental aspects of kinetic plasmas and for applications mainly to astrophysical problems. AstroGK is an Eulerian slab code that solves the electromagnetic Gyrokinetic-Maxwell equations in five-dimensional phase space, and is derived from the existing gyrokinetics code GS2 by removing magnetic geometry effects. Algorithms used in the code are described. The code is benchmarked using linear and nonlinear problems. Serial and parallel performance scalings are also presented.

preprint2009arXiv

Resolving velocity space dynamics in continuum gyrokinetics

Many plasmas of interest to the astrophysical and fusion communities are weakly collisional. In such plasmas, small scales can develop in the distribution of particle velocities, potentially affecting observable quantities such as turbulent fluxes. Consequently, it is necessary to monitor velocity space resolution in gyrokinetic simulations. In this paper, we present a set of computationally efficient diagnostics for measuring velocity space resolution in gyrokinetic simulations and apply them to a range of plasma physics phenomena using the continuum gyrokinetic code GS2. For the cases considered here, it is found that the use of a collisionality at or below experimental values allows for the resolution of plasma dynamics with relatively few velocity space grid points. Additionally, we describe implementation of an adaptive collision frequency which can be used to improve velocity space resolution in the collisionless regime, where results are expected to be independent of collision frequency.

preprint2008arXiv

A Model of Turbulence in Magnetized Plasmas: Implications for the Dissipation Range in the Solar Wind

This paper studies the turbulent cascade of magnetic energy in weakly collisional magnetized plasmas. A cascade model is presented, based on the assumptions of local nonlinear energy transfer in wavenumber space, critical balance between linear propagation and nonlinear interaction times, and the applicability of linear dissipation rates for the nonlinearly turbulent plasma. The model follows the nonlinear cascade of energy from the driving scale in the MHD regime, through the transition at the ion Larmor radius into the kinetic Alfven wave regime, in which the turbulence is dissipated by kinetic processes. The turbulent fluctuations remain at frequencies below the ion cyclotron frequency due to the strong anisotropy of the turbulent fluctuations, k_parallel << k_perp (implied by critical balance). In this limit, the turbulence is optimally described by gyrokinetics; it is shown that the gyrokinetic approximation is well satisfied for typical slow solar wind parameters. Wave phase velocity measurements are consistent with a kinetic Alfven wave cascade and not the onset of ion cyclotron damping. The conditions under which the gyrokinetic cascade reaches the ion cyclotron frequency are established. Cascade model solutions imply that collisionless damping provides a natural explanation for the observed range of spectral indices in the dissipation range of the solar wind. The dissipation range spectrum is predicted to be an exponential fall off; the power-law behavior apparent in observations may be an artifact of limited instrumental sensitivity. The cascade model is motivated by a programme of gyrokinetic simulations of turbulence and particle heating in the solar wind.