Scribd
Upload
46 views2 pages

Numerical Methods of Thermo-Fluid Dynamics I: Deliverable Task Ii: Numerical Solution of Boundary Layer Equation

This document outlines the second deliverable task for the course Numerical Methods of Thermo-Fluid Dynamics. Students are asked to solve the boundary layer equations over a flat plate numerically using MATLAB. This involves: 1) Discretizing the governing equations, writing iterative formulas, and commenting on the schemes. 2) Developing a MATLAB program to solve the equations over the domain. 3) Generating plots of the numerical solutions and explaining the results. 4) Comparing profiles to analytical Blasius solutions and stating conclusions. The report and code are due on January 8th, 2018.

Uploaded by

Jishnu Jayaraj
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
46 views2 pages

Numerical Methods of Thermo-Fluid Dynamics I: Deliverable Task Ii: Numerical Solution of Boundary Layer Equation

This document outlines the second deliverable task for the course Numerical Methods of Thermo-Fluid Dynamics. Students are asked to solve the boundary layer equations over a flat plate numerically using MATLAB. This involves: 1) Discretizing the governing equations, writing iterative formulas, and commenting on the schemes. 2) Developing a MATLAB program to solve the equations over the domain. 3) Generating plots of the numerical solutions and explaining the results. 4) Comparing profiles to analytical Blasius solutions and stating conclusions. The report and code are due on January 8th, 2018.

Uploaded by

Jishnu Jayaraj
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

Numerical Methods of Thermo-Fluid Dynamics I

Winter Semester 2017-2018

DELIVERABLE TASK II:


Numerical Solution of Boundary Layer Equation
Given: Monday, 27/11/2017
Deadline: 8/1/2018

Chair of Fluid Mechanics


Department of Biochemical Engineering, Technical Faculty
Friedrich-Alexander University Erlangen-Nuremberg

Prof. Dr. Marc Avila, M.Sc. Nan Chen


2 Numerical Methods of Thermo-Fluid Dynamics

Deliverable Task II

The goal of the second deliverable task is to solve the boundary layer equation over a flat plate:

u v
+ =0 (1)
x y
u u 1 2 u
u +v = . (2)
x y Re y2

The Reynolds number is defined as Re = u L/, where is the kinematic viscosity, L is the plate
length and u is the free-stream velocity of the outer flow. Here the velocity components u and
v are dimensionless, u = u/u and v = v/u , where u and v are the dimensional velocities.
Lengths have been made dimensionless by using the plate length, e.g. x = x/L. The boundary
conditions are as follows:

y=0: u=v=0 no-slip condition (3)


y: u 1 (u u ) free outer flow. (4)

(5p) Discretise equation (1). One possible scheme for u


x is as follows:

1 ui,j ui1,j ui,j1 ui1,j1


   
u
= + (5)
x i,j 2 x x

Discretise y ,
v
write the iterative formula for v based on your schemes, show the order of
accuracy. Note that you could choose other schemes to approximate x .
u

(10p) Discretise equation (2), write the iterative formula for u, give comments on your
scheme (explicit/implicit, order of accuracy, stability constraint, etc.)

(25p) Write a MATLAB/OCTAVE program that solves the boundary layer equation (0 x
1 and 0 y 2, where = 5 is the boundary layer thickness at x = 1).
Re

(30p) Calculate the solution for Re = 10000, generate 2 figures for the numerical results of
u and v of the whole domain, explain your results.

(30p) Consider u = 10[m/s] and L = 1[m]. For Re = 10000, compare your u and v
velocity profiles obtained above at x = 0.0005 and x = 0.5 with the results you obtained
from the Blasius equation in Exercise Sheet. Is there a good agreement between the two?
Why? Explain the results and state your conclusions.

The Deliverable Task II contains one report and one matlab code file, which must be submitted to
[email protected] before January 8th, 2018, 10:00am.

You might also like