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

Darryl D. Holm

Darryl D. Holm appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math-phmath.MPnlin.CDphysics.flu-dynnlin.SImath.DSmath.APmath.DGmath.PRnlin.PSquant-phcond-mat.stat-mechmath.NAphysics.class-phphysics.plasm-phApplicationsComputer Visioncond-mat.mes-hallcond-mat.otherhep-thmath.CAmath.OCnlin.AOphysics.chem-ph

Catalog footprint

What is connected

46works
24topics
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

46 published item(s)

preprint2022arXiv

Solution properties of the incompressible Euler system with rough path advection

The present paper aims to establish the local well-posedness of Euler's fluid equations on geometric rough paths. In particular, we consider the Euler equations for the incompressible flow of an ideal fluid whose Lagrangian transport velocity possesses an additional rough-in-time, divergence-free vector field. In recent work, we have demonstrated that this system can be derived from Clebsch and Hamilton-Pontryagin variational principles that possess a perturbative geometric rough path Lie-advection constraint. In this paper, we prove the local well-posedness of the system in $L^2$-Sobolev spaces $H^m$ with integer regularity $m\ge \lfloor d/2\rfloor+2$ and establish a Beale-Kato-Majda (BKM) blow-up criterion in terms of the $L^1_tL^\infty_x$-norm of the vorticity. In dimension two, we show that the $L^p$-norms of the vorticity are conserved, which yields global well-posedness and a Wong-Zakai approximation theorem for the stochastic version of the equation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Variational principles for fluid dynamics on rough paths

In this paper, we introduce a new framework for parametrization schemes (PS) in GFD. Using the theory of controlled rough paths, we derive a class of rough geophysical fluid dynamics (RGFD) models as critical points of rough action functionals. These RGFD models characterize Lagrangian trajectories in fluid dynamics as geometric rough paths (GRP) on the manifold of diffeomorphic maps. Three constrained variational approaches are formulated for the derivation of these models. The first is the Clebsch formulation, in which the constraints are imposed as rough advection laws. The second is the Hamilton-Pontryagin formulation, in which the constraints are imposed as right-invariant rough vector fields. The third is the Euler--Poincaré formulation in which the variations are constrained. These variational principles lead directly to the Lie--Poisson Hamiltonian formulation of fluid dynamics on geometric rough paths. The GRP framework preserves the geometric structure of fluid dynamics obtained by using Lie group reduction to pass from Lagrangian to Eulerian variational principles, thereby yielding a rough formulation of the Kelvin circulation theorem. The rough-path variational approach includes non-Markovian perturbations of the Lagrangian fluid trajectories. In particular, memory effects can be introduced through this formulation through a judicious choice of the rough path (e.g. a realization of a fractional Brownian motion). In the special case when the rough path is a realization of a semimartingale, we recover the SGFD models in Holm (2015). However, by eliminating the need for stochastic variational tools, we retain a pathwise interpretation of the Lagrangian trajectories. In contrast, the Lagrangian trajectories in the stochastic framework are described by stochastic integrals which do not have a pathwise interpretation. Thus, the rough path formulation restores this property.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A Particle Filter for Stochastic Advection by Lie Transport (SALT): A case study for the damped and forced incompressible 2D Euler equation

