lec5_notes

lec5_notes - 16.920J/SMA 5212 Numerical Methods for Partial...

Info iconThis preview shows pages 1–6. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: 16.920J/SMA 5212 Numerical Methods for Partial Differential Equations Lecture 5 Finite Differences: Parabolic Problems B. C. Khoo Thanks to Franklin Tan 19 February 2003 16.920J/SMA 5212 Numerical Methods for PDEs 2 OUTLINE Governing Equation Stability Analysis 3 Examples Relationship between and h Implicit Time-Marching Scheme Summary Slide 2 GOVERNING EQUATION Consider the Parabolic PDE in 1-D If viscosity Diffusion Equation If thermal conductivity Heat Conduction Equation Slide 3 STABILITY ANALYSIS Discretization Keeping time continuous, we carry out a spatial discretization of the RHS of [ ] 2 2 0, u u x t x = subject to at 0, at u u x u u x = = = = x = x = u u ( 29 , ? u x t = 2 2 u u t x = x = x = x 1 x 2 x 1 N x- N x 16.920J/SMA 5212 Numerical Methods for PDEs 3 Slide 4 STABILITY ANALYSIS Discretization which is second-order accurate. Schemes of other orders of accuracy may be constructed. Slide 5 Construction of Spatial Difference Scheme of Any Order p The idea of constructing a spatial difference operator is to represent the spatial differential operator at a location by the neighboring nodal points, each with its own weightage. The order of accuracy, p of a spatial difference scheme is represented as ( 29 p O x . Generally, to represent the spatial operator to a higher order of accuracy, more nodal points must be used. Consider the following procedure of determining the spatial operator j du dx up to the order of accuracy ( 29 2 O x : There is a total of 1 grid points such that , 0,1, 2,...., j N x j x j N + = = 2 2 Use the Central Difference Scheme for u x 2 1 1 2 2 2 2 ( ) j j j j u u u u O x x x +-- + = + j- 2 j- 1 j j +1 j+2 j du dx 16.920J/SMA 5212 Numerical Methods for PDEs 4 1. Let j du dx be represented by u at the nodes j- 1, j , and j +1 with 1 - , and 1 being the coefficients to be determined, i.e. ( 29 1 1 1 1 p j j j j du u u u O x dx -- + + + + = 2. Seek Taylor Expansions for 1 j u- , j u and 1 j u + about j u and present them in a table as shown below. (Note that p is not known a priori but is determined at the end of the analysis when the s are made known.) u j u j u j u j j u 0 1 0 0 1 1 j u -- 1 - 1 x -- 2 1 1 2 x - 3 1 1 6 x -- j u 0 0 0 1 1 j u + 1 1 x 2 1 1 2 x 3 1 1 6 x 1 1 k j k j k k u u = + =- + 1 S 2 S 3 S 4 S ( 1 ) This column consists of all the terms on the LHS of (1). Each cell in this row comprises the sum of its corresponding column. 16.920J/SMA 5212 Numerical Methods for PDEs 5 where 1 1 2 3 4 1 .......
View Full Document

This note was uploaded on 11/08/2011 for the course AERO 16.872 taught by Professor Danielhastings during the Fall '03 term at MIT.

Page1 / 31

lec5_notes - 16.920J/SMA 5212 Numerical Methods for Partial...

This preview shows document pages 1 - 6. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online