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

J. W. Burby

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

ResearcherUnclaimed source record
physics.plasm-phmath-phmath.MPphysics.comp-phphysics.flu-dynmath.NAastro-ph.SRcond-mat.mes-hallnlin.CDphysics.acc-ph

Catalog footprint

What is connected

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

20 published item(s)

preprint2020arXiv

General formulas for adiabatic invariants in nearly-periodic Hamiltonian systems

While it is well-known that every nearly-periodic Hamiltonian system possesses an adiabatic invariant, extant methods for computing terms in the adiabatic invariant series are inefficient. The most popular method involves the heavy intermediate calculation of a non-unique near-identity coordinate transformation, even though the adiabatic invariant itself is a uniquely-defined scalar. A less well-known method, developed by S. Omohundro, avoids calculating intermediate sequences of coordinate transformations but is also inefficient as it involves its own sequence of complex intermediate calculations. In order to improve the efficiency of future calculations of adiabatic invariants, we derive generally-applicable, readily computable formulas for the first several terms in the adiabatic invariant series. To demonstrate the utility of these formulas, we apply them to charged particle dynamics in a strong magnetic field and magnetic field-line dynamics when the field lines are nearly closed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Nontrivial topology in the continuous spectrum of a magnetized plasma

Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can coexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the Alfvén continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear Alfvén eigenmode (RSAE) is well known within magnetic fusion as it has been observed in several tokamaks. In examining how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology. More broadly, this work raises the possibility that fusion devices could be carefully tailored to produce other topological states with potentially useful behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Slow manifold reduction for plasma science

The classical Chapman-Enskog procedure admits a substantial geometrical generalization known as slow manifold reduction. This generalization provides a paradigm for deriving and understanding most reduced models in plasma physics that are based on controlled approximations applied to problems with multiple timescales. In this Review we develop the theory of slow manifold reduction with a plasma physics audience in mind. In particular we illustrate (a) how the slow manifold concept may be used to understand \emph{breakdown} of a reduced model over sufficiently-long time intervals, and (b) how a discrete-time analogue of slow manifold theory provides a useful framework for developing implicit integrators for temporally-stiff plasma models. For readers with more advanced mathematical training we also use slow manifold reduction to explain the phenomenon of inheritance of Hamiltonian structure in dissipation-free reduced plasma models. Various facets of the theory are illustrated in the context of the Abraham-Lorentz model of a single charged particle experiencing its own radiation drag. As a culminating example we derive the slow manifold underlying kinetic quasineutral plasma dynamics up to first-order in perturbation theory. This first-order result incorporates several physical effects associated with small deviations from exact charge neutrality that lead to slow drift away from predictions based on the leading-order approximation $n_e = Z_i \,n_i$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-dependent relaxed magnetohydrodynamics -- inclusion of cross helicity constraint using phase-space action

A new formulation of time-dependent Relaxed Magnetohydrodynamics (RxMHD) is derived variationally from Hamilton's Action Principle using microscopic conservation of mass, and macroscopic conservation of total magnetic helicity, cross helicity and entropy, as the only constraints on variations of density, pressure, fluid velocity, and magnetic vector potential over a relaxation domain. A novel phase-space version of the MHD Lagrangian is derived, which gives Euler--Lagrange equations consistent with previous work on exact ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether $\vec{E} + \vec{u}\times\vec{B} = 0$ should be a constraint is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Guiding center dynamics as motion on a formal slow manifold in loop space

Since the late 1950's, the dynamics of a charged particle's ``guiding center" in a strong, inhomogeneous magnetic field have been understood in terms of near-identity coordinate transformations. The basic idea has been to approximately transform away the coupling between the fast gyration around magnetic fields lines and the remaining slow dynamics. This basic understanding now serves as a foundation for describing the kinetic theory of strongly magnetized plasmas. I present a new way to understand guiding center dynamics that does not involve complicated coordinate transformations. Starting from a dynamical systems formulation of the motion of parameterized loops in a charged particle's phase space, I identify a formal slow manifold in loop space. Dynamics on this formal slow manifold are equivalent to guiding center dynamics to all orders in perturbation theory. After demonstrating that loop space dynamics comprises an infinite-dimensional noncanonical Hamiltonian system, I recover the well-known Hamiltonian formulation of guiding center motion by restricting the (pre-) symplectic structure on loop space to the finite-dimensional guiding center formal slow manifold.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally-plane-wave solutions of nonlinear partial differential equations (PDE). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small, and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether's theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by $S^1$. As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Lifting of the Vlasov-Maxwell Bracket by Lie-transform Method

