SINGULAR SOLITONS, MOMENTUM MAPS & COMPUTATIONAL ANATOMY
|
|
- Geraldine Hamilton
- 6 years ago
- Views:
Transcription
1 SINGULAR SOLITONS, MOMENTUM MAPS & COMPUTATIONAL ANATOMY Darryl D. Holm Los Alamos National Laboratory Mathematics Imperial College London dholm@lanl.gov, d.holm@imperial.ac.uk IMA Shape Spaces Workshop April 6, 2006 I shall speak of things... so singular in their oddity as in some manner to instruct, or at least entertain, without wearying. Lorenzo da Ponte IMA Shape Spaces Workshop, April 6,
2 Collaborators & References for this talk DDH, J. T. Rananather, A. Trouvé and L. Younes, Soliton Dynamics in Computational Anatomy. NeuroImage 23 (2004) S DDH & J. E. Marsden, Momentum maps & measure valued solutions of the Euler-Poincar e equations for the diffeomorphism group. Progr. Math. 232 (2004) DDH, J. E. Marsden and T. S. Ratiu, The Euler-Poincar e equations and semidirect products with applications to continuum theories. Adv. in Math., 137 (1998) IMA Shape Spaces Workshop, April 6,
3 More collaborators & references for this talk C. J. Cotter and DDH, Discrete momentum maps for lattice EPDiff, DDH and M. F. Staley, Interaction Dynamics of Singular Wave Fronts. In preparation, see staley/ R. Camassa and DDH, An Integrable Shallow Water Equation with Peaked Solitons. Phys. Rev. Lett. 71 (1993) Mumford D Pattern theory and vision. In Questions Mathématiques En Traitement Du Signal et de L Image, Chapter 3, pp Paris: Institut Henri Poincaré IMA Shape Spaces Workshop, April 6,
4 Problem & Approach for Computational Anatomy In a very large part of morphology, our essential task lies in the comparison of related forms rather than in the precise definition of each.... This process of comparison, of recognizing in one form a definite permutation or deformation of another,... lies within the immediate province of mathematics and finds its solution in... the Theory of Transformations. D Arcy Thompson, On Growth and Form (1917) IMA Shape Spaces Workshop, April 6,
5 Summary of our discussion today Computational Anatomy (CA) compares shapes (graphical structures) by making a geodesic deformation from one shape to the another. Among these graphical structures, landmarks and image outlines in CA are found to be singular solutions of the geodesic EPDiff equation. A momentum map for singular solutions of EPDiff yields their canonical Hamiltonian formulation, which provides a complete parameterization of the landmarks and image outlines by their canonical positions and momenta. IMA Shape Spaces Workshop, April 6,
6 Summary of discussion (cont) The momentum map provides an isomorphism between landmarks (and outlines) for images and singular (weak) solutions of EPDiff. (These are solitons in 1D.) This isomorphism provides for CA: (1) a complete and non-redundant data representation (2) a dynamical paradigm in which image outlines interact by exchange of momentum (3) methods for numerical simulation & data assimilation Euler-Poincaré theory also provides a framework for unifying and extending the various approaches in CA IMA Shape Spaces Workshop, April 6,
7 The importance of momentum for CA Completes the representation of images (momentum of cartoons) Informs template matching of the possibility of soliton-like collisions and momentum exchange in image outline interactions Encodes the subsequent deformation into the initial locus and momentum of an image outline Provides numerical simulation methods using the momentum map for right action as a data structure. Accomplishes matching and data assimilation via the adjoint linear problem for template matching, using the initial momentum as a control variable. All these momentum properties flow from the EPDiff equation IMA Shape Spaces Workshop, April 6,
8 Outline of the rest of the talk Describe the template matching variational problems of computational anatomy Motivate the EPDiff evolution equation. Describe the singular solutions for the EPDiff equation with diffeomorphism group Diff. Explain why these singular solutions exist (singular momentum map) Discuss the consequences of EPDiff for CA. Give numerical examples (by Colin J Cotter, Imperial College London) IMA Shape Spaces Workshop, April 6,
9 Cost Darryl D. Holm, Imperial College London & LANL Template matching assigns the cost for a given comparison Cost(t ϕ t ) = 1 0 l(u t ) dt as a functional defined on curves ϕ t in a Lie group with tangents dϕ t dt = u t ϕ t, I t = ϕ t I 0 (1) In the following, the function u t l(u t ) = u t 2 g is taken as a squared functional norm on the space of velocity vectors. The vector space of right invariant instantaneous velocities, u t = (dϕ t /dt) ϕ 1 t, forms the tangent space at the identity of the considered Lie group, and is isomorphic to the group s Lie algebra, denoted g. IMA Shape Spaces Workshop, April 6,
10 Problem statement Darryl D. Holm, Imperial College London & LANL Given the cost of a time-dependent deformation process, defined by Cost(t ϕ t ) = 1 most problems in CA can be formulated as: 0 u t 2 g dt, (2) Find the deformation path with minimal cost, under the constraint that it carries the template to the target. Such problems have a remarkable analogy with fluid dynamics. IMA Shape Spaces Workshop, April 6,
11 Mathematical analogy: template matching vs fluid dynamics (I) The frameworks in CA and fluid dynamics both involve a right-invariant stationary principle with action, or cost function 1 A = l(u t ) dt } 0{{} Action The main differences are:, l(u t ) = 1 2 u t 2 g }{{} Lagrangian, u t = (dϕ t /dt) ϕ 1 t }{{} Right invariant Template matching Cost function is designed for the application Optimal control problem Fluid dynamics Cost function is Kinetic energy Initial value problem IMA Shape Spaces Workshop, April 6,
12 Mathematical analogy: CA and fluids (cont1) (II) The geodesic evolution for both template matching and fluid dynamics is governed by EPDiff [HoMaRa1998, Mumford1998, Younes1998], ( t + u ) m + ( u) T m + m(div u) = 0. (3) Here u = G m, where G denotes convolution with the Green s kernel G for the operator L, when m = δl =: Lu δu Here, the operator L is symmetric and positive definite for the Cost(t ϕ t ) = 1 0 l(u t ) dt = 1 0 u t 2 g dt = with L 2 pairing, whenever u t 2 g is a norm. 1 0 u t, Lu t dt IMA Shape Spaces Workshop, April 6,
13 Mathematical analogy: CA and fluids (cont2) (III) The flows in CA and fluid dynamics both evolve under a left group action on a linear representation space, I t = ϕ t I 0 They differ in the roles of their advected quantities, a t = a 0 ϕ 1 t : Template matching Image properties are passive Fluid dynamics Advected quantities may affect fluid flows Images are only swept along Waves may form in fluids for l(u, a) for l(u) since restoring forcing depend on a (IV) Some features of CA are a bit more like fluid dynamics, e.g., metamorphoses have a semidirect-product structure. Can metamorphoses have waves? IMA Shape Spaces Workshop, April 6,
14 How EPDiff emerges in CA Darryl D. Holm, Imperial College London & LANL Choose the cost function for continuously morphing I 0 into I 1 as Cost(t ϕ t ) = 1 0 l(u t ) dt = 1 0 u t 2 g dt, where u t is the velocity of the fluid deformation at time t and u t 2 g = u t, Lu t, and L is a positive symmetric linear operator. Then, the momentum governing the process, m t = Lu t, with Green s function G : u t = G m t satisfies EPDiff. Namely, ( t + u ) m + ( u) T m + m(div u) = 0 This equation arises in both template matching and fluid dynamics, and it informs both fields of endeavor. IMA Shape Spaces Workshop, April 6,
15 More about EPDiff Darryl D. Holm, Imperial College London & LANL EPDiff may be written in a Lagrangian sense as d ( ) m dx d n dx x = 0 along dt dt = u = G m That is, EPDiff convects the one-form density of momentum as ( ) ( ) t + u m t dx t d n x t = 0, with Lie derivative u wrt velocity u. Consequently, the solution of EPDiff is, ) m t dx t d n x t = (m 0 dx 0 d n x 0 for velocity u = G m = ϕ t ϕ 1 t and ϕ t Diff. ϕ 1 t Recognizing this as coadjoint motion will be a useful step. IMA Shape Spaces Workshop, April 6,
16 Deriving EPDiff from Hamilton s principle. Euler-Poincaré Reduction starts with a right (or left) G invariant Lagrangian L : T G R on the tangent bundle of a Lie group G. Right invariance of the Lagrangian may be written as L(g(t), ġ(t)) = L(g(t)h, ġ(t)h), for all h G A G invariant Lagrangian defined on T G possesses a symmetryreduced Hamilton s principle defined on the Lie algebra g T G T G/G g. Stationarity of the symmetry-reduced Hamilton s principle yields the Euler-Poincaré equations on the dual Lie algebra g. For G = Diff, this equation is EPDiff. IMA Shape Spaces Workshop, April 6,
17 Theorem 0.1 (Euler-Poincaré Reduction) Let G be a Lie group, L : T G R a right-invariant Lagrangian, and l := L g : g R be its restriction to the Lie algebra g. For a curve g(t) G, let ξ(t) = ġ(t) g(t) 1 g. Then the following four statements are equivalent: IMA Shape Spaces Workshop, April 6,
18 (i) g(t) satisfies the Euler-Lagrange equations for Lagrangian L defined on G. (ii) The variational principle holds, for variations with fixed endpoints, δ b a L(g(t), ġ(t))dt = 0 (iii) The (right invariant) Euler-Poincaré equation holds: µ = ad δl ξ µ with µ = δξ and where ξ = ġg 1 (t). (If G = Diff, this is EPDiff.) l(ξ) := L(e, ġg 1 (t)) (iv) For an arbitrary path η(t) in g which vanishes at the endpoints, i.e., η(a) = η(b) = 0, the reduced variational principle holds on g, δ b a l(ξ(t))dt = 0, for δξ = η [ξ, η] =: η ad ξ η IMA Shape Spaces Workshop, April 6,
19 Example: Outlines & momentum measures Outline matching: Given two collections of curves c 1,..., c N and C 1,..., C N in Ω, find a time-dependent diffeomorphic process (t ϕ t ) of minimal action (or cost) such that ϕ 0 = id and ϕ 1 (c i ) = C i for i = 1,..., N. The matching problem for the image outlines seeks singular momentum solutions which naturally emerge in the computation of geodesics. IMA Shape Spaces Workshop, April 6,
20 Image outlines as Singular Momentum Solutions of EPDiff For example, in the 2D plane, EPDiff has weak singular momentum solutions that are expressed as [CaHo1993, HoSt2003, HoMa2004] N m(x, t) = P a (t, s)δ ( x Q a (t, s) ) ds, (4) a=1 s where s is a Lagrangian coordinate defined along a set of N curves in the plane moving with the flow by the equations x = Q a (t, s) and supported on the delta functions in the EPDiff solution (4). = PPT Thus, the singular momentum solutions of EPDiff represent evolving wavefronts supported on delta functions defined along curves by (4). These solutions exist in any dimension & provide CA matching for points (landmarks), curves, surfaces, in any combination. IMA Shape Spaces Workshop, April 6,
21 Lagrangian representation of the singular solutions of EPDiff Substituting the singular momentum solution formula (4) for s S into EPDiff (3), then integrating against a smooth test function implies canonical Hamiltonian equations for Lagrangian wavefronts N t Q a(s, t) = P b (s, t) G(Q a (s, t), Q b (s, t) ) ds b=1 = u(q a (s, t), t) = δh N, δp a N t P a(s, t) = P a (s, t) P b (s, t) (5) b=1 Q a (s, t) G( Q a (s, t), Q b (s, t) ) ds = δh N δq a. The s S are Lagrangian coordinates; Q a (s, t) moves with the flow. IMA Shape Spaces Workshop, April 6,
22 Canonical Hamiltonian dynamics The singular momentum solutions (4) satisfy geodesic Hamiltonian phase space dynamics (5) for the canonically conjugate vector parameters Q a (s, t) and P a (s, t) with a = 1, 2... N. The Hamiltonian is H N = 1 2 N a, b=1 ( Pa (s, t) P b (s, t) ) G ( Q a (s, t), Q b (s, t) ) ds ds IMA Shape Spaces Workshop, April 6,
23 Here is the Geometry Leading to the Numerics Basic observation that ties everything together in n dimensions: Theorem (Holm and Marsden, 2004): The singular momentum solutions T Emb(S, R n ) g : (P, Q) m define an equivariant momentum map. The embedded manifold S is the support set of the P s and Q s. The momentum map is for left action of the diffeos on S. The whole system is right invariant. Its momentum map for right action is conserved. These constructions persist for a certain class of numerical schemes They apply for every choice of norm for template matching. IMA Shape Spaces Workshop, April 6,
24 Wait a second, what is a momentum map? Give an example A momentum map J : T Q g is a Hamiltonian for canonical action of a Lie group G on phase space T Q. It is expressed in terms of the pairing, : g g R as J, ξ = p, ξ q =: q p, ξ, where (q, p) Tq Q and ξ q is the infinitesimal generator of the action of the Lie algebra element ξ g on q in the manifold Q. The standard example is ξ q = ξ q for R 3 R 3 R 3, with pairing, given by scalar product of vectors. The momentum map is then J ξ = p ξ q = q p ξ J = q p This is angular momentum, the Hamiltonian for phase-space rotations. IMA Shape Spaces Workshop, April 6,
25 How is the singular solution ansatz a momentum map? A momentum map J : T Q g satisfies the defining relation, J(α q ), ξ }{{}, : g g R On LHS, α q T q Q and ξ g. = α q, ξ Q (q) }{{}, : T q Q T q Q R On RHS, ξ Q is the infinitesimal generator of the action of Lie group G on Q (associated to Lie algebra element ξ g, the Lie algebra g of G.) G acts on T Q from the left by g Q = g Q (composition of functions) Singular solution ansatz does satisfy the momentum map formula, by ( Pi (s) δ (x Q(s)) d k s ) ξ i (x)d n x = (Q, P), ξ Q R n S which verifies that both sides are equal = S P i (s)ξ i (Q(s))d k s IMA Shape Spaces Workshop, April 6,
26 What about the momentum map for relabeling symmetry? Relabeling η : S S is given by the right action Diff(S) : Q η = Q η. The infinitesimal generator of this right action is the vector field X Emb(S,R n )(Q) = d Q η ɛ = T Q ξ dɛ ɛ=0 where ξ X is tangent to the curve η ɛ Diff at ɛ = 0. Thus for right action the infinitesimal generator needed in the momentum map formula is ξ Q (s) = T Q ξ(s) IMA Shape Spaces Workshop, April 6,
27 For right action, the momentum map formula yields: J S (Q, P), ξ = (Q, P), T Q ξ = P i (s) Qi (s) S s m ξm (s) d k s ( ) = ξ P(s) dq(s) d k s S ( = P(s) dq(s) d k s, ξ(s) ) S = P dq, ξ for the pairing of the one-form density P dq with the vector field ξ. Symmetry under relabeling implies the momentum map for right action (relabeling momentum) is conserved on the embedded surface J S (Q, P) = P(s) dq(s) IMA Shape Spaces Workshop, April 6,
28 EPDiff dynamics informs optimal control for CA CA must compare two geometric objects, and thus it is concerned with an optimal control problem. However, the initial value problem for EPDiff also has important consequences for CA applications. When matching two geometric structures, the momentum at time t=0 contains all required information for reconstructing the target from the template. This is done via Hamiltonian geodesic flow. Being canonically conjugate, the momentum has exactly the same dimension as the matched structures, so there is no redundancy. Right invariance mods out the relabeling motions from the optimal solution. This symmetry also yields a conserved momentum map. IMA Shape Spaces Workshop, April 6,
29 Besides being one-to-one, the momentum representation is defined on a linear space, being dual to the velocity vectors. This means one may, for example,: study linear instability of CA processes, take averages and apply statistics to the space of image contours. The advantage is the ease of building, sampling and estimating statistical models on a linear space. IMA Shape Spaces Workshop, April 6,
30 Conclusions so far We have: (1) identified momentum as a key concept in the representation of image data for CA and (2) discussed analogies with fluid dynamics. The fundamental idea transferring from fluid dynamics to CA is the idea of momentum maps corresponding to group actions. Relabeling by right action Emb(S, R n ) G Emb(S, R n ) is a symmetry of the template Hamiltonian, so its momentum map is conserved. (Right action also generates the steady EPDiff solutions.) Left action G Emb(S, R n ) Emb(S, R n ) is not a symmetry. Momentum map J represents the weak solution of EPDiff supported on Emb(S, R n ) whose motion is a coadjoint orbit of the left action. IMA Shape Spaces Workshop, April 6,
31 Examples Let Q 0 and Q 1 be two embeddings of S 1 in R 2 which represent two shapes, each a closed planar curve. We seek a 1-parameter family of embeddings Q(t) : S 1 [0, 1] R 2 so that Q(0) = Q 0 and Q(1) matches Q 1 (up to relabeling). Q(t) is found by minimizing the constrained norm of its velocity. To find the equation for Q we require extremal values of the action A = 0 2 u(t) 2 g d t + P(s, t) ( Q(s, t) u(q(s, t))) d t, 0 S 1 i.e. we seek time-series of vector fields u(t) which are minimized in some norm subject to the constraint that Q is advected by the flow using the Lagrange multipliers P (which we call momentum). IMA Shape Spaces Workshop, April 6,
32 The minimizing solutions are δl δu = P(s, t)δ(x Q(s, t)) d s, S 1 Ṗ(s, t) = P(s, t) u(q(s, t), t), Q(s, t) = u(q(s, t), t), subject to Q(s, 0) = Q 0 (s). Now one must seek initial momentum P(s, 0) which takes shape Q 0 (s) to shape Q 1 (s). This last problem can be solved by using a gradient algorithm, where the gradient of the residual error with respect to P(s, 0) is calculated using the adjoint. Refrain from discussing the adjoint dynamics on T g and T g IMA Shape Spaces Workshop, April 6,
33 In the movie examples we use the Variational Particle-Mesh (VPM) method (Cotter, 2005) to discretize the equations, as follows: Discretize the velocity on an Eulerian grid and approximate u there. Replace S 1 by representing the shape by a finite set of Lagrangian particles {Q β } n p β=1. Interpolate from the grid to the particles using basis functions n g n g u(q β ) = u k ψ k (Q β ), with ψ k (x) = 1, x k=1 The action for the continuous time motion on the grid then becomes 1 1 A = 0 2 u(t) 2 grid + P β ( Q β u k ψ k (Q β )) d t, β k and one can obtain a fully discrete method by discretizing the action in time. The discrete adjoint is then applied in computing the inversion. k=1 IMA Shape Spaces Workshop, April 6,
34 Acknowledgments This work was partially supported by US DOE, under contract W-7405-ENG-36 for Los Alamos National Laboratory, and Office of Science ASCAR/AMS/MICS. IMA Shape Spaces Workshop, April 6,
35 Relation to Metamorphosis Consider the action for Metamorphosis, A = 1 0 l dt = u(t) 2 g + 1 Darryl D. Holm, Imperial College London & LANL σ 2 X t X + u X 2 dt where X is an element of any vector representation space of G, u denotes Lie derivative wrt velocity u and we shall take the norms ( ) 2 to define appropriate pairings, ( ) 2 = ( ), ( ). Stationarity with respect to variations of the velocity u yields δu : Lu = Z X X g where Z X = 1 ( ) σx 2 t X + u X = δl δ( t X) and the operation is defined with any chosen pairing by Z X X, η := Z X, η X =: Z X, Xη IMA Shape Spaces Workshop, April 6,
36 Relation to Metamorphosis (cont) The diamond operation is skew symmetric X Z X, η = Z X X, η and it satisfies the chain rule under Lie derivative, Darryl D. Holm, Imperial College London & LANL ξ (Z X X), η = ( ξ Z X ) X, η + Z X ξ X, η Consequently, one finds ( t + ad u)lu = ( t + u )Lu = ( t + u )(Z X X) = 0 IMA Shape Spaces Workshop, April 6,
37 Euler Poincaré equation via Clebsch relations Stationary variations of the constrained Eulerian action, S = l(ξ, a) + v, a t + ξa dt yield the following Clebsch relations, δξ : δl δξ v a = 0 (momentum map) δa : δl δa v t ξ v = 0 δv : a t + ξ a = 0 (advection relation) The Clebsch relations will recover the Euler Poincaré equations, ( ) δl t + ad ξ δξ = δl δa a and a t + ξa = 0. by using the properties of the diamond operation. (6) IMA Shape Spaces Workshop, April 6,
38 The diamond operation is defined by v a, η v, η a = v, a η. This operation is antisymmetric, v a, η = a v, η, and satisfies the chain rule under the Lie derivative, ξ (v a), η = ( ξ v) a, η + v ( ξ a), η. These two properties of the operation and the Clebsch relations (6) together imply ( ) (v t + ) δl ξ a = δa a This manipulation recovers the EP motion equation, since δl ξ δξ = δl ad ξ δξ, for one-form densities such as δl/δξ = v a. IMA Shape Spaces Workshop, April 6,
Metamorphoses for Pattern Matching
Metamorphoses for Pattern Matching Laurent Younes Joint work with Darryl Holm and Alain Trouvé Johns Hopkins University July 23 2007 Happy Birthday Darryl ( ) Metamorphoses for Pattern Matching July 23
More informationHAPPY BIRTHDAY TONY! Chicago Blochfest June 30, 2105
HAPPY BIRTHDAY TONY! 1 CLEBSCH OPTIMAL CONTROL Tudor S. Ratiu Section de Mathématiques Ecole Polytechnique Fédérale de Lausanne, Switzerland tudor.ratiu@epfl.ch 2 PLAN OF THE PRESENTATION Symmetric representation
More informationA CRASH COURSE IN EULER-POINCARÉ REDUCTION
A CRASH COURSE IN EULER-POINCARÉ REDUCTION HENRY O. JACOBS Abstract. The following are lecture notes from lectures given at the Fields Institute during the Legacy of Jerrold E. Marsden workshop. In these
More informationThe momentum map representation of images
The momentum map representation of images M. Bruveris 1, F. Gay-Balmaz, D. D. Holm 1, and T. S. Ratiu 3 Revised version August 1, 1 arxiv:91.99v [nlin.cd] 8 Apr 15 Abstract This paper discusses the mathematical
More informationGauge Fixing and Constrained Dynamics in Numerical Relativity
Gauge Fixing and Constrained Dynamics in Numerical Relativity Jon Allen The Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation is reviewed. Gauge freedom is discussed and
More informationM3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011
M3-4-5 A16 Notes for Geometric Mechanics: Oct Nov 2011 Text for the course: Professor Darryl D Holm 25 October 2011 Imperial College London d.holm@ic.ac.uk http://www.ma.ic.ac.uk/~dholm/ Geometric Mechanics
More informationThe Averaged Fluid and EPDiff Equations: (A Deja Vu Lecture) Jerrold E. Marsden
C A L T E C H Control & Dynamical Systems The Averaged Fluid and EPDiff Equations: (A Deja Vu Lecture) Jerrold E. Marsden Control and Dynamical Systems, Caltech http://www.cds.caltech.edu/ marsden/ Boulder,
More informationGeodesic shooting on shape spaces
Geodesic shooting on shape spaces Alain Trouvé CMLA, Ecole Normale Supérieure de Cachan GDR MIA Paris, November 21 2014 Outline Riemannian manifolds (finite dimensional) Spaces spaces (intrinsic metrics)
More informationReview of Lagrangian Mechanics and Reduction
Review of Lagrangian Mechanics and Reduction Joel W. Burdick and Patricio Vela California Institute of Technology Mechanical Engineering, BioEngineering Pasadena, CA 91125, USA Verona Short Course, August
More informationarxiv: v3 [math.na] 7 Jul 2009
The variational particle-mesh method for matching curves arxiv:7.v [math.na] 7 Jul 9 C J Cotter Department of Aeronautics, Imperial College London, SW7 AZ E-mail: colin.cotter@imperial.ac.u Abstract. Diffeomorphic
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More informationThe Euler Poincaré Equations in Geophysical Fluid Dynamics
The Euler Poincaré Equations in Geophysical Fluid Dynamics Darryl D. Holm 1 Jerrold E. Marsden 2 Tudor S. Ratiu 3 1 Theoretical Division and Center for Nonlinear Studies Los Alamos National Laboratory,
More informationSeptember 21, :43pm Holm Vol 2 WSPC/Book Trim Size for 9in by 6in
1 GALILEO Contents 1.1 Principle of Galilean relativity 2 1.2 Galilean transformations 3 1.2.1 Admissible force laws for an N-particle system 6 1.3 Subgroups of the Galilean transformations 8 1.3.1 Matrix
More informationRotating Concentric Circular Peakons arxiv:nlin/ v1 [nlin.si] 5 Dec 2003
Rotating Concentric Circular Peakons arxiv:nlin/0312012v1 [nlin.si] 5 Dec 2003 Darryl D. Holm Theoretical Division and Center for Nonlinear Studies Los Alamos National Laboratory, MS B284 Los Alamos, NM
More information[#1] R 3 bracket for the spherical pendulum
.. Holm Tuesday 11 January 2011 Solutions to MSc Enhanced Coursework for MA16 1 M3/4A16 MSc Enhanced Coursework arryl Holm Solutions Tuesday 11 January 2011 [#1] R 3 bracket for the spherical pendulum
More informationBACKGROUND IN SYMPLECTIC GEOMETRY
BACKGROUND IN SYMPLECTIC GEOMETRY NILAY KUMAR Today I want to introduce some of the symplectic structure underlying classical mechanics. The key idea is actually quite old and in its various formulations
More informationDISCRETE VARIATIONAL OPTIMAL CONTROL
DISCRETE VARIATIONAL OPTIMAL CONTROL FERNANDO JIMÉNEZ, MARIN KOBILAROV, AND DAVID MARTÍN DE DIEGO Abstract. This paper develops numerical methods for optimal control of mechanical systems in the Lagrangian
More informationChapter 1. Principles of Motion in Invariantive Mechanics
Chapter 1 Principles of Motion in Invariantive Mechanics 1.1. The Euler-Lagrange and Hamilton s equations obtained by means of exterior forms Let L = L(q 1,q 2,...,q n, q 1, q 2,..., q n,t) L(q, q,t) (1.1)
More informationIntegrable evolution equations on spaces of tensor densities
Integrable evolution equations on spaces of tensor densities J. Lenells, G. Misiolek and F. Tiglay* April 11, 21 Lenells, Misiolek, Tiglay* AMS Meeting, Minnesota, April 11, 21 slide 1/22 A family of equations
More informationDiscrete Dirac Mechanics and Discrete Dirac Geometry
Discrete Dirac Mechanics and Discrete Dirac Geometry Melvin Leok Mathematics, University of California, San Diego Joint work with Anthony Bloch and Tomoki Ohsawa Geometric Numerical Integration Workshop,
More informationDirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems
Dirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems Hiroaki Yoshimura Mechanical Engineering, Waseda University Tokyo, Japan Joint Work with Jerrold E. Marsden
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationSolutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations
D. D. Holm Solutions to M3-4A16 Assessed Problems # 3 15 Dec 2010 1 Solutions of M3-4A16 Assessed Problems # 3 [#1] Exercises in exterior calculus operations Vector notation for differential basis elements:
More informationInvariant Lagrangian Systems on Lie Groups
Invariant Lagrangian Systems on Lie Groups Dennis Barrett Geometry and Geometric Control (GGC) Research Group Department of Mathematics (Pure and Applied) Rhodes University, Grahamstown 6140 Eastern Cape
More informationarxiv: v3 [cs.cv] 20 Oct 2018
STRING METHODS FOR STOCHASTIC IMAGE AND SHAPE MATCHING ALEXIS ARNAUDON, DARRYL D HOLM, AND STEFAN SOMMER arxiv:1805.06038v3 [cs.cv] 20 Oct 2018 Abstract. Matching of images and analysis of shape differences
More informationLecture I: Constrained Hamiltonian systems
Lecture I: Constrained Hamiltonian systems (Courses in canonical gravity) Yaser Tavakoli December 15, 2014 1 Introduction In canonical formulation of general relativity, geometry of space-time is given
More informationarxiv: v5 [math-ph] 1 Oct 2014
VARIATIONAL INTEGRATORS FOR UNDERACTUATED MECHANICAL CONTROL SYSTEMS WITH SYMMETRIES LEONARDO COLOMBO, FERNANDO JIMÉNEZ, AND DAVID MARTÍN DE DIEGO arxiv:129.6315v5 [math-ph] 1 Oct 214 Abstract. Optimal
More informationCocycles and stream functions in quasigeostrophic motion
Journal of Nonlinear Mathematical Physics Volume 15, Number 2 (2008), 140 146 Letter Cocycles and stream functions in quasigeostrophic motion Cornelia VIZMAN West University of Timişoara, Romania E-mail:
More informationThe Geometry of Euler s equation. Introduction
The Geometry of Euler s equation Introduction Part 1 Mechanical systems with constraints, symmetries flexible joint fixed length In principle can be dealt with by applying F=ma, but this can become complicated
More informationarxiv: v1 [cs.cv] 20 Nov 2017
STOCHASTIC METAMORPHOSIS WITH TEMPLATE UNCERTAINTIES ALEXIS ARNAUDON, DARRYL D HOLM, AND STEFAN SOMMER arxiv:1711.7231v1 [cs.cv] 2 Nov 217 Abstract. In this paper, we investigate two stochastic perturbations
More informationDIRAC COTANGENT BUNDLE REDUCTION HIROAKI YOSHIMURA JERROLD E. MARSDEN. (Communicated by Juan-Pablo Ortega)
JOURNAL OF GEOMETRIC MECHANICS doi:10.3934/jgm.2009.1.87 c American Institute of Mathematical Sciences Volume 1, Number 1, March 2009 pp. 87 158 DIRAC COTANGENT BUNDLE REDUCTION HIROAKI YOSHIMURA Applied
More informationEuler-Poincaré reduction in principal bundles by a subgroup of the structure group
Euler-Poincaré reduction in principal bundles by a subgroup of the structure group M. Castrillón López ICMAT(CSIC-UAM-UC3M-UCM) Universidad Complutense de Madrid (Spain) Joint work with Pedro L. García,
More informationA DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS
International Journal of Scientific and Research Publications Volume 5 Issue 9 September 2015 1 A DIFFERENTIAL GEOMETRIC APPROACH TO FLUID MECHANICS Mansour Hassan Mansour * M A Bashir ** * Department
More informationSub-Riemannian geometry in groups of diffeomorphisms and shape spaces
Sub-Riemannian geometry in groups of diffeomorphisms and shape spaces Sylvain Arguillère, Emmanuel Trélat (Paris 6), Alain Trouvé (ENS Cachan), Laurent May 2013 Plan 1 Sub-Riemannian geometry 2 Right-invariant
More informationNonlinear MHD Stability and Dynamical Accessibility
Nonlinear MHD Stability and Dynamical Accessibility Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University Philip J. Morrison Department of Physics and Institute
More informationarxiv:chao-dyn/ v1 31 May 1999
Lagrangian Reduction, the Euler Poincaré Equations, and Semidirect Products arxiv:chao-dyn/9906004v1 31 May 1999 Hernán Cendra Universidad Nacional del Sur 8000 Bahia Blanca, Argentina uscendra@criba.edu.ar
More informationHamiltonian flows, cotangent lifts, and momentum maps
Hamiltonian flows, cotangent lifts, and momentum maps Jordan Bell jordan.bell@gmail.com Department of Mathematics, University of Toronto April 3, 2014 1 Symplectic manifolds Let (M, ω) and (N, η) be symplectic
More informationGeometric, Variational Discretization of Continuum Theories
Geometric, Variational Discretization of Continuum Theories E. S. Gawli P. Mullen D. Pavlov J. E. Marsden M. Desbrun Physica D: Nonlinear Phenomena (40(1), pp. 174-1760, 011) Abstract This study derives
More informationInvariants and Labels in Lie Poisson Systems
Invariants and Labels in Lie Poisson Systems JEAN-LUC THIFFEAULT a AN P.J. MORRISON b Institute for Fusion Studies and epartment of Physics,University of Texas at Austin, Austin, Texas, 78712-1060 ABSTRACT:
More informationLECTURE 3 MATH 261A. Office hours are now settled to be after class on Thursdays from 12 : 30 2 in Evans 815, or still by appointment.
LECTURE 3 MATH 261A LECTURES BY: PROFESSOR DAVID NADLER PROFESSOR NOTES BY: JACKSON VAN DYKE Office hours are now settled to be after class on Thursdays from 12 : 30 2 in Evans 815, or still by appointment.
More informationBRST 2006 (jmf) 7. g X (M) X ξ X. X η = [ξ X,η]. (X θ)(η) := X θ(η) θ(x η) = ξ X θ(η) θ([ξ X,η]).
BRST 2006 (jmf) 7 Lecture 2: Symplectic reduction In this lecture we discuss group actions on symplectic manifolds and symplectic reduction. We start with some generalities about group actions on manifolds.
More informationGeometric, Variational Discretization of Continuum Theories. Evan S. Gawlik
Geometric, Variational Discretization of Continuum Theories Evan S. Gawli June 010 A Senior Thesis Presented to the CALIFORNIA INSTITUTE OF TECHNOLOGY In Partial Fulfillment of the Requirements for a Minor
More informationRESEARCH REPORT. Stochastic metamorphosis with template uncertainties. Alexis Arnaudon, Darryl D. Holm and Stefan Sommer
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING 217 www.csgb.dk RESEARCH REPORT Alexis Arnaudon, Darryl D. Holm and Stefan Sommer Stochastic metamorphosis with template uncertainties No. 13, November
More informationVariational principles and Hamiltonian Mechanics
A Primer on Geometric Mechanics Variational principles and Hamiltonian Mechanics Alex L. Castro, PUC Rio de Janeiro Henry O. Jacobs, CMS, Caltech Christian Lessig, CMS, Caltech Alex L. Castro (PUC-Rio)
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationHamiltonian Dynamics In The Theory of Abstraction
Hamiltonian Dynamics In The Theory of Abstraction Subhajit Ganguly. email: gangulysubhajit@indiatimes.com Abstract: This paper deals with fluid flow dynamics which may be Hamiltonian in nature and yet
More informationLegendre Transforms, Calculus of Varations, and Mechanics Principles
page 437 Appendix C Legendre Transforms, Calculus of Varations, and Mechanics Principles C.1 Legendre Transforms Legendre transforms map functions in a vector space to functions in the dual space. From
More informationM3/4A16 Assessed Coursework 1 Darryl Holm Due in class Thursday November 6, 2008 #1 Eikonal equation from Fermat s principle
D. D. Holm November 6, 2008 M3/416 Geom Mech Part 1 1 M3/416 ssessed Coursework 1 Darryl Holm Due in class Thursday November 6, 2008 #1 Eikonal equation from Fermat s principle #1a Prove that the 3D eikonal
More informationCommon Hamiltonian structure of the shallow water equations with horizontal temperature gradients and magnetic fields
Common Hamiltonian structure of the shallow water equations with horizontal temperature gradients and magnetic fields Paul J. Dellar OCIAM, Mathematical Institute, 24 29 St Giles, Oxford, OX1 3LB, United
More informationarxiv: v1 [nlin.cd] 17 Mar 2011
Waltzing peakons and compacton pairs in a cross-coupled Camassa-Holm equation arxiv:1103.336v1 [nlin.cd] 17 Mar 011 Colin J Cotter 1, Darryl D Holm, Rossen I Ivanov,3 and James R Percival 1 Department
More informationSymmetries in Semiclassical Mechanics
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 2, 923 930 Symmetries in Semiclassical Mechanics Oleg Yu. SHVEDOV Sub-Department of Quantum Statistics and Field Theory, Department
More informationLarge-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods
Large-scale atmospheric circulation, semi-geostrophic motion and Lagrangian particle methods Colin Cotter (Imperial College London) & Sebastian Reich (Universität Potsdam) Outline 1. Hydrostatic and semi-geostrophic
More informationDeformations of coisotropic submanifolds in symplectic geometry
Deformations of coisotropic submanifolds in symplectic geometry Marco Zambon IAP annual meeting 2015 Symplectic manifolds Definition Let M be a manifold. A symplectic form is a two-form ω Ω 2 (M) which
More informationSYMPLECTIC GEOMETRY: LECTURE 5
SYMPLECTIC GEOMETRY: LECTURE 5 LIAT KESSLER Let (M, ω) be a connected compact symplectic manifold, T a torus, T M M a Hamiltonian action of T on M, and Φ: M t the assoaciated moment map. Theorem 0.1 (The
More informationTHE GEOMETRY AND DYNAMICS OF INTERACTING RIGID BODIES AND POINT VORTICES. Joris Vankerschaver. Eva Kanso. Jerrold Marsden
JOURNAL OF GEOMETRIC MECHANICS doi:10.3934/jgm.2009.1.227 c American Institute of Mathematical Sciences Volume 1, Number 2, June 2009 pp. 227 270 THE GEOMETRY AND DYNAMICS OF INTERACTING RIGID BODIES AND
More informationLecture 1: Water waves and Hamiltonian partial differential equations
Lecture : Water waves and Hamiltonian partial differential equations Walter Craig Department of Mathematics & Statistics Waves in Flows Prague Summer School 208 Czech Academy of Science August 30 208 Outline...
More informationRegularization by noise in infinite dimensions
Regularization by noise in infinite dimensions Franco Flandoli, University of Pisa King s College 2017 Franco Flandoli, University of Pisa () Regularization by noise King s College 2017 1 / 33 Plan of
More informationMetrics on the space of shapes
Metrics on the space of shapes IPAM, July 4 David Mumford Division of Applied Math Brown University Collaborators Peter Michor Eitan Sharon What is the space of shapes? S = set of all smooth connected
More informationA LITTLE TASTE OF SYMPLECTIC GEOMETRY: THE SCHUR-HORN THEOREM CONTENTS
A LITTLE TASTE OF SYMPLECTIC GEOMETRY: THE SCHUR-HORN THEOREM TIMOTHY E. GOLDBERG ABSTRACT. This is a handout for a talk given at Bard College on Tuesday, 1 May 2007 by the author. It gives careful versions
More informationThe Geometry and Dynamics of Interacting Rigid Bodies and Point Vortices
The Geometry and Dynamics of Interacting Rigid Bodies and Point Vortices Joris Vankerschaver a,b, Eva Kanso c, Jerrold E. Marsden a a Control and Dynamical Systems California Institute of Technology MC
More informationOn the coupling between an ideal fluid and immersed particles
On the coupling between an ideal fluid and immersed particles Henry O. Jacobs d,, Tudor S. Ratiu e, Mathieu Desbrun f a Imperial College London, Mathematics Department, 18 Queen s Gate Rd., London SW72AZ,
More informationNoether Symmetries and Conserved Momenta of Dirac Equation in Presymplectic Dynamics
International Mathematical Forum, 2, 2007, no. 45, 2207-2220 Noether Symmetries and Conserved Momenta of Dirac Equation in Presymplectic Dynamics Renato Grassini Department of Mathematics and Applications
More informationDiffraction by Edges. András Vasy (with Richard Melrose and Jared Wunsch)
Diffraction by Edges András Vasy (with Richard Melrose and Jared Wunsch) Cambridge, July 2006 Consider the wave equation Pu = 0, Pu = D 2 t u gu, on manifolds with corners M; here g 0 the Laplacian, D
More informationOn some relations between Optimal Transport and Stochastic Geometric Mechanics
Title On some relations between Optimal Transport and Stochastic Geometric Mechanics Banff, December 218 Ana Bela Cruzeiro Dep. Mathematics IST and Grupo de Física-Matemática Univ. Lisboa 1 / 2 Title Based
More informationMultivariate Statistical Analysis of Deformation Momenta Relating Anatomical Shape to Neuropsychological Measures
Multivariate Statistical Analysis of Deformation Momenta Relating Anatomical Shape to Neuropsychological Measures Nikhil Singh, Tom Fletcher, Sam Preston, Linh Ha, J. Stephen Marron, Michael Wiener, and
More informationSynchro-Betatron Motion in Circular Accelerators
Outlines Synchro-Betatron Motion in Circular Accelerators Kevin Li March 30, 2011 Kevin Li Synchro-Betatron Motion 1/ 70 Outline of Part I Outlines Part I: and Model Introduction Part II: The Transverse
More informationA Lattice Approximation of Dirac Equation
A Lattice Approximation of Dirac Equation Jerzy Kijowsi and Artur Thielmann Institute for Theoretical Physics Polish Academy of Sciences Al. Lotniów 3/46, 0-668 Warsaw, Poland ABSTRACT A different point
More informationarxiv:math-ph/ v1 22 May 2003
On the global version of Euler-Lagrange equations. arxiv:math-ph/0305045v1 22 May 2003 R. E. Gamboa Saraví and J. E. Solomin Departamento de Física, Facultad de Ciencias Exactas, Universidad Nacional de
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationSYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS. 1. Introduction
SYMPLECTIC MANIFOLDS, GEOMETRIC QUANTIZATION, AND UNITARY REPRESENTATIONS OF LIE GROUPS CRAIG JACKSON 1. Introduction Generally speaking, geometric quantization is a scheme for associating Hilbert spaces
More informationOptimal Collision Avoidance and Formation Switching on Riemannian Manifolds 1
Optimal Collision Avoidance and Formation Switching on Riemannian Manifolds Jianghai Hu and Shankar Sastry {jianghai,sastry}@eecs.berkeley.edu Department of Electrical Eng. & Comp. Science University of
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationHamiltonian Field Theory
Hamiltonian Field Theory August 31, 016 1 Introduction So far we have treated classical field theory using Lagrangian and an action principle for Lagrangian. This approach is called Lagrangian field theory
More informationBalanced models in Geophysical Fluid Dynamics: Hamiltonian formulation, constraints and formal stability
Balanced models in Geophysical Fluid Dynamics: Hamiltonian formulation, constraints and formal stability Onno Bokhove 1 Introduction Most fluid systems, such as the three-dimensional compressible Euler
More informationHamiltonian aspects of fluid dynamics
Hamiltonian aspects of fluid dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS 01/29/08, 01/31/08 Joris Vankerschaver (CDS) Hamiltonian aspects of fluid dynamics 01/29/08, 01/31/08 1 / 34 Outline
More informationQuantum Theory and Group Representations
Quantum Theory and Group Representations Peter Woit Columbia University LaGuardia Community College, November 1, 2017 Queensborough Community College, November 15, 2017 Peter Woit (Columbia University)
More informationControl Systems on Lie Groups
Control Systems on Lie Groups P.S. Krishnaprasad Department of Electrical and Computer Engineering and Institute for Systems Research University of Maryland, College Park, MD 20742 1 Control Systems on
More informationThe Principle of Least Action
The Principle of Least Action Anders Svensson Abstract In this paper, the principle of least action in classical mechanics is studied. The term is used in several different contexts, mainly for Hamilton
More informationSolutions to the Hamilton-Jacobi equation as Lagrangian submanifolds
Solutions to the Hamilton-Jacobi equation as Lagrangian submanifolds Matias Dahl January 2004 1 Introduction In this essay we shall study the following problem: Suppose is a smooth -manifold, is a function,
More informationMATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS
MATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS Poisson Systems and complete integrability with applications from Fluid Dynamics E. van Groesen Dept. of Applied Mathematics University oftwente
More informationThe Toda Lattice. Chris Elliott. April 9 th, 2014
The Toda Lattice Chris Elliott April 9 th, 2014 In this talk I ll introduce classical integrable systems, and explain how they can arise from the data of solutions to the classical Yang-Baxter equation.
More informationHamiltonian Systems of Negative Curvature are Hyperbolic
Hamiltonian Systems of Negative Curvature are Hyperbolic A. A. Agrachev N. N. Chtcherbakova Abstract The curvature and the reduced curvature are basic differential invariants of the pair: Hamiltonian system,
More informationHamilton-Jacobi theory on Lie algebroids: Applications to nonholonomic mechanics. Manuel de León Institute of Mathematical Sciences CSIC, Spain
Hamilton-Jacobi theory on Lie algebroids: Applications to nonholonomic mechanics Manuel de León Institute of Mathematical Sciences CSIC, Spain joint work with J.C. Marrero (University of La Laguna) D.
More informationL 2 Geometry of the Symplectomorphism Group
University of Notre Dame Workshop on Innite Dimensional Geometry, Vienna 2015 Outline 1 The Exponential Map on D s ω(m) 2 Existence of Multiplicity of Outline 1 The Exponential Map on D s ω(m) 2 Existence
More informationSOLUTIONS TO THE GINZBURG LANDAU EQUATIONS FOR PLANAR TEXTURES IN SUPERFLUID 3 He
SOLUTIONS TO THE GINZBURG LANDAU EQUATIONS FOR PLANAR TEXTURES IN SUPERFLUID 3 He V. L. GOLO, M. I. MONASTYRSKY, AND S. P. NOVIKOV Abstract. The Ginzburg Landau equations for planar textures of superfluid
More informationConstrained shape spaces of infinite dimension
Constrained shape spaces of infinite dimension Sylvain Arguillère, Emmanuel Trélat (Paris 6), Alain Trouvé (ENS Cachan), Laurent May 2013 Goal: Adding constraints to shapes in order to better fit the modeled
More informationFrom physical assumptions to classical Hamiltonian and Lagrangian particle mechanics
From physical assumptions to classical Hamiltonian and Lagrangian particle mechanics Gabriele Carcassi, Christine A. Aidala, David John Baker and Lydia Bieri University of Michigan Foundations of Physics
More informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationThe geometry of hydrodynamic integrability
1 The geometry of hydrodynamic integrability David M. J. Calderbank University of Bath October 2017 2 What is hydrodynamic integrability? A test for integrability of dispersionless systems of PDEs. Introduced
More information1 Hamiltonian formalism
1 Hamiltonian formalism 1.1 Hamilton s principle of stationary action A dynamical system with a finite number n degrees of freedom can be described by real functions of time q i (t) (i =1, 2,..., n) which,
More informationA little taste of symplectic geometry
A little taste of symplectic geometry Timothy Goldberg Thursday, October 4, 2007 Olivetti Club Talk Cornell University 1 2 What is symplectic geometry? Symplectic geometry is the study of the geometry
More informationTHE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS
ASIAN J. MATH. c 2009 International Press Vol. 13, No. 1, pp. 007 014, March 2009 002 THE POINCARÉ RECURRENCE PROBLEM OF INVISCID INCOMPRESSIBLE FLUIDS Y. CHARLES LI Abstract. Nadirashvili presented a
More informationSymmetries, Conservation Laws and Hamiltonian Structures in Geophysical Fluid Dynamics
Symmetries, Conservation Laws and Hamiltonian Structures in Geophysical Fluid Dynamics Miguel A. Jiménez-Urias 1 Department of Oceanography University of Washington AMATH 573 1 Kinematics Canonical vs
More informationDynamics and Stability application to submerged bodies, vortex streets and vortex-body systems
Dynamics and Stability application to submerged bodies, vortex streets and vortex-body systems Eva Kanso University of Southern California CDS 140B Introduction to Dynamics February 5 and 7, 2008 Fish
More informationInstability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems
Instability of Periodic Orbits of Some Rhombus and Parallelogram Four Body Problems by Abdalla Mansur A thesis submitted to the Department of Mathematics & Statistics in conformity with the requirements
More informationWe begin our discussion of special relativity with a power point presentation, available on the website.
Special Relativity We begin our discussion of special relativity with a power point presentation, available on the website.. Spacetime From the power point presentation, you know that spacetime is a four
More informationGlobal Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds
1 Global Formulations of Lagrangian and Hamiltonian Dynamics on Embedded Manifolds By Taeyoung Lee, Melvin Leok, and N. Harris McClamroch Mechanical and Aerospace Engineering, George Washington University,
More informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
More information1 M3-4-5A16 Assessed Problems # 1: Do all three problems
D. D. Holm M3-4-5A34 Assessed Problems # 1 Due 1 Feb 2013 1 1 M3-4-5A16 Assessed Problems # 1: Do all three problems Exercise 1.1 (Quaternions in Cayley-Klein (CK) parameters). Express all of your answers
More informationLecture II: Hamiltonian formulation of general relativity
Lecture II: Hamiltonian formulation of general relativity (Courses in canonical gravity) Yaser Tavakoli December 16, 2014 1 Space-time foliation The Hamiltonian formulation of ordinary mechanics is given
More information