BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

M. Barnes

M. Barnes contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.plasm-phastro-ph.HEphysics.comp-phphysics.space-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

8 published item(s)

preprint2022arXiv

A novel approach to radially global gyrokinetic simulation using the flux-tube code $\texttt{stella}$

A novel approach to global gyrokinetic simulation is implemented in the flux-tube code $\texttt{stella}$. This is done by using a subsidiary expansion of the gyrokinetic equation in the perpendicular scale length of the turbulence, originally derived by Parra and Barnes [Plasma Phys. Controlled Fusion, $\textbf{57}$ 054003, 2015], which allows the use of Fourier basis functions while enabling the effect of radial profile variation to be included in a perturbative way. Radial variation of the magnetic geometry is included by utilizing a global extension of the Grad-Shafranov equation and the Miller equilibrium equations which is obtained through Taylor expansion. Radial boundary conditions that employ multiple flux-tube simulations are also developed, serving as a more physically motivated replacement to the conventional Dirichlet radial boundary conditions that are used in global simulation. It is shown that these new boundary conditions eliminate much of the numerical artefacts generated near the radial boundary when expressing a non-periodic function using a spectral basis. We then benchmark the new approach both linearly and nonlinearly using a number of standard test cases.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Extended electron tails in electrostatic microinstabilities and the nonadiabatic response of passing electrons

Ion-gyroradius-scale microinstabilities typically have a frequency comparable to the ion transit frequency. Due to the small electron-to-ion mass ratio and the large electron transit frequency, it is conventionally assumed that passing electrons respond adiabatically in ion-gyroradius-scale modes. However, in gyrokinetic simulations of ion-gyroradius-scale modes in axisymmetric toroidal magnetic fields, the nonadiabatic response of passing electrons can drive the mode, and generate fluctuations with narrow radial layers, which may have consequences for turbulent transport in a variety of circumstances. In flux tube simulations, in the ballooning representation, these instabilities reveal themselves as modes with extended tails. The small electron-to-ion mass ratio limit of linear gyrokinetics for electrostatic instabilities is presented, in axisymmetric toroidal magnetic geometry, including the nonadiabatic response of passing electrons and associated narrow radial layers. This theory reveals the existence of ion-gyroradius-scale modes driven solely by the nonadiabatic passing electron response, and recovers the usual ion-gyroradius-scale modes driven by the response of ions and trapped electrons, where the nonadiabatic response of passing electrons is small. The collisionless and collisional limits of the theory are considered, demonstrating parallels in structure and physical processes to neoclassical transport theory. By examining initial-value simulations of fastest-growing eigenmodes, the predictions for mass-ratio scaling are tested and verified numerically for a range of collision frequencies. Insights from the small electron-to-ion mass ratio theory may lead to a computationally efficient treatment of extended modes.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Interpreting Radial Correlation Doppler Reflectometry using Gyrokinetic Simulations

A linear response, local model for the DBS amplitude applied to gyrokinetic simulations shows that radial correlation Doppler reflectometry measurements (RCDR, Schirmer et al., Plasma Phys. Control. Fusion 49 1019 (2007)) are not sensitive to the average turbulence radial correlation length, but to a correlation length that depends on the binormal wavenumber $k_\perp$ selected by the Doppler backscattering (DBS) signal. Nonlinear gyrokinetic simulations show that the turbulence naturally exhibits a non-separable power law spectrum in wavenumber space, leading to a power law dependence of the radial correlation length with binormal wavenumber $l_r \sim C k_\perp^{-α} (α\approx 1)$ which agrees with the inverse proportionality relationship between the measured $l_r$ and $k_\perp $ in experiments (Fernandez-Marina et al., Nucl. Fusion 54 072001 (2014)). This offers the possibility of characterizing the eddy aspect ratio in the perpendicular plane to the magnetic field and motivates future use of a non-separable turbulent spectrum to quantitatively interpret RCDR and potentially other turbulence diagnostics. The radial correlation length is only measurable when the radial resolution at the cutoff location $W_n$ satisfies $W_n \ll l_r$, while the measurement becomes dominated by $W_n$ for $W_n \gg l_r$. This suggests that $l_r$ is likely inaccessible for electron-scale DBS measurements ($k_\perpρ_s > 1$). The effect of $W_n$ on ion-scale radial correlation lengths could be non-negligible.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Electrostatic gyrokinetic simulations in Wendelstein 7-X geometry: benchmark between the codes stella and GENE

The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and for the preparation of the next campaigns. At the same time, new stellarator gyrokinetic codes have just been or are being developed. It is therefore desirable to have a sufficiently complete, documented and verified set of gyrokinetic simulations in W7-X geometry against which new codes or upgrades of existing codes can be tested and benchmarked. This paper attemps to provide such a set of simulations in the form of a comprehensive benchmark between the recently developed code stella and the well-established code GENE. The benchmark consists of electrostatic gyrokinetic simulations in W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode (TEM)} stability analyses, computation of linear zonal flow responses and calculation of ITG-driven heat fluxes.

0 citations0 reviewsNo code linkedOpen access
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*.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ITG, TEM and ETG modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. It is also confirmed that thermo-diffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

0 citations0 reviewsNo code linkedOpen access