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

Philip J. Morrison

Philip 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-phphysics.flu-dynmath-phmath.MPnlin.CDmath.NAmath.DSphysics.class-phphysics.comp-phastro-ph.IMastro-ph.SRmath.APmath.SGnlin.PSNumerical Analysis

Catalog footprint

What is connected

30works
15topics
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

30 published item(s)

preprint2022arXiv

A Hamiltonian model for the macroscopic Maxwell equations using exterior calculus

A Hamiltonian field theory for the macroscopic Maxwell equations with fully general polarization and magnetization is stated in the language of differential forms. The precise procedure for translating the vector calculus formulation into differential forms is discussed in detail. We choose to distinguish between straight and twisted differential forms so that all integrals be taken over densities (i.e. twisted top forms). This ensures that the duality pairings, which are stated as integrals over densities, are orientation independent. The relationship between functional differentiation with respect to vector fields and with respect to differential forms is established using the chain rule. The theory is developed such that the Poisson bracket is metric and orientation independent with all metric dependence contained in the Hamiltonian. As is typically seen in the exterior calculus formulation of Maxwell's equations, the Hodge star operator plays a key role in modeling the constitutive relations. As a demonstration of the kind of constitutive models this theory accommodates, the paper concludes with several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A mimetic discretization of the macroscopic Maxwell equations in Hamiltonian form

A mimetic spectral element discretization, utilizing a novel Galerkin projection Hodge star operator, of the macroscopic Maxwell equations in Hamiltonian form is presented. The idea of splitting purely topological and metric dependent quantities is natural in the Hamiltonian modeling framework as the Poisson bracket is metric free with the Hamiltonian containing all metric information. This idea may be incorporated into the mimetic spectral element method by directly discretizing the Poincaré duality structure. This "split exterior calculus mimetic spectral element method" yields spatially discretized Maxwell's equations which are Hamiltonian and exactly and strongly conserve Gauss's laws. Moreover, the new discrete Hodge star operator is itself of interest as a partition of the purely topological and metric dependent portions of the Hodge star operator. As a simple test case, the numerical results of applying this method to a one-dimensional version of Maxwell's equations are given.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dual PIC: a structure preserving method for discretizing Lie-Poisson brackets

We consider a general discretization strategy for Hamiltonian field theories generated by Lie-Poisson brackets which we call dual PIC (DPIC). This method involves prescribing two different discrete representations of the dynamical variable which are constrained as a Casimir invariant of the flow to coincide with one another via an L2 projection throughout the entire simulation. This allows one to leverage the relative advantages of each discrete representation. We begin by describing DPIC as applied to a general Lie-Poisson system and then provide illustrative examples: the discretization of the two-dimensional vorticity equations and the Vlasov-Poisson equation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Fluid Flow Induced Deformation of a Boundary Hair

The deformation of a dense carpet of hair due to Stokes flow in a channel can be described by a nonlinear integro-differential equation for the shape of a single hair, which possesses several solutions for a given choice of parameters. While being posed in a previous study and bearing resemblance to the pendulum problem from mechanics, this equation has not been analytically solved until now. Despite the presence on an integral with a nonlinear functional dependence on the dependent variable, the system is integrable. We compare the analytically obtained solution to a finite-difference numerical approach, identify the physically realizable solution branch, and briefly study the solution structure through a conserved energy-like quantity. Time-dependent fluid-structure interactions are a rich and complex subject to investigate and we argue that the solution discussed herein can be used as a basis for understanding these systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The kinetic origin of the fluid helicity -- a symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g. that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov's Poisson algebra). The implication is that some "kinetic effect" can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a "sub-algebra" of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved, if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry, and show that, for a special class of flows (so-called epi-2 dimensional flows), the helicity symmetry is written as $\partial_γ=0$ for a coordinate $γ$ of the configuration space.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A class of three-dimensional gyroviscous magnetohydrodynamic models

