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
|
|
- Bernice Gallagher
- 6 years ago
- Views:
Transcription
1 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 without the permission of the Author.
2 Travelling Wave Solutions of Multisymplectic Discretisations of Wave Equations A thesis presented in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Mathematics at Massey University, Manawatu, New Zealand. Fleur Cordelia McDonald 2013
3
4 Abstract Symplectic integrators for Hamiltonian ODEs have been well studied over the years and a lot is known about these integrators. They preserve the symplecticity of the system which automatically ensures the preservation of other geometric properties of the system, such as a nearby Hamiltonian and periodic and quasiperiodic orbits. It is then natural to ask how this situation generalises to Hamiltonian PDEs, which leads us to the concept of multisymplectic integration. In this thesis we study the question of how well multisymplectic integrators capture the long-time dynamics of multi-hamiltonian PDEs. We approach this question in two ways numerically and through backward error analysis (BEA). As multi-hamiltonian PDEs possess travelling wave solutions, we wish to see how well multisymplectic integrators preserve these types of solutions. We mainly use the leapfrog method applied to the nonlinear wave equation as our test problem and look for the preservation of periodic travelling waves. We call the resulting equation the discrete travelling wave equation. It cannot be solved exactly. Therefore, our analysis begins with numerically solving the discrete travelling wave equation for simplified nonlinearities. Next, we mov on to analysing periodic solution for a smooth nonlinearity. This results in the presence of resonances in the solutions for certain combinations of the parameters. Finally, we use backward error analysis to compare and back up our results from numerical analysis.
5 Acknowledgements My biggest thanks goes to my supervisor Professor Robert McLachlan. Thank you for taking me on for the past few years and giving me the chance to learn something new. You have been great at steering me in the right direction throughout my research. Thanks for your enthusiasm and perseverance with me. Not only have you helped me with my research but you have helped me grow in other ways by encouraging me to attend many conferences and workshops within and outside of New Zealand. I would not have gone to any of these without your encouragement, but I am very glad that I did. They gave me confidence in presenting, the chance to work with others, and the opportunity to meet many people within the mathematical community. I appreciate all the opportunities you have given to me. Thanks also to the Marsden Fund of the Royal Society of New Zealand for funding my PhD. I would have not been able to continue my studies without this funding. I largely have Associate Professor Brian Moore to thank for the topic of my research. Brian came to visit Massey University in early 2010 from Florida. He was working on travelling waves and while he was here we looked at travelling waves for a semi-discrete semi-linear wave equation using the McKean caricature. We moved on to looking at the fully discrete system with Brian and when he left Robert and I continued on with our study of discrete travelling waves by varying the nonlinearity. Thanks Brian more for giving us the idea for such an interesting topic. Next, I would like to thank Professor Reinout Quispel. It has been good to meet you at various conferences in New Zealand and Australia throughout the years. You always seemed quite interested in our research and were very happy to give your ideas on these. I remember you visiting Massey University soon after Robert found our equation was like a multistep method. It was then that you gave us the idea of applying backward error analysis (BEA) to our equation. I am very grateful for your input over the last few years. Next, I would like to thank my two co-supervisors, Associate Professor ii
6 Matt Perlmutter and Associate Professor Bruce van Brunt. Thanks Matt for being my co-supervisor for the first part of my PhD until you had to move overseas. I know we did not get a chance to talk much while you were my co-supervisor, but thanks for the discussions we had when we met again at a conference down south and then back in Palmerston North. I enjoyed getting to know you during this time. Thank you Bruce for being my co-supervisor throughout the remainder of my studies. You are always very approachable and happy to discuss anything with me. Next, I would like to thank Professor Brynjulf Owren, Professor Elena Celledoni and their students. It was so good to meet you all and spend time with you when you all came to New Zealand in Thank you all for all the research discussions and social catch-ups we had. Finally, I would you like to thank my family, close friends, and my treasured animal friends Coco, Tammy, Tortie, Stripes, Muzbie, Toby, Jake, Ferghus and Flynn, for always supporting me and bringing happiness to my life throughout my many years of university study. iii
7 Contents 1 Geometric Integration of Ordinary Differential Equations Introduction to Ordinary Differential Equations (ODEs) Flow Map Numerical Solutions of ODEs Errors in Numerical Methods Other One-Step Methods Runge Kutta Methods Errors for Runge Kutta Methods Splitting and Composition Methods Hamiltonian Systems N-Body Newtonian Gravitational Problem Molecular Dynamics Poisson Systems Symplecticity Conserved Quantities Symmetries Reversible Differential Equations Geometric Integrators Symplectic Integrators Preserving First Integrals Symmetric and Time-Reversible Integrators N-Body Simulations Using a Symplectic Integrator iv
8 1.9.1 The Leapfrog Method Applied to the Lennard Jones Potential The Outer Solar System Conservation Laws Backward Error Analysis (BEA) Introduction The Modified Differential Equation (MDE) Geometric Properties and the Modified Differential Equation Long-Time Energy Conservation Multisymplectic Integration of Partial Differential Equations Introduction to Partial Differential Equations (PDEs) Hamiltonian PDEs Multi-Hamiltonian PDEs Conservation Laws Solving PDEs Numerically Multisymplectic Integrators The Preissman Box Scheme Constructions of Multisymplectic integrators Discrete Lagrangian Approach Gaussian Runge Kutta Methods Partitioned Runge Kutta methods Behaviour of Multisymplectic Integrators The Dispersion Relation Numerical Experiments Case Study: Numerical Solution of Burgers Equation Burgers Equation A Multisymplectic Scheme for Burgers Equation Solutions to the 8-Point Burgers Scheme (2.26) v
9 3 Introduction to Travelling Waves Travelling Waves Travelling Wave Solutions for the Multi-Hamiltonian Formulation Travelling Wave Solutions for the Nonlinear Wave Equation Travelling Wave Solutions in the Literature Lattice Dynamical Systems (LDSs) Coupled Map Lattices (CMLs) Travelling Wave Solutions for Multisymplectic Discretisations of Wave Equations with Simplified Nonlinearities A Multisymplectic Discretisation of the Nonlinear Wave Equation (3.4) The Discrete Travelling Wave Equation Discrete Travelling Wave Equation with McKean Nonlinearity The Continuous Case Symmetric PDEs and Symmetric Methods: An Example The Semi-Discrete Case The Discrete Case Summary of Results Periodic Travelling Wave Solutions Discrete Travelling Wave Equation with Sawtooth Nonlinearity Analytic Solution Discrete Solution Numerical Solution of the Discrete Travelling Wave Equation with Smooth Nonlinearity The sine-gordon equation Numerical Method for Periodic Travelling Wave Solutions Checking our Numerical Solution Resonance vi
10 5.4 Resonances in the Discrete Travelling Wave Equation Initial Results Measuring the Size of the Resonances Resonance Results Backward Error Analysis of the Discrete Travelling Wave Equation Multistep Methods Linear Multistep Methods for First Order Systems Example of the Behaviour of Multistep methods Multistep Methods for Second Order Equations The Underlying One-Step Method Non-Symplecticity of the Underlying One-Step Method Backward Error Analysis for Multistep Methods for Second Order Differential Equations The (Principal) Modified Equation The Parasitic Modified Equation Conservation Properties from the Underlying One-Step Method The Discrete Travelling Wave Equation as a Multistep Method Backward Error Analysis of the Discrete Travelling Wave Equation Numerical Validation of the BEA The Modified Hamiltonian of the Discrete Travelling Wave Equation Backward Error Analysis for Multisymplectic Integrations of Multi-Hamiltonian PDEs Summary and Conclusions 227 vii
11 List of Figures 1.1 First step in the formation of the rooted tree for the elementary differential f f (f, f) Second step in the formation of the rooted tree for the elementary differential f f (f, f) Final step in the formation of the rooted tree for the elementary differential f f (f, f) Illustration of lower and upper vertices in a rooted tree Illustration showing the leaves of rooted trees Labeling of a rooted tree τ for the related polynomial φ(τ) Labeling of a rooted tree τ for the related polynomial φ(τ), with factors assigned to each label Labeling of a rooted tree τ for the related polynomial φ(τ), with factors assigned to each label and leaf Labeling of the leaves of a rooted tree τ for the related polynomial γ(τ) Labeling of the leaves and other vertices of a rooted tree τ for the related polynomial γ(τ) Flow of a Hamiltonian system, (the Pendulum), showing preservation of areas. This figure is reproduced from [53] The 2-form ω for d = viii
12 1.13 Exact solution of the harmonic oscillator (left), compared with two numerical solutions, the Euler method - solution spirals out clockwise from the centre (middle) and a symplectic integrator (right) Energy error of the Euler method (left) and a symplectic integrator (right) when applied to the harmonic oscillator The attractive and repulsive forces of the Lennard Jones potential Results of the simulation: Snapshots of the dynamics of the two-dimensional Lennard Jones system Energy error of the two-dimensional Lennard Jones system The initial position of the outer planets relative to the sun Orbits of the outer solar system The energy error of the outer solar system after time steps The energy error of the outer solar system after time steps The MDE as the exact solution of a perturbed problem. This diagram is reproduced from [37] Comparison of truncated solutions of the modified differential equation (1.34) and the numerical solution from Euler s method (red) of the equation ẏ = y 2, y(0) = 1. The blue curve gives the truncation after the 1st term (or the exact solution), the green curve the truncation after the 2nd term, and the yellow curve the truncation after the 3rd term Solution of the pendulum and its modified Hamiltonian for the symplectic Euler method Snapshots in time of the leapfrog method applied to the semidiscretisation of the sine-gordon equation Energy error for time steps of the leapfrog method applied to the sine-gordon equation ix
13 2.3 Uniform grid for multisymplectic discretisations The dispersion relation of the linearised wave equation Solution to Burgers equation at different points in time for ν = 1, x = 0.01, and t = Solution to Burgers equation at time t = 0.3 for ν = , x = 0.01, and t = Solitary wave corresponding to a homoclinic orbit Solitary waves corresponding to heteroclinic orbits A homoclinic orbit A heteroclinic orbit FitzHugh-Nagumo travelling wave solutions reproduced from [48] Combustion travelling wave solutions Periodic travelling wave Periodic orbit Phase portrait of the travelling wave equation of the sine-gordon equation Propagation failure for different values of parameters: a versus c for different values of σ, where a (0, 1) and c is the wave speed. This figure is reproduced from [25] Left: A single pulse solution. Right: Zoom on the oscillatory tail from the solution on the left. This figure is reproduced from [49] The five-point stenci used in the leapfrog method Grid points of the discrete travelling wave equation (4.4) with spacing κ (red) and σ (blue) McKean nonlinearity: V (u) = u h(u a) Phase portrait of the travelling wave equation of the nonlinear wave equation with McKean nonlinearity x
14 4.5 Expected travelling wave solution of the nonlinear wave equation with McKean nonlinearity The singularities of equation (4.18) Contour in the upper half plane Travelling wave solution (4.19) of (4.6) with c = 0.1 and ξ = The phase portrait of the symmetric PDE (4.20) Exact solution of the symmetric PDE (4.20) Numerical solution of the symmetric PDE (4.20) with a symmetric method The zeros of R Semi (s) with a 4th order cosine approximation on the left and an 8th order cosine approximation on the right The solution of (4.25) for c = 0: κ = 0.1 left, κ = 0.2 right The solution of (4.25) for σ = κ: c = 0.2 left, c = 0.9 right Zeros of R Disc (s) for σ = 2κ and their periodic images Solution from a real pole for κ = 1 and κ = 0.2, and c = σ Solution from an imaginary pole for κ = 1 and κ = 0.2, and σ 2 c = Solution of the discrete travelling wave equation with McKean nonlinearity for σ = 2κ Expected periodic travelling wave solution of the nonlinear wave equation with McKean nonlinearity Periodic travelling wave solution of the discrete differentialdifference equation (4.25) for parameters τ = 30, σ = 2κ, κ = 0.01, and c = Zoom of Figure 4.20 showing the occurrence of wiggles Zoom of solution for irrational r = 1 2 showing the occurrence of wiggles Wiggles in the periodic travelling wave solution of the nonlinear wave equation with McKean nonlinearity Sawtooth function (4.30) xi
15 4.25 Phase portrait of the travelling wave equation of the nonlinear wave equation with sawtooth nonlinearity Phase portrait and expected solution of the nonlinear wave equation with sawtooth nonlinearity Exact solutions of the nonlinear wave equation with sawtooth nonlinearity Comparison of the analytic solution and the discrete solution of the nonlinear wave equation with sawtooth nonlinearity, for fixed periods and fixed c, for rational values of σ κ Comparison of the analytic solution and the discrete solution of the nonlinear wave equation with sawtooth nonlinearity, for fixed periods and fixed c, for irrational values of σ κ Phase portrait of the travelling wave equation of the nonlinear wave equation with smooth nonlinearity Numerical solution ϕ(ξ) of the discrete travelling wave equation (5.3). ξ is on the x-axis and ϕ on the y-axis Two cases for numerical solutions of the discrete travelling wave equation. ξ is on the x-axis and ϕ on the y-axis. (5.3) Left: An orbit for the leapfrog method applied to the reduced ODE (5.2), and, Right: the corresponding solution in red, with the numerical solution of the discrete travelling wave equation (5.3) for σ = κ = 0.5 superimposed in blue Left: An orbit for the leapfrog method applied to the reduced ODE (5.2), and, Right: the corresponding solution in red, with the numerical solution of the discrete travelling wave equation (5.3) for σ = κ = 0.8 superimposed in blue Continuation in T starting with T = 2π for σ = 1, κ = 1 2, c = Continuation in T starting with T = 2π for σ = 1, κ = 1 2, c = A chosen section of the solution ϕ xii
16 5.9 A short continuation simulation showing one major resonance peak. The y-axis gives the size of the resonance defined by (5.6) Comparison of the resonances for a rational ratio of σ and an κ irrational one. The y-axis gives the size of the resonance defined by (5.6) Continuation in T starting with T = 2π for different rational values of σ. The y-axis gives the size of the resonance defined κ by (5.6) Continuation simulation for σ = 2κ: resonances plotted against T and T. The y-axis gives the size of the resonance defined by 2κ (5.6) Continuation in T starting with T = 2π for different rational values of σ. The y-axis gives the size of the resonance defined κ by (5.6) The left plot shows the resonances for σ = m for the case n = 1 κ n and m = 1, 2, 3, 4 and the plot on the right shows the case n = 2 and m = 3, 5, 7. The y-axis gives the size of the resonance defined by (5.6) Continuation in T starting with T = 2π for fixed rational values of σ κ for different values of σ. The y-axis gives the size of the resonance defined by (5.6) Continuation in T tending to infinity starting with T = 2π. The y-axis gives the size of the resonance defined by (5.6) Continuation in T tending to infinity starting with T = 2π. The y-axis gives the size of the resonance defined by (5.6) The explicit Adams method with step size h = 0.5, initial value (q 0, p 0 ) = (0.7, 0) and the explicit midpoint rule with h = 0.4, initial value (q 0, p 0 ) = (0, 1.1), applied to the pendulum Energy error for three different linear multistep methods applied to the Kepler problem, reproduced from [37]. The x-axis gives the number of periods, and the y-axis gives the energy error xiii
17 6.3 The solution of the 2-step symmetric explicit midpoint rule plotted with the MDE. The blue solution curve corresponds to the same orbit of the numerical solution (also in blue) Relabeling of the five grid points in (6.15) to match that of a multistep method (6.6) for σ = 2κ. The grid points are given as red dots with the spacing between each equal to κ. The labels of the grid points of (6.15) are in blue and those of (6.6) are in green Roots of the characteristic polynomial of the discrete travelling wave equation for c = 1.3 Case Roots of the characteristic polynomial of the discrete travelling wave equation for c = 0.5 Case 3. and Roots of the characteristic polynomial of the discrete travelling wave equation for c = 0 Case Phase portraits of the solution of the modified differential equation of the discrete travelling wave equation (6.23) Left: Orbits of fixed period T = 2π, for increasing κ = σ. The blue orbit has very small value of σ = 0 and the black orbit has a large value of σ = 1.5. The curves in between these have increasing values of σ from the blue orbit to the black orbit. Right: Orbits of fixed amplitude equal to 1.6, for increasing κ = σ. The blue orbit has very small value of σ = 0 and the orange orbit has a large value of σ = 1.8. The curves in between these have increasing values of σ from the blue orbit to the orange orbit Truncated solutions of the modified differential equation for the discrete travelling wave equation: y 1 is the green curve, y 2, the red curve, and y 3 the orange curve. The numerical solution is given by the blue curve xiv
18 6.11 On the left we show the error (the difference between the numerical solution and the solution of the modified equation) of the three truncations of the modified equation for the discrete travelling wave equation: y 1 -green, y 2 -red, y 3 -orange. On the right is a zoom of the graph on the left Resonances of the numerical solution from Chapter 5 given by the blue curve. Resonances of the modified differential equation for the discrete travelling wave equation given by the red curve Contour plot of the resonances for T = 2π and c = 1.3. The lighter the shading the bigger the resonances are On the left is the modified Hamiltonian with small values of σ and κ and on the right the modified Hamiltonian with large values of σ and κ Modified Hamiltonian of the 1st order nonlinear wave system with change of coordinates Modified Hamiltonians of the leapfrog method applied to the nonlinear wave equation as a 1st order system for small σ and κ on the left and large σ and κ on the right Modified Hamiltonians of the Preissman box scheme applied to the nonlinear wave equation as a 1st order system for small σ and κ on the left and large σ and κ on the right xv
19 List of Tables 1.1 Trees of order 3 with their elementary differentials and coefficients Conservation laws for ODEs and PDEs xvi
Energy Preserving Numerical Integration Methods
Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim)
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 informationANALYTICAL MECHANICS. LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS
ANALYTICAL MECHANICS LOUIS N. HAND and JANET D. FINCH CAMBRIDGE UNIVERSITY PRESS Preface xi 1 LAGRANGIAN MECHANICS l 1.1 Example and Review of Newton's Mechanics: A Block Sliding on an Inclined Plane 1
More informationSummer Review Packet. for students entering. IB Math SL
Summer Review Packet for students entering IB Math SL The problems in this packet are designed to help you review topics that are important to your success in IB Math SL. Please attempt the problems on
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 informationSummer Review Packet. for students entering. AP Calculus BC
Summer Review Packet for students entering AP Calculus BC The problems in this packet are designed to help you review topics that are important to your success in AP Calculus. Please attempt the problems
More informationENERGY PRESERVING INTEGRATION OF BI-HAMILTONIAN PARTIAL DIFFERENTIAL EQUATIONS
TWMS J. App. Eng. Math. V., N.,, pp. 75-86 ENERGY PRESERVING INTEGRATION OF BI-HAMILTONIAN PARTIAL DIFFERENTIAL EQUATIONS B. KARASÖZEN, G. ŞİMŞEK Abstract. The energy preserving average vector field (AVF
More informationHigh performance computing and numerical modeling
High performance computing and numerical modeling Volker Springel Plan for my lectures Lecture 1: Collisional and collisionless N-body dynamics Lecture 2: Gravitational force calculation Lecture 3: Basic
More informationSeparatrix Map Analysis for Fractal Scatterings in Weak Interactions of Solitary Waves
Separatrix Map Analysis for Fractal Scatterings in Weak Interactions of Solitary Waves By Yi Zhu, Richard Haberman, and Jianke Yang Previous studies have shown that fractal scatterings in weak interactions
More informationDifferential Equations with Boundary Value Problems
Differential Equations with Boundary Value Problems John Polking Rice University Albert Boggess Texas A&M University David Arnold College of the Redwoods Pearson Education, Inc. Upper Saddle River, New
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 Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationAM 205 Final Project The N-Body Problem
AM 205 Final Project The N-Body Problem Leah Birch Elizabeth Finn Karen Yu December 14, 2012 Abstract The N-Body Problem can be solved using a variety of numeric integrators. Newton s Law of Universal
More informationDoctor of Philosophy
MAINTAINING A COMMON ARBITRARY UNIT IN SOCIAL MEASUREMENT STEPHEN HUMPHRY 2005 Submitted in fulfillment of the requirements of the degree of Doctor of Philosophy School of Education, Murdoch University,
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 informationCLASSICAL MECHANICS. The author
CLASSICAL MECHANICS Gregory s Classical Mechanics is a major new textbook for undergraduates in mathematics and physics. It is a thorough, self-contained and highly readable account of a subject many students
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 informationFree surface flows over submerged obstructions
Free surface flows over submerged obstructions A thesis submitted to the School of Computing Sciences of the University of East Anglia in partial fulfilment of the requirements for the degree of Doctor
More informationAPPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems
APPLIED PARTIM DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fourth Edition Richard Haberman Department of Mathematics Southern Methodist University PEARSON Prentice Hall PEARSON
More informationPARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICAL METHODS PARTIAL DIFFERENTIAL EQUATIONS I YEAR B.Tech By Mr. Y. Prabhaker Reddy Asst. Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad. SYLLABUS OF MATHEMATICAL
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
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 informationDaba Meshesha Gusu and O.Chandra Sekhara Reddy 1
International Journal of Basic and Applied Sciences Vol. 4. No. 1 2015. Pp.22-27 Copyright by CRDEEP. All Rights Reserved. Full Length Research Paper Solutions of Non Linear Ordinary Differential Equations
More informationIndustrial Rotating Kiln Simulation
Industrial Rotating Kiln Simulation This thesis is presented for the degree of Doctor of Philosophy Faculty of Science University of Technology, Sydney 1999 Submitted by Dennis Van Puyvelde, B. Chem. Eng.
More informationMATHEMATICAL METHODS INTERPOLATION
MATHEMATICAL METHODS INTERPOLATION I YEAR BTech By Mr Y Prabhaker Reddy Asst Professor of Mathematics Guru Nanak Engineering College Ibrahimpatnam, Hyderabad SYLLABUS OF MATHEMATICAL METHODS (as per JNTU
More informationOn the construction of discrete gradients
On the construction of discrete gradients Elizabeth L. Mansfield School of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, U.K. G. Reinout W. Quispel Department
More informationChapter 23. Predicting Chaos The Shift Map and Symbolic Dynamics
Chapter 23 Predicting Chaos We have discussed methods for diagnosing chaos, but what about predicting the existence of chaos in a dynamical system. This is a much harder problem, and it seems that the
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationThe Transition to Chaos
Linda E. Reichl The Transition to Chaos Conservative Classical Systems and Quantum Manifestations Second Edition With 180 Illustrations v I.,,-,,t,...,* ', Springer Dedication Acknowledgements v vii 1
More informationDynamics and Time Integration. Computer Graphics CMU /15-662, Fall 2016
Dynamics and Time Integration Computer Graphics CMU 15-462/15-662, Fall 2016 Last Quiz: IK (and demo) Last few classes Keyframing and motion capture - Q: what can one NOT do with these techniques? The
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
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 informationIntroduction to Applied Nonlinear Dynamical Systems and Chaos
Stephen Wiggins Introduction to Applied Nonlinear Dynamical Systems and Chaos Second Edition With 250 Figures 4jj Springer I Series Preface v L I Preface to the Second Edition vii Introduction 1 1 Equilibrium
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
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 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 informationENERGY PRESERVATION (DISSIPATION) IN NUMERICAL ODEs and PDEs
ENERGY PRESERVATION (DISSIPATION) IN NUMERICAL ODEs and PDEs E. CELLEDONI Department of Mathematical Sciences, NTNU 749 Trondheim, Norway. email: elenac@math.ntnu.no Abstract This note considers energy-preservation
More informationIntroduction to Hamiltonian systems
Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008 Examples Mathematical biology: Lotka-Volterra model First numerical methods
More informationTowards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity
Towards Discrete Exterior Calculus and Discrete Mechanics for Numerical Relativity Melvin Leok Mathematics, University of Michigan, Ann Arbor. Joint work with Mathieu Desbrun, Anil Hirani, and Jerrold
More informationBoundary. DIFFERENTIAL EQUATIONS with Fourier Series and. Value Problems APPLIED PARTIAL. Fifth Edition. Richard Haberman PEARSON
APPLIED PARTIAL DIFFERENTIAL EQUATIONS with Fourier Series and Boundary Value Problems Fifth Edition Richard Haberman Southern Methodist University PEARSON Boston Columbus Indianapolis New York San Francisco
More informationSymplectic Runge-Kutta methods satisfying effective order conditions
/30 Symplectic Runge-Kutta methods satisfying effective order conditions Gulshad Imran John Butcher (supervisor) University of Auckland 5 July, 20 International Conference on Scientific Computation and
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
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 informationEstimation for state space models: quasi-likelihood and asymptotic quasi-likelihood approaches
University of Wollongong Research Online University of Wollongong Thesis Collection 1954-2016 University of Wollongong Thesis Collections 2008 Estimation for state space models: quasi-likelihood and asymptotic
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationOrbital and Celestial Mechanics
Orbital and Celestial Mechanics John P. Vinti Edited by Gim J. Der TRW Los Angeles, California Nino L. Bonavito NASA Goddard Space Flight Center Greenbelt, Maryland Volume 177 PROGRESS IN ASTRONAUTICS
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 informationUncommon (in)stability properties in the numerical method of characteristic applied to hyperbolic equations
University of Vermont ScholarWorks @ UVM Graduate College Dissertations and Theses Dissertations and Theses 2016 Uncommon (in)stability properties in the numerical method of characteristic applied to hyperbolic
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationImmaculate Heart Academy Summer Math Assignment for Algebra II Honors and Algebra II/Pre- Calculus- STEM COURSE CODES (5330 and STEM)
Immaculate Heart Academy Summer Math Assignment for Algebra II Honors and Algebra II/Pre- Calculus- STEM COURSE CODES (50 and 551- STEM) LEARN PRACTICE EXCEL You are taking Algebra II Honors or Algebra
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationApplication of a Helmholtz resonator excited by grazing flow for manipulation of a turbulent boundary layer
Application of a Helmholtz resonator excited by grazing flow for manipulation of a turbulent boundary layer Farzin Ghanadi School of Mechanical Engineering The University of Adelaide South Australia, 5005
More informationContents. I Basic Methods 13
Preface xiii 1 Introduction 1 I Basic Methods 13 2 Convergent and Divergent Series 15 2.1 Introduction... 15 2.1.1 Power series: First steps... 15 2.1.2 Further practical aspects... 17 2.2 Differential
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More information5. Ordinary Differential Equations. Indispensable for many technical applications!
Indispensable for many technical applications! Numerisches Programmieren, Hans-Joachim Bungartz page 1 of 30 5.1. Introduction Differential Equations One of the most important fields of application of
More informationTyn Myint-U Lokenath Debnath. Linear Partial Differential Equations for Scientists and Engineers. Fourth Edition. Birkhauser Boston Basel Berlin
Tyn Myint-U Lokenath Debnath Linear Partial Differential Equations for Scientists and Engineers Fourth Edition Birkhauser Boston Basel Berlin Preface to the Fourth Edition Preface to the Third Edition
More informationADVANCED ENGINEERING MATHEMATICS
ADVANCED ENGINEERING MATHEMATICS DENNIS G. ZILL Loyola Marymount University MICHAEL R. CULLEN Loyola Marymount University PWS-KENT O I^7 3 PUBLISHING COMPANY E 9 U Boston CONTENTS Preface xiii Parti ORDINARY
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 informationSOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS
SOLVING ODE s NUMERICALLY WHILE PRESERVING ALL FIRST INTEGRALS G.R.W. QUISPEL 1,2 and H.W. CAPEL 3 Abstract. Using Discrete Gradient Methods (Quispel & Turner, J. Phys. A29 (1996) L341-L349) we construct
More informationGeometric Integrators for ODEs
Geometric Integrators for ODEs Robert I McLachlan and G Reinout W Quispel Massey University, Palmerston North, New Zealand La Trobe University, Bundoora, VIC 3083, Australia Abstract. Geometric integration
More informationKink, singular soliton and periodic solutions to class of nonlinear equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10 Issue 1 (June 015 pp. 1 - Applications and Applied Mathematics: An International Journal (AAM Kink singular soliton and periodic
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 informationxvi xxiii xxvi Construction of the Real Line 2 Is Every Real Number Rational? 3 Problems Algebra of the Real Numbers 7
About the Author v Preface to the Instructor xvi WileyPLUS xxii Acknowledgments xxiii Preface to the Student xxvi 1 The Real Numbers 1 1.1 The Real Line 2 Construction of the Real Line 2 Is Every Real
More informationAS 1 Math Structure for BSc (Ed) (Primary 2 CS Track) AS 1 Math Structure for BSc (Ed) (Secondary)
ACADEMIC SUBJECT: MATHEMATICS Table 1: AS 1 Math Structure for BSc (Ed) (Primary 2 CS Track) AS 1 Math Structure for BSc (Ed) (Secondary) Year 1 2 3 4 Course Code Title Course Category No. of AUs Prerequisites
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 informationBD1. Reciprocal space and Fourier transforms. IB Mineral Sciences Module B: Reciprocal Space, Symmetry and Crystallography
Reciprocal space and Fourier transforms The aim of this practical is to give you a chance to develop your intuition about the relationship between real space and reciprocal space. There are two halves
More informationTRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM
TRANSFER TO THE COLLINEAR LIBRATION POINT L 3 IN THE SUN-EARTH+MOON SYSTEM HOU Xi-yun,2 TANG Jing-shi,2 LIU Lin,2. Astronomy Department, Nanjing University, Nanjing 20093, China 2. Institute of Space Environment
More informationAP Calculus AB. Syllabus. Course Overview and Philosophy. Course Planner:
AP Calculus AB Syllabus Course Overview and Philosophy We cover all topics in the Calculus AB topic outline in the AP Calculus Course Description, and we also include integration by parts, shell method,
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationMR. YATES. Vocabulary. Quadratic Cubic Monomial Binomial Trinomial Term Leading Term Leading Coefficient
ALGEBRA II WITH TRIGONOMETRY COURSE OUTLINE SPRING 2009. MR. YATES Vocabulary Unit 1: Polynomials Scientific Notation Exponent Base Polynomial Degree (of a polynomial) Constant Linear Quadratic Cubic Monomial
More informationMODELING BY NONLINEAR DIFFERENTIAL EQUATIONS
MODELING BY NONLINEAR DIFFERENTIAL EQUATIONS Dissipative and Conservative Processes WORLD SCIENTIFIC SERIES ON NONLINEAR SCIENCE Editor: Leon O. Chua University of California, Berkeley Series A. Volume
More informationHigh order multisymplectic Runge Kutta methods
High order multisymplectic Runge Kutta methods Robert I. McLachlan, Brett N. Ryland, and Yajuan Sun Abstract We study the spatial semidiscretizations obtained by applying Runge Kutta (RK) and partitioned
More informationRiemann Solvers and Numerical Methods for Fluid Dynamics
Eleuterio R Toro Riemann Solvers and Numerical Methods for Fluid Dynamics A Practical Introduction With 223 Figures Springer Table of Contents Preface V 1. The Equations of Fluid Dynamics 1 1.1 The Euler
More informationTRAVELLING WAVE SOLUTIONS
TRAVELLING WAVE SOLUTIONS A dissertation submitted to the University of Manchester for the degree of Master of Science in the Faculty of Science and Engineering 2016 Rami Achouri School of Mathematics
More informationCompacton-like solutions in some nonlocal hydrodynamic-type models
Compacton-like solutions in some nonlocal hydrodynamic-type models Vsevolod Vladimirov AGH University of Science and technology, Faculty of Applied Mathematics Protaras, October 26, 2008 WMS AGH Compactons
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 informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationFrom Last Time. Gravitational forces are apparent at a wide range of scales. Obeys
From Last Time Gravitational forces are apparent at a wide range of scales. Obeys F gravity (Mass of object 1) (Mass of object 2) square of distance between them F = 6.7 10-11 m 1 m 2 d 2 Gravitational
More informationIn this lesson, students model filling a rectangular
NATIONAL MATH + SCIENCE INITIATIVE Mathematics Fill It Up, Please Part III Level Algebra or Math at the end of a unit on linear functions Geometry or Math as part of a unit on volume to spiral concepts
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More information5.12 The Aerodynamic Assist Trajectories of Vehicles Propelled by Solar Radiation Pressure References...
1 The Two-Body Problem... 1 1.1 Position of the Problem... 1 1.2 The Conic Sections and Their Geometrical Properties... 12 1.3 The Elliptic Orbits... 20 1.4 The Hyperbolic and Parabolic Trajectories...
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 informationx 2e e 3x 1. Find the equation of the line that passes through the two points 3,7 and 5, 2 slope-intercept form. . Write your final answer in
Algebra / Trigonometry Review (Notes for MAT0) NOTE: For more review on any of these topics just navigate to my MAT187 Precalculus page and check in the Help section for the topic(s) you wish to review!
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 informationModule 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 informationInstructors Manual Algebra and Trigonometry, 2e Cynthia Y. Young
Dear Instructor, Instructors Manual Algebra and Trigonometry, 2e Cynthia Y. Young I hope that when you read my Algebra and Trigonometry book it will be seamless with how you teach. My goal when writing
More informationCAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS
CAM Ph.D. Qualifying Exam in Numerical Analysis CONTENTS Preliminaries Round-off errors and computer arithmetic, algorithms and convergence Solutions of Equations in One Variable Bisection method, fixed-point
More informationMean-Motion Resonance and Formation of Kirkwood Gaps
Yan Wang Project 1 PHYS 527 October 13, 2008 Mean-Motion Resonance and Formation of Kirkwood Gaps Introduction A histogram of the number of asteroids versus their distance from the Sun shows some distinct
More informationComputational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras
Computational Techniques Prof. Dr. Niket Kaisare Department of Chemical Engineering Indian Institute of Technology, Madras Module No. # 07 Lecture No. # 05 Ordinary Differential Equations (Refer Slide
More informationSouthington High School 720 Pleasant Street Southington, CT 06489
BLUE KNIGHTS Southington High School 720 Pleasant Street Southington, CT 06489 Phone: (860) 628-3229 Fax: (860) 628-3397 Home Page: www.southingtonschools.org Principal Brian Stranieri Assistant Principals
More informationMATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS
MATHEMATICAL STRUCTURES IN CONTINUOUS DYNAMICAL SYSTEMS Poisson Systems and complete integrability with applications from Fluid Dynamics E. van Groesen Dept. of Applied Mathematics University oftwente
More informationA SEMI-ANALYTICAL ORBIT PROPAGATOR PROGRAM FOR HIGHLY ELLIPTICAL ORBITS
A SEMI-ANALYTICAL ORBIT PROPAGATOR PROGRAM FOR HIGHLY ELLIPTICAL ORBITS M. Lara, J. F. San Juan and D. Hautesserres Scientific Computing Group and Centre National d Études Spatiales 6th International Conference
More informationTHE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS
THE NUMERICAL TREATMENT OF DIFFERENTIAL EQUATIONS 2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. BY DR. LOTHAR COLLATZ
More informationThis content has been downloaded from IOPscience. Please scroll down to see the full text.
This content has been downloaded from IOPscience. Please scroll down to see the full text. Download details: IP Address: 46.3.203.124 This content was downloaded on 30/12/2017 at 22:16 Please note that
More informationPRINCIPLES OF PHYSICS. \Hp. Ni Jun TSINGHUA. Physics. From Quantum Field Theory. to Classical Mechanics. World Scientific. Vol.2. Report and Review in
LONDON BEIJING HONG TSINGHUA Report and Review in Physics Vol2 PRINCIPLES OF PHYSICS From Quantum Field Theory to Classical Mechanics Ni Jun Tsinghua University, China NEW JERSEY \Hp SINGAPORE World Scientific
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
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 informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationIntegrating-Factor-Based 2-Additive Runge Kutta methods for Advection-Reaction-Diffusion Equations
Integrating-Factor-Based 2-Additive Runge Kutta methods for Advection-Reaction-Diffusion Equations A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements
More information