Lecture 4.5 Schemes for Parabolic Type Equations
|
|
- Matthew Hawkins
- 6 years ago
- Views:
Transcription
1 Lecture 4.5 Schemes for Parabolic Type Equations 1
2 Difference Schemes for Parabolic Equations One-dimensional problems: Consider the unsteady diffusion problem (parabolic in nature) in a thin wire governed by the differential equation u t u k, ( a, b ), t 0 (4.5.1) Assume that the initial conditions, the distribution of u at t = 0 and the boundary conditions, u at = a and b are given.
3 Forward time central space (FTCS) scheme A simple and easiest scheme to compute the numerical solution of (4.5.1) is the FTCS (forward time and central space) scheme which is an eplicit method. An eplicit scheme uses a stencil in which only one unknown is written in terms of the remaining known values at other stencil points. The FTCS approimation of Eq. (4.5.1) is 1 n 1 n 1 n n n ui ui k u i 1 ui ui 1 t n 1 n n n n t ui ui r ui 1 ui ui 1, r k u ru (1 r) u ru i 1,,, n 1, n 0,1, n 1 n n n i i 1 i i 1 (4.5.) where, the superscript n represents the time level. The discretization and the stencil of the FTCS Eq. (4.5.) is shown in the Fig
4 Fig Discretization and the stencil of FTCS scheme In the Fig , there is only one stencil point at n+1 th time level which is the unknown and three points at n th time level which anyway are known. Therefore using Eq. (4.5.), one unknown at a time at n+1 time level can be computed by varying i = 1,,..., n -1. 4
5 Taylor series epansion of Eq. (4.5.) demonstrates that the FTCS scheme is first order accurate in time and second order in space. Backward time central space (BTCS) scheme The BTCS approimation of Eq. (4.5.1) gives 1 n n 1 1 n n n ui ui k u i 1 ui ui 1 t n n n n 1 t rui 1 (1 r) ui rui 1 ui, r k, i 1,,, n 1, n 1,, (4.5.3) Equation (4.5.3) is an implicit scheme which has more than one unknown at n th time level, therefore one equation alone can t be solved unless it is clubbed with more number of equations to close the system. 5
6 This is done by grouping all the discretized equations at a particular time level and solving them in one step using say, Thomas algorithm if the resultant system is a tri-diagonal one. The same is repeated by incrementing the value of n until the required time level is reached. This type of procedure is called a time marching scheme. 6
7 The computational stencil for the BTCS scheme is as shown in Fig Fig Fig Computational stencil for BTCS scheme At each time step, the scheme Eq. (4.5.3) also can be written as n 1 n 1 n 1 n t rui 1 (1 r) ui rui 1 ui, r k, i 1,,, n 1, 0,1, n (4.5.4) 7
8 Weighted average scheme Weighted average of Eqs. (4.5.1) and (4.5.3), gives 1 n 1 n k n n n n 1 n 1 n 1 ui ui (1 ) u i 1 ui ui 1 ui 1 ui ui 1 t r u (1 r) u r u (1 ) r u (1 (1 ) r) u (1 ) r u n 1 n 1 n 1 n n n i 1 i i 1 i 1 i i 1 t r k, 0 1, i 1,,, n, 0,1, n (4.5.4) For θ equals to 0 and 1, Eq. (4.5.4) gives FTCS and BTCS schemes, respectively. The Taylor series epansion of (4.5.4) gives u u 1 u t u k u k t u k t t t 1 t n 3 n 4 n 4 n i i i i t 3 4 (4.5.4) Therefore, weighted average scheme Eq. (4.5.4) is second order 1 accurate in space and first order in time if. The scheme is 1 also second order accurate in time if. The weighted average scheme with is known as Crank-Nicolson scheme. 8 1
9 Numerical Illustration u Consider, (0,1), t with initial conditions u(,0) = sin t 0 and boundary conditions zero at = 0 and 1. Use step sizes 0. and 0.01 in and t directions, respectively, Compare, after ten time steps, the numerical solutions obtained with FTCS, BTCS and Crank-Nicolson schemes with the t analytical solution e sin. For the step size 0., we have u r k t 1* *0. With r = 0.3, the FTCS, BTCS and Crank-Nicolson schemes are given by 0.3 9
10 u 0.3u 0.4u 0.3u n 1 n n n i i 1 i i 1 0.3u 1.6u 0.3u u n 1 n 1 n 1 n i 1 i i 1 i 0.15u 1.3u 0.15u 0.15u 0.7u 0.15u n 1 n 1 n 1 n n n i 1 i i 1 i 1 i i 1 for i 1,,3, 4, n 0,1,,,9 (4.5.6) Equation (4.5.6) depends only on the value of r and is independent of the step sizes. The solution and the percentage errors generated by the three schemes (FTCS, BTCS and Crank_Nicolson), after marching 10 times in the time direction using Eq. (4.5.6), are compared in the Table
11 X Analytical FTCS BTCS Crank-Nicolson Solution Error Solution Error Solution Error % % % % % % Table Comparison of the analytical and numerical solutions and their errors The initial and boundary conditions in the above computations are taken from the eact solution. 11
12 Two-dimensional problems Consider the unsteady diffusion over a flat plate governed by the differential equation u u u k, (, y) ( a, b) X ( c, d), t 0 t y (4.5.7) Assume that the initial and boundary conditions on u are known. Eplicit scheme Approimating the time derivative with forward difference and space derivatives with central differences gives a scheme u u k u u u u u u t y n 1 n n n n n n n i, j i, j i 1, j i, j i 1, j i, j 1 i, j i, j 1 u u r u u u r u u u n 1 n n n n n n n i, j i, j 1 i 1, j i, j i 1, j i, j 1 i, j i, j 1 n 1 n n n n n ui, j r ui, j 1 ru 1 i 1, j (1 r1 r ) u i, j ru 1 i 1, j r ui, j 1 t r k, r k t y 1 i 1,,, n 1, j 1,,, n 1, n 0,1, y (4.5.8) 1
13 Here, y is the step length in y direction. Taylor series epansion of Eq. (4.5.8) shows that the eplicit scheme is first order accurate in time and second order in space (both in and y directions). Weighted average or Crank-Nicolson type of approimation to (4.5.7) gives a penta-diagonal system like Eq. (4.4.) at the n+1 th time level, solving such a system is very epensive computationally, therefore, alternatively ADI method can be developed as follows: n n n n n n n n i, j i, j i 1, j i, j i 1, j i, j 1 i, j i, j 1 k u u u u u u u u t / y (4.5.9) n n n n n 1 n 1 n 1 n 1 i, j i, j i 1, j i, j i 1, j i, j 1 i, j i, j 1 k u u u u u u u u t / y 13
14 First on j is constant lines, using the first tri-diagonal part of Eq. (4.5.9), solution at n+1/ time level is obtained. In the second step, using the solution at the n+1/ time level over the i is constant lines and using second tri-diagonal part of Eq. (4.5.9), the solution is marched to the n+1 time level. Therefore, for each time level, Eq. (4.5.9) gives (n -1 + n y -1) tridiagonal systems in and y directions which are mush easier to solve using Thomas algorithm than the penta-diagonal system (of size L X L, where L = (n -1) * (n y -1)) appears in the weighted average or Crank-Nicolson schemes for two dimensional problems. 14
15 Convergence Hear, convergence means, the convergence of the numerical solution to the analytical solution For parabolic and also for all time dependent problems, the convergence of the numerical solution to the corresponding analytical solution is carried out through the testing for consistency and stability since, according to La equivalence theorem, Consistency and stability are necessary and sufficient conditions for convergence of the finite difference solutions of any time dependent problem. 15
16 Consistency of a Numerical Scheme Under the limiting case of step lengths tending to zero, if a finite difference scheme converge to the corresponding differential equation then such a scheme is called consistent. Mathematically, it is tested by looking at the truncation error of the scheme as the step lengths tend to zero. If the truncation error tends to zero as the step lengths tend zero then the numerical scheme is said to be consistent. 16
17 Numerical Illustration Consider the Taylor series epansion Eq. (4.5.5) of the Weighted average scheme Eq. (4.5.4) r u (1 r) u r u (1 ) r u (1 (1 ) r) u (1 ) r u given by n 1 n 1 n 1 n n n i 1 i i 1 i 1 i i 1 u u 1 u t u k u k t u k t t t 1 t n 3 n 4 n 4 n i i i i t 3 4 Taking step lengths and t tending to zero, the series epansion converges to the governing equation k, therefore, the weighted average scheme is consistent. u t u 17
18 Stability of a Numerical Scheme Stability of a numerical scheme deals with the growth of the rounding errors during the time marching process. Let us illustrate, the stability through the following observation Repeat the computations of Eq. (4.5.6) once again with time steps and 0.0 (that is, for the value of r as 0.45 and 0.55, respectively) for FTCS and Crank-Nicolson (CN) schemes The corresponding solution with r =.45 and.55 are presented in Tables 4.5. and 4.5.3, respectively 18
19 X Analytical FTCS Crank-Nicolson Solution Error Solution Error % % % % Table 4.5. Comparison of the solution and errors with r = X Analytical FTCS Crank-Nicolson Solution Error Solution Error % % % % Table Comparison of the solution and errors with r =
20 Its clear from the Tables 4.5. and that, with r = 0.45, the FTCS scheme produces solutions with errors less than.% and CN scheme with errors less than 1.3%. However, if the value of r is increased to 0.55, the CN scheme continues to give solutions with similar errors while the errors in FTCS scheme increased by many folds. The behavior of increasing errors with FTCS becomes worse and completely dominated by these errors if we still continue the computations for higher time levels. This is due to the instable nature of the FTCS scheme when the value of r is greater than 0.5. Mathematically this can be understood using the following analysis: 0
21 n Let (u) = 0 be a linear difference scheme and is its numerical n n solution, U is the eact solution and is its error at n th i E i time level at the nodal point i then we have n n n u (4.5.10) i Ui Ei and n n n n ( ui ) = ( U E ) = ( ) + ( E n ) = ( E n ) = 0 (4.5.11) i i i i U i That is, the error also satisfies the same difference equation which numerical solution satisfies. Therefore, one can study the behavior of the error by studying the numerical solution itself. u i 1
22 Further, if the numerical solution is assumed to be periodic, achieved by reflecting the solution in the region (0, L) in to (-L, L) and epressible it in terms of finite Fourier series (since the domain is of finite length which is discretized with finite step length) then it can be written as where j i j u K e K e K e (4.5.1) N is the number of points in the discretization, (N=L/ ) I n K j 1 is the amplitude of the j th harmonic, k j is the wave number n N n Ik N n Ik i N n Ii i j j j j N j N j N is the phase angle given by k j.
23 Note: k j varies from N to N instead of - to, because the maimum and minimum resolvable wavelengths (λ) are only L and, respectively and the maimum wavelength L is discretized with N+1 points, that is from N to N. In the actual computation, due to the linear nature of the difference scheme, it is enough to use one Fourier mode, instead of (4.5.1) and looking at the raise or damping of the amplification factor G which is defined as the ratio of the amplitude at n+1 and n th time levels, that G is defined as n 1 Ki G (4.5.13) K n Any scheme (linear) is said to be stable if Otherwise is said to be unstable. i G < 1 (4.5.1) 3
24 Numerical Eample: Discuss the stability criteria of the weighted average scheme r u (1 r) u r u (1 ) r u (1 (1 ) r) u (1 ) r u n 1 n 1 n 1 n n n i 1 i i 1 i 1 i i 1 (4.5.15) n n Ii Substituting u K e in Eq. (4.5.15) and simplifying for i amplification factor G, gives G n 1 I I K (1 ) re (1 (1 ) r) (1 ) re 1 4(1 ) r sin n I I K re (1 r) re 1 4 r sin (4.5.16) For θ = 0, that is, for FTCS scheme, G <1 implies which is true whenever the value of r < 1/. 1 4r sin 1 4
25 Therefore, FTCS scheme is only conditionally stable. However, for θ greater or equal to ½, Eq. (4.5.16) is unconditionally stable. Once again looking at the Tables 4.5. and 4.5.3, it is clear that for r = 0.45, FTCS scheme is able to produce accurate solutions but the round of errors are dominated when the value of r is raised to 0.55 because FTCS scheme is not stable at r = On the other hand due its stable nature of CN scheme for all values of r, the round of errors are under control at both r = 0.45 and
26 Eample: Analyze the stability of the scheme t u u r u u u r u u u, r k, r k n 1 n n n n n n n i, j i, j 1 i 1, j i, j i 1, j i, j 1 i, j i, j 1 1 t y n Ii Ii n n y Substituting either Ke or Ke for u ij, depending on variation of u in or y directions we get n 1 K i i i y G 1 r1 e e r e e n K 1 r cos 1 r cos r1sin 4rsin (4.5.17) Then the condition G <1 is satisfied whenever r 1 +r < ½. If r 1 = r then condition for the stability of the two dimensional diffusion problem is r < ¼ which is even stronger than the 6 condition of the corresponding 1D case. y y i y
27 Eercise Problems Discuss the consistency and stability criteria for the following schemes: 1. Crank-Nicolson (CN) scheme for the two-dimensional diffusion problem u u r u u 3u u u n 1 n n n n n n i, j i, j 1 i, j i 1, j i, j i 1, j i, j r u u 3u u u n n n n n i, j 1 i, j 1 i, j i, j 1 i, j 3. For the scheme obtained by discretizing the time derivative with forward difference approimation and space derivatives with second order central difference approimations of the equation, u a u where a, a y and μ are a u u u y t y y 7 constants.
28 Summary of Lecture 4.5 Various finite difference approimations of parabolic equations like unsteady diffusion equations are introduced in this lecture. END OF LECTURE 4.5 8
Finite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
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 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 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 informationMulti-Factor Finite Differences
February 17, 2017 Aims and outline Finite differences for more than one direction The θ-method, explicit, implicit, Crank-Nicolson Iterative solution of discretised equations Alternating directions implicit
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 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 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 and Its Difference Approximations 1-D Initial Boundary Value
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 informationBeam Propagation Method Solution to the Seminar Tasks
Beam Propagation Method Solution to the Seminar Tasks Matthias Zilk The task was to implement a 1D beam propagation method (BPM) that solves the equation z v(xz) = i 2 [ 2k x 2 + (x) k 2 ik2 v(x, z) =
More informationChapter 5. Formulation of FEM for Unsteady Problems
Chapter 5 Formulation of FEM for Unsteady Problems Two alternatives for formulating time dependent problems are called coupled space-time formulation and semi-discrete formulation. The first one treats
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Packet The Group Velocity Phys 375 Overview and Motivation: Last time we looked at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
More informationCharacteristic finite-difference solution Stability of C C (CDS in time/space, explicit): Example: Effective numerical wave numbers and dispersion
Spring 015 Lecture 14 REVIEW Lecture 13: Stability: Von Neumann Ex.: 1st order linear convection/wave eqn., F-B scheme Hyperbolic PDEs and Stability nd order wave equation and waves on a string Characteristic
More informationFUNDAMENTALS OF FINITE DIFFERENCE METHODS
FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur Acknowledgements It has been a pleasure
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 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 informationq t = F q x. (1) is a flux of q due to diffusion. Although very complex parameterizations for F q
! Revised Tuesday, December 8, 015! 1 Chapter 7: Diffusion Copyright 015, David A. Randall 7.1! Introduction Diffusion is a macroscopic statistical description of microscopic advection. Here microscopic
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 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 informationBTCS Solution to the Heat Equation
BTCS Solution to the Heat Equation ME 448/548 Notes Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@mepdxedu ME 448/548: BTCS Solution to the Heat Equation Overview
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 information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
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 informationErrors Intensive Computation
Errors Intensive Computation Annalisa Massini - 2015/2016 OVERVIEW Sources of Approimation Before computation modeling empirical measurements previous computations During computation truncation or discretization
More informationA Propagating Wave Packet The Group Velocity
Lecture 7 A Propagating Wave Pacet The Group Velocity Phys 375 Overview and Motivation: Last time we looed at a solution to the Schrödinger equation (SE) with an initial condition (,) that corresponds
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 informationAn Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems
An Exponential High-Order Compact ADI Method for 3D Unsteady Convection Diffusion Problems Yongbin Ge, 1 Zhen F. Tian, 2 Jun Zhang 3 1 Institute of Applied Mathematics and Mechanics, Ningxia University,
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationChapter 4. Nonlinear Hyperbolic Problems
Chapter 4. Nonlinear Hyperbolic Problems 4.1. Introduction Reading: Durran sections 3.5-3.6. Mesinger and Arakawa (1976) Chapter 3 sections 6-7. Supplementary reading: Tannehill et al sections 4.4 and
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 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 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 informationComputational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras
Computational Techniques Prof. Sreenivas Jayanthi. Department of Chemical Engineering Indian institute of Technology, Madras Module No. # 05 Lecture No. # 24 Gauss-Jordan method L U decomposition method
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationThe Fundamental Theorem of Calculus Part 3
The Fundamental Theorem of Calculus Part FTC Part Worksheet 5: Basic Rules, Initial Value Problems, Rewriting Integrands A. It s time to find anti-derivatives algebraically. Instead of saying the anti-derivative
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 7:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture7/7_1.htm 1 of 1 6/20/2012 12:26 PM The Lecture deals with: Errors and Stability Analysis file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture7/7_2.htm
More informationPlatzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Introduction to PDEs and Numerical Methods Lecture 6: Numerical solution of the heat equation with FD method: method of lines, Euler
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 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 informationModule 3: BASICS OF CFD. Part A: Finite Difference Methods
Module 3: BASICS OF CFD Part A: Finite Difference Methods THE CFD APPROACH Assembling the governing equations Identifying flow domain and boundary conditions Geometrical discretization of flow domain Discretization
More informationPartial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction
Page of 8 Partial Differential Equations (PDEs) and the Finite Difference Method (FDM). An introduction FILE:Chap 3 Partial Differential Equations-V6. Original: May 7, 05 Revised: Dec 9, 06, Feb 0, 07,
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 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 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 informationBasics of Discretization Methods
Basics of Discretization Methods In the finite difference approach, the continuous problem domain is discretized, so that the dependent variables are considered to exist only at discrete points. Derivatives
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 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 informationConditional stability of Larkin methods with non-uniform grids
Theoret. Appl. Mech., Vol.37, No., pp.139-159, Belgrade 010 Conditional stability of Larkin methods with non-uniform grids Kazuhiro Fukuyo Abstract Stability analysis based on the von Neumann method showed
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 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 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 information2.13 Linearization and Differentials
Linearization and Differentials Section Notes Page Sometimes we can approimate more complicated functions with simpler ones These would give us enough accuracy for certain situations and are easier to
More informationExact and Approximate Numbers:
Eact and Approimate Numbers: The numbers that arise in technical applications are better described as eact numbers because there is not the sort of uncertainty in their values that was described above.
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 Solution Techniques in Mechanical and Aerospace Engineering
Numerical Solution Techniques in Mechanical and Aerospace Engineering Chunlei Liang LECTURE 3 Solvers of linear algebraic equations 3.1. Outline of Lecture Finite-difference method for a 2D elliptic PDE
More informationA CCD-ADI method for unsteady convection-diffusion equations
A CCD-ADI method for unsteady convection-diffusion equations Hai-Wei Sun, Leonard Z. Li Department of Mathematics, University of Macau, Macao Abstract With a combined compact difference scheme for the
More informationChapter 2 Finite-Difference Discretization of the Advection-Diffusion Equation
Chapter Finite-Difference Discretization of the Advection-Diffusion Equation. Introduction Finite-difference methods are numerical methods that find solutions to differential equations using approximate
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 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 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 informationTENSOR TRANSFORMATION OF STRESSES
GG303 Lecture 18 9/4/01 1 TENSOR TRANSFORMATION OF STRESSES Transformation of stresses between planes of arbitrar orientation In the 2-D eample of lecture 16, the normal and shear stresses (tractions)
More informationThe Finite Difference Method
Chapter 5. The Finite Difference Method This chapter derives the finite difference equations that are used in the conduction analyses in the next chapter and the techniques that are used to overcome computational
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 information+ y = 1 : the polynomial
Notes on Basic Ideas of Spherical Harmonics In the representation of wavefields (solutions of the wave equation) one of the natural considerations that arise along the lines of Huygens Principle is the
More informationEuler-Maclaurin summation formula
Physics 4 Spring 6 Euler-Maclaurin summation formula Lecture notes by M. G. Rozman Last modified: March 9, 6 Euler-Maclaurin summation formula gives an estimation of the sum N in f i) in terms of the integral
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationThis theorem guarantees solutions to many problems you will encounter. exists, then f ( c)
Maimum and Minimum Values Etreme Value Theorem If f () is continuous on the closed interval [a, b], then f () achieves both a global (absolute) maimum and global minimum at some numbers c and d in [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 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 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 informationName of the Student: Unit I (Solution of Equations and Eigenvalue Problems)
Engineering Mathematics 8 SUBJECT NAME : Numerical Methods SUBJECT CODE : MA6459 MATERIAL NAME : University Questions REGULATION : R3 UPDATED ON : November 7 (Upto N/D 7 Q.P) (Scan the above Q.R code for
More information1 Finite difference example: 1D implicit heat equation
1 Finite difference example: 1D implicit heat equation 1.1 Boundary conditions Neumann and Dirichlet We solve the transient heat equation ρc p t = ( k ) (1) on the domain L/2 x L/2 subject to the following
More information(a) Show that there is a root α of f (x) = 0 in the interval [1.2, 1.3]. (2)
. f() = 4 cosec 4 +, where is in radians. (a) Show that there is a root α of f () = 0 in the interval [.,.3]. Show that the equation f() = 0 can be written in the form = + sin 4 Use the iterative formula
More informationlecture 7: Trigonometric Interpolation
lecture : Trigonometric Interpolation 9 Trigonometric interpolation for periodic functions Thus far all our interpolation schemes have been based on polynomials However, if the function f is periodic,
More informationComputer Aided Design of Thermal Systems (ME648)
Computer Aided Design of Thermal Systems (ME648) PG/Open Elective Credits: 3-0-0-9 Updated Syallabus: Introduction. Basic Considerations in Design. Modelling of Thermal Systems. Numerical Modelling and
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 information2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13
2.29 Numerical Fluid Mechanics Fall 2009 Lecture 13 REVIEW Lecture 12: Classification of Partial Differential Equations (PDEs) and eamples with finite difference discretizations Parabolic PDEs Elliptic
More informationIntroduction to Differential Equations
Math0 Lecture # Introduction to Differential Equations Basic definitions Definition : (What is a DE?) A differential equation (DE) is an equation that involves some of the derivatives (or differentials)
More informationChapter 6. Nonlinear Equations. 6.1 The Problem of Nonlinear Root-finding. 6.2 Rate of Convergence
Chapter 6 Nonlinear Equations 6. The Problem of Nonlinear Root-finding In this module we consider the problem of using numerical techniques to find the roots of nonlinear equations, f () =. Initially we
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 informationCS205b/CME306. Lecture 4. x v. + t
CS05b/CME306 Lecture 4 Time Integration We now consider seeral popular approaches for integrating an ODE Forward Euler Forward Euler eolution takes on the form n+ = n + Because forward Euler is unstable
More informationA Padé approximation to the scalar wavefield extrapolator for inhomogeneous media
A Padé approimation A Padé approimation to the scalar wavefield etrapolator for inhomogeneous media Yanpeng Mi, Zhengsheng Yao, and Gary F. Margrave ABSTRACT A seismic wavefield at depth z can be obtained
More informationNewton's Laws You should be able to state these laws using both words and equations.
Review before first test Physical Mechanics Fall 000 Newton's Laws You should be able to state these laws using both words and equations. The nd law most important for meteorology. Second law: net force
More informationPart E1. Transient Fields: Leapfrog Integration. Prof. Dr.-Ing. Rolf Schuhmann
Part E1 Transient Fields: Leapfrog Integration Prof. Dr.-Ing. Rolf Schuhmann MAXWELL Grid Equations in time domain d 1 h() t MC e( t) dt d 1 e() t M Ch() t j( t) dt Transient Fields system of 1 st order
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationQuantum Dynamics. March 10, 2017
Quantum Dynamics March 0, 07 As in classical mechanics, time is a parameter in quantum mechanics. It is distinct from space in the sense that, while we have Hermitian operators, X, for position and therefore
More informationSPS Mathematical Methods
SPS 2281 - Mathematical Methods Assignment No. 2 Deadline: 11th March 2015, before 4:45 p.m. INSTRUCTIONS: Answer the following questions. Check our answer for odd number questions at the back of the tetbook.
More informationChapter Parabolic Partial Differential Equations
hapter. Parabolic Partial Differential Equations After reading this chapter, you should be able to:. Use numerical methods to solve parabolic partial differential equations by explicit, implicit, and rank-nicolson
More informationα x x 0 α x x f(x) α x x α x ( 1) f(x) x f(x) x f(x) α x = α x x 2
Quadratic speedup for unstructured search - Grover s Al- CS 94- gorithm /8/07 Spring 007 Lecture 11 01 Unstructured Search Here s the problem: You are given an efficient boolean function f : {1,,} {0,1},
More informationFurther factorising, simplifying, completing the square and algebraic proof
Further factorising, simplifying, completing the square and algebraic proof 8 CHAPTER 8. Further factorising Quadratic epressions of the form b c were factorised in Section 8. by finding two numbers whose
More informationSOLUTIONS BY SUBSTITUTIONS
25 SOLUTIONS BY SUBSTITUTIONS 71 25 SOLUTIONS BY SUBSTITUTIONS REVIEW MATERIAL Techniques of integration Separation of variables Solution of linear DEs INTRODUCTION We usually solve a differential equation
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationSeismic Waves Propagation in Complex Media
H4.SMR/1586-1 "7th Workshop on Three-Dimensional Modelling of Seismic Waves Generation and their Propagation" 5 October - 5 November 004 Seismic Waves Propagation in Comple Media Fabio ROMANELLI Dept.
More informationIntroduction to Differential Equations. National Chiao Tung University Chun-Jen Tsai 9/14/2011
Introduction to Differential Equations National Chiao Tung Universit Chun-Jen Tsai 9/14/011 Differential Equations Definition: An equation containing the derivatives of one or more dependent variables,
More informationA WAVELET-TAYLOR GALERKIN METHOD FOR PARABOLIC AND HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
International Journal of Computational Methods Vol. 2, No. (25) 75 97 c World Scientific Publishing Company A WAVELET-TAYLOR GALERKIN METHOD FOR PARABOLIC AND HYPERBOLIC PARTIAL DIFFERENTIAL EQUATIONS
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 informationPartial Differential Equations
Next: Using Matlab Up: Numerical Analysis for Chemical Previous: Ordinary Differential Equations Subsections Finite Difference: Elliptic Equations The Laplace Equations Solution Techniques Boundary Conditions
More informationUniversity of Alberta ENGM 541: Modeling and Simulation of Engineering Systems Laboratory #5
University of Alberta ENGM 54: Modeling and Simulation of Engineering Systems Laboratory #5 M.G. Lipsett, Updated 00 Integration Methods with Higher-Order Truncation Errors with MATLAB MATLAB is capable
More informationY m = y n e 2πi(m 1)(n 1)/N (1) Y m e 2πi(m 1)(n 1)/N (2) m=1
The Discrete Fourier Transform (Bretherton notes): 1 Definition Let y n, n = 1,..., N be a sequence of N possibly comple values. The discrete Fourier transform () of this sequence is the sequence Y m,
More information