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

P. J. Morrison

P. J. Morrison 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.flu-dynnlin.CDchao-dyncond-mat.stat-mechmath.DSastro-phastro-ph.IMastro-ph.SRgr-qcmath.DGmath.NAnlin.PSphysics.class-phphysics.opticsphysics.space-phquant-ph

Catalog footprint

What is connected

44works
19topics
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

44 published item(s)

preprint2026arXiv

Thermodynamics of driven systems via the Kuramoto-Sivashinsky equation

We examine the differences between the driven turbulence described by the Kuramoto-Sivashinsky (KS) equation and the second law of thermodynamics. A general velocity and entropy density system is analyzed with the unified thermodynamic algorithm of metriplectic dynamics, and we show that the positive spectra of the KS equation due to an external energy source prevent its metriplectic description. A variant of the KS equation is produced that monotonically generates an entropy, but the only equilibria of this variant system are spatially constant. Numerical experiments are performed comparing the evolution of the KS equation and its thermodynamic variant. The entropy of this thermodynamic system is increased further by the driving effects of the KS equation, reconciling the generation of entropy with the energy source of the KS equation. Further numerical experiments restrict the positive spectra in the KS equation to determine the effect on the system time evolution. While rescaling the growth rates of instabilities reproduces similar behavior on a slower time scale, introduction of individual positive spectra reproduces the formation of equilibria, relative equilibria, and a transition to chaos.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Energy-Casimir, dynamically accessible, and Lagrangian stability of extended magnetohydrodynamic equilibria

The formal stability analysis of Eulerian extended magnetohydrodynamics (XMHD) equilibria is considered within the noncanonical Hamiltonian framework by means of the energy-Casimir variational principle and the dynamically accessible stability method. Specifically, we find explicit sufficient stability conditions for axisymmetric XMHD and Hall MHD (HMHD) equilibria with toroidal flow and for equilibria with arbitrary flow under constrained perturbations. The dynamically accessible, second-order variation of the Hamiltonian, which can potentially provide explicit stability criteria for generic equilibria, is also obtained. Moreover, we examine the Lagrangian stability of the general quasineutral two-fluid model written in terms of MHD-like variables, by finding the action and the Hamiltonian functionals of the linearized dynamics, working within a mixed Lagrangian-Eulerian framework. Upon neglecting electron mass, we derive a HMHD energy principle, and in addition, the perturbed induction equation arises from Hamilton's equations of motion in view of a consistency condition for the relation between the perturbed magnetic potential and the canonical variables.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lagrangian and Dirac constraints for the ideal incompressible fluid and magnetohydrodynamics

The incompressibility constraint for fluid flow was imposed by Lagrange in the so-called Lagrangian variable description using his method of multipliers in the Lagrangian (variational) formulation. An alternative is the imposition of incompressibility in the Eulerian variable description by a generalization of Dirac's constraint method using noncanonical Poisson brackets. Here it is shown how to impose the incompressibility constraint using Dirac's method in terms of both the canonical Poisson brackets in the Lagrangian variable description and the noncanonical Poisson brackets in the Eulerian description, allowing for the advection of density. Both cases give dynamics of infinite-dimensional geodesic flow on the group of volume preserving diffeomorphisms and explicit expressions for this dynamics in terms of the constraints and original variables is given. Because Lagrangian and Eulerian conservation laws are not identical, comparison of the various methods is made.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Derivation of the Hall and Extended Magnetohydrodynamics Brackets

There are several plasma models intermediate in complexity between ideal magnetohydrodynamics (MHD) and two-fluid theory, with Hall and Extended MHD being two important examples. In this paper we investigate several aspects of these theories, with the ultimate goal of deriving the noncanonical Poisson brackets used in their Hamiltonian formulations. We present fully Lagrangian actions for each, as opposed to the fully Eulerian, or mixed Eulerian-Lagrangian, actions that have appeared previously. As an important step in this process we exhibit each theory's two advected fluxes (in analogy to ideal MHD's advected magnetic flux), discovering also that with the correct choice of gauge they have corresponding Lie-dragged potentials resembling the electromagnetic vector potential, and associated conserved helicities. Finally, using the Euler-Lagrange map, we show how to derive the noncanonical Eulerian brackets from canonical Lagrangian ones.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hamiltonian Magnetohydrodynamics: Lagrangian, Eulerian, and Dynamically Accessible Stability -- Examples with Translation Symmetry