In this work, we combine a stochastic model reduction with a particle filter augmented with tempering and jittering, and apply the combined algorithm to a damped and forced incompressible 2D Euler dynamics defined on a simply connected bounded domain. We show that using the combined algorithm, we are able to assimilate data from a reference system state (the ``truth") modelled by a highly resolved numerical solution of the flow that has roughly $3.1\times10^6$ degrees of freedom, into a stochastic system having two orders of magnitude less degrees of freedom, which is able to approximate the true state reasonably accurately for $5$ large scale eddy turnover times, using modest computational hardware. The model reduction is performed through the introduction of a stochastic advection by Lie transport (SALT) model as the signal on a coarser resolution. The SALT approach was introduced as a general theory using a geometric mechanics framework from Holm, Proc. Roy. Soc. A (2015). This work follows on the numerical implementation for SALT presented by Cotter et al, SIAM Multiscale Model. Sim. (2019) for the flow in consideration. The model reduction is substantial: The reduced SALT model has $4.9\times 10^4$ degrees of freedom. Results from reliability tests on the assimilated system are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Modelling uncertainty using stochastic transport noise in a 2-layer quasi-geostrophic model

The stochastic variational approach for geophysical fluid dynamics was introduced by Holm (Proc Roy Soc A, 2015) as a framework for deriving stochastic parameterisations for unresolved scales. This paper applies the variational stochastic parameterisation in a two-layer quasi-geostrophic model for a beta-plane channel flow configuration. We present a new method for estimating the stochastic forcing (used in the parameterisation) to approximate unresolved components using data from the high resolution deterministic simulation, and describe a procedure for computing physically-consistent initial conditions for the stochastic model. We also quantify uncertainty of coarse grid simulations relative to the fine grid ones in homogeneous (teamed with small-scale vortices) and heterogeneous (featuring horizontally elongated large-scale jets) flows, and analyse how the spread of stochastic solutions depends on different parameters of the model. The parameterisation is tested by comparing it with the true eddy-resolving solution that has reached some statistical equilibrium and the deterministic solution modelled on a low-resolution grid. The results show that the proposed parameterisation significantly depends on the resolution of the stochastic model and gives good ensemble performance for both homogeneous and heterogeneous flows, and the parameterisation lays solid foundations for data assimilation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A Geometric Diffuse-Interface Method for Droplet Spreading

This paper exploits the theory of geometric gradient flows to introduce an alternative regularization of the thin-film equation. The solution properties of this regularization are investigated via a sequence of numerical simulations whose results lead to new perspectives on thin-film behavior. The new perspectives in large-scale droplet-spreading dynamics are elucidated by comparing numerical-simulation results for the solution properties of the current model with corresponding known properties of three different alternative models. The three specific comparisons in solution behavior are made with the slip model, the precursor-film method and the diffuse-interface model.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Lagrangian averaged stochastic advection by Lie transport for fluids

We formulate a class of stochastic partial differential equations based on Kelvin's circulation theorem for ideal fluids. In these models, the velocity field is randomly transported by white-noise vector fields, as well as by its own average over realizations of this noise. We call these systems the Lagrangian averaged stochastic advection by Lie transport (LA SALT) equations. These equations are nonlinear and non-local, in both physical and probability space. Before taking this average, the equations recover the Stochastic Advection by Lie Transport (SALT) fluid equations introduced by Holm (2015). Remarkably, the introduction of the non-locality in probability space in the form of momentum transported by its own mean velocity gives rise to a closed equation for the expectation field which comprises Navier--Stokes equations with Lie--Laplacian "dissipation". As such, this form of non-locality provides a regularization mechanism. The formalism we develop is closely connected to the stochastic Weber velocity framework of Constantin and Iyer (2008) in the case when the noise correlates are taken to be the constant basis vectors in $\mathbb{R}^3$ and, thus, the Lie--Laplacian reduces to the usual Laplacian. We extend this class of equations to allow for advected quantities to be present and affect the flow through exchange of kinetic and potential energies. The statistics of the solutions for the LA SALT fluid equations are found to be changing dynamically due to an array of intricate correlations among the physical variables. The statistical properties of the LA SALT physical variables propagate as local evolutionary equations which when spatially integrated become dynamical equations for the variances of the fluctuations. Essentially, the LA SALT theory is a non-equilibrium stochastic linear response theory for fluctuations in SALT fluids with advected quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Stochastic Large Deformation Model for Computational Anatomy

In the study of shapes of human organs using computational anatomy, variations are found to arise from inter-subject anatomical differences, disease-specific effects, and measurement noise. This paper introduces a stochastic model for incorporating random variations into the Large Deformation Diffeomorphic Metric Mapping (LDDMM) framework. By accounting for randomness in a particular setup which is crafted to fit the geometrical properties of LDDMM, we formulate the template estimation problem for landmarks with noise and give two methods for efficiently estimating the parameters of the noise fields from a prescribed data set. One method directly approximates the time evolution of the variance of each landmark by a finite set of differential equations, and the other is based on an Expectation-Maximisation algorithm. In the second method, the evaluation of the data likelihood is achieved without registering the landmarks, by applying bridge sampling using a stochastically perturbed version of the large deformation gradient flow algorithm. The method and the estimation algorithms are experimentally validated on synthetic examples and shape data of human corpora callosa.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Noise and dissipation in rigid body motion

Using the rigid body as an example, we illustrate some features of stochastic geometric mechanics. These features include: i) a geometric variational motivation for the noise structure involving Lie-Poisson brackets and momentum maps, ii) stochastic coadjoint motion with double bracket dissipation, iii) the Lie-Poisson Fokker-Planck description and its stationary solutions, iv) random dynamical systems, random attractors and SRB measures connected to statistical physics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Covariant un-reduction for curve matching

The process of un-reduction, a sort of reversal of reduction by the Lie group symmetries of a variational problem, is explored in the setting of field theories. This process is applied to the problem of curve matching in the plane, when the curves depend on more than one independent variable. This situation occurs in a variety of instances such as matching of surfaces or comparison of evolution between species. A discussion of the appropriate Lagrangian involved in the variational principle is given, as well as some initial numerical investigations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Un-reduction

