blob: d62873e680b2e127de99388336ad7b3d172f5e52 (
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
|
clc
clear
printf("example 8.7 page number 367\n\n")
//to find the flow rate in an orifice
density_of_water = 1000; //in kg/m3
viscosity = 1*10^-3; //in Pa-s
pipe_diameter = 250; //in mm
orifice_diameter = 50; // in mm
density_of_mercury = 13600; // in mm
manometer_height = 242; //in mm
//calculation
height_water_equivalent = (density_of_mercury-density_of_water)*(manometer_height*10^-3)/(density_of_water) //in m
//assuming Re>30000
Co = 0.61;
velocity = Co*(2*9.8*height_water_equivalent/(1-(orifice_diameter/pipe_diameter)^4))^0.5; //in m/s
//checking Reynold's number
Re = (orifice_diameter*10^-3*velocity*density_of_water)/viscosity;
printf("reynolds number = %f\nwhich is greater than 30000",Re)
if Re>30000 then printf("\n\nvelocity of water = %f m/s",velocity)
end
rate_of_flow = (3.14*(orifice_diameter*10^-3)^2/4)*velocity*density_of_water;
printf("\n\nrate of flow = %f litre/s",rate_of_flow)
|
