blob: 11c8becc25d8f968391bd545ce64b3b72c9aa5f5 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
clc;
clear all;
disp("displacement thickness")
disp("velocity distribution = u/U=y/d")
disp("i) the displacement thickness, d*:")
disp("d*=intergrate((1-u/U,U,0,d)")
disp("d*=d/2")
disp("ii) the momentum thickness, th")
disp("th = integrate(u(1-u/U)/U,u,0,d)")
disp("th=d/6")
disp("iii) The energy thickness de")
disp("de=integrate (u/U*(1-(u/U)^2),U,0,d")
disp("de=d/4")
disp("iv) The value of d*/th")
disp("d*/th =3")
|