This paper provides a full geometric development of a new technique called un-reduction, for dealing with dynamics and optimal control problems posed on spaces that are unwieldy for numerical implementation. The technique, which was originally concieved for an application to image dynamics, uses Lagrangian reduction by symmetry in reverse. A deeper understanding of un-reduction leads to new developments in image matching which serve to illustrate the mathematical power of the technique.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variational Principles for Stochastic Fluid Dynamics

This paper derives stochastic partial differential equations (SPDEs) for fluid dynamics from a stochastic variational principle (SVP). The Legendre transform of the Lagrangian formulation of these SPDEs yields their Lie-Poisson Hamiltonian form. The paper proceeds by: taking variations in the SVP to derive stochastic Stratonovich fluid equations; writing their Itô representation; and then investigating the properties of these stochastic fluid models in comparison with each other, and with the corresponding deterministic fluid models. The circulation properties of the stochastic Stratonovich fluid equations are found to closely mimic those of the deterministic ideal fluid models. As with deterministic ideal flows, motion along the stochastic Stratonovich paths also preserves the helicity of the vortex field lines in incompressible stochastic flows. However, these Stratonovich properties are not apparent in the equivalent Itô representation, because they are disguised by the quadratic covariation drift term arising in the Stratonovich to Itô transformation. This term is a geometric generalisation of the quadratic covariation drift term already found for scalar densities in Stratonovich's famous 1966 paper. The paper also derives motion equations for two examples of stochastic geophysical fluid dynamics (SGFD); namely, the Euler-Boussinesq and quasigeostropic approximations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Higher order Lagrange-Poincaré and Hamilton-Poincaré reductions

Motivated by the problem of longitudinal data assimilation, e.g., in the registration of a sequence of images, we develop the higher-order framework for Lagrangian and Hamiltonian reduction by symmetry in geometric mechanics. In particular, we obtain the reduced variational principles and the associated Poisson brackets. The special case of higher order Euler-Poincaré and Lie-Poisson reduction is also studied in detail.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Matrix G-Strands

We discuss three examples in which one may extend integrable Euler--Poincaré ODEs to integrable Euler--Poincaré PDEs in the matrix G-Strand context. After describing matrix G-Strand examples for $SO(3)$ and $SO(4)$ we turn our attention to $SE(3)$ where the matrix G-Strand equations recover the exact rod theory in the convective representation. We then find a zero curvature representation (ZCR) of these equations and establish the conditions under which they are completely integrable. Thus, the G-Strand equations turn out to be a rich source of integrable systems. The treatment is meant to be expository and most concepts are explained in examples in the language of vectors in $\mathbb{R}^3$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

$G$-Strands and Peakon Collisions on ${\rm Diff}(\mathbb{R})$

A $G$-strand is a map $g:\mathbb{R}\times\mathbb{R}\to G$ for a Lie group $G$ that follows from Hamilton's principle for a certain class of $G$-invariant Lagrangians. Some $G$-strands on finite-dimensional groups satisfy 1+1 space-time evolutionary equations that admit soliton solutions as completely integrable Hamiltonian systems. For example, the ${\rm SO}(3)$-strand equations may be regarded physically as integrable dynamics for solitons on a continuous spin chain. Previous work has shown that $G$-strands for diffeomorphisms on the real line possess solutions with singular support (e.g. peakons). This paper studies collisions of such singular solutions of $G$-strands when $G={\rm Diff}(\mathbb{R})$ is the group of diffeomorphisms of the real line $\mathbb{R}$, for which the group product is composition of smooth invertible functions. In the case of peakon-antipeakon collisions, the solution reduces to solving either Laplace's equation or the wave equation (depending on a sign in the Lagrangian) and is written in terms of their solutions. We also consider the complexified systems of $G$-strand equations for $G={\rm Diff}(\mathbb{R})$ corresponding to a harmonic map $g: \mathbb{C}\to{\rm Diff}(\mathbb{R})$ and find explicit expressions for its peakon-antipeakon solutions, as well.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Euler-Poincaré equations for $G$-Strands

The $G$-strand equations for a map $\mathbb{R}\times \mathbb{R}$ into a Lie group $G$ are associated to a $G$-invariant Lagrangian. The Lie group manifold is also the configuration space for the Lagrangian. The $G$-strand itself is the map $g(t,s): \mathbb{R}\times \mathbb{R}\to G$, where $t$ and $s$ are the independent variables of the $G$-strand equations. The Euler-Poincaré reduction of the variational principle leads to a formulation where the dependent variables of the $G$-strand equations take values in the corresponding Lie algebra $\mathfrak{g}$ and its co-algebra, $\mathfrak{g}^*$ with respect to the pairing provided by the variational derivatives of the Lagrangian. We review examples of different $G$-strand constructions, including matrix Lie groups and diffeomorphism group. In some cases the $G$-strand equations are completely integrable 1+1 Hamiltonian systems that admit soliton solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Inexact trajectory planning and inverse problems in the Hamilton--Pontryagin framework

