Exponential integrators for semilinear parabolic problems
|
|
- Kathleen Parrish
- 5 years ago
- Views:
Transcription
1 Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p.
2 Outline Exponential integrators general class of methods order conditions for stiff problems explicit exponential Runge-Kutta methods of order 4 summary Implementation preconditioned Krylov subspace methods for matrix functions summary and outlook p.2
3 Exponential integrators recent interest in exponential integrators for semilinear problems u + Au = g(t, u) Cox and Matthews, JCP 2002 Trefethen and Kassam, to appear in SISC Krogstad, Thesis and Preprint 2003 Minchev, Thesis 2004 Berland, Owren, and Skaflestad, Preprint 2004 Hochbruck and Ostermann, APNUM 2004, Preprint 2004 p.3
4 Exponential integrators, cont d numerical experiments in these references: show that exponential integrators outperform standard methods for a number of examples implementation only for small problems or problems with periodic boundary conditions in this talk: present error bounds discuss implementation p.4
5 PART I Exponential Integrators joint work with Alexander Ostermann University of Innsbruck, Austria p.5
6 General methods u + Au = g(u) motivated by the construction of collocation methods u n+ = χ( ha)u n + h U n,i = χ i ( ha)u n + h s b i ( ha)g n,i, G n,i = g(u n,i ) i= s a ij ( ha)g n,j j= (Friedli 978; Strehmel, Weiner 987) A 0: Runge-Kutta method with b i = b i (0), a ij = a ij (0) consistency: χ(0) = χ i (0) = explicit methods: χ (z), a ij (z) 0, i j s p.6
7 General properties simplifying assumptions (preservation of equilibria) s j= a ij (z) = χ i(z) z, s i= b i (z) = χ(z) z equivalent numerical scheme s u n+ = u n + h U n,i = u n + h G n,i = g(u n,i ) i= s j= b i ( ha) ( G n,i Au n ), a ij ( ha) ( G n,j Au n ), p.7
8 for small A 0 Friedli 978 Strehmel, Weiner 987 Order conditions bi-colored trees, same trees as for W-methods (Berland, Owren, Skaflestad, 2004) up to desired order χ(z) = e z, χ i (z) = e c iz b i, a ij linear combination of ϕ,..., ϕ s ϕ j ( ta) = (j )!t j t 0 e (t τ)a τ j dτ, j p.8
9 for small A 0 Friedli 978 Strehmel, Weiner 987 Order conditions bi-colored trees, same trees as for W-methods (Berland, Owren, Skaflestad, 2004) up to desired order χ(z) = e z, χ i (z) = e c iz b i, a ij linear combination of ϕ,..., ϕ s ϕ j ( ta) = (j )!t j t 0 e (t τ)a τ j dτ, j these order conditions are not sufficient for error bounds p.8
10 Analytical framework u (t) + Au(t) = g(u) (X, ) Banach space, D(A) domain of A in X Assumption on A (Henry, Pazy): A : D(A) X X sectorial: A densely defined, closed linear operator on X, satisfying resolvent condition λ (λi A) M λ a a ϑ for à = A + ωi, ω > a, fractional powers Ãα well defined p.9
11 Analytical framework, II Assumption on g: V = D(Ãα ), v V = Ãα v X g(v) g(w) X L v w V framework = stability φ( ta) X X + φ( ta) V V + t γ Ã γ φ( ta) X X C, for 0 γ, φ = exp, ϕ i, b i, a ij in this talk: X = C(Ω), =, D(A) = C 2 (Ω) C 0 (Ω) here α = 0 p.0
12 Error analysis insert exact solution into numerical scheme u(t n + c i h) = e ciha u(t n ) + h a ij ( ha)f(t n + c j h) + n,i, u(t n+ ) = e ha u(t n ) + h b i ( ha)f(t n + c i h) + δ n+ defects n,i = j h j ψ j,i ( ha)f (j ) (t n ) +... ψ j,i ( ha) = ϕ j ( c i ha) c }{{} =: ϕ j,i j i s k= a ik ( ha) c j k (j )! and analogous for δ n+ and ψ j p.
13 Stiff order conditions No. order order condition ψ ( ha) = ψ 2 ( ha) = ψ i ( ha) = ψ 3 ( ha) = i b i( ha) J ψ 2,i ( ha) = ψ 4 ( ha) = i b i( ha) J ψ 3,i ( ha) = i b i( ha) J j a ij( ha) J ψ 2,j ( ha) = i b i( ha)c i K ψ 2,i ( ha) = 0 p.2
14 Main result Theorem (H, Ostermann 2004) Assume stiff order conditions are satisfied up to order p, 2 p 4 and that ψ p (0) = 0. Further assume that the remaining conditions of order p are satisfied in the weaker form with Then b i (0) instead of b i ( ha), 2 i s. u n u(t n ) Ch p, where C is independent of n and h. p.3
15 Main result Theorem (H, Ostermann 2004) Assume stiff order conditions are satisfied up to order p, 2 p 4 and that ψ p (0) = 0. Further assume that the remaining conditions of order p are satisfied in the weaker form with Then b i (0) instead of b i ( ha), 2 i s. u n u(t n ) Ch p, where C is independent of n and h. Replacing A by 0 elsewhere in the order conditions leads to order reductions in general p.3
16 Second order methods order conditions for order two: b ( ha) + b 2 ( ha) = ϕ ( ha) order b 2 ( ha)c 2 = ϕ 2 ( ha) order 2 a 2 ( ha) = c 2 ϕ ( c 2 ha) order 2 yields one-parameter family of second order methods: 0 c 2 c 2 ϕ ( c 2 ha) ϕ ( ha) c 2 ϕ 2 ( ha) c2 ϕ 2 ( ha) 2nd condition can be weakened to b 2 (0)c 2 = ϕ 2 (0) = 2 p.4
17 Second order methods, cont d one-parameter family using only ϕ : 0 c 2 c 2 ϕ ( c 2 ha) ( 2c 2 )ϕ ( ha) 2c2 ϕ ( ha) most attractive choice: c 2 = 2 which yields b = 0 exponential Runge methods p.5
18 Higher order methods order three conditions (no. 5) for s = 3: b 2 ( ha) J c 2 2ϕ 2 ( c 2 ha)+b 3 ( ha) J ψ 2,3 ( ha) = 0 for all J where ψ 2,3 ( ha) = c 2 3ϕ 2 ( c 3 ha) c 2 a 32 ( ha) p.6
19 Higher order methods order three conditions (no. 5) for s = 3: b 2 ( ha) J c 2 2ϕ 2 ( c 2 ha)+b 3 ( ha) J ψ 2,3 ( ha) = 0 for all J where ψ 2,3 ( ha) = c 2 3ϕ 2 ( c 3 ha) c 2 a 32 ( ha) possible solutions b 2 = 0, then ψ 2,3 ( ha) = 0 b 2 = γb 3, then γc 2 2ϕ 2 ( c 2 ha) + ψ 2,3 ( ha) = 0 p.6
20 4-stage explicit exponential methods consider only methods which reduce to classical Runge-Kutta method for A = 0: several options for designing exponential methods we study methods known to perform well in experiments method of Cox and Matthews (2002) method of Krogstad (2003) p.7
21 Order conditions observations taken from order conditions any fourth order method has to involve ϕ, ϕ 2, and ϕ 3 n,2 = O(h 2 ) for any explicit method since b i (0) 0 for all i we have to duplicate nodes to get cancellation of defects in U n,2, U n,3, U n,4 p.8
22 Cox and Matthews method (ETDRK4) ϕ, ϕ,3 2 ϕ,3 ( ϕ0,3 ) 0 ϕ,3 ϕ 3ϕ 2 + 4ϕ 3 2ϕ 2 4ϕ 3 2ϕ 2 4ϕ 3 4ϕ 3 ϕ 2 satisfies (from total of 9 conditions, 6 9 are for order 4) conditions 4, and 6 (ψ 4 (0) = 0) only weakened form of conditions 5 and 9 (arguments of b i evaluated for A = 0) very weak form of conditions 7 and 8 (all arguments evaluated for A = 0) p.9
23 Cox and Matthews method (ETDRK4) Theorem for periodic boundary conditions: order 4 in the worst case: order 2 if A γ J, J = g is bounded: order 2 + γ u p.20
24 Krogstad s method (ETDRK4-B) ϕ,2 2 ϕ,3 ϕ 2,3 ϕ 2,3 ϕ,4 2ϕ 2,4 0 2ϕ 2,4 ϕ 3ϕ 2 + 4ϕ 3 2ϕ 2 4ϕ 3 2ϕ 2 4ϕ 3 ϕ 2 + 4ϕ 3 satisfies (from total of 9 conditions, 6 9 are for order 4) conditions 5, 9, and 6 (ψ 4 (0) = 0) very weak form of conditions 7 and 8 (all arguments evaluated for A = 0) p.2
25 Krogstad s method (ETDRK4-B) Theorem for periodic boundary conditions: order 4 in the worst case: order 3 if A γ J, J = g is bounded: order 3 + γ u p.22
26 Krogstad s method (ETDRK4-B) Theorem for periodic boundary conditions: order 4 in the worst case: order 3 if A γ J, J = g is bounded: order 3 + γ u is it possible to obtain full order 4? p.22
27 Full order four Runge-Kutta-type methods no 4-stage method of order 4 exists p.23
28 Full order four Runge-Kutta-type methods no 4-stage method of order 4 exists new 4th order method with s = 5, c 5 = / ϕ,2 2 ϕ,3 ϕ 2,3 ϕ 2,3 ϕ,4 2ϕ 2,4 ϕ 2,4 ϕ 2,4 2 2 ϕ,5 2a 5,2 a 5,4 a 5,2 a 5,2 4 ϕ 2,5 a 5,2 ϕ 3ϕ 2 + 4ϕ ϕ 2 + 4ϕ 3 4ϕ 2 8ϕ 3 with a 5,2 = 2 ϕ 2,5 ϕ 3,4 + 4 ϕ 2,4 2 ϕ 3,5 p.23
29 Numerical experiment I test problem: u u xx = +u + g(t), Dirichlet b.c. 2 with exact solution x( x)e t, x = /00, L -norm errors 0 0 Runge Heun 0 2 Krogstad Cox Matthews in this case, γ = p.24
30 Numerical experiment I test problem: u u xx = +u + g(t), Dirichlet b.c. 2 with exact solution x( x)e t, x = /200, L -norm errors 0 0 Runge Heun 0 2 Krogstad Cox Matthews in this case, γ = p.24
31 Numerical experiment II test problem: u u xx = 0 u(x, t)dx, Dirichlet b.c. x = /200, L -norm errors 0 0 new method Krogstad 0 2 CoxMatthews h 2.5 h h in this case, γ = / p.25
32 Summary presented error bounds for explicit exponential Runge-Kutta methods for abstract odes u + Au = g(t, u) new order conditions for stiff problems verified sharpness of these bounds numerically References: M. Hochbruck, A. Ostermann, Exponential Runge-Kutta methods for parabolic problems, to appear in APNUM 2004 M. Hochbruck, A. Ostermann, Explicit exponential Runge-Kutta methods for semilinear parabolic problems, Preprint 2004 p.26
33 PART II Implementation joint work with Jasper van den Eshof University of Düsseldorf, Germany p.27
34 Approximation of matrix operators assume A symmetric, positive semidefinite (but A sectorial works as well) yesterday: showed that Krylov approximations of exp( ha)b and ϕ j ( ha)b in each step: one matrix-vector multiplication Ax always converge superlinearly superlinear convergence starts at m ha /2 p.28
35 Approximation of matrix operators assume A symmetric, positive semidefinite (but A sectorial works as well) yesterday: showed that Krylov approximations of exp( ha)b and ϕ j ( ha)b in each step: one matrix-vector multiplication Ax always converge superlinearly superlinear convergence starts at m ha /2 want to apply preconditioning in order to get mesh independent convergence (as for multigrid) p.28
36 A is discretization of elliptic PDE eigenvalues of D Laplacian not of importance for result important eigenvalues exponential function is rapidly increasing only smallest eigenvalues of A are important but these are hard to find by Lanczos (Kuijlaars 2000) p.29
37 A is discretization of elliptic PDE eigenvalues of D Laplacian not of importance for result important eigenvalues exponential function is rapidly increasing only smallest eigenvalues of A are important but these are hard to find by Lanczos (Kuijlaars 2000) Idea: spectral transformation (I + γa) p.29
38 Preconditioning Lanczos approximations Idea: Lanczos process w.r.t. (I + γa) (instead of A) in each step: one linear system solve (I + γa)x = b obtain basis V m, tridiagonal T m approximation y m (τ) = V m exp( τ T m )e, Tm = γ (T m I). Related work by Moret and Novati, 2002 p.30
39 Questions optimal choice of shift parameter γ a posteriori error estimates iterative solution of linear systems as inner iteration of Lanczos process (preconditioned conjugate gradient iteration) Convergence: A priori: mesh independent but sublinear convergence Table to determine optimal γ for desired accuracy (computed with Remez algorithm) p.3
40 Example, D periodic Poisson, dimension 0 5 y m (τ) = V m exp( τ T m )e, γ = τ/0 0 0 Error Iteration Error for τ = /2, τ = /20, τ = /50, τ = /2000 and upper bound (dotted). p.32
41 Numerical experiment: mesh independence d 2 dx 2 + d2 dy 2 discretized on uniform grids, ɛ = 0 8, γ = τ/0 shifted systems solved by SAMG Fortran package dimension Fixed τ = /0 65/3 70/3 78/3 78/3 90/3 9/3 τ = /00 04/5 48/6 89/8 90/8 08/8 08/8 Relaxed τ = /0 43/3 48/3 52/3 53/3 58/3 60/3 τ = /00 62/5 32/6 48/8 55/8 6/8 64/8 Table reports V-cycles/Lanczos iterations p.33
42 Summary Generalized concept for preconditioning matrix functions spectral transformation Worst case convergence bounded by restricted rational approximations Optimal linear solver optimal method Relaxation of inner iteration results in gain up to 40 percent Oscillatory problems? Reference: J. van den Eshof, M. Hochbruck Preconditioning Lanczos approximations to the matrix exponential, submitted to SISC, March 2004 p.34
EXPLICIT 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 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 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
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 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 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 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 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 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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
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 informationFast Iterative Solution of Saddle Point Problems
Michele Benzi Department of Mathematics and Computer Science Emory University Atlanta, GA Acknowledgments NSF (Computational Mathematics) Maxim Olshanskii (Mech-Math, Moscow State U.) Zhen Wang (PhD student,
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 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 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 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 informationBackground. Background. C. T. Kelley NC State University tim C. T. Kelley Background NCSU, Spring / 58
Background C. T. Kelley NC State University tim kelley@ncsu.edu C. T. Kelley Background NCSU, Spring 2012 1 / 58 Notation vectors, matrices, norms l 1 : max col sum... spectral radius scaled integral norms
More informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
More informationNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Solving the nonlinear Schödinger equation using exponential integrators PREPRINT NUMERICS NO.
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Solving the nonlinear Schödinger equation using exponential integrators by Håvard Berland and Bård Skaflestad PREPRINT NUMERICS NO. 5/2005 NORWEGIAN UNIVERSITY
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 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 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 informationMarlis Hochbruck 1, Michael Hönig 1 and Alexander Ostermann 2
ESAIM: M2AN 43 (29) 79 72 DOI:.5/m2an/292 ESAIM: Mathematical Modelling and Numerical Analysis www.esaim-m2an.org REGULARIZATION OF NONLINEAR ILL-POSED PROBLEMS BY EXPONENTIAL INTEGRATORS Marlis Hochbruck,
More informationComputation Fluid Dynamics
Computation Fluid Dynamics CFD I Jitesh Gajjar Maths Dept Manchester University Computation Fluid Dynamics p.1/189 Garbage In, Garbage Out We will begin with a discussion of errors. Useful to understand
More informationRecent advances in approximation using Krylov subspaces. V. Simoncini. Dipartimento di Matematica, Università di Bologna.
Recent advances in approximation using Krylov subspaces V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The framework It is given an operator
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 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 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 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 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 informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
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 informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
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 information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationRitz Value Bounds That Exploit Quasi-Sparsity
Banff p.1 Ritz Value Bounds That Exploit Quasi-Sparsity Ilse Ipsen North Carolina State University Banff p.2 Motivation Quantum Physics Small eigenvalues of large Hamiltonians Quasi-Sparse Eigenvector
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 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 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 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 informationMotivation: Sparse matrices and numerical PDE's
Lecture 20: Numerical Linear Algebra #4 Iterative methods and Eigenproblems Outline 1) Motivation: beyond LU for Ax=b A little PDE's and sparse matrices A) Temperature Equation B) Poisson Equation 2) Splitting
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 informationBoundary Value Problems - Solving 3-D Finite-Difference problems Jacob White
Introduction to Simulation - Lecture 2 Boundary Value Problems - Solving 3-D Finite-Difference problems Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Reminder about
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.4: Krylov Subspace Methods 2 / 51 Krylov Subspace Methods In this chapter we give
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 informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationApproximation of functions of large matrices. Part I. Computational aspects. V. Simoncini
Approximation of functions of large matrices. Part I. Computational aspects V. Simoncini Dipartimento di Matematica, Università di Bologna and CIRSA, Ravenna, Italy valeria@dm.unibo.it 1 The problem Given
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 informationNumerical Study of Oscillatory Regimes in the KP equation
Numerical Study of Oscillatory Regimes in the KP equation C. Klein, MPI for Mathematics in the Sciences, Leipzig, with C. Sparber, P. Markowich, Vienna, math-ph/"#"$"%& C. Sparber (generalized KP), personal-homepages.mis.mpg.de/klein/
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 informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
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 informationThe amount of work to construct each new guess from the previous one should be a small multiple of the number of nonzeros in A.
AMSC/CMSC 661 Scientific Computing II Spring 2005 Solution of Sparse Linear Systems Part 2: Iterative methods Dianne P. O Leary c 2005 Solving Sparse Linear Systems: Iterative methods The plan: Iterative
More informationOn solving linear systems arising from Shishkin mesh discretizations
On solving linear systems arising from Shishkin mesh discretizations Petr Tichý Faculty of Mathematics and Physics, Charles University joint work with Carlos Echeverría, Jörg Liesen, and Daniel Szyld October
More informationSOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS. Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA
1 SOLVING SPARSE LINEAR SYSTEMS OF EQUATIONS Chao Yang Computational Research Division Lawrence Berkeley National Laboratory Berkeley, CA, USA 2 OUTLINE Sparse matrix storage format Basic factorization
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
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 informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationSolving Ill-Posed Cauchy Problems in Three Space Dimensions using Krylov Methods
Solving Ill-Posed Cauchy Problems in Three Space Dimensions using Krylov Methods Lars Eldén Department of Mathematics Linköping University, Sweden Joint work with Valeria Simoncini February 21 Lars Eldén
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
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 informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
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 informationDiscretization of PDEs and Tools for the Parallel Solution of the Resulting Systems
Discretization of PDEs and Tools for the Parallel Solution of the Resulting Systems Stan Tomov Innovative Computing Laboratory Computer Science Department The University of Tennessee Wednesday April 4,
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationMultigrid and Iterative Strategies for Optimal Control Problems
Multigrid and Iterative Strategies for Optimal Control Problems John Pearson 1, Stefan Takacs 1 1 Mathematical Institute, 24 29 St. Giles, Oxford, OX1 3LB e-mail: john.pearson@worc.ox.ac.uk, takacs@maths.ox.ac.uk
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
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 informationMultigrid Methods for Discretized PDE Problems
Towards Metods for Discretized PDE Problems Institute for Applied Matematics University of Heidelberg Feb 1-5, 2010 Towards Outline A model problem Solution of very large linear systems Iterative Metods
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationThe Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge
The Dirichlet problem for non-divergence parabolic equations with discontinuous in time coefficients in a wedge Vladimir Kozlov (Linköping University, Sweden) 2010 joint work with A.Nazarov Lu t u a ij
More informationKrylov single-step implicit integration factor WENO methods for advection-diffusion-reaction equations
Accepted Manuscript Krylov single-step implicit integration factor WENO methods for advection diffusion reaction equations Tian Jiang, Yong-Tao Zhang PII: S0021-9991(16)00029-2 DOI: http://dx.doi.org/10.1016/j.jcp.2016.01.021
More informationFinite Element Methods
Solving Operator Equations Via Minimization We start with several definitions. Definition. Let V be an inner product space. A linear operator L: D V V is said to be positive definite if v, Lv > for every
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationSpectral Methods for Reaction Diffusion Systems
WDS'13 Proceedings of Contributed Papers, Part I, 97 101, 2013. ISBN 978-80-7378-250-4 MATFYZPRESS Spectral Methods for Reaction Diffusion Systems V. Rybář Institute of Mathematics of the Academy of Sciences
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 11 Partial Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002.
More informationLecture 42 Determining Internal Node Values
Lecture 42 Determining Internal Node Values As seen in the previous section, a finite element solution of a boundary value problem boils down to finding the best values of the constants {C j } n, which
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
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 informationKey words. preconditioned conjugate gradient method, saddle point problems, optimal control of PDEs, control and state constraints, multigrid method
PRECONDITIONED CONJUGATE GRADIENT METHOD FOR OPTIMAL CONTROL PROBLEMS WITH CONTROL AND STATE CONSTRAINTS ROLAND HERZOG AND EKKEHARD SACHS Abstract. Optimality systems and their linearizations arising in
More informationHigh order splitting methods for analytic semigroups exist
BIT manuscript No. (will be inserted by the editor High order splitting methods for analytic semigroups exist Eskil Hansen Alexander Ostermann Received: date / Accepted: date Abstract In this paper, we
More informationPOD for Parametric PDEs and for Optimality Systems
POD for Parametric PDEs and for Optimality Systems M. Kahlbacher, K. Kunisch, H. Müller and S. Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria DMV-Jahrestagung 26,
More informationTime-dependent Dirichlet Boundary Conditions in Finite Element Discretizations
Time-dependent Dirichlet Boundary Conditions in Finite Element Discretizations Peter Benner and Jan Heiland November 5, 2015 Seminar Talk at Uni Konstanz Introduction Motivation A controlled physical processes
More information8.1 Concentration inequality for Gaussian random matrix (cont d)
MGMT 69: Topics in High-dimensional Data Analysis Falll 26 Lecture 8: Spectral clustering and Laplacian matrices Lecturer: Jiaming Xu Scribe: Hyun-Ju Oh and Taotao He, October 4, 26 Outline Concentration
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
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 informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationNumerical Solution of the Generalized Burgers-Huxley Equation by Exponential Time Differencing Scheme
International Journal of Bioinformatics and Biomedical Engineering Vol., No., 5, pp. 3-5 http://www.aiscience.org/journal/ijbbe Numerical Solution of the Generalized Burgers-Huxley Equation by Exponential
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem in a 2D Box Region Let us consider a model problem of parabolic
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationMatrix square root and interpolation spaces
Matrix square root and interpolation spaces Mario Arioli and Daniel Loghin m.arioli@rl.ac.uk STFC-Rutherford Appleton Laboratory, and University of Birmingham Sparse Days,CERFACS, Toulouse, 2008 p.1/30
More informationNumerical Schemes from the Perspective of Consensus
Numerical Schemes from the Perspective of Consensus Exploring Connections between Agreement Problems and PDEs Department of Electrical Engineering and Computer Sciences University of California, Berkeley
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More information