Source author record

D. D. Holm

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

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

16 published item(s)

preprint2016arXiv

Weak dual pairs and jetlet methods for ideal incompressible fluid models in $n\geq 2$ dimensions

We review the role of dual pairs in mechanics and use them to derive particle-like solutions to regularized incompressible fluid systems. In our case we have a dual pair resulting from the action of diffeomorphisms on point particles (essentially by moving the points). We then augment our dual pair by considering the action of diffeomorphisms on Taylor series, also known as jets. The augmented weak dual pairs induce a hierarchy of particle-like solutions and conservation laws with particles carrying a copy of a jet group. We call these augmented particles jetlets. The jet groups serve as finite-dimensional models of the diffeomorphism group itself, and so the jetlet particles serve as a finite-dimensional model of the self-similarity exhibited by ideal incompressible fluids. The conservation law associated to jetlet solutions is shown to be a shadow of Kelvin's circulation theorem. Finally, we study the dynamics of infinite time particle mergers. We prove that two merging particles at the zeroth level in the hierarchy yield dynamics which asymptotically approach that of a single particle in the first level in the hierarchy. This merging behavior is then verified numerically as well as the exchange of angular momentum which must occur during a near collision of two particles. The resulting particle-like solutions suggest a new class of meshless methods which work in dimensions $n \geq 2$ and which exhibit a shadow of Kelvin's circulation theorem. More broadly, this provides one of the first finite-dimensional models of self-similarity in ideal fluids.

preprint2015arXiv

The momentum map representation of images

This paper discusses the mathematical framework for designing methods of large deformation matching (LDM) for image registration in computational anatomy. After reviewing the geometrical framework of LDM image registration methods, a theorem is proved showing that these methods may be designed by using the actions of diffeomorphisms on the image data structure to define their associated momentum representations as (cotangent lift) momentum maps. To illustrate its use, the momentum map theorem is shown to recover the known algorithms for matching landmarks, scalar images and vector fields. After briefly discussing the use of this approach for Diffusion Tensor (DT) images, we explain how to use momentum maps in the design of registration algorithms for more general data structures. For example, we extend our methods to determine the corresponding momentum map for registration using semidirect product groups, for the purpose of matching images at two different length scales. Finally, we discuss the use of momentum maps in the design of image registration algorithms when the image data is defined on manifolds instead of vector spaces.

preprint2014arXiv

A jetlet hierarchy for ideal fluid dynamics

Truncated Taylor expansions of smooth flow maps are used in Hamilton's principle to derive a multiscale Lagrangian particle representation of ideal fluid dynamics. Numerical simulations for scattering of solutions at one level of truncation are found to produce solutions at higher levels. These scattering events to higher levels in the Taylor expansion are interpreted as modeling a cascade to smaller scales.

preprint2013arXiv

A variational formulation of vertical slice models

A variational framework is defined for vertical slice models with three dimensional velocity depending only on x and z. The models that result from this framework are Hamiltonian, and have a Kelvin-Noether circulation theorem that results in a conserved potential vorticity in the slice geometry. These results are demonstrated for the incompressible Euler--Boussinesq equations with a constant temperature gradient in the $y$-direction (the Eady--Boussinesq model), which is an idealised problem used to study the formation and subsequent evolution of weather fronts. We then introduce a new compressible extension of this model. Unlike the incompressible model, the compressible model does not produce solutions that are also solutions of the three-dimensional equations, but it does reduce to the Eady--Boussinesq model in the low Mach number limit. This means that this new model can be used in asymptotic limit error testing for compressible weather models running in a vertical slice configuration.

preprint2013arXiv

Stretching and folding processes in the 3D Euler and Navier-Stokes equations

Stretching and folding dynamics in the incompressible, stratified 3D Euler and Navier-Stokes equations are reviewed in the context of the vector $\bdB = \nabla q\times\nablaθ$ where $q=\bom\cdot\nablaθ$. The variable $θ$ is the temperature and $\bdB$ satisfies $\partial_{t}\bdB = \mbox{curl}\,(\bu\times\bdB)$. These ideas are then discussed in the context of the full compressible Navier-Stokes equations where $q$ takes the two forms $q = \bom\cdot\nablaρ$ and $q = \bom\cdot\nabla(\lnρ)$.