The Vlasov-Maxwell equations possess a Hamiltonian structure expressed in terms of a Hamiltonian functional and a functional bracket. In the present paper, the transformation ("lift") of the Vlasov-Maxwell bracket induced by the dynamical reduction of single-particle dynamics is investigated when the reduction is carried out by Lie-transform perturbation methods. The ultimate goal of this work is to derive explicit Hamiltonian formulations for the guiding-center and gyrokinetic Vlasov-Maxwell equations that have important applications in our understanding of turbulent magnetized plasmas. Here, it is shown that the general form of the reduced Vlasov-Maxwell equations possesses a Hamiltonian structure defined in terms of a reduced Hamiltonian functional and a reduced bracket that automatically satisfies the standard bracket properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Chasing Hamiltonian structure in gyrokinetic theory

Hamiltonian structure is pursued and uncovered in collisional and collisionless gyrokinetic theory. A new Hamiltonian formulation of collisionless electromagnetic theory is presented that is ideally suited to implementation on modern supercomputers. The method used to uncover this structure is described in detail and applied to a number of examples, where several well-known plasma models are endowed with a Hamiltonian structure for the first time. The first energy- and momentum-conserving formulation of full-F collisional gyrokinetics is presented. In an effort to understand the theoretical underpinnings of this result at a deeper level, a \emph{stochastic} Hamiltonian modeling approach is presented and applied to pitch angle scattering. Interestingly, the collision operator produced by the Hamiltonian approach is equal to the Lorentz operator plus higher-order terms, but does not exactly conserve energy. Conversely, the classical Lorentz collision operator is provably not Hamiltonian in the stochastic sense.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Comment on "Symplectic integration of magnetic systems": a proof that the Boris algorithm is not variational

The Boris algorithm for integrating charged particle trajectories in electric and magnetic fields is popular due to its simple implementation, rapid iteration, and observed long-term numerical fidelity. The underlying cause of this long-term fidelity has become a matter of controversy, with one article claiming the method to be symplectic [S. D. Webb, J. Comput. Phys. 270 (2014) 570], and others claiming the method to be volume preserving but not symplectic [e.g. H. Qin et al., Phys. Plasmas 20 (2013) 084503]. To resolve the discrepancy, this letter leverages a discrete Helmholtz condition to demonstrate that no variational formulation of the Boris algorithm exists, indicating that the long-term fidelity should be attributed to the volume-preserving properties of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Equivalence of two independent calculations of the higher order guiding center Lagrangian

