Geometric Numerical Integration

Size: px
Start display at page:

Download "Geometric Numerical Integration"

Transcription

1 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 (1855), multistep methods, problem of Bashforth (1883) Runge (1895) and Kutta (1901), one-step methods. Stiff differential equations (since about 1950), see [5, 2, 1] Dahlquist (1963), A-stability of multistep methods Gear (1971), backward differentiation code. Geometric numerical integration (since 1986), see [3, 7, 6] structure-preserving integration of differential equations, Hamiltonian, reversible, divergence-free, Poisson systems. Provisional contents of the lectures Hamiltonian systems symplectic transformations, theorem of Poincaré generating functions Symplectic numerical integrators Störmer Verlet method symplectic Runge Kutta methods composition and splitting methods Backward error analysis long-time energy conservation perturbed integrable Hamiltonian systems Hamiltonian systems on manifolds constrained mechanical systems (DAE s) numerical integrator RATTLE Differential equations with highly oscillatory solutions Fermi Pasta Ulam type problems modulated Fourier expansion 1

2 1 Hamiltonian system (double well potential) Consider a differential equation q = U(q) which can also be written as q = p ṗ = U(q) or q = p H(p, q) ṗ = q H(p, q) with H(p, q) = 1 2 pt p + U(q). We have energy conservation along exact solutions: H ( p(t), q(t) ) = Const Numerical experiment with U(q) = (q 2 1) 2, taken from vilmart/ A quartic Hamiltonian system (double well) 11/22/09 2:18 PM der 1: Symplectic Euler sleep= 30 ms speed= x 1 Start stepsize h= 0.05 q(0)= 0.01 H? p(0)= 0.0 init file:///users/hairer/cours/geo_int/javaanim/martinet.html Page 1 of 2 Explicit Euler Symplectic Euler q n+1 = q n + h p n p n+1 = p n h U(q n ) q n+1 = q n + h p n p n+1 = p n h U(q n+1 ) or q n+1 = q n + h p n+1 p n+1 = p n h U(q n ) 2

3 2 N-body problem (planetary motion) For p i, q i R 3 we consider the Hamiltonian (kinetic plus potential energy) H(p, q) = 1 2 N i=1 with the gravitational potential 1 m i p T i p i + 1 i<j N U ij ( qi q j ) U ij (r) = G m i m j r Illustration: sun jupiter saturn ( vilmart/) The three body problem Sun-Jupiter-Saturn 11/22/09 5:46 PM sleep= 40 ms speed= x 50 restart Start order 1: Explicit Euler stepsize h= 3 file:///users/hairer/cours/geo_int/javaanim/sunjupitersaturn.html Page 1 of 2 Conserved quantities (first integrals): Hamiltonian (total energy) H(p, q) linear momentum P (p, q) = N i=1 N angular momentum L(p, q) = q i p i i=1 3 p i

4 3 N-body problem (molecular dynamics) N-body problem as before, but with Lennard Jones potential ( (σij ) 12 ( σij ) ) 6 U ij (r) = 4 ε ij r r Numerical simulation with initial configuration (left picture) and solution after sufficiently long time (right picture), see vilmart/ Molecular dynamics simulation: Lennard-Jones potential Molecular dynamics simulation: Lennard-Jones potential 11/22/09 6:20 PM 11/22/09 5:47 PM sleep= 40 ms speed= x 1 restart Start sleep= 40 ms speed= x 1 restart Start order 2: Stormer-Verlet stepsize h= 0.05 order 2: Stormer-Verlet stepsize h= 0.05 file:///users/hairer/cours/geo_int/javaanim/lennardjones.html file:///users/hairer/cours/geo_int/javaanim/lennardjones.html Page 1 of 2 Page 1 of 2 4 Idea of backward error analysis For a differential equation ẏ = f(y) consider a numerical solution obtained by a one-step method y n+1 = Φ h (y n ). Find a modified differential equation ẏ = f h (y) of the form ẏ = f(y) + hf 2 (y) + h 2 f 3 (y) +..., such that its solution ỹ(t) satisfies (formally) y n = ỹ(nh). 4

