BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. Barnes

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

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-phastro-ph.HEphysics.flu-dynphysics.space-phastro-ph.SRphysics.acc-ph

Catalog footprint

What is connected

25works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

25 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
preprint2015arXiv

Dependence of intrinsic rotation reversals on collisionality in MAST

Tokamak plasmas rotate even without external injection of momentum. A Doppler backscattering system installed at MAST has allowed this intrinsic rotation to be studied in Ohmic L-mode and H-mode plasmas, including the first observation of intrinsic rotation reversals in a spherical tokamak. Experimental results are compared to a novel 1D model, which captures the collisionality dependence of the radial transport of toroidal angular momentum due to the effect of neoclassical flows on turbulent fluctuations. The model is able to accurately reproduce the change in sign of core toroidal rotation, using experimental density and temperature profiles from shots with rotation reversals as inputs and no free parameters fit to experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-m Kink/Tearing Modes in Cylindrical Geometry

The global ideal kink equation, for cylindrical geometry and zero beta, is simplified in the high poloidal mode number limit and used to determine the tearing stability parameter, $Δ^\prime$. In the presence of a steep monotonic current gradient, $Δ^\prime$ becomes a function of a parameter, $σ_0$, characterising the ratio of the maximum current gradient to magnetic shear, and $x_s$, characterising the separation of the resonant surface from the maximum of the current gradient. In equilibria containing a current "spike", so that there is a non-monotonic current profile, $Δ^\prime$ also depends on two parameters: $κ$, related to the ratio of the curvature of the current density at its maximum to the magnetic shear, and $x_s$, which now represents the separation of the resonance from the point of maximum current density. The relation of our results to earlier studies of tearing modes and to recent gyro-kinetic calculations of current driven instabilities, is discussed, together with potential implications for the stability of the tokamak pedestal.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Intrinsic rotation driven by non-Maxwellian equilibria in tokamak plasmas

The effect of small deviations from a Maxwellian equilibrium on turbulent momentum transport in tokamak plasmas is considered. These non-Maxwellian features, arising from diamagnetic effects, introduce a strong dependence of the radial flux of co-current toroidal angular momentum on collisionality: As the plasma goes from nearly collisionless to weakly collisional, the flux reverses direction from radially inward to outward. This indicates a collisionality-dependent transition from peaked to hollow rotation profiles, consistent with experimental observations of intrinsic rotation.

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Experimental Signatures of Critically Balanced Turbulence in MAST

Beam Emission Spectroscopy (BES) measurements of ion-scale density fluctuations in the MAST tokamak are used to show that the turbulence correlation time, the drift time associated with ion temperature or density gradients, the particle (ion) streaming time along the magnetic field and the magnetic drift time are consistently comparable, suggesting a "critically balanced" turbulence determined by the local equilibrium. The resulting scalings of the poloidal and radial correlation lengths are derived and tested. The nonlinear time inferred from the density fluctuations is longer than the other times; its ratio to the correlation time scales as $ν_{*i}^{-0.8\pm0.1}$, where $ν_{*i}=$ ion collision rate/streaming rate. This is consistent with turbulent decorrelation being controlled by a zonal component, invisible to the BES, with an amplitude exceeding the drift waves' by $\sim ν_{*i}^{-0.8}$.

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Critically balanced ion temperature gradient turbulence in fusion plasmas

Scaling laws for ion temperature gradient driven turbulence in magnetized toroidal plasmas are derived and compared with direct numerical simulations. Predicted dependences of turbulence fluctuation amplitudes, spatial scales, and resulting heat fluxes on temperature gradient and magnetic field line pitch are found to agree with numerical results in both the driving and inertial ranges. Evidence is provided to support the critical balance conjecture that parallel streaming and nonlinear perpendicular decorrelation times are comparable at all spatial scales, leading to a scaling relationship between parallel and perpendicular spatial scales. This indicates that even strongly magnetized plasma turbulence is intrinsically three-dimensional.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Injection and dump considerations for a 16.5 TeV HE-LHC

Injection and beam dumping is considered for a 16.5 TeV hadron accelerator in the current LHC tunnel, with an injection energy in the range 1 - 1.3 TeV. The present systems are described and the possible upgrade scenarios investigated for higher beam rigidity. In addition to the required equipment performance, the machine protection related aspects are explored. The expected constraints on the machine layout are also given. The technological challenges for the different equipment subsystems are detailed, and areas where R&D is necessary are highlighted.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Momentum injection in tokamak plasmas and transitions to reduced transport

The effect of momentum injection on the temperature gradient in tokamak plasmas is studied. A plausible scenario for transitions to reduced transport regimes is proposed. The transition happens when there is sufficient momentum input so that the velocity shear can suppress or reduce the turbulence. However, it is possible to drive too much velocity shear and rekindle the turbulent transport. The optimal level of momentum injection is determined. The reduction in transport is maximized in the regions of low or zero magnetic shear.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transport Bifurcation Induced by Sheared Toroidal Flow in Tokamak Plasmas

First-principles numerical simulations are used to describe a transport bifurcation in a differentially rotating tokamak plasma. Such a bifurcation is more probable in a region of zero magnetic shear than one of finite magnetic shear because in the former case the component of the sheared toroidal flow that is perpendicular to the magnetic field has the strongest suppressing effect on the turbulence. In the zero-magnetic-shear regime, there are no growing linear eigenmodes at any finite value of flow shear. However, subcritical turbulence can be sustained, owing to the transient growth of modes driven by the ion temperature gradient (ITG) and the parallel velocity gradient (PVG). Nonetheless, in a parameter space containing a wide range of temperature gradients and velocity shears, there is a sizeable window where all turbulence is suppressed. Combined with the relatively low transport of momentum by collisional (neoclassical) mechanisms, this produces the conditions for a bifurcation from low to high temperature and velocity gradients. The path of this bifurcation is mapped out using interpolation from a large number of simulations. Numerical simulations are also used to construct a parametric model which accurately describes the combined effect of the temperature gradient and the flow gradient over a wide range of their values. Using this parametric model, it is shown that in this reduced-transport state, heat is transported almost neoclassically, while momentum transport is dominated by subcritical PVG turbulence. It is further shown that for any given input of torque, there is an optimum input of heat which maximises the temperature gradient. The parametric model describes both the behaviour of the subcritical turbulence and the complicated effect of the flow shear on the transport stiffness. It may prove useful for transport modelling of tokamaks with sheared flows.

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Transport Bifurcation in a Rotating Tokamak Plasma

The effect of flow shear on turbulent transport in tokamaks is studied numerically in the experimentally relevant limit of zero magnetic shear. It is found that the plasma is linearly stable for all non-zero flow shear values, but that subcritical turbulence can be sustained nonlinearly at a wide range of temperature gradients. Flow shear increases the nonlinear temperature gradient threshold for turbulence but also increases the sensitivity of the heat flux to changes in the temperature gradient, except over a small range near the threshold where the sensitivity is decreased. A bifurcation in the equilibrium gradients is found: for a given input of heat, it is possible, by varying the applied torque, to trigger a transition to significantly higher temperature and flow gradients.

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

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

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

0 citations0 reviewsNo code linkedOpen access