We study a trajectory-planning problem whose solution path evolves by means of a Lie group action and passes near a designated set of target positions at particular times. This is a higher-order variational problem in optimal control, motivated by potential applications in computational anatomy and quantum control. Reduction by symmetry in such problems naturally summons methods from Lie group theory and Riemannian geometry. A geometrically illuminating form of the Euler-Lagrange equations is obtained from a higher-order Hamilton-Pontryagin variational formulation. In this context, the previously known node equations are recovered with a new interpretation as Legendre-Ostrogradsky momenta possessing certain conservation properties. Three example applications are discussed as well as a numerical integration scheme that follows naturally from the Hamilton-Pontryagin principle and preserves the geometric properties of the continuous-time solution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Integrable G-Strands on semisimple Lie groups

The present paper derives systems of partial differential equations that admit a quadratic zero curvature representation for an arbitrary real semisimple Lie algebra. It also determines the general form of Hamilton's principles and Hamiltonians for these systems and analyzes the linear stability of their equilibrium solutions in the examples of $\mathfrak{so}(3)$ and $\mathfrak{sl}(2,\mathbb{R})$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Interaction Dynamics of Singular Wave Fronts

Some of the most impressive singular wave fronts seen in Nature are the transbasin oceanic internal waves, which may be observed from the Space Shuttle as they propagate and interact with each other, for example, in the South China Sea. The characteristic feature of these strongly nonlinear wavefronts is that they reconnect when two of them collide transversely. We derive the EPDiff equation, and use it to model this phenomenon as elastic collisions between singular wave fronts (solitons) whose momentum is distributed along curves moving in the plane. Numerical methods for EPDiff based on compatible differencing algorithms (CDAs) are used for simulating these collisions among curves. The numerical results show the same nonlinear behavior of wavefront reconnections as that observed for internal waves in the South China Sea. We generalize the singular solutions of EPDiff for other applications, in computational anatomy and in imaging science, where the singular wavefronts are evolving image outlines, whose momentum may be distributed on surfaces moving though space in three dimensions. The key idea is always momentum exchange during collisions of the wavefronts. A suite of 2d and 3d numerical simulations provide collision rules for the wavefront reconnection phenomenon in a variety of scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the persistence properties of the cross-coupled Camassa-Holm system

In this paper we examine the evolution of solutions, that initially have compact support, of a recently-derived system of cross-coupled Camassa-Holm equations. The analytical methods which we employ provide a full picture for the persistence of compact support for the momenta. For solutions of the system itself, the answer is more convoluted, and we determine when the compactness of the support is lost, replaced instead by an exponential decay rate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Relative Geodesics in the Special Euclidean Group

We propose a notion of distance between two parametrized planar curves, called their discrepancy, and defined intuitively as the minimal amount of deformation needed to deform the source curve into the target curve. A precise definition of discrepancy is given as follows. A curve of transformations in the special Euclidean group SE(2) is said to be admissible if it maps the source curve to the target curve under the point-wise action of SE(2) on the plane. After endowing the group SE(2) with a left-invariant metric, we define a relative geodesic in SE(2) to be a critical point of the energy functional associated to the metric, over all admissible curves. The discrepancy is then defined as the value of the energy of the minimizing relative geodesic. In the first part of the paper, we derive a scalar ODE which is a necessary condition for a curve in SE(2) to be a relative geodesic, and we discuss some of the properties of the discrepancy. In the second part of the paper, we consider discrete curves, and by means of a variational principle, we derive a system of discrete equations for the relative geodesics. We finish with several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Toda lattice G-Strands

Hamilton's principle is used to extend for the Toda lattice ODEs to systems of PDEs called the Toda lattice strand equations (T-Strands). The T-Strands in the $n$-particle Toda case comprise $4n-2$ quadratically nonlinear PDEs in one space and one time variable. T-Strands form a symmetric hyperbolic Lie-Poisson Hamiltonian system of quadratically nonlinear PDEs with constant characteristic velocities. The travelling wave solutions for the two-particle T-Strand equations are solved geometrically, and their Lax pair is given to show how nonlinearity affects the solution. The three-particle T-Strands equations are also derived from Hamilton's principle. For both the two-particle and three-particle T-Strand PDEs the determining conditions for the existence of a quadratic zero-curvature relation (ZCR) exactly cancel the nonlinear terms in the PDEs. Thus, the two-particle and three-particle T-Strand PDEs do not pass the ZCR test for integrability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

$G$-Strands