5 Numerical experiment with double well potential. The following pictures show the exact solution (in phase space (q, p)) and solutions of the modified equation for various numerical integrators. exact solution explicit Euler symplectic Euler explicit in p, implicit in q symplectic Euler explicit in q, implicit in p Modified differential equation for explicit Euler: ( ) ( ) q p = ṗ U + h ( ) ( ) U (q) (q) 2 U + h2 2 U (q)p (q)p 4 2 U U U p Modified differential equation for symplectic Euler (expl. in q, impl. in p): ( ) ( ) q p = ṗ U + h ( ) ( ) U (q) (q) 2 U + h2 2 U (q)p (q)p 12 2 U U U p this modified equation is Hamiltonian with H(p, q) = 1 2 p2 + U(q) h ) 2 U (q)p + (U h2 (q) 2 + U (q)p

6 5 Constrained mechanical system Position of corners q i R 3, i = 0,..., 5 (6 3 = 18 variables) Constraints motion of one piece (red) is prescribed, i.e., q 0, q 1, q 2 are fixed (9 conditions), distances between neighbor corners is unity (4 conditions), orthogonality between neighbor edges (q n 1 q n ) (q n+1 q n ) (5 conditions). Graphics: J.-P. Eckmann & M. Hairer Dynamics (red piece is fixed, corners have equal mass one, massless edges) H(p, q) = p T i p i + i=3 5 q iz, i=3 where q iz is the vertical component of q i. We obtain a differential-algebraic equation (DAE) q i = pi H(p, q) ṗ i = qi H(p, q) G T (q)λ 0 = g(q) Here, q is the vector composed by q 3, q 4, q 5, p = q is the velocity, g(q) = 0 represents the 9 constraints for q, and G(q) = g (q) (matrix of dimension 9 9). 6

7 6 Euler s equations of motion for a free rigid body The angular momentum vector y = (y 1, y 2, y 3 ) T satisfies ẏ 1 = (I3 1 I2 1 ) y 3 y 2 ẏ 1 0 y 3 y 2 y 1 /I 1 ẏ 2 = (I1 1 I3 1 ) y 1 y 3 or ẏ 2 = y 3 0 y 1 y 2 /I 2 ẏ 3 = (I2 1 I3 1 ) y 2 y 1 ẏ 3 y 2 y 1 0 y 3 /I 3 Hamiltonian: H(y) = 1 2 ( y y2 2 + y2 3 ), Casimir: C(y) = 1 ( ) y1 2 +y2 2 +y3 2 I 1 I 2 I 3 2 1st picture: integration with explicit Euler 2nd picture: trapezoidal rule with projection onto the sphere 3rd picture: implicit midpoint rule 7

8 References [1] J. C. BUTCHER, Numerical Methods for Ordinary Differential Equations, John Wiley & Sons Ltd., Chichester, 2nd ed., [2] P. DEUFLHARD AND F. BORNEMANN, Numerische Mathematik 2, de Gruyter Lehrbuch. [de Gruyter Textbook], Walter de Gruyter & Co., Berlin, revised ed., Gewöhnliche Differentialgleichungen. [Ordinary differential equations]. [3] E. HAIRER, C. LUBICH, AND G. WANNER, Geometric Numerical Integration. Structure-Preserving Algorithms for Ordinary Differential Equations, Springer Series in Computational Mathematics 31, Springer-Verlag, Berlin, 2nd ed., [4] E. HAIRER, S. P. NØRSETT, AND G. WANNER, Solving Ordinary Differential Equations I. Nonstiff Problems, Springer Series in Computational Mathematics 8, Springer, Berlin, 2nd ed., [5] E. HAIRER AND G. WANNER, Solving Ordinary Differential Equations II. Stiff and Differential-Algebraic Problems, Springer Series in Computational Mathematics 14, Springer-Verlag, Berlin, 2nd ed., [6] B. LEIMKUHLER AND S. REICH, Simulating Hamiltonian Dynamics, Cambridge Monographs on Applied and Computational Mathematics 14, Cambridge University Press, Cambridge, [7] J. M. SANZ-SERNA AND M. P. CALVO, Numerical Hamiltonian problems, vol. 7 of Applied Mathematics and Mathematical Computation, Chapman & Hall, London,

SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS

SYMMETRIC 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 information

GEOMETRIC INTEGRATION OF ORDINARY DIFFERENTIAL EQUATIONS ON MANIFOLDS

GEOMETRIC 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 information

Backward error analysis

Backward 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 information

Introduction to Hamiltonian systems

Introduction 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 information

Symplectic integration with Runge-Kutta methods, AARMS summer school 2015

Symplectic 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 information

Reducing round-off errors in symmetric multistep methods

Reducing 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 information

1 Ordinary differential equations

1 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 information

Efficiency of Runge-Kutta Methods in Solving Simple Harmonic Oscillators

Efficiency 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 information

Reducing round-off errors in rigid body dynamics

Reducing 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 information

Lecture 4: High oscillations

Lecture 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 information

REVERSIBLE LONG-TERM INTEGRATION WITH VARIABLE STEPSIZES. 1. Introduction. We consider the numerical treatment of systems

REVERSIBLE LONG-TERM INTEGRATION WITH VARIABLE STEPSIZES. 1. Introduction. We consider the numerical treatment of systems SIAM J. SCI. COMPUT. c 1997 Society for Industrial and Applied Mathematics Vol. 18, No. 1, pp. 257 269, January 1997 014 REVERSIBLE LONG-TERM INTEGRATION WITH VARIABLE STEPSIZES ERNST HAIRER AND DANIEL

More information

Chapter 1 Symplectic Integrator and Beam Dynamics Simulations

Chapter 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 information

Algebraic structures of B-series

Algebraic structures of B-series Algebraic structures of B-series Philippe Chartier, Ernst Hairer, and Gilles Vilmart November 10, 2009 Abstract B-series are a fundamental tool in practical and theoretical aspects of numerical integrators

More information

Symmetric multistep methods for charged-particle dynamics

Symmetric 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 information

SOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS

SOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS SOME PROPERTIES OF SYMPLECTIC RUNGE-KUTTA METHODS ERNST HAIRER AND PIERRE LEONE Abstract. We prove that to every rational function R(z) satisfying R( z)r(z) = 1, there exists a symplectic Runge-Kutta method

More information

GEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS

GEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS BIT Numerical Mathematics 25, Vol. 45, No., pp. 79 723 c Springer GEOMETRIC INTEGRATION METHODS THAT PRESERVE LYAPUNOV FUNCTIONS V. GRIMM and G.R.W. QUISPEL Department of Mathematics, and Centre for Mathematics

More information

IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM

IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 8 (1), 17-30 (2002) IMPLICIT INTERVAL MULTISTEP METHODS FOR SOLVING THE INITIAL VALUE PROBLEM MAŁGORZATA JANKOWSKA 1, ANDRZEJ MARCINIAK 1,2 1 Poznań University

More information

Numerical Integration of Equations of Motion

Numerical Integration of Equations of Motion GraSMech course 2009-2010 Computer-aided analysis of rigid and flexible multibody systems Numerical Integration of Equations of Motion Prof. Olivier Verlinden (FPMs) Olivier.Verlinden@fpms.ac.be Prof.

More information

Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations

Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang

More information

Study the Numerical Methods for Solving System of Equation

Study the Numerical Methods for Solving System of Equation Study the Numerical Methods for Solving System of Equation Ravi Kumar 1, Mr. Raj Kumar Duhan 2 1 M. Tech. (M.E), 4 th Semester, UIET MDU Rohtak 2 Assistant Professor, Dept. of Mechanical Engg., UIET MDU

More information

Solving Orthogonal Matrix Differential Systems in Mathematica

Solving Orthogonal Matrix Differential Systems in Mathematica Solving Orthogonal Matrix Differential Systems in Mathematica Mark Sofroniou 1 and Giulia Spaletta 2 1 Wolfram Research, Champaign, Illinois, USA. marks@wolfram.com 2 Mathematics Department, Bologna University,

More information

Mini-course on Geometric Numerical Integration

Mini-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 information

MODIFIED DIFFERENTIAL EQUATIONS. Dedicated to Prof. Michel Crouzeix

MODIFIED DIFFERENTIAL EQUATIONS. Dedicated to Prof. Michel Crouzeix ESAIM: PROCEEDINGS, September 2007, Vol.21, 16-20 Gabriel Caloz & Monique Dauge, Editors MODIFIED DIFFERENTIAL EQUATIONS Pilippe Cartier 1, Ernst Hairer 2 and Gilles Vilmart 1,2 Dedicated to Prof. Micel

More information

ON EXPLICIT INTERVAL METHODS OF ADAMS-BASHFORTH TYPE

ON EXPLICIT INTERVAL METHODS OF ADAMS-BASHFORTH TYPE COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 8 (2), 46-57 (2002) ON EXPLICIT INTERVAL METHODS OF ADAMS-BASHFORTH TYPE MAŁGORZATA JANKOWSKA 1, ANDRZEJ MARCINIAK 1, 2 1 Poznań University of Technology,

More information

High-order symmetric schemes for the energy conservation of polynomial Hamiltonian problems.

High-order symmetric schemes for the energy conservation of polynomial Hamiltonian problems. High-order symmetric schemes for the energy conservation of polynomial Hamiltonian problems. Felice Iavernaro and Donato Trigiante Abstract We define a class of arbitrary high order symmetric one-step

More information

Variable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems

Variable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems Variable Step Runge-Kutta-Nyström Methods for the Numerical Solution of Reversible Systems J. R. Cash and S. Girdlestone Department of Mathematics, Imperial College London South Kensington London SW7 2AZ,

More information

Implicit-explicit variational integration of highly oscillatory problems

Implicit-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 information

Fourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review

Fourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review Fourth-order symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: a review G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan

More information

ONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE

ONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 5,53-65 (1999) ONE- AND TWO-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE ANDRZEJ MARCINIAK 1, BARBARA SZYSZKA 2 1 Institute of Computing Science,

More information

High-order Symmetric Schemes for the Energy Conservation of Polynomial Hamiltonian Problems 1 2

High-order Symmetric Schemes for the Energy Conservation of Polynomial Hamiltonian Problems 1 2 European Society of Computational Methods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Industrial and Applied Mathematics (JNAIAM) vol 4, no -2, 29, pp 87- ISSN 79 84 High-order

More information

THREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE

THREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE COMPUTATIONAL METHODS IN SCIENCE AND TECHNOLOGY 6, 41-59 (2000) THREE- AND FOUR-STAGE IMPLICIT INTERVAL METHODS OF RUNGE-KUTTA TYPE KAROL GAJDA 1, ANDRZEJ MARCINIAK 2, BARBARA SZYSZKA 1 1 Institute of

More information

Draft TAYLOR SERIES METHOD FOR SYSTEM OF PARTICLES INTERACTING VIA LENNARD-JONES POTENTIAL. Nikolai Shegunov, Ivan Hristov

Draft 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 information

Modified Equations for Variational Integrators

Modified 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 information

Numerical Methods for Differential Equations

Numerical Methods for Differential Equations Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods

More information

Mathematics for chemical engineers. Numerical solution of ordinary differential equations

Mathematics for chemical engineers. Numerical solution of ordinary differential equations Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler

More information

Ernst Hairer 1, Robert I. McLachlan 2 and Robert D. Skeel 3

Ernst 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 information

Scientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of

Scientific Computing with Case Studies SIAM Press, Lecture Notes for Unit V Solution of Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 1 Dianne P. O Leary c 2008 1 The

More information

A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations

A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations A New Block Method and Their Application to Numerical Solution of Ordinary Differential Equations Rei-Wei Song and Ming-Gong Lee* d09440@chu.edu.tw, mglee@chu.edu.tw * Department of Applied Mathematics/

More information

Initial value problems for ordinary differential equations

Initial value problems for ordinary differential equations AMSC/CMSC 660 Scientific Computing I Fall 2008 UNIT 5: Numerical Solution of Ordinary Differential Equations Part 1 Dianne P. O Leary c 2008 The Plan Initial value problems (ivps) for ordinary differential

More information

Butcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients

Butcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge

More information

Lecture 2: Symplectic integrators

Lecture 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 information

A NOTE ON EXPLICIT THREE-DERIVATIVE RUNGE-KUTTA METHODS (ThDRK)

A NOTE ON EXPLICIT THREE-DERIVATIVE RUNGE-KUTTA METHODS (ThDRK) BULLETIN OF THE INTERNATIONAL MATHEMATICAL VIRTUAL INSTITUTE ISSN 303-4874 (p), ISSN (o) 303-4955 www.imvibl.org / JOURNALS / BULLETIN Vol. 5(015), 65-7 Former BULLETIN OF THE SOCIETY OF MATHEMATICIANS

More information

Numerical Methods for Differential Equations

Numerical Methods for Differential Equations Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the

More information

Geometric Integrators with Application To Hamiltonian Systems

Geometric 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 information

Numerical Methods for the Solution of Differential Equations

Numerical Methods for the Solution of Differential Equations Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential

More information

the solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t

the solution of ()-(), an ecient numerical treatment requires variable steps. An alternative approach is to apply a time transformation of the form t Asymptotic Error Analysis of the Aaptive Verlet Metho Stephane Cirilli, Ernst Hairer Beneict Leimkuhler y May 3, 999 Abstract The Aaptive Verlet metho [7] an variants [6] are time-reversible schemes for

More information

Parallel-in-time integrators for Hamiltonian systems

Parallel-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 information

Chapter VII. Non-Canonical Hamiltonian Systems

Chapter 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 information

Symplectic integration. Yichao Jing

Symplectic 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 information

Energy-Preserving Runge-Kutta methods

Energy-Preserving Runge-Kutta methods Energy-Preserving Runge-Kutta methods Fasma Diele, Brigida Pace Istituto per le Applicazioni del Calcolo M. Picone, CNR, Via Amendola 122, 70126 Bari, Italy f.diele@ba.iac.cnr.it b.pace@ba.iac.cnr.it SDS2010,

More information

Discontinuous Collocation Methods for DAEs in Mechanics

Discontinuous 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 information

Review Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester

Review Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE

More information

THE Hamilton equations of motion constitute a system

THE 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 information

Lecture IV: Time Discretization

Lecture 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 information

Introduction to numerical methods for Ordinary Differential Equations

Introduction 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 information

The family of Runge Kutta methods with two intermediate evaluations is defined by

The family of Runge Kutta methods with two intermediate evaluations is defined by AM 205: lecture 13 Last time: Numerical solution of ordinary differential equations Today: Additional ODE methods, boundary value problems Thursday s lecture will be given by Thomas Fai Assignment 3 will

More information

Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations

Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao

More information

CHAPTER 10: Numerical Methods for DAEs

CHAPTER 10: Numerical Methods for DAEs CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct

More information

Symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited

Symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited Symplectic exponentially-fitted modified Runge-Kutta methods of the Gauss type: revisited G. Vanden Berghe and M. Van Daele Vakgroep Toegepaste Wiskunde en Informatica, Universiteit Gent, Krijgslaan 8

More information

Solving Ordinary Differential equations

Solving Ordinary Differential equations Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation

More information

A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS

A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 3, Fall 2009 A SYMBOLIC-NUMERIC APPROACH TO THE SOLUTION OF THE BUTCHER EQUATIONS SERGEY KHASHIN ABSTRACT. A new approach based on the use of new

More information

1. Introduction. We consider canonical Hamiltonian systems in the form 0 I y(t 0 ) = y 0 R 2m J = I 0

1. Introduction. We consider canonical Hamiltonian systems in the form 0 I y(t 0 ) = y 0 R 2m J = I 0 ENERGY AND QUADRATIC INVARIANTS PRESERVING INTEGRATORS BASED UPON GAUSS COLLOCATION FORMULAE LUIGI BRUGNANO, FELICE IAVERNARO, AND DONATO TRIGIANTE Abstract. We introduce a new family of symplectic integrators

More information

Module 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008

Module 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008 Outlines Module 4: for ODE Part I: Basic Part II: Advanced Lehrstuhl Informatik V Winter 2007/2008 Part I: Basic 1 Direction Fields 2 Euler s Method Outlines Part I: Basic Part II: Advanced 3 Discretized

More information

On the stability regions of implicit linear multistep methods

On the stability regions of implicit linear multistep methods On the stability regions of implicit linear multistep methods arxiv:1404.6934v1 [math.na] 8 Apr 014 Lajos Lóczi April 9, 014 Abstract If we apply the accepted definition to determine the stability region

More information

Application of Geometric Integration to some Mechanical Problems

Application of Geometric Integration to some Mechanical Problems Application of Geometric Integration to some Mechanical Problems Kenth Engø y and Arne Marthinsen z Department of Informatics, University of Bergen, N-5 Bergen, Norway Department of Mathematical Sciences,

More information

Estimating Accuracy in Classical Molecular Simulation

Estimating 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 information

369 Nigerian Research Journal of Engineering and Environmental Sciences 2(2) 2017 pp

369 Nigerian Research Journal of Engineering and Environmental Sciences 2(2) 2017 pp 369 Nigerian Research Journal of Engineering and Environmental Sciences 2(2) 27 pp. 369-374 Original Research Article THIRD DERIVATIVE MULTISTEP METHODS WITH OPTIMIZED REGIONS OF ABSOLUTE STABILITY FOR

More information

Exponential integrators for oscillatory second-order differential equations

Exponential 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 information

A NEW CLASS OF VARIABLE STEP-SIZE METHODS FOR MULTIBODY DYNAMICS

A NEW CLASS OF VARIABLE STEP-SIZE METHODS FOR MULTIBODY DYNAMICS MULTIBODY DYNAMICS 2005, ECCOMAS Thematic Conference J.M. Goicolea, J. Cuadrado, J.C. García Orden (eds.) Madrid, Spain, 21 24 June 2005 A NEW CLASS OF VARIABLE STEP-SIZE METHODS FOR MULTIBODY DYNAMICS

More information

DYNAMICS, NUMERICAL ANALYSIS, AND SOME GEOMETRY

DYNAMICS, NUMERICAL ANALYSIS, AND SOME GEOMETRY DYNAMICS, NUMERICAL ANALYSIS, AND SOME GEOMETRY LUDWIG GAUCKLER, ERNST HAIRER, AND CHRISTIAN LUBICH Abstract. Geometric aspects play an important role in the construction and analysis of structure-preserving

More information

Validated Explicit and Implicit Runge-Kutta Methods

Validated Explicit and Implicit Runge-Kutta Methods Validated Explicit and Implicit Runge-Kutta Methods Alexandre Chapoutot joint work with Julien Alexandre dit Sandretto and Olivier Mullier U2IS, ENSTA ParisTech 8th Small Workshop on Interval Methods,

More information

Geometric numerical integration illustrated by the Störmer Verlet method

Geometric 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 information

An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation

An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation An Implicit Runge Kutta Solver adapted to Flexible Multibody System Simulation Johannes Gerstmayr 7. WORKSHOP ÜBER DESKRIPTORSYSTEME 15. - 18. March 2005, Liborianum, Paderborn, Germany Austrian Academy

More information

Stability Ordinates of Adams Predictor-Corrector Methods

Stability Ordinates of Adams Predictor-Corrector Methods BIT manuscript No. will be inserted by the editor) Stability Ordinates of Adams Predictor-Corrector Methods Michelle L. Ghrist Bengt Fornberg Jonah A. Reeger Received: date / Accepted: date Abstract How