A Hamiltonian and action principle formalism for deriving three-dimensional gyroviscous magnetohydrodynamic models is presented. The uniqueness of the approach in constructing the gyroviscous tensor from first principles and its ability to explain the origin of the gyromap and the gyroviscous terms are highlighted. The procedure allows for the specification of free functions, which can be used to generate a wide range of gyroviscous models. Through the process of reduction, the noncanonical Hamiltonian bracket is obtained and briefly analysed.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deformation of Lie-Poisson algebras and chirality

Linearization of a Hamiltonian system around an equilibrium point yields a set of Hamiltonian-symmetric spectra: If $λ$ is an eigenvalue of the linearized generator, $-λ$ and $\barλ$ (hence, $-\barλ$) are also eigenvalues -- the former implies a time-reversal symmetry, while the latter guarantees the reality of the solution. However, linearization around a singular equilibrium point (which commonly exists in noncanonical Hamiltonian systems) works out differently, resulting in breaking of the Hamiltonian symmetry of spectra; time-reversal asymmetry causes chirality. This interesting phenomenon was first found in analyzing the chiral motion of the rattleback, a boat-shaped top having misaligned axes of inertia and geometry [Phys. Lett. A 381 (2017), 2772--2777]. To elucidate how chiral spectra are generated, we study the 3-dimensional Lie-Poisson systems, and classify the prototypes of singularities that cause symmetry breaking. The central idea is the deformation of the underlying Lie algebra; invoking Bianchi's list of all 3-dimensional Lie algebras, we show that the so-called class-B algebras, which are produced by asymmetric deformations of the simple algebra so(3), yield chiral spectra when linearized around their singularities. The theory of deformation is generalized to higher dimensions, including the infinite-dimensional Poisson manifolds relevant to fluid mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A method for Hamiltonian truncation: A four-wave example

A method for extracting finite-dimensional Hamiltonian systems from a class of 2+1 Hamiltonian mean field theories is presented. These theories possess noncanonical Poisson brackets, which normally resist Hamiltonian truncation, but a process of beatification by coordinate transformation near a reference state is described in order to perturbatively overcome this difficulty. Two examples of four-wave truncation of Euler's equation for scalar vortex dynamics are given and compared: one a direct non-Hamiltonian truncation of the equations of motion, the other obtained by beatifying the Poisson bracket and then truncating.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Concomitant Hamiltonian and topological structures of extended magnetohydrodynamics

The paper describes the unique geometric properties of ideal magnetohydrodynamics (MHD), and demonstrates how such features are inherited by extended MHD, viz. models that incorporate two-fluid effects (the Hall term and electron inertia). The generalized helicities, and other geometric expressions for these models are presented in a topological context, emphasizing their universal facets. Some of the results presented include: the generalized Kelvin circulation theorems; the existence of two Lie-dragged 2-forms; and two concomitant helicities that can be studied via the Jones polynomial, which is widely utilized in Chern-Simons theory. The ensuing commonality is traced to the existence of an underlying Hamiltonian structure for all the extended MHD models, exemplified by the presence of a unique noncanonical Poisson bracket, and its associated energy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Explicit high-order noncanonical symplectic algorithms for ideal two-fluid systems

An explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as discrete differential form fields on a fixed mesh. With the assistance of Whitney interpolating forms, this scheme preserves the gauge symmetry of the electromagnetic field, and the pressure field is naturally derived from the discrete internal energy. The whole system is solved using the Hamiltonian splitting method discovered by He et al., which was been successfully adopted in constructing symplectic particle-in-cell schemes. Because of its structure preserving and explicit nature, this algorithm is especially suitable for large-scale simulations for physics problems that are multi-scale and require long-term fidelity and accuracy. The algorithm is verified via two tests: studies of the dispersion relation of waves in a two-fluid plasma system and the oscillating two-stream instability.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Internal wave pressure, velocity, and energy flux from density perturbations

