Exponential integrators for oscillatory second-order differential equations
|
|
- Ruth Watts
- 5 years ago
- Views:
Transcription
1 Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007
2 Outline Motivation Gautschi-type exponential integrators One- and two-step formulations, properties Main result Numerical example: Fermi-Pasta-Ulam problem
3 Motivation: Sine-Gordon equation u = u xx sinu, t [0, T], x Ω = (0, 1), b.c. + i.c. abstract framework: u = Au g(u), A = u xx, g(u) = sinu u V = D(A 1/2) ) = H0(Ω), 1 u L 2 (Ω) energy: H(u, u ) = 1 2 u A1/2 u 2 + G(u) where T T G(u) = g(tu), u dt = sin(tu(x))u(x)dxdt 0 0 Ω = L 2,, =, L 2
4 Assumptions gradient operator g: u 2 C m G(u) C M + A 1/2 u 2 with moderate constancs C m, C M finite energy: H e (u, u ) = 1 2 u A1/2 u 2 = H(u, u ) G(u) moderately bounded whenever H(u, u ) is
5 Some facts A unbouded operator u highly oscillatory in time (temporal derivaties not bounded) u satisfies finite energy assumptions high frequencies cause small time steps energy conservation
6 Spatial discretization: oscillatory ode second-order ode y + Ω 2 N y = g(y), y(0) = y 0, y (0) = y 0 with finite energy Ω N = A 1/2 N, H(y, y ) = 1 2 y yt Ω 2 N y 1 2 K Ω 2 = A N sym. pos. semidef.
7 Gautschi-type exponential integrator y + Ω 2 (t, y)y = g(y), y(0) = y 0, y (0) = y 0 for constant g and Ω, exact solution satisfies y(t+h) 2y(t)+y(t h) = h 2 ψ(hω) ( g Ω 2 y(t) ), ψ(x) = sinc 2 x 2 (variation-of-constants formula) Gautschi-type exponential integrator y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 ) n y n, Ωn = Ω(t n, y n ) (Verlet for Ω = 0)
8 Gautschi-type exponential integrator y + Ω 2 (t, y)y = g(y), y(0) = y 0, y (0) = y 0 for constant g and Ω, exact solution satisfies y(t+h) 2y(t)+y(t h) = h 2 ψ(hω) ( g Ω 2 y(t) ), ψ(x) = sinc 2 x 2 (variation-of-constants formula) Gautschi-type exponential integrator y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 ) n y n, Ωn = Ω(t n, y n ) (Verlet for Ω = 0) choice of g n?
9 Gautschi-type exponential integrator y n+1 2y n + y n 1 = h 2 ψ(hω n ) ( g n Ω 2 n y n ) obvious choice: g n = g(y n ) Gautschi 61 resonance problems for hω k jπ, ω k eigenvalue of Ω n better choice: g n = g(φ(hω n )y n ), φ filter function φ(0) = 1, φ(kπ) = 0, k = 1, 2, 3,... convergence result: (H., Lubich, 99, Grimm 02, 05) Assumptions: g smooth, bounded energy: y n y(t n ) h 2 C(t n ), C(t n ) e tnl
10 Effect of filter function 10 1 without φ with φ
11 Numerical example: time step h = 0.02 Verlet scheme 128 Fourier modes 512 Fourier modes 2048 Fourier modes Gautschi-type exponential integrator
12 One-step formulation rewrite 2nd order ode as system of 1st order odes, apply variation-of-constants formula: exact solution satisfies [ ] y(t + h) y (t + h) =R(hΩ) [ ] y(t) y (t) motivates numerical scheme [ ] [ ] [ yn+1 yn = R(hΩ) + y n+1 y n + t+h t [ Ω 1 ] sin(t + h s)ω g ( y(s) ) ds cos(t + h s)ω 1 2 h2 Ψg(Φy n ) 1 2 h (Ψ 0 g(φy n ) + Ψ 1 g(φy n+1 ) where [ cos hω Ω R(hΩ) := 1 ] sinhω, Φ = φ(hω), Ψ = ψ(hω),... ΩsinhΩ cos hω ] ),
13 Family of exponential integrators where [ yn+1 y n+1 ] = R(hΩ) [ yn y n ] [ h2 Ψg(Φy n ) 1 2 h (Ψ 0 g(φy n ) + Ψ 1 g(φy n+1 ) Φ = φ(hω), Ψ = ψ(hω), Ψ 0 = ψ 0 (hω), Ψ 1 = ψ 1 (hω) assumptions on ψ, ψ, ψ 0, ψ 1 even analytic functions φ(0) = ψ(0) = ψ 0 (0) = ψ 1 (0) = 1 bounded on non-negative real axis ] ),
14 Properties of one-step scheme symmetric if and only if ψ(ξ) = sinc(ξ)ψ 1 (ξ), ψ 0 (ξ) = cos(ξ)ψ 1 (ξ), symmetric methods can be cast into equivalent two-step formulation with starting values y n+1 2 cos hω y n + y n 1 = h 2 Ψg(Φy n ), y 0, y 1 = cos hω y 0 + hsinchω y h2 Ψg(Φy 0 ).
15 One- and two-step formulations, order mollified impulse method (García-Archilla, Sanz-Serna, Skeel, 1998) order 2 as one-step method with particular starting value y0, y 1 as above order 1 as two-step method with exact starting values Gautschi-type integrator (H., Lubich, 1998) order 2 as two-step method for arbitrary starting values close enough to exact solution symmetric one-step formulation leads to ψ1 with singularities at integer multiples of π
16 Symplecticity necessary and sufficient condition for one-step methods being symplecitic (Hairer, Lubich, Wanner, GNI, 2002) ψ(ξ) = sinc(ξ)φ(ξ), however (Hairer, Lubich, 2000) for Ω = ω > 0 and linear problems, i.e. g(y) = By, energy is conserved up to O(h) for all values of hω if and only if ψ(ξ) = sinc 2 (ξ)φ(ξ)
17 Symplecticity necessary and sufficient condition for one-step methods being symplecitic (Hairer, Lubich, Wanner, GNI, 2002) ψ(ξ) = sinc(ξ)φ(ξ), however (Hairer, Lubich, 2000) for Ω = ω > 0 and linear problems, i.e. g(y) = By, energy is conserved up to O(h) for all values of hω if and only if ψ(ξ) = sinc 2 (ξ)φ(ξ) indicates that methods satisfying the latter condition are preferable to symplectic ones
18 Assumptions on filter functions max χ(ξ) M 1, χ = φ, ψ, ψ 0, ψ 1 ξ 0 max φ(ξ) 1 ξ 0 ξ M 2. ( 1 max ξ 0 sin ξ sinc 2ξ ) 2 ψ(ξ) M max ξ 0 ξ sin ξ (sincξ χ(ξ)) M 4, χ = φ, ψ 0, ψ 1 2 [ yn+1 y n+1 ] = R(hΩ) [ yn y n ] [ h2 Ψg(Φy n ) 1 2 h (Ψ 0 g(φy n ) + Ψ 1 g(φy n+1 ) ] ),
19 Theorem (Grimm, H., 2006) Assumptions; exact solution y satisfies finite-energy condition conditions on filter functions are satisfied then y(t n ) y n h 2 C, t 0 t n = t 0 + nh t 0 + T where the constant C depends on T, K, M 1,...,M 4, g, g y, and g yy. additional conditions on filter functions y (t n ) y n hc, t 0 t n = t 0 + nh t 0 + T
20 Numerical examples Fermi-Pasta-Ulam problem stiff harmonic soft nonlinear new choice of filter function (H., Grimm, 06) ψ(ξ) = sinc 3 (ξ), φ(ξ) = sinc(ξ) order two energy conserved up to O(h) for linear problems
21 Gautschi-type methods ψ(ξ) φ(ξ) A sinc 2 ( 1 2ξ) 1 Gautschi, 61 B sinc(ξ) 1 Deuflhard, 79 C sinc(ξ)φ(ξ) sinc(ξ) García-Archilla et al., 96 D sinc 2 ( 1 2 ξ) sinc(ξ)( sin2( 1 2 ξ)) H., Lubich, 98 E sinc 2 (ξ) 1 Hairer, Lubich, 00 G sinc 3 (ξ) sinc(ξ) Grimm, H., 06
22 Maximum error of total energy (A) (B) (C) π 2π 3π 4π π 2π 3π 4π π 2π 3π 4π 0.2 (D) 0.2 (E) 0.2 (G) π 2π 3π 4π π 2π 3π 4π π 2π 3π 4π interval [0, 1000], h = 0.02, error vs. hω
23 Global error at t = (A) 10 0 (B) 10 0 (C) (D) 10 0 (E) 10 0 (G) error vs. step size, ω = 1000
24 Maximum deviation of oscillatory energy (A) (B) (C) π 2π 3π 4π π 2π 3π 4π π 2π 3π 4π 0.2 (D) 0.2 (E) 0.2 (G) π 2π 3π 4π π 2π 3π 4π π 2π 3π 4π interval [0, 1000], h = 0.02, error vs. hω
25 Comments on new method order 2 (according to theorem) nearly conserved energy (for linear problems according to Hairer, Lubich, 00) no resonances for oscillatory energy (surprise, because there is no method which uniformely conserves oscillatory energy on interval of length > 2π for linear problems, Hairer, Lubich 00)
26 Sketch of proof substitute exact solution into numerical scheme defects derive expressions for defects substract numerical solution, obtain error recursion (discrete variation-of-constants formula) use explicit expression for defects to bound all sums arising apply Gronwall Lemma
27 Summary nonsmooth error bounds for family of exponential integrators characterized second order methods in terms of properties of filter functions accuracy in time independent on spatial discretization results valid for abstract ode s suggest new choice of filter function with favorable properties on fpu example
Exponential integrators
Exponential integrators Marlis Hochbruck Heinrich-Heine University Düsseldorf, Germany Alexander Ostermann University Innsbruck, Austria Helsinki, May 25 p.1 Outline Time dependent Schrödinger equations
More informationLecture 4: High oscillations
Geometric Numerical Integration TU München Ernst Hairer January February Lecture 4: High oscillations Table of contents A Fermi Pasta Ulam type problem Application of classical integrators 4 3 Exponential
More informationNumerical Integrators for Highly Oscillatory Hamiltonian Systems: A Review
Numerical Integrators for Highly Oscillatory Hamiltonian Systems: A Review David Cohen 1, Tobias Jahnke 2, Katina Lorenz 1, and Christian Lubich 1 1 Mathematisches Institut, Univiversität Tübingen, 7276
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 informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
More informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationFinal: Solutions Math 118A, Fall 2013
Final: Solutions Math 118A, Fall 2013 1. [20 pts] For each of the following PDEs for u(x, y), give their order and say if they are nonlinear or linear. If they are linear, say if they are homogeneous or
More informationImplicit-explicit variational integration of highly oscillatory problems
Implicit-explicit variational integration of igly oscillatory problems Ari Stern Structured Integrators Worksop April 9, 9 Stern, A., and E. Grinspun. Multiscale Model. Simul., to appear. arxiv:88.39 [mat.na].
More informationConservation of energy, momentum and actions in numerical discretizations of nonlinear wave equations
Conservation of energy, momentum and actions in numerical discretizations of nonlinear wave equations D. COHEN 1, E. HAIRER 2 and CH. LUBICH 3 April 28, 2007 1 Department of Mathematical Sciences, NTNU,
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 informationMini-course on Geometric Numerical Integration
Mini-course on Geometric Numerical Integration Umeå University and University of Innsbruck david.cohen@umu.se http://snovit.math.umu.se/~david/ This mini-course is supported by an initiation grant from
More informationExponential Integration for Hamiltonian Monte Carlo
Wei-Lun Chao WEILUNC@USC.EDU Justin Solomon 2 JUSTIN.SOLOMON@STANFORD.EDU Dominik L. Michels 2 MICHELS@CS.STANFORD.EDU Fei Sha FEISHA@USC.EDU Department of Computer Science, University of Southern California,
More informationA proof for the full Fourier series on [ π, π] is given here.
niform convergence of Fourier series A smooth function on an interval [a, b] may be represented by a full, sine, or cosine Fourier series, and pointwise convergence can be achieved, except possibly at
More informationExponential Runge Kutta methods for the Schrödinger equation
Exponential Runge Kutta methods for the Schrödinger equation Guillaume Dujardin Department of Applied Mathematics and Theoretical Physics, Center for Mathematical Sciences, Wilberforce road, CAMBRIDGE
More informationstiff harmonic soft nonlinear
LONG-TME ENERGY CONSERVATON OF NUMERCAL METHODS FOR OSCLLATORY DFFERENTAL EQUATONS ERNST HARER AND CHRSTAN LUBCH y Abstract. We consider second-order dierential systems where high-frequency oscillations
More informationVibrating-string problem
EE-2020, Spring 2009 p. 1/30 Vibrating-string problem Newton s equation of motion, m u tt = applied forces to the segment (x, x, + x), Net force due to the tension of the string, T Sinθ 2 T Sinθ 1 T[u
More informationA parallel implementation of a two-dimensional fluid laser-plasma integrator for stratified plasma-vacuum systems
A parallel implementation of a two-dimensional fluid laser-plasma integrator for stratified plasma-vacuum systems Ch. Karle a, J. Schweitzer b, M. Hochbruck b, K.H. Spatschek a, a Institut für Theoretische
More informationLecture 9 Infinite Impulse Response Filters
Lecture 9 Infinite Impulse Response Filters Outline 9 Infinite Impulse Response Filters 9 First-Order Low-Pass Filter 93 IIR Filter Design 5 93 CT Butterworth filter design 5 93 Bilinear transform 7 9
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 informationQuasi-stroboscopic averaging: The non-autonomous case
Outline Quasi-stroboscopic averaging: The non-autonomous case Philippe Chartier INRIA and ENS Cachan Bruz Joint work with Ander Murua and Jesus-Maria Sanz-Serna Geometric Numerical Integration, Oberwolfach
More informationWORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS. Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid
WORD SERIES FOR THE ANALYSIS OF SPLITTING SDE INTEGRATORS Alfonso Álamo/J. M. Sanz-Serna Universidad de Valladolid/Universidad Carlos III de Madrid 1 I. OVERVIEW 2 The importance of splitting integrators
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 informationQuasi-stroboscopic averaging: from B-series to numerical methods?
Outline Quasi-stroboscopic averaging: from B-series to numerical methods? Philippe Chartier INRIA, ENS Cachan Bruz, IRMAR, University of Rennes I Joint work with Ander Murua and Jesus-Maria Sanz-Serna
More informationarxiv: v1 [math.na] 23 Jun 2010
Structure preserving Stochastic Impulse Methods for stiff Langevin systems with a uniform global error of order 1 or 1/2 on position Molei Tao, Houman Owhadi, Jerrold E. Marsden arxiv:16.4657v1 math.na
More informationExponential Runge-Kutta methods for parabolic problems
Exponential Runge-Kutta methods for parabolic problems Marlis Hochbruck a,, Alexander Ostermann b a Mathematisches Institut, Heinrich-Heine Universität Düsseldorf, Universitätsstraße 1, D-4225 Düsseldorf,
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationFourier mode dynamics for NLS Synchronization in fiber lasers arrays
Fourier mode dynamics for NLS Synchronization in fiber lasers arrays Nonlinear Schrodinger workshop Heraklio, May 21 2013 Jean-guy Caputo Laboratoire de Mathématiques INSA de Rouen, France A. Aceves, N.
More informationUniformly accurate averaging numerical schemes for oscillatory evolution equations
Uniformly accurate averaging numerical schemes for oscillatory evolution equations Philippe Chartier University of Rennes, INRIA Joint work with M. Lemou (University of Rennes-CNRS), F. Méhats (University
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 informationAccuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations
COMMUICATIOS I COMPUTATIOAL PHYSICS Vol. 5, o. -4, pp. 779-79 Commun. Comput. Phys. February 009 Accuracy Enhancement Using Spectral Postprocessing for Differential Equations and Integral Equations Tao
More informationFirst Order Differential Equations Lecture 3
First Order Differential Equations Lecture 3 Dibyajyoti Deb 3.1. Outline of Lecture Differences Between Linear and Nonlinear Equations Exact Equations and Integrating Factors 3.. Differences between Linear
More informationUnderstand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.
Review Outline To review for the final, look over the following outline and look at problems from the book and on the old exam s and exam reviews to find problems about each of the following topics.. Basics
More informationA note on the uniform perturbation index 1
Rostock. Math. Kolloq. 52, 33 46 (998) Subject Classification (AMS) 65L5 M. Arnold A note on the uniform perturbation index ABSTRACT. For a given differential-algebraic equation (DAE) the perturbation
More informationUNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH
UNIVERSITY of LIMERICK OLLSCOIL LUIMNIGH Faculty of Science and Engineering Department of Mathematics and Statistics END OF SEMESTER ASSESSMENT PAPER MODULE CODE: MA4006 SEMESTER: Spring 2011 MODULE TITLE:
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 informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 5. 18.305 MIT Fall 2005). D. Margetis and R. Rosales MIT, Math. Dept., Cambridge, MA 02139). November 23, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.011: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 2010 ANSWER BOOKLET
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 informationFourier Series. Spectral Analysis of Periodic Signals
Fourier Series. Spectral Analysis of Periodic Signals he response of continuous-time linear invariant systems to the complex exponential with unitary magnitude response of a continuous-time LI system at
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.
More informationThree approaches for the design of adaptive time-splitting methods
Three approaches for the design of adaptive time-splitting methods Mechthild Thalhammer Leopold Franzens Universität Innsbruck, Austria Workshop on Geometric Integration and Computational Mechanics Organisers:
More informationQuestion Paper Code : AEC11T02
Hall Ticket No Question Paper Code : AEC11T02 VARDHAMAN COLLEGE OF ENGINEERING (AUTONOMOUS) Affiliated to JNTUH, Hyderabad Four Year B. Tech III Semester Tutorial Question Bank 2013-14 (Regulations: VCE-R11)
More informationORDINARY DIFFERENTIAL EQUATIONS
ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We
More 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 informationn 1 f n 1 c 1 n+1 = c 1 n $ c 1 n 1. After taking logs, this becomes
Root finding: 1 a The points {x n+1, }, {x n, f n }, {x n 1, f n 1 } should be co-linear Say they lie on the line x + y = This gives the relations x n+1 + = x n +f n = x n 1 +f n 1 = Eliminating α and
More informationLecture No 1 Introduction to Diffusion equations The heat equat
Lecture No 1 Introduction to Diffusion equations The heat equation Columbia University IAS summer program June, 2009 Outline of the lectures We will discuss some basic models of diffusion equations and
More informationFind the Fourier series of the odd-periodic extension of the function f (x) = 1 for x ( 1, 0). Solution: The Fourier series is.
Review for Final Exam. Monday /09, :45-:45pm in CC-403. Exam is cumulative, -4 problems. 5 grading attempts per problem. Problems similar to homeworks. Integration and LT tables provided. No notes, no
More informationEXPONENTIAL ROSENBROCK-TYPE METHODS
EXPONENTIAL ROSENBROCK-TYPE METHODS MARLIS HOCHBRUCK, ALEXANDER OSTERMANN, AND JULIA SCHWEITZER Abstract. We introduce a new class of exponential integrators for the numerical integration of large-scale
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationLECTURE 5: THE METHOD OF STATIONARY PHASE
LECTURE 5: THE METHOD OF STATIONARY PHASE Some notions.. A crash course on Fourier transform For j =,, n, j = x j. D j = i j. For any multi-index α = (α,, α n ) N n. α = α + + α n. α! = α! α n!. x α =
More informationFrequency-Domain C/S of LTI Systems
Frequency-Domain C/S of LTI Systems x(n) LTI y(n) LTI: Linear Time-Invariant system h(n), the impulse response of an LTI systems describes the time domain c/s. H(ω), the frequency response describes the
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationMassachusetts Institute of Technology
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.11: Introduction to Communication, Control and Signal Processing QUIZ 1, March 16, 21 QUESTION BOOKLET
More informationNumerical methods for eigenvalue problems
Numerical methods for eigenvalue problems D. Löchel Supervisors: M. Hochbruck und M. Tokar Mathematisches Institut Heinrich-Heine-Universität Düsseldorf GRK 1203 seminar february 2008 Outline Introduction
More informationRepresenting a Signal
The Fourier Series Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity and timeinvariance of the system and represents the
More informationExponential integration of large systems of ODEs
Exponential integration of large systems of ODEs Jitse Niesen (University of Leeds) in collaboration with Will Wright (Melbourne University) 23rd Biennial Conference on Numerical Analysis, June 2009 Plan
More informationSolutions to Problems in Chapter 4
Solutions to Problems in Chapter 4 Problems with Solutions Problem 4. Fourier Series of the Output Voltage of an Ideal Full-Wave Diode Bridge Rectifier he nonlinear circuit in Figure 4. is a full-wave
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationA Composite Runge Kutta Method for the Spectral Solution of Semilinear PDEs
Journal of Computational Physics 8, 57 67 (00) doi:0.006/jcph.00.77 A Composite Runge Kutta Method for the Spectral Solution of Semilinear PDEs Tobin A. Driscoll Department of Mathematical Sciences, University
More informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationStochastic Particle Dynamics for Unresolved Degrees of Freedom. Sebastian Reich in collaborative work with Colin Cotter (Imperial College)
Stochastic Particle Dynamics for Unresolved Degrees of Freedom Sebastian Reich in collaborative work with Colin Cotter (Imperial College) 1. The Motivation Classical molecular dynamics (MD) as well as
More informationThe Plan. Initial value problems (ivps) for ordinary differential equations (odes) Review of basic methods You are here: Hamiltonian systems
Scientific Computing with Case Studies SIAM Press, 2009 http://www.cs.umd.edu/users/oleary/sccswebpage Lecture Notes for Unit V Solution of Differential Equations Part 2 Dianne P. O Leary c 2008 The Plan
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 informationBessel Functions Michael Taylor. Lecture Notes for Math 524
Bessel Functions Michael Taylor Lecture Notes for Math 54 Contents 1. Introduction. Conversion to first order systems 3. The Bessel functions J ν 4. The Bessel functions Y ν 5. Relations between J ν and
More informationKeble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II
Tomi Johnson 1 Keble College - Hilary 2014 CP3&4: Mathematical methods I&II Tutorial 5 - Waves and normal modes II Prepare full solutions to the problems with a self assessment of your progress on a cover
More informationBOUNDARY PROBLEMS IN HIGHER DIMENSIONS. kt = X X = λ, and the series solutions have the form (for λ n 0):
BOUNDARY PROBLEMS IN HIGHER DIMENSIONS Time-space separation To solve the wave and diffusion equations u tt = c u or u t = k u in a bounded domain D with one of the three classical BCs (Dirichlet, Neumann,
More informationProf. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait
Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use
More informationExponential integrators
Acta Numerica (2), pp. 29 286 c Cambridge University Press, 2 doi:.7/s962492948 Printed in the United Kingdom Exponential integrators Marlis Hochbruck Karlsruher Institut für Technologie, Institut für
More informationHETEROGENEOUS MULTISCALE METHODS FOR MECHANICAL SYSTEMS WITH VIBRATIONS
HETEROGENEOUS MULTISCALE METHODS FOR MECHANICAL SYSTEMS WITH VIBRATIONS M. P. CALVO AND J. M. SANZ-SERNA Abstract. We are concerned with heterogeneous multiscale methods for the numerical integration of
More informationNotes on the Periodically Forced Harmonic Oscillator
Notes on the Periodically orced Harmonic Oscillator Warren Weckesser Math 38 - Differential Equations 1 The Periodically orced Harmonic Oscillator. By periodically forced harmonic oscillator, we mean the
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 informationArchimedes Center for Modeling, Analysis & Computation. Singular solutions in elastodynamics
Archimedes Center for Modeling, Analysis & Computation Singular solutions in elastodynamics Jan Giesselmann joint work with A. Tzavaras (University of Crete and FORTH) Supported by the ACMAC project -
More informationSolutions Preliminary Examination in Numerical Analysis January, 2017
Solutions Preliminary Examination in Numerical Analysis January, 07 Root Finding The roots are -,0, a) First consider x 0 > Let x n+ = + ε and x n = + δ with δ > 0 The iteration gives 0 < ε δ < 3, which
More informationMAXIMAL AVERAGE ALONG VARIABLE LINES. 1. Introduction
MAXIMAL AVERAGE ALONG VARIABLE LINES JOONIL KIM Abstract. We prove the L p boundedness of the maximal operator associated with a family of lines l x = {(x, x 2) t(, a(x )) : t [0, )} when a is a positive
More informationOrdinary Differential Equations
Ordinary Differential Equations Introduction: first order ODE We are given a function f(t,y) which describes a direction field in the (t,y) plane an initial point (t 0,y 0 ) We want to find a function
More informationwe get y 2 5y = x + e x + C: From the initial condition y(0) = 1, we get 1 5 = 0+1+C; so that C = 5. Completing the square to solve y 2 5y = x + e x 5
Math 24 Final Exam Solution 17 December 1999 1. Find the general solution to the differential equation ty 0 +2y = sin t. Solution: Rewriting the equation in the form (for t 6= 0),we find that y 0 + 2 t
More informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationNEW NUMERICAL INTEGRATORS BASED ON SOLVABILITY AND SPLITTING Fernando Casas Universitat Jaume I, Castellón, Spain Fernando.Casas@uji.es (on sabbatical leave at DAMTP, University of Cambridge) Edinburgh,
More informationThree-dimensional relativistic particle-in-cell hybrid code based on an exponential integrator
Three-dimensional relativistic particle-in-cell hybrid code based on an exponential integrator T. Tückmantel, A. Pukhov Institut für Theoretische Physik I Heinrich-Heine Universität Düsseldorf 45, Düsseldorf,
More informationLECTURE 12 Sections Introduction to the Fourier series of periodic signals
Signals and Systems I Wednesday, February 11, 29 LECURE 12 Sections 3.1-3.3 Introduction to the Fourier series of periodic signals Chapter 3: Fourier Series of periodic signals 3. Introduction 3.1 Historical
More informationMAS 315 Waves 1 of 7 Answers to Example Sheet 3. NB Questions 1 and 2 are relevant to resonance - see S2 Q7
MAS 35 Waves o 7 Answers to Example Sheet 3 NB Questions and are relevant to resonance - see S Q7. The CF is A cosnt + B sin nt (i ω n: Try PI y = C cosωt. OK provided C( ω + n = C = GS is y = A cosnt
More informationDIEUDONNE AGBOR AND JAN BOMAN
ON THE MODULUS OF CONTINUITY OF MAPPINGS BETWEEN EUCLIDEAN SPACES DIEUDONNE AGBOR AND JAN BOMAN Abstract Let f be a function from R p to R q and let Λ be a finite set of pairs (θ, η) R p R q. Assume that
More informationA stroboscopic averaging technique for highly-oscillatory Schrödinger equation
A stroboscopic averaging technique for highly-oscillatory Schrödinger equation F. Castella IRMAR and INRIA-Rennes Joint work with P. Chartier, F. Méhats and A. Murua Saint-Malo Workshop, January 26-28,
More informationLagrangian submanifolds and generating functions
Chapter 4 Lagrangian submanifolds and generating functions Motivated by theorem 3.9 we will now study properties of the manifold Λ φ X (R n \{0}) for a clean phase function φ. As shown in section 3.3 Λ
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 informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
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 informationCircular Membranes. Farlow, Lesson 30. November 21, Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431
The Problem Polar coordinates Solving the Problem by Separation of Variables Circular Membranes Farlow, Lesson 30 Department of Mathematical Sciences Florida Atlantic University Boca Raton, FL 33431 November
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
More informationSuperconvergence analysis of multistep collocation method for delay Volterra integral equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 216, pp. 25-216 Superconvergence analysis of multistep collocation method for delay Volterra integral equations
More informationChapter 6: The Laplace Transform. Chih-Wei Liu
Chapter 6: The Laplace Transform Chih-Wei Liu Outline Introduction The Laplace Transform The Unilateral Laplace Transform Properties of the Unilateral Laplace Transform Inversion of the Unilateral Laplace
More informationHomework 6 Math 309 Spring 2016
Homework 6 Math 309 Spring 2016 Due May 18th Name: Solution: KEY: Do not distribute! Directions: No late homework will be accepted. The homework can be turned in during class or in the math lounge in Pedelford
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 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 information