Runge Kutta Chebyshev methods for parabolic problems
|
|
- Erick Rice
- 5 years ago
- Views:
Transcription
1 Runge Kutta Chebyshev methods for parabolic problems Xueyu Zhu Division of Appied Mathematics, Brown University December 2, 2009 Xueyu Zhu 1/18
2 Outline Introdution Consistency condition Stability Properties Integration formula Numerical experiments and Performance issues Xueyu Zhu 2/18
3 Introduction Target Problem (semi-discrete, multi-space): U(t) = F (t, U(t)), 0 < t T, U(0) given F (t, U(t)) is originated from spatial discretization of parabolic equaiton, such as u t = u u t = ɛ u + f(u, x, t) Stiff system: explicit method : easy but time-step restrictive implicit method : larger time steps but expensive What s remedy? Xueyu Zhu 3/18
4 RKC formula S-stage RKC scheme: Y 0 = U n, Y 1 = Y 0 + ũ 1 τf 0, Y j = µ j Y j 1 + ν j Y j 2 + (1 µ j ν j )Y 0 + ũ j τf j 1 + γ j τf 0, (2 j s) U n+1 = Y s, n = 0, 1,.. which can be rewritten in the general form of Runge Kutta methods: j 1 Y j = U n + τ a jl F (t + c l τ, Y l ) l=0 where τ = t n t n 1 to specify a method, we need to specify s and c l, a jl Xueyu Zhu 4/18
5 Consistency condition Order 1: Suppose U(t) is sufficiently smooth, all Y j statisfy the following expansion: Y j = U(t n ) + c j τ U(t n ) + X j τ 2 U (2) (t n )) + O(τ 3 ), consistent of order 1 if c s = 1 based on the Taylor s expansion (Y j approximates U(t n + c j τ)). Subsititute this expression into RKC formula, we get X 0 = X 1 = 0 X j = µ j X j 1 + ν j X j 2 + µ j c j 1 Xueyu Zhu 5/18
6 Order 2: consistent of order requires X j = 1 2 c2 j, which leads to c 2 2 = 2ũ 2 c 1, c 2 3 = µ 3 c ũ 3 c 2, c 2 j = µ j c 2 j 1 + ν j c 2 j 2 + 2ũ j c j 1, 4 j s Xueyu Zhu 6/18
7 Stabiliy region and Stability function Scalar test problem: U = λu(t) For each stage: Y j = P j (z)y 0, z = λτ where P j is a polynomial of degree j At the final stage, U n+1 = P s (z)u n Put it back into the RKC formula, we get the recursion foruma for P j (z) without any assumption on P j (z): P 0 (z) = 1, P 1 (z) = 1 + µ 1 z P j (z) = (1 µ j ν j ) + γ j z + (µ j + µ j z)p j 1 (z) + ν j P j 2 (z) Xueyu Zhu 7/18
8 Stability Optimal stability per step: Recalling the consistency condition, each stability function approximates e c jz for z 0 as : P j (z) = 1 + c j z + X j z 2 + O(z 3 ) The stability region requires that S = {z C : P s 1} we want the largest boundary value on negative real axis: β(s) = max{ z : z 0, P s 1} Ease of implementation: need analytical form for integration coefficients Xueyu Zhu 8/18
9 Shifted Chebyshev polynomial Largest negative real stability boundary β(s) leads to Chebyshev polynomial of the first kind: T s (x) = cos(s arccos x), 1 x 1 for 1st order, the stability function : P s (z) = T s (1 + z ), β(s) z 0 s2 gives β(s) = 2s 2. Recall three-terms recursion formula for Chebyshev polynomial: P 0 (z) = 1, P 1 (z) = 1+ z s 2, P j(z) = 2(1+ z s 2 )P j 1 P j 2 (z), j 2 compare to the recursive formua, it is easy to get the analytical form of µ j, ν j, µ j, γ j : µ 1 = 1 s 2, µ j = 2, µ j = 2 s 2, ν j = 1, γ j = 0 Xueyu Zhu 9/18
10 First order forumla for 1st and 2nd order RKC, the stability functions fit in the general form: 1st order: P j (z) = a j + b j T j (w 0 + w 1 z), 0 j s a j = 0, b j = T 1 (w 0 ), w 0 = 1 + ɛ s 2, w 1 = T s(w 0 ) T s(w 0 ), 0 j s leads to β(s) = (2 4 3 ɛ)s2 as ɛ 0. ɛ is a damping factor choose ɛ = 0.05, then β(s) 1.9s 2. Xueyu Zhu 10/18
11 with this set of parameters, we can match them to three-terms of recursion formula to get the integration coeffcients: µ 1 = w 1 w 0, µ j = 2w 0 b j b j 1, ν j = b j b j 2 b j µ j = 2w 1, γ j = 0, 2 j s b j 1 c j = T s(w 0 ) T j (w 0) T s(w 0 ) T j (w 0 ) Xueyu Zhu 11/18
12 Second order formula 2nd order: a j = 1 b j T j (w 0 ), b j = T j (w 0) (T, (2 j s) (w 0 )) 2 w 0 = 1 + ɛ s 2, w 1 = T s(w 0 ) T s (w 0 ) a 0 = 1 b 0, a 1 = 1 b 1 w 0, b 0 = b 1 = b 2 for this set of parameters, we get β(s) (w 0 + 1)T s (w 0 ) 2 T s (w 0 ) 3 (s2 1)( ɛ) as ɛ 0. Xueyu Zhu 12/18
13 with this set of parameters, we can match them to three-terms of recursion formula to get the integration coeffcients: µ 1 = b 1 w 1, µ j = 2w 0 b j b j 1, ν j = b j b j 2 µ j = 2w 1 b j b j 1, γ j = (1 b j 1 T j 1 (w 0 )) µ j, 2 j s c 1 = c 2 T 2 (w 0), c j = T s(w 0 ) T j (w 0) T s (w 0 ) T j (w 0) Xueyu Zhu 13/18
14 Checking Consistency Expanding P s (z): P s (z) = 1 + z + z s2 10( 1 + s 2 ) z3... which is indeed 2nd order consistency Xueyu Zhu 14/18
15 Example 1: linear heat problem u t = u + f(x, y, z, t), 0 < x, y, z < 1, t > 0 take h = 0.025, 39 3 = equations Table: result of RKC2 and VODPK tol err steps F-evals CPU 1e e-02/0.99e-00 6/7 402/46 186/35 1e e-02/0.83e-01 15/16 729/ /122 1e e-03/0.10e-01 27/34 786/ /185 1e e-04/0.12e-02 57/ / /371 1e e-05/0.13e / / /770 1e e-06/0.19e / / /913 more accurate; approximate same amout of CPU time with relatively strict tolerence. Xueyu Zhu 15/18
16 Example 2: Combustion problem c t = c Dce δ T, LTt = T + αdce δ T, 0 < x, y, z < 1, t > 0 take N = 40, 10 6 equations Table: result of RKC2 and VODPK tol err steps F-evals CPU 1e e-00/0.87e-00 51/33 525/ /412 1e e-00/0.76e /91 781/ /957 1e e-01/0.12e / / /1702 1e e-02/0.12e / / /2376 lower accruacy is expected from the local instablity of the problem. RKC2 is more accurate and faster in general. Xueyu Zhu 16/18
17 Conclusion Desgined for modestly stiff, semi-discrete problems Construction of the extened stability region by shifted Chebyshev polynomial Qudratic increase of β(s) with number of stages Not good for eigenvalues of Jacobian that is far from negative real axis Xueyu Zhu 17/18
18 Reference J.G. Verwer, W.H. Hundsdorfer and B.P. Sommeijer, Convergence properties of the Runge-Kutta-Chebyshev method, Numer. Math. 57, , B.P. Sommeijer, L.F. Shampine and J.G. Verwer, RKC: An explicit solver for parabolic PDEs, J. Comp. Appl. Math. 88, , J.G. Verwer, B.P. Sommeijer and W. Hundsdorfer, RKC time-stepping for Advection-Diffusion-Reaction Problems, J. Comput. Phys. 201, 61-79, J.G. Verwer and B.P. Sommeijer, An Implicit-Explicit Runge-Kutta-Chebyshev Scheme for Diffusion-Reaction Equations, SIAM J. Scientific Computing 25, , J.G. Verwer, Explicit Runge-Kutta methods for parabolic partial differential equations, Applied Numerical Mathematics 22 (1996) Xueyu Zhu 18/18
Runge-Kutta-Chebyshev Projection Method
Runge-Kutta-Chebyshev Projection Method Zheming Zheng Linda Petzold Department of Mechanical Engineering, University of California Santa Barbara, Santa Barbara, CA 93106, USA June 8, 2006 This work was
More informationIRKC: an IMEX Solver for Stiff Diffusion-Reaction PDEs
IRKC: an IMEX Solver for Stiff Diffusion-Reaction PDEs L.F. Shampine Mathematics Department, Southern Methodist University Dallas, TX 75275-156, USA B.P. Sommeijer & J.G. Verwer Center for Mathematics
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 informationFlexible Stability Domains for Explicit Runge-Kutta Methods
Flexible Stability Domains for Explicit Runge-Kutta Methods Rolf Jeltsch and Manuel Torrilhon ETH Zurich, Seminar for Applied Mathematics, 8092 Zurich, Switzerland email: {jeltsch,matorril}@math.ethz.ch
More informationModule 4: Numerical Methods for ODE. Michael Bader. Winter 2007/2008
Outlines Module 4: for ODE Part I: Basic Part II: Advanced Lehrstuhl Informatik V Winter 2007/2008 Part I: Basic 1 Direction Fields 2 Euler s Method Outlines Part I: Basic Part II: Advanced 3 Discretized
More 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 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 informationFaster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs
Faster Kinetics: Accelerate Your Finite-Rate Combustion Simulation with GPUs Christopher P. Stone, Ph.D. Computational Science and Engineering, LLC Kyle Niemeyer, Ph.D. Oregon State University 2 Outline
More informationRKC time-stepping for advection diffusion reaction problems
Journal of Computational Physics 21 (24) 61 79 www.elsevier.com/locate/jcp RKC time-stepping for advection diffusion reaction problems J.G. Verwer *, B.P. Sommeijer, W. Hundsdorfer CWI, P.O. Box 9479,
More informationNumerische Mathematik
Numer. Math. (2007) 106:303 334 DOI 10.1007/s00211-006-0059-5 Numerische Mathematik Essentially optimal explicit Runge Kutta methods with application to hyperbolic parabolic equations Manuel Torrilhon
More informationMath 660 Lecture 4: FDM for evolutionary equations: ODE solvers
Math 660 Lecture 4: FDM for evolutionary equations: ODE solvers Consider the ODE u (t) = f(t, u(t)), u(0) = u 0, where u could be a vector valued function. Any ODE can be reduced to a first order system,
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 informationNumerical techniques for stiff and multiscale differential equations
2009 SIAM Annual Meeting July 6-10, 2009 Numerical techniques for stiff and multiscale differential equations Assyr Abdulle - Mathematics Section Swiss Federal Institute of Technology, Lausanne (EPFL)
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 informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationTime-adaptive methods for the incompressible Navier-Stokes equations
Time-adaptive methods for the incompressible Navier-Stokes equations Joachim Rang, Thorsten Grahs, Justin Wiegmann, 29.09.2016 Contents Introduction Diagonally implicit Runge Kutta methods Results with
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationHigh-Order Multi-Implicit Spectral Deferred Correction Methods for Problems of Reactive Flow
High-Order Multi-Implicit Spectral Deferred Correction Methods for Problems of Reactive Flow Anne Bourlioux Département de Mathématiques et Statistique, Université de Montréal, CP 6128 succ. Centre-Ville,
More informationarxiv: v2 [math.na] 22 Mar 2011
Application of Operator Splitting to Solve Reaction Diffusion Equations arxiv:1011.4810v2 [math.na] 22 Mar 2011 Abstract Tamás Ladics Szent István University, Ybl Miklós College of Building Approximate
More informationIntegrating-Factor-Based 2-Additive Runge Kutta methods for Advection-Reaction-Diffusion Equations
Integrating-Factor-Based 2-Additive Runge Kutta methods for Advection-Reaction-Diffusion Equations A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements
More 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 informationReview Higher Order methods Multistep methods Summary HIGHER ORDER METHODS. P.V. Johnson. School of Mathematics. Semester
HIGHER ORDER METHODS School of Mathematics Semester 1 2008 OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE 1 REVIEW 2 HIGHER ORDER METHODS 3 MULTISTEP METHODS 4 SUMMARY OUTLINE
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 informationarxiv: v1 [math.na] 15 Nov 2008
arxiv:0811.2481v1 [math.na] 15 Nov 2008 A Phase-Fitted Runge-Kutta-Nyström method for the Numerical Solution of Initial Value Problems with Oscillating Solutions Abstract D. F. Papadopoulos a, Z. A. Anastassi
More informationNumerical Oscillations and how to avoid them
Numerical Oscillations and how to avoid them Willem Hundsdorfer Talk for CWI Scientific Meeting, based on work with Anna Mozartova (CWI, RBS) & Marc Spijker (Leiden Univ.) For details: see thesis of A.
More 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 informationMartin Luther Universität Halle Wittenberg Institut für Mathematik
Martin Luther Universität Halle Wittenberg Institut für Mathematik Implicit peer methods for large stiff ODE systems St. Beck, R. Weiner, H. Podhaisky and B. A. Schmitt Report No. 07 (2010) Editors: Professors
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 informationSTABILITY FOR PARABOLIC SOLVERS
Review STABILITY FOR PARABOLIC SOLVERS School of Mathematics Semester 1 2008 OUTLINE Review 1 REVIEW 2 STABILITY: EXPLICIT METHOD Explicit Method as a Matrix Equation Growing Errors Stability Constraint
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 informationarxiv: v2 [math.na] 1 Apr 2016
High order schemes based on operator splitting and deferred corrections for stiff time dependent PDE s Max Duarte Matthew Emmett arxiv:1407.0195v2 [math.na] 1 Apr 2016 May 14, 2018 Abstract We consider
More informationOptimal Runge Kutta Stability Regions
Optimal Runge Kutta Stability Regions David I. Ketcheson Aron J. Ahmadia arxiv:11.335v1 [math.na] 1 Jan 1 March, 1 Abstract The stable step size for numerical integration of an initial value problem depends
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 informationModule 6: Implicit Runge-Kutta Methods Lecture 17: Derivation of Implicit Runge-Kutta Methods(Contd.) The Lecture Contains:
The Lecture Contains: We continue with the details about the derivation of the two stage implicit Runge- Kutta methods. A brief description of semi-explicit Runge-Kutta methods is also given. Finally,
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 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 informationStepsize Restrictions for Boundedness and Monotonicity of Multistep Methods
Stepsize Restrictions for Boundedness and Monotonicity of Multistep Methods W. Hundsdorfer, A. Mozartova, M.N. Spijker Abstract In this paper nonlinear monotonicity and boundedness properties are analyzed
More informationOn the mesh relaxation time in the moving mesh method
On the mesh relaxation time in the moving mesh method John Stockie http://www.math.sfu.ca/~stockie/ Department of Mathematics Simon Fraser University August 3, 2005 Adaptivity and Beyond Vancouver, August
More informationAdaptive Time Space Discretization for Combustion Problems
Konrad-Zuse-Zentrum fu r Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany JENS LANG 1,BODO ERDMANN,RAINER ROITZSCH Adaptive Time Space Discretization for Combustion Problems 1 Talk
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 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 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 informationFinite difference methods for the diffusion equation
Finite difference methods for the diffusion equation D150, Tillämpade numeriska metoder II Olof Runborg May 0, 003 These notes summarize a part of the material in Chapter 13 of Iserles. They are based
More informationNumerical Methods in Physics and Astrophysics
Kostas Kokkotas 2 October 20, 2014 2 http://www.tat.physik.uni-tuebingen.de/ kokkotas Kostas Kokkotas 3 TOPICS 1. Solving nonlinear equations 2. Solving linear systems of equations 3. Interpolation, approximation
More informationUNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH 2 IMPLICIT-EXPLICIT NUMERICAL METHOD
UNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH IMPLICIT-EXPLICIT NUMERICAL METHOD ANDREW JORGENSON Abstract. Systems of nonlinear partial differential equations modeling turbulent fluid flow
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 informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Partial differential equations arise in a number of physical problems, such as fluid flow, heat transfer, solid mechanics and biological processes. These
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 informationThe RAMSES code and related techniques I. Hydro solvers
The RAMSES code and related techniques I. Hydro solvers Outline - The Euler equations - Systems of conservation laws - The Riemann problem - The Godunov Method - Riemann solvers - 2D Godunov schemes -
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 informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
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 informationA New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws
A New Fourth-Order Non-Oscillatory Central Scheme For Hyperbolic Conservation Laws A. A. I. Peer a,, A. Gopaul a, M. Z. Dauhoo a, M. Bhuruth a, a Department of Mathematics, University of Mauritius, Reduit,
More informationarxiv: v2 [math.oc] 4 Oct 2018
Explicit Stabilised Gradient Descent for Faster Strongly Convex Optimisation arxiv:185.7199v2 [math.oc] 4 Oct 218 Armin Eftekhari Alan Turing Institute, London, UK aeftekhari@turing.ac.uk Gilles Vilmart
More informationNumerical Analysis. A Comprehensive Introduction. H. R. Schwarz University of Zürich Switzerland. with a contribution by
Numerical Analysis A Comprehensive Introduction H. R. Schwarz University of Zürich Switzerland with a contribution by J. Waldvogel Swiss Federal Institute of Technology, Zürich JOHN WILEY & SONS Chichester
More informationDifferential Equations
Pysics-based simulation xi Differential Equations xi+1 xi xi+1 xi + x x Pysics-based simulation xi Wat is a differential equation? Differential equations describe te relation between an unknown function
More informationSelected HW Solutions
Selected HW Solutions HW1 1 & See web page notes Derivative Approximations. For example: df f i+1 f i 1 = dx h i 1 f i + hf i + h h f i + h3 6 f i + f i + h 6 f i + 3 a realmax 17 1.7014 10 38 b realmin
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationNumerical Treatment of IVP ODEs... in a Nutshell
Numerical Treatment of IVP ODEs... in a Nutshell AT Patera, M Yano October 31, 2014 Draft V1.3 MIT 2014. From Math, Numerics, & Programming for Mechanical Engineers... in a Nutshell by AT Patera and M
More informationQuarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations
MATEMATIKA, 2011, Volume 27, Number 2, 199 208 c Department of Mathematical Sciences, UTM Quarter-Sweep Gauss-Seidel Method for Solving First Order Linear Fredholm Integro-differential Equations 1 E. Aruchunan
More informationSaul Abarbanel; Half a century of scientific work. Bertil Gustafsson, Uppsala University
Saul Abarbanel; Half a century of scientific work Bertil Gustafsson, Uppsala University Grew up in Tel Aviv Served in Israeli Army during the War of Independence 1948 1950 MIT 1952 1959 Ph.D 1959, Theoretical
More informationQuadratic SDIRK pair for treating chemical reaction problems.
Quadratic SDIRK pair for treating chemical reaction problems. Ch. Tsitouras TEI of Chalkis, Dept. of Applied Sciences, GR 34400 Psahna, Greece. I. Th. Famelis TEI of Athens, Dept. of Mathematics, GR 12210
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
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 informationChapter 10 Exercises
Chapter 10 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rl/fdmbook Exercise 10.1 (One-sided and
More informationRicharson Extrapolation for Runge-Kutta Methods
Richarson Extrapolation for Runge-Kutta Methods Zahari Zlatevᵃ, Ivan Dimovᵇ and Krassimir Georgievᵇ ᵃ Department of Environmental Science, Aarhus University, Frederiksborgvej 399, P. O. 358, 4000 Roskilde,
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 informationA self adjusting multirate algorithm based on the TR-BDF2 method
A self adjusting multirate algorithm based on the TR-BDF2 method arxiv:8.98v [math.na] 27 Jan 28 Luca Bonaventura (), Francesco Casella (2) Ludovica Delpopolo (), Akshay Ranade (3) September, 28 () MOX
More informationButcher tableau Can summarize an s + 1 stage Runge Kutta method using a triangular grid of coefficients
AM 205: lecture 13 Last time: ODE convergence and stability, Runge Kutta methods Today: the Butcher tableau, multi-step methods, boundary value problems Butcher tableau Can summarize an s + 1 stage Runge
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 informationNewton s Method and Efficient, Robust Variants
Newton s Method and Efficient, Robust Variants Philipp Birken University of Kassel (SFB/TRR 30) Soon: University of Lund October 7th 2013 Efficient solution of large systems of non-linear PDEs in science
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 informationPart 1. The diffusion equation
Differential Equations FMNN10 Graded Project #3 c G Söderlind 2016 2017 Published 2017-11-27. Instruction in computer lab 2017-11-30/2017-12-06/07. Project due date: Monday 2017-12-11 at 12:00:00. Goals.
More informationNumerical Smoothness and PDE Error Analysis
Numerical Smoothness and PDE Error Analysis GLSIAM 2013 Central Michigan University Tong Sun Bowling Green State University April 20, 2013 Table of contents Why consider Numerical Smoothness? What is Numerical
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 informationFour Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential Equations
International journal of scientific and technical research in engineering (IJSTRE) www.ijstre.com Volume Issue ǁ July 206. Four Point Gauss Quadrature Runge Kuta Method Of Order 8 For Ordinary Differential
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 informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationTime stepping methods
Time stepping methods ATHENS course: Introduction into Finite Elements Delft Institute of Applied Mathematics, TU Delft Matthias Möller (m.moller@tudelft.nl) 19 November 2014 M. Möller (DIAM@TUDelft) Time
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationUNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH 2 NUMERICAL METHOD
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING, SERIES B Volume 5, Number 3, Pages 7 87 c 04 Institute for Scientific Computing and Information UNCONDITIONAL STABILITY OF A CRANK-NICOLSON ADAMS-BASHFORTH
More informationCommunications in Applied Mathematics and Computational Science
Communications in Applied Mathematics and Computational Science OPTIMAL STABILITY POLYNOMIALS FOR NUMERICAL INTEGRATION OF INITIAL VALUE PROBLEMS DAVID I. KETCHESON AND ARON J. AHMADIA vol. 7 no. 1 msp
More informationAn Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations
An Accurate Fourier-Spectral Solver for Variable Coefficient Elliptic Equations Moshe Israeli Computer Science Department, Technion-Israel Institute of Technology, Technion city, Haifa 32000, ISRAEL Alexander
More informationAdaptive Integration of Stiff ODE
INFORMATICA, 2000, Vol. 11, No. 4, 371 380 371 2000 Institute of Mathematics and Informatics, Vilnius Adaptive Integration of Stiff ODE Raimondas ČIEGIS, Olga SUBOČ, Vilnius Gediminas Technical University
More informationOn the Diagonal Approximation of Full Matrices
On the Diagonal Approximation of Full Matrices Walter M. Lioen CWI P.O. Box 94079, 090 GB Amsterdam, The Netherlands ABSTRACT In this paper the construction of diagonal matrices, in some
More informationHigher Order Taylor Methods
Higher Order Taylor Methods Marcelo Julio Alvisio & Lisa Marie Danz May 6, 2007 Introduction Differential equations are one of the building blocks in science or engineering. Scientists aim to obtain numerical
More informationExponentially Fitted Error Correction Methods for Solving Initial Value Problems
KYUNGPOOK Math. J. 52(2012), 167-177 http://dx.doi.org/10.5666/kmj.2012.52.2.167 Exponentially Fitted Error Correction Methods for Solving Initial Value Problems Sangdong Kim and Philsu Kim Department
More informationLinearly implicit Runge Kutta methods for advection reaction diffusion equations
Applied Numerical Mathematics 37 (2001) 535 549 Linearly implicit Runge Kutta methods for advection reaction diffusion equations M.P. Calvo, J. de Frutos, J. Novo Departamento de Matemática Aplicada y
More informationarxiv: v3 [math.na] 12 Jul 2012
Optimal stability polynomials for numerical integration of initial value problems David I. Ketcheson Aron J. Ahmadia arxiv:11.335v3 [math.na] 1 Jul 1 July 13, 1 Abstract We consider the problem of finding
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 06-15 A Note on the Numerical Simulation of Kleijn s Benchmark Problem S. van Veldhuizen, C. Vuik, C.R. Kleijn ISSN 1389-6520 Reports of the Department of Applied
More informationODEs. PHY 688: Numerical Methods for (Astro)Physics
ODEs ODEs ODEs arise in many physics problems Classifications: As with the other topics, there are a large number of different methods Initial value problems Boundary value problems Eigenvalue problems
More information4 Stability analysis of finite-difference methods for ODEs
MATH 337, by T. Lakoba, University of Vermont 36 4 Stability analysis of finite-difference methods for ODEs 4.1 Consistency, stability, and convergence of a numerical method; Main Theorem In this Lecture
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 informationBlanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS
Opuscula Matematica Vol. 26 No. 3 26 Blanca Bujanda, Juan Carlos Jorge NEW EFFICIENT TIME INTEGRATORS FOR NON-LINEAR PARABOLIC PROBLEMS Abstract. In tis work a new numerical metod is constructed for time-integrating
More informationarxiv: v2 [math-ph] 29 Nov 2010
An Efficient NRxx Method for Boltzmann-BGK Equation Zhenning Cai, Ruo Li and Yanli Wang arxiv:1011.5789v2 [math-ph] 29 Nov 2010 November 6, 2018 Abstract In [3], we proposed a numerical regularized moment
More informationApplication of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms
Future Generation Computer Systems () 65 7 Application of the Kurganov Levy semi-discrete numerical scheme to hyperbolic problems with nonlinear source terms R. Naidoo a,b, S. Baboolal b, a Department
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationMath 273, Final Exam Solutions
Math 273, Final Exam Solutions 1. Find the solution of the differential equation y = y +e x that satisfies the condition y(x) 0 as x +. SOLUTION: y = y H + y P where y H = ce x is a solution of the homogeneous
More information