blob: f0c445893139194d89cca33ca7f49d64c3a8adc6 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
//exapple 3.1
clc; funcprot(0);
// Initialization of Variable
rho=998;
mu=1.002/1000;
x=48/100;
u=19.6/100;
x1=30/100;
b=2.6;
//calculation
//part1
disp("fluid in boundary layer would be entirely in streamline motion ");
Re=rho*x*u/mu;
printf("reynolds no is %.2e",Re);
//part 2
Re1=rho*x1*u/mu;
delta=x1*4.64*Re1^-.5;
disp(delta*1000,"boundary layer width in (mm):");
//part3
y=0.5*delta;//middle of boundary layer
ux=3/2*u*y/delta-.5*u*(y/delta)^3;
disp(ux*100,"velocity of water in (cm/s):");
//part4
R=0.323*rho*u^2*Re1^-0.5;
disp(R,"shear stress at 30cm in (N/m^2):");
//part5
Rms=0.646*rho*u^2*Re^-0.5;
disp(Rms,"mean shear stress experienced over whole plate in (N/m^2)");
//part6
F=Rms*x*b;
disp(F,"total force experienced by the plate in (N)")
|
