Symplectic Runge-Kutta methods satisfying effective order conditions
|
|
- Blaise Briggs
- 5 years ago
- Views:
Transcription
1 /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 Differential equations (SciCADE 20) University of Toronto, Canada
2 2/30 Introduction The idea of effective order was first introduced by Butcher in 969 for explicit Runge Kutta methods as a mean of overcoming the 5 th order 5 stage barrier. 2 This was extended to Singly Implicit Runge Kutta methods by Butcher and Chartier in The idea was later used for Diagonally Extended Singly Implicit Runge Kutta methods by Butcher and Diamantakis in 998 and by Butcher and Chen in The accuracy of symplectic integrators for Hamiltonian systems was enhanced using effective order by M A Lopéz-Marcos, J M Sanz-Serna and R D Skeel in 996.
3 Introduction The idea of effective order was first introduced by Butcher in 969 for explicit Runge Kutta methods as a mean of overcoming the 5 th order 5 stage barrier. 2 This was extended to Singly Implicit Runge Kutta methods by Butcher and Chartier in The idea was later used for Diagonally Extended Singly Implicit Runge Kutta methods by Butcher and Diamantakis in 998 and by Butcher and Chen in The accuracy of symplectic integrators for Hamiltonian systems was enhanced using effective order by M A Lopéz-Marcos, J M Sanz-Serna and R D Skeel in 996. The purpose of this presentation is to develop some special symplectic effective order methods with low implementation costs. 2/30
4 3/30 Differential Equations with Invariants Consider an initial value problem y (x) = f(y(x)), y(x 0 ) = y 0. () Suppose Q(y) = y T My is a quadratic invariant, that is or Q (y)f(y) = 0, y T Mf(y) = 0. Methods which conserve quadratic invariants are said to be Canonical or Symplectic.
5 4/30 Outline Effective order Algebraic interpretation Computational interpretation
6 4/30 Outline Effective order Algebraic interpretation Computational interpretation 2 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Order conditions for Symplectic methods
7 4/30 Outline Effective order Algebraic interpretation Computational interpretation 2 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Order conditions for Symplectic methods 3 Effective order with symplectic integrator New Implicit Methods Cheap implementation Transformation Starting method
8 4/30 Outline Effective order Algebraic interpretation Computational interpretation 2 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Order conditions for Symplectic methods 3 Effective order with symplectic integrator New Implicit Methods Cheap implementation Transformation Starting method 4 Numerical Experiments
9 4/30 Outline Effective order Algebraic interpretation Computational interpretation 2 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Order conditions for Symplectic methods 3 Effective order with symplectic integrator New Implicit Methods Cheap implementation Transformation Starting method 4 Numerical Experiments 5 Conclusions
10 Outline Effective order Algebraic interpretation Computational interpretation 2 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Order conditions for Symplectic methods 3 Effective order with symplectic integrator New Implicit Methods Cheap implementation Transformation Starting method 4 Numerical Experiments 5 Conclusions 6 Future work 4/30
11 5/30 Effective order Algebraic interpretation Algebraic interpretation We introduce an algebraic system which represents indvidual Runge Kutta methods and also composition of methods. 2 This centres on the meaning of order for Runge Kutta methods and leads to the possible generalisation to the effective order. 3 We introduce a group G whose elements are mapping from T (rooted trees) to real numbers and where the group operation is defined according to the algebraic theory of Runge Kutta methods or to the theory of B series. 4 Members of G represents Runge Kutta methods with E representing the exact solution. That is, E : T R is defined by E(t) = γ(t), t
12 Effective order Algebraic interpretation Table: Group Operation t r(t) α(t) β(t) (αβ)(t) E(t) α β α β 0 +β 2 α 2 β 2 α 2 β 0 +β 2 +α β 2 3 α 3 β 3 α 3 β 0 +β 3 +α 2 β +2α β α 4 β 4 α 4 β 0 +β 4 +α β 2 +α 2 β 4 α 5 β 5 α 5 β 0 +β 5 +3α β 3 +α 2 β +α 3 β 4 α 6 β 6 α 6 β 0 +β 6 +α β 4 +α β α 2 β 2 +α 2 β 2 +α α 2 β α 7 β 7 α 7 β 0 +β 7 +2α β 4 +α 3 β +α 2 β α 8 β 8 α 8 β 0 +β 8 +α 2 β 2 +α β 4 +α 4 β 24 6/30
13 7/30 We introduce N p as a normal subgroup, which is defined by N p = {α G : α(t) = 0,whenever r(t) p} A Runge- Kutta method with group element α is of order p, if it is in the same coset as EN p, that is αn p = EN p A Runge- Kutta method has an effective order p if there exist another Runge - Kutta method with corresponding group element β, such that βαn p = EβN p
14 8/30 Computational interpretation Computational interpretation The conjugacy concept in group theory provides a computational interpretation of the effective order. This means that, every input value for effective order method α is perturbed by a method β. Therefore the starting method β offers some freedom of the order conditions of the effective order method α. Every output value could also be perturbed back to the origional trajectory using method β. α n β β x 0 E n x n
15 9/30 Symplectic Runge Kutta methods Symplectic Runge Kutta methods A Runge Kutta method is said to be canonical or symplectic if the numerical solution y n also has the quadratic invariant Q(y n ) i.e. y n,y n = y n,y n A method has this property if and only if, for all i,j. b i a ij +b j a ji b i b j = 0
16 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Superfluous and Non superfluous trees It is a consequence of the symplectic condition that if τ and τ 2 are rooted trees corressponding to the same tree τ then φ(τ ) = γ(τ ), φ(τ 2) = γ(τ 2 ) For Symplectic Runge Kutta methods, we distinguish trees in two ways. Superfluous trees, Non Superfluous trees. 2 τ τ τ 2 0/30
17 /30 Symplectic Runge Kutta methods Superfluous and Non superfluous trees Order conditions corresponding to non superfluous trees are transformed into one order condition. a τ b τ a τ b
18 2/30 Symplectic Runge Kutta methods Order conditions for Symplectic methods Order Conditions for Symplectic Runge Kutta methods The number of order conditions for symplectic Runge Kutta methods is less than the number of order conditions for a general Runge Kutta method. Case : Suppose the method is of order at least, ( b i = ), i b i a ij + i,j i,j b j a ji i,j b i b j = 0 i,j b i a ij = 2 Therefore for symplectic Runge Kutta method the second order condition is automatically satisfied and hence not required.
19 Symplectic Runge Kutta methods Order conditions for Symplectic methods Order Conditions for Symplectic Runge Kutta methods Case 2 : Suppose the method is of order at least 2, ( i b i c i = ), Consider the symplectic condition, 2 b i a ij +b j a ji b i b j = 0 Multiply with c j and take summation, b i a ij c j + b j c j a ji b i b j c j = 0 i,j i,j i,j ( b i a ij c j 6 )+( b j cj 2 3 ) = 0 i,j j
20 3/30 Symplectic Runge Kutta methods Order conditions for Symplectic methods Order Conditions for Symplectic Runge Kutta methods Case 2 : Suppose the method is of order at least 2, ( i b i c i = ), Consider the symplectic condition, 2 b i a ij +b j a ji b i b j = 0 Multiply with c j and take summation, b i a ij c j + b j c j a ji b i b j c j = 0 i,j i,j i,j ( b i a ij c j 6 )+( b j cj 2 3 ) = 0 i,j j
21 4/30 Symplectic Runge Kutta methods Order conditions for Symplectic methods Order General RK method Symplectic RK method Table: Order conditions for general and symplectic Runge Kutta methods up to order 4.
22 Symplectic Runge Kutta methods Order conditions for Symplectic methods Gauss method For example, consider applying Gauss method to the harmonic oscillator problem given below: The energy is given by, q = p, p = q. H = p2 2 + q2 2. The exact solution is, [ ] [ ][ ] p(x) cos(x) sin(x) p(0) =. q(x) sin(x) cos(x) q(0) This problem describes the motion of a unit mass attached to a spring with momentum p and position co-ordinates q defines a Hamiltonian system. 5/30
23 6/30 Symplectic Runge Kutta methods Order conditions for Symplectic methods we consider the two stage order four Gauss method, (2) 2 2 Gauss method approximately conserve the total energy of the himiltonian system.
24 7/30 Effective order with symplectic integrator Effective order with symplectic integrator We show a way of analyzing methods of effective order 4 having symplectic condition with three stages. (βα)( ) = β( )+α( ) (βα)( ) = β( )+2β( )α( )+β 2 ( )α( )+α( ) (βα)( ) = α( )+3β( )α( )+3β 2 ( )α( ) +β 3 ( )α( )+β( )
25 Effective order with symplectic integrator New Implicit Methods New Implicit methods We present two examples of symplectic effective order methods: Method A - has real and distinct eigenvalues Method B - has complex eigenvalues /30
26 9/30 Effective order with symplectic integrator Cheap implementation Cheap implementation Since method A have real eigenvalues, it is therefore of interest.this is because we can obtain a cheaper implementation. Here we are considering only method A. The general form of an s-stage implicit method is y n+ = y n +h Y i = y n +h s b i f(x n +hc i,y i ), i= s a ij f(x n +hc j,y j ) j= The stage equations can be written in the form Y = e y n +h(a I m )F(Y)
27 20/30 Effective order with symplectic integrator Cheap implementation We use modified Newton Raphson iteration scheme to solve the above equation. This can be defined as where M Y [k] = G(Y [k] ) Y [k+] = Y [k] + Y [k] M = I s I m h(a J) G(Y [k] ) = Y [k] +e y n +h(a I m )F(Y) The total computational cost in this scheme include the evaluation of F and G, the evaluation of J, the evaluation of M, LU factorization of the iteration matrix, M, back substitution to get the Newton update vector, Y [k].
28 Effective order with symplectic integrator Transformation Transformation To reduce computational cost of fully implicit RK method, we use transformation. The transformation matrix T for method A is given by T = The transformation matrix has the property, T AT = where , , and are three distinct real eigenvalues 2/30
29 22/30 Effective order with symplectic integrator Starting method Starting method Solution of these equations give the starting method (α)( ) = (α)( ) = 2β( )+ 3 (α)( ) = 3β( )+3β( )+ 4 which is given by
30 23/30 Numerical Experiments Numerical Experiments The Kepler s problem x = y, x 2 = y 2, x y =, y (x 2 2 = +x2 2 )3 2 (x 2 +x2 2 )3 2 x 2 where (x,x 2 ) are the position coordinates and (y,y 2 ) are the velocity components of the body. (x,x 2,y,y 2 ) = ( e,0,0, (+e)/( e)), H = 2 (y2 +y2) 2. x 2 +x2 2
31 24/30 Numerical Experiments Kepler problem (e=0,h=0.0, n = 0 6 ) Graph for Hamiltonian: error VS time 20 x
32 25/30 Numerical Experiments Kepler problem (e=0.5,h=0.0, n = 0 6 ) For Hamiltonian: error VS time 20 x
33 26/30 Numerical Experiments The simple Pendulum p = sin(q), q = p, (p,q) = (0,2.3). H = p2 2 cos(q).
34 27/30 Numerical Experiments Simple Pendulum (h=0.05, n = 0 6 ) For Hamiltonian: error VS time 6 x x 0 4
35 28/30 Conclusions Conclusions For problems that conserve some sort of invariant structure, it is a good idea to use numerical methods which mimic this behaviour. 2 Symplectic Runge Kutta methods have this role for many important problems. 3 Because of greater flexibilty,effective order methods can provide greater efficiency as compared with methods with classical order. 4 It is possible to obtain cheap implementation cost if A has real eigenvalues. 5 These methods are suited for parallel computers which have very large number of processors.
36 29/30 Future work Future work Error estimates 2 Working on implicit methods with optimal choices of parameters. 3 Construct general linear methods with closely related properties. 4 Generalization of effective order on partioned Runge Kutta methods for separable Hamiltonian.
37 30/30 Future work THANK YOU
Parallel Methods for ODEs
Parallel Methods for ODEs Levels of parallelism There are a number of levels of parallelism that are possible within a program to numerically solve ODEs. An obvious place to start is with manual code restructuring
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 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 informationRunge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot. Department U2IS ENSTA ParisTech SCAN Uppsala
Runge-Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto Alexandre Chapoutot Department U2IS ENSTA ParisTech SCAN 2016 - Uppsala Contents Numerical integration Runge-Kutta with interval
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
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 informationInvestigation on the Most Efficient Ways to Solve the Implicit Equations for Gauss Methods in the Constant Stepsize Setting
Applied Mathematical Sciences, Vol. 12, 2018, no. 2, 93-103 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2018.711340 Investigation on the Most Efficient Ways to Solve the Implicit Equations
More informationSOME 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 informationConjugate Gradient (CG) Method
Conjugate Gradient (CG) Method by K. Ozawa 1 Introduction In the series of this lecture, I will introduce the conjugate gradient method, which solves efficiently large scale sparse linear simultaneous
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 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 informationNumerical solution of stiff ODEs using second derivative general linear methods
Numerical solution of stiff ODEs using second derivative general linear methods A. Abdi and G. Hojjati University of Tabriz, Tabriz, Iran SciCADE 2011, Toronto 0 5 10 15 20 25 30 35 40 1 Contents 1 Introduction
More informationEnergy-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 informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
Lab 1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationExponentially Fitted Symplectic Runge-Kutta-Nyström methods
Appl. Math. Inf. Sci. 7, No. 1, 81-85 (2013) 81 Applied Mathematics & Information Sciences An International Journal Exponentially Fitted Symplectic Runge-Kutta-Nyström methods Th. Monovasilis 1, Z. Kalogiratou
More informationReview 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 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 informationRUNGE-KUTTA SCHEMES FOR HAMILTONIAN SYSTEMS
BIT28 (1988),877-8.83 RUNGE-KUTTA SCHEMES FOR HAMILTONIAN SYSTEMS J. M. SANZ-SERNA Departamento de Materndtica Apiicada y Computaci6n, Facuttad de Ciencias, Universidad de Valladolid, Valladolid, Spain
More informationAN OVERVIEW. Numerical Methods for ODE Initial Value Problems. 1. One-step methods (Taylor series, Runge-Kutta)
AN OVERVIEW Numerical Methods for ODE Initial Value Problems 1. One-step methods (Taylor series, Runge-Kutta) 2. Multistep methods (Predictor-Corrector, Adams methods) Both of these types of methods are
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 information(f(x) P 3 (x)) dx. (a) The Lagrange formula for the error is given by
1. QUESTION (a) Given a nth degree Taylor polynomial P n (x) of a function f(x), expanded about x = x 0, write down the Lagrange formula for the truncation error, carefully defining all its elements. How
More informationEXAMPLE OF ONE-STEP METHOD
EXAMPLE OF ONE-STEP METHOD Consider solving y = y cos x, y(0) = 1 Imagine writing a Taylor series for the solution Y (x), say initially about x = 0. Then Y (h) = Y (0) + hy (0) + h2 2 Y (0) + h3 6 Y (0)
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationNumerical Methods for Initial Value Problems; Harmonic Oscillators
1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.
More informationP-stable Exponentially fitted Obrechkoff methods
P-stable Exponentially fitted Obrechkoff methods Marnix Van Daele, G. Vanden Berghe Marnix.VanDaele@UGent.be Vakgroep Toegepaste Wiskunde en Informatica Universiteit Gent SciCADE 7 p. /4 Outline Introduction
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 informationOn the reliability and stability of direct explicit Runge-Kutta integrators
Global Journal of Pure and Applied Mathematics. ISSN 973-78 Volume, Number 4 (), pp. 3959-3975 Research India Publications http://www.ripublication.com/gjpam.htm On the reliability and stability of direct
More informationORBIT 14 The propagator. A. Milani et al. Dipmat, Pisa, Italy
ORBIT 14 The propagator A. Milani et al. Dipmat, Pisa, Italy Orbit 14 Propagator manual 1 Contents 1 Propagator - Interpolator 2 1.1 Propagation and interpolation............................ 2 1.1.1 Introduction..................................
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationImproved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3
Improved Starting Methods for Two-Step Runge Kutta Methods of Stage-Order p 3 J.H. Verner August 3, 2004 Abstract. In [5], Jackiewicz and Verner derived formulas for, and tested the implementation of two-step
More informationCHAPTER 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 informationUniversity of Utrecht
University of Utrecht Bachelor thesis The bifurcation diagram of the second nontrivial normal form of an axially symmetric perturbation of the isotropic harmonic oscillator Author: T. Welker Supervisor:
More informationAn 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 informationReview. Numerical Methods Lecture 22. Prof. Jinbo Bi CSE, UConn
Review Taylor Series and Error Analysis Roots of Equations Linear Algebraic Equations Optimization Numerical Differentiation and Integration Ordinary Differential Equations Partial Differential Equations
More informationCoupled differential equations
Coupled differential equations Example: dy 1 dy 11 1 1 1 1 1 a y a y b x a y a y b x Consider the case with b1 b d y1 a11 a1 y1 y a1 ay dy y y e y 4/1/18 One way to address this sort of problem, is to
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 informationSolution of Matrix Eigenvalue Problem
Outlines October 12, 2004 Outlines Part I: Review of Previous Lecture Part II: Review of Previous Lecture Outlines Part I: Review of Previous Lecture Part II: Standard Matrix Eigenvalue Problem Other Forms
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 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 informationNumerical 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 informationMB4018 Differential equations
MB4018 Differential equations Part II http://www.staff.ul.ie/natalia/mb4018.html Prof. Natalia Kopteva Spring 2015 MB4018 (Spring 2015) Differential equations Part II 0 / 69 Section 1 Second-Order Linear
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationValidated 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 informationDepartment of Applied Mathematics and Theoretical Physics. AMA 204 Numerical analysis. Exam Winter 2004
Department of Applied Mathematics and Theoretical Physics AMA 204 Numerical analysis Exam Winter 2004 The best six answers will be credited All questions carry equal marks Answer all parts of each question
More informationRunge Kutta Theory and Constraint Programming
Runge Kutta Theory and Constraint Programming Julien Alexandre dit Sandretto julien.alexandre-dit-sandretto@ensta-paristech.fr U2IS, ENSTA ParisTech, Université Paris-Saclay, 828 bd des Maréchaux, 9762
More informationPart 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 informationNumerical Algorithms 0 (1998)?{? 1. Position Versus Momentum Projections for. Constrained Hamiltonian Systems. Werner M. Seiler
Numerical Algorithms 0 (1998)?{? 1 Position Versus Momentum Projections for Constrained Hamiltonian Systems Werner M. Seiler Lehrstuhl I fur Mathematik, Universitat Mannheim, 68131 Mannheim, Germany E-mail:
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 informationBefore you begin read these instructions carefully.
MATHEMATICAL TRIPOS Part IB Tuesday, 5 June, 2012 9:00 am to 12:00 pm PAPER 1 Before you begin read these instructions carefully. Each question in Section II carries twice the number of marks of each question
More informationLOWELL WEEKLY JOURNAL
Y G y G Y 87 y Y 8 Y - $ X ; ; y y q 8 y $8 $ $ $ G 8 q < 8 6 4 y 8 7 4 8 8 < < y 6 $ q - - y G y G - Y y y 8 y y y Y Y 7-7- G - y y y ) y - y y y y - - y - y 87 7-7- G G < G y G y y 6 X y G y y y 87 G
More informationSolving 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 informationAdvanced methods for ODEs and DAEs
Advanced methods for ODEs and DAEs Lecture 13: Overview Bojana Rosic, 13. Juli 2016 Dynamical system is a system that evolves over time possibly under external excitations: river flow, car, etc. The dynamics
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationFourth-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 informationNumerical 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 informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationStructure-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 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 informationPrevious Year Questions & Detailed Solutions
Previous Year Questions & Detailed Solutions. The rate of convergence in the Gauss-Seidal method is as fast as in Gauss Jacobi smethod ) thrice ) half-times ) twice 4) three by two times. In application
More informationMath Matrix Algebra
Math 44 - Matrix Algebra Review notes - (Alberto Bressan, Spring 7) sec: Orthogonal diagonalization of symmetric matrices When we seek to diagonalize a general n n matrix A, two difficulties may arise:
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
More informationA NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN FILM FLOW PROBLEM
Indian J. Pure Appl. Math., 491): 151-167, March 218 c Indian National Science Academy DOI: 1.17/s13226-18-259-6 A NEW INTEGRATOR FOR SPECIAL THIRD ORDER DIFFERENTIAL EQUATIONS WITH APPLICATION TO THIN
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationOrdinary Differential Equations (ODEs)
Ordinary Differential Equations (ODEs) 1 Computer Simulations Why is computation becoming so important in physics? One reason is that most of our analytical tools such as differential calculus are best
More informationON THE IMPLEMENTATION OF SINGLY IMPLICIT RUNGE-KUTTA METHODS
MATHEMATICS OF COMPUTATION VOLUME 57, NUMBER 196 OCTOBER 1991, PAGES 663-672 ON THE IMPLEMENTATION OF SINGLY IMPLICIT RUNGE-KUTTA METHODS G. J. COOPER Abstract. A modified Newton method is often used to
More informationB-Series and their Cousins - A. Contemporary Geometric View
B-Series and their Cousins - A Contemporary Geometric View Hans Munthe-Kaas 1 Department of Mathematics University of Bergen Norway url: hans.munthe-kaas.no mail: hans.munthe-kaas@uib.no SciCADE Bath,
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationRecent Progress in the Integration of Poisson Systems via the Mid Point Rule and Runge Kutta Algorithm
Recent Progress in te Integration of Poisson Systems via te Mid Point Rule and Runge Kutta Algoritm Klaus Bucner, Mircea Craioveanu and Mircea Puta Abstract Some recent progress in te integration of Poisson
More informationParamodular theta series
Paramodular theta series Rainer Schulze-Pillot Universität des Saarlandes, Saarbrücken 10. 5. 2011 Rainer Schulze-Pillot (Univ. d. Saarlandes) Paramodular theta series 10. 5. 2011 1 / 16 Outline 1 The
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 informationNumerical Analysis: Interpolation Part 1
Numerical Analysis: Interpolation Part 1 Computer Science, Ben-Gurion University (slides based mostly on Prof. Ben-Shahar s notes) 2018/2019, Fall Semester BGU CS Interpolation (ver. 1.00) AY 2018/2019,
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationDifferential Equations
Differential Equations Theory, Technique, and Practice George F. Simmons and Steven G. Krantz Higher Education Boston Burr Ridge, IL Dubuque, IA Madison, Wl New York San Francisco St. Louis Bangkok Bogota
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 informationVariable 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 informationM.SC. PHYSICS - II YEAR
MANONMANIAM SUNDARANAR UNIVERSITY DIRECTORATE OF DISTANCE & CONTINUING EDUCATION TIRUNELVELI 627012, TAMIL NADU M.SC. PHYSICS - II YEAR DKP26 - NUMERICAL METHODS (From the academic year 2016-17) Most Student
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 informationDeferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods
Deferred correction based on exponentially-fitted mono-implicit Runge-Kutta methods M. Van Daele Department of Applied Mathematics, Computer Science and Statistics Ghent University 3th International Conference
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 informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationA 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 information6.3 Partial Fractions
6.3 Partial Fractions Mark Woodard Furman U Fall 2009 Mark Woodard (Furman U) 6.3 Partial Fractions Fall 2009 1 / 11 Outline 1 The method illustrated 2 Terminology 3 Factoring Polynomials 4 Partial fraction
More information1. Consider the initial value problem: find y(t) such that. y = y 2 t, y(0) = 1.
Engineering Mathematics CHEN30101 solutions to sheet 3 1. Consider the initial value problem: find y(t) such that y = y 2 t, y(0) = 1. Take a step size h = 0.1 and verify that the forward Euler approximation
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 informationA Study on Linear and Nonlinear Stiff Problems. Using Single-Term Haar Wavelet Series Technique
Int. Journal of Math. Analysis, Vol. 7, 3, no. 53, 65-636 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.988/ijma.3.3894 A Study on Linear and Nonlinear Stiff Problems Using Single-Term Haar Wavelet Series
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 informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More informationMini project ODE, TANA22
Mini project ODE, TANA22 Filip Berglund (filbe882) Linh Nguyen (linng299) Amanda Åkesson (amaak531) October 2018 1 1 Introduction Carl David Tohmé Runge (1856 1927) was a German mathematician and a prominent
More informationLorentzian elasticity arxiv:
Lorentzian elasticity arxiv:1805.01303 Matteo Capoferri and Dmitri Vassiliev University College London 14 July 2018 Abstract formulation of elasticity theory Consider a manifold M equipped with non-degenerate
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationMathematical Analysis II, 2018/19 First semester
Mathematical Analysis II, 208/9 First semester Yoh Tanimoto Dipartimento di Matematica, Università di Roma Tor Vergata Via della Ricerca Scientifica, I-0033 Roma, Italy email: hoyt@mat.uniroma2.it We basically
More informationLecture 44. Better and successive approximations x2, x3,, xn to the root are obtained from
Lecture 44 Solution of Non-Linear Equations Regula-Falsi Method Method of iteration Newton - Raphson Method Muller s Method Graeffe s Root Squaring Method Newton -Raphson Method An approximation to the
More informationPolytechnic Institute of NYU MA 2132 Final Practice Answers Fall 2012
Polytechnic Institute of NYU MA Final Practice Answers Fall Studying from past or sample exams is NOT recommended. If you do, it should be only AFTER you know how to do all of the homework and worksheet
More informationLINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM
LINEAR ALGEBRA BOOT CAMP WEEK 4: THE SPECTRAL THEOREM Unless otherwise stated, all vector spaces in this worksheet are finite dimensional and the scalar field F is R or C. Definition 1. A linear operator
More informationCDS 101 Precourse Phase Plane Analysis and Stability
CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/
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 information