Energy Preserving Numerical Integration Methods
|
|
- Francis Baker
- 5 years ago
- Views:
Transcription
1 Energy Preserving Numerical Integration Methods Department of Mathematics and Statistics, La Trobe University Supported by Australian Research Council Professorial Fellowship Collaborators: Celledoni (Trondheim) Grimm (Erlangen) McLachlan (Massey) McLaren (La Trobe) O Neale (La Trobe) Owren (Trondheim)
2 Outline of the talk Energy-preserving integrators for ODEs
3 Outline of the talk Energy-preserving integrators for ODEs Energy-preserving integrators for PDEs
4 PART I Energy-preserving integrators for Ordinary Differential Equations
5 Consider a Poisson ODE: dx dt = S(x) H(x). Here S(x) is a skew matrix, and H is the gradient of the Hamiltonian energy function.
6 Consider a Poisson ODE: dx dt = S(x) H(x). Here S(x) is a skew matrix, and H is the gradient of the Hamiltonian energy function. Definition (Skew matrix) S(x) is skew if (a S(x)b) = (b S(x)a) a, b
7 ENERGY-PRESERVING DISCRETE GRADIENT METHOD Definition (Discrete gradient) A discrete gradient H(x n, x n+1 ) is defined by (i) (x n+1 x n ) H(x n, x n+1 ) H(x n+1 ) H(x n ) and (ii) lim xn+1 x n H(x n, x n+1 ) = H(x n )
8 Theorem Let H be a discrete gradient. Then the discretization x n+1 x n t = S(x n ) H(x n, x n+1 ) is energy-preserving.
9 Theorem Let H be a discrete gradient. Then the discretization x n+1 x n t = S(x n ) H(x n, x n+1 ) is energy-preserving. Proof. H(x n+1 ) H(x n ) = H(x n, x n+1 ) (x n+1 x n ) = t H(x n, x n+1 ) S(x n ) H(x n, x n+1 ) = 0
10 TWO EXAMPLES OF DISCRETE GRADIENTS Remember (x n+1 x n ) H(x n, x n+1 ) H(x n+1 ) H(x n ) Example 1. (Itoh-Abe discrete gradient): H := ( H(xn+1,y n) H(x n,y n) x n+1 x n H(x n+1,y n+1 ) H(x n+1,y n) y n+1 y n ). (This can be generalised to any dimension)
11 TWO EXAMPLES OF DISCRETE GRADIENTS Remember (x n+1 x n ) H(x n, x n+1 ) H(x n+1 ) H(x n ) Example 1. (Itoh-Abe discrete gradient): H := ( H(xn+1,y n) H(x n,y n) x n+1 x n H(x n+1,y n+1 ) H(x n+1,y n) y n+1 y n ). (This can be generalised to any dimension) Example 2. ( Average" discrete gradient): H := 1 0 H(ξx n+1 + (1 ξ)x n ) dξ
12 Proof. 1 (x n+1 x n ) H = (x n+1 x n ) H(ξx n+1 + (1 ξ)x n ) dξ 0 1 d = 0 dξ H(ξx n+1 + (1 ξ)x n ) dξ = H(ξx n+1 + (1 ξ)x n ) = H(x n+1 ) H(x n ) ξ=1 ξ=0
13 So x n+1 x n t = S(x) 1 0 H(ξx n+1 + (1 ξ)x n ) dξ is an energy preserving integrator for the ODE dx dt = S(x) H(x)
14 So x n+1 x n t = S(x) 1 0 H(ξx n+1 + (1 ξ)x n ) dξ is an energy preserving integrator for the ODE dx dt = S(x) H(x) From now on we take S(x) = S (constant).
15 So x n+1 x n t = S(x) 1 0 H(ξx n+1 + (1 ξ)x n ) dξ is an energy preserving integrator for the ODE dx dt = S(x) H(x) From now on we take S(x) = S (constant). Then x n+1 x n t x n+1 x n t = S = H(ξx n+1 + (1 ξ)x n ) dξ S H(ξx n+1 + (1 ξ)x n ) dξ
16 So x n+1 x n t = S(x) 1 0 H(ξx n+1 + (1 ξ)x n ) dξ is an energy preserving integrator for the ODE dx dt = S(x) H(x) From now on we take S(x) = S (constant). Then x n+1 x n t x n+1 x n t = S = H(ξx n+1 + (1 ξ)x n ) dξ S H(ξx n+1 + (1 ξ)x n ) dξ But note that S H is the vector field! We get:
17 Theorem (The average vector field method") The numerical integration method x n+1 x n t = 1 0 f(ξx n+1 + (1 ξ)x n ) dξ preserves the energy H(x) exactly for any Hamiltonian ODE with constant symplectic structure, i.e. for dx dt = f(x) with f(x) = S H(x)
18 Example (DOUBLE WELL POTENTIAL) dx dt dy dt = y = 2x(1 x 2 )
19 Example (DOUBLE WELL POTENTIAL) AVF integrator dx dt dy dt = y = 2x(1 x 2 ) ( 1 xn+1 x n t y n+1 y n ) = ( 1 2 (y n+1 + y n ) x n+1 + x n 1 2 (x3 n + x 2 nx n+1 + x n x 2 n+1 + x3 n+1 ) )
20 Energy-preserving integrators for Partial Differential Equations
21 Hamiltonian PDEs Definition The PDE ( u t = f u, u ) x, 2 u x 2,... is Hamiltonian if it is of the form ( f u, u ) x, 2 u x 2,... = Ŝ δh δu δh where Ŝ is a constant skew operator, and δu derivative of H. is the variational
22 Hamiltonian PDEs Definition The PDE ( u t = f u, u ) x, 2 u x 2,... is Hamiltonian if it is of the form ( f u, u ) x, 2 u x 2,... = Ŝ δh δu δh where Ŝ is a constant skew operator, and is the variational δu derivative of H. (Note: f, u, and x can also be taken to be vectors. Although u is usually real, it may also be complex).
23 Hamiltonian PDEs So a Hamiltonian PDE is written u t = Ŝ δh δu
24 Hamiltonian PDEs So a Hamiltonian PDE is written u t = Ŝ δh δu Definition (Variational derivative) Let then δh δu H(u) := H ( u, u ) x, 2 u x 2,... dx H := u ( ) ( ) H + 2 H x u x x 2 u xx
25 Hamiltonian PDEs So a Hamiltonian PDE is written u t = Ŝ δh δu Definition (Variational derivative) Let then δh δu H(u) := H ( u, u ) x, 2 u x 2,... dx H := u ( ) ( ) H + 2 H x u x x 2 u xx Definition (skew operator) Ŝ is skew if a.ŝb = b.ŝa, a, b where a.c := a(x)b(x)dx
26 Example (1. sine-gordon equation (a)) 2 ϕ t 2 = 2 ϕ x 2 + α sin ϕ
27 Example (1. sine-gordon equation (a)) Continuous: H = Ŝ = 2 ϕ t 2 = 2 ϕ x 2 + α sin ϕ [ 1 2 π2 + 1 ( ) ϕ 2 + α(1 cos(ϕ))] dx 2 x ( )
28 Example (1. sine-gordon equation (a)) Continuous: 2 ϕ t 2 = 2 ϕ x 2 + α sin ϕ H = Ŝ = Semi-discrete [ 1 2 π2 + 1 ( ) ϕ 2 + α(1 cos(ϕ))] dx 2 x ( ) ( ϕ x ϕ ) n ϕ n 1 x [ 1 H = x 2 π2 n + n ( 0 id S = id 0 1 2( x) 2 (ϕ n ϕ n 1 ) 2 + α (1 cos(ϕ n )) ), denotes standard gradient ]
29 SG equation (B.D.), Nspace= 101, t =0.01 avf2 Midpoint exact SG equation (B.D.), Nspace= 101, t =0.01 avf2 Midpoint Energy Global Error time time Figure: sine-gordon equation, Backward Differences for spatial discretizations : Energy (left) and global error (right) vs time, for AVF and implicit midpoint integrators.
30 Example (2. sine-gordon equation (b)) 2 ϕ t 2 = 2 ϕ x 2 + α sin ϕ
31 Example (2. sine-gordon equation (b)) Continuous: H = Ŝ = 2 ϕ t 2 = 2 ϕ x 2 + α sin ϕ [ 1 2 π2 + 1 ( ) ϕ 2 + α(1 cos(ϕ))] dx 2 x ( )
32 Example (2. sine-gordon equation (b)) Continuous: 2 ϕ t 2 = 2 ϕ x 2 + α sin ϕ H = Ŝ = Semi-discrete [ 1 2 π2 + 1 ( ) ϕ 2 + α(1 cos(ϕ))] dx 2 x ( ) ( ϕ x ϕ ) n+1 ϕ n 1 2 x [ 1 H = x 2 π2 n + n ( ) 0 id S =. id 0 1 8( x) 2 (ϕ n ϕ n 1 ) 2 + α (1 cos(ϕ n )) ]
33 SG equation (C.D.), Nspace= 101, t =0.01 avf2 Midpoint exact SG equation (C.D.), Nspace= 101, t =0.01 avf2 Midpoint Energy Global Error time time Figure: sine-gordon equation, Central Differences for spatial discretizations : Energy (left) and global error (right) vs time, for AVF and implicit midpoint integrators.
34 Example (3. Korteweg-deVries equation (a)) u t = 6u u x 3 u x 3
35 Example (3. Korteweg-deVries equation (a)) Continuous: u t = 6u u x 3 u x 3 H = Ŝ = x [ ] 1 2 u2 x u 3 dx
36 Example (3. Korteweg-deVries equation (a)) Continuous: u t = 6u u x 3 u x 3 H = Ŝ = x [ ] 1 2 u2 x u 3 dx Semi-discrete: (u 3 u 3 n) H = x [ ] 1 2( x) 2 (u n u n 1 ) 2 u 3 n n S =... 2 x
37 KdV equation, Nspace= 401, t =0.001 avf2 Midpoint exact Energy time Figure: KdV equation: Exact energy, and energy vs time given by AVF and implicit midpoint methods, using discretization (a).
38 Example (4. Korteweg-deVries equation (b)) u t = 6u u x 3 u x 3
39 Example (4. Korteweg-deVries equation (b)) Continuous: u t = 6u u x 3 u x 3 H = Ŝ = x [ ] 1 2 u2 x u 3 dx
40 Example (4. Korteweg-deVries equation (b)) Continuous: u t = 6u u x 3 u x 3 H = Ŝ = x [ ] 1 2 u2 x u 3 dx Semi-discrete: (u 3 u n u n+1 u n+2 ) (for illustrative purposes only) H = x [ ] 1 2( x) 2 (u n u n 1 ) 2 u n u n+1 u n+2 n S =... 2 x
41 Example (4. Korteweg-deVries equation (b)) Continuous: u t = 6u u x 3 u x 3 H = Ŝ = x [ ] 1 2 u2 x u 3 dx Semi-discrete: (u 3 u n u n+1 u n+2 ) (for illustrative purposes only) H = x [ ] 1 2( x) 2 (u n u n 1 ) 2 u n u n+1 u n+2 n S =... 2 x Note: Large freedom in semi-discretisation H without destroying energy preservation. But S must be skew!
42 KdV equation, Nspace= 401, t =0.001 avf2 Midpoint Energy time Figure: KdV equation: Energy vs time given by AVF and implicit midpoint methods, using discretization (b).
43 Example (5. NLS equation) Continuous: ( ) ( u 0 i t u = i 0 ) ( δh δu δh δu ), (1) where u denotes the complex conjugate of u. [ ] H = u 2 x + γ 2 u 4 dx, (2) ( ) 0 i Ŝ =. i 0 (3)
44 Example (5. NLS equation) Continuous: ( ) ( u 0 i t u = i 0 ) ( δh δu δh δu ), (1) where u denotes the complex conjugate of u. [ ] H = u 2 x + γ 2 u 4 dx, (2) ( ) 0 i Ŝ =. i 0 (3) Semi-discrete: H = j [ 1 ( x) 2 u j+1 u j 2 + γ ] 2 u j 4, (4) ( ) 0 id S = i. (5) id 0
45 16 x NLS equation, Nspace= 201, Δt = NLS equation, Nspace= 201, Δt =0.05 avf2 Midpoint Energy Error avf2 Midpoint Global Error time time Figure: Nonlinear Schrödinger equation: Energy error (left) and global error (right) vs time, for AVF and implicit midpoint integrators.
46 14 x 10 6 NLS equation, Nspace= 201, Δt =0.05 Probability Error avf2 Midpt time Figure: Nonlinear Schrödinger equation: Total probability error vs time, for AVF and implicit midpoint integrators.
47 Example (6. Nonlinear Wave Equation) Continuous: H = Ŝ = ϕ t 2 = ( 2 x + 2 y)ϕ ϕ 3 1 ( [ 1 2 (π2 + ϕ 2 x + ϕ 2 y) + 1 ] 4 ϕ4 dx dy )
48 Example (6. Nonlinear Wave Equation) Continuous: H = Ŝ = ϕ t 2 = ( 2 x + 2 y)ϕ ϕ 3 1 ( [ 1 2 (π2 + ϕ 2 x + ϕ 2 y) + 1 ] 4 ϕ4 dx dy ) Semi-discrete: We discretize the Hamiltonian in space with a tensor product Lagrange quadrature formula based on p + 1 Gauss-Lobatto-Legendre (GLL) quadrature nodes in each space direction:
49 Nonlinear Wave Equation II H = 1 2 p j 1 =0 j 2 =0 where d j1,k = dl k(x) dx p w j1 w j2 (π j 2 1,j ( p ) 2 d j1,kϕ k,j2 k=0 ( p m=0 d j2,mϕ j1,m ) ϕ4 j 1,j 2, and l k (x) is the k-th Lagrange basis x=xj1 function based on the GLL quadrature nodes x 0,..., x p, and with w 0,..., w p the corresponding quadrature weights. )
50 T = 0 T = x 10 3 x T = T = x 10 3 x Figure: Snapshots of the solution of the 2D wave equation at different times. AVF method with step-size t = Space discretization with 6 Gauss Lobatto nodes in each space direction. Numerical solution interpolated on a equidistant grid of 21 nodes in each space direction.
51 0.5 x D Wave equation avf and ode15s Energy error time Figure: The 2D wave equation. MATLAB routine ode15s with absolute and relative tolerance (dashed line), and AVF method with step size t = 10/(2 5 ) (solid line). Energy error versus time. Time interval [0, 10]. Space discretization with 6 Gauss Lobatto nodes in each space direction.
52 Leapfrog
53 Exponential Integrator
54 GENERALIZATIONS AND FURTHER REMARKS Instead of using finite differences for the (semi-)discretization of the spatial derivatives one may use e.g. a spectral discretization.
55 GENERALIZATIONS AND FURTHER REMARKS Instead of using finite differences for the (semi-)discretization of the spatial derivatives one may use e.g. a spectral discretization. The method presented also applies to linear PDEs, e.g. the (Linear) Time-dependent Schrödinger equation, 3D Maxwell s equation.
56 GENERALIZATIONS AND FURTHER REMARKS Instead of using finite differences for the (semi-)discretization of the spatial derivatives one may use e.g. a spectral discretization. The method presented also applies to linear PDEs, e.g. the (Linear) Time-dependent Schrödinger equation, 3D Maxwell s equation. The method presented can be generalised to u t = ˆN δh δu where ˆN is a constant negative (semi) definite operator, and where H H is a Lyapunov function, i.e. t 0. (e.g. Allen-Cahn eq, Cahn-Hilliard eq, Ginzburg-Landau eq.)
57 GENERALIZATIONS AND FURTHER REMARKS The method can also be generalised to preserving any number of integrals (not necessarily energy).
58 GENERALIZATIONS AND FURTHER REMARKS The method can also be generalised to preserving any number of integrals (not necessarily energy). If one replaces the integral in the Average Vector Field by quadrature of a certain order, the resulting method will preserve energy to a certain order.
59 GENERALIZATIONS AND FURTHER REMARKS The method can also be generalised to preserving any number of integrals (not necessarily energy). If one replaces the integral in the Average Vector Field by quadrature of a certain order, the resulting method will preserve energy to a certain order. The Average Vector Field Method is a B-series method.
60 Some references to our work ODEs 1. McLachlan, Quispel & Robidoux, A unified approach to Hamiltonian systems, Poisson systems, gradient systems and systems with Lyapunov functions and/or first integrals, Phys. Rev. Lett. 81(1998) McLachlan, Quispel & Robidoux, Geometric integration using discrete gradients, Phil. Trans. Roy. Soc. A357(1999) Quispel & McLaren, A new class of energy-preserving numerical integration methods, J. Phys A41(2008) (7pp). PDEs 4. Celledoni, Grimm, McLachlan, McLaren, O Neale & Quispel, Preserving energy resp. dissipation in numerical PDEs, using the Average Vector Field method, Submitted for publication. Of course there are also many important publications by Furihata, Matsuo, Yaguchi, and collaborators.
norges teknisk-naturvitenskapelige universitet Preserving energy resp. dissipation in numerical PDEs, using the Average Vector Field method
norges teknisk-naturvitenskapelige universitet Preserving energy resp dissipation in numerical PDEs, using the Average Vector Field method by E Celledoni,V Grimm, RI McLachlan, DI McLaren, DRJ O Neale,
More informationarxiv: v1 [math.na] 21 Feb 2012
Preserving energy resp. dissipation in numerical PDEs using the Average Vector Field method E. Celledoni a, V. Grimm b, R. I. McLachlan c, D. I. McLaren d, D. O Neale d,e, B. Owren a, G. R. W. Quispel
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 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 informationENERGY 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 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 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 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 informationDiscrete Variational Derivative Method
.... Discrete Variational Derivative Method one of structure preserving methods for PDEs Daisue Furihata Osaa Univ. 211.7.12 furihata@cmc.osaa-u.ac.jp (Osaa Univ.)Discrete Variational Derivative Method
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationA Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning
Int. J. Contemp. Math. Sciences, Vol. 2, 2007, no. 22, 1097-1106 A Numerical Solution of the Lax s 7 th -order KdV Equation by Pseudospectral Method and Darvishi s Preconditioning M. T. Darvishi a,, S.
More informationEnergy Stable Model Order Reduction for the Allen-Cahn Equation
Energy Stable Model Order Reduction for the Allen-Cahn Equation Bülent Karasözen joint work with Murat Uzunca & Tugba Küçükseyhan Institute of Applied Mathematics & Department of Mathematics, METU, Ankara
More informationFixed point iteration and root finding
Fixed point iteration and root finding The sign function is defined as x > 0 sign(x) = 0 x = 0 x < 0. It can be evaluated via an iteration which is useful for some problems. One such iteration is given
More informationGEOMETRIC 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 informationThe Spectral-Element Method: Introduction
The Spectral-Element Method: Introduction Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Computational Seismology 1 / 59 Outline 1 Introduction 2 Lagrange
More informationThe elliptic sinh-gordon equation in the half plane
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 8 25), 63 73 Research Article The elliptic sinh-gordon equation in the half plane Guenbo Hwang Department of Mathematics, Daegu University, Gyeongsan
More informationA SIXTH ORDER AVERAGED VECTOR FIELD METHOD *
Journal of Computational Mathematics Vol.34, No.5, 6, 479 498. http://www.global-sci.org/jcm doi:.48/jcm.6-m5-65 A SIXTH ORDER AVERAGED VECTOR FIELD METHOD * Haochen Li Jiangsu Key Laboratory for NSLSCS,
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 informationImplicit-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 informationCONSTRUCTING OF CONSTRAINT PRESERVING SCHEME FOR EINSTEIN EQUATIONS
CONSTRUCTING OF CONSTRAINT PRESERVING SCHEME FOR EINSTEIN EQUATIONS TAKUYA TSUCHIYA AND GEN YONEDA arxiv:1610.04543v3 [gr-qc] 14 Jul 2017 Abstract. We propose a new numerical scheme of evolution for the
More informationInfinite-Dimensional Dynamical Systems in Mechanics and Physics
Roger Temam Infinite-Dimensional Dynamical Systems in Mechanics and Physics Second Edition With 13 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vii ix GENERAL
More informationJim Lambers MAT 460/560 Fall Semester Practice Final Exam
Jim Lambers MAT 460/560 Fall Semester 2009-10 Practice Final Exam 1. Let f(x) = sin 2x + cos 2x. (a) Write down the 2nd Taylor polynomial P 2 (x) of f(x) centered around x 0 = 0. (b) Write down the corresponding
More informationKrylov integrators for Hamiltonian systems. Eirola, Timo
https://helda.helsinki.fi Krylov integrators for Hamiltonian systems Eirola, Timo 2018-10-10 Eirola, T & Koskela, A 2018, ' Krylov integrators for Hamiltonian systems ' BIT Numerical Mathematics. https://doi.org/10.1007/s10543-018-0732-y
More informationMA2501 Numerical Methods Spring 2015
Norwegian University of Science and Technology Department of Mathematics MA5 Numerical Methods Spring 5 Solutions to exercise set 9 Find approximate values of the following integrals using the adaptive
More informationMath Numerical Analysis
Math 541 - Numerical Analysis Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center San Diego State University
More informationThe integrating factor method (Sect. 1.1)
The integrating factor method (Sect. 1.1) Overview of differential equations. Linear Ordinary Differential Equations. The integrating factor method. Constant coefficients. The Initial Value Problem. Overview
More informationEnergy-preserving Pseudo-spectral Methods for Klein-Gordon-Schrödinger Equations
.. Energy-preserving Pseudo-spectral Methods for Klein-Gordon-Schrödinger Equations Jingjing Zhang Henan Polytechnic University December 9, 211 Outline.1 Background.2 Aim, Methodology And Motivation.3
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
More informationIntroduction to the Hirota bilinear method
Introduction to the Hirota bilinear method arxiv:solv-int/9708006v1 14 Aug 1997 J. Hietarinta Department of Physics, University of Turku FIN-20014 Turku, Finland e-mail: hietarin@utu.fi Abstract We give
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationSurprising Computations
.... Surprising Computations Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca www.cs.ubc.ca/ ascher/ Uri Ascher (UBC) Surprising Computations Fall 2012 1 / 67 Motivation.
More informationseismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model
seismic wave propagation through spherical sections including anisotropy and visco-elastic dissipation 3D heterogeneities of all Earth model properties inversion of seismic observables for 3D Earth structure
More informationIntroduction of Partial Differential Equations and Boundary Value Problems
Introduction of Partial Differential Equations and Boundary Value Problems 2009 Outline Definition Classification Where PDEs come from? Well-posed problem, solutions Initial Conditions and Boundary Conditions
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationAlexei F. Cheviakov. University of Saskatchewan, Saskatoon, Canada. INPL seminar June 09, 2011
Direct Method of Construction of Conservation Laws for Nonlinear Differential Equations, its Relation with Noether s Theorem, Applications, and Symbolic Software Alexei F. Cheviakov University of Saskatchewan,
More informationHigh-order time-splitting methods for Schrödinger equations
High-order time-splitting methods for Schrödinger equations and M. Thalhammer Department of Mathematics, University of Innsbruck Conference on Scientific Computing 2009 Geneva, Switzerland Nonlinear Schrödinger
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Higher Order
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 informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More informationGeometric-integration tools for the simulation of musical sounds
Japan J. Indust. Appl. Math. https://doi.org/10.1007/s13160-017-09-6 ORIGINAL PAPER Area Geometric-integration tools for the simulation of musical sounds Ai Ishikawa 1 Dominik L. Michels Takaharu Yaguchi
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 informationWRT in 2D: Poisson Example
WRT in 2D: Poisson Example Consider 2 u f on [, L x [, L y with u. WRT: For all v X N, find u X N a(v, u) such that v u dv v f dv. Follows from strong form plus integration by parts: ( ) 2 u v + 2 u dx
More informationHamiltonian partial differential equations and Painlevé transcendents
The 6th TIMS-OCAMI-WASEDA Joint International Workshop on Integrable Systems and Mathematical Physics March 22-26, 2014 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN
More informationWell-balanced DG scheme for Euler equations with gravity
Well-balanced DG scheme for Euler equations with gravity Praveen Chandrashekar praveen@tifrbng.res.in Center for Applicable Mathematics Tata Institute of Fundamental Research Bangalore 560065 Dept. of
More informationContinuous and Discrete Homotopy Operators with Applications in Integrability Testing. Willy Hereman
Continuous and Discrete Homotopy Operators with Applications in Integrability Testing Willy Hereman Department of Mathematical and Computer Sciences Colorado School of Mines Golden, Colorado, USA http://www.mines.edu/fs
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 informationIntroduction and some preliminaries
1 Partial differential equations Introduction and some preliminaries A partial differential equation (PDE) is a relationship among partial derivatives of a function (or functions) of more than one variable.
More information2. Examples of Integrable Equations
Integrable equations A.V.Mikhailov and V.V.Sokolov 1. Introduction 2. Examples of Integrable Equations 3. Examples of Lax pairs 4. Structure of Lax pairs 5. Local Symmetries, conservation laws and the
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More informationLecture 10. (2) Functions of two variables. Partial derivatives. Dan Nichols February 27, 2018
Lecture 10 Partial derivatives Dan Nichols nichols@math.umass.edu MATH 233, Spring 2018 University of Massachusetts February 27, 2018 Last time: functions of two variables f(x, y) x and y are the independent
More informationGalerkin method for the numerical solution of the RLW equation using quintic B-splines
Journal of Computational and Applied Mathematics 19 (26) 532 547 www.elsevier.com/locate/cam Galerkin method for the numerical solution of the RLW equation using quintic B-splines İdris Dağ a,, Bülent
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationDispersion relations, stability and linearization
Dispersion relations, stability and linearization 1 Dispersion relations Suppose that u(x, t) is a function with domain { < x 0}, and it satisfies a linear, constant coefficient partial differential
More informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
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 informationSoliton Solutions of a General Rosenau-Kawahara-RLW Equation
Soliton Solutions of a General Rosenau-Kawahara-RLW Equation Jin-ming Zuo 1 1 School of Science, Shandong University of Technology, Zibo 255049, PR China Journal of Mathematics Research; Vol. 7, No. 2;
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationLecture 1: Water waves and Hamiltonian partial differential equations
Lecture : Water waves and Hamiltonian partial differential equations Walter Craig Department of Mathematics & Statistics Waves in Flows Prague Summer School 208 Czech Academy of Science August 30 208 Outline...
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationNumerical Computational Techniques for Nonlinear Optimal Control
Seminar at LAAS Numerical Computational echniques for Nonlinear Optimal Control Yasuaki Oishi (Nanzan University) July 2, 2018 * Joint work with N. Sakamoto and. Nakamura 1. Introduction Optimal control
More informationA Very General Hamilton Jacobi Theorem
University of Salerno 9th AIMS ICDSDEA Orlando (Florida) USA, July 1 5, 2012 Introduction: Generalized Hamilton-Jacobi Theorem Hamilton-Jacobi (HJ) Theorem: Consider a Hamiltonian system with Hamiltonian
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 informationOutline. 1 Boundary Value Problems. 2 Numerical Methods for BVPs. Boundary Value Problems Numerical Methods for BVPs
Boundary Value Problems Numerical Methods for BVPs Outline Boundary Value Problems 2 Numerical Methods for BVPs Michael T. Heath Scientific Computing 2 / 45 Boundary Value Problems Numerical Methods for
More informationNumerical Analysis Preliminary Exam 10.00am 1.00pm, January 19, 2018
Numerical Analysis Preliminary Exam 0.00am.00pm, January 9, 208 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationSINC PACK, and Separation of Variables
SINC PACK, and Separation of Variables Frank Stenger Abstract This talk consists of a proof of part of Stenger s SINC-PACK computer package (an approx. 400-page tutorial + about 250 Matlab programs) that
More informationHamiltonian partial differential equations and Painlevé transcendents
Winter School on PDEs St Etienne de Tinée February 2-6, 2015 Hamiltonian partial differential equations and Painlevé transcendents Boris DUBROVIN SISSA, Trieste Cauchy problem for evolutionary PDEs with
More informationAP Calculus 2004 AB FRQ Solutions
AP Calculus 4 AB FRQ Solutions Louis A. Talman, Ph. D. Emeritus Professor of Mathematics Metropolitan State University of Denver July, 7 Problem. Part a The function F (t) = 8 + 4 sin(t/) gives the rate,
More information10.34 Numerical Methods Applied to Chemical Engineering. Quiz 2
10.34 Numerical Methods Applied to Chemical Engineering Quiz 2 This quiz consists of three problems worth 35, 35, and 30 points respectively. There are 4 pages in this quiz (including this cover page).
More informationInterpolation Theory
Numerical Analysis Massoud Malek Interpolation Theory The concept of interpolation is to select a function P (x) from a given class of functions in such a way that the graph of y P (x) passes through the
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 informationIBVPs for linear and integrable nonlinear evolution PDEs
IBVPs for linear and integrable nonlinear evolution PDEs Dionyssis Mantzavinos Department of Applied Mathematics and Theoretical Physics, University of Cambridge. Edinburgh, May 31, 212. Dionyssis Mantzavinos
More informationLectures 9-10: Polynomial and piecewise polynomial interpolation
Lectures 9-1: Polynomial and piecewise polynomial interpolation Let f be a function, which is only known at the nodes x 1, x,, x n, ie, all we know about the function f are its values y j = f(x j ), j
More informationIntroduction. J.M. Burgers Center Graduate Course CFD I January Least-Squares Spectral Element Methods
Introduction In this workshop we will introduce you to the least-squares spectral element method. As you can see from the lecture notes, this method is a combination of the weak formulation derived from
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationDispersion in Shallow Water
Seattle University in collaboration with Harvey Segur University of Colorado at Boulder David George U.S. Geological Survey Diane Henderson Penn State University Outline I. Experiments II. St. Venant equations
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 informationSpectral element schemes for the. Korteweg-de Vries and Saint-Venant equations
Spectral element schemes for the Korteweg-de Vries and Saint-Venant equations R. PASQUETTI a a. Université Côte d Azur, CNRS, Inria, LJAD, France, richard.pasquetti@unice.fr... Résumé : Les sytèmes hyperboliques
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More informationA MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC PROBLEMS IN SEMI-INFINITE MEDIA
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 A MODIFIED DECOUPLED SCALED BOUNDARY-FINITE ELEMENT METHOD FOR MODELING 2D IN-PLANE-MOTION TRANSIENT ELASTODYNAMIC
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt =
More informationCOURSE Numerical integration of functions
COURSE 6 3. Numerical integration of functions The need: for evaluating definite integrals of functions that has no explicit antiderivatives or whose antiderivatives are not easy to obtain. Let f : [a,
More informationQualitative behaviour of numerical methods for SDEs and application to homogenization
Qualitative behaviour of numerical methods for SDEs and application to homogenization K. C. Zygalakis Oxford Centre For Collaborative Applied Mathematics, University of Oxford. Center for Nonlinear Analysis,
More informationChapter 1: Introduction
Chapter 1: Introduction Definition: A differential equation is an equation involving the derivative of a function. If the function depends on a single variable, then only ordinary derivatives appear and
More informationNumerical simulation for a nonlinear partial differential equation with variable coefficients by means of the discrete variational derivative method
Journal of Computational and Applied Mathematics 94 (6) 45 459 www.elsevier.com/locate/cam Numerical simulation for a nonlinear partial differential equation with variable coefficients by means of the
More informationCOURSE Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method
COURSE 7 3. Numerical integration of functions (continuation) 3.3. The Romberg s iterative generation method The presence of derivatives in the remainder difficulties in applicability to practical problems
More informationSection 5.8. Taylor Series
Difference Equations to Differential Equations Section 5.8 Taylor Series In this section we will put together much of the work of Sections 5.-5.7 in the context of a discussion of Taylor series. We begin
More informationAM205: Assignment 3 (due 5 PM, October 20)
AM25: Assignment 3 (due 5 PM, October 2) For this assignment, first complete problems 1, 2, 3, and 4, and then complete either problem 5 (on theory) or problem 6 (on an application). If you submit answers
More informationWeak Convergence of Numerical Methods for Dynamical Systems and Optimal Control, and a relation with Large Deviations for Stochastic Equations
Weak Convergence of Numerical Methods for Dynamical Systems and, and a relation with Large Deviations for Stochastic Equations Mattias Sandberg KTH CSC 2010-10-21 Outline The error representation for weak
More informationNATIONAL UNIVERSITY OF SINGAPORE PC5215 NUMERICAL RECIPES WITH APPLICATIONS. (Semester I: AY ) Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE PC55 NUMERICAL RECIPES WITH APPLICATIONS (Semester I: AY 0-) Time Allowed: Hours INSTRUCTIONS TO CANDIDATES. This examination paper contains FIVE questions and comprises
More informationSpectral Elements. Introduction recalling the elastic wave equation
Spectral Elements Introduction recalling the elastic wave equation The spectral-element method: General concept domain mapping from space-continuous to space-discrete time extrapolation Gauss-Lobatto-Legendre
More informationNumerical integration and differentiation. Unit IV. Numerical Integration and Differentiation. Plan of attack. Numerical integration.
Unit IV Numerical Integration and Differentiation Numerical integration and differentiation quadrature classical formulas for equally spaced nodes improper integrals Gaussian quadrature and orthogonal
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS 1. Basic Terminology A differential equation is an equation that contains an unknown function together with one or more of its derivatives. 1 Examples: 1. y = 2x + cos x 2. dy dt
More informationMATH 251 Final Examination December 16, 2015 FORM A. Name: Student Number: Section:
MATH 5 Final Examination December 6, 5 FORM A Name: Student Number: Section: This exam has 7 questions for a total of 5 points. In order to obtain full credit for partial credit problems, all work must
More informationMath 112 Rahman. Week Taylor Series Suppose the function f has the following power series:
Math Rahman Week 0.8-0.0 Taylor Series Suppose the function f has the following power series: fx) c 0 + c x a) + c x a) + c 3 x a) 3 + c n x a) n. ) Can we figure out what the coefficients are? Yes, yes
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More informationDmitry Pelinovsky Department of Mathematics, McMaster University, Canada
Spectrum of the linearized NLS problem Dmitry Pelinovsky Department of Mathematics, McMaster University, Canada Collaboration: Marina Chugunova (McMaster, Canada) Scipio Cuccagna (Modena and Reggio Emilia,
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More information