A $G$-strand is a map $g(t,{s}):\,\mathbb{R}\times\mathbb{R}\to G$ for a Lie group $G$ that follows from Hamilton's principle for a certain class of $G$-invariant Lagrangians. The SO(3)-strand is the $G$-strand version of the rigid body equation and it may be regarded physically as a continuous spin chain. Here, $SO(3)_K$-strand dynamics for ellipsoidal rotations is derived as an Euler-Poincaré system for a certain class of variations and recast as a Lie-Poisson system for coadjoint flow with the same Hamiltonian structure as for a perfect complex fluid. For a special Hamiltonian, the $SO(3)_K$-strand is mapped into a completely integrable generalization of the classical chiral model for the SO(3)-strand. Analogous results are obtained for the $Sp(2)$-strand. The $Sp(2)$-strand is the $G$-strand version of the $Sp(2)$ Bloch-Iserles ordinary differential equation, whose solutions exhibit dynamical sorting. Numerical solutions show nonlinear interactions of coherent wave-like solutions in both cases. ${\rm Diff}(\mathbb{R})$-strand equations on the diffeomorphism group $G={\rm Diff}(\mathbb{R})$ are also introduced and shown to admit solutions with singular support (e.g., peakons).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collisionless kinetic theory of rolling molecules

We derive a collisionless kinetic theory for an ensemble of molecules undergoing nonholonomic rolling dynamics. We demonstrate that the existence of nonholonomic constraints leads to problems in generalizing the standard methods of statistical physics. In particular, we show that even though the energy of the system is conserved, and the system is closed in the thermodynamic sense, some fundamental features of statistical physics such as invariant measure do not hold for such nonholonomic systems. Nevertheless, we are able to construct a consistent kinetic theory using Hamilton's variational principle in Lagrangian variables, by regarding the kinetic solution as being concentrated on the constraint distribution. A cold fluid closure for the kinetic system is also presented, along with a particular class of exact solutions of the kinetic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Euler-Poincaré equations for anelastic fluid flows

We show that the ideal (nondissipative) form of the dynamical equations for the Lipps-Hemler formulation of the anelastic fluid model follow as Euler-Poincaré equations, obtained from a constrained Hamilton's principle expressed in the Eulerian fluid description. This establishes the mathematical framework for the following properties of these anelastic equations: the Kelvin-Noether circulation theorem, conservation of potential vorticity on fluid parcels, and the Lie-Poisson Hamiltonian formulation possessing conserved Casimirs, conserved domain integrated energy and an associated variational principle satisfied by the equilibrium solutions. We then introduce a modified set of anelastic equations that represent the mean anelastic motion, averaged over subgrid scale rapid fluctuations, while preserving the mathematical properties of the Euler-Poincaré framework.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Multiscale Turbulence Models Based on Convected Fluid Microstructure

The Euler-Poincaré approach to complex fluids is used to derive multiscale equations for computationally modelling Euler flows as a basis for modelling turbulence. The model is based on a \emph{kinematic sweeping ansatz} (KSA) which assumes that the mean fluid flow serves as a Lagrangian frame of motion for the fluctuation dynamics. Thus, we regard the motion of a fluid parcel on the computationally resolvable length scales as a moving Lagrange coordinate for the fluctuating (zero-mean) motion of fluid parcels at the unresolved scales. Even in the simplest 2-scale version on which we concentrate here, the contributions of the fluctuating motion under the KSA to the mean motion yields a system of equations that extends known results and appears to be suitable for modelling nonlinear backscatter (energy transfer from smaller to larger scales) in turbulence using multiscale methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantum splines

A quantum spline is a smooth curve parameterised by time in the space of unitary transformations, whose associated orbit on the space of pure states traverses a designated set of quantum states at designated times, such that the trace norm of the time rate of change of the associated Hamiltonian is minimised. The solution to the quantum spline problem is obtained, and is applied in an example that illustrates quantum control of coherent states. An efficient numerical scheme for computing quantum splines is discussed and implemented in the examples.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Euler-Poincaré formulation of hybrid plasma models

Three different hybrid Vlasov-fluid systems are derived by applying reduction by symmetry to Hamilton's variational principle. In particular, the discussion focuses on the Euler-Poincaré formulation of three major hybrid MHD models, which are compared in the same framework. These are the current-coupling scheme and two different variants of the pressure-coupling scheme. The Kelvin-Noether theorem is presented explicitly for each scheme, together with the Poincaré invariants for its hot particle trajectories. Extensions of Ertel's relation for the potential vorticity and for its gradient are also found in each case, as well as new expressions of cross helicity invariants.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Invariant higher-order variational problems II

Motivated by applications in computational anatomy, we consider a second-order problem in the calculus of variations on object manifolds that are acted upon by Lie groups of smooth invertible transformations. This problem leads to solution curves known as Riemannian cubics on object manifolds that are endowed with normal metrics. The prime examples of such object manifolds are the symmetric spaces. We characterize the class of cubics on object manifolds that can be lifted horizontally to cubics on the group of transformations. Conversely, we show that certain types of non-horizontal geodesics on the group of transformations project to cubics. Finally, we apply second-order Lagrange--Poincaré reduction to the problem of Riemannian cubics on the group of transformations. This leads to a reduced form of the equations that reveals the obstruction for the projection of a cubic on a transformation group to again be a cubic on its object manifold.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The n:m resonance dual pair