preprint2013arXiv

Variational formulations of sound-proof models

We derive a family of ideal (nondissipative) 3D sound-proof fluid models that includes both the Lipps-Hemler anelastic approximation (AA) and the Durran pseudo-incompressible approximation (PIA). This family of models arises in the Euler-Poincaré framework involving a constrained Hamilton's principle expressed in the Eulerian fluid description. The derivation in this framework establishes the following properties of each member of the entire family: the Kelvin-Noether circulation theorem, conservation of potential vorticity on fluid parcels, a Lie-Poisson Hamiltonian formulation possessing conserved Casimirs, a conserved domain integrated energy and an associated variational principle satisfied by the equilibrium solutions. \smallskip Having set the stage with the derivations of 3D models using the constrained Hamilton's principle, we then derive the corresponding 2D vertical slice models for these sound-proof theories.

preprint2012arXiv

Quasi-conservation laws for compressible 3D Navier-Stokes flow

We formulate the quasi-Lagrangian fluid transport dynamics of mass density $ρ$ and the projection $q=\bom\cdot\nablaρ$ of the vorticity $\bom$ onto the density gradient, as determined by the 3D compressible Navier-Stokes equations for an ideal gas, although the results apply for an arbitrary equation of state. It turns out that the quasi-Lagrangian transport of $q$ cannot cross a level set of $ρ$. That is, in this formulation, level sets of $ρ$ (isopychnals) are impermeable to the transport of the projection $q$.

preprint2011arXiv

Enstrophy bounds and the range of space-time scales in the hydrostatic primitive equations

The hydrostatic primitive equations (HPE) form the basis of most numerical weather, climate and global ocean circulation models. Analytical (not statistical) methods are used to find a scaling proportional to $(Nu\,Ra\,Re)^{1/4}$ for the range of horizontal spatial sizes in HPE solutions, which is much broader than currently achievable computationally. The range of scales for the HPE is determined from an analytical bound on the time-averaged enstrophy of the horizontal circulation. This bound allows the formation of very small spatial scales, whose existence would excite unphysically large linear oscillation frequencies and gravity wave speeds.

preprint2011arXiv

Geometric dynamics of optimization

This paper investigates a family of dynamical systems arising from an evolutionary re-interpretation of certain optimal control and optimization problems. We focus particularly on the application in image registration of the theory of \emph{metamorphosis}. Metamorphosis is a means of tracking the optimal changes of shape that are necessary for registration of images with various types of data structures, without requiring that the transformations of shape be diffeomorphisms. This is a rich field whose possibilities are just beginning to be developed. In particular, metamorphosis and its related variants in the geometric approach to control and optimization can be expected to produce many exciting opportunities for new applications and analysis in geometric dynamics.

preprint2010arXiv

Invariant higher-order variational problems

We investigate higher-order geometric $k$-splines for template matching on Lie groups. This is motivated by the need to apply diffeomorphic template matching to a series of images, e.g., in longitudinal studies of Computational Anatomy. Our approach formulates Euler-Poincaré theory in higher-order tangent spaces on Lie groups. In particular, we develop the Euler-Poincaré formalism for higher-order variational problems that are invariant under Lie group transformations. The theory is then applied to higher-order template matching and the corresponding curves on the Lie group of transformations are shown to satisfy higher-order Euler-Poincaré equations. The example of SO(3) for template matching on the sphere is presented explicitly. Various cotangent bundle momentum maps emerge naturally that help organize the formulas. We also present Hamiltonian and Hamilton-Ostrogradsky Lie-Poisson formulations of the higher-order Euler-Poincaré theory for applications on the Hamiltonian side.

preprint2010arXiv

Multi-component generalizations of the {CH} equation: Geometrical Aspects, Peakons and Numerical Examples

