Researcher profile

P. J. Morrison

P. J. Morrison contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
13topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

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

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.

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.

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

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.

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.

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.

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.

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.

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.

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)].

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.

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.

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.

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.