7 Hyperbolic Differential Equations
|
|
- Reginald Grant
- 5 years ago
- Views:
Transcription
1 Numerical Analysis of Differential Equations Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena. Important application areas include Acoustics Electromagnetics Seismics Optics Fluid Mechanics Modeling the propagation of small perturbations of an equilibrium typically leads to linear equations. 7 Hyperbolic Differential Equations TU Bergakademie Freiberg, SS 2012
2 Numerical Analysis of Differential Equations The Linear Advection Equation We consider the simplest hyperbolic PDE, the linear advection equation, sometimes called the one-way wave equation. The Cauchy problem for this PDE reads with a constant c. u t + cu x = 0, x R, t > 0, (7.1a) u(x, 0) = u 0 (x), (7.1b) The solution of (7.1) poses little difficulty, as the exact solution is easily verified to be u(x, t) = u 0 (x ct). Nonetheless, this simple problem serves to illustrate typical difficulties encountered in the numerical solution of hyperbolic PDEs. 7.1 The Linear Advection Equation TU Bergakademie Freiberg, SS 2012
3 Numerical Analysis of Differential Equations 245 We again introduce a space-time grid {(x j, t n ) : x j = j x, j Z; t n = n t, n N 0 }. Natural explicit scheme: central difference approximation of u x combined with explicit Euler method in time, results in or, in explicit form, U n+1 j t U n j + c U n j+1 U n j 1 2 x = 0 U n+1 j = U n j c t 2 x (U n j+1 U n j 1). (7.2) As we will later see, this is a useless method. However, replacing Uj n by the average 1 2 (U j 1 n +U j+1 n ) leads to the eminently useful Lax-Friedrichs method U n+1 j = 1 2 (U n j 1 + U n j+1) c t 2 x (U n j+1 U n j 1). (7.3) 7.1 The Linear Advection Equation TU Bergakademie Freiberg, SS 2012
4 Numerical Analysis of Differential Equations 246 We will see that (7.3) is Lax-Richtmyer stable if c t x 1. In particular, this permits choosing a timestep of t = O( x). 7.1 The Linear Advection Equation TU Bergakademie Freiberg, SS 2012
5 Numerical Analysis of Differential Equations Method of Lines Discretisation Formulating an IBVP for (7.1a) on the spatial domain x (0, 1) requires a boundary condition on one of the two endpoints. When c > 0, a boundary condition at x = 0 leads to a well-posed problem, e.g. u(0, t) = g 0 (t), t 0. In this case the left boundary is called the inflow boundary, and the right the outflow boundary. Their roles are reversed when c < 0. For the two schemes considered so far, some numerical boundary condition must be chosen to determine the values at the outflow boundary, which complicates the stability analysis. We therefore simplify matters by choosing periodic boundary conditions and obtain the IBVP... u(0, t) = u(1, t), t 0, 7.2 Method of Lines Discretisation TU Bergakademie Freiberg, SS 2012
6 Numerical Analysis of Differential Equations 248 u t + cu x = 0, x (0, 1), t > 0, (7.4a) u(0, t) = u(1, t), t 0, (7.4b) u(x, 0) = u 0 (x). The vector of approximations at the spatial gridpoints (semidiscretization!) U 1 (t) U (t) =., U j(t) u(x j, t), U J+1 (t) must then solve the system of ODEs U 1(t) = c [ U2 (t) U J+1 (t) ], 2 x U j(t) = c [ Uj+1 (t) U j 1 (t) ], 2 j J, 2 x U J+1(t) = c [ U1 (t) U J (t) ]. 2 x (7.4c) 7.2 Method of Lines Discretisation TU Bergakademie Freiberg, SS 2012
7 Numerical Analysis of Differential Equations 249 In matrix notation: U (t) = AU (t), U (0) = U 0, where [U 0 ] j = u 0 (x j ), and A = c 2 x R (J+1) (J+1). (7.5) The matrix A is skew symmetric (A = A), and its eigenvalues are λ j = ic x with associated eigenvectors sin(2πj x), j = 1, 2,..., J + 1, [v j ] k = e 2πikj x, j = 1, 2,..., J Method of Lines Discretisation TU Bergakademie Freiberg, SS 2012
8 Numerical Analysis of Differential Equations 250 The eigenvalues of A are thus contained in the interval [ ic x, ic x] on the imaginary axis. For absolute stability in the method of lines, a time-stepping scheme is required whose region of absolute stability contains this interval. For the explicit Euler method (7.2) this is not the case for any value of t/ x. This method is therefore not absolutely stable for any value of this grid ratio. A different choice, however, satisfies the weaker variant (6.11) of Lax- Richtmyer stability: for the eigenvalues of the matrix B( t) = I + ta we have 1 + tλ j 2 = 1 + ( ) 2 c t sin 2 (2πj x) 1 + x ( ) 2 c t, x so that choosing t = x 2 results in 1 + tλ j c 2 t and thus to the bound B n = (I + ta) n (1 + c 2 t) n/2 e c2 T/2, n t T. 7.2 Method of Lines Discretisation TU Bergakademie Freiberg, SS 2012
9 Numerical Analysis of Differential Equations The Leapfrog Method If, in (7.1a) or (7.4a), using the same central difference formlula for u x, one uses the midpoint rule y n+1 = y n t f(t n, y n ) in place of the explicit Euler method, one obtains the so-called leapfrog scheme U n+1 j = U n 1 j c t x (U j+1 n Uj 1). n (7.6) Since the region of absolute stability of the midpoint rule is the interval {iy : 1 < y < 1} on the imaginary axis, the leapfrog scheme is absolutely stable if c t x < The Leapfrog Method TU Bergakademie Freiberg, SS 2012
10 Numerical Analysis of Differential Equations The Lax-Friedrichs Scheme The replacement 1 2 (U n j 1 + U n j+1 ) = U n j (U n j 1 2U n j + U n j+1 ) modifies the Lax-Friedrichs scheme (7.3) to U n+1 j = U n j c t 2 x (U n j+1 U n j 1) (U n j 1 2U n j and, upon rearranging the terms to + U n j+1), U n+1 j t U n j + c U n j+1 U n j 1 2 x = x2 2 t U n j 1 2U n j + U n j+1 x 2. (7.7) Consistency analysis reveals this scheme to be consistent with the PDE u t + cu x = 0. However, it looks more like a discretization of the advectiondiffusion equation u t + cu x = ɛu xx, ɛ = x2 2 t. 7.4 The Lax-Friedrichs Scheme TU Bergakademie Freiberg, SS 2012
11 Numerical Analysis of Differential Equations 253 The scheme (7.7) thus consists of the explicit Euler method applied to A ɛ = c 2 x U (t) = A ɛ U (t), with (7.8) ɛ x The matrix (7.8) is obtained from the matrix (7.5) by adding a small multiple of a symmetric, negative definite difference operator, causing its eigelvalues to move away from the imaginary axis into the left half plane. Due to the periodic boundary conditions both matrices are circulant and therefore have the same (orthogonal) eigenvectors The eigenvalues are λ ɛ j = ic x sin(2πj x) 2ɛ (1 cos(2πj x), j = 1, 2,..., J. (7.9) x2 7.4 The Lax-Friedrichs Scheme TU Bergakademie Freiberg, SS 2012
12 Numerical Analysis of Differential Equations Eigenvalues of ta ɛ for c = 1, x = 1/50 t = 0.8 x for different values of ɛ expl. Euler (ɛ = 0) ɛ = ɛ = Lax-Friedrichs (ɛ = ) 7.4 The Lax-Friedrichs Scheme TU Bergakademie Freiberg, SS 2012
13 Numerical Analysis of Differential Equations The Lax-Wendroff Scheme Possibilities for increasing consistency order w.r.t. t: Use 2nd order time stepping method such as midpoint rule, leapfrog. Disadvantages: involves three consecutive time levels (large memory requirements, especially in 3D), BC more difficult to implement. Leapfrog nondissipative (neutral stability), problems with variable coefficients/nonlinear terms Trapezoidal rule, but this is implicit. Two-stage Runge-Kutta method Taylor series approach applied to u (t) = Au(t), taking account of u (t) = Au (t) = A 2 u(t): u n+1 = u n + tau n + t2 2 A2 u n. 7.5 The Lax-Wendroff Scheme TU Bergakademie Freiberg, SS 2012
14 Numerical Analysis of Differential Equations 256 If A is the matrix from (7.5), this method written out is U n+1 j = U n j c t 2 x (U n j+1 U n j 1) + c2 t 2 8 x 2 (U n j 2 2U n j Bad: 5-point stencil in x, difficult near boundary. + U n j+2). Remedy: note that last term represents approximation to 1 2 c2 t 2 u xx. Replace this term by usual 3-point formula for u xx. This results in the Lax- Wendroff scheme U n+1 j = U n j c t 2 x (U n j+1 U n j 1) + c2 t 2 2 x 2 (U n j 1 2U n j + U n j+1). (7.10) 7.5 The Lax-Wendroff Scheme TU Bergakademie Freiberg, SS 2012
15 Numerical Analysis of Differential Equations Asymptotic Stability Note: The Lax-Wendroff scheme corresponds to the explicit Euler method applied to U (t) = A ɛ U (t) with A ɛ as in (7.8), except that now ɛ = c 2 t/2 in place of ɛ = x 2 /(2 t) as for the Lax-Friedrichs method. Just as in (7.9), we investigate the eigenvalues of ta ɛ : tλ ɛ j = i c t x sin(jπ x) + ( ) 2 c t (cos(jπ x 1), j = 1, 2,..., J. x These are located on the boundary of an ellipse with midpoint (c t/ x) 2 and semi-axes of length (c c/ x) 2 and c t/ x. This ellipse is contained in the region of absolute stability of the Euler scheme if c t x The Lax-Wendroff Scheme TU Bergakademie Freiberg, SS 2012
16 Numerical Analysis of Differential Equations 258 Based on the same parameters as in the previous section this figure shows the eigenvalues of the Lax-Wendroff method Lax-Wendroff (ɛ = 0.005) 7.5 The Lax-Wendroff Scheme TU Bergakademie Freiberg, SS 2012
17 Numerical Analysis of Differential Equations Upwind Schemes Besides the central difference formulas we have used up to now for approximating u x in the linear advection equation (7.1a), one could also consider using the one-sided difference formulas u x (x j, t) U j U j 1 x or u x (x j, t) U j+1 U j x Combined with a forward difference in time these lead to the schemes U n+1 j = U n j c t x (U j U j 1 ), U n+1 j = U n j c t x (U j+1 U j ). Both are consistent with order one in x and t.. (7.11a) (7.11b) 7.6 Upwind Schemes TU Bergakademie Freiberg, SS 2012
18 Numerical Analysis of Differential Equations 260 Despite their low order of accuracy these two schemes mimic a qualitative property of the solution of (7.1a), namely the propagation of information at speed c from left to right (c > 0) or right to left (c < 0). More precisely, for the exact solution we have u(x j, t + t) = u 0 (x j c(t + t)) = u(x j c t, t), i.e., the value of the solution in the next time step is determined by data which, at the current time step, lies to the left of x j (c > 0) or to the right of x j (c < 0). This suggests using (7.11a) when c > 0 and (7.11b) when c < Upwind Schemes TU Bergakademie Freiberg, SS 2012
19 Numerical Analysis of Differential Equations 261 To determine the stability properties of these schemes, note that (7.11a) may also be rearranged as U n+1 j = U n j c t 2 x (U n j+1 U n j 1) + c t 2 x (U n j+1 2U n j + U n j 1) (7.12) which is a scheme of the form (7.8) with ɛ = c x/2. Already established: such a scheme is stable if c t x 1 and 2 < 2ɛ t x 2 < Upwind Schemes TU Bergakademie Freiberg, SS 2012
20 Numerical Analysis of Differential Equations Numerical Example We consider the IVP (7.11a) with c = 1 and initial data u 0 (x) = e 20(x 2)2 + e (x 5)2 in the time range t [0, 17] and compare the exact solution with the approximations provided by the upwind, Lax-Wendroff and leapfrog schemes at the end of the time interval. The initial data consists of two Gaussian peaks, one of which is less well resolved by the spatial grid than the other. The spatial mesh size is x = 0.05, the timestep t = 0.8 x; we use homogeneous Dirichlet conditions as both (true and numerical) BC. The solid black line shows the extact solution, the red dots their numerical approximations at the grid points. 7.7 Numerical Example TU Bergakademie Freiberg, SS 2012
21 Numerical Analysis of Differential Equations Exakte Lösung t=0 1.5 Upwind, t= Numerical Example TU Bergakademie Freiberg, SS 2012
22 Numerical Analysis of Differential Equations Lax Wendroff, t= Leapfrog, t= Numerical Example TU Bergakademie Freiberg, SS 2012
23 Numerical Analysis of Differential Equations 265 Observations: The upwind scheme is strongly diffusive and approximates both peaks less well than the remaining methods. This is partly due to its low order of accuracy and can be improved by refining the mesh. Lax-Wendroff and leapfrog display a strong dispersion, i.e., different wave components in the numerical approximation are propagated with different speeds. 7.7 Numerical Example TU Bergakademie Freiberg, SS 2012
24 Numerical Analysis of Differential Equations 266 Goals of Section 7. You should be able to state the linear advection equation and discuss its solution properties; be familiar with the construction principles of the leapfrog, Lax-Friedrichs, Lax-Wendroff and upwind schemes; know the difference between absolute stability and Lax-Richtmyer stability; be able to carry out a von Neumann stability analysis for numerical schemes approximating the linear advection equation; be able to implement simple difference scheme for the linear advection equation. 7.7 Numerical Example TU Bergakademie Freiberg, SS 2012
Numerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
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 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 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 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 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 informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
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 informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationIntroduction to PDEs and Numerical Methods: Exam 1
Prof Dr Thomas Sonar, Institute of Analysis Winter Semester 2003/4 17122003 Introduction to PDEs and Numerical Methods: Exam 1 To obtain full points explain your solutions thoroughly and self-consistently
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationLecture 17: Initial value problems
Lecture 17: Initial value problems Let s start with initial value problems, and consider numerical solution to the simplest PDE we can think of u/ t + c u/ x = 0 (with u a scalar) for which the solution
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 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 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 informationApplied Mathematics 205. Unit III: Numerical Calculus. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit III: Numerical Calculus Lecturer: Dr. David Knezevic Unit III: Numerical Calculus Chapter III.3: Boundary Value Problems and PDEs 2 / 96 ODE Boundary Value Problems 3 / 96
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationCS205B / CME306 Homework 3. R n+1 = R n + tω R n. (b) Show that the updated rotation matrix computed from this update is not orthogonal.
CS205B / CME306 Homework 3 Rotation Matrices 1. The ODE that describes rigid body evolution is given by R = ω R. (a) Write down the forward Euler update for this equation. R n+1 = R n + tω R n (b) Show
More informationNumerical Methods for PDEs
Numerical Methods for PDEs Problems 1. Numerical Differentiation. Find the best approximation to the second drivative d 2 f(x)/dx 2 at x = x you can of a function f(x) using (a) the Taylor series approach
More 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 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 informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationPartial differential equations
Partial differential equations Many problems in science involve the evolution of quantities not only in time but also in space (this is the most common situation)! We will call partial differential equation
More informationDevelopment and stability analysis of the inverse Lax-Wendroff boundary. treatment for central compact schemes 1
Development and stability analysis of the inverse Lax-Wendroff boundary treatment for central compact schemes François Vilar 2 and Chi-Wang Shu 3 Division of Applied Mathematics, Brown University, Providence,
More informationINTRODUCTION TO PDEs
INTRODUCTION TO PDEs In this course we are interested in the numerical approximation of PDEs using finite difference methods (FDM). We will use some simple prototype boundary value problems (BVP) and initial
More informationNumerical Methods for Engineers and Scientists
Numerical Methods for Engineers and Scientists Second Edition Revised and Expanded Joe D. Hoffman Department of Mechanical Engineering Purdue University West Lafayette, Indiana m MARCEL D E К К E R MARCEL
More information12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis
ATH 337, by T. Lakoba, University of Vermont 113 12 The Heat equation in one spatial dimension: Simple explicit method and Stability analysis 12.1 Formulation of the IBVP and the minimax property of its
More informationAM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs
AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationPartial Differential Equations II
Partial Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Partial Differential Equations II 1 / 28 Almost Done! Homework
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 informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationNumerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer C2 b 2
Numerical Methods for Conservation Laws WPI, January 2006 C. Ringhofer ringhofer@asu.edu, C2 b 2 2 h2 x u http://math.la.asu.edu/ chris Last update: Jan 24, 2006 1 LITERATURE 1. Numerical Methods for Conservation
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationLecture Introduction
Lecture 1 1.1 Introduction The theory of Partial Differential Equations (PDEs) is central to mathematics, both pure and applied. The main difference between the theory of PDEs and the theory of Ordinary
More informationNumerical methods Revised March 2001
Revised March 00 By R. W. Riddaway (revised by M. Hortal) Table of contents. Some introductory ideas. Introduction. Classification of PDE's.3 Existence and uniqueness.4 Discretization.5 Convergence, consistency
More informationClassification of partial differential equations and their solution characteristics
9 TH INDO GERMAN WINTER ACADEMY 2010 Classification of partial differential equations and their solution characteristics By Ankita Bhutani IIT Roorkee Tutors: Prof. V. Buwa Prof. S. V. R. Rao Prof. U.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS NUMERICAL FLUID MECHANICS FALL 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 02139 2.29 NUMERICAL FLUID MECHANICS FALL 2011 QUIZ 2 The goals of this quiz 2 are to: (i) ask some general
More informationFinite Differences: Consistency, Stability and Convergence
Finite Differences: Consistency, Stability and Convergence Varun Shankar March, 06 Introduction Now that we have tackled our first space-time PDE, we will take a quick detour from presenting new FD methods,
More informationA Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations
A Bound-Preserving Fourth Order Compact Finite Difference Scheme for Scalar Convection Diffusion Equations Hao Li Math Dept, Purdue Univeristy Ocean University of China, December, 2017 Joint work with
More informationDo not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE712B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.
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 informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More information1 Upwind scheme for advection equation with variable. 2 Modified equations: numerical dissipation and dispersion
1 Upwind sceme for advection equation wit variable coefficient Consider te equation u t + a(x)u x Applying te upwind sceme, we ave u n 1 = a (un u n 1), a 0 u n 1 = a (un +1 u n ) a < 0. CFL condition
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationENO and WENO schemes. Further topics and time Integration
ENO and WENO schemes. Further topics and time Integration Tefa Kaisara CASA Seminar 29 November, 2006 Outline 1 Short review ENO/WENO 2 Further topics Subcell resolution Other building blocks 3 Time Integration
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationNumerical Solution of partial differential equations
G. D. SMITH Brunei University Numerical Solution of partial differential equations FINITE DIFFERENCE METHODS THIRD EDITION CLARENDON PRESS OXFORD Contents NOTATION 1. INTRODUCTION AND FINITE-DIFFERENCE
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More informationNumerical Solutions to PDE s
Introduction Numerical Solutions to PDE s Mathematical Modelling Week 5 Kurt Bryan Let s start by recalling a simple numerical scheme for solving ODE s. Suppose we have an ODE u (t) = f(t, u(t)) for some
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationNotes: Outline. Shock formation. Notes: Notes: Shocks in traffic flow
Outline Scalar nonlinear conservation laws Traffic flow Shocks and rarefaction waves Burgers equation Rankine-Hugoniot conditions Importance of conservation form Weak solutions Reading: Chapter, 2 R.J.
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 informationProblem Set 4 Issued: Wednesday, March 18, 2015 Due: Wednesday, April 8, 2015
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0139.9 NUMERICAL FLUID MECHANICS SPRING 015 Problem Set 4 Issued: Wednesday, March 18, 015 Due: Wednesday,
More informationNumerical Algorithms for Visual Computing II 2010/11 Example Solutions for Assignment 6
Numerical Algorithms for Visual Computing II 00/ Example Solutions for Assignment 6 Problem (Matrix Stability Infusion). The matrix A of the arising matrix notation U n+ = AU n takes the following form,
More informationOutline. 1 Partial Differential Equations. 2 Numerical Methods for PDEs. 3 Sparse Linear Systems
Outline Partial Differential Equations Numerical Methods for PDEs Sparse Linear Systems 1 Partial Differential Equations 2 Numerical Methods for PDEs 3 Sparse Linear Systems Michael T Heath Scientific
More informationHIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 HIGH ORDER LAX-WENDROFF-TYPE SCHEMES FOR LINEAR WAVE PROPAGATION
More informationRelevant self-assessment exercises: [LIST SELF-ASSESSMENT EXERCISES HERE]
Chapter 6 Finite Volume Methods In the previous chapter we have discussed finite difference methods for the discretization of PDEs. In developing finite difference methods we started from the differential
More informationThe Lattice Boltzmann method for hyperbolic systems. Benjamin Graille. October 19, 2016
The Lattice Boltzmann method for hyperbolic systems Benjamin Graille October 19, 2016 Framework The Lattice Boltzmann method 1 Description of the lattice Boltzmann method Link with the kinetic theory Classical
More informationREVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS
REVIEW OF SUMMATION-BY-PARTS SCHEMES FOR INITIAL-BOUNDARY-VALUE PROBLEMS MAGNUS SVÄRD AND JAN NORDSTRÖM Abstract. High-order finite difference methods are efficient, easy to program, scale well in multiple
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationDissipation and Dispersion
Consider the problem with periodic boundary conditions Dissipation and Dispersion u t = au x 0 < x < 1, t > 0 u 0 = sin 40 πx u(0, t) = u(1, t) t > 0 If a > 0 then the wave is moving to the left and if
More informationECE539 - Advanced Theory of Semiconductors and Semiconductor Devices. Numerical Methods and Simulation / Umberto Ravaioli
ECE539 - Advanced Theory of Semiconductors and Semiconductor Devices 1 General concepts Numerical Methods and Simulation / Umberto Ravaioli Introduction to the Numerical Solution of Partial Differential
More informationNumerical Integration of Linear and Nonlinear Wave Equations
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Dissertations, Theses, and Student Research Papers in Mathematics Mathematics, Department of 12-2004 Numerical Integration
More informationNumerical Methods for Partial Differential Equations CAAM 452. Spring 2005
Numerical Methods for Partial Differential Equations Instructor: Tim Warburton Class Location: Duncan Hall 1046 Class Time: 9:5am to 10:40am Office Hours: 10:45am to noon in DH 301 CAAM 45 Spring 005 Homeworks
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Modified Equation of a Difference Scheme What is a Modified Equation of a Difference
More information3.3. Phase and Amplitude Errors of 1-D Advection Equation
3.3. Phase and Amplitude Errors of 1-D Advection Equation Reading: Duran section 2.4.2. Tannehill et al section 4.1.2. The following example F.D. solutions of a 1D advection equation show errors in both
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationDepartment of Mathematics California State University, Los Angeles Master s Degree Comprehensive Examination in. NUMERICAL ANALYSIS Spring 2015
Department of Mathematics California State University, Los Angeles Master s Degree Comprehensive Examination in NUMERICAL ANALYSIS Spring 2015 Instructions: Do exactly two problems from Part A AND two
More informationVISCOUS FLUX LIMITERS
VISCOUS FLUX LIMITERS E. F. Toro Department of Aerospace Science College of Aeronautics Cranfield Institute of Technology Cranfield, Beds MK43 OAL England. Abstract We present Numerical Viscosity Functions,
More informationDiagonal-norm upwind SBP operators
Diagonal-norm upwind SBP operators Ken Mattsson June 8, 16 Abstract High-order accurate first derivative finite difference operators are derived that naturally introduce artificial dissipation. The boundary
More informationCharacteristics for IBVP. Notes: Notes: Periodic boundary conditions. Boundary conditions. Notes: In x t plane for the case u > 0: Solution:
AMath 574 January 3, 20 Today: Boundary conditions Multi-dimensional Wednesday and Friday: More multi-dimensional Reading: Chapters 8, 9, 20 R.J. LeVeque, University of Washington AMath 574, January 3,
More informationThe Initial Value Problem for Ordinary Differential Equations
Chapter 5 The Initial Value Problem for Ordinary Differential Equations In this chapter we begin a study of time-dependent differential equations, beginning with the initial value problem (IVP) for a time-dependent
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 9
Spring 2015 Lecture 9 REVIEW Lecture 8: Direct Methods for solving (linear) algebraic equations Gauss Elimination LU decomposition/factorization Error Analysis for Linear Systems and Condition Numbers
More informationSurprising Computations
.... Surprising Computations Uri Ascher Department of Computer Science University of British Columbia ascher@cs.ubc.ca www.cs.ubc.ca/ ascher/ Uri Ascher (UBC) Surprising Computations Fall 2012 1 / 67 Motivation.
More informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationNumerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement
Numerical Solutions for Hyperbolic Systems of Conservation Laws: from Godunov Method to Adaptive Mesh Refinement Romain Teyssier CEA Saclay Romain Teyssier 1 Outline - Euler equations, MHD, waves, hyperbolic
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
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 informationWhat we ll do: Lecture 21. Ordinary Differential Equations (ODEs) Differential Equations. Ordinary Differential Equations
What we ll do: Lecture Ordinary Differential Equations J. Chaudhry Department of Mathematics and Statistics University of New Mexico Review ODEs Single Step Methods Euler s method (st order accurate) Runge-Kutta
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationKINGS COLLEGE OF ENGINEERING DEPARTMENT OF MATHEMATICS ACADEMIC YEAR / EVEN SEMESTER QUESTION BANK
KINGS COLLEGE OF ENGINEERING MA5-NUMERICAL METHODS DEPARTMENT OF MATHEMATICS ACADEMIC YEAR 00-0 / EVEN SEMESTER QUESTION BANK SUBJECT NAME: NUMERICAL METHODS YEAR/SEM: II / IV UNIT - I SOLUTION OF EQUATIONS
More informationx n+1 = x n f(x n) f (x n ), n 0.
1. Nonlinear Equations Given scalar equation, f(x) = 0, (a) Describe I) Newtons Method, II) Secant Method for approximating the solution. (b) State sufficient conditions for Newton and Secant to converge.
More informationInverse Lax-Wendroff Procedure for Numerical Boundary Conditions of. Conservation Laws 1. Abstract
Inverse Lax-Wendroff Procedure for Numerical Boundary Conditions of Conservation Laws Sirui Tan and Chi-Wang Shu 3 Abstract We develop a high order finite difference numerical boundary condition for solving
More informationFinite Difference Methods (FDMs) 2
Finite Difference Methods (FDMs) 2 Time- dependent PDEs A partial differential equation of the form (15.1) where A, B, and C are constants, is called quasilinear. There are three types of quasilinear equations:
More informationConservation Laws and Finite Volume Methods
Conservation Laws and Finite Volume Methods AMath 574 Winter Quarter, 2011 Randall J. LeVeque Applied Mathematics University of Washington January 3, 2011 R.J. LeVeque, University of Washington AMath 574,
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationFluid Dynamics. Part 2. Massimo Ricotti. University of Maryland. Fluid Dynamics p.1/17
Fluid Dynamics p.1/17 Fluid Dynamics Part 2 Massimo Ricotti ricotti@astro.umd.edu University of Maryland Fluid Dynamics p.2/17 Schemes Based on Flux-conservative Form By their very nature, the fluid equations
More informationarxiv: v1 [physics.comp-ph] 22 Feb 2013
Numerical Methods and Causality in Physics Muhammad Adeel Ajaib 1 University of Delaware, Newark, DE 19716, USA arxiv:1302.5601v1 [physics.comp-ph] 22 Feb 2013 Abstract We discuss physical implications
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationNumerical resolution of discontinuous Galerkin methods for time dependent. wave equations 1. Abstract
Numerical resolution of discontinuous Galerkin methods for time dependent wave equations Xinghui Zhong 2 and Chi-Wang Shu Abstract The discontinuous Galerkin DG method is known to provide good wave resolution
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 informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More information