The Lax pair formulation of the two-component Camassa-Holm equation (CH2) is generalized to produce an integrable multi-component family, CH(n,k), of equations with $n$ components and $1\le |k|\le n$ velocities. All of the members of the CH(n,k) family show fluid-dynamics properties with coherent solitons following particle characteristics. We determine their Lie-Poisson Hamiltonian structures and give numerical examples of their soliton solution behaviour. We concentrate on the CH(2,k) family with one or two velocities, including the CH(2,-1) equation in the Dym position of the CH2 hierarchy. A brief discussion of the CH(3,1) system reveals the underlying graded Lie-algebraic structure of the Hamiltonian formulation for CH(n,k) when $n\ge3$.

preprint2010arXiv

Stretching and folding diagnostics in solutions of the three-dimensional Euler and Navier-Stokes equations

Two possible diagnostics of stretching and folding (S&F) in fluid flows are discussed, based on the dynamics of the gradient of potential vorticity ($q = \bom\cdot\nablaθ$) associated with solutions of the three-dimensional Euler and Navier-Stokes equations. The vector $\bdB = \nabla q \times \nablaθ$ satisfies the same type of stretching and folding equation as that for the vorticity field $\bom $ in the incompressible Euler equations (Gibbon & Holm, 2010). The quantity $θ$ may be chosen as the potential temperature for the stratified, rotating Euler/Navier-Stokes equations, or it may play the role of a seeded passive scalar for the Euler equations alone. The first discussion of these S&F-flow diagnostics concerns a numerical test for Euler codes and also includes a connection with the two-dimensional surface quasi-geostrophic equations. The second S&F-flow diagnostic concerns the evolution of the Lamb vector $\bsD = \bom\times\bu$, which is the nonlinearity for Euler's equations apart from the pressure. The curl of the Lamb vector ($\boldsymbol{\varpi} := \bsD$) turns out to possess similar stretching and folding properties to that of the $\bdB$-vector.

preprint2010arXiv

The dynamics of the gradient of potential vorticity

The transport of the potential vorticity gradient $\bnabla{q}$ along surfaces of constant temperature $θ$ is investigated for the stratified Euler, Navier-Stokes and hydrostatic primitive equations of the oceans and atmosphere using the divergenceless flux vector $\bdB = \bnabla Q(q)\times\bnablaθ$, for any smooth function $Q(q)$. The flux $\bdB$ is shown to satisfy $$ \partial_t\bdB - {curl} (\bU\times\bdB) = - \bnabla\big[qQ'(q) {div} \bU\big]\times\bnablaθ, $$ where $\bU$ is a formal transport velocity of PV flux. While the left hand side of this expression is reminiscent of the frozen-in magnetic field flux in magnetohydrodynamics, the non-zero right hand side means that $\bdB$ is not frozen into the flow of $\bU$ when ${div} \bU \neq 0$. The result may apply to measurements of potential vorticity and potential temperature at the tropopause.

preprint2010arXiv

The gradient of potential vorticity, quaternions and an orthonormal frame for fluid particles

The gradient of potential vorticity (PV) is an important quantity because of the way PV (denoted as $q$) tends to accumulate locally in the oceans and atmospheres. Recent analysis by the authors has shown that the vector quantity $\bdB = \bnabla q\times \bnablaθ$ for the three-dimensional incompressible rotating Euler equations evolves according to the same stretching equation as for $\bom$ the vorticity and $\bB$, the magnetic field in magnetohydrodynamics (MHD). The $\bdB$-vector therefore acts like the vorticity $\bom$ in Euler's equations and the $\bB$-field in MHD. For example, it allows various analogies, such as stretching dynamics, helicity, superhelicity and cross helicity. In addition, using quaternionic analysis, the dynamics of the $\bdB$-vector naturally allow the construction of an orthonormal frame attached to fluid particles\,; this is designated as a quaternion frame. The alignment dynamics of this frame are particularly relevant to the three-axis rotations that particles undergo as they traverse regions of a flow when the PV gradient $\bnabla q$ is large.