summaryrefslogtreecommitdiff
path: root/1727/CH7
diff options
context:
space:
mode:
Diffstat (limited to '1727/CH7')
-rwxr-xr-x1727/CH7/EX7.1/7_1.sce20
-rwxr-xr-x1727/CH7/EX7.10/7_10.sce15
-rwxr-xr-x1727/CH7/EX7.11/7_11.sce15
-rwxr-xr-x1727/CH7/EX7.12/7_12.sce17
-rwxr-xr-x1727/CH7/EX7.13/7_13.sce12
-rwxr-xr-x1727/CH7/EX7.14/7_14.sce21
-rwxr-xr-x1727/CH7/EX7.15/7_15.sce20
-rwxr-xr-x1727/CH7/EX7.16/7_16.sce11
-rwxr-xr-x1727/CH7/EX7.17/7_17.sce9
-rwxr-xr-x1727/CH7/EX7.18/7_18.sce12
-rwxr-xr-x1727/CH7/EX7.19/7_19.sce30
-rwxr-xr-x1727/CH7/EX7.2/7_2.sce14
-rwxr-xr-x1727/CH7/EX7.20/7_20.sce28
-rwxr-xr-x1727/CH7/EX7.21/7_21.sce19
-rwxr-xr-x1727/CH7/EX7.22/7_22.sce19
-rwxr-xr-x1727/CH7/EX7.3/7_3.sce20
-rwxr-xr-x1727/CH7/EX7.4/7_4.sce20
-rwxr-xr-x1727/CH7/EX7.5/7_5.sce17
-rwxr-xr-x1727/CH7/EX7.6/7_6.sce15
-rwxr-xr-x1727/CH7/EX7.7/7_7.sce16
-rwxr-xr-x1727/CH7/EX7.8/7_8.sce15
-rwxr-xr-x1727/CH7/EX7.9/7_9.sce21
22 files changed, 386 insertions, 0 deletions
diff --git a/1727/CH7/EX7.1/7_1.sce b/1727/CH7/EX7.1/7_1.sce
new file mode 100755
index 000000000..fb517f441
--- /dev/null
+++ b/1727/CH7/EX7.1/7_1.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+h1=4 //m
+muw=0.001 //Ns/m^2
+l=1.5 //m
+B=0.15/1000 //m
+len=11.2 //m
+//calculations
+P1=g*rho*h1
+V=P1*B^2 /(12*muw*l)
+A=B*len
+Q=A*V
+Q=7112.4
+tau= B/2 *(P1)/l
+//results
+printf("Average velocity through the crack = %.3f m/s",V)
+printf("\n rate of leakage = %.1f l/hr",Q)
+printf("\n Shear stress = %.3f N/m^2",tau)
diff --git a/1727/CH7/EX7.10/7_10.sce b/1727/CH7/EX7.10/7_10.sce
new file mode 100755
index 000000000..1bf2d36f2
--- /dev/null
+++ b/1727/CH7/EX7.10/7_10.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+l=1 //m
+b=0.3 //m
+Q=4.2 //m^3/s
+//calculations
+A=l*b
+R=A/(2*(l+b))
+d5=1.62/24.15
+d=d5^(1/5)
+Pr=2*(l+b)/(%pi*d)
+//results
+printf("The rectangular cross section will cost %.2f times that of a circular cross section",Pr)
diff --git a/1727/CH7/EX7.11/7_11.sce b/1727/CH7/EX7.11/7_11.sce
new file mode 100755
index 000000000..d1e4c4428
--- /dev/null
+++ b/1727/CH7/EX7.11/7_11.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+d1=2.5*10^-2 //m
+d2=7.2*10^-2 //m
+Q=100*10^-3 //m^3/hr
+//calculations
+V1=Q/(60*%pi/4*d1^2)
+V2=(d1/d2)^2 *V1
+dp= -(V2^2 -V1^2 + (V1-V2)^2)/(2*g)
+Pdiff=dp*g*rho
+//results
+printf("pressure difference = %.2f kN/m^2",Pdiff/1000)
+//The answers are a bit different due to rounding off error
diff --git a/1727/CH7/EX7.12/7_12.sce b/1727/CH7/EX7.12/7_12.sce
new file mode 100755
index 000000000..210fb8efb
--- /dev/null
+++ b/1727/CH7/EX7.12/7_12.sce
@@ -0,0 +1,17 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+d2=30/100 //cm
+d1=60/100 //cm
+Pu=105 //kN/m^2
+Pd=75 //kN/m^2
+Cc=0.65
+//calculations
+V22=(2*g/(1 - (d2/d1)^4 + (1/Cc -1)^2)) *(Pu-Pd)*10^3 /(rho*g)
+V2=sqrt(V22)
+Q=%pi/4 *V2 *d2^2
+hl=(1/Cc -1)^2 *V2^2 /(2*g)
+//results
+printf("Flow rate = %.3f m^3/s",Q)
+printf("\n Head loss = %.3f m",hl)
diff --git a/1727/CH7/EX7.13/7_13.sce b/1727/CH7/EX7.13/7_13.sce
new file mode 100755
index 000000000..6c38e1f49
--- /dev/null
+++ b/1727/CH7/EX7.13/7_13.sce
@@ -0,0 +1,12 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+d=9 //m
+dia=0.3 //m
+//calculations
+V302= 2*g*d/(0.5 + 20 + 2.53+101+0.66+41.47+2.07)
+V30=sqrt(V302)
+Q=%pi/4 *dia^2 *V30
+//results
+printf("Flow rate = %.3f m^3/s",Q)
diff --git a/1727/CH7/EX7.14/7_14.sce b/1727/CH7/EX7.14/7_14.sce
new file mode 100755
index 000000000..82b9151ce
--- /dev/null
+++ b/1727/CH7/EX7.14/7_14.sce
@@ -0,0 +1,21 @@
+clc
+clear
+//Initialization of variables
+h=6 //m
+rho=930 //kg/m^3
+Q=3/60 //m^3/s
+d=0.15 //m
+L=20 //m
+mu=0.006
+g=9.81 //m/s^2
+//calculations
+V=Q/(%pi/4 *d^2)
+RN=V*d*rho/mu
+f=0.316/RN^0.25
+hl=f*L/d *V^2 /(2*g)
+Hp=h+hl
+gam=rho*g
+W=gam*Q
+Power= W*Hp
+//results
+printf("Power required = %.3f kW",Power/1000)
diff --git a/1727/CH7/EX7.15/7_15.sce b/1727/CH7/EX7.15/7_15.sce
new file mode 100755
index 000000000..ad6dae2cd
--- /dev/null
+++ b/1727/CH7/EX7.15/7_15.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+d=0.02 //m
+d2=1.2 //m
+f=0.01
+L=250
+ken=0.5
+g=9.81
+h1=8 //m
+h2=4 //m
+//calculations
+V2=2*g/(1+ken+ f*L/d)
+V=sqrt(V2)
+Q=%pi/4 *d^2 *V
+function t=time(h)
+ t=-%pi/4 *d2^2 /Q /sqrt(h)
+endfunction
+ti=intg(h1,h2,time)
+//results
+printf("Time required = %d sec",ti)
diff --git a/1727/CH7/EX7.16/7_16.sce b/1727/CH7/EX7.16/7_16.sce
new file mode 100755
index 000000000..09d452925
--- /dev/null
+++ b/1727/CH7/EX7.16/7_16.sce
@@ -0,0 +1,11 @@
+clc
+//Initialization of variables
+d1=0.1 //m
+d2=0.05 //m
+l1=20 //m
+l2=20 //m
+f=0.02
+//calculations
+Kl=(f*l2/d2 *(d1/d2)^4 - f*l1/d1)
+//results
+printf("Loss coefficient = %d ",Kl)
diff --git a/1727/CH7/EX7.17/7_17.sce b/1727/CH7/EX7.17/7_17.sce
new file mode 100755
index 000000000..c70ce8ff6
--- /dev/null
+++ b/1727/CH7/EX7.17/7_17.sce
@@ -0,0 +1,9 @@
+clc
+clear
+//Initialization of variables
+g=9.81
+rati=1.265
+//calculations
+percent = (rati-1)*100
+//results
+printf("Increase in discharge = %.1f",percent)
diff --git a/1727/CH7/EX7.18/7_18.sce b/1727/CH7/EX7.18/7_18.sce
new file mode 100755
index 000000000..21c35fbc1
--- /dev/null
+++ b/1727/CH7/EX7.18/7_18.sce
@@ -0,0 +1,12 @@
+clc
+//Initialization of variables
+Q=0.6 //m^3/s
+l1=1200 //m
+l2=800 //m
+d1=0.3 //m
+//calculations
+V1=1.02 //m/s
+d5= d1*l2*4^2 *Q^2 /(l1*%pi^2 *V1^2)
+d=d5^(1/5)
+//results
+printf("diameter of the single pipe = %.2f m",d)
diff --git a/1727/CH7/EX7.19/7_19.sce b/1727/CH7/EX7.19/7_19.sce
new file mode 100755
index 000000000..2d0a661f2
--- /dev/null
+++ b/1727/CH7/EX7.19/7_19.sce
@@ -0,0 +1,30 @@
+clc
+//Initialization of variables
+g=9.81
+Q=0.18 //m^3/s
+d3=0.3//m
+f=0.032
+L3=360 //m
+z=25.5 //m
+z2=30 //m
+L2=450 //m
+d2=0.45//m
+L1=950 //m
+d1=0.45 //m
+zn=18 //m
+rho=1000
+//calculations
+V3=Q/(%pi/4 *d3^2)
+hl3=f*L3/d3 *(V3^2 /(2*g))
+Z2=z+hl3
+hl2=Z2-z2
+V2= sqrt(2*g*d2*hl2/(f*L2))
+Q2=%pi/4 *d2^2 *V2
+V1=V2+ (d3/d2)^2 *V3
+hl1=f*L1/d1*V1^2 /(2*g)
+Hp= hl1+ Z2-zn
+gam=rho*g
+P=gam*Hp
+//results
+printf("Discharge into the reservoir = %.3f m^3/s",Q2)
+printf("\n Pressure maintained by the pump = %.2f kN/m^2",P/1000)
diff --git a/1727/CH7/EX7.2/7_2.sce b/1727/CH7/EX7.2/7_2.sce
new file mode 100755
index 000000000..f2a7b653e
--- /dev/null
+++ b/1727/CH7/EX7.2/7_2.sce
@@ -0,0 +1,14 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=1200 //kg/m^3
+mu=0.005 //Ns/m^2
+d=0.006 //m
+Re=2000
+V=0.15 //m/s
+//calculations
+Vc=Re*mu/(d*rho)
+Vr=V/Vc
+T0=8*mu*V/d
+//results
+printf("Shear stress = %d N/m^2",T0)
diff --git a/1727/CH7/EX7.20/7_20.sce b/1727/CH7/EX7.20/7_20.sce
new file mode 100755
index 000000000..a40bba9a0
--- /dev/null
+++ b/1727/CH7/EX7.20/7_20.sce
@@ -0,0 +1,28 @@
+clc
+//Initialization of variables
+h=[1 2 1.9 1.96]
+z1=10 //m
+z2=5 //m
+z3=7.5 //m
+f=0.04
+l1=100 //m
+l2=50 //m
+l3=70 //m
+d1=0.1 //m
+d2=0.075 //m
+d3=0.06 //m
+g=9.81 //m/s^2
+//calculations
+Q1=sqrt(d1^5 *(%pi/4)^2 *2*g/(f*l1)) .*sqrt(z1-h)
+Q2=sqrt(d2^5 *(%pi/4)^2 *2*g/(f*l2)) .*sqrt(h+z2)
+Q3=sqrt(d3^5 *(%pi/4)^2 *2*g/(f*l3)) .*sqrt(h+z3)
+len=length(h)
+for i=1:len
+ Q=Q2(i)+Q3(i)
+ if (Q1(i) == Q) then
+ break;
+ end
+end
+printf("height h = %.2f m",h(i))
+printf("\nDischarge in BC Q2 = %.2f lps",Q2(i)*1000)
+printf("\nDischarge in BD Q3 = %.2f lps",Q3(i)*1000)
diff --git a/1727/CH7/EX7.21/7_21.sce b/1727/CH7/EX7.21/7_21.sce
new file mode 100755
index 000000000..9d8579ea5
--- /dev/null
+++ b/1727/CH7/EX7.21/7_21.sce
@@ -0,0 +1,19 @@
+clc
+//Initialization of variables
+e=0.8
+output=400 //kW
+H=150 //m
+rho=1000
+g=9.81
+f=0.028
+l=1250 //m
+//calculations
+gam=rho*g
+inpu=output/e
+Q=inpu*10^3 /(2/3 *gam*H)
+hl=1/3 *H
+d5= f*l*Q^2 /(2*g* %pi/4 * %pi/4 *hl)
+d=d5^(1/5)
+//results
+printf("Smallest diameter of pen stock = %d cm",d*100)
+
diff --git a/1727/CH7/EX7.22/7_22.sce b/1727/CH7/EX7.22/7_22.sce
new file mode 100755
index 000000000..5359f63c9
--- /dev/null
+++ b/1727/CH7/EX7.22/7_22.sce
@@ -0,0 +1,19 @@
+clc
+//Initialization of variables
+f=0.04
+H=30 //m
+l=200 //m
+d=0.075 //m
+g=9.81
+rho=1000
+gam=rho*g
+//calculations
+h=2/3 *H
+vj=sqrt(2*g*h)
+hl= 1/3 *H
+V= sqrt(hl*d*2*g/(f*l))
+dj= d*(sqrt(V/vj))
+Power= 2/3 *gam*%pi/4 *d^2 *V*H
+//results
+printf("Size of nozzle = %.1f cm",dj*100)
+printf("\n Max power = %.2f kW",Power/1000)
diff --git a/1727/CH7/EX7.3/7_3.sce b/1727/CH7/EX7.3/7_3.sce
new file mode 100755
index 000000000..15789dc40
--- /dev/null
+++ b/1727/CH7/EX7.3/7_3.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+Q=0.45/(60*1000) //m^3/s
+d=0.003 //m
+depth=0.95 //m
+alpha=2
+len=1.25 //m
+//calculations
+A=%pi/4 *d^2
+V=Q/A
+nu= (depth - alpha*V^2 /(2*g))*g*d^2 /(32*V*len)
+Re=V*d/nu
+//results
+if Re<2000 then
+ printf("Flow is laminar")
+else
+ printf("Flow is not laminar")
+end
diff --git a/1727/CH7/EX7.4/7_4.sce b/1727/CH7/EX7.4/7_4.sce
new file mode 100755
index 000000000..e02cac4b4
--- /dev/null
+++ b/1727/CH7/EX7.4/7_4.sce
@@ -0,0 +1,20 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=787 //kg/m^3
+Q=90*10^-3 //m^3/hr
+d=0.015 //m
+k=0.0045*10^-2 //m
+nu=1.6e-6
+l=5 //m
+//calculations
+V=Q/(60*%pi/4 *d^2)
+Rn=V*d/nu
+e=k/d
+disp("From moody diagram, f=0.028")
+f=0.028
+hl=f*l/d *V^2 /(2*g)
+Power=rho*g*Q/60 *hl
+//result
+printf("Head loss = %.2f m",hl)
+printf("\n power required =%.3f kW",Power/1000)
diff --git a/1727/CH7/EX7.5/7_5.sce b/1727/CH7/EX7.5/7_5.sce
new file mode 100755
index 000000000..e3565a172
--- /dev/null
+++ b/1727/CH7/EX7.5/7_5.sce
@@ -0,0 +1,17 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=870 //kg/m^3
+Q=2*10^-3 //m^3/s
+d=0.03 //m
+mu=5*10^-4
+l=50 //m
+//calculations
+V=Q/(%pi/4 *d^2)
+RN=rho*V*d/mu
+f=0.017
+hl=f*l/d *V^2/(2*g)
+Ploss=rho*g*hl
+//results
+printf("Loss of pressure = %.1f kN/m^2",Ploss/1000)
+//The answers are a bit different due to rounding off error in textbook
diff --git a/1727/CH7/EX7.6/7_6.sce b/1727/CH7/EX7.6/7_6.sce
new file mode 100755
index 000000000..716ee26f7
--- /dev/null
+++ b/1727/CH7/EX7.6/7_6.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=813 //kg/m^3
+Q=0.007 //m^3/hr
+d=0.01//m
+mu=0.002 //Ns/m^2
+l=30 //m
+//calculations
+V=Q/(60*%pi/4*d^2)
+RN=V*d*rho/mu
+f=0.316/RN^(0.25)
+h=(1+f*l/d)*V^2 /(2*g)
+//result
+printf("Height required = %.2f m",h)
diff --git a/1727/CH7/EX7.7/7_7.sce b/1727/CH7/EX7.7/7_7.sce
new file mode 100755
index 000000000..f8048975b
--- /dev/null
+++ b/1727/CH7/EX7.7/7_7.sce
@@ -0,0 +1,16 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+hl=0.02
+d=1.2 //m
+l=1 //m
+k=0.5 *10^-2 //m
+//calculations
+v2f=hl*(2*g*d)/l
+e=k/d
+f=0.028
+V=sqrt(v2f/f)
+Q=%pi/4 *d^2 *V
+//results
+printf("Rate of flow = %.2f m^3/s",Q)
diff --git a/1727/CH7/EX7.8/7_8.sce b/1727/CH7/EX7.8/7_8.sce
new file mode 100755
index 000000000..1e21a893a
--- /dev/null
+++ b/1727/CH7/EX7.8/7_8.sce
@@ -0,0 +1,15 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+e=0.03*10^-2 //m
+l=3000 //m
+Q=300*10^-3 //m^3/s
+nu=10^-5 //m^2/s
+hl=24 //m
+//calculations
+d5f= l*Q/(%pi/4) * Q/(%pi/4) /(hl*2*g)
+f=0.022
+d=(d5f*f)^(1/5)
+//results
+printf("Size of the required pipe = %d cm",d*100)
diff --git a/1727/CH7/EX7.9/7_9.sce b/1727/CH7/EX7.9/7_9.sce
new file mode 100755
index 000000000..55778e060
--- /dev/null
+++ b/1727/CH7/EX7.9/7_9.sce
@@ -0,0 +1,21 @@
+clc
+//Initialization of variables
+g=9.81 //m/s^2
+rho=10^3 //kg/m^3
+d=0.3 //m
+per=25/100
+Q=0.1 //m^3/s
+k0=0.025*10^-2 //m
+nu=0.000001
+year=10
+//calculations
+V=Q/(%pi/4 *d^2)
+RN=V*d/nu
+e1=k0/d
+f1=0.019
+f2=(1+per)*f1
+e2=0.002
+k2=e2*d
+rate = (k2-k0)*100/year
+//results
+printf("Rate of increase =%.4f cm/year",rate)
generated by cgit (git 2.34.1) at 2024-12-20 14:44:17 +0000