The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {\bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {\bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated explicitly.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Field theory and weak Euler-Lagrange equation for classical particle-field systems

It is commonly believed as a fundamental principle that energy-momentum conservation of a physical system is the result of space-time symmetry. However, for classical particle-field systems, e.g., charged particles interacting through self-consistent electromagnetic or electrostatic fields, such a connection has only been cautiously suggested. It has not been formally established. The difficulty is due to the fact that the dynamics of particles and the electromagnetic fields reside on different manifolds. We show how to overcome this difficulty and establish the connection by generalizing the Euler-Lagrange equation, the central component of a field theory, to a so-called weak form. The weak Euler-Lagrange equation induces a new type of flux, called the weak Euler-Lagrange current, which enters conservation laws. Using field theory together with the weak Euler-Lagrange equation developed here, energy-momentum conservation laws that are difficult to find otherwise can be systematically derived from the underlying space-time symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The initial value problem in Lagrangian drift kinetic theory

Existing high-order variational drift kinetic theories contain unphysical rapidly varying modes that are not seen at low-orders. These unphysical modes, which may be rapidly oscillating, damped, or growing, are ushered in by a failure of conventional high-order drift kinetic theory to preserve the structure of its parent model's initial value problem (Vlasov-Poisson for electrostatics, Vlasov-Darwin or Vlasov-Maxwell for electromagnetics.) In short, the system phase space is unphysically enlarged in conventional high-order variational drift kinetic theory. I present an alternative, "renormalized" variational approach to drift kinetic theory that manifestly respects the parent model's initial value problem. The basic philosophy underlying this alternate approach is that high-order drift kinetic theory ought to be derived by truncating the all-orders system phase space Lagrangian instead of the usual "field+particle" Lagrangian. For the sake of clarity, this story is told first through the lens of a finite-dimensional toy model of high-order variational drift kinetics, the analogous full-on drift kinetic results are discussed subsequently. The renormalized drift kinetic system, while just as accurate as conventional formulations, does not support the troublesome rapidly varying mode.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hamiltonian mechanics of generalized eikonal waves

In accordance with the Keller-Maslov global WKB theory, a semiclassical scalar wave field is best encoded as a triple consisting of (i) a Lagrangian submanifold $Λ$ in the ray phase space, (ii) a density $μ$ on $Λ$, and (iii) an overall phase factor $ϕ$. We present the Hamiltonian structure of the Cauchy problem for such a "geometric semiclassical state" in the special case where the wave operator is Hermetian. Variational, symplectic, and Poisson formulations of the time evolution equations for $(Λ,μ,ϕ)$ are identitfied. Because we work in terms of the Keller-Maslov global WKB ansatz, as opposed to the more restrictive $ψ=a \exp(i S/ε)$, all of our results are insensitive to the presence of caustics. In particular, because the variational principle is insensitive to caustics, the latter may be used to construct structure-perserving numerical integrators for scalar wave equations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Initializing and stabilizing variational multistep algorithms for modeling dynamical systems

Backward error initialization and parasitic mode control are well-suited for use in algorithms that arise from a discrete variational principle on phase-space dynamics. Dynamical systems described by degenerate Lagrangians, such as those occurring in phase-space action principles, lead to variational multistep algorithms for the integration of first-order differential equations. As multistep algorithms, an initialization procedure must be chosen and the stability of parasitic modes assessed. The conventional selection of initial conditions using accurate one-step methods does not yield the best numerical performance for smoothness and stability. Instead, backward error initialization identifies a set of initial conditions that minimize the amplitude of undesirable parasitic modes. This issue is especially important in the context of structure-preserving multistep algorithms where numerical damping of the parasitic modes would violate the conservation properties. In the presence of growing parasitic modes, the algorithm may also be periodically re-initialized to prevent the undesired mode from reaching large amplitude. Numerical examples of variational multistep algorithms are presented in which the backward error initialized trajectories outperform those initialized using highly accurate approximations of the true solution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variational integration for ideal magnetohydrodynamics with built-in advection equations

Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variational integrators for perturbed non-canonical Hamiltonian systems

Finite-dimensional non-canonical Hamiltonian systems arise naturally from Hamilton's principle in phase space. We present a method for deriving variational integrators that can be applied to perturbed non-canonical Hamiltonian systems on manifolds based on discretizing this phase-space variational principle. Relative to the perturbation parameter $ε$, this type of integrator can take $O(1)$ time steps with arbitrary accuracy in $ε$ by leveraging the unperturbed dynamics. Moreover, these integrators are coordinate independent in the sense that their time-advance rules transform correctly when passing from one phase space coordinate system to another.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Lagrangian perspective on the stability of ideal MHD equilibria with flow

We take a careful look at two approaches to deriving stability criteria for ideal MHD equilibria. One is based on a tedious analysis of the linearized equations of motion, while the other examines the second variation of the MHD Hamiltonian computed with proper variational constraints. For equilibria without flow, the two approaches are known to be fully consistent. However, for equilibria with flow, the stability criterion obtained from the constrained variation approach was claimed to be stronger than that derived using the linearized equations of motion. We show this claim is incorrect by deriving and comparing both criteria within the same framework. It turns out that the criterion obtained from the constrained variation approach has stricter requirements on the initial perturbations than the other. Such requirements naturally emerge in our new treatment of the constrained variation approach using the Euler-Poincaré structure of ideal MHD, which is more direct and simple than the previous derivation from the Poisson perspective.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Automation of The Guiding Center Expansion

We report on the use of the recently-developed Mathematica package \emph{VEST} (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the transformation order-by-order. This procedure has several novel features. (1) It is designed to allow the user to easily explore the guiding center transformation's numerous non-unique forms or representations. (2) The procedure proceeds entirely in cartesian position and velocity coordinates, thereby producing manifestly gyrogauge invariant results; the commonly-used perpendicular unit vector fields $e_1,e_2$ are never even introduced. (3) It is easy to apply in the derivation of higher-order contributions to the guiding center transformation without fear of human error. Our code therefore stands as a useful tool for exploring subtle issues related to the physics of toroidal momentum conservation in tokamaks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gyrosymmetry: Global Considerations

In the guiding center theory, smooth unit vectors perpendicular to the magnetic field are required to define the gyrophase. The question of global existence of these vectors is addressed using a general result from the theory of characteristic classes. It is found that there is, in certain cases, an obstruction to global existence. In these cases, the gyrophase cannot be defined globally. The implications of this fact on the basic structure of the guiding center theory are discussed. In particular it is demonstrated that the guiding center asymptotic expansion of the equations of motion can still be performed in a globally consistent manner when a single global convention for measuring gyrophase is unavailable. The latter fact is demonstrated directly by deriving a new expression for the guiding-center Poincaré-Cartan form exhibiting no dependence on the choice of perpendicular unit vectors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Hamiltonian Mechanics of Stochastic Acceleration

We show how to find the physical Langevin equation describing the trajectories of particles undergoing collisionless stochastic acceleration. These stochastic differential equations retain not only one-, but two-particle statistics, and inherit the Hamiltonian nature of the underlying microscopic equations. This opens the door to using stochastic variational integrators to perform simulations of stochastic interactions such as Fermi acceleration. We illustrate the theory by applying it to two example problems.

0 citations0 reviewsNo code linkedOpen access