More information

ECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations

ECE257 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 information

ENERGY PRESERVING METHODS FOR VOLTERRA LATTICE EQUATION

ENERGY PRESERVING METHODS FOR VOLTERRA LATTICE EQUATION TWMS J. App. Eng. Math. V., N.,, pp. 9- ENERGY PRESERVING METHODS FOR VOLTERRA LATTICE EQUATION B. KARASÖZEN, Ö. ERDEM, Abstract. We investigate linear energy preserving methods for the Volterra lattice

More information

THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands

THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS. Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands THE θ-methods IN THE NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS Karel J. in t Hout, Marc N. Spijker Leiden, The Netherlands 1. Introduction This paper deals with initial value problems for delay

More information

16.7 Multistep, Multivalue, and Predictor-Corrector Methods

16.7 Multistep, Multivalue, and Predictor-Corrector Methods 16.7 Multistep, Multivalue, and Predictor-Corrector Methods 747 } free_vector(ysav,1,nv); free_vector(yerr,1,nv); free_vector(x,1,kmaxx); free_vector(err,1,kmaxx); free_matrix(dfdy,1,nv,1,nv); free_vector(dfdx,1,nv);

More information

Technische Universität Berlin

Technische Universität Berlin Technische Universität Berlin Institut für Mathematik M7 - A Skateboard(v1.) Andreas Steinbrecher Preprint 216/8 Preprint-Reihe des Instituts für Mathematik Technische Universität Berlin http://www.math.tu-berlin.de/preprints

