N-Particle Dynamics of the Euler Equations for Planar Diffeomorphisms
Total Page:16
File Type:pdf, Size:1020Kb
N-Particle dynamics of the Euler equations for planar diffeomorphisms Robert I McLachlan1 and Stephen Marsland2 1 IFS, Massey University, Palmerston North, New Zealand 2 IIS&T, Massey University, Palmerston North, New Zealand Email: r.mclachlan,s.r.marsland @massey.ac.nz { } Abstract The Euler equations associated with diffeomorphism groups have received much recent study because of their links with fluid dynamics, computer vision, and mechanics. In this paper, we consider the dynamics of N point particles or ‘blobs’ moving under the action of the Euler equations associated with the group of diffeomorphisms of the plane in a variety of different metrics. The 2 body problem is always integrable, and we analyze its phase portrait under different metrics. In particular, we show that 2-body capturing orbits (in which the distances between the particles tend to 0 as t ) can occur when the kernel is sufficiently smooth and the relative initial velocity of the→ particles ∞ is sufficiently large. We compute the dynamics of these ‘dipoles’ with respect to other test particles, and supplement the calculations with simulations for larger N that illustrate the different regimes. 1 Introduction The Euler equations for planar diffeomorphisms are m˙ + u m + uT m + m(div u)=0, (1) · ∇ ∇ · wherem ˙ denotes differentiation with respect to time, u(x, t) (u, x Rn, t R) is a velocity field, and m(x, t) its associated momentum. The velocity u and momentum∈ m∈ are related by m = u (2) A arXiv:nlin/0512030v1 [nlin.CD] 13 Dec 2005 where is an elliptic operator (e.g. = (1 2)k) called the inertia operator. ThereA has been much recent interestA in the− ∇ Euler equations (1) because they arise in several different fields, including computer vision and fluid dynamics. In computer vision they appear in two and three dimensions in the field of image registration, such as in the averaged template matching equations [13, 15] and the geodesic interpolating clamped-plate spline [12], while in fluid dynamics they are the limiting case of the shallow water wave Camassa-Holm equations in one and two dimensions [7], and the motion of an ideal incompressible fluid corresponds to geodesics on the group of volume-preserving diffeomorphisms [1, 5] Euler equations such as Eq. (1) have a natural geometric origin, which we sketch here although these details are not needed in the paper. Let G be a Lie or diffeomorphism group with Lie algebra g. Let G be equipped with a left- or right-invariant metric that restricts to a metric , on g. Typically, this metric is defined by a linear inertia operator : g g∗ via hh ii A → u, v := u, v . (3) hh ii h A i 1 The geodesic equation on T G can be reduced to give a non-canonical Hamiltonian evolution equation on g called the Euler–Poincar´e equation, or transferred (via the Legendre transform, g g∗, which in this case is u m := u) to a Lie-Poisson system on g∗ called the Euler equations→ → A ∗ m˙ = ad −1 m, (4) ± A m where the sign is + for left- and for right-invariant metrics. The variable m g∗ is called the momentum and u g the velocity−. ∈ The most famous∈ Euler equation on a Lie group is the equation of motion of a free rigid body (G = SO(3), m = body angular momentum). An infinite-dimensional example is the Landau-Lifshitz equation on the loop group G = C∞(R, SO(3)). Examples on diffeomorphism n groups include the Euler fluid equations on the group G = Diff vol(R ) of volume-preserving diffeomorphisms with respect to the L2 metric (m = vorticity) [5], the Camassa-Holm equation 1 1 n 1 (G = Diff(S ), H metric) [7], and the second-grade fluid equations (G = Diffvol(R ), H metric) [16]. There has been recent interest in the use of diffeomorphic deformations (‘warps’) to align groups of medical images. This alignment allows the range of variation in the appearance of biological structures to be measured, providing an approach to automated disease diagnosis. One approach to this problem, known as Computational Anatomy, is to introduce group actions as deformable templates that are warped via the actions of a group onto other images [6]. This work has been fundamental to a large amount of research on aligning images through landmark matching, where corresponding points are defined on a set of images, and diffeomorphic warps used to align them. The corresponding group for image analysis is the full diffeomorphism group, not the volume-preserving subgroup. Under a right-invariant Riemannian metric, it can be shown that the geodesics of the motion of a set of landmarks can be computed as an optimisation problem. See [12, 13] for an overview. For G = Diff(Rn) the Euler equations are given by Eqs. (1,2) (see [8, 15] for further details). The inverse of the inertia operator is given by convolution with the Green’s function G of , i.e., u = G m, where denotes convolutionA and G(x, x′)= δ(x x′) for x, x′ Rn. We shallA only consider∗ rotationally∗ invariant and diagonalA ; in this case −G(x, x′) = G(∈x x′ ) for a scalar function G, which we call the kernel of G.A Examples in the literature includek − Gaussiank [9] and Bessel function kernels [16] associated with the H1 metric. In this paper, we consider a family of kernels associated with Hk metrics for various k, which includes the Gaussian kernel in the limit as k . A striking feature→ ∞ of Euler equations on diffeomorphism groups is that they admit (formally, at least) exact solutions in which the momentum is concentrated at a finite set of points that we call particles. For fluid equations these are point vortices, which are widely studied both in their own right and as a means of approximating the evolution of smooth or other vorticities. At first sight it is remarkable that a PDE should have such a finite-dimensional reduction, unrelated to integrability. (Such a reduction, however, is not unique to Euler equations. On the other hand, there are many particle-based numerical methods, collectively known as Smoothed Particle Hydrodynamics [14], which are not exact solutions of the PDE they approximate.) For the 2D and 3D Euler fluid equations, convergence of the point vortex solutions to solu- tions for smooth initial data has been established [2]. The speed of convergence can be improved by smoothing out the point vortices to vortex blobs [4], even though the (e.g., Gaussian) blobs are no longer an exact solution of the Euler equations. Instead, their evolution can be regarded as that of delta-functions under a slightly different inertia operator. Altering the shape of the blob (which corresponds to the kernel) is equivalent to altering the metric. In this way we are lead to consider a wide class of metrics. For an introduction to many-body problems, particu- larly arising in Newtonian and relativistic mechanics, see [3]. Point vortices defined over Dirac delta-functions, and the regularised version, vortex blob methods, are introduced in [4]; see [10] for a review of their use. The methods used for image registration and template matching of images in [13, 12], and references therein, correspond to computing the motion of particles in Diff(R2). Various metrics 2 have been used, including Gaussians, H1, and Hk. Diffeomorphic image registration can be considered as a boundary value problem, since the start and end points of the particles are defined by the images under consideration. However, it can equivalently be treated as an initial value problem under the definition of the starting point and initial velocity of the particles [17, 12]. This paper provides an initial study of the dynamics of particles for the group Diff(R2) of all diffeomorphisms of the plane, as are used in the image registration framework. Following the approach taken in fluid dynamics, we consider a set of N particles moving under the action of a kernel G. We begin the paper by explicitly computing the Hamiltonian for the case of N particles, and computing the equations of motion in terms of the Green’s function. Following this, in section 3 we use delta-function particles and Helmholtz-style inertia operators to find the kernels Gk corresponding to Hk metrics. The limit of these, H∞, corresponds to a Gaussian kernel for delta-function particles, where the inertia operator tends to exp( ε2 2). We show how even when blobs are required, rather than particles, we can still conside− r∇ delta-function particles through a change in the inertia operator. We then shift our attention to the case N = 2, and examine the simplest possible interaction of particles. This two body problem is always integrable. We derive the reduced Hamiltonian and consider how the particle dynamics change for different metrics. The most striking feature of the dynamics is the phenomenon of particle capture, in which particles cohere or stick together. In a typical capture orbit, the distance between two particles is (e−at) as t , but there is no actual collision. We identify the regimes of capture and (moreO traditional)→ scattering ∞ orbits. For different metrics, we also show the phase portraits of the dynamics, and the appearance of typical scattering orbits. The analysis shows that the phase portrait depends sensitively on the behaviour of G close to r = 0. Only when the kernel G is smooth at r = 0 and the relative initial velocities of the particles are sufficiently large will the particles capture each other, otherwise they scatter. For the cases k 2 the Green’s function is not sufficiently smooth to typically allow capture orbits.