In this paper we build dual pairs of Poisson maps \[ \RR\stackrel{R_\pm}{\longleftarrow}(D,\om_\pm) \stackrel{Π_\pm}{\longrightarrow} B \] associated to $n:m$ resonance, as well as to $n:-m$ resonance. Except for the above mentioned cases $1:\pm1$, these are not pairs of momentum maps. Here $D$ is an open subset of $\CC^2$ with the above mentioned symplectic forms $\om_\pm$, and $B$ an open subset of $\RR^3$. The Poisson structure on $B$, which depends on the natural numbers $n$ and $m$, is not Lie-Poisson. Instead, its symplectic leaves are the Kummer shapes: bounded surfaces for $n:m$ resonance, and unbounded surfaces for $n:-m$ resonance

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Double bracket dissipation in kinetic theory for particles with anisotropic interactions

We derive equations of motion for the dynamics of anisotropic particles directly from the dissipative Vlasov kinetic equations, with the dissipation given by the double bracket approach (Double Bracket Vlasov, or DBV). The moments of the DBV equation lead to a nonlocal form of Darcy's law for the mass density. Next, kinetic equations for particles with anisotropic interaction are considered and also cast into the DBV form. The moment dynamics for these double bracket kinetic equations is expressed as Lie-Darcy continuum equations for densities of mass and orientation. We also show how to obtain a Smoluchowski model from a cold plasma-like moment closure of DBV. Thus, the double bracket kinetic framework serves as a unifying method for deriving different types of dynamics, from density--orientation to Smoluchowski equations. Extensions for more general physical systems are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exact geometric theory of dendronized polymer dynamics

Dendronized polymers consist of an elastic backbone with a set of iterated branch structures (dendrimers)attached at every base point of the backbone. The conformations of such molecules depend on the elastic deformation of the backbone and the branches, as well as on nonlocal (e.g., electrostatic, or Lennard-Jones) interactions between the elementary molecular units comprising the dendrimers and/or backbone. We develop a geometrically exact theory for the dynamics of such polymers, taking into account both local (elastic) and nonlocal interactions. The theory is based on applying symmetry reduction of Hamilton's principle for a Lagrangian defined on the tangent bundle of iterated semidirect products of the rotation groups that represent the relative orientations of the dendritic branches of the polymer. The resulting symmetry-reduced equations of motion are written in conservative form.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extreme events in solutions of hydrostatic and non-hydrostatic climate models

Initially this paper reviews the mathematical issues surrounding the hydrostatic (HPE) and non-hydrostatic (NPE) primitive equations that have been used extensively in numerical weather prediction and climate modelling. Cao and Titi (2005, 2007) have provided a new impetus to this by proving existence and uniqueness of solutions of viscous HPE on a cylinder with Neumann-like boundary conditions on the top and bottom. In contrast, the regularity of solutions of NPE remains an open question. With this HPE regularity result in mind, the second issue examined in this paper is whether extreme events are allowed to arise spontaneously in their solutions. Such events could include, for example, the sudden appearance and disappearance of locally intense fronts that do not involve deep convection. Analytical methods are used to show that for viscous HPE, the creation of small-scale structures is allowed locally in space and time at sizes that scale inversely with the Reynolds number.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Helical states of nonlocally interacting molecules and their linear stability: geometric approach

The equations for strands of rigid charge configurations interacting nonlocally are formulated on the special Euclidean group, SE(3), which naturally generates helical conformations. Helical stationary shapes are found by minimizing the energy for rigid charge configurations positioned along an infinitely long molecule with charges that are off-axis. The classical energy landscape for such a molecule is complex with a large number of energy minima, even when limited to helical shapes. The question of linear stability and selection of stationary shapes is studied using an SE(3) method that naturally accounts for the helical geometry. We investigate the linear stability of a general helical polymer that possesses torque-inducing non-local self-interactions and find the exact dispersion relation for the stability of the helical shapes with an arbitrary interaction potential. We explicitly determine the linearization operators and compute the numerical stability for the particular example of a linear polymer comprising a flexible rod with a repeated configuration of two equal and opposite off-axis charges, thereby showing that even in this simple case the non-local terms can induce instability that leads to the rod assuming helical shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Smooth and Peaked Solitons of the CH equation

