31 computational grid used for the nite di erence

Info iconThis preview shows page 1. Sign up to view the full content.

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

Unformatted text preview: re r = (2.3.3b) t 2: x The initial and boundary conditions (2.3.1b, 2.3.1c) are 0 Uj =( j ) x j =0 1 ::: J (2.3.3c) 18 Finite Di erence Methods n 0 U = n UJ =0 n> 0 (2.3.3d) : The computational stencil for (2.3.3a) is shown in Figure 2.3.2. The parameter analogous to the Courant number (2.2.3) for the kinematic wave equation (2.2.1). 11 00 11 00 is n+1 11 00 j-1 r 11 00 j n j+1 Figure 2.3.2: Computational stencil of the forward time-centered space scheme (2.3.3a) for the heat conduction equation (2.3.1a). The solution of (2.3.3a) is obtained in the same manner as the nite di erence solutions of (2.2.2b) and (2.2.4) for the kinematic wave equation (2.2.1). Thus, using the initial data (2.3.3c), we calculate a solution j1 at the interior mesh points, =12 ; 1, of time level 1 using (2.3.3a) with ranging from 1 to ; 1. The 1 boundary conditions (2.3.3d) with = 1 determine 01 and J . Knowing the discrete solution at time level 1, we proceed to determine it at time level 2, etc. in the sam...
View Full Document

This document was uploaded on 03/16/2014 for the course CSCI 6840 at Rensselaer Polytechnic Institute.

Ask a homework question - tutors are online