Determination of energy transport is crucial for understanding the energy budget and fluid circulation in density varying fluids such as the ocean and the atmosphere. However, it is rarely possible to determine the energy flux field $\mathbf{J} = p \mathbf{u}$, which requires simultaneous measurements of the pressure and velocity perturbation fields, $p$ and $\mathbf{u}$. We present a method for obtaining the instantaneous $\mathbf{J}(x,z,t)$ from density perturbations alone: a Green's function-based calculation yields $p$, and $\mathbf{u}$ is obtained by integrating the continuity equation and the incompressibility condition. We validate our method with results from Navier-Stokes simulations: the Green's function method is applied to the density perturbation field from the simulations, and the result for $\mathbf{J}$ is found to agree typically to within $1\%$ with $\mathbf{J}$ computed directly using $p$ and $ \mathbf{u}$ from the Navier-Stokes simulation. We also apply the Green's function method to density perturbation data from laboratory schlieren measurements of internal waves in a stratified fluid, and the result for $\mathbf{J}$ agrees to within $6\%$ with results from Navier-Stokes simulations. Our method for determining the instantaneous velocity, pressure, and energy flux fields applies to any system described by a linear approximation of the density perturbation field, e.g., to small amplitude lee waves and propagating vertical modes. The method can be applied using our Matlab graphical user interface EnergyFlux.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability Boundaries and Sufficient Stability Conditions for Stably Stratified, Monotonic Shear Flows

Linear stability of inviscid, parallel, and stably stratified shear flow is studied under the assumption of smooth strictly monotonic profiles of shear flow and density, so that the local Richardson number is positive everywhere. The marginally unstable modes are systematically found by solving a one-parameter family of regular Sturm-Liouville problems, which can determine the stability boundaries more efficiently than solving the Taylor-Goldstein equation directly. By arguing for the non-existence of a marginally unstable mode, we derive new sufficient conditions for stability, which generalize the Rayleigh-Fjørtoft criterion for unstratified shear flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Explosive magnetic reconnection caused by an X-shaped current-vortex layer in a collisionless plasma

A mechanism for explosive magnetic reconnection is investigated by analyzing the nonlinear evolution of a collisionless tearing mode in a two-fluid model that includes the effects of electron inertia and temperature. These effects cooperatively enable a fast reconnection by forming an X-shaped current-vortex layer centered at the reconnection point. A high-resolution simulation of this model for an unprecedentedly small electron skin depth $d_e$ and ion-sound gyroradius $ρ_s$, satisfying $d_e=ρ_s$, shows an explosive tendency for nonlinear growth of the tearing mode, where it is newly found that the explosive widening of the X-shaped layer occurs locally around the reconnection point with the length of the X shape being shorter than the domain length and the wavelength of the linear tearing mode. The reason for the onset of this locally enhanced reconnection is explained theoretically by developing a novel nonlinear and nonequilibrium inner solution that models the local X-shaped layer, and then matching it to an outer solution that is approximated by a linear tearing eigenmode with a shorter wavelength than the domain length. This theoretical model proves that the local reconnection can release the magnetic energy more efficiently than the global one and the estimated scaling of the explosive growth rate agrees well with the simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Energy-Casimir stability of hybrid Vlasov-MHD models

Different variants of hybrid kinetic-fluid models are considered for describing the interaction of a bulk fluid plasma obeying MHD and an energetic component obeying a kinetic theory. Upon using the Vlasov kinetic theory for energetic particles, two planar Vlasov-MHD models are compared in terms of their stability properties. This is made possible by the Hamiltonian structures underlying the considered hybrid systems, whose infinite number of invariants makes the energy-Casimir method effective for determining stability. Equilibrium equations for the models are obtained from a variational principle and in particular a generalized hybrid Grad-Shafranov equation follows for one of the considered models. The stability conditions are then derived and discussed with particular emphasis on kinetic particle effects on classical MHD stability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Hybrid Vlasov-MHD models: Hamiltonian vs. non-Hamiltonian