The relations between smooth and peaked soliton solutions are reviewed for the Camassa-Holm (CH) shallow water wave equation in one spatial dimension. The canonical Hamiltonian formulation of the CH equation in action-angle variables is expressed for solitons by using the scattering data for its associated isospectral eigenvalue problem, rephrased as a Riemann-Hilbert problem. The momentum map from the action-angle scattering variables $T^*({\mathbb{T}^N})$ to the flow momentum ($\mathfrak{X}^*$) provides the Eulerian representation of the $N$-soliton solution of CH in terms of the scattering data and squared eigenfunctions of its isospectral eigenvalue problem. The dispersionless limit of the CH equation and its resulting peakon solutions are examined by using an asymptotic expansion in the dispersion parameter. The peakon solutions of the dispersionless CH equation in one dimension are shown to generalize in higher dimensions to peakon wave-front solutions of the EPDiff equation whose associated momentum is supported on smoothly embedded subspaces. The Eulerian representations of the singular solutions of both CH and EPDiff are given by the (cotangent-lift) momentum maps arising from the left action of the diffeomorphisms on smoothly embedded subspaces.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of subfilter-scale physics on regularization models

The subfilter-scale (SFS) physics of regularization models are investigated to understand the regularizations' performance as SFS models. The strong suppression of spectrally local SFS interactions and the conservation of small-scale circulation in the Lagrangian-averaged Navier-Stokes alpha-model (LANS-alpha) is found to lead to the formation of rigid bodies. These contaminate the superfilter-scale energy spectrum with a scaling that approaches k^1 as the SFS spectra is resolved. The Clark-alpha and Leray-alpha models, truncations of LANS-alpha, do not conserve small-scale circulation and do not develop rigid bodies. LANS-alpha, however, is closest to Navier-Stokes in intermittency properties. All three models are found to be stable at high Reynolds number. Differences between L^2 and H^1 norm models are clarified. For magnetohydrodynamics (MHD), the presence of the Lorentz force as a source (or sink) for circulation and as a facilitator of both spectrally nonlocal large to small scale interactions as well as local SFS interactions prevents the formation of rigid bodies in Lagrangian-averaged MHD (LAMHD-alpha). We find LAMHD-alpha performs well as a predictor of superfilter-scale energy spectra and of intermittent current sheets at high Reynolds numbers. We expect it may prove to be a generally applicable MHD-LES.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Euler's fluid equations: Optimal Control vs Optimization

An optimization method used in image-processing (metamorphosis) is found to imply Euler's equations for incompressible flow of an inviscid fluid, without requiring that the Lagrangian particle labels exactly follow the flow lines of the Eulerian velocity vector field. Thus, an optimal control problem and an optimization problem for incompressible ideal fluid flow both yield the \emph {same} Euler fluid equations, although their Lagrangian parcel dynamics are \emph{different}. This is a result of the \emph{gauge freedom} in the definition of the fluid pressure for an incompressible flow, in combination with the symmetry of fluid dynamics under relabeling of their Lagrangian coordinates. Similar ideas are also illustrated for SO(N) rigid body motion.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lagrange-Poincare field equations

The Lagrange-Poincare equations of classical mechanics are cast into a field theoretic context together with their associated constrained variational principle. An integrability/reconstruction condition is established that relates solutions of the original problem with those of the reduced problem. The Kelvin-Noether theorem is formulated in this context. Applications to the isoperimetric problem, the Skyrme model for meson interaction, metamorphosis image dynamics, and molecular strands illustrate various aspects of the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Random Hamiltonian in thermal equilibrium

