diff options
Diffstat (limited to '1268/CH11')
| -rw-r--r-- | 1268/CH11/EX11.1/a_1.sce | 13 | ||||
| -rw-r--r-- | 1268/CH11/EX11.10/a_10.sce | 9 | ||||
| -rw-r--r-- | 1268/CH11/EX11.14/a_14.sce | 11 | ||||
| -rw-r--r-- | 1268/CH11/EX11.15/a_15.sce | 5 | ||||
| -rw-r--r-- | 1268/CH11/EX11.17/a_17.sce | 10 | ||||
| -rw-r--r-- | 1268/CH11/EX11.19/a_19.sce | 10 | ||||
| -rw-r--r-- | 1268/CH11/EX11.2/a_2.sce | 7 | ||||
| -rw-r--r-- | 1268/CH11/EX11.3/a_3.sce | 22 | ||||
| -rw-r--r-- | 1268/CH11/EX11.6/a_6.sce | 10 | ||||
| -rw-r--r-- | 1268/CH11/EX11.7/a_7.sce | 9 | ||||
| -rw-r--r-- | 1268/CH11/EX11.8/a_8.sce | 10 | ||||
| -rw-r--r-- | 1268/CH11/EX11.9/a_9.sce | 6 |
12 files changed, 122 insertions, 0 deletions
diff --git a/1268/CH11/EX11.1/a_1.sce b/1268/CH11/EX11.1/a_1.sce new file mode 100644 index 000000000..b13af849d --- /dev/null +++ b/1268/CH11/EX11.1/a_1.sce @@ -0,0 +1,13 @@ +clc;
+disp("Example A.1")
+g=9.81
+density=1000 // of water in kg/m^3
+densitym=13600 // of mercury in kg/m^3
+h=0.1 // in m
+p1=density*g*h
+p2=p1+(densitym*g*h)
+waterhead=p2/(density*g)
+hghead=p2/(g*densitym)
+
+disp(waterhead,"Head of water is ")
+disp(hghead,"Head of mercury is ")
diff --git a/1268/CH11/EX11.10/a_10.sce b/1268/CH11/EX11.10/a_10.sce new file mode 100644 index 000000000..558103a91 --- /dev/null +++ b/1268/CH11/EX11.10/a_10.sce @@ -0,0 +1,9 @@ +clc;
+disp("Example A.10")
+density=1000 // in kg/m^3
+Q=0.1/60 // flow rate in m^3/s
+mew=0.001 // viscosity in kg/ms
+tau=Q*density/(%pi*2)
+Re=4*tau/mew
+disp(Re,"Reynolds number is ")
+disp("It indicates viscous flow!")
diff --git a/1268/CH11/EX11.14/a_14.sce b/1268/CH11/EX11.14/a_14.sce new file mode 100644 index 000000000..1633e2c0f --- /dev/null +++ b/1268/CH11/EX11.14/a_14.sce @@ -0,0 +1,11 @@ +clc;
+disp("Example A.14")
+d=0.025 // diameter in m
+density= 1000 // in kg/m^3
+G=4*40/(60*%pi*d*d) // flow rate
+U=G/density // velocity in m/s
+mew=0.025
+Re=d*G/mew
+Q=40/60000
+delP=128*Q*mew/(%pi*d*d*d*d)
+disp(delP,"Pressure gradient for the liquid is ")
diff --git a/1268/CH11/EX11.15/a_15.sce b/1268/CH11/EX11.15/a_15.sce new file mode 100644 index 000000000..a74ddcbd9 --- /dev/null +++ b/1268/CH11/EX11.15/a_15.sce @@ -0,0 +1,5 @@ +clc;
+disp("Example A.15")
+dratio=4/5
+pratio=dratio^-5
+disp(pratio,"Pressure drop is ")
diff --git a/1268/CH11/EX11.17/a_17.sce b/1268/CH11/EX11.17/a_17.sce new file mode 100644 index 000000000..c1445ae25 --- /dev/null +++ b/1268/CH11/EX11.17/a_17.sce @@ -0,0 +1,10 @@ +clc;
+disp("Example A.17")
+d=14e-3 // diameter in m
+l=100 // length in m
+v=1e-5 // kinematic viscosity in kg/ms
+Re=2100
+U=Re*v/d
+h=32*v*l*U/(9.81*d*d)
+disp(h,"Head loss is ")
+
diff --git a/1268/CH11/EX11.19/a_19.sce b/1268/CH11/EX11.19/a_19.sce new file mode 100644 index 000000000..3b6fed749 --- /dev/null +++ b/1268/CH11/EX11.19/a_19.sce @@ -0,0 +1,10 @@ +clc;
+disp("Example A.19")
+n=100 // in rpm
+omega=2*%pi*n/60
+r=0.05 // radius in m
+u=r*omega // velocity in m/s
+gap=0.001 // in m
+mew=0.5 // in kg/ms
+tau=mew*u/gap
+disp(tau,"Shear stress is ")
diff --git a/1268/CH11/EX11.2/a_2.sce b/1268/CH11/EX11.2/a_2.sce new file mode 100644 index 000000000..315aaa478 --- /dev/null +++ b/1268/CH11/EX11.2/a_2.sce @@ -0,0 +1,7 @@ +clc;
+disp("Example A.2")
+x=1
+y=5
+z=x/y
+theta=asin(z)*180/%pi
+disp(theta, "The angle of inclination is ")
diff --git a/1268/CH11/EX11.3/a_3.sce b/1268/CH11/EX11.3/a_3.sce new file mode 100644 index 000000000..940fa7e52 --- /dev/null +++ b/1268/CH11/EX11.3/a_3.sce @@ -0,0 +1,22 @@ +clear;
+disp("Example A.3")
+d=0.097 // diameter in m
+gradp= 16 // pressure gradient in N/m
+density=1000 // in kg/m^3
+tau=gradp*d/4
+u=(tau/density)^0.5
+y=0.02 // in m
+v=1e-6 // kinematic viscosity iin m^2/s
+ydash=y*u/v
+
+if(ydash>30) then
+ udash=2.5*log(ydash)+5.5
+ ugrad=2.5/ydash
+ ratio=2*y/d // ratio of ydash/rdash=y/r
+ x=(1-ratio)/ugrad-1
+ disp(x,"Turbulent viscosity to molecular viscosity at 2m is ")
+end
+rdash=d*u/(2*v)
+ydash=rdash/2
+x=(0.5*ydash/2.5)-1
+disp(x,"Turbulent viscosity to molecular viscosity at the point of maximum viscosity is ")
diff --git a/1268/CH11/EX11.6/a_6.sce b/1268/CH11/EX11.6/a_6.sce new file mode 100644 index 000000000..61a321f0f --- /dev/null +++ b/1268/CH11/EX11.6/a_6.sce @@ -0,0 +1,10 @@ +clear;
+disp("Example A.6")
+d=1e-4 // diameter in m
+mew=1e-3 // viscosity in kg/ms
+densityp=1200 // of particle in kg/m^3
+density= 1000 // of water in kg/m^3
+t=0.256*densityp*d*d/mew
+U=densityp*d*d*9.81*(1-(density/densityp))/(18*mew)
+Re=d*U*density/mew
+disp(t,"Time is ")
diff --git a/1268/CH11/EX11.7/a_7.sce b/1268/CH11/EX11.7/a_7.sce new file mode 100644 index 000000000..34fd29f06 --- /dev/null +++ b/1268/CH11/EX11.7/a_7.sce @@ -0,0 +1,9 @@ +clc;
+disp("Example A.7")
+Q=700000/3600 // in m^3/s
+d=0.2 // diameter in m
+v=1.2e-5 // kinematic viscosity
+Re=4*Q/(%pi*d*v)
+ratio=0.0013
+f=0.3313/((log((ratio/3.7)+(5.74/(Re^0.9))))^2)
+disp(f,"Friction factor is ")
diff --git a/1268/CH11/EX11.8/a_8.sce b/1268/CH11/EX11.8/a_8.sce new file mode 100644 index 000000000..d8ebbf63e --- /dev/null +++ b/1268/CH11/EX11.8/a_8.sce @@ -0,0 +1,10 @@ +clc;
+disp("Example A.8")
+Re=90000
+f1=0.3313/((log(5.74/(Re^0.9)))^2)
+f2=0.079/(Re^0.25)
+if((f1-f2)<0.001) then
+ disp("Excellent agreement in friction factor is seen")
+end
+
+
diff --git a/1268/CH11/EX11.9/a_9.sce b/1268/CH11/EX11.9/a_9.sce new file mode 100644 index 000000000..8a9f92fa3 --- /dev/null +++ b/1268/CH11/EX11.9/a_9.sce @@ -0,0 +1,6 @@ +clc;
+disp("Example A.9")
+dratio=10
+uratio=1/dratio
+fratio=uratio*uratio*dratio*dratio
+disp(fratio,"The ratio of model force to prototype force is ")
|