Parallel-in-time integrators for Hamiltonian systems
|
|
- Rosalyn Cunningham
- 6 years ago
- Views:
Transcription
1 Parallel-in-time integrators for Hamiltonian systems Claude Le Bris ENPC and INRIA Visiting Professor, The University of Chicago joint work with X. Dai (Paris 6 and Chinese Academy of Sciences), F. Legoll (ENPC and INRIA) and Y. Maday (Paris 6 and Brown) SSC Seminar, October 2012 The University of Chicago p. 1
2 Parallel in time algorithm for ODEs ẋ = f(x), x R d Most algorithms are sequential in nature: x n+1 = Φ T (x n ). The parareal algorithm (Lions, Maday and Turinici, 2001) is the first efficient version of a time integrator that is parallel in time. It is based upon two integrators to propagate the system over a time T : a fine integrator F T, to advance of T by many tiny timesteps; a coarse integrator G T, to advance of T by a few large timesteps (hence much cheaper). F T = (Φ δt ) T/δt and G T = (Φ dt ) T/dT with δt dt Φ δt a single step of any one-step integrator. The University of Chicago p. 2
3 The parareal paradigm - 1 The parareal algorithm combines the two integrators as follows: first, coarse propagation that yields { } x k=0 n Then, iterate over k 0: compute jumps (in parallel): n : n, x k=0 n+1 = G T (x k=0 n ) J k n = F T (x k n) G T (x k n) sequential update to obtain { } x k+1 n : n n, x k+1 n+1 = G T(x k+1 n ) + Jn k The fine solver is called only in the parallel part of the algorithm. The University of Chicago p. 3
4 The parareal paradigm - 2 x k+1 n+1 = G T(x k+1 n ) + F T (x k n) G T (x k n) If F G, then x k+1 n+1 = F T(x k+1 n ): similar to a fine sequential integration. It turns out that, at iteration k, x k k = (F T ) (k) (x 0 ), as if the fine scheme was used sequentially all the way from x 0 up to time k T. In practice, convergence occurs much faster: only a few iterations are needed to obtain a good accuracy on a large range. The approach has been successfully tested, in the past few years, on many equations (ODEs or PDEs, e.g. heat equation). The University of Chicago p. 4
5 A simple test-case: ẋ = x (coarse and fine propagators: forward Euler) 1 1e-05 fine integrator (sequential) parareal iteration 0 (coarse integrator) parareal iteration 3 parareal iteration 6 parareal iteration 10 1e-10 1e-15 1e-20 1e-25 1e-30 1e Trajectory x k (t) ( T = 0.5, dt = T, δt = 0.005) t At iteration k = 10, trajectory is exact up to time 10 T = 5, and is accurate over a much longer range. The University of Chicago p. 5
6 Hamiltonian dynamics Consider a Hamiltonian function H(q, p), and the Hamiltonian dynamics q = H (q, p), p ṗ = H q (q, p), q Rd, p R d. Introducing x = (q, p) R 2d, recast the dynamics as ẋ = f(x) = J H(x), J = 0 Id Id 0. Separated case: H(q, p) = 1 2 pt p + V (q). Then q = p, ṗ = V (q). Verlet algorithm (explicit, second order, symmetric, symplectic): p n+1/2 = p n δt 2 V (q n), q n+1 = q n + δt p n+1/2, p n+1 = p n+1/2 δt 2 V (q n+1). The University of Chicago p. 6
7 Reversibility ẋ = f(x) = J H(x) Let ϕ t (x) be the flow of the equation, namely the solution at time t with I.C. x: dϕ t (x) = J H(ϕ t (x)), ϕ t=0 (x) = x. dt we will focus on the case H(q, p) = 1 2 pt p + V (q), for which the flow is reversible: t, where ρ(x) = ρ(q, p) = (q, p). ρ ϕ t = (ϕ t ) 1 ρ The University of Chicago p. 7
8 Reversibility ẋ = f(x) = J H(x) Let ϕ t (x) be the flow of the equation, namely the solution at time t with I.C. x: dϕ t (x) = J H(ϕ t (x)), ϕ t=0 (x) = x. dt we will focus on the case H(q, p) = 1 2 pt p + V (q), for which the flow is reversible: t, where ρ(x) = ρ(q, p) = (q, p). ρ ϕ t = (ϕ t ) 1 ρ a subset of such systems is the case of integrable reversible systems: up to a reversible change of variables (q, p) (a, θ), the dynamics reads ȧ = 0, θ = ω(a), a R d, θ T d. Thus a(q, p) R d is left invariant by the dynamics. The University of Chicago p. 7
9 Backward error analysis for ODEs exact flow ϕ δt (x) ẋ = f(x) numerical scheme exact flow Φ δt (x) ẋ = f δt (x) Exact flow of the modified equation at time δt Flow Φ δt (x) given by the numerical scheme after one-time step. The University of Chicago p. 8
10 Classical analysis for Hamiltonian systems Symmetric integrators on reversible systems: δt, Φ δt Φ δt = Id Symmetric integrator Reversible modified equation Reversible modified equation + integrable reversible system + non resonant condition Conservation of invariants (hence of energy). Another class of efficient schemes: Symplectic integrators The University of Chicago p. 9
11 Parareal scheme for Hamiltonian dynamics ẋ = f(x) = J H(x), x = (q, p) R 2d Standard parareal algorithm: x k+1 n+1 = G T(x k+1 n ) + F T (x k n) G T (x k n) with F T = (Φ δt ) T/δt and G T = (Φ dt ) T/dT with δt dt, where Φ δt one step of (say) Verlet algorithm with δt time-step. Symplectic or symmetric algorithms are known to be efficient for Hamiltonian dynamics. PROBLEM: Even if F T and G T are symplectic (or symmetric), the resulting parareal algorithm is NOT symplectic (or symmetric)... Long-term (and geometric) properties of the approach are to be clarified! The University of Chicago p. 10
12 Harmonic oscillator case: q = p, ṗ = q 1e Sequential coarse scheme Parareal it. #1 Parareal it. #4 Parareal it. #10 Sequential fine scheme e-04 1e-06 1e-08 1e t Error on the energy 1 2 q p2 The University of Chicago p. 11
13 Outline Adapt the parareal algorithm to the Hamiltonian framework by using symmetrization projection on the constant energy manifold M = { (q, p) R 2d ; H(q, p) = H 0 } The University of Chicago p. 12
14 Symmetric parareal algorithm Symmetrizing by hand a given integrator is easy: consider a numerical flow Φ T, that reads x n+1 = Φ T (x n ) consider the adjoint of Φ T, defined by Φ T = Φ 1 T then the scheme x n+1 = Φ T/2 Φ T/2 (x n) is symmetric: x n+1/2 = Φ T/2 (x n), x n+1 = Φ T/2 (x n+1/2 ) e.g. Φ T/2 (x n+1/2 ) = x n, x n+1 = Φ T/2 (x n+1/2 ) Use this observation to make a symmetric version of parareal: Φ T the standard parareal algorithm; we can write explicitly what is Φ T. Remark: it is much more difficult to make symplectic a given integrator... The University of Chicago p. 13
15 Symmetric parareal The standard form of the parareal algorithm writes x k+1 n+1 = G T(x k+1 n ) + F T (x k n) G T (x k n) ( ) It is NOT a one-step integrator on x n! The University of Chicago p. 14
16 Symmetric parareal The standard form of the parareal algorithm writes x k+1 n+1 = G T(x k+1 n ) + F T (x k n) G T (x k n) ( ) It is NOT a one-step integrator on x n! Define U n = (x 0 n, x 1 n,...,x K n ) Knowing U n, the parareal algorithm defines U n+1 : we write U n+1 = Φ T (U n ). The symmetrized form of Φ, that is Φ T/2 Φ T/2, reads U n = Φ T/2 (U n+1/2 ), U n+1 = Φ T/2 (U n+1/2 ). The University of Chicago p. 14
17 Symmetric parareal The standard form of the parareal algorithm writes x k+1 n+1 = G T(x k+1 n ) + F T (x k n) G T (x k n) ( ) It is NOT a one-step integrator on x n! Define U n = (x 0 n, x 1 n,...,x K n ) Knowing U n, the parareal algorithm defines U n+1 : we write U n+1 = Φ T (U n ). The symmetrized form of Φ, that is Φ T/2 Φ T/2, reads U n = Φ T/2 (U n+1/2 ), U n+1 = Φ T/2 (U n+1/2 ). In the specific context of the parareal algorithm ( ), this two-step iteration yields the symmetric parareal integrator. The University of Chicago p. 14
18 Symmetric parareal: some basic remarks - 1 At the initial iteration (k = 0), we use the symmetric version of G: x k=0 n+1/2 = G 1 T/2 (xk=0 n ), x k=0 n+1 = G T/2 (x k=0 n+1/2 ). we next iterate over k: [ x k+1 n+1/2 = G 1 T/2 x k+1 n F T/2 (x k n+1/2 ) + G T/2(x k n+1/2 ], ) x k+1 n+1 = G T/2 (x k+1 n+1/2 ) + F T/2(x k n+1/2 ) G T/2(x k n+1/2 ) The flow is symmetric in the following sense: if, for any k and any n, (x 0 n,...,x k n, x k+1 n ) (x 0 n+1,...,x k n+1, x k+1 n+1 ) by the previous integrator, then (x 0 n+1,...,x k n+1, x k+1 n+1 ) (x0 n,...,x k n, x k+1 n ) by the exact same algorithm reversing the time. The University of Chicago p. 15
19 Symmetric parareal: some basic remarks - 2 [ x k+1 n+1/2 = G 1 T/2 x k+1 n F T/2 (x k n+1/2 ) + G T/2(x k n+1/2 ], ) x k+1 n+1 = G T/2 (x k+1 n+1/2 ) + F T/2(x k n+1/2 ) G T/2(x k n+1/2 ) If F G, then the symmetrized form reads x k+1 n+1 = F T/2 F 1 T/2 (xk+1 n ) standard symmetrized version of the fine (= coarse) propagator. Formally taking the limit k + yields x n+1 = F T/2 F 1 T/2 (x n ). Limit of the algorithm in terms of parareal iterations: standard symmetrized form of F T. Note also that we never call F 1. All the expensive computations (involving F) are performed in parallel. The University of Chicago p. 16
20 Harmonic oscillator case: q = p, ṗ = q, symmetric algorithm Sequential coarse scheme Parareal it. #1 Parareal it. #4 Parareal it. #10 Sequential fine scheme e-04 1e-06 1e t Error on the energy 1 2 q p2 The University of Chicago p. 17
21 Parareal integrator interpretation A parareal integrator may be seen, at parareal iteration k, as an integrator of a system consisting of k + 1 identical replicas of the original system. Consider the symmetric parareal integrator, used on an integrable reversible Hamiltonian system: The first replica (k = 0) is integrated by the symmetric algorithm G T/2 G 1 T/2 : energy preservation. Replicas are noninteracting: the system of replicas is an integrable reversible system, with energy = k j=0 H(q j, p j ) At iteration k + 1, we have integrated an integrable reversible system with a symmetric algorithm: its energy is preserved. Hence, for all k, k j=0 H(q j, p j ) is preserved: the energy of each replica H(q j, p j ) is preserved. However, the non resonant condition is NOT satisfied: replicas...! The University of Chicago p. 18
22 Resonances Consider a set of k identical replicas of the original system: before time discretization, they are resonant (identical!), but uncoupled: energy preservation,... OK the parareal time discretization introduces coupling. these resonances are present in the standard version as well as in the symmetric version of the algorithm. even if the symmetric algorithm is used, these resonances may impede preservation of the invariants (energy,... ) in the long-time limit. One possibility to avoid resonances: use a timestep dt (for G T ) that depends on k. x k=0 n+1 = G (k=0) T (x k=0 n ) x k=1 n+1 = G (k=1) T (x k=1 n ) + F T (x k=0 n ) G (k=1) T (x k=0 n ) + symmetrization. The University of Chicago p. 19
23 Fix resonances by using dt k Error on energy, symetric parareal (shift), Verlet e-04 1e-05 1e-06 1e-07 1e-08 1e-09 1e $t$ Kepler, symmetric parareal based on Verlet, dt k Error on the energy H(q, p) = pt p 2 1 q The University of Chicago p. 20
24 Going from dt k to a smaller dt with no k dependence Error on energy, symetric parareal (shift), Verlet e-04 1e-05 1e-06 1e-07 1e-08 1e-09 Error on energy, symetric parareal (without shift), Verlet e-04 1e-05 1e-06 1e-07 1e-08 1e-09 Initial coarse scheme Parareal iteration #1 Parareal iteration #4 Parareal iteration #10 At convergence 1e $t$ 1e $t$ Kepler, symmetric parareal based on Verlet Error on the energy Left: dt k Right: dt smaller Both pictures with equal cost The University of Chicago p. 21
25 Another idea On top of symmetrizing, let us project on the constant energy manifold M = { } x = (q, p) R 2d ; H(q, p) = H 0, in a way that preserves symmetry. The University of Chicago p. 22
26 Another idea On top of symmetrizing, let us project on the constant energy manifold in a way that preserves symmetry. M = { x = (q, p) R 2d ; H(q, p) = H 0 }, Take any symmetric integrator Ψ T, and consider x n x n+1 defined by x n = x n + µ H(x n ), x n+1 = Ψ T ( x n ), with µ chosen such that x n+1 M. x n+1 = x n+1 + µ H(x n+1 ), As the same Lagrange multiplier µ is used in the first and third lines, and as Ψ T is symmetric, the integrator x n x n+1 is symmetric. The University of Chicago p. 22
27 Symmetric projection x n = x n + µ H(x n ), x n+1 = Ψ T ( x n ), x n+1 = x n+1 + µ H(x n+1 ), u 3 M u 2 u 1 u 0 Need to solve a nonlinear problem to compute x n+1 and µ. In practice: use a Newton-like algorithm, and stop after a few iterations, or when the residu is small enough. In the following examples, stop after 2 iterations. The University of Chicago p. 23
28 Application to the parareal setting - 1 We need to start from a symmetric algorithm. We again introduce U n := (x 0 n, x 1 n,, x K n ), and consider the symmetric parareal algorithm previously obtained: U n+1 = Ψ sym T (U n), where the map Ψ sym T is symmetric in the classical sense. Introduce H k (U n ) := H(x k n), 1 k K and consider a symmetric algorithm, based on the symmetric map Ψ sym T and on a symmetric projection on the manifold where all the energies H k are preserved. The University of Chicago p. 24
29 Application to the parareal setting - 2 Ũ n = U n + K µ k H k (U n ) k=1 Ũ n+1 = Ψ sym T (Ũn) U n+1 = K Ũn+1 + µ k H k (U n+1 ), with the Lagrange multipliers µ k such that H k (U n+1 ) = H 0 for any 1 k K. Resulting algorithm: k=1 Initialization: at k = 0, use the symmetrized version of G: Set x 0 n+1/2 = x0 n+1/2. x 0 n+1/2 = G 1 T/2 (x0 n), x 0 n+1 = G T/2 (x 0 n+1/2 ). The University of Chicago p. 25
30 Symmetric parareal integrator with projection Assume iteration k is completed. We compute the solution at iteration k + 1 by: Set x k+1 0 = u 0 ; For 0 n N 1, compute x k+1 n+1 x k+1 n = x k+1 x k+1 n+1/2 = G 1 T/2 from xk+1 n n ) n + µ k+1 H(x k+1 as [ x k+1 n F T/2 ( x k n+1/2 ) + G T/2( x k n+1/2 ) ] x k+1 n+1 = G T/2 ( x k+1 n+1/2 ) + F T/2( x k n+1/2 ) G T/2( x k n+1/2 ) x k+1 n+1 = x k+1 n+1 + µ k+1 H(x k+1 n+1 ) where µ k+1 is such that H(x k+1 n+1 ) = H 0. As before, µ k+1 and x k+1 n+1 are determined iteratively. Jumps are precomputed in parallel, before the nonlinear iterations at step k + 1. All the expensive computations (fine propagator) are performed in parallel. The University of Chicago p. 26
31 Kepler problem (N max iter = 2): energy preservation e-04 1e-06 1e-08 1e-10 1e-12 Sequential coarse scheme 1e-14 Parareal it. #1 Parareal it. #4 Parareal it. #10 Sequential fine scheme 1e Relative error on the energy H(q, p) (δt = 10 4, dt = 0.01, T = 0.2). t Energy is well preserved at any parareal iteration k 1, although we limit the number of iterations to solve the nonlinear projection step. The University of Chicago p. 27
32 Kepler problem (N max iter = 2): angular momentum preservation e-04 Sequential coarse scheme Parareal it. #1 Parareal it. #4 Parareal it. #7 Parareal it. #10 Sequential fine scheme 1e-06 1e-08 1e-10 1e-12 1e-14 1e Relative error on the angular momentum q p t For all parareal iterations k 7, the angular momentum is preserved with a relative accuracy of 10 4 over the complete time range. The University of Chicago p. 28
33 Kepler problem (N max iter = 2): trajectory accuracy Sequential coarse scheme Parareal it. #1 Parareal it. #3 Parareal it. #5 Parareal it. #6 Sequential fine scheme e-04 1e-05 1e-06 1e-07 1e t Error on the trajectory After only 6 iterations, errors on the trajectory are comparable to the errors on the trajectory of the fine scheme used sequentially on the complete time range. The University of Chicago p. 29
34 Kepler problem: summary of the results energy is well preserved at any parareal iteration k 1 good preservation of the angular momentum (although we project on the manifold where only the energy is preserved) only 6 parareal iterations are needed to reach the best possible accuracy on the trajectory. These results are better than those obtained with the standard parareal algorithm or its symmetrized version: energy was badly preserved the standard parareal algorithm coupled with a standard, not symmetrized, projection step: better preservation of energy (remember N max iter < ), of the angular momentum, and trajectory convergence is reached for fewer parareal iterations. Both symmetrization and projection are useful. The University of Chicago p. 30
35 The outer solar system H(q, p) = 1 2 pt M 1 p + V (q), q R 18, p R 18, with V (q) = 1 i<j 6 Gm i m j q i q j, M = diag(m i). First particle Sun, with mass m m i for any planet (i 2). We show here that the symmetric parareal algorithm with symmetric projection on the constant energy manifold again behaves well. that the coarse solver can be driven by a simplified dynamics without damaging the good properties of the algorithm. The University of Chicago p. 31
36 A simplified coarse integrator Consider the simplified dynamics H(q, p) = 1 2 pt M 1 p + V simp (q), q R 18, p R 18, with V simp (q) = 2 j 6 Gm 1 m j q 1 q j. In V simp, we only take into account the gravitational interaction between the Sun (at position q 1 ) and the other planets (at position q j, 2 j 6), and we ignore the interaction between pairs of planets. m 1 m i The exact dynamics is a perturbation of the simplified dynamics. The simplified system is less expensive to simulate than the exact one, since V simp is a sum of 5 terms, whereas V is a sum of 15 terms. The University of Chicago p. 32
37 Numerical results (N max iter = 2): energy preservation 1e-04 Sequential coarse scheme Parareal it. #1 Parareal it. #4 Parareal it. #8 Parareal it. #10 Sequential fine scheme 1e-06 1e-08 1e-10 1e-12 1e t Relative error on the energy (δt = 10 2, dt = 50, T = 200). For k 8, energy is preserved up to the prescribed tolerance. The University of Chicago p. 33
38 Numerical results: trajectory accuracy Sequential coarse scheme Parareal it. #1 Parareal it. #5 Parareal it. #8 Parareal it. #10 Parareal it. #14 Sequential fine scheme e-04 1e-06 1e-08 1e t Error on the trajectory For k 5, large trajectory error. At k = 15, the trajectory error is comparable to the error made by the fine propagator used sequentially. The University of Chicago p. 34
39 Numerical results: angular momentum preservation e-04 Sequential coarse scheme Parareal it. #1 Parareal it. #5 Parareal it. #10 Parareal it. #15 Sequential fine scheme 1e-06 1e-08 1e-10 1e-12 1e-14 1e t Relative error on the angular momentum Good preservation of the angular momentum (error < 1 %), although the trajectory may be wrong. The University of Chicago p. 35
40 Numerical results: qualitative trajectory Sun Jupiter Saturn Uranus Neptune Pluto Sun Jupiter Saturn Uranus Neptune Pluto Sun Jupiter Saturn Uranus Neptune Pluto Sun Jupiter Saturn Uranus Neptune Pluto k = 1 (top left), k = 5 (top right), k = 10 (bottom left), k = 15 (bottom right). For k = 5, the trajectory is quantitatively wrong but qualitatively correct! For k 10, the trajectory is quantitatively correct. The University of Chicago p. 36
41 Computational speed-up With our parameters, the computation cost is mostly due to the fine-scale integrations. Hence, the speed-up is G = T K T = number of processors number of parareal iterations For the outer-solar system test-case, provided we have 1000 processors. G 66 The University of Chicago p. 37
42 Symplectic parareal algorithm Unclear (no standard way to make a given algorithm symplectic). A first possible approach by G. Bal and Q. Wu (2008): build a generating function S by interpolating parareal results: in practice, for some chosen φ i (q, P), S(q, P) = i a i φ i (q, P), a i to determine. derive from this generating function a symplectic algorithm: p = P + S q (q, P), Q = q + S P (q, P). How to go further? The University of Chicago p. 38
43 Conclusions the standard parareal algorithm is neither symplectic nor symmetric, even if the building blocks are: this creates issues in the long-time behaviour. easy to symmetrize, not enough due to resonances... symmetrize + symmetric projection good results! Both ingredients are needed. Some open questions: Symplectic parareal algorithm? Understand better the behaviour of the integrators through their modified equations (work in progress). X. Dai, CLB, F. Legoll, Y. Maday, Symmetric parareal algorithms for Hamiltonian systems, arxiv preprint , to appear in M2AN. Related work: M. Gander, E. Hairer, Analysis for parareal algorithms applied to Hamiltonian differential equations, preprint The University of Chicago p. 39
Chaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationPARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM
PARAREAL TIME DISCRETIZATION FOR PARABOLIC CONTROL PROBLEM Daoud S Daoud Dept. of Mathematics Eastern Mediterranean University Famagusta-North Cyprus via Mersin 10-Turkey. e mail: daoud.daoud@emu.edu.tr
More informationIntroduction to Hamiltonian systems
Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008 Examples Mathematical biology: Lotka-Volterra model First numerical methods
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationParareal Methods. Scott Field. December 2, 2009
December 2, 2009 Outline Matlab example : Advantages, disadvantages, and survey of usage in literature Overview Overview Notation Algorithm Model equation example Proposed by Lions, Maday, Turinici in
More informationParallel in Time Algorithms for Multiscale Dynamical Systems using Interpolation and Neural Networks
Parallel in Time Algorithms for Multiscale Dynamical Systems using Interpolation and Neural Networks Gopal Yalla Björn Engquist University of Texas at Austin Institute for Computational Engineering and
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationProper initial conditions for long-term integrations of the solar system
Proper initial conditions for long-term integrations of the solar system M. Arminjon Laboratoire Sols, Solides, Structures [CNRS / Université Joseph Fourier / Institut National Polytechnique de Grenoble]
More informationResonance In the Solar System
Resonance In the Solar System Steve Bache UNC Wilmington Dept. of Physics and Physical Oceanography Advisor : Dr. Russ Herman Spring 2012 Goal numerically investigate the dynamics of the asteroid belt
More informationNumerical methods in molecular dynamics and multiscale problems
Numerical methods in molecular dynamics and multiscale problems Two examples T. Lelièvre CERMICS - Ecole des Ponts ParisTech & MicMac project-team - INRIA Horizon Maths December 2012 Introduction The aim
More informationGeometric Numerical Integration
Geometric Numerical Integration (Ernst Hairer, TU München, winter 2009/10) Development of numerical ordinary differential equations Nonstiff differential equations (since about 1850), see [4, 2, 1] Adams
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 informationNumerical Algorithms as Dynamical Systems
A Study on Numerical Algorithms as Dynamical Systems Moody Chu North Carolina State University What This Study Is About? To recast many numerical algorithms as special dynamical systems, whence to derive
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationDraft TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL. Nikolai Shegunov, Ivan Hristov
TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL Nikolai Shegunov, Ivan Hristov March 9, Seminar in IMI,BAS,Soa 1/25 Mathematical model We consider the Hamiltonian :
More informationStability of the Parareal Algorithm
Stability of the Parareal Algorithm Gunnar Andreas Staff and Einar M. Rønquist Norwegian University of Science and Technology Department of Mathematical Sciences Summary. We discuss the stability of the
More informationNumerical Optimal Control Overview. Moritz Diehl
Numerical Optimal Control Overview Moritz Diehl Simplified Optimal Control Problem in ODE path constraints h(x, u) 0 initial value x0 states x(t) terminal constraint r(x(t )) 0 controls u(t) 0 t T minimize
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationThe first order formalism and the transition to the
Problem 1. Hamiltonian The first order formalism and the transition to the This problem uses the notion of Lagrange multipliers and Legendre transforms to understand the action in the Hamiltonian formalism.
More informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationPhysics 115/242 Comparison of methods for integrating the simple harmonic oscillator.
Physics 115/4 Comparison of methods for integrating the simple harmonic oscillator. Peter Young I. THE SIMPLE HARMONIC OSCILLATOR The energy (sometimes called the Hamiltonian ) of the simple harmonic oscillator
More informationOrbital Motion in Schwarzschild Geometry
Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation
More information1 Ordinary differential equations
Numerical Analysis Seminar Frühjahrssemester 08 Lecturers: Prof. M. Torrilhon, Prof. D. Kressner The Störmer-Verlet method F. Crivelli (flcrivel@student.ethz.ch May 8, 2008 Introduction During this talk
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
More informationComputational Geometric Uncertainty Propagation for Hamiltonian Systems on a Lie Group
Computational Geometric Uncertainty Propagation for Hamiltonian Systems on a Lie Group Melvin Leok Mathematics, University of California, San Diego Foundations of Dynamics Session, CDS@20 Workshop Caltech,
More informationRunga-Kutta Schemes. Exact evolution over a small time step: Expand both sides in a small time increment: d(δt) F (x(t + δt),t+ δt) Ft + FF x ( t)
Runga-Kutta Schemes Exact evolution over a small time step: x(t + t) =x(t)+ t 0 d(δt) F (x(t + δt),t+ δt) Expand both sides in a small time increment: x(t + t) =x(t)+x (t) t + 1 2 x (t)( t) 2 + 1 6 x (t)+
More informationTime-Parallel Algorithms for Weather Prediction and Climate Simulation
Time-Parallel Algorithms for Weather Prediction and Climate Simulation Jean Côté Adjunct Prof. ESCER Centre (Étude et la Simulation du Climat à l Échelle Régionale) Sciences de la Terre & de l Atmosphère
More informationNumerical Methods of Applied Mathematics -- II Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 18.336 Numerical Methods of Applied Mathematics -- II Spring 2009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationRATIO BASED PARALLEL TIME INTEGRATION (RAPTI) FOR INITIAL VALUE PROBLEMS
RATIO BASED PARALLEL TIME INTEGRATION (RAPTI) FOR INITIAL VALUE PROBLEMS Nabil R. Nassif Noha Makhoul Karam Yeran Soukiassian Mathematics Department, American University of Beirut IRISA, Campus Beaulieu,
More informationAM 205 Final Project The N-Body Problem
AM 205 Final Project The N-Body Problem Leah Birch Elizabeth Finn Karen Yu December 14, 2012 Abstract The N-Body Problem can be solved using a variety of numeric integrators. Newton s Law of Universal
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Numerics Michael Bader SCCS Technical University of Munich Summer 018 Recall: Molecular Dynamics System of ODEs resulting force acting on a molecule: F i = j
More informationJoint work with Nguyen Hoang (Univ. Concepción, Chile) Padova, Italy, May 2018
EXTENDED EULER-LAGRANGE AND HAMILTONIAN CONDITIONS IN OPTIMAL CONTROL OF SWEEPING PROCESSES WITH CONTROLLED MOVING SETS BORIS MORDUKHOVICH Wayne State University Talk given at the conference Optimization,
More informationLecture Notes of EE 714
Lecture Notes of EE 714 Lecture 1 Motivation Systems theory that we have studied so far deals with the notion of specified input and output spaces. But there are systems which do not have a clear demarcation
More informationUnder evolution for a small time δt the area A(t) = q p evolves into an area
Physics 106a, Caltech 6 November, 2018 Lecture 11: Hamiltonian Mechanics II Towards statistical mechanics Phase space volumes are conserved by Hamiltonian dynamics We can use many nearby initial conditions
More informationLEGENDRE POLYNOMIALS AND APPLICATIONS. We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.
LEGENDRE POLYNOMIALS AND APPLICATIONS We construct Legendre polynomials and apply them to solve Dirichlet problems in spherical coordinates.. Legendre equation: series solutions The Legendre equation is
More informationEstimating Accuracy in Classical Molecular Simulation
Estimating Accuracy in Classical Molecular Simulation University of Illinois Urbana-Champaign Department of Computer Science Institute for Mathematics and its Applications July 2007 Acknowledgements Ben
More informationMA/CS 615 Spring 2019 Homework #2
MA/CS 615 Spring 019 Homework # Due before class starts on Feb 1. Late homework will not be given any credit. Collaboration is OK but not encouraged. Indicate on your report whether you have collaborated
More informationHamiltonian Dynamics
Hamiltonian Dynamics CDS 140b Joris Vankerschaver jv@caltech.edu CDS Feb. 10, 2009 Joris Vankerschaver (CDS) Hamiltonian Dynamics Feb. 10, 2009 1 / 31 Outline 1. Introductory concepts; 2. Poisson brackets;
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 informationLong simulations of the Solar System: Brouwer s Law and chaos
ANZIAM J. 46 (E) pp.c1086 C1103, 2005 C1086 Long simulations of the Solar System: Brouwer s Law and chaos K. R. Grazier W. I. Newman James M. Hyman Philip. W. Sharp (received 28 October 2004, revised 18
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationHamiltonian Lecture notes Part 3
Hamiltonian Lecture notes Part 3 Alice Quillen March 1, 017 Contents 1 What is a resonance? 1 1.1 Dangers of low order approximations...................... 1. A resonance is a commensurability.......................
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationLiouville Equation. q s = H p s
Liouville Equation In this section we will build a bridge from Classical Mechanics to Statistical Physics. The bridge is Liouville equation. We start with the Hamiltonian formalism of the Classical Mechanics,
More informationDynamics and Time Integration. Computer Graphics CMU /15-662, Fall 2016
Dynamics and Time Integration Computer Graphics CMU 15-462/15-662, Fall 2016 Last Quiz: IK (and demo) Last few classes Keyframing and motion capture - Q: what can one NOT do with these techniques? The
More informationMean-Motion Resonance and Formation of Kirkwood Gaps
Yan Wang Project 1 PHYS 527 October 13, 2008 Mean-Motion Resonance and Formation of Kirkwood Gaps Introduction A histogram of the number of asteroids versus their distance from the Sun shows some distinct
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 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 informationClassical Equations of Motion
3 Classical Equations of Motion Several formulations are in use Newtonian Lagrangian Hamiltonian Advantages of non-newtonian formulations more general, no need for fictitious forces better suited for multiparticle
More informationSymplectic Runge-Kutta methods satisfying effective order conditions
/30 Symplectic Runge-Kutta methods satisfying effective order conditions Gulshad Imran John Butcher (supervisor) University of Auckland 5 July, 20 International Conference on Scientific Computation and
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More information15 Nonlinear Equations and Zero-Finders
15 Nonlinear Equations and Zero-Finders This lecture describes several methods for the solution of nonlinear equations. In particular, we will discuss the computation of zeros of nonlinear functions f(x).
More informationTHE Hamilton equations of motion constitute a system
Proceedings of the World Congress on Engineering 0 Vol I WCE 0, July 4-6, 0, London, U.K. Systematic Improvement of Splitting Methods for the Hamilton Equations Asif Mushtaq, Anne Kværnø, and Kåre Olaussen
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationThe Outer Planets (pages )
The Outer Planets (pages 720 727) Gas Giants and Pluto (page 721) Key Concept: The first four outer planets Jupiter, Saturn, Uranus, and Neptune are much larger and more massive than Earth, and they do
More informationTheory of mean motion resonances.
Theory of mean motion resonances. Mean motion resonances are ubiquitous in space. They can be found between planets and asteroids, planets and rings in gaseous disks or satellites and planetary rings.
More informationICES REPORT Parareal methods for highly oscillatory ordinary differential equations
ICES REPORT 15-6 March 215 Parareal methods for highly oscillatory ordinary differential equations by G. Ariel, S. Kim, and R. Tsai The Institute for Computational Engineering and Sciences The University
More informationMidterm EXAM PHYS 401 (Spring 2012), 03/20/12
Midterm EXAM PHYS 401 (Spring 2012), 03/20/12 Name: Signature: Duration: 75 minutes Show all your work for full/partial credit! In taking this exam you confirm to adhere to the Aggie Honor Code: An Aggie
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationScalable solution of non-linear time-dependent systems
Scalable solution of non-linear time-dependent systems Mark Maienschein-Cline and L.R. Scott April 13, 211 Abstract We study parallel solution methods for time-dependent problems where the domain decomposition
More informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationStochastic Variational Integrators
Stochastic Variational Integrators Jeremy Schmitt Variational integrators are numerical geometric integrator s derived from discretizing Hamilton s principle. They are symplectic integrators that exhibit
More informationIntroduction to the Solar System
Introduction to the Solar System Sep. 11, 2002 1) Introduction 2) Angular Momentum 3) Formation of the Solar System 4) Cowboy Astronomer Review Kepler s Laws empirical description of planetary motion Newton
More informationNumerical Methods for Embedded Optimization and Optimal Control. Exercises
Summer Course Numerical Methods for Embedded Optimization and Optimal Control Exercises Moritz Diehl, Daniel Axehill and Lars Eriksson June 2011 Introduction This collection of exercises is intended to
More information15. Hamiltonian Mechanics
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 15. Hamiltonian Mechanics Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationANALYTICAL MECHANICS. LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS
ANALYTICAL MECHANICS LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS Preface xi 1 LAGRANGIAN MECHANICS l 1.1 Example and Review of Newton's Mechanics: A Block Sliding on an Inclined Plane 1
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 informationDynamics and Time Integration
Lecture 17: Dynamics and Time Integration Computer Graphics CMU 15-462/15-662, Fall 2015 Last time: animation Added motion to our model Interpolate keyframes Still a lot of work! Today: physically-based
More informationProblem Set 5 Solution
Problem Set 5 Solution Friday, 14 October 011 Physics 111 Proble The Setup For a pair of point masses, and m, interacting via a conservative force directed along the line joining the masses, show that
More informationA645/A445: Exercise #1. The Kepler Problem
A645/A445: Exercise #1 The Kepler Problem Due: 2017 Sep 26 1 Numerical solution to the one-degree-of-freedom implicit equation The one-dimensional implicit solution is: t = t o + x x o 2(E U(x)). (1) The
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationNumerical Linear and Multilinear Algebra in Quantum Tensor Networks
Numerical Linear and Multilinear Algebra in Quantum Tensor Networks Konrad Waldherr October 20, 2013 Joint work with Thomas Huckle QCCC 2013, Prien, October 20, 2013 1 Outline Numerical (Multi-) Linear
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 17, 23) I. INTRODUCTION One frequently obtains detailed dynamical information
More informationLecture Outlines. Chapter 6. Astronomy Today 7th Edition Chaisson/McMillan Pearson Education, Inc.
Lecture Outlines Chapter 6 Astronomy Today 7th Edition Chaisson/McMillan Chapter 6 The Solar System Units of Chapter 6 6.1 An Inventory of the Solar System 6.2 Measuring the Planets 6.3 The Overall Layout
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationPhysics Unit 7: Circular Motion, Universal Gravitation, and Satellite Orbits. Planetary Motion
Physics Unit 7: Circular Motion, Universal Gravitation, and Satellite Orbits Planetary Motion Geocentric Models --Many people prior to the 1500 s viewed the! Earth and the solar system using a! geocentric
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 informationReducing round-off errors in rigid body dynamics
Reducing round-off errors in rigid body dynamics Gilles Vilmart INRIA Rennes, ENS Cachan Bretagne, avenue Robert Schuman, 357 Bruz, France Université de Genève, Section de mathématiques, -4 rue du Lièvre,
More informationPhysics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion
Physics 115/242 The leapfrog method and other symplectic algorithms for integrating Newton s laws of motion Peter Young (Dated: April 14, 2009) I. INTRODUCTION One frequently obtains detailed dynamical
More informationThe Big Picture. Discuss Examples of unpredictability. Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986)
The Big Picture Discuss Examples of unpredictability Odds, Stanisław Lem, The New Yorker (1974) Chaos, Scientific American (1986) Lecture 2: Natural Computation & Self-Organization, Physics 256A (Winter
More informationDirect Methods. Moritz Diehl. Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U.
Direct Methods Moritz Diehl Optimization in Engineering Center (OPTEC) and Electrical Engineering Department (ESAT) K.U. Leuven Belgium Overview Direct Single Shooting Direct Collocation Direct Multiple
More informationAPPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS
LECTURE 10 APPLICATIONS OF FD APPROXIMATIONS FOR SOLVING ORDINARY DIFFERENTIAL EQUATIONS Ordinary Differential Equations Initial Value Problems For Initial Value problems (IVP s), conditions are specified
More informationAsynchronous Parareal in time discretization for partial differential equations
Asynchronous Parareal in time discretization for partial differential equations Frédéric Magoulès, Guillaume Gbikpi-Benissan April 7, 2016 CentraleSupélec IRT SystemX Outline of the presentation 01 Introduction
More informationSolutions to homework assignment #3 Math 119B UC Davis, Spring = 12x. The Euler-Lagrange equation is. 2q (x) = 12x. q(x) = x 3 + x.
1. Find the stationary points of the following functionals: (a) 1 (q (x) + 1xq(x)) dx, q() =, q(1) = Solution. The functional is 1 L(q, q, x) where L(q, q, x) = q + 1xq. We have q = q, q = 1x. The Euler-Lagrange
More informationLeap Frog Solar System
Leap Frog Solar System David-Alexander Robinson Sch. 08332461 20th November 2011 Contents 1 Introduction & Theory 2 1.1 The Leap Frog Integrator......................... 2 1.2 Class.....................................
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More informationPARALLELIZATION IN TIME OF (STOCHASTIC) ORDINARY DIFFERENTIAL EQUATIONS
PARALLELIZATION IN TIME OF (STOCHASTIC) ORDINARY DIFFERENTIAL EQUATIONS GUILLAUME BAL Abstract. This paper analyzes some properties of the parareal algorithm, which can be used to parallelize the time
More informationElectrostatic Imaging via Conformal Mapping. R. Kress. joint work with I. Akduman, Istanbul and H. Haddar, Paris
Electrostatic Imaging via Conformal Mapping R. Kress Göttingen joint work with I. Akduman, Istanbul and H. Haddar, Paris Or: A new solution method for inverse boundary value problems for the Laplace equation
More informationInvariant Manifolds and Transport in the Three-Body Problem
Dynamical S C C A L T E C H Invariant Manifolds and Transport in the Three-Body Problem Shane D. Ross Control and Dynamical Systems California Institute of Technology Classical N-Body Systems and Applications
More informationSome Collision solutions of the rectilinear periodically forced Kepler problem
Advanced Nonlinear Studies 1 (2001), xxx xxx Some Collision solutions of the rectilinear periodically forced Kepler problem Lei Zhao Johann Bernoulli Institute for Mathematics and Computer Science University
More informationConservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations.
Conservation Laws: Systematic Construction, Noether s Theorem, Applications, and Symbolic Computations. Alexey Shevyakov Department of Mathematics and Statistics, University of Saskatchewan, Saskatoon,
More informationAstronomy 241: Foundations of Astrophysics I. The Solar System
Astronomy 241: Foundations of Astrophysics I. The Solar System Astronomy 241 is the first part of a year-long introduction to astrophysics. It uses basic classical mechanics and thermodynamics to analyze
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationQuest Chapter 14. Think back to circular motion. What keeps something moving around an object instead of to it?
1 Since the moon is gravitationally attracted to the Earth, why doesn t it simply crash into the Earth? 1.When the moon moves close to the Earth, the air on the Earth repels it. 2. The moon s tangential
More informationNoether s Theorem: Uses and Abuses
Noether s Theorem: Uses and Abuses Ryan Browne December 15, 2011 Contents 1 Introduction 1 2 Formulation 2 2.1 Definitions............................. 2 2.2 Formal Statement of Noether s Theorem............
More informationNon-Differentiable Embedding of Lagrangian structures
Non-Differentiable Embedding of Lagrangian structures Isabelle Greff Joint work with J. Cresson Université de Pau et des Pays de l Adour CNAM, Paris, April, 22nd 2010 Position of the problem 1. Example
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
More information