readme.txt

Computtional Fluid Dynamics (CFD) steady 1-Dimensional convection and diffusion using central diff

Computtional Fluid Dynamics (CFD)
steady 1-Dimensional convection and diffusion
using central differencing, upwind differencing,
hybrid differencing and power law schemes,
programmed in MATLAB

9th March 2005
-------------------------------------------------------
Prepared by:
Mohd Anuar Abd Aziz
Student Me 071642
Mechanical Department
Tenaga Nasional University.
(UNITEN)
Malaysia.


e-mail: [e-mail]
Reference:-

Versteeg, H.K., Malalasekera,W, An introduction to computational
fluid dynamics, The finite volume methods,Chapter 5.

Lectures and Notes by Dr Zamri, vice dean of College of Engineering

------------------------------------------------------

MATLAB Version: 6.5 R13.

------------------------------------------------------
MEHB463: COMPUTATIONAL FLUID DYNAMICS
Assignment No 3

The above three methods of iteration is used as in accordance 
to completion of assignment 3 in the CFD class of this term 
under the supervision of Dr Zamri.

The problem describes a diffusion and convection terms which
is to be solved and discretised using the above mentioned schemes
Given boundary conditions, phia=4,x=0 and phib=1,x=L=1.

Considering two cases to demonstrate the advantage and disadvantage
of 4 properties of discretization schemes;(1)conservativeness,
(2)transportiveness, (3) boundedness, (4)accuracy.

case1: u=0.2 m/s, rho=1,r=1,mesh,N=10
case2: u=3.0 m/s, rho=1,r=1,mesh,N=10

 

Program Listing
The encoding is done using Matlab Ver 6.5, release R13.
Below describes the functions and files in the running of the program. 
Please see attached program listings for the full program.

cfd3.m	Main function. Will call tdma function
tdma.m	Function that receives velocity, number of mesh , and S=types of scheme
		and return the value of x value and the phi that is solved using
		tridiagonal matrices (thomas method).

Please run cfd3.m to run the comparison program of the schemes.
>>cfd1 <enter>

Note: Initials conditions of rho,r,L, etc, should be done in the tdma.m mfile.
      I only made the fuction to receive types, mesh, and velocity of the scheme.
      The program is actually not done in general form as to cater for the needs of students.
      Enjoy guys.