Introduction to Hamiltonian systems
|
|
- Laurel Holt
- 6 years ago
- Views:
Transcription
1 Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008
2 Examples Mathematical biology: Lotka-Volterra model First numerical methods Mathematical pendulum Kepler problem Outer solar system Molecular dynamics First integrals Energy, linear invariants Quadratic and polynomial invariants Reversible differential equations Symmetric methods
3 Lotka-Volterra model I u(t) number of predators v(t) number of prey u = u(v 2) v = v(1 u) general autonomous system of odes ẏ = f (y) y point in phase space f (y) vector field (velocity in y) flow: ϕ t : y 0 y(t) if y(0) = y 0
4 Lotka-Volterra model II v v 5 v u u u u number of predators, v number of prey
5 Invariant of Lotka-Volterra model equations u = u(v 2), v = v(1 u) divide by each other and separation of variables with invariant 0 = 1 u u u v 2 v v = d I(u, v) dt I(u, v) = lnu u + 2 lnv v every solution lies on level curve of I level curves are closed thus all solutions are periodic
6 First numerical methods autonomous problem y = f (y) explicit Euler method: y n+1 = y n + hf (y n ) implicit Euler method: y n+1 = y n + hf (y n+1 ) implicit midpoint rule y n+1 = y n + hf ( ) yn + y n+1 2 discrete or numerical flow: Φ h : y n y n+1
7 Partitioned systems symplectic Euler partitioned system u = f (u, v), v = g(u, v) combine explicit and implicit Euler: symplectic Euler u n+1 = u n + hf (u n, v n+1 ) v n+1 = v n + hg(u n, v n+1 ) (SE1) or u n+1 = u n + hf (u n+1, v n ) v n+1 = v n + hg(u n+1, v n ) (SE2) SE1 becomes explicit if f (u, v) = f (u), g(u, v) = g(v) SE2 becomes explicit if f (u, v) = f (v), g(u, v) = g(u)
8 Lotka-Volterra model experiment explicit Euler implicit Euler symplectic Euler v v y 0 v 6 6 y 49 y y y 0 y y 0 y u 2 4 u 2 4 u
9 Hamiltonian problem Hamiltonian H(p, q) = H(p 1,...,p d, q 1,...,q d ) (total energy) q 1,...,q d positions p 1,...,p d momenta Hamiltonian equations of motion ṗ = H q, H q = q H = q = H p ( ) H T q energy conservation: H(p(t), q(t)) = const for all t
10 Mathematical pendulum mass m = 1, massless rod of length l = 1, gravitational acceleration g = 1 Hamiltonian H(p, q) = 1 2 p2 cos q cos q q l equations of motion ṗ = H q, q = H p m ṗ = sinq, q = p or q = sinq vector field 2π-periodic in q = phase space cylinder R S 1 flow ϕ t (p, q) is an area preserving mapping
11 Area preservation
12 Pendulum numerical experiment explicit Euler h = 0.2 symplectic Euler h = 0.3 Störmer-Verlet h = 0.6
13 Kepler problem two-body problem 1st body as center of coordinate system (p, q) coordinates of second body Hamiltonian H(p 1, p 2, q 1, q 2 ) = 1 2 (p2 1 + p 2 2) (q q 2 2) 1/2 equations of motion: first integrals q i = p i, ṗ i = H qi = q i (q q 2 2) 3/2 total energy H(p, q) angular momentum L(p 1, p 2, q 1, q 2 ) = q 1 p 2 q 2 p 1 (Kepler s second law)
14 Numerical example Kepler problem steps h = implicit midpoint 1 explicit Euler steps h = 0.05 symplectic Euler steps h = Störmer Verlet steps h = 0.05
15 Numerical example Kepler problem II.02 conservation of energy explicit Euler, h = symplectic Euler, h = global error of solution explicit Euler, h = symplectic Euler, h =
16 Qualitative long-time behavior Kepler problem method error in H error in L global error explicit Euler O(th) O(th) O(t 2 h) symplectic Euler O(h) 0 O(th) implicit midpoint O(h 2 ) 0 O(th 2 ) Störmer-Verlet O(h 2 ) 0 O(th 2 )
17 Outer solar system Hamiltonian H(p, q) = i=0 1 m i p T i p i g 5 i 1 i=1 j=0 m i m j q i q j astronomical units (1 A.U. = km) masses relative to mass of sun m 0 = (account for inner planets) g = gravitational constant initial positions and initial velocity from Sept. 5, 1994, 0h00
18 Outer solar system numerical example explicit Euler, h = 10 implicit Euler, h = 10 symplectic Euler, h = 100 Störmer Verlet, h = 200
19 Molecular dynamics Hamiltonian H(p, q) = 1 2 N i=1 1 m i p T i p i + N i 1 ) V ij ( q i q j i=2 j=1 V ij (r) potential function q i, p i positions and momenta of atoms m i atomic mass of ith atom in molecular dynamics: V ij Lennard-Jones potential.2 V ij (r) = 4ε ij ( (σij r ) 12 ( σij ) ) 6 r
20 Numerical experiment frozen argon crystal N = 7 argon atoms in a plane temperature T = 1 Nk B N m i q i 2 i=1
21 Numerical experiment argon crystal explicit Euler, h = 0.5[fs] symplectic Euler, h = 10[fs] total energy explicit Euler, h = 10[fs] symplectic Euler, h = 10[fs] temperature total energy temperature Verlet, h = 40[fs] Verlet, h = 80[fs] Verlet, h = 10[fs] Verlet, h = 20[fs]
22 First integrals Definition. A non-constant function I(y) is called a first integral of ẏ = f (y) if I (y)f (y) = 0 for all y. synonyms: invariant, conserved quantity, constant of motion
23 Examples of first integrals total energy H(p, q) in Hamiltonian systems total linear and angular momentum of N-body systems H(p, q) = 1 2 N i=1 1 m i p T i p i + N i 1 V ij (r ij ), i=2 j=1 r ij = q i q j equations of motion q i = 1 m i p i, ṗ i = N j=1 ν ij (q i q j ), ν ij = V ij(r ij )/r ij linear invariants I(y) = d T y, d constant, s.t. d T f (y) = 0
24 Quadratic and polynomial invariants consider Ẏ = A(Y )Y, A(Y ) skew symmetric for all Y where Y is a vector or a matrix Theorem. The quadratic function I(Y ) = Y T Y is invariant. In particular, orthogonality of Y 0 is conserved. Lemma. Let Y, A(Y ) R n,n. If tracea(y ) = 0 for all Y, then dety is an invariant. det Y represents volume of parallelepiped generated by columns of Y volume convervation for trace A(Y ) = 0
25 Reversible differential equations Definition. Let ρ be an invertible linear transformation in the phase space of ẏ = f (y). The differential equation and the vector field f (y) are called ρ-reversible if ρf (y) = f (ρy) for all y v y f (y) v y 0 ϕ t y 1 ρf (y) f (ρy) ρ u ρ ρ u ρy ρf (y) ρy 0 ϕ t ρy 1
26 Reversible vector fields examples partitioned system where u = f (u, v), v = g(u, v) f (u, v) = f (u, v), g(u, v) = g(u, v) is (ρ)-reversible for ρ(u, v) = (u, v) second order differential equations are (ρ)-reversible ü = g(u) u = v, v = g(u) Do numerical methods produce a reversible numerical flow when applied to a reversible differential equation?
27 Symmetric methods Definition. A numerical one-step method Φ h is symmetric or time reversible if Φ h Φ h = id. y 1 = Φ h (y 0 ) is symmetric if exchanging leaves the method unaltered y 0 y 1 and h h Examples: implicit midpoint rule, Störmer-Verlet method Theorem. If a numerical method applied to a ρ-reversible differential equations satisfies ρ Φ h = Φ h ρ then Φ h is ρ-reversible if and only if Φ h is a symmetric method.
Geometric 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 informationGeometric Integrators with Application To Hamiltonian Systems
United Arab Emirates University Scholarworks@UAEU Theses Electronic Theses and Dissertations 9-2015 Geometric Integrators with Application To Hamiltonian Systems Hebatallah Jamil Al Sakaji Follow this
More informationParallel-in-time integrators for Hamiltonian systems
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
More informationIntroduction to numerical methods for Ordinary Differential Equations
Introduction to numerical methods for Ordinary Differential Equations Charles-Edouard Bréhier Abstract The aims of these lecture notes are the following. We introduce Euler numerical schemes, and prove
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 informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
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 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 informationMini-course on Geometric Numerical Integration
Mini-course on Geometric Numerical Integration Umeå University and University of Innsbruck david.cohen@umu.se http://snovit.math.umu.se/~david/ This mini-course is supported by an initiation grant from
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 informationSymplectic integration with Runge-Kutta methods, AARMS summer school 2015
Symplectic integration with Runge-Kutta methods, AARMS summer school 2015 Elena Celledoni July 13, 2015 1 Hamiltonian systems and their properties We consider a Hamiltonian system in the form ẏ = J H(y)
More informationsecond order Runge-Kutta time scheme is a good compromise for solving ODEs unstable for oscillators
ODE Examples We have seen so far that the second order Runge-Kutta time scheme is a good compromise for solving ODEs with a good precision, without making the calculations too heavy! It is unstable for
More informationPartitioned Runge-Kutta Methods for Semi-explicit Differential-Algebraic Systems of Index 2
Partitioned Runge-Kutta Methods for Semi-explicit Differential-Algebraic Systems of Index 2 A. Murua Konputazio Zientziak eta A. A. Saila, EHU/UPV, Donostia-San Sebastián email ander@si.ehu.es December
More informationEnergy-stepping integrators in
Energy-stepping integrators in Lagrangian mechanics M. Gonzalez, B. Schmidt, M. Ortiz California Institute of Technology Technische Universität München Structured Integrators Workshop Caltech, May 8, 2009
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationGEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/01/4105-0996 $16.00 2001, Vol. 41, No. 5, pp. 996 1007 c Swets & Zeitlinger GEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationChapter 1 Symplectic Integrator and Beam Dynamics Simulations
Chapter 1 and Beam Accelerator Physics Group, Journal Club November 9, 2010 and Beam NSLS-II Brookhaven National Laboratory 1.1 (70) Big Picture for Numerical Accelerator Physics and Beam For single particle
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationErnst Hairer 1, Robert I. McLachlan 2 and Robert D. Skeel 3
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher ON ENERGY CONSERVATION OF THE SIMPLIFIED TAKAHASHI IMADA METHOD, Ernst Hairer 1,
More informationu n 2 4 u n 36 u n 1, n 1.
Exercise 1 Let (u n ) be the sequence defined by Set v n = u n 1 x+ u n and f (x) = 4 x. 1. Solve the equations f (x) = 1 and f (x) =. u 0 = 0, n Z +, u n+1 = u n + 4 u n.. Prove that if u n < 1, then
More informationSYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS
BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université
More informationExponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
More informationChapter I. Examples and Numerical Experiments
Chapter I Examples and Numerical Experiments This chapter introduces some interesting examples of differential equations and illustrates different types of qualitative behaviour of numerical methods We
More informationChapter 4. Systems of ODEs. Phase Plane. Qualitative Methods
Chapter 4 Systems of ODEs. Phase Plane. Qualitative Methods Contents 4.0 Basics of Matrices and Vectors 4.1 Systems of ODEs as Models 4.2 Basic Theory of Systems of ODEs 4.3 Constant-Coefficient Systems.
More informationCS520: numerical ODEs (Ch.2)
.. CS520: numerical ODEs (Ch.2) Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca people.cs.ubc.ca/ ascher/520.html Uri Ascher (UBC) CPSC 520: ODEs (Ch. 2) Fall
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
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 informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationGeometric numerical integration illustrated by the Störmer Verlet method
Acta Numerica (003), pp. 399 450 c Cambridge University Press, 003 DOI: 0.07/S096499000044 Printed in the United Kingdom Geometric numerical integration illustrated by the Störmer Verlet method Ernst Hairer
More informationProblem Solving Session 10 Simple Harmonic Oscillator Solutions
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics 8.01 Problem Solving Session 10 Simple Harmonic Oscillator Solutions W13D3-0 Group Problem Gravitational Simple Harmonic Oscillator Two identical
More informationCOMPARISON OF GEOMETRIC INTEGRATOR METHODS FOR HAMILTON SYSTEMS
COMPARISON OF GEOMETRIC INTEGRATOR METHODS FOR HAMILTON SYSTEMS A Thesis Submitted to the Graduate School of Engineering and Sciences of İzmir Institute of Technology in Partial Fulfillment of the Requirements
More informationChapter 4 Molecular Dynamics and Other Dynamics
Chapter 4 Molecular Dynamics and Other Dynamics Molecular dynamics is a method in which the motion of each individual atom or molecule is computed according to Newton s second law. It usually involves
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 informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationNumerical methods for eigenvalue problems
Numerical methods for eigenvalue problems D. Löchel Supervisors: M. Hochbruck und M. Tokar Mathematisches Institut Heinrich-Heine-Universität Düsseldorf GRK 1203 seminar february 2008 Outline Introduction
More informationA little taste of symplectic geometry
A little taste of symplectic geometry Mathematics Seminar The Richard Stockton College of New Jersey Pomona, New Jersey Timothy E. Goldberg Cornell University Ithaca, New York goldberg@math.cornell.edu
More informationSymmetric multistep methods for charged-particle dynamics
Version of 12 September 2017 Symmetric multistep methods for charged-particle dynamics Ernst Hairer 1, Christian Lubich 2 Abstract A class of explicit symmetric multistep methods is proposed for integrating
More informationLecture 1: Oscillatory motions in the restricted three body problem
Lecture 1: Oscillatory motions in the restricted three body problem Marcel Guardia Universitat Politècnica de Catalunya February 6, 2017 M. Guardia (UPC) Lecture 1 February 6, 2017 1 / 31 Outline of the
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationSpringer Series in Computational Mathematics
Springer Series in Computational Mathematics 31 Editorial Board R. Bank R.L. Graham J. Stoer R. Varga H. Yserentant Ernst Hairer Christian Lubich Gerhard Wanner Geometric Numerical Integration Structure-Preserving
More informationEfficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators
MATEMATIKA, 8, Volume 3, Number, c Penerbit UTM Press. All rights reserved Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators Annie Gorgey and Nor Azian Aini Mat Department of Mathematics,
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 32 Theme of Last Few
More information1 2 predators competing for 1 prey
1 2 predators competing for 1 prey I consider here the equations for two predator species competing for 1 prey species The equations of the system are H (t) = rh(1 H K ) a1hp1 1+a a 2HP 2 1T 1H 1 + a 2
More informationOrdinary Differential Equations I
Ordinary Differential Equations I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations I 1 / 27 Theme of Last Three
More informationLecture 4: High oscillations
Geometric Numerical Integration TU München Ernst Hairer January February Lecture 4: High oscillations Table of contents A Fermi Pasta Ulam type problem Application of classical integrators 4 3 Exponential
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 informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationMath 128A Spring 2003 Week 12 Solutions
Math 128A Spring 2003 Week 12 Solutions Burden & Faires 5.9: 1b, 2b, 3, 5, 6, 7 Burden & Faires 5.10: 4, 5, 8 Burden & Faires 5.11: 1c, 2, 5, 6, 8 Burden & Faires 5.9. Higher-Order Equations and Systems
More informationMath 692A: Geometric Numerical Integration
Math 692A: Geometric Numerical Integration Melvin Leok Mathematics, Purdue University. mleok@math.purdue.edu http://www.math.purdue.edu/ mleok/ Mathematics, Purdue University NSF DMS-54747, 714223 and
More informationLecture 2: Symplectic integrators
Geometric Numerical Integration TU Müncen Ernst Hairer January February 010 Lecture : Symplectic integrators Table of contents 1 Basic symplectic integration scemes 1 Symplectic Runge Kutta metods 4 3
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 informationSuperintegrability and exactly solvable problems in classical and quantum mechanics
Superintegrability and exactly solvable problems in classical and quantum mechanics Willard Miller Jr. University of Minnesota W. Miller (University of Minnesota) Superintegrability Penn State Talk 1 /
More informationReduced models and numerical analysis in molecular quantum dynamics I. Variational approximation
Reduced models and numerical analysis in molecular quantum dynamics I. Variational approximation Christian Lubich Univ. Tübingen Bressanone/Brixen, 14 February 2011 Computational challenges and buzzwords
More informationEULER-LAGRANGE TO HAMILTON. The goal of these notes is to give one way of getting from the Euler-Lagrange equations to Hamilton s equations.
EULER-LAGRANGE TO HAMILTON LANCE D. DRAGER The goal of these notes is to give one way of getting from the Euler-Lagrange equations to Hamilton s equations. 1. Euler-Lagrange to Hamilton We will often write
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 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 informationSymplectic integration. Yichao Jing
Yichao Jing Hamiltonian & symplecticness Outline Numerical integrator and symplectic integration Application to accelerator beam dynamics Accuracy and integration order Hamiltonian dynamics In accelerator,
More informationDiscontinuous Collocation Methods for DAEs in Mechanics
Discontinuous Collocation Methods for DAEs in Mechanics Scott Small Laurent O. Jay The University of Iowa Iowa City, Iowa, USA July 11, 2011 Outline 1 Introduction of the DAEs 2 SPARK and EMPRK Methods
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 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 informationBIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs
BIFURCATION PHENOMENA Lecture 1: Qualitative theory of planar ODEs Yuri A. Kuznetsov August, 2010 Contents 1. Solutions and orbits. 2. Equilibria. 3. Periodic orbits and limit cycles. 4. Homoclinic orbits.
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
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 informationChapter VII. Non-Canonical Hamiltonian Systems
Chapter VII. Non-Canonical Hamiltonian Systems We discuss theoretical properties and the structure-preserving numerical treatment of Hamiltonian systems on manifolds and of the closely related class of
More informationANALYTISK MEKANIK I HT 2016
Karlstads Universitet Fysik ANALYTISK MEKANIK I HT 2016 Kursens kod: FYGB08 Kursansvarig lärare: Jürgen Fuchs rum 21F 316 tel. 054-700 1817 el.mail: juerfuch@kau.se FYGB08 HT 2016 Exercises 1 2016-12-14
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 informationAPPLICATION OF CANONICAL TRANSFORMATION TO GENERATE INVARIANTS OF NON-CONSERVATIVE SYSTEM
Indian J pure appl Math, 39(4: 353-368, August 2008 c Printed in India APPLICATION OF CANONICAL TRANSFORMATION TO GENERATE INVARIANTS OF NON-CONSERVATIVE SYSTEM ASHWINI SAKALKAR AND SARITA THAKAR 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 informationTheory of Ordinary Differential Equations. Stability and Bifurcation I. John A. Burns
Theory of Ordinary Differential Equations Stability and Bifurcation I John A. Burns Center for Optimal Design And Control Interdisciplinary Center for Applied Mathematics Virginia Polytechnic Institute
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationCopyright is owned by the Author of the thesis. Permission is given for a copy to be downloaded by an individual for the purpose of research and
Copyright is owned by the Author of the thesis. Permission is given for a copy to be downloaded by an individual for the purpose of research and private study only. The thesis may not be reproduced elsewhere
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 informationOn linear and non-linear equations.(sect. 2.4).
On linear and non-linear equations.sect. 2.4). Review: Linear differential equations. Non-linear differential equations. Properties of solutions to non-linear ODE. The Bernoulli equation. Review: Linear
More informationPOPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL. If they co-exist in the same environment:
POPULATION DYNAMICS: TWO SPECIES MODELS; Susceptible Infected Recovered (SIR) MODEL Next logical step: consider dynamics of more than one species. We start with models of 2 interacting species. We consider,
More information= 0. = q i., q i = E
Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More information18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More informationA Brief Outline of Math 355
A Brief Outline of Math 355 Lecture 1 The geometry of linear equations; elimination with matrices A system of m linear equations with n unknowns can be thought of geometrically as m hyperplanes intersecting
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
More informationCPS 5310: Parameter Estimation. Natasha Sharma, Ph.D.
Example Suppose our task is to determine the net income for year 2019 based on the net incomes given below Year Net Income 2016 48.3 million 2017 90.4 million 2018 249.9 million Last lecture we tried to
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 informationCopyright is owned by the Author of the thesis. Permission is given for a copy to be downloaded by an individual for the purpose of research and
Copyright is owned by the Author of the thesis. Permission is given for a copy to be downloaded by an individual for the purpose of research and private study only. The thesis may not be reproduced elsewhere
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 informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationEXTENSIONS OF THE HHT-α METHOD TO DIFFERENTIAL-ALGEBRAIC EQUATIONS IN MECHANICS
EXTENSIONS OF THE HHT-α METHOD TO DIFFERENTIAL-ALGEBRAIC EQUATIONS IN MECHANICS LAURENT O. JAY AND DAN NEGRUT Abstract. We present second-order extensions of the Hilber-Hughes-Taylor-α HHT-α) method for
More informationChaotic 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 informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
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 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 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 informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More informationFirst Order Systems of Linear Equations. or ODEs of Arbitrary Order
First Order Systems of Linear Equations or ODEs of Arbitrary Order Systems of Equations Relate Quantities Examples Predator-Prey Relationships r 0 = r (100 f) f 0 = f (r 50) (Lokta-Volterra Model) Systems
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationVariational construction of periodic and connecting orbits in the planar Sitnikov problem. Mitsuru Shibayama(Kyoto University)
Variational construction of periodic and connecting orbits in the planar Sitnikov problem Mitsuru Shibayama(Kyoto University) 1 The three-body problem Consider the planar three-body problem which is governed
More informationAVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS
AVERAGING AND RECONSTRUCTION IN HAMILTONIAN SYSTEMS Kenneth R. Meyer 1 Jesús F. Palacián 2 Patricia Yanguas 2 1 Department of Mathematical Sciences University of Cincinnati, Cincinnati, Ohio (USA) 2 Departamento
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 informationLecture Tutorial: Angular Momentum and Kepler s Second Law
2017 Eclipse: Research-Based Teaching Resources Lecture Tutorial: Angular Momentum and Kepler s Second Law Description: This guided inquiry paper-and-pencil activity helps students to describe angular
More informationLectures on Dynamical Systems. Anatoly Neishtadt
Lectures on Dynamical Systems Anatoly Neishtadt Lectures for Mathematics Access Grid Instruction and Collaboration (MAGIC) consortium, Loughborough University, 2007 Part 3 LECTURE 14 NORMAL FORMS Resonances
More informationVariational Integrators for Electrical Circuits
Variational Integrators for Electrical Circuits Sina Ober-Blöbaum California Institute of Technology Joint work with Jerrold E. Marsden, Houman Owhadi, Molei Tao, and Mulin Cheng Structured Integrators
More information