Momentum Maps and Measure-valued Solutions (Peakons, Filaments and Sheets) for the EPDiff Equation

We study the dynamics of measure-valued solutions of what we call the EPDiff equations, standing for the {\it Euler-Poincar\'e equations associated with the diffeomorphism group (of $\mathbb{R}^n$ or an $n$-dimensional manifold $M$)}. Our main focus will be on the case of quadratic Lagrangians; that is, on geodesic motion on the diffeomorphism group with respect to the right invariant Sobolev $H^1$ metric. The corresponding Euler-Poincar\'e (EP) equations are the EPDiff equations, which coincide with the averaged template matching equations (ATME) from computer vision and agree with the Camassa-Holm (CH) equations in one dimension. The corresponding equations for the volume preserving diffeomorphism group are the well-known LAE (Lagrangian averaged Euler) equations for incompressible fluids. We first show that the EPDiff equations are generated by a smooth vector field on the diffeomorphism group for sufficiently smooth solutions. This is analogous to known results for incompressible fluids--both the Euler equations and the LAE equations--and it shows that for sufficiently smooth solutions, the equations are well-posed for short time. In fact, numerical evidence suggests that, as time progresses, these smooth solutions break up into singular solutions which, at least in one dimension, exhibit soliton behavior. With regard to these non-smooth solutions, we study measure-valued solutions that generalize to higher dimensions the peakon solutions of the (CH) equation in one dimension. One of the main purposes of this paper is to show that many of the properties of these measure-valued solutions may be understood through the fact that their solution ansatz is a momentum map. Some additional geometry is also pointed out, for example, that this momentum map is one leg of a natural dual pair.Comment: 27 pages, 2 figures, To Alan Weinstein on the occasion of his 60th Birthda

A simulation comparison of methods for new product location

Includes bibliographical references (p. 29-31)

Discrete Mechanics and Optimal Control Applied to the Compass Gait Biped

This paper presents a methodology for generating locally optimal control policies for simple hybrid mechanical systems, and illustrates the method on the compass gait biped. Principles from discrete mechanics are utilized to generate optimal control policies as solutions of constrained nonlinear optimization problems. In the context of bipedal walking, this procedure provides a comparative measure of the suboptimality of existing control policies. Furthermore, our methodology can be used as a control design tool; to demonstrate this, we minimize the specific cost of transport of periodic orbits for the compass gait biped, both in the fully actuated and underactuated case

A Discrete Theory of Connections on Principal Bundles

Connections on principal bundles play a fundamental role in expressing the equations of motion for mechanical systems with symmetry in an intrinsic fashion. A discrete theory of connections on principal bundles is constructed by introducing the discrete analogue of the Atiyah sequence, with a connection corresponding to the choice of a splitting of the short exact sequence. Equivalent representations of a discrete connection are considered, and an extension of the pair groupoid composition, that takes into account the principal bundle structure, is introduced. Computational issues, such as the order of approximation, are also addressed. Discrete connections provide an intrinsic method for introducing coordinates on the reduced space for discrete mechanics, and provide the necessary discrete geometry to introduce more general discrete symmetry reduction. In addition, discrete analogues of the Levi-Civita connection, and its curvature, are introduced by using the machinery of discrete exterior calculus, and discrete connections.Comment: 38 pages, 11 figures. Fixed labels in figure

The Euler-Poincaré Equations in Geophysical Fluid Dynamics

Recent theoretical work has developed the Hamilton's-principle analog of Lie-Poisson Hamiltonian systems defined on semidirect products. The main theoretical results are twofold: 1. Euler-Poincaré equations (the Lagrangian analog of Lie-Poisson Hamiltonian equations) are derived for a parameter dependent Lagrangian from a general variational principle of Lagrange d'Alembert type in which variations are constrained; 2. an abstract Kelvin-Noether theorem is derived for such systems. By imposing suitable constraints on the variations and by using invariance properties of the Lagrangian, as one does for the Euler equations for the rigid body and ideal fluids, we cast several standard Eulerian models of geophysical fluid dynamics (GFD) at various levels of approximation into Euler-Poincaré form and discuss their corresponding Kelvin-Noether theorems and potential vorticity conservation laws. The various levels of GFD approximation are related by substituting a sequence of velocity decompositions and asymptotic expansions into Hamilton's principle for the Euler equations of a rotating stratified ideal incompressible fluid. We emphasize that the shared properties of this sequence of approximate ideal GFD models follow directly from their Euler-Poincaré formulations. New modifications of the Euler-Boussinesq equations and primitive equations are also proposed in which nonlinear dispersion adaptively filters high wavenumbers and thereby enhances stability and regularity without compromising either low wavenumber behavior or geophysical balances

Heteroclinic connections between periodic orbits and resonance transitions in celestial mechanics

In this paper we apply dynamical systems techniques to the problem of heteroclinic connections and resonance transitions in the planar circular restricted three-body problem. These related phenomena have been of concern for some time in topics such as the capture of comets and asteroids and with the design of trajectories for space missions such as the Genesis Discovery Mission. The main new technical result in this paper is the numerical demonstration of the existence of a heteroclinic connection between pairs of periodic orbits: one around the libration point L1 and the other around L2, with the two periodic orbits having the same energy. This result is applied to the resonance transition problem and to the explicit numerical construction of interesting orbits with prescribed itineraries. The point of view developed in this paper is that the invariant manifold structures associated to L1 and L2 as well as the aforementioned heteroclinic connection are fundamental tools that can aid in understanding dynamical channels throughout the solar system as well as transport between the "interior" and "exterior" Hill's regions and other resonant phenomena

The Hamiltonian structure of a two-dimensional rigid circular cylinder interacting dynamically with N point vortices

This paper studies the dynamical fluid plus rigid-body system consisting of a two-dimensional rigid cylinder of general cross-sectional shape interacting with N point vortices. We derive the equations of motion for this system and show that, in particular, if the vortex strengths sum to zero and the rigid-body has a circular shape, the equations are Hamiltonian with respect to a Poisson bracket structure that is the sum of the rigid body Lie–Poisson bracket on Se(2)*, the dual of the Lie algebra of the Euclidean group on the plane, and the canonical Poisson bracket for the dynamics of N point vortices in an unbounded plane. We then use this Hamiltonian structure to study the linear and nonlinear stability of the moving Föppl equilibrium solutions using the energy-Casimir method

A Nonlinear Analysis of the Averaged Euler Equations

This paper develops the geometry and analysis of the averaged Euler equations for ideal incompressible flow in domains in Euclidean space and on Riemannian manifolds, possibly with boundary. The averaged Euler equations involve a parameter $\alpha$; one interpretation is that they are obtained by ensemble averaging the Euler equations in Lagrangian representation over rapid fluctuations whose amplitudes are of order $\alpha$. The particle flows associated with these equations are shown to be geodesics on a suitable group of volume preserving diffeomorphisms, just as with the Euler equations themselves (according to Arnold's theorem), but with respect to a right invariant $H^1$ metric instead of the $L^2$ metric. The equations are also equivalent to those for a certain second grade fluid. Additional properties of the Euler equations, such as smoothness of the geodesic spray (the Ebin-Marsden theorem) are also shown to hold. Using this nonlinear analysis framework, the limit of zero viscosity for the corresponding viscous equations is shown to be a regular limit, {\it even in the presence of boundaries}.Comment: 25 pages, no figures, Dedicated to Vladimir Arnold on the occasion of his 60th birthday, Arnold Festschrift Volume 2 (in press

Statistical Theory of Asteroid Escape Rates

Transition states in phase space are identified and shown to regulate the rate of escape of asteroids temporarily captured in circumplanetary orbits. The transition states, similar to those occurring in chemical reaction dynamics, are then used to develop a statistical semianalytical theory for the rate of escape of asteroids temporarily captured by Mars. Theory and numerical simulations are found to agree to better than 1%. These calculations suggest that further development of transition state theory in celestial mechanics, as an alternative to large-scale numerical simulations, will be a fruitful approach to mass transport calculations