6 files changed, 104 insertions, 0 deletions
diff --git a/1052/CH15/EX15.2/152.sce b/1052/CH15/EX15.2/152.sce
new file mode 100755
index 000000000..684693fd9
--- /dev/null
+++ b/1052/CH15/EX15.2/152.sce
@@ -0,0 +1,14 @@
+clc;
+//Example 15.2
+//page no 169
+printf(" Example 15.2 page no 169\n\n");
+//nitrogen gas is flowing in a duct,neglect compressibility effects
+T=293//temperature,k
+R=8314.4//gas constant
+k=1.4//for nitrogen
+M=28//molecular weight of nitrogen
+c=sqrt(k*R*T/M)//speed of sound in nitrogen 
+printf("\n speed of sound on nitrogen c=%f m/s",c);
+v=82//flow velocity 
+M_a=v/c//mach no.
+printf("\n mach no. M_a=%f ",M_a);
diff --git a/1052/CH15/EX15.3/153.sce b/1052/CH15/EX15.3/153.sce
new file mode 100755
index 000000000..0c3e3b5b0
--- /dev/null
+++ b/1052/CH15/EX15.3/153.sce
@@ -0,0 +1,21 @@
+clc;
+//Example 15.3
+//page no 170
+printf("Example 15.3 page no 170\n\n");
+//propane is flowing in a tube
+k=1.3//degree of freedom for propane
+T=290//temperature,k
+M=44//mol. weight 
+R=8314.4//gas constant
+c=sqrt((k*R*T)/M)//speed of sound in propane
+printf("\n speed of sound in propane c=%f m/s",c);
+v=43//average velocity
+M_a=v/c//mach no.
+printf("\n M_a mach no=%f ",M_a);
+//mach no is< 0.3,that's why flow is incompressible
+rho=6.39//density,kg/m^3
+meu=8e-6//viscosity ,m^2/s
+D=0.0254//inside diameter of tube
+R_e=D*rho*v/meu//reynolds no.
+printf("\n reynolds no R_e=%f ",R_e);
+//because R_e is >4000,flow is turbulent
diff --git a/1052/CH15/EX15.6/156.sce b/1052/CH15/EX15.6/156.sce
new file mode 100755
index 000000000..4e242676d
--- /dev/null
+++ b/1052/CH15/EX15.6/156.sce
@@ -0,0 +1,23 @@
+clc;
+//Example 15.6
+//page no 173
+printf("Example 15.6 page no 173\n\n");
+//methane is flowing through a horizontal steel pipe
+m_dot=10//mass flow rate, lb/s
+D=1//diameter of pipe,ft
+G=m_dot/((%pi/4)*D^2)//mass velocity flux
+P=89.7//inlet pressure
+T=530//temprature,k
+MW=16//mol. weight
+R=10.73//gas constant
+//applying eq 15.7
+rho=P*MW/(R*T)//density
+f=0.008//friction factor
+L=15840//length of pipe,ft
+g_c=32.2//gravitational constant
+P_drop=(2*f*L*(G^2))/(g_c*rho*D)//pressure drop
+P1=89.7//inlet pressure,psia
+P2=P1-(P_drop/144)
+P2=54.7//corrected value
+P_drop=P1-P2//updated value of P_drop
+printf("\n pressure drop P_drop=%f psia",P_drop);
diff --git a/1052/CH15/EX15.7/157.sce b/1052/CH15/EX15.7/157.sce
new file mode 100755
index 000000000..848c97cca
--- /dev/null
+++ b/1052/CH15/EX15.7/157.sce
@@ -0,0 +1,10 @@
+clc;
+//Example 15.7
+//page no 174
+printf("Example 15.7 page no 174\n\n");
+//refr to example 15.6
+D=1//diameter of pipe
+G=12.7//mass velocity flux
+meu=7.39e-6//viscosity,lb/ft.s
+R_e=(D*G)/(meu)//reynolds no
+printf("\n reynolds no R_e=%f ",R_e);
diff --git a/1052/CH15/EX15.8/158.sce b/1052/CH15/EX15.8/158.sce
new file mode 100755
index 000000000..80c8b4a3b
--- /dev/null
+++ b/1052/CH15/EX15.8/158.sce
@@ -0,0 +1,26 @@
+clc;
+//Example 15.8
+//page no 174
+printf("Example no page no 174\n\n");
+//air flowing through a steel pipe 
+P_1=2.7//pressure,atm
+T=288//temperature,k
+v=30//velocity at the entrance of the pipe ,m/s
+Mw=29//mol. weight of air
+V=22.4//standard volume
+T_s=273//st. temp
+P_s=1//st. pressure
+rho=(Mw*P_1*T_s)/(V*T*P_s)//density 
+printf("\ density rho =%f kg/m^3",rho);
+G=v*rho//mass veocity flux
+printf("\n G mass velocity flux =%f kg/m^2.s",G);
+f=0.004//friction factor
+D=0.085//diameter ,m
+L=65//length of pipe,m
+//gravitational constant
+P_2=P_1-2*f*L*G^2/(rho*D*101325)//pressure drop across the line
+//factor 101325 for atm
+printf("\n pressure drop P__2=%f atm",P_2);
+P_drop=P_1-P_2//pressure drop
+printf("\n P_drop pressure=%f atm",P_drop);
+
diff --git a/1052/CH15/EX15.9/159.sce b/1052/CH15/EX15.9/159.sce
new file mode 100755
index 000000000..1f8381188
--- /dev/null
+++ b/1052/CH15/EX15.9/159.sce
@@ -0,0 +1,10 @@
+clc;
+//Example 15.9
+//page no 175
+printf(" Example 15.9 page no 175\n\n");
+//refer to Example 15.9
+meu=1.74e-5//viscosity,kg/m.s
+D=0.085//diameter of pipe
+G=99.3//mass velocity flux
+R_e=D*G/meu//reynolds no.
+printf("\n reynolds no R_e=%f ",R_e);