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

G. W. Hammett

G. 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-phphysics.comp-phphysics.space-phastro-phastro-ph.HEastro-ph.SRmath.NAnlin.CD

Catalog footprint

What is connected

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

17 published item(s)

preprint2022arXiv

Kinetic modeling of neutral transport for a continuum gyrokinetic code

We present the first-of-its-kind coupling of a continuum full-f gyrokinetic turbulence model with a 6D continuum model for kinetic neutrals, carried out using the Gkeyll code. Our objective is to improve the first-principles understanding of the role of neutrals in plasma fueling, detachment, and their interaction with edge plasma profiles and turbulence statistics. Our model includes only atomic hydrogen and incorporates electron-impact ionization, charge exchange, and wall recycling. These features have been successfully verified with analytical predictions and benchmarked with the DEGAS2 Monte Carlo neutral code. We carry out simulations for a scrape-off layer (SOL) with simplified geometry and NSTX parameters. We compare these results to a baseline simulation without neutrals and find that neutral interactions reduce the normalized density fluctuation levels and associated skewness and kurtosis, while increasing auto-correlation times. A flatter density profile is also observed, similar to the SOL density shoulder formation in experimental scenarios with high fueling.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent broadening of electron heat-flux width in electromagnetic gyrokinetic simulations of a helical scrape-off layer model

We demonstrate that cross-field transport in the scrape-off layer (SOL) can be moderately increased by electromagnetic effects in high-beta regimes, resulting in a broader electron heat-flux width on the endplates. This conclusion is taken from full-$f$ electromagnetic gyrokinetic simulations of a helical SOL model that roughly approximates the SOL of the National Spherical Torus Experiment (NSTX). The simulations have been performed with the Gkeyll code, which recently became the first code to demonstrate the capability to simulate electromagnetic gyrokinetic turbulence on open magnetic field lines with sheath boundary conditions. We scan the source rate and thus $β$ so that the normalized pressure gradient (the MHD ballooning parameter $α\propto \partial β/ \partial r \propto β/ L_p$) is scanned over an experimentally-relevant range, $α= 0.3-1.5$. While there is little change in the pressure gradient scale length $L_p$ near the midplane as beta is increased, a 10% increase in cross-field transport near the midplane results in an increase in the electron heat-flux width $λ_q$ and a 25% reduction of the peak electron heat flux to the endplates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Reduction of transport due to magnetic shear in gyrokinetic simulations of the scrape-off layer

The effect of varying magnetic shear on scrape-off layer turbulence and profiles is studied via electromagnetic gyrokinetic simulations of a helical scrape-off layer model. We develop a model helical geometry with magnetic shear and a corresponding field-aligned coordinate system, which is used for simulations with the Gkeyll code. We find that perpendicular transport is reduced in cases with stronger shear, resulting in higher peak particle and heat fluxes to the endplates. Electromagnetic effects slightly increase transport in strong shear cases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Electromagnetic full-$f$ gyrokinetics in the tokamak edge with discontinuous Galerkin methods

We present an energy-conserving discontinuous Galerkin scheme for the full-$f$ electromagnetic gyrokinetic system in the long-wavelength limit. We use the symplectic formulation and solve directly for $\partial A_\parallel/\partial t$, the inductive component of the parallel electric field, using a generalized Ohm's law derived directly from the gyrokinetic equation. Linear benchmarks are performed to verify the implementation and show that the scheme avoids the Ampère cancellation problem. We perform a nonlinear electromagnetic simulation in a helical open-field-line system as a rough model of the tokamak scrape-off layer using parameters from the National Spherical Torus Experiment (NSTX). This is the first published nonlinear electromagnetic gyrokinetic simulation on open field lines. Comparisons are made to a corresponding electrostatic simulation.

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scrape-Off Layer Turbulence in Tokamaks Simulated with a Continuum Gyrokinetic Code

We are developing a new continuum gyrokinetic code, Gkeyll, for use in edge plasma simulations, and here present initial simulations of turbulence on open field lines with model sheath boundary conditions. The code implements an energy conserving discontinuous Galerkin scheme, applicable to a general class of Hamiltonian equations. Several applications to test problems have been done, including a calculation of the parallel heat-flux on divertor plates resulting from an ELM crash in JET, for a 1x/1v SOL scenario explored previously, where the ELM is modeled as a time-dependent intense upstream source. Here we present initial simulations of turbulence on open field lines in the LAPD linear plasma device. We have also done simulations in a helical open-field-line geometry. While various simplifications have been made at present, this still includes some of the key physics of SOL turbulence, such as bad-curvature drive for instabilities and rapid parallel losses with sheath boundary conditions. This is useful for demonstrating the overall feasibility of this approach and for initial physics studies of SOL turbulence. We developed a novel version of DG that uses Maxwellian-weighted basis functions while still preserving exact particle and energy conservation. The Maxwellian-weighted DG method achieves the same error with 4 times less computational cost in 1v, or 16 times lower cost in the 2 velocity dimensions of gyrokinetics (assuming memory bandwidth is the limiting factor).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Gyrokinetic 1D Scrape-Off Layer Model of an ELM Heat Pulse