This paper investigates hybrid kinetic-MHD models, where a hot plasma (governed by a kinetic theory) interacts with a fluid bulk (governed by MHD). Different nonlinear coupling schemes are reviewed, including the pressure-coupling scheme (PCS) used in modern hybrid simulations. This latter scheme suffers from being non-Hamiltonian and to not exactly conserve total energy. Upon adopting the Vlasov description for the hot component, the non-Hamiltonian PCS and a Hamiltonian variant are compared. Special emphasis is given to the linear stability of Alfvén waves, for which it is shown that a spurious instability appears at high frequency in the non-Hamiltonian version. This instability is removed in the Hamiltonian version.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Variational Necessary and Sufficient Stability Conditions for Inviscid Shear Flow

A necessary and sufficient condition for linear stability of inviscid parallel shear flow is formulated by developing a novel variational principle, where the velocity profile is assumed to be monotonic and analytic. It is shown that unstable eigenvalues of Rayleigh's equation (which is a non-selfadjoint eigenvalue problem) can be associated with positive eigenvalues of a certain selfadjoint operator. The stability is therefore determined by maximizing a quadratic form, which is theoretically and numerically more tractable than directly solving Rayleigh's equation. This variational stability criterion is based on the understandings of Krein signature for continuous spectrum and is applicable to other stability problems of infinite-dimensional Hamiltonian systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Continuum Hamiltonian Hopf Bifurcation I

Hamiltonian bifurcations in the context of noncanonical Hamiltonian matter models are described. First, a large class of 1 + 1 Hamiltonian multi-fluid models is considered. These models have linear dynamics with discrete spectra, when linearized about homogeneous equilibria, and these spectra have counterparts to the steady state and Hamiltonian Hopf bifurcations when equilibrium parameters are varied. Examples of fluid sound waves and plasma and gravitational streaming are treated in detail. Next, using these 1 + 1 examples as a guide, a large class of 2 + 1 Hamiltonian systems is introduced, and Hamiltonian bifurcations with continuous spectra are examined. It is shown how to attach a signature to such continuous spectra, which facilitates the description of the continuous Hamiltonian Hopf bifurcation. This chapter lays the groundwork for Kreĭn-like theorems associated with the CHH bifurcation that are more rigorously discussed in our companion chapter \cite{chaptII}.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Discontinuous Galerkin Methods for the Vlasov-Maxwell Equations

Discontinuous Galerkin methods are developed for solving the Vlasov-Maxwell system, methods that are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in general are hard to achieve within other numerical method frameworks for simulating the Vlasov-Maxwell system. The proposed scheme employs discontinuous Galerkin discretizations for both the Vlasov and the Maxwell equations, resulting in a consistent description of the distribution function and electromagnetic fields. It is proven, up to some boundary effects, that charge is conserved and the total energy can be preserved with suitable choices of the numerical flux for the Maxwell equations and the underlying approximation spaces. Error estimates are established for several flux choices. The scheme is tested on the streaming Weibel instability: the order of accuracy and conservation properties of the proposed method are verified.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy stability analysis for a hybrid fluid-kinetic plasma model

In plasma physics, a hybrid fluid-kinetic model is composed of a magnetohydrodynamics (MHD) part that describes a bulk fluid component and a Vlasov kinetic theory part that describes an energetic plasma component. While most hybrid models in the plasma literature are non-Hamiltonian, this paper investigates a recent Hamiltonian variant in its two-dimensional configuration. The corresponding Hamiltonian structure is described along with its Casimir invariants. Then, the energy-Casimir method is used to derive explicit sufficient stability conditions, which imply a stable spectrum and suggest nonlinear stability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Singular Casimir Elements of the Euler Equation and Equilibrium Points

