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

Gregory W. Hammett

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

ResearcherUnclaimed source record
physics.plasm-phastro-phphysics.comp-phastro-ph.IMastro-ph.SRnlin.CDphysics.space-ph

Catalog footprint

What is connected

10works
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

10 published item(s)

preprint2026arXiv

Conservative velocity mappings for discontinuous Galerkin kinetics

Continuum computational kinetic plasma models evolve the distribution function of a plasma species $f_s$ on a phase-space grid over time. In many problems of interest the distribution function has limited extent in velocity space; hence, using a uniform, highly refined mesh would be costly and slow. Nonuniform velocity grids can reduce the computational cost by placing more degrees of freedom where $f_s$ is appreciable and fewer where it is not. In this work we introduce a first-of-its kind discontinuous Galerkin approach to nonuniform velocity-space discretization using mapped velocity coordinates. This new method is presented in the context of a gyrokinetic model used to study magnetized plasmas. We create discretizations of collisionless and collisional terms using mappings in a way that exactly conserves particles and energy. Numerical tests of such properties are presented, and we show that this new discretization can reproduce earlier gyrokinetic simulations using grids with up to 6-60 times fewer cells and 22X-60X speed-ups depending on dimensionality, geometry and plasma parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species interactions are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, satisfies the H-Theorem unlike analogous Bhatnagar-Gross-Krook operators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator

A conservative discontinuous Galerkin scheme for a nonlinear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it has a non-decreasing entropy functional, and conserves particles, momentum and energy. Discretely these conservative properties are maintained exactly as well, independent of numerical resolution. In this work the phase space discretization is performed using a novel version of the discontinuous Galerkin scheme, carefully constructed using concepts of weak equality and recovery. Discrete time advancement is carried out with an explicit time-stepping algorithm, whose stability limits we explore. The formulation and implementation within the long-wavelength gyrokinetic solver of Gkeyll are validated with relaxation tests, collisional Landau-damping benchmarks and the study of 5D gyrokinetic turbulence on helical, open field lines.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluid & Gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas (T. N. Bernard, et. al., Phys. of Plasmas 26, 042301 (2019)). Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs which can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, as well as improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Investigating shear flow through continuum gyrokinetic simulations of limiter biasing in the Texas Helimak

Previous limiter-biasing experiments on the Texas Helimak, a simple magnetized torus, have been inconclusive on the effect of flow shear on turbulence levels. To investigate this, the first gyrokinetic simulations of limiter biasing in the Helimak using the plasma physics code Gkeyll have been carried out, and results are presented here. For the scenarios considered, turbulence is mostly driven by the interchange instability, which depends on gradients of equilibrium density profiles. An analysis of both experimental and simulation data demonstrates that shear rates are mostly less than than local linear growth rates, and not all requirements for shear stabilization are met. Rather, the mostly vertical shear flow has an important effect on bulk transport and experimental equilibrium density profiles, and changes to the gradients correspond to changes in turbulence levels.

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Fast Semi-implicit Method for Anisotropic Diffusion

Simple finite differencing of the anisotropic diffusion equation, where diffusion is only along a given direction, does not ensure that the numerically calculated heat fluxes are in the correct direction. This can lead to negative temperatures for the anisotropic thermal diffusion equation. In a previous paper we proposed a monotonicity-preserving explicit method which uses limiters (analogous to those used in the solution of hyperbolic equations) to interpolate the temperature gradients at cell faces. However, being explicit, this method was limited by a restrictive Courant-Friedrichs-Lewy (CFL) stability timestep. Here we propose a fast, conservative, directionally-split, semi-implicit method which is second order accurate in space, is stable for large timesteps, and is easy to implement in parallel. Although not strictly monotonicity-preserving, our method gives only small amplitude temperature oscillations at large temperature gradients, and the oscillations are damped in time. With numerical experiments we show that our semi-implicit method can achieve large speed-ups compared to the explicit method, without seriously violating the monotonicity constraint. This method can also be applied to isotropic diffusion, both on regular and distorted meshes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gyrokinetic statistical absolute equilibrium and turbulence

A paradigm based on the absolute equilibrium of Galerkin-truncated inviscid systems to aid in understanding turbulence [T.-D. Lee, "On some statistical properties of hydrodynamical and magnetohydrodynamical fields," Q. Appl. Math. 10, 69 (1952)] is taken to study gyrokinetic plasma turbulence: A finite set of Fourier modes of the collisionless gyrokinetic equations are kept and the statistical equilibria are calculated; possible implications for plasma turbulence in various situations are discussed. For the case of two spatial and one velocity dimension, in the calculation with discretization also of velocity $v$ with $N$ grid points (where $N+1$ quantities are conserved, corresponding to an energy invariant and $N$ entropy-related invariants), the negative temperature states, corresponding to the condensation of the generalized energy into the lowest modes, are found. This indicates a generic feature of inverse energy cascade. Comparisons are made with some classical results, such as those of Charney-Hasegawa-Mima in the cold-ion limit. There is a universal shape for statistical equilibrium of gyrokinetics in three spatial and two velocity dimensions with just one conserved quantity. Possible physical relevance to turbulence, such as ITG zonal flows, and to a critical balance hypothesis are also discussed.

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

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Electron Heating in Hot Accretion Flows

Local (shearing box) simulations of the nonlinear evolution of the magnetorotational instability in a collisionless plasma show that angular momentum transport by pressure anisotropy ($p_\perp \ne p_\parallel$, where the directions are defined with respect to the local magnetic field) is comparable to that due to the Maxwell and Reynolds stresses. Pressure anisotropy, which is effectively a large-scale viscosity, arises because of adiabatic invariants related to $p_\perp$ and $p_\parallel$ in a fluctuating magnetic field. In a collisionless plasma, the magnitude of the pressure anisotropy, and thus the viscosity, is determined by kinetic instabilities at the cyclotron frequency. Our simulations show that $\sim 50$ % of the gravitational potential energy is directly converted into heat at large scales by the viscous stress (the remaining energy is lost to grid-scale numerical dissipation of kinetic and magnetic energy). We show that electrons receive a significant fraction ($\sim [T_e/T_i]^{1/2}$) of this dissipated energy. Employing this heating by an anisotropic viscous stress in one dimensional models of radiatively inefficient accretion flows, we find that the radiative efficiency of the flow is greater than 0.5% for $\dot{M} \gtrsim 10^{-4} \dot{M}_{Edd}$. Thus a low accretion rate, rather than just a low radiative efficiency, is necessary to explain the low luminosity of many accreting black holes. For Sgr A* in the Galactic Center, our predicted radiative efficiencies imply an accretion rate of $\approx 3 \times 10^{-8} M_\odot {\rm yr^{-1}}$ and an electron temperature of $\approx 3 \times 10^{10}$ K at $\approx 10$ Schwarzschild radii; the latter is consistent with the brightness temperature inferred from VLBI observations.

0 citations0 reviewsNo code linkedOpen access