More information

Ordinary Differential Equations II

Ordinary 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 information

Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations

Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential

More information

arxiv: v1 [math.na] 31 Oct 2016

arxiv: v1 [math.na] 31 Oct 2016 RKFD Methods - a short review Maciej Jaromin November, 206 arxiv:60.09739v [math.na] 3 Oct 206 Abstract In this paper, a recently published method [Hussain, Ismail, Senua, Solving directly special fourthorder

More information

The Geometry of Euler s equation. Introduction

The 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 information

About the numerical conservation of first integral by a class of symmetric methods

About the numerical conservation of first integral by a class of symmetric methods ÁÈ ÊÌÁÅ ÆÌÇ Á Å Ì Å ÌÁ Ö Ò Ó Ö Ó ÈÇÄÁÌ ÆÁ Ç Á ÅÁÄ ÆÇ ÓÙØ Ø ÒÙÑ Ö Ð ÓÒ ÖÚ Ø ÓÒ Ó Ö Ø ÒØ Ö Ð Ý Ð Ó ÝÑÑ ØÖ Ñ Ø Ó Ê «ÐÐ È Ú Ò ÓÐÐ Þ ÓÒ ÉÙ ÖÒ Ô ÖØ Ñ ÒØÓ ÒÙÑ ÖÓ É ÁÒ Ö ØÓ Ò Ð Ö Ú Ø Ð Ô ÖØ Ñ ÒØÓ Ò Ø ¹½½¹¾¼½¼

More information

Ordinary Differential Equations II

Ordinary 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 information

Finite Element Decompositions for Stable Time Integration of Flow Equations

Finite Element Decompositions for Stable Time Integration of Flow Equations MAX PLANCK INSTITUT August 13, 2015 Finite Element Decompositions for Stable Time Integration of Flow Equations Jan Heiland, Robert Altmann (TU Berlin) ICIAM 2015 Beijing DYNAMIK KOMPLEXER TECHNISCHER

More information

Chapter 4 Molecular Dynamics and Other Dynamics

Chapter 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 information

Applications of Mathematics

Applications of Mathematics Applications of Mathematics Zdzisław Jackiewicz; Rossana Vermiglio Order conditions for partitioned Runge-Kutta methods Applications of Mathematics Vol 45 2000) No 4 30--36 Persistent URL: http://dmlcz/dmlcz/3444