A framework for the investigation of disordered quantum systems in thermal equilibrium is proposed. The approach is based on a dynamical model--which consists of a combination of a double-bracket gradient flow and a uniform Brownian fluctuation--that `equilibrates' the Hamiltonian into a canonical distribution. The resulting equilibrium state is used to calculate quenched and annealed averages of quantum observables.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Square Root Depth Wave Equations

We introduce a set of coupled equations for multilayer water waves that removes the ill-posedness of the multilayer Green-Naghdi (MGN) equations in the presence of shear. The new well-posed equations are Hamiltonian and in the absence of imposed background shear they retain the same travelling wave solutions as MGN. We call the new model the Square Root Depth equations, from the modified form of their kinetic energy of vertical motion. Our numerical results show how the Square Root Depth equations model the effects of multilayer wave propagation and interaction, with and without shear.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Hamiltonian statistical mechanics

A framework for statistical-mechanical analysis of quantum Hamiltonians is introduced. The approach is based upon a gradient flow equation in the space of Hamiltonians such that the eigenvectors of the initial Hamiltonian evolve toward those of the reference Hamiltonian. The nonlinear double-bracket equation governing the flow is such that the eigenvalues of the initial Hamiltonian remain unperturbed. The space of Hamiltonians is foliated by compact invariant subspaces, which permits the construction of statistical distributions over the Hamiltonians. In two dimensions, an explicit dynamical model is introduced, wherein the density function on the space of Hamiltonians approaches an equilibrium state characterised by the canonical ensemble. This is used to compute quenched and annealed averages of quantum observables.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Applications of Poisson Geometry to Physical Problems

We consider Lagrangians in Hamilton's principle defined on the tangent space $TG$ of a Lie group $G$. Invariance of such a Lagrangian under the action of $G$ leads to the symmetry-reduced Euler-Lagrange equations called the Euler-Poincaré equations. In this case, the invariant Lagrangian is defined on the Lie algebra of the group and its Euler-Poincaré equations are defined on the dual Lie algebra, where dual is defined by the operation of taking variational derivative. On the Hamiltonian side, the Euler-Poincaré equations are Lie-Poisson and they possess accompanying momentum maps, which encode both their conservation laws and the geometry of their solution space. The standard Euler-Poincaré examples are treated, including particle dynamics, the rigid body, the heavy top and geodesic motion on Lie groups. Additional topics deal with Fermat's principle, the $\mathbb{R}^3$ Poisson bracket, polarized optical traveling waves, deformable bodies (Riemann ellipsoids) and shallow water waves, including the integrable shallow water wave systems associated with geodesic motion on the diffeomorphisms. The lectures end with the semidirect-product Euler-Poincaré reduction theorem for ideal fluid dynamics. This theorem introduces the Euler--Poincaré variational principle for incompressible and compressible motions of ideal fluids, with applications to geophysical fluids. It also leads to their Lie-Poisson Hamiltonian formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Complex Trajectories of a Simple Pendulum

The motion of a classical pendulum in a gravitational field of strength g is explored. The complex trajectories as well as the real ones are determined. If g is taken to be imaginary, the Hamiltonian that describes the pendulum becomes PT-symmetric. The classical motion for this PT-symmetric Hamiltonian is examined in detail. The complex motion of this pendulum in the presence of an external periodic forcing term is also studied.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Soliton Dynamics in Computational Anatomy

Computational anatomy (CA) has introduced the idea of anatomical structures being transformed by geodesic deformations on groups of diffeomorphisms. Among these geometric structures, landmarks and image outlines in CA are shown to be singular solutions of a partial differential equation that is called the geodesic EPDiff equation. A recently discovered momentum map for singular solutions of EPDiff yields their canonical Hamiltonian formulation, which in turn provides a complete parameterization of the landmarks by their canonical positions and momenta. The momentum map provides an isomorphism between landmarks (and outlines) for images and singular soliton solutions of the EPDiff equation. This isomorphism suggests a new dynamical paradigm for CA, as well as new data representation.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

The Camassa-Holm hierarchy, related N-dimensional integrable systems, and algebro-geometric solution on a symplectic submanifold

This paper shows that the Camassa-Holm (CH) spectral problem yields two different integrable hierarchies of nonlinear evolution equations (NLEEs), one is of negative order CH hierachy while the other one is of positive order CH hierarchy. The two CH hierarchies possess the zero curvature representations through solving a key matrix equation. We find that the well-known CH equation is included in the negative order CH hierarchy while a Dym type equation is included in the positive order CH hierarchy. Furthermore, under two constraint conditions between the eigenfunctions and the potentials, the CH spectral problem is cast in: (enumerate) a new Neumann-like N-dimensional system when it is restricted into a symplectic submanifold of $\R^{2N}$ which is proven to be integrable by using the Dirac-Poisson bracket and the r-matrix process; and a new Bargmann-like N-dimensional system when it is considered in the whole $\R^{2N}$ which is proven to be integrable by using the standard Poisson bracket and the r-matrix process. (enumerate) In the paper, we present two $4\times4$ instead of $N\times N$ r-matrix structures. One is for the Neumann-like CH system (not the peaked CH system), while the other one is for the Bargmann-like CH system. The whole CH hierarchy (both positive and negative order) is shown to have the parametric solution which obey the constraint relation. In particular, the CH equation constrained to some symplectic submanifold, and the Dym type equation have the parametric solutions. Moreover, we see that the kind of parametric solution of the CH equation is not gauge equivalent to the peakons. Solving the parametric representation of solution on the symplectic submanifold gives a class of new algebro-geometric solution of the CH equation.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

An Integrable Shallow Water Equation with Linear and Nonlinear Dispersion

We study a class of 1+1 quadratically nonlinear water wave equations that combines the linear dispersion of the Korteweg-deVries (KdV) equation with the nonlinear/nonlocal dispersion of the Camassa-Holm (CH) equation, yet still preserves integrability via the inverse scattering transform (IST) method. This IST-integrable class of equations contains both the KdV equation and the CH equation as limiting cases. It arises as the compatibility condition for a second order isospectral eigenvalue problem and a first order equation for the evolution of its eigenfunctions. This integrable equation is shown to be a shallow water wave equation derived by asymptotic expansion at one order higher approximation than KdV. We compare its traveling wave solutions to KdV solitons.

0 citations0 reviewsNo code linkedOpen access