diff options
Diffstat (limited to '587/CH8/EX8.4')
-rwxr-xr-x | 587/CH8/EX8.4/example8_4.sce | 26 |
1 files changed, 26 insertions, 0 deletions
diff --git a/587/CH8/EX8.4/example8_4.sce b/587/CH8/EX8.4/example8_4.sce new file mode 100755 index 000000000..9c91fedd1 --- /dev/null +++ b/587/CH8/EX8.4/example8_4.sce @@ -0,0 +1,26 @@ +clear;
+clc;
+
+//Example8.4[Pressure Drop in a Water tube]
+Tw=15;//temp of water while entering[degree Celcius]
+rho=999.1;//[kg/m^3]
+mu=1.138*10^(-3);//Viscosity[kg/m.s]
+id=0.05;//Internal diameter[m]
+V=5.5*10^(-3);//Flow rate[m^3/s]
+l=60;//length of tube[m]
+e=0.002*10^(-3);//[m]
+//Solution:-
+v=V/(%pi*(id^2)*(1/4));//Mean Velocity[m/s]
+Re=rho*v*id/mu;
+disp(Re,"Reynolds Number is")
+//Flow is turbulent
+r=e/id;//Relative roughness of the tube
+function[Func]=fric(fac)
+ Func(1)=(1/(fac(1)^(1/2)))+(2*log((0.00004/3.7)+(2.51/(122900*fac(1)^(1/2)))));
+ deff('[Func]=fric(fac)',['fric_1=(1/(fac(1)^(1/2)))+(2*log((0.00004/3.7)+(2.51/(122900*fac(1)^(1/2)))))'])
+endfunction
+disp(xs,"Friction Factor is")
+del_P=xs*l*rho*(v^2)/(2*id);//[kPa]
+disp("Pa",del_P,"The pressure drop is")
+W_pump=V*del_P;//[W]
+disp("W",W_pump,"The required poer input tp overcome the frictional losses in the tube is")
\ No newline at end of file |