The problem of the nonequivalence of the sets of equilibrium points and energy-Casimir extremal points, which occurs in the noncanonical Hamiltonian formulation of equations describing ideal fluid and plasma dynamics, is addressed in the context of the Euler equation for an incompressible inviscid fluid. The problem is traced to a Casimir deficit, where Casimir elements constitute the center of the Lie-Poisson algebra underlying the Hamiltonian formulation, and this leads to a study of the symplectic operator defining the Poisson bracket. The kernel of the symplectic operator, for this typical example of an infinite-dimensional Hamiltonian system for media in terms of Eulerian variables, is analyzed. For two-dimensional flows, a rigorously solvable system is formulated. The nonlinearity of the Euler equation makes the symplectic operator inhomogeneous on phase space (the function space of the state variable), and it is seen that this creates a singularity where the nullity of the symplectic operator (the "dimension" of the center) changes. Singular Casimir elements stemming from this singularity are unearthed using a generalization of the functional derivative that occurs in the Poisson bracket.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stability of compressible reduced magnetohydrodynamic equilibria - analogy with magnetorotational instability

Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared with those obtained for ideal magnetohydrodynamics (MHD) from the classical δW criterion. An identification of the terms in the second variation of the free energy functional for CRMHD with those of δW is made: two destabilizing effects present for CRMHD turn out to correspond to the kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending and compressibility are also identified in the reduced model. Also, using the EC method, stability conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the condition for magnetorotational instability (MRI) for the latter model, corresponds to the condition for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of the curvature term and pressure gradient, respectively. Using the EC method, stability conditions for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, the Rayleigh hydrodynamical stability condition is recovered.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Variational Approach to Necessary and Sufficient Stability Conditions for Inviscid Shear Flow

A necessary and sufficient condition for linear stability of inviscid parallel shear flow is formulated by a novel variational method, where the velocity profile is assumed to be monotonic and analytic. Unstable eigenvalues of the Rayleigh equation are shown to be associated with positive eigenvalues of a certain selfadjoint operator. The stability is therefore simply determined by maximizing a quadratic form, which is theoretically and numerically more tractable than directly solving the Rayleigh equation. This variational approach is based on the Hamiltonian nature of the inviscid fluid and will be applicable to other hydrodynamic stability problems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gradient flows in the normal and Kähler metrics and triple bracket generated metriplectic systems

The dynamics of gradient and Hamiltonian flows with particular application to flows on adjoint orbits of a Lie group and the extension of this setting to flows on a loop group are discussed. Different types of gradient flows that arise from different metrics including the so-called normal metric on adjoint orbits of a Lie group and the Kähler metric are compared. It is discussed how a Kähler metric can arise from a complex structure induced by the Hilbert transform. Hybrid and metriplectic flows that arise when one has both Hamiltonian and gradient components are examined. A class of metriplectic systems that is generated by completely antisymmetric triple brackets is described and for finite-dimensional systems given a Lie algebraic interpretation. A variety of explicit examples of the several types of flows are given.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hamiltonian magnetohydrodynamics: symmetric formulation, Casimir invariants, and equilibrium variational principles

The noncanonical Hamiltonian formulation of magnetohydrodynamics (MHD) is used to construct variational principles for symmetric equilibrium configurations of magnetized plasma including flow. In particular, helical symmetry is considered and results on axial and translational symmetries are retrieved as special cases of the helical configurations. The symmetry condition, which allows the description in terms of a magnetic flux function, is exploited to deduce a symmetric form of the noncanonical Poisson bracket of MHD. Casimir invariants are then obtained directly from the Poisson bracket. Equilibria are obtained from an energy-Casimir principle and reduced forms of this variational principle are obtained by the elimination of algebraic constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Study of conservation and recurrence of Runge-Kutta discontinuous Galerkin schemes for Vlasov-Poisson systems