More information

CS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations

CS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80

More information

Adebayo O. Adeniran1, Saheed O. Akindeinde2 and Babatunde S. Ogundare2 * Contents. 1. Introduction

Adebayo O. Adeniran1, Saheed O. Akindeinde2 and Babatunde S. Ogundare2 * Contents. 1. Introduction Malaya Journal of Matematik Vol. No. 73-73 8 https://doi.org/.37/mjm/ An accurate five-step trigonometrically-fitted numerical scheme for approximating solutions of second order ordinary differential equations

More information

Part IB Numerical Analysis

Part IB Numerical Analysis Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after

More information

Numerical Differential Equations: IVP

Numerical Differential Equations: IVP Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential

More information

THEWAY in which time evolution acts on observables may be rephrased in

THEWAY in which time evolution acts on observables may be rephrased in Appendix F Implementing evolution F.1 Koopmania THEWAY in which time evolution acts on observables may be rephrased in the language of functional analysis, by introducing the Koopman operator, whose action

More information

An algebraic approach to invariant preserving integators: The case of quadratic and Hamiltonian invariants

An algebraic approach to invariant preserving integators: The case of quadratic and Hamiltonian invariants Numerische Mathematik manuscript No. (will be inserted by the editor) An algebraic approach to invariant preserving integators: The case of quadratic and Hamiltonian invariants Philippe Chartier 1, Erwan

More information

Structure-preserving General Linear Methods

Structure-preserving General Linear Methods Structure-preserving General Linear Methods submitted by Terence James Taylor Norton for the degree of Doctor of Philosophy of the University of Bath Department of Mathematical Sciences July 205 COPYRIGHT

More information

Numerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li

Numerical 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 information

Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations

Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep

More information

Implicit-explicit exponential integrators

Implicit-explicit exponential integrators Implicit-explicit exponential integrators Bawfeh Kingsley Kometa joint work with Elena Celledoni MaGIC 2011 Finse, March 1-4 1 Introduction Motivation 2 semi-lagrangian Runge-Kutta exponential integrators

More information

Long simulations of the Solar System: Brouwer s Law and chaos

Long 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 information

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Singly diagonally implicit Runge-Kutta methods with an explicit first stage

NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Singly diagonally implicit Runge-Kutta methods with an explicit first stage NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Singly diagonally implicit Runge-Kutta methods with an explicit first stage by Anne Kværnø PREPRINT NUMERICS NO. 1/2004 NORWEGIAN UNIVERSITY OF SCIENCE AND

More information