Exponential multistep methods of Adams-type
|
|
- Kory Little
- 5 years ago
- Views:
Transcription
1 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 of the Helmholtz Association
2 Outline Motivation Exponential Adams methods Linearized exponential multistep methods Implementation Application: Simulation of optical resonators 1
3 Motivation general initial value problem u (t) = F ( t, u(t) ), u(t 0 ) = u 0. (in this talk: autonomous problems only) convergence analysis for semilinear initial value problems variation of constants formula u (t) = Au(t) + g ( t, u(t) ), u(0) = u 0 h u(t n+1 ) = e ha u(t n ) + e (h τ)a g ( u(t n + τ) ) dτ 0 t n = t 0 + nh, n = 0, 1,... idea of multistep methods: replace g by interpolation polynomial 2
4 Exponential Adams methods Given approximations u j u(t j ), consider interpolation polynomial p n through ( tn k+1, g(u n k+1 ) ),..., ( t n, g(u n ) ), given by k 1 p n (t n + θh) = G n + j=1 ( 1) j ( θ j where j G n denotes jth backward difference ) j G n, G j = g(u j ) 0 G n = G n, j G n = j 1 G n j 1 G n 1, j = 1, 2,.... replace nonlinearity g by p n in v.o.c formula 3
5 Exponential Adams methods, cont d numerical scheme h u n+1 = e ha u n + e (h τ)a p n (t n + τ)dτ 0 k 1 = u n + hϕ 1 ( ha)f(u n ) + h γ j ( ha) j G n j=1 with weights γ j (z) = ( 1) j 1 0 ( ) θ e (1 θ)z dθ, j 0. j Certaine, 1960, Nørsett 1969, Cox and Matthews 2002 rational variants: Lambert and Sigurdsson 1972, Verwer
6 Exponential Adams methods, cont d in terms of ϕ-functions the weights γ j are ϕ j (z) = γ 1 = ϕ 2 1 e (1 θ)z 0 γ 2 = ϕ ϕ 2 θ j 1 γ 3 = ϕ 4 + ϕ ϕ 2 dθ, j 1, (j 1)! γ 4 = ϕ ϕ ϕ ϕ 2 γ 5 = ϕ 6 + 2ϕ ϕ ϕ ϕ 2 5
7 Examples of exponential Adams methods k = 1: exponential Euler method k = 2: u n+1 = u n + hϕ 1 ( ha)f(u n ) u n+1 = u n + hϕ 1 ( ha)f(u n ) + hϕ 2 ( ha) ( g(u n ) g(u n 1 ) ) interpretation as corrected exponential Euler step implementation can take advantage of j G n =O(h j ) H., Lubich, Selhofer 1998: exp4 methods are not invariant under linearization 6
8 Starting values use interpolation polynomial p in ( t0, g(u 0 ) ),..., ( t k 1, g(u k 1 ) ) within v.o.c. formula over interval of length mh for m = 1,..., k 1 approximate k 1 u m = u 0 + mhϕ 1 ( mha)f(u 0 ) + h σ m,l ( ha) l G 0, Calvo and Palencia, 2006 solve nonlinear system for u 1,..., u k 1 using fixed point iteration l=1 7
9 Linearized exponential multistep methods construction involves two steps linearize u = F(u) in each step at u n u(t n ) to get u (t) = J n u(t) + g n ( u(t) ) with J n = F u (u n), g n (u) = F(u) J n u apply explicit exponential integrator to linearized problem generalization to non-autonomous problems (not in this talk) H., Ostermann, Schweitzer (2006, 2009), Tokman (2006) 8
10 Examples k = 1: linearized exponential Euler method or exponential Rosenbrock Euler method u n+1 = u n + hϕ 1 (hj n )F(u n ) second-order convergent (H., Ostermann, Schweitzer, 2009) k = 2: Tokman 2006 u n+1 = u n + h ϕ 1 (hj n )F(u n ) 2h 3 ϕ 2(hJ n ) ( g n (u n ) g n (u n 1 ) ) third-order convergent (H., Ostermann, 2010) k = 2: u n+1 = u n + hϕ 1 (hj n )F(u n ) 2hϕ 3 (hj n ) ( g n (u n ) g n (u n 1 ) ) third-order convergent (H., Ostermann, 2010) 9
11 Linearized exponential multistep methods general construction of higher order methods v.o.c. formula applied to linearized ode yields exploit u (t) = J n u(t) + g n ( u(t)) h u(t n+1 ) = e hj n u(t n ) + e (h τ)j ( n g n u(tn + τ) ) dτ 0 g n u (u n) = 0, by approximating g n by Hermite interpolation polynomial p n of degree k satisfying p n (t n ) = 0 10
12 Linearized exponential Adams method, cont d Hermite interpolation polynomial p n of degree k interpolating in ( tn k+1, g n (u n k+1 ) ),..., ( t n, g n (u n ) ) and satisfying p n(t n ) = 0 k 1 p n (t n + θh) = G n,n + j=1 ( θ ( 1) j+1 θ j ) j 1 l l G n,n, l=1 where G n,m = g n (u m ) and j G n,m denotes the jth backward difference (w.r.t. m) 11
13 Linearized exponential Adams method, cont d inserting Hermite interpolation polynomial into v.o.c. formula with weights 1 u n+1 = e hj n u n + h e h(1 θ)j n p n (t n + θh)dθ 0 k 1 = u n + hϕ 1 (hj n )F(u n ) + h j=1 γ 2 = 2ϕ 3 γ 3 = 3ϕ 4 ϕ 3 γ 4 = 4ϕ 5 3ϕ ϕ 3 γ j+1 (hj n ) γ 5 = 5ϕ 6 6ϕ ϕ ϕ 3 j 1 l l G n,n l=1 12
14 Convergence assumptions on semilinear initial value problem u (t) = F ( u(t) ) = Au(t) + g ( u(t) ), u(0) = u 0 X Banach space, A : X X e ta X X + t γ A γ e ta X X C γ, γ, t 0. g locally Lipschitz-continuous in a strip along the exact solution f (t) = g ( u(t) ) sufficiently smooth 13
15 Convergence theorem u (t) = Au(t) + g(u(t)) exponential (k + 1)-step Adams method f (t) = g ( u(t) ) satisfies f C k+1 ([0, T ], X ) Then, if the error bound u j u(t j ) V c 0 h k+1, j = 1,..., k, u n u(t n ) V C h k+1 sup f (k+1) (t) 0 t t n holds uniformly in 0 nh T. C = C(T ), independent of n and h H., Ostermann,
16 Convergence theorem u (t) = Au(t) + g(u(t)) linearized exponential k-step Adams method f (t) = g ( u(t) ) satisfies f C k+1 ([0, T ], X ) Then, if the error bound u j u(t j ) V c 0 h k+1, j = 1,..., k 1, u n u(t n ) V C h k+1 sup 0 t t n ( f (k+1) (t) + u (k+1) (t) V ) holds uniformly in 0 nh T. C = C(T ), independent of n and h H., Ostermann,
17 Outline of proof using interpolation in exact data, exact solution satisfies with defect 1 u(t n+1 ) = e hj n u(t n ) + h e h(1 θ)j n p n (t n + θh) + δ n δ n+1 = h e h(1 θ)j ( n f n (t n + θh) p n (t n + θh) ) dθ 0 estimating interpolation error yields the bounds f n (t n + θh) p n (t n + θh) δ n+1 Ch k+2, δ n+1 V Ch k+2 α 15
18 Outline of proof, cont. error recursion 1 e n+1 = e hj n e n + h e h(1 θ)j ( n p n (t n + θh) p n (t n + θh) ) dθ 0 δ n+1 stability is not trivial (H., Ostermann, Schweitzer 2009) employ Lipschitz condition g(u j ) g ( u(t j ) ) L e j V e n V C max e j V + Ch j=1,...,k 1 n 1 j=0 1 t α n j ( ej V + h k+1) stated error bound follows from a discrete Gronwall lemma 16
19 Numerical example U t 2 U x 2 = 1 + Φ(x, t), x, t [0, 1] 1 + U2 hom. Dirichlet b.c., Φ s.t. U(x, t) = x(1 x) e t, N = 200 grid points error h h 2 h 3 h 4 h 5 h step size 10 1 Adams k = 1,..., 6 linearized Adams k = 1,..., 5 17
20 Implementation issues error analysis shows that j G n = O(h j ) for sufficiently smooth solutions Krylov subspace methods become cheaper with increasing j (Tokman, 2006) higher order does not cost very much can be exploited by reformulation of linearized exponential Adams methods k 1 u n+1 = u n + hϕ 1 (hj n )F(u n ) + h j=1 γ j+1 (hj n ) j 1 l l G n,n l=1 k 1 = u n + hϕ 1 (hj n )F(u n ) + h β k,l (hj n ) l G n,n l=1 18
21 Multiple time stepping construction of exponential Adams methods was based on replacing nonlinearity g by local interpolation polynomial p n in v.o.c. formula h u(t n+1 ) = e ha u(t n ) + 0 e (h τ)a g ( u(t n + τ) ) dτ h u n+1 = e ha u n + e (h τ)a p n (t n + τ)dτ 0 different interpretation: u n+1 = y n (h) exact solution of y n(τ) = Ay n (τ) + p n (t n + τ), y n (0) = u n, multiple time stepping: solve this ode by smaller time steps Grote, Mitkova, 2009, 2010 analogously for linearized exponential Adams methods 19
22 Numerical example, 2d U t U = U 2, x [0, 1]2, t [0, 0.2] hom. Dirichlet b.c., N = 75 grid points in each direction 10 4 j = 1 j = j = 3 j = norms of jth backward differences Krylov dimensions j = 0 20
23 Partitioned problems partitioned semilinear problem [ ] [ ] [ v L Z v = w Y B w ] [ a(t, v, w) + b(t, v, w) ] [ v = A w ] + g(t, v, w) with L corresponding to stiff components, B to nonstiff components [ ] [ ] [ ] [ ] L 0 0 Z v a(t, v, w) A =, g(t, v, w) = + Y 0 0 B w b(t, v, w) matrix functions [ φ( ha) = φ(hl) 0 hy φ [1] (hl) φ(0)i ], φ [1] φ(z) φ(0) (z) = z attractive, if number of stiff components small 21
24 Application: Simulation of optical resonators optical discs and ring resonators have become essential building blocks of integrated all-optical circuits filtering multiplexing dispersion compensation advanced sensing challenges: long time simulation necessary to extract full spectral behavior of resonator complex geometries joint work with Kurt Busch, Jens Niegemann, KIT Photonics group and Abdullah Demirel 22
25 1d test problem Maxwell s equation discontinuous Galerkin discretization, order 4 for E x and B z locally refined grid 10 3 p= p= p= error versus coarse grid step size, x fine = x coarse /p 23
26 Mathematical model of ring resonator modeling and spatial discretization by Kurt Busch and Jens Niegemann Maxwell s equation in 2d PML for transparent boundary conditions locally refined grid discontinuous Galerkin discretization 10 d.o.f. per triangle for each of E x, E y, B z total of d.o.f. 300 components are considered as stiff small matrix functions computed via diagonalization work in progress... only preliminary results 24
27 Histogram of radii of triangles incircles Count Radius of incircle 25
28 Eigenvalues
29 Spatial grid gain about 30% of computational time against explicit RK method 27
30 28 Simulation snapshots
31 Concluding remarks Summary new class of exponential multistep methods convergence analysis of two classes of exponential multistep methods link to multiple time stepping or explicit local time stepping by Mitkova and Grote application to optical resonators Outlook optimize methods for ring resonator simulation different implementations more comparisons with existing methods multiple time stepping rational Krylov methods preconditioning 29
Exponential 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 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 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 informationEXPLICIT EXPONENTIAL RUNGE-KUTTA METHODS FOR SEMILINEAR PARABOLIC PROBLEMS
EXPLICIT EXPONENTIAL RUNGE-KUTTA METHODS FOR SEMILINEAR PARABOLIC PROBLEMS MARLIS HOCHBRUCK AND ALEXANDER OSTERMANN, Abstract. The aim of this paper is to analyze explicit exponential Runge-Kutta methods
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 informationLinear algebra for exponential integrators
Linear algebra for exponential integrators Antti Koskela KTH Royal Institute of Technology Beräkningsmatematikcirkus, 28 May 2014 The problem Considered: time integration of stiff semilinear initial value
More informationExponential 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 informationExponential integrators and functions of the matrix exponential
Exponential integrators and functions of the matrix exponential Paul Matthews, Stephen Cox, Hala Ashi and Linda Cummings School of Mathematical Sciences, University of Nottingham, UK Introduction to exponential
More informationImplementation of exponential Rosenbrock-type integrators
Implementation of exponential Rosenbrock-type integrators Marco Caliari a,b, Alexander Ostermann b, a Department of Pure and Applied Mathematics, University of Padua, Via Trieste 63, I-35121 Padova, Italy
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 informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationExponential integrators for oscillatory second-order differential equations
Exponential integrators for oscillatory second-order differential equations Marlis Hochbruck and Volker Grimm Mathematisches Institut Heinrich Heine Universität Düsseldorf Cambridge, March 2007 Outline
More informationMeshfree Exponential Integrators
Meshfree joint work with A. Ostermann (Innsbruck) M. Caliari (Verona) Leopold Franzens Universität Innsbruck Innovative Integrators 3 October 2 Meshfree Problem class: Goal: Time-dependent PDE s with dominating
More informationApproximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method
Approximating the matrix exponential of an advection-diffusion operator using the incomplete orthogonalization method Antti Koskela KTH Royal Institute of Technology, Lindstedtvägen 25, 10044 Stockholm,
More informationDIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS. Eskil Hansen and Alexander Ostermann
Manuscript submitted to Website: http://aimsciences.org AIMS Journals Volume X, Number 0X, XX 00X pp. 1 11 DIMENSION SPLITTING FOR TIME DEPENDENT OPERATORS Eskil Hansen and Alexander Ostermann Institut
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationMcGill University Department of Mathematics and Statistics. Ph.D. preliminary examination, PART A. PURE AND APPLIED MATHEMATICS Paper BETA
McGill University Department of Mathematics and Statistics Ph.D. preliminary examination, PART A PURE AND APPLIED MATHEMATICS Paper BETA 17 August, 2018 1:00 p.m. - 5:00 p.m. INSTRUCTIONS: (i) This paper
More informationEquivalence Theorems and Their Applications
Equivalence Theorems and Their Applications Tan Bui-Thanh, Center for Computational Geosciences and Optimization Institute for Computational Engineering and Sciences (ICES) The University of Texas at Austin,
More informationError analysis of implicit Euler methods for quasilinear hyperbolic evolution equations
Error analysis of implicit Euler methods for quasilinear hyperbolic evolution equations Marlis Hochbruck Tomislav Pažur CRC Preprint 2016/9 March 2016 KARLSRUHE INSTITUTE OF TECHNOLOGY KIT The Research
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationTakens embedding theorem for infinite-dimensional dynamical systems
Takens embedding theorem for infinite-dimensional dynamical systems James C. Robinson Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. E-mail: jcr@maths.warwick.ac.uk Abstract. Takens
More informationExponential Integrators
Exponential Integrators John C. Bowman (University of Alberta) May 22, 2007 www.math.ualberta.ca/ bowman/talks 1 Exponential Integrators Outline Exponential Euler History Generalizations Stationary Green
More informationSolving Ordinary Differential Equations
Solving Ordinary Differential Equations Sanzheng Qiao Department of Computing and Software McMaster University March, 2014 Outline 1 Initial Value Problem Euler s Method Runge-Kutta Methods Multistep Methods
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 informationKrylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations
mathematics Article Krylov Implicit Integration Factor Methods for Semilinear Fourth-Order Equations Michael Machen and Yong-Tao Zhang * Department of Applied and Computational Mathematics and Statistics,
More informationInterval Methods and Taylor Model Methods for ODEs
Interval Methods and Taylor Model Methods for ODEs Markus Neher, Dept. of Mathematics KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 TM VII, Key West KIT University of the State of Baden-Wuerttemberg and Interval
More information1 Ordinary Differential Equations
Ordinary Differential Equations.0 Mathematical Background.0. Smoothness Definition. A function f defined on [a, b] is continuous at ξ [a, b] if lim x ξ f(x) = f(ξ). Remark Note that this implies existence
More informationFDM for wave equations
FDM for wave equations Consider the second order wave equation Some properties Existence & Uniqueness Wave speed finite!!! Dependence region Analytical solution in 1D Finite difference discretization Finite
More informationCo-sponsored School/Workshop on Integrable Systems and Scientific Computing June Exponential integrators for stiff systems
2042-2 Co-sponsored School/Workshop on Integrable Systems and Scientific Computing 15-20 June 2009 Exponential integrators for stiff systems Paul Matthews University of Nottingham U.K. e-mail: Paul.Matthews@nottingham.ac.uk
More informationKrylov methods for the computation of matrix functions
Krylov methods for the computation of matrix functions Jitse Niesen (University of Leeds) in collaboration with Will Wright (Melbourne University) Heriot-Watt University, March 2010 Outline Definition
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationRUNGE-KUTTA TIME DISCRETIZATIONS OF NONLINEAR DISSIPATIVE EVOLUTION EQUATIONS
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 631 640 S 0025-571805)01866-1 Article electronically published on December 19, 2005 RUNGE-KUTTA TIME DISCRETIZATIONS OF NONLINEAR DISSIPATIVE EVOLUTION
More informationExponential Integrators
Exponential Integrators John C. Bowman and Malcolm Roberts (University of Alberta) June 11, 2009 www.math.ualberta.ca/ bowman/talks 1 Outline Exponential Integrators Exponential Euler History Generalizations
More informationEfficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques
1 Efficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques From Resonators to Open Domains Rob Remis Delft University of Technology November 3, 2017 ICERM Brown University 2 Thanks
More informationComputing the Action of the Matrix Exponential
Computing the Action of the Matrix Exponential Nick Higham School of Mathematics The University of Manchester higham@ma.man.ac.uk http://www.ma.man.ac.uk/~higham/ Joint work with Awad H. Al-Mohy 16th ILAS
More informationTest you method and verify the results with an appropriate function, e.g. f(x) = cos 1 + sin
Advanced methods for ODEs - Lab exercises Verona, 20-30 May 2015 1. FFT and IFFT in MATLAB: Use the Fast Fourier Transform to compute compute the k-th derivative of an appropriate function. In MATLAB the
More informationPreliminary Examination, Numerical Analysis, August 2016
Preliminary Examination, Numerical Analysis, August 2016 Instructions: This exam is closed books and notes. The time allowed is three hours and you need to work on any three out of questions 1-4 and any
More informationOrdinary Differential Equations
Ordinary Differential Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro September 19, 2014 1 / 55 Motivation
More informationConvexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls
1 1 Convexity of the Reachable Set of Nonlinear Systems under L 2 Bounded Controls B.T.Polyak Institute for Control Science, Moscow, Russia e-mail boris@ipu.rssi.ru Abstract Recently [1, 2] the new convexity
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
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 informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
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 informationA STUDY OF MULTIGRID SMOOTHERS USED IN COMPRESSIBLE CFD BASED ON THE CONVECTION DIFFUSION EQUATION
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationDynamical Low-Rank Approximation to the Solution of Wave Equations
Dynamical Low-Rank Approximation to the Solution of Wave Equations Julia Schweitzer joint work with Marlis Hochbruck INSTITUT FÜR ANGEWANDTE UND NUMERISCHE MATHEMATIK 1 KIT Universität des Landes Baden-Württemberg
More informationMarlis Hochbruck 1, Michael Hönig 1 and Alexander Ostermann 2
Mathematical Modelling and Numerical Analysis Modélisation Mathématique et Analyse Numérique Will be set by the publisher REGULARIZATION OF NONLINEAR ILL-POSED PROBLEMS BY EXPONENTIAL INTEGRATORS Marlis
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
More informationSemi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations
Semi-implicit Krylov Deferred Correction Methods for Ordinary Differential Equations Sunyoung Bu University of North Carolina Department of Mathematics CB # 325, Chapel Hill USA agatha@email.unc.edu Jingfang
More 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 informationOn the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods
On the efficiency of the Peaceman-Rachford ADI-dG method for wave-type methods Marlis Hochbruck, Jonas Köhler CRC Preprint 2017/34 (revised), March 2018 KARLSRUHE INSTITUTE OF TECHNOLOGY KIT The Research
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
More informationLecture 10: Finite Differences for ODEs & Nonlinear Equations
Lecture 10: Finite Differences for ODEs & Nonlinear Equations J.K. Ryan@tudelft.nl WI3097TU Delft Institute of Applied Mathematics Delft University of Technology 21 November 2012 () Finite Differences
More informationarxiv: v1 [math.na] 6 Nov 2017
Efficient boundary corrected Strang splitting Lukas Einkemmer Martina Moccaldi Alexander Ostermann arxiv:1711.02193v1 [math.na] 6 Nov 2017 Version of November 6, 2017 Abstract Strang splitting is a well
More informationSouthern Methodist University.
Title: Continuous extensions Name: Lawrence F. Shampine 1, Laurent O. Jay 2 Affil./Addr. 1: Department of Mathematics Southern Methodist University Dallas, TX 75275 USA Phone: +1 (972) 690-8439 E-mail:
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 informationPreconditioning for Nonsymmetry and Time-dependence
Preconditioning for Nonsymmetry and Time-dependence Andy Wathen Oxford University, UK joint work with Jen Pestana and Elle McDonald Jeju, Korea, 2015 p.1/24 Iterative methods For self-adjoint problems/symmetric
More information1. Introduction. In the present paper, we consider a nonautonomous differential equation involving a time-dependent linear operator A
SIAM J. NUMER. ANAL. Vol. 44, No. 2, pp. 851 864 c 2006 Society for Industrial and Applied Mathematics A FOURTH-ORDER COMMUTATOR-FREE EXPONENTIAL INTEGRATOR FOR NONAUTONOMOUS DIFFERENTIAL EQUATIONS MECHTHILD
More informationNon-smooth data error estimates for linearly implicit Runge Kutta methods
IMA Journal of Numerical Analysis 2000 20, 167 184 Non-smooth data error estimates for linearly implicit Runge Kutta methods ALEXANDER OSTERMANN AND MECHTHILD THALHAMMER Institut für Mathematik und Geometrie,
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 informationAn Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations.
An Efficient Algorithm Based on Quadratic Spline Collocation and Finite Difference Methods for Parabolic Partial Differential Equations by Tong Chen A thesis submitted in conformity with the requirements
More informationPartial regularity for fully nonlinear PDE
Partial regularity for fully nonlinear PDE Luis Silvestre University of Chicago Joint work with Scott Armstrong and Charles Smart Outline Introduction Intro Review of fully nonlinear elliptic PDE Our result
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationExponential Time Differencing for Stiff Systems
Journal of Computational Physics 176, 43 455 (22) doi:1.16/jcph.22.6995, available online at http://www.idealibrary.com on Exponential Time Differencing for Stiff Systems S. M. Cox and P. C. Matthews School
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
More informationLecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.
Lecture 4 Chapter 4: Lyapunov Stability Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 4 p. 1/86 Autonomous Systems Consider the autonomous system ẋ
More informationu = (A + F )u, u(0) = η, (1.1)
A CONVERGENCE ANALYSIS OF THE PEACEMAN RACHFORD SCHEME FOR SEMILINEAR EVOLUTION EQUATIONS ESKIL HANSEN AND ERIK HENNINGSSON Abstract. The Peaceman Rachford scheme is a commonly used splitting method for
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationDispersive numerical schemes for Schrödinger equations
Dispersive numerical schemes for Schrödinger equations Enrique Zuazua joint work with L. Ignat zuazua@bcamath.org Basque Center for Applied Mathematics (BCAM), Bilbao, Basque Country, Spain IMA Workshop:
More informationA Sobolev trust-region method for numerical solution of the Ginz
A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012
More informationComparing Leja and Krylov approximations of large scale matrix exponentials
Comparing Leja and Krylov approximations of large scale matrix exponentials L. Bergamaschi 1, M. Caliari 2, A. Martínez 2, and M. Vianello 2 1 Dept. of Math. Methods and Models, University of Padova, berga@dmsa.unipd.it
More informationITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS
ITERATIVE METHODS FOR NONLINEAR ELLIPTIC EQUATIONS LONG CHEN In this chapter we discuss iterative methods for solving the finite element discretization of semi-linear elliptic equations of the form: find
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Multistep Methods Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 1. Runge Kutta methods 2. Embedded RK methods
More informationNumerical Methods for Differential Equations
Numerical Methods for Differential Equations Chapter 2: Runge Kutta and Linear Multistep methods Gustaf Söderlind and Carmen Arévalo Numerical Analysis, Lund University Textbooks: A First Course in the
More informationExam in TMA4215 December 7th 2012
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 9 Contact during the exam: Elena Celledoni, tlf. 7359354, cell phone 48238584 Exam in TMA425 December 7th 22 Allowed
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
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 informationNOTES ON LINEAR ODES
NOTES ON LINEAR ODES JONATHAN LUK We can now use all the discussions we had on linear algebra to study linear ODEs Most of this material appears in the textbook in 21, 22, 23, 26 As always, this is a preliminary
More informationDesign of optimal Runge-Kutta methods
Design of optimal Runge-Kutta methods David I. Ketcheson King Abdullah University of Science & Technology (KAUST) D. Ketcheson (KAUST) 1 / 36 Acknowledgments Some parts of this are joint work with: Aron
More information2 Numerical Methods for Initial Value Problems
Numerical Analysis of Differential Equations 44 2 Numerical Methods for Initial Value Problems Contents 2.1 Some Simple Methods 2.2 One-Step Methods Definition and Properties 2.3 Runge-Kutta-Methods 2.4
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationIEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 57, NO. 2, FEBRUARY X/$ IEEE
IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 57, NO. 2, FEBRUARY 2009 475 Time-Domain Simulations of the Nonlinear Maxwell Equations Using Operator-Exponential Methods Martin Pototschnig, Jens Niegemann,
More informationFavourable time integration methods for. Non-autonomous evolution equations
Favourable time integration methods for non-autonomous evolution equations International Conference on Scientific Computation and Differential Equations (SciCADE) Session on Numerical integration of non-autonomous
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationA space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationProjection Methods. Michal Kejak CERGE CERGE-EI ( ) 1 / 29
Projection Methods Michal Kejak CERGE CERGE-EI ( ) 1 / 29 Introduction numerical methods for dynamic economies nite-di erence methods initial value problems (Euler method) two-point boundary value problems
More informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationThe LEM exponential integrator for advection-diffusion-reaction equations
The LEM exponential integrator for advection-diffusion-reaction equations Marco Caliari a, Marco Vianello a,, Luca Bergamaschi b a Dept. of Pure and Applied Mathematics, University of Padova. b Dept. of
More informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
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 informationThe ReLPM Exponential Integrator for FE Discretizations of Advection-Diffusion Equations
The ReLPM Exponential Integrator for FE Discretizations of Advection-Diffusion Equations Luca Bergamaschi 1, Marco Caliari 2, and Marco Vianello 3 1 Dip.to di Met. e Mod. Mat. per le Scienze Appl., Università
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More information