10 files changed, 191 insertions, 0 deletions
diff --git a/2504/CH7/EX7.1/7_1.sce b/2504/CH7/EX7.1/7_1.sce
new file mode 100755
index 000000000..e353c6823
--- /dev/null
+++ b/2504/CH7/EX7.1/7_1.sce
@@ -0,0 +1,9 @@
+clc
+//initialisation of variables
+clear
+r= 1.5
+f= 0.025
+//CALCULATIONS
+r1= (2/f)*(r^2-1)
+//RESULTS
+printf ('ratio L/D2 = %.f',r1)
diff --git a/2504/CH7/EX7.10/7_10.sce b/2504/CH7/EX7.10/7_10.sce
new file mode 100755
index 000000000..79f22c8de
--- /dev/null
+++ b/2504/CH7/EX7.10/7_10.sce
@@ -0,0 +1,25 @@
+clc
+//initialisation of variables
+clear
+p2= 67.2 //lbf/in^2
+p1= 63 //lbf/in62
+r= 1.4
+n= 0.6
+T1= 870 //R
+ma1= 0.8 //ft/sec
+//CALCULATIONS
+dt= (p2/p1)^((r-1)/r)-1
+dt1= dt/n
+T2= T1*(1+dt1)
+Ts1= T1*(1+0.5*(r-1)*ma1^2)
+ps1= p1*(1+0.5*(r-1)*ma1^2)^(r/(r-1))
+ps2= p2*(Ts1/T2)^(r/(r-1))
+dp= ps1-ps2
+//RESULTS
+ printf ('dT = %.5f ',dt)
+ printf ('\n dT1 = %.5f ',dt1)
+ printf ('\n Temperature = %.f R',T2)
+ printf ('\n Temperature = %.1f R',Ts1)
+ printf ('\n Pressure = %.1f lbf/in^2',ps1)
+  printf ('\n Pressure = %.1f lbf/in^2',ps2)
+   printf ('\n pressure difference = %.1f lbf/in^2',dp)
diff --git a/2504/CH7/EX7.11/7_11.sce b/2504/CH7/EX7.11/7_11.sce
new file mode 100755
index 000000000..eef7c66cd
--- /dev/null
+++ b/2504/CH7/EX7.11/7_11.sce
@@ -0,0 +1,29 @@
+clc
+//initialisation of variables
+clear
+r= 1.4
+ma3= 3 //ft/sec
+ps= 80 //lbf/ft^2
+Ts= 840 //R
+r1= 53.3 //ft lbm/ft^3
+A3= 2 //in^2
+g= 32.2 //ft/sec^2
+ma1= 1.6
+//CALCULATIONS
+R= (1+(r-1)*0.5*ma3^2)^(r/(r-1))
+p3= ps/R
+R1= 1+0.5*(r-1)*ma3^2
+T3= Ts/R1
+w3= p3*144/(r1*T3)
+V3= ma3*sqrt(r*r1*g*T3)
+m= w3*V3*A3/144
+ra= ((r+1)/(2*r*ma1^2-(r-1)))^(1/(r-1))*(0.5*(r+1)*ma1^2/(1+0.5*(r-1)*ma1^2))^(r/(r-1))
+ps2= ps*ra
+dp= ps-ps2
+//RESULTS
+ printf ('outlet pressure = %.2f lbf/in^2',p3)
+  printf ('\n outlet temperature = %.f R',T3)
+  printf ('\n mass flow rate = %.3f lbm/sec',m)
+  printf ('\n mass flow rate = %.3f lbm/sec',m)
+   printf ('\n ps2 = %.1f lbf/in^2',ps2)
+   printf ('\n pressure difference = %.1f lbf/in^2',dp)
diff --git a/2504/CH7/EX7.2/7_2.sce b/2504/CH7/EX7.2/7_2.sce
new file mode 100755
index 000000000..d10e67017
--- /dev/null
+++ b/2504/CH7/EX7.2/7_2.sce
@@ -0,0 +1,22 @@
+clc
+//initialisation of variables
+clear
+a= 6 //degrees
+r= 1.5
+l= 100 //ft
+f= 0.025
+K= 0.15
+//CALCULATIONS
+R= r^4-1
+R1= cotd(a/2)*(1-(1/r))
+p1= f*l
+p2= 2.5*(l-p1)/l
+p3= (1-r^2)^2
+p4= K*p3
+pt= p4+p2
+//RESULTS
+printf ('lowest ratio = %.2f',R)
+printf ('\n contribtuion of friction in pipe = %.1f lbf/ft^2',p1)
+printf ('\n contribtuion of diffuser in pipe = %.3f lbf/ft^2',p2)
+printf ('\n stagnant pressure drop = %.3f lbf/ft^2',p3)
+printf ('\n contribtuion of friction in pipe after reduction = %.3f lbf/ft^2',pt)
diff --git a/2504/CH7/EX7.3/7_3.sce b/2504/CH7/EX7.3/7_3.sce
new file mode 100755
index 000000000..69cf756f9
--- /dev/null
+++ b/2504/CH7/EX7.3/7_3.sce
@@ -0,0 +1,21 @@
+clc
+//initialisation of variables
+clear
+d= 4 //in
+q= 0.5 //ft^3/sec
+w= 62.4 //lb/ft^3
+u= 2.7*10^-5 //lbf sec/ft^2
+e= 0.0005 //ft
+g= 32.1 //ft/sec^2
+f= 0.0235
+lt= 400 //ft
+//CALCULATIONS
+V= 4*q/(%pi*(d/12)^2)
+Re= w*V*(d/12)/(u*g)
+r= e/(d/12)
+dz= (V^2/(2*g))*(1.7+f*lt/(d/12))
+//RESULTS
+printf ('mean flow velocity = %.2f ft/sec',V)
+printf ('\n Reynolds number = %.2e',Re)
+printf ('\n Relative roughness = %.4f',r)
+printf ('\n difference in the levels of water = %.1f ft',dz)
diff --git a/2504/CH7/EX7.4/7_4.sce b/2504/CH7/EX7.4/7_4.sce
new file mode 100755
index 000000000..8f59b08ef
--- /dev/null
+++ b/2504/CH7/EX7.4/7_4.sce
@@ -0,0 +1,9 @@
+clc
+//initialisation of variables
+clear
+d= 4 //in
+v= 6.64 //ft/sec
+//CALCULATIONS
+Q= %pi*0.25*(d/12)^2*v
+//RESULTS
+printf ('Flow rate= %.3f ft^3/sec',Q)
diff --git a/2504/CH7/EX7.5/7_5.sce b/2504/CH7/EX7.5/7_5.sce
new file mode 100755
index 000000000..eafb56f6d
--- /dev/null
+++ b/2504/CH7/EX7.5/7_5.sce
@@ -0,0 +1,9 @@
+clc
+//initialisation of variables
+clear
+d= 0.366 //ft
+i= 12
+//CALCULATIONS
+pd= d*i
+//RESULTS
+printf ('Required pipe diameter = %.2f in',pd)
diff --git a/2504/CH7/EX7.6/7_6.sce b/2504/CH7/EX7.6/7_6.sce
new file mode 100755
index 000000000..9d2ed5b51
--- /dev/null
+++ b/2504/CH7/EX7.6/7_6.sce
@@ -0,0 +1,20 @@
+clc
+//initialisation of variables
+clear
+Ps1= 1050 //lbf/ft^2
+fr= 10.7
+p= 36.6 //lbf/ft^2
+p1= 195 //lbf/ft^2
+fr1= 16
+fr2= 1.8
+//CALCULATIONS
+p2= fr*p
+dp= Ps1-p2
+lc= dp/p
+sp= Ps1+p1-p*(fr1+fr2)
+lc1= sp/p
+//RESULTS
+printf ('Pressure = %.f lbf/ft^2',p1)
+printf ('\n pressure difference = %.f lbf/ft^2',dp)
+printf ('\n Loss coefficient = %.f ',lc)
+printf ('\n Loss coefficient = %.1f ',lc1)
diff --git a/2504/CH7/EX7.7/7_7.sce b/2504/CH7/EX7.7/7_7.sce
new file mode 100755
index 000000000..299851d2f
--- /dev/null
+++ b/2504/CH7/EX7.7/7_7.sce
@@ -0,0 +1,26 @@
+clc
+//initialisation of variables
+clear
+p1= 50 //lbf/in^2
+R= 96.3 //ft lbf/lbm R
+T= 80 //F
+p2= 20 //lbf/in^2
+r= 1.31
+u= 2.34*10^-7 //lbf sec/ft^2
+e= 0.00005 //ft
+m= 5*10^4 //lbm/sec
+d= 1.5 //ft
+g= 32.2 //ft/sec^2
+f= 0.113
+//CALCULATIONS
+w1= p1*144/(R*(460+T))
+V1= 4*(m/3600)/(%pi*w1*d^2)
+Ma1= V1/(r*R*g*(460+T))^0.5
+Re= w1*V1*d/(u*g)
+dx= (((1/(r*Ma1^2))*10*(1-(p2/p1)^2))+log(p2/p1))*d/f
+//RESULTS
+printf ('density = %.3f lbm/ft^3',w1)
+printf ('\n mean flow velocity = %.1f ft/sec',V1)
+printf ('\n Match number = %.4f ',Ma1)
+printf ('\n Reynolds number = %.2e ',Re)
+printf ('\n Length of pipe = %.2e ft',dx)
diff --git a/2504/CH7/EX7.9/7_9.sce b/2504/CH7/EX7.9/7_9.sce
new file mode 100755
index 000000000..9dfa64453
--- /dev/null
+++ b/2504/CH7/EX7.9/7_9.sce
@@ -0,0 +1,21 @@
+clc
+//initialisation of variables
+clear
+r= 1.4
+R= 53.3 //ft lbf/lbm R
+g= 32.2 //ft/sec^2
+T1= 410 //R
+v= 2500 //ft/sec
+P1= 628 //lbf/in^2
+//CALCULATIONS
+v1= sqrt(r*g*R*T1)
+Ma1= v/v1
+Ts1= T1*(1+0.5*(r-1)*Ma1^2)
+Ps1= P1*(1+0.5*(r-1)*Ma1^2)^(r/(r-1))
+Ps2= Ps1*((r+1)/(2*r*Ma1^2-r+1))^(1/(r-1))*(0.5*(r+1)*Ma1^2/(1+0.5*(r-1)*Ma1^2))^(r/(r-1))
+//RESULTS
+printf ('acoustic velocity = %.f ft/sec',v1)
+printf ('\n Match number = %.2f ',Ma1)
+printf ('\n Stagnition temperature = %.f R',Ts1)
+printf ('\n Stagnition pressure = %.f lbf/ft^2',Ps1)
+printf ('\n Stagnition pressure = %.f lbf/ft^2',Ps2)