In this paper we consider Runge-Kutta discontinuous Galerkin (RKDG) schemes for Vlasov-Poisson systems that model collisionless plasmas. One-dimensional systems are emphasized. The RKDG method, originally devised to solve conservation laws, is seen to have excellent conservation properties, be readily designed for arbitrary order of accuracy, and capable of being used with a positivity-preserving limiter that guarantees positivity of the distribution functions. The RKDG solver for the Vlasov equation is the main focus, while the electric field is obtained through the classical representation by Green's function for the Poisson equation. A rigorous study of recurrence of the DG methods is presented by Fourier analysis, and the impact of different polynomial spaces and the positivity-preserving limiters on the quality of the solutions is ascertained. Several benchmark test problems, such as Landau damping, two-stream instability and the KEEN (Kinetic Electrostatic Electron Nonlinear) wave, are given.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the Hamiltonian formulation of incompressible ideal fluids and magnetohydrodynamics via Dirac's theory of constraints

The Hamiltonian structures of the incompressible ideal fluid, including entropy advection, and magnetohydrodynamics are investigated by making use of Dirac's theory of constrained Hamiltonian systems. A Dirac bracket for these systems is constructed by assuming a primary constraint of constant density. The resulting bracket is seen to naturally project onto solenoidal velocity fields.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Derivation of reduced two-dimensional fluid models via Dirac's theory of constrained Hamiltonian systems

We present a Hamiltonian derivation of a class of reduced plasma two-dimensional fluid models, an example being the Charney-Hasegawa-Mima equation. These models are obtained from the same parent Hamiltonian model, which consists of the ion momentum equation coupled to the continuity equation, by imposing dynamical constraints. It is shown that the Poisson bracket associated with these reduced models is the Dirac bracket obtained from the Poisson bracket of the parent model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mode signature and stability for a Hamiltonian model of electron temperature gradient turbulence

Stability properties and mode signature for equilibria of a model of electron temperature gradient (ETG) driven turbulence are investigated by Hamiltonian techniques. After deriving the infinite families of Casimir invariants, associated with the noncanonical Poisson bracket of the model, a sufficient condition for stability is obtained by means of the Energy-Casimir method. Mode signature is then investigated for linear motions about homogeneous equilibria. Depending on the sign of the equilibrium "translated" pressure gradient, stable equilibria can either be energy stable, i.e.\ possess definite linearized perturbation energy (Hamiltonian), or spectrally stable with the existence of negative energy modes (NEMs). The ETG instability is then shown to arise through a Kre\uın-type bifurcation, due to the merging of a positive and a negative energy mode, corresponding to two modified drift waves admitted by the system. The Hamiltonian of the linearized system is then explicitly transformed into normal form, which unambiguously defines mode signature. In particular, the fast mode turns out to always be a positive energy mode (PEM), whereas the energy of the slow mode can have either positive or negative sign.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On Krein-like theorems for noncanonical Hamiltonian systems with continuous spectra: application to Vlasov-Poisson

The notions of spectral stability and the spectrum for the Vlasov-Poisson system linearized about homogeneous equilibria, f_0(v), are reviewed. Structural stability is reviewed and applied to perturbations of the linearized Vlasov operator through perturbations of f_0. We prove that for each f_0 there is an arbitrarily small delta f_0' in W^{1,1}(R) such that f_0+delta f_0$ is unstable. When $f_0$ is perturbed by an area preserving rearrangement, f_0 will always be stable if the continuous spectrum is only of positive signature, where the signature of the continuous spectrum is defined as in previous work. If there is a signature change, then there is a rearrangement of f_0 that is unstable and arbitrarily close to f_0 with f_0' in W^{1,1}. This result is analogous to Krein's theorem for the continuous spectrum. We prove that if a discrete mode embedded in the continuous spectrum is surrounded by the opposite signature there is an infinitesimal perturbation in C^n norm that makes f_0 unstable. If f_0 is stable we prove that the signature of every discrete mode is the opposite of the continuum surrounding it.

0 citations0 reviewsNo code linkedOpen access