Because different constraints are imposed, stability conditions for dissipationless fluids and magnetofluids may take different forms when derived within the Lagrangian, Eulerian (energy-Casimir), or dynamical accessible frameworks. This is in particular the case when flows are present. These differences are explored explicitly by working out in detail two magnetohydrodynamic examples: convection against gravity in a stratified fluid and translationally invariant perturbations of a rotating magnetized plasma pinch. In this second example we show in explicit form how to perform the time-dependent relabeling introduced in Andreussi {\it et al.}\ [Phys.\ Plasmas {\bf20}, 092104 (2013)] that makes it possible to reformulate Eulerian equilibria with flows as Lagrangian equilibria in the relabeled variables. The procedures detailed in the present article provide a paradigm that can be applied to more general plasma configurations and in addition extended to more general plasma descriptions where dissipation is absent.

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
preprint2016arXiv

Multiscale Equatorial Electrojet Turbulence: Energy Conservation, Coupling, and Cascades in a Baseline 2-D Fluid Model

Progress in understanding the coupling between plasma instabilities in the equatorial electrojet based on a unified fluid model is reported. A deeper understanding of the linear and nonlinear evolution and the coupling of the gradient-drift and Farley-Buneman instabilities is achieved by studying the effect of different combinations of the density-gradient scale-lengths (Ln) and cross-field (E?B) drifts on the plasma turbulence. Mechanisms and channels of energy transfer are illucidated for these multiscale instabilities. Energy for the unified model is examined, including the injected, conservative redistribution (between fields and scales), and ultimate dissipation. Various physical mechanisms involved in the energetics are categorized as sources, sinks, nonlinear transfer, and coupling to show that the system satisfies the fundamental law of energy Oonservation. The physics of the nonlinear transfer terms is studied to identify their roles in producing energy cascades { the transference of energy from the dominant unstable wavenumbers and driving irregularities of different scale-lengths. The theory of two-step energy cascading to generate the 3-meter plasma irregularities in the equatorial electrojet is verified. In addition, the nonlinearity of the system allows the possibility for a reverse energy cascading, responsible for generating large-scale plasma structures at the top of the electrojet.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Hamiltonian Five-Field Gyrofluid Model

A Lie-Poisson bracket is presented for a five-field gyrofluid model, thereby showing the model to be Hamiltonian. The model includes the effects of magnetic field curvature and describes the evolution of the electron and ion gyro-center densities, the parallel component of the ion and electron velocities, and the ion temperature. The quasineutrality property and Ampere's law determine respectively the electrostatic potential and magnetic flux. The Casimir invariants are presented, and shown to be associated to five Lagrangian invariants advected by distinct velocity fields. A linear, local study of the model is conducted both with and without Landau and diamagnetic resonant damping terms. Stability criteria and dispersion relations for the electrostatic and the electromagnetic cases are derived and compared with their analogs for fluid and kinetic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hamiltonian closures for fluid models with four moments by dimensional analysis

Fluid reductions of the Vlasov-Amp{è}re equations that preserve the Hamiltonian structure of the parent kinetic model are investigated. Hamiltonian closures using the first four moments of the Vlasov distribution are obtained, and all closures provided by a dimensional analysis procedure for satisfying the Jacobi identity are identified. Two Hamiltonian models emerge, for which the explicit closures are given, along with their Poisson brackets and Casimir invariants.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hamiltonian fluid closures of the Vlasov-Amp{è}re equations: from water-bags to N moment models

Moment closures of the Vlasov-Amp{è}re system, whereby higher moments are represented as functions of lower moments with the constraint that the resulting fluid system remains Hamiltonian, are investigated by using water-bag theory. The link between the water-bag formalism and fluid models that involve density, fluid velocity, pressure and higher moments is established by introducing suitable thermodynamic variables. The cases of one, two and three water-bags are treated and their Hamiltonian structures are provided. In each case, we give the associated fluid closures and we discuss their Casimir invariants. We show how the method can be extended to an arbitrary number of fields, i.e., an arbitrary number of water-bags and associated moments. The thermodynamic interpretation of the resulting models is discussed. Finally, a general procedure to derive Hamiltonian N-field fluid models is proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Remarkable connections between extended magnetohydrodynamics models

Through the use of suitable variable transformations, the commonality of all extended magnetohydrodynamics (MHD) models is established. Remarkable correspondences between the Poisson brackets of inertialess Hall MHD and inertial MHD (which has electron inertia, but not the Hall drift) and extended MHD (which has both effects), are established. The helicities (two in all) for each of these models are obtained through these correspondences. The commonality of all the extended MHD models is traced to the existence of two Lie-dragged 2-forms, which are closely associated with the canonical momenta of the two underlying species. The Lie-dragging of these 2-forms by suitable velocities also leads to the correct equations of motion. The Hall MHD Poisson bracket is analyzed in detail, and the Jacobi identity is verified through a detailed proof and this proof ensures the Jacobi identity for the Poisson brackets of all the models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theory and applications of the Vlasov equation

Forty articles have been recently published in EPJD as contributions to the topical issue "Theory and applications of the Vlasov equation". The aim of this topical issue was to provide a forum for the presentation of a broad variety of scientific results involving the Vlasov equation. In this editorial, after some introductory notes, a brief account is given of the main points addressed in these papers and of the perspectives they open.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Weakly nonlinear dynamics in noncanonical Hamiltonian systems with applications to fluids and plasmas

A method, called beatification, is presented for rapidly extracting weakly nonlinear Hamiltonian systems that describe the dynamics near equilibria for systems possessing Hamiltonian form in terms of noncanonical Poisson brackets. The procedure applies to systems like fluids and plasmas in terms of Eulerian variables that have such noncanonical Poisson brackets, i.e., brackets with nonstandard and possibly degenerate form. A collection of examples of both finite and infinite dimensions is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A hierarchy of noncanonical Hamiltonian systems: circulation laws in an extended phase space

The dynamics of an ideal fluid or plasma is constrained by topological invariants such as the circulation of (canonical) momentum or, equivalently, the flux of the vorticity or magnetic fields. In the Hamiltonian formalism, topological invariants restrict the orbits to submanifolds of the phase space. While the coadjoint orbits have a natural symplectic structure, the global geometry of the degenerate (constrained) Poisson manifold can be very complex. Some invariants are represented by the center of the Poisson algebra (i.e., the Casimir elements such as the helicities), and then, the global structure of phase space is delineated by Casimir leaves. However, a general constraint is not necessarily integrable, which precludes the existence of an appropriate Casimir element; the circulation is an example of such an invariant. In this work, we formulate a systematic method to embed a Hamiltonian system in an extended phase space; we introduce mock fields and extend the Poisson algebra so that the mock fields are Lie-dragged by the flow vector. A mock field defines a new Casimir element, a cross helicity, which represents topological constraints including the circulation. Unearthing a Casimir element brings about immense advantage in the study of dynamics and equilibria -- the so-called energy-Casimir method becomes ready available. Yet, a mock field does not a priori have a physical meaning. Here we proffer an interpretation of a Casimir element obtained, e.g., by such a construction as an adiabatic invariant associated with a hidden "microscopic" angle variable, and in this way give the mock field a physical meaning. We proceed further and consider a perturbation of the Hamiltonian by a canonical pair, composed of the Casimir element and the angle, that causes the topological constraint to be unfrozen. The theory is applied to the tearing modes of MHD.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Action Principles for Extended MHD Models

The general, non-dissipative, two-fluid model in plasma physics is Hamiltonian, but this property is sometimes lost or obscured in the process of deriving simplified (or reduced) two-fluid or one-fluid models from the two-fluid equations of motion. To ensure that the reduced models are Hamiltonian, we start with the general two-fluid action functional, and make all the approximations, changes of variables, and expansions directly within the action context. The resulting equations are then mapped to the Eulerian fluid variables using a novel nonlocal Lagrange-Euler map. Using this method, we recover Lüst's general two-fluid model, extended MHD, Hall MHD, and electron MHD from a unified framework. The variational formulation allows us to use Noether's theorem to derive conserved quantities for each symmetry of the action.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Experimental determination of radiated internal wave power without pressure field data

We present a method to determine, using only velocity field data, the time-averaged energy flux $\left<\boldsymbol{J}\right>$ and total radiated power $P$ for two-dimensional internal gravity waves. Both $\left<\boldsymbol{J}\right>$ and $P$ are determined from expressions involving only a scalar function, the stream function $ψ$. We test the method using data from a direct numerical simulation for tidal flow of a stratified fluid past a knife edge. The results for the radiated internal wave power given by the stream function method agree to within 0.5% with results obtained using pressure and velocity data from the numerical simulation. The results for the radiated power computed from the stream function agree well with power computed from the velocity and pressure if the starting point for the stream function computation is on a solid boundary, but if a boundary point is not available, care must be taken to choose an appropriate starting point. We also test the stream function method by applying it to laboratory data for tidal flow past a knife edge, and the results are found to agree with the direct numerical simulation. Supplementary Material includes a Matlab code with a graphical user interface (GUI) that can be used to compute the energy flux and power from any two-dimensional velocity field data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hamiltonian and action formalisms for two-dimensional gyroviscous MHD

A general procedure for constructing action principles for continuum models via a generalization of Hamilton's principle of mechanics is described. Through the procedure, an action principle for a gyroviscous magnetohydrodynamics (MHD) model is constructed. The model is shown to agree with a reduced version of Braginskii's fluid equations. The construction reveals the origin of the gyromap, a device used to derive previous gyrofluid models. Also, a systematic reduction procedure is presented for obtaining the Hamiltonian structure in terms of the noncanonical Poisson bracket. The construction procedure yields a class of Casimir invariants, which are then used to variational principles for equilibrium equations with flow and gyroviscosity. The procedure for obtaining reduced fluid models with gyroviscosity is also described.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Inertial magnetohydrodynamics

A version of extended magnetohydrodynamics (MHD) that incorporates electron inertia is obtained by constructing an action principle. Unlike MHD which freezes in magnetic flux, the present theory freezes in an alternative flux related to the electron canonical momentum. The associated Hamiltonian formulation is derived and reduced models that have previously been used to describe collisionless reconnection are obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On energy conservation in extended magnetohydrodynamics

A systematic study of energy conservation for extended magnetohydrodynamic (MHD) models that include Hall terms and electron inertia is performed. It is observed that commonly used models do not conserve energy in the ideal limit, i.e., when viscosity and resistivity are neglected. In particular, a term in the momentum equation that is often neglected is seen to be needed for conservation of energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The action principle for generalized fluid motion including gyroviscosity

A general set of fluid equations that allow for energy-conserving momentum transport by gyroscopic motion of fluid elements is obtained. The equations are produced by a class of action principles that yield a large subset of the known fluid and magnetofluid models, including gyroviscosity. Analysis of the action principle yields broad, model-independent results regarding the conservation laws of energy and linear and angular momenta. The formalism is illustrated by studying fluid models with intrinsic angular momentum that may appear in the contexts of condensed matter, biological, and other areas of physics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Continuum Hamiltonian Hopf Bifurcation II

Building on the development of [MOR13], bifurcation of unstable modes that emerge from continuous spectra in a class of infinite-dimensional noncanonical Hamiltonian systems is investigated. Of main interest is a bifurcation termed the continuum Hamiltonian Hopf (CHH) bifurcation, which is an infinite-dimensional analog of the usual Hamiltonian Hopf (HH) bifurcation. Necessary notions pertaining to spectra, structural stability, signature of the continuous spectra, and normal forms are described. The theory developed is applicable to a wide class of 2+1 noncanonical Hamiltonian matter models, but the specific example of the Vlasov-Poisson system linearized about homogeneous (spatially independent) equilibria is treated in detail. For this example, structural (in)stability is established in an appropriate functional analytic setting, and two kinds of bifurcations are considered, one at infinite and one at finite wavenumber. After defining and describing the notion of dynamical accessibility, Kreĭn-like theorems regarding signature and bifurcation are proven. In addition, a canonical Hamiltonian example, composed of a negative energy oscillator coupled to a heat bath, is treated and our development is compared to pervious work in this context. A careful counting of eigenvalues, with and without symmetry, is undertaken, leading to the definition of a degenerate continuum steady state (CSS) bifurcation. It is described how the CHH and CSS bifurcations can be viewed as linear normal forms associated with the nonlinear single-wave model described in [BAL13], which is a companion to the present work and that of [MOR13].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Deformation of vortex patches by boundaries

The deformation of two-dimensional vortex patches in the vicinity of fluid boundaries is investigated. The presence of a boundary causes an initially circular patch of uniform vorticity to deform. Sufficiently far away from the boundary, the deformed shape is well approximated by an ellipse. This leading order elliptical deformation is investigated via the elliptic moment model of Melander, Zabusky & Styczek [M. V. Melander, N. J. Zabusky & A. S. Styczek, J. Fluid. Mech., 167, 95 (1986)]. When the boundary is straight, the centre of the elliptic patch remains at a constant distance from the boundary, and the motion is integrable. Furthermore, since the straining flow acting on the patch is constant in time, the problem is that of an elliptic vortex patch in constant strain, which was analysed by Kida [S. Kida, J. Phys. Soc. Japan, 50, 3517 (1981)]. For more complicated boundary shapes, such as a square corner, the motion is no longer integrable. Instead, there is an adiabatic invariant for the motion. This adiabatic invariant arises due to the separation in times scales between the relatively rapid time scale associated with the rotation of the patch and the slower time scale associated with the self-advection of the patch along the boundary. The interaction of a vortex patch with a circular island is also considered. Without a background flow, conservation of angular impulse implies that the motion is again integrable. The addition of an irrotational flow past the island can drive the patch towards the boundary, leading to the possibility of large deformations and breakup.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hamiltonian Magnetohydrodynamics: Lagrangian, Eulerian, and Dynamically Accessible Stability - Theory

Stability conditions of magnetized plasma flows are obtained by exploiting the Hamiltonian structure of the magnetohydrodynamics (MHD) equations and, in particular, by using three kinds of energy principles. First, the Lagrangian variable energy principle is described and sufficient stability conditions are presented. Next, plasma flows are described in terms of Eulerian variables and the noncanonical Hamiltonian formulation of MHD is exploited. For symmetric equilibria, the energy-Casimir principle is expanded to second order and sufficient conditions for stability to symmetric perturbation are obtained. Then, dynamically accessible variations, i.e. variations that explicitly preserve invariants of the system, are introduced and the respective energy principle is considered. General criteria for stability are obtained, along with comparisons between the three different approaches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local thermodynamics of a magnetized, anisotropic plasma

An expression for the internal energy of a fluid element in a weakly coupled, magnetized, anisotropic plasma is derived from first principles. The result is a function of entropy, particle density and magnetic field, and as such plays the role of a thermodynamic potential: it determines in principle all thermodynamic properties of the fluid element. In particular it provides equations of state for the magnetized plasma. The derivation uses familiar fluid equations, a few elements of kinetic theory, the MHD version of Faraday's law, and certain familiar stability and regularity conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Nonlinear variational method for predicting fast collisionless magnetic reconnection

A mechanism for fast magnetic reconnection in collisionless plasma is studied for understanding sawtooth collapse in tokamak discharges by using a two-fluid model for cold ions and electrons. Explosive growth of the tearing mode enabled by electron inertia is analytically estimated by using an energy principle with a nonlinear displacement map. Decrease of the potential energy in the nonlinear regime (where the island width exceeds the electron skin depth) is found to be steeper than in the linear regime, resulting in accelerated reconnection. Release of potential energy by such a fluid displacement leads to unsteady and strong convective flow, which is not damped by the small dissipation effects in high-temperature tokamak plasmas. Direct numerical simulation in slab geometry substantiates the theoretical prediction of the nonlinear growth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pattern formation in Hamiltonian systems with continuous spectra; a normal-form single-wave model

Pattern formation in biological, chemical and physical problems has received considerable attention, with much attention paid to dissipative systems. For example, the Ginzburg--Landau equation is a normal form that describes pattern formation due to the appearance of a single mode of instability in a wide variety of dissipative problems. In a similar vein, a certain "single-wave model" arises in many physical contexts that share common pattern forming behavior. These systems have Hamiltonian structure, and the single-wave model is a kind of Hamiltonian mean-field theory describing the patterns that form in phase space. The single-wave model was originally derived in the context of nonlinear plasma theory, where it describes the behavior near threshold and subsequent nonlinear evolution of unstable plasma waves. However, the single-wave model also arises in fluid mechanics, specifically shear-flow and vortex dynamics, galactic dynamics, the XY and Potts models of condensed matter physics, and other Hamiltonian theories characterized by mean field interaction. We demonstrate, by a suitable asymptotic analysis, how the single-wave model emerges from a large class of nonlinear advection-transport theories. An essential ingredient for the reduction is that the Hamiltonian system has a continuous spectrum in the linear stability problem, arising not from an infinite spatial domain but from singular resonances along curves in phase space whereat wavespeeds match material speeds (wave-particle resonances in the plasma problem, or critical levels in fluid problems). The dynamics of the continuous spectrum is manifest as the phenomenon of Landau damping when the system is ... Such dynamical phenomena have been rediscovered in different contexts, which is unsurprising in view of the normal-form character of the single-wave model.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Rank change in Poisson dynamical systems

It is shown in this paper how a connection may be made between the symmetry generators of the Hamiltonian (or potential) invariant under a symmetry group $G$, and the subcasimirs that come about when the rank of the Poisson structure of a dynamical system drops by an even integer. This {\em kinematics-dynamics} connection is made by using the algebraic geometry of the orbit space in the vicinity of rank change, and the extra null eigenvectors of the mass matrix (Hessian with respect to symmetry generators) of the Hamiltonian (or potential). Some physical interpretations of this point of view include a control-theoretic prescription to study stability on various symplectic leaves of the Poisson structure. Methods of Invariant Theory are utilized to provide parametrization for the leaves of a Poisson dynamical system for the case where a compact Lie group acts properly on the phase space, which is assumed to be modeled by Poisson geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Response to Comment on `Undamped electrostatic plasma waves' [Phys. Plasmas 19, 092103 (2012)]

Numerical and experimental evidence is given for the occurrence of the plateau states and concomitant corner modes proposed in \cite{valentini12}. It is argued that these states provide a better description of reality for small amplitude off-dispersion disturbances than the conventional Bernstein-Greene-Kruskal or cnoidal states such as those proposed in \cite{comment}

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Unfreezing Casimir invariants: singular perturbations giving rise to forbidden instabilities

The infinite-dimensional mechanics of fluids and plasmas can be formulated as "noncanonical" Hamiltonian systems on a phase space of Eulerian variables. Singularities of the Poisson bracket operator produce singular Casimir elements that foliate the phase space, imposing topological constraints on the dynamics. Here we proffer a physical interpretation of Casimir elements as \emph{adiabatic invariants} ---upon coarse graining microscopic angle variables, we obtain a macroscopic hierarchy on which the separated action variables become adiabatic invariants. On reflection, a Casimir element may be \emph{unfrozen} by recovering a corresponding angle variable; such an increase in the number of degrees of freedom is, then, formulated as a \emph{singular perturbation}. As an example, we propose a canonization of the resonant-singularity of the Poisson bracket operator of the linearized magnetohydrodynamics equations, by which the ideal obstacle (resonant Casimir element) constraining the dynamics is unfrozen, giving rise to a tearing-mode instability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A general theory for gauge-free lifting

A theory for lifting equations of motion for charged particle dynamics, subject to given electromagnetic like forces, up to a gauge-free system of coupled Hamiltonian Vlasov-Maxwell like equations is given. The theory provides very general expressions for the polarization and magnetization vector fields in terms of the particle dynamics description of matter. Thus, as is common in plasma physics, the particle dynamics replaces conventional constitutive relations for matter. Several examples are considered including the usual Vlasov-Maxwell theory, a guiding center kinetic theory, Vlasov-Maxwell theory with the inclusion of spin, and a Vlasov-Maxwell theory with the inclusion of Dirac's magnetic monopoles. All are shown to be Hamiltonian field theories and the Jacobi identity is proven directly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effective transport barriers in nontwist systems

In fluids and plasmas with zonal flow reversed shear, a peculiar kind of transport barrier appears in the shearless region, one that is associated with a proper route of transition to chaos. Using a symplectic nontwist maps, which model such zonal flows, we investigate these barriers. In particular, the standard nontwist map, a paradigm of nontwist systems, is used to analyze the parameter dependence of the transport through a broken shearless barrier. Varying a proper control parameter, we identify the onset of structures with high stickiness that give rise to an effective barrier near the broken shearless curve. Several diagnostics are used: stickiness is measured by escape time and transmissivity plots, while barriers that separate different regions of stickiness are identified by employing two indicators in the phase space, the finite-time rotation number and the finite-time Lyapunov exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Magnetic moment type of lifting from particle dynamics to Vlasov-Maxwell dynamics

Techniques for coordinate changes that depend on both dependent and independent variables are developed and applied to the Maxwell-Vlasov Hamiltonian theory. Particle coordinate changes with a new velocity variable dependent on the magnetic field, with spatial coordinates unchanged, are lifted to transform the noncanonical Poisson bracket and, thus, the field Hamiltonian structure of the Vlasov-Maxwell equation. Several examples are given including magnetic coordinates, where the velocity is decomposed into components parallel and perpendicular to the local magnetic field, and the case of spherical velocity coordinates. An example of the lifting procedure is performed to obtain a simplified version of gyrokinetics, where the magnetic moment is used as a coordinate and the dynamics is reduced by elimination of the electric field energy in the Hamiltonian.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonlinear Acceleration Mechanism of Collisionless Magnetic Reconnection

A mechanism for fast magnetic reconnection in collisionless plasma is studied for understanding sawtooth collapse in tokamak discharges. Nonlinear growth of the tearing mode driven by electron inertia is analytically estimated by invoking the energy principle for the first time. Decrease of potential energy in the nonlinear regime (where the island width exceeds the electron skin depth) is found to be steeper than in the linear regime, resulting in acceleration of the reconnection. Release of free energy by such ideal fluid motion leads to unsteady and strong convective flow, which theoretically corroborates the inertia-driven collapse model of the sawtooth crash [D. Biskamp and J. F. Drake, Phys. Rev. Lett. 73, 971 (1994)].

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Undamped electrostatic plasma waves

Electrostatic waves in a collision-free unmagnetized plasma of electrons with fixed ions are investigated for electron equilibrium velocity distribution functions that deviate slightly from Maxwellian. Of interest are undamped waves that are the small amplitude limit of nonlinear excitations, such as electron acoustic waves (EAWs). A deviation consisting of a small plateau, a region with zero velocity derivative over a width that is a very small fraction of the electron thermal speed, is shown to give rise to new undamped modes, which here are named {\it corner modes}. The presence of the plateau turns off Landau damping and allows oscillations with phase speeds within the plateau. These undamped waves are obtained in a wide region of the $(k,ω_{_R})$ plane ($ω_{_R}$ being the real part of the wave frequency and $k$ the wavenumber), away from the well-known `thumb curve' for Langmuir waves and EAWs based on the Maxwellian. Results of nonlinear Vlasov-Poisson simulations that corroborate the existence of these modes are described. It is also shown that deviations caused by fattening the tail of the distribution shift roots off of the thumb curve toward lower $k$-values and chopping the tail shifts them toward higher $k$-values. In addition, a rule of thumb is obtained for assessing how the existence of a plateau shifts roots off of the thumb curve. Suggestions are made for interpreting experimental observations of electrostatic waves, such as recent ones in nonneutral plasmas.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A discontinuous Galerkin method for the Vlasov-Poisson system

A discontinuous Galerkin method for approximating the Vlasov-Poisson system of equations describing the time evolution of a collisionless plasma is proposed. The method is mass conservative and, in the case that piecewise constant functions are used as a basis, the method preserves the positivity of the electron distribution function and weakly enforces continuity of the electric field through mesh interfaces and boundary conditions. The performance of the method is investigated by computing several examples and error estimates associated system's approximation are stated. In particular, computed results are benchmarked against established theoretical results for linear advection and the phenomenon of linear Landau damping for both the Maxwell and Lorentz distributions. Moreover, two nonlinear problems are considered: nonlinear Landau damping and a version of the two-stream instability are computed. For the latter, fine scale details of the resulting long-time BGK-like state are presented. Conservation laws are examined and various comparisons to theory are made. The results obtained demonstrate that the discontinuous Galerkin method is a viable option for integrating the Vlasov-Poisson system.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Finite-time rotation number: a fast indicator for chaotic dynamical structures

Lagrangian coherent structures are effective barriers, sticky regions, that separate phase space regions of different dynamical behavior. The usual way to detect such structures is via finite-time Lyapunov exponents. We show that similar results can be obtained for single-frequency systems from finite-time rotation numbers, which are much faster to compute. We illustrate our claim by considering examples of continuous and discrete-time dynamical systems of physical interest.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gauge-Free Hamiltonian Structure of the Spin Maxwell-Vlasov Equations

We derive the gauge-free Hamiltonian structure of an extended kinetic theory, for which the intrinsic spin of the particles is taken into account. Such a semi-classical theory can be of interest for describing, e.g., strongly magnetized plasma systems. We find that it is possible to construct generalized brackets for the extended phase space, and discuss the implications of our findings.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hamiltonian Dynamics of Spatially-Homogeneous Vlasov-Einstein Systems

We introduce a new matter action principle, with a wide range of applicability, for the Vlasov equation in terms of a conjugate pair of functions. Here we apply this action principle to the study of matter in Bianchi cosmological models in general relativity. The Bianchi models are spatially-homogeneous solutions to the Einstein field equations, classified by the three-dimensional Lie algebra that describes the symmetry group of the model. The Einstein equations for these models reduce to a set of coupled ordinary differential equations. The class A Bianchi models admit a Hamiltonian formulation in which the components of the metric tensor and their time derivatives yield the canonical coordinates. The evolution of anisotropy in the vacuum Bianchi models is determined by a potential due to the curvature of the model, according to its symmetry. For illustrative purposes, we examine the evolution of anisotropy in models with Vlasov matter. The Vlasov content is further simplified by the assumption of cold, counter-streaming matter, a kind of matter that is far from thermal equilibrium and is not describable by an ordinary fluid model nor other more simplistic matter models. Qualitative differences and similarities are found in the dynamics of certain vacuum class A Bianchi models and Bianchi Type I models with cold, counter-streaming Vlasov matter potentials analogous to the curvature potentials of corresponding vacuum models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Caldeira-Leggett Model, Landau Damping, and the Vlasov-Poisson System

The Caldeira-Leggett Hamiltonian (Eq. (1) below) describes the interaction of a discrete harmonic oscillator with a continuous bath of harmonic oscillators. This system is a standard model of dissipation in macroscopic low temperature physics, and has applications to superconductors, quantum computing, and macroscopic quantum tunneling. The similarities between the Caldeira-Leggett model and the linearized Vlasov-Poisson equation are analyzed, and it is shown that the damping in the Caldeira-Leggett model is analogous to that of Landau damping in plasmas [1]. An invertible linear transformation [2, 3] is presented that converts solutions of the Caldeira-Leggett model into solutions of the linearized Vlasov-Poisson system.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Hamiltonian description of incompressible fluid ellipsoids

We construct the noncanonical Poisson bracket associated with the phase space of first order moments of the velocity field and quadratic moments of the density of a fluid with a free- boundary, constrained by the condition of incompressibility. Two methods are used to obtain the bracket, both based on Dirac's procedure for incorporating constraints. First, the Poisson bracket of moments of the unconstrained Euler equations is used to construct a Dirac bracket, with Casimir invariants corresponding to volume preservation and incompressibility. Second, the Dirac procedure is applied directly to the continuum, noncanonical Poisson bracket that describes the compressible Euler equations, and the moment reduction is applied to this bracket. When the Hamiltonian can be expressed exactly in terms of these moments, a closure is achieved and the resulting finite-dimensional Hamiltonian system provides exact solutions of Euler's equations. This is shown to be the case for the classical, incompressible Riemann ellipsoids, which have velocities that vary linearly with position and have constant density within an ellipsoidal boundary. The incompressible, noncanonical Poisson bracket differs from its counterpart for the compressible case in that it is not of Lie-Poisson form.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

A necessary and sufficient instability condition for inviscid shear flow

The linear stability of inviscid, incompressible, two-dimensional, plane parallel, shear flow was considered over a century ago by Rayleigh, Kelvin, and others. A principal result on the subject is Rayleigh's celebrated inflection point theorem {R80}, which states that for an equilibrium flow to be unstable, the equilibrium velocity profile must contain an inflection point. That is, if the velocity profile is given by $U(y)$, where $y$ is the cross-stream coordinate, then there must be a point, $y=y_I$, for which $U''(y_I)=0$. Much later, in 1950, Fjørtoft {F50} generalized the theorem by showing that, moreover, if there is one inflection point, then $U'''(y_I)/U'(y_I)<0$ is required for instability (see {Bar} for further extensions). Both Rayleigh's Theorem and Fjørtoft's subsequent generalization are necessary conditions for instability, but they are not sufficient. That is, even though an equilibrium profile may contain a vorticity minimum, it is not necessarily unstable. The point of this paper is to derive, for a large class of equilibrium velocity profiles, a condition that is necessary and sufficient for instability.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Invariants and Labels in Lie-Poisson Systems

Reduction is a process that uses symmetry to lower the order of a Hamiltonian system. The new variables in the reduced picture are often not canonical: there are no clear variables representing positions and momenta, and the Poisson bracket obtained is not of the canonical type. Specifically, we give two examples that give rise to brackets of the noncanonical Lie-Poisson form: the rigid body and the two-dimensional ideal fluid. From these simple cases, we then use the semidirect product extension of algebras to describe more complex physical systems. The Casimir invariants in these systems are examined, and some are shown to be linked to the recovery of information about the configuration of the system. We discuss a case in which the extension is not a semidirect product, namely compressible reduced MHD, and find for this case that the Casimir invariants lend partial information about the configuration of the system.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Negative Energy Modes and Gravitational Instability of Interpenetrating Fluids

We study the longitudinal instabilities of two interpenetrating fluids interacting only through gravity. When one of the constituents is of relatively low density, it is possible to have a band of unstable wave numbers well separated from those involved in the usual Jeans instability. If the initial streaming is large enough, and there is no linear instability, the indefinite sign of the free energy has the possible consequence of explosive interactions between positive and negative energy modes in the nonlinear regime. The effect of dissipation on the negative energy modes is also examined.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Exactly Conservative Integrators

Traditional numerical discretizations of conservative systems generically yield an artificial secular drift of any nonlinear invariants. In this work we present an explicit nontraditional algorithm that exactly conserves these invariants. We illustrate the general method by applying it to the three-wave truncation of the Euler equations, the Lotka--Volterra predator--prey model, and the Kepler problem. This method is discussed in the context of symplectic (phase space conserving) integration methods as well as nonsymplectic conservative methods. We comment on the application of our method to general conservative systems.

0 citations0 reviewsNo code linkedOpen access