An electrostatic gyrokinetic-based model is applied to simulate parallel plasma transport in the scrape-off layer to a divertor plate. The authors focus on a test problem that has been studied previously, using parameters chosen to model a heat pulse driven by an edge-localized mode (ELM) in JET. Previous work has used direct particle-in-cell equations with full dynamics, or Vlasov or fluid equations with only parallel dynamics. With the use of the gyrokinetic quasineutrality equation and logical sheath boundary conditions, spatial and temporal resolution requirements are no longer set by the electron Debye length and plasma frequency, respectively. This test problem also helps illustrate some of the physics contained in the Hamiltonian form of the gyrokinetic equations and some of the numerical challenges in developing an edge gyrokinetic code.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale nature of the dissipation range in gyrokinetic simulations of Alfvénic turbulence

Nonlinear energy transfer and dissipation in Alfvén wave turbulence are analyzed in the first gyrokinetic simulation spanning all scales from the tail of the MHD range to the electron gyroradius scale. For typical solar wind parameters at 1 AU, about 30% of the nonlinear energy transfer close to the electron gyroradius scale is mediated by modes in the tail of the MHD cascade. Collisional dissipation occurs across the entire kinetic range $k_\perpρ_i\gtrsim 1$. Both mechanisms thus act on multiple coupled scales, which have to be retained for a comprehensive picture of the dissipation range in Alfvénic turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On discontinuous Galerkin discretizations of second-order derivatives

Some properties of a Local discontinuous Galerkin (LDG) algorithm are demonstrated for the problem of evaluting a second derivative $g = f_{xx}$ for a given $f$. (This is a somewhat unusual problem, but it is useful for understanding the initial transient response of an algorithm for diffusion equations.) LDG uses an auxiliary variable to break this up into two first order equations and then applies techniques by analogy to DG algorithms for advection algorithms. This introduces an asymmetry into the solution that depends on the choice of upwind directions for these two first order equations. When using piecewise linear basis functions, this LDG solution $g_h$ is shown not to converge in an $L_2$ norm because the slopes in each cell diverge. However, when LDG is used in a time-dependent diffusion problem, this error in the second derivative term is transient and rapidly decays away, so that the overall error is bounded. I.e., the LDG approximation $f_h(x,t)$ for a diffusion equation $\partial f / \partial t = f_{xx}$ converges to the proper solution (as has been shown before), even though the initial rate of change $\partial f_h / \partial t$ does not converge. We also show results from the Recovery discontinuous Galerkin (RDG) approach, which gives symmetric solutions that can have higher rates of convergence for a stencil that couples the same number of cells.

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Applications of Large Eddy Simulation methods to gyrokinetic turbulence

The Large Eddy Simulation (LES) approach - solving numerically the large scales of a turbulent system and accounting for the small-scale influence through a model - is applied to nonlinear gyrokinetic systems that are driven by a number of different microinstabilities. Comparisons between modeled, lower resolution, and higher resolution simulations are performed for an experimental measurable quantity, the electron density fluctuation spectrum. Moreover, the validation and applicability of LES is demonstrated through a series of diagnostics based on the free energetics of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparing linear ion-temperature-gradient-driven mode stability of the National Compact Stellarator Experiment and a shaped tokamak

One metric for comparing confinement properties of different magnetic fusion energy configurations is the linear critical gradient of drift wave modes. The critical gradient scale length determines the ratio of the core to pedestal temperature when a plasma is limited to marginal stability in the plasma core. The gyrokinetic turbulence code GS2 was used to calculate critical temperature gradients for the linear, collisionless ion temperature gradient (ITG) mode in the National Compact Stellarator Experiment (NCSX) and a prototypical shaped tokamak, based on the profiles of a JET H-mode shot and the stronger shaping of ARIES-AT. While a concern was that the narrow cross section of NCSX at some toroidal locations would result in steep gradients that drive instabilities more easily, it is found that other stabilizing effects of the stellarator configuration offset this so that the normalized critical gradients for NCSX are competitive with or even better than for the tokamak. For the adiabatic ITG mode, NCSX and the tokamak had similar critical gradients, though beyond marginal stability, NCSX had larger growth rates. However, for the kinetic ITG mode, NCSX had a higher critical gradient and lower growth rates until a/L_T is approximately 1.5 times a/L_{T,crit}, when it surpassed the tokamak's. A discussion of the results presented with respect to a/L_T vs. R/L_T is included.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gyrokinetic studies of the effect of beta on drift-wave stability in NCSX

The gyrokinetic turbulence code GS2 was used to investigate the effects of plasma beta on linear, collisionless ion temperature gradient (ITG) modes and trapped electron modes (TEM) in National Compact Stellarator Experiment (NCSX) geometry. Plasma beta affects stability in two ways: through the equilibrium and through magnetic fluctuations. The first was studied here by comparing ITG and TEM stability in two NCSX equilibria of differing beta values, revealing that the high beta equilibrium was marginally more stable than the low beta equilibrium in the adiabatic-electron ITG mode case. However, the high beta case had a lower kinetic-electron ITG mode critical gradient. Electrostatic and electromagnetic ITG and TEM mode growth rate dependencies on temperature gradient and density gradient were qualitatively similar. The second beta effect is demonstrated via electromagnetic ITG growth rates' dependency on GS2's beta input parameter. A linear benchmark with gyrokinetic codes GENE and GKV-X is also presented.

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

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

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