initial commit / add all books

17 files changed, 327 insertions, 0 deletions

diff --git a/51/CH11/EX11.1/11_1.sce b/51/CH11/EX11.1/11_1.sce
new file mode 100755
index 000000000..c64990b38
--- /dev/null
+++ b/51/CH11/EX11.1/11_1.sce
@@ -0,0 +1,17 @@
+clc;

+clear;

+D=4;//in

+T1=540;//degree R

+p1=100;//psia

+T2=453;//degree R

+p2=18.4;//psia

+k=1.4;

+R=1716/32.174;//ft*lb/(lbm*(degree R))

+cv=R/(k-1);//ft*lb/(lbm*(degree R))

+udiff=cv*(T2-T1);//ft*lb/lbm; change in internal energy

+disp("ft*lb/lbm",udiff,"a)The change in internal energy between (1) and (2)=")

+cp=k*cv;//ft*lb/(lbm*(degree R))

+hdiff=cp*(T2-T1);//ft*lb/lbm; change in enthalpy

+disp("ft*lb/lbm",hdiff,"b)The  change in enthalpl energy between (1) and (2)=")

+ddiff=(1/R)*((p2*144/T2)-(p1*144/T1));//lbm/(ft^3); change in density

+disp("lbm/(ft^3)",ddiff,"The change in density betwenn (1) and (2)=")

diff --git a/51/CH11/EX11.11/11_11.sce b/51/CH11/EX11.11/11_11.sce
new file mode 100755
index 000000000..84dc3b739
--- /dev/null
+++ b/51/CH11/EX11.11/11_11.sce
@@ -0,0 +1,29 @@
+clc;

+clear;

+k=1.4;   

+T0=518.67;//degree R

+T1=514.55;//degree R

+p1=14.3;//psia

+R=53.3;//(ft*lb)/(lbm* degree R)

+cp=R*k/(k-1);//(ft*lb/(lbm* degree R))

+Tratio=T1/T0;

+Ma=(((1/Tratio)-1)/((k-1)/2))^0.5;

+x=(R*T1*k*32.2)^0.5;//ft/sec; where x=(R*T1*k)^0.5

+y=p1*144/(R*T1)*(Ma*x);//lbm/((ft^2)*sec); where y=d*V

+//for p=7 psia 

+p=7;//psia

+fn=poly([(-T0) 1 ((y*y/(2*cp*p*p*144*144/(R^2)))/32.2)],"T","c");

+r=roots(fn);

+T=r(1);//K

+sdif=(cp*log(T/T1))-(R*log(p/p1));//(ft*lb)/(lbm* degree R)

+disp("K",T,"The corrosponding value of temperature for Fanno for downstream pressure of 7psia=")

+disp("(ft*lb)/(lbm* degree R)",sdif,"The corrosponding value of entropy change for Fanno for downstream pressure of 7psia=")

+count=1;

+for i=1.4:0.1:7

+    root=roots(poly([(-T0) 1 ((y*y/(2*cp*i*i*144*144/(R^2)))/32.2)],"T","c"));

+    temp(count)=root(1);

+    s(count)=(cp*log(temp(count)/T1))-(R*log(i/p1));

+    count=count+1;

+end

+plot2d(s,temp)

+xtitle("T vs s-s1","s-s1, ((ft*lb)/(lbm* degree R))","T, Degree R")
\ No newline at end of file
diff --git a/51/CH11/EX11.11/11_11graph.jpg b/51/CH11/EX11.11/11_11graph.jpg
new file mode 100755
index 000000000..24fcddb4c
--- /dev/null
+++ b/51/CH11/EX11.11/11_11graph.jpg
diff --git a/51/CH11/EX11.12/11_12.sce b/51/CH11/EX11.12/11_12.sce
new file mode 100755
index 000000000..070326747
--- /dev/null
+++ b/51/CH11/EX11.12/11_12.sce
@@ -0,0 +1,42 @@
+clc;

+clear;

+T0=288;//K

+p0=101;//kPa(abs)

+l=2;//m

+D=0.1;//m

+f=0.02;

+k=1.4;

+x=f*l/D;

+Tratio=2/(k+1);//where Tratio is Tcritical/T0

+Tcritical=Tratio*T0;//K = T2

+Vcritical=(286.9*Tcritical*k)^0.5;//m/sec = V2

+//from value of x, the following are found

+Ma=0.63;

+Trat=1.1;//where Trat=T1/Tcritical

+Vrat=0.66;//where Vrat=V1/Vcritical

+prat=1.7;//where prat=p1/pcritical

+pratio=1.16;//where pratio=p0,1/p0critical

+//from value of Ma, the following are found

+Tfraction=0.93;//whereTfraction=T1/T0

+pfraction=0.76;//where pfraction=p1/p0,1

+dfraction=0.83;//where dfraction=d1/d0,1

+//hence,

+V1=Vrat*Vcritical;//m/sec

+d1=dfraction*(1.23);//kg/(m^3)

+m=d1*%pi*(D^2)*V1/4;//kg/sec

+T1=Tfraction*T0;//K

+p1=pfraction*p0;//kPa(abs)

+T01=T0;//K and T01=T02

+p01=p0;//kPa(abs)

+p2=(1/prat)*(pfraction)*p01;//kpa(abs)

+p02=(1/pratio)*p01;//kPa(abs)

+disp("K",Tcritical,"Critical temperature=")

+disp("m/sec",Vcritical,"Critical velocity=")

+disp("m/sec",V1,"Velocity at inlet=")

+disp("kg/sec",m,"Maximum mass flowrate=")

+disp("K",T1,"Temperature at inlet=")

+disp("kPa(abs)",p1,"Pressure at inlet=")

+disp("K",T01,"stagnation temperature at inlet and exit=")

+disp("kPa(abs)",p01,"The stagnation pressure at inlet=")

+disp("kPa(abs)",p2,"Pressure at exit=")

+disp("kPa(abs)",p02,"The stagnation pressure at exit=")

diff --git a/51/CH11/EX11.13/11_13.sce b/51/CH11/EX11.13/11_13.sce
new file mode 100755
index 000000000..d17ab92ea
--- /dev/null
+++ b/51/CH11/EX11.13/11_13.sce
@@ -0,0 +1,28 @@
+clc;

+clear;

+T0=288;//K

+p0=101;//kPa(abs)

+l=2;//m

+D=0.1;//m

+f=0.02;

+pd=45;//kPa(abs)

+f=0.02;

+k=1.4;

+lnew=(50/100)*l;

+x=lnew*f/D;

+//from this value of x, following are found

+Ma=0.7;

+prat=1.5;//where prat=p1/pcritical

+//from this value of Ma, following are found

+pratio=0.72;//where pratio=p1/p0

+dratio=0.79;//where dratio=d1/d0,1

+Vratio=0.73;//where Vratio=V1/Vcritical

+//hence,

+p2=(1/prat)*pratio*p0;//kPa(abs)

+pcritical=p2; 

+//we find that pd<pcritical

+d1=dratio*(1.23);//kg/(m^3)

+Vcritical=(286.9*Tcritical*k)^0.5;//m/sec = V2

+V1=Vratio*Vcritical;//m/sec

+m=d1*%pi*(D^2)*V1/4;//kg/sec

+disp("kg/sec",1.65,"is less than the flowrate for the longer tube =","kg/sec,",m,"The flowrate for the smaller tube=")
\ No newline at end of file
diff --git a/51/CH11/EX11.14/11_14.sce b/51/CH11/EX11.14/11_14.sce
new file mode 100755
index 000000000..4e6540df9
--- /dev/null
+++ b/51/CH11/EX11.14/11_14.sce
@@ -0,0 +1,24 @@
+clc;

+clear;

+T0=288;//K

+p0=101;//kPa(abs)

+l=2;//m

+D=0.1;//m

+f=0.02;

+pd=45;//kPa(abs)

+f=0.02;

+m=1.65;//kg/sec

+lnew=l/2;//m

+

+x=f*l/D;

+//from this value of x, Ma at exit is found as

+Ma=0.7;

+//and p2/pcritical is found as

+pratio=1.5;

+//and, from example 11.12,

+prat=1.7;//where prat=p1/pcritical

+pfraction=0.76;//where pfraction=p1/p0,1

+//Hence,

+p2=pratio*(1/prat)*pfraction*p0;//kPa(abs)

+disp(Ma,"The Mach number at the exit=")

+disp("kPa(abs)",p2,"The back pressure required=")
\ No newline at end of file
diff --git a/51/CH11/EX11.15/11_15.sce b/51/CH11/EX11.15/11_15.sce
new file mode 100755
index 000000000..de06b1d06
--- /dev/null
+++ b/51/CH11/EX11.15/11_15.sce
@@ -0,0 +1,31 @@
+clc;

+clear;

+k=1.4;   

+T0=518.67;//degree R

+T1=514.55;//degree R

+p1=14.3;//psia

+

+R=53.3;//(ft*lb)/(lbm*degree R)

+cp=R*k/(k-1);//(ft*lb/(lbm* degree R))

+Tratio=T1/T0;

+Ma=(((1/Tratio)-1)/((k-1)/2))^0.5;

+x=(R*T1*k*32.2)^0.5;//ft/sec; where x=(R*T1*k)^0.5

+y=p1*144/(R*T1)*(Ma*x);//lbm/((ft^2)*sec); where y=d*V

+z=R*T1/(p1*144);//(ft^3)/lbm

+c=(p1)+(y*y*z/(32.2*144));//psia; =constant

+//when downstream pressure p=13.5 psia

+p=13.5;//psia 

+a=(y^2)*R/(p*144*32.2*144);//(lb/(in^2))/degree R

+fn=poly([(p-c) a],"T","c");

+T=roots(fn);//degree R

+sdif=(cp*log(T/T1))-(R*log(p/p1));//ft*lb/(lbm*degree R)

+disp("degree R",T,"The corrosponding value of temperature for the downstream pressure of 13.5 psia=")

+disp("ft*lb/(lbm*degree R)",sdif,"The corrosponding value of change in entropy for the downstream pressure of 13.5 psia=")

+count=1;

+for i=1:0.1:13.5

+    temp(count)=roots(poly([(i-c) ((y^2)*R/(i*144*32.2*144))],"T","c"));

+    s(count)=(cp*log(temp(count)/T1))-(R*log(i/p1));

+    count=count+1;

+end

+plot2d(s,temp,rect=[100,500,400,3000])

+xtitle("T vs s-s1","s-s1, ((ft*lb)/(lbm* degree R))","T, Degree R")
\ No newline at end of file
diff --git a/51/CH11/EX11.15/11_15graph.jpg b/51/CH11/EX11.15/11_15graph.jpg
new file mode 100755
index 000000000..acf036b99
--- /dev/null
+++ b/51/CH11/EX11.15/11_15graph.jpg
diff --git a/51/CH11/EX11.18/11_18.sce b/51/CH11/EX11.18/11_18.sce
new file mode 100755
index 000000000..2aa3f5de2
--- /dev/null
+++ b/51/CH11/EX11.18/11_18.sce
@@ -0,0 +1,18 @@
+clc;

+clear;

+p=60;//psia

+T=1000;//degree R

+px=12;//psia

+k=1.4;

+R=53.3;//ft*lb/(lbm*degree R)

+pratio=p/px;

+//for this value of pratio, Max is calculated as

+Max=1.9;

+//using this value of Max,  Tx/T0,x is found as

+Tratio=0.59;

+//T=T0,x=T0,y

+Tx=Tratio*T;//degree R

+cx=(R*Tx*k)^0.5;//ft/sec

+Vx=1.87*cx*(32.2^0.5);//ft/sec

+disp(Max,"The Mach number for the flow=")

+disp("ft/sec",Vx,"The velocity of the flow=")

diff --git a/51/CH11/EX11.19/11_19.sce b/51/CH11/EX11.19/11_19.sce
new file mode 100755
index 000000000..6a03c6881
--- /dev/null
+++ b/51/CH11/EX11.19/11_19.sce
@@ -0,0 +1,31 @@
+clc;

+clear;

+x1=0.5;//m

+x2=0.3;//m

+Acritical=0.1;//m^2

+//at x1, Max1 is found as

+Max1=2.8;

+//and px/p0,x is found as

+pratio1=0.04;

+//For this value of Max, py/px is found as

+prat1=9;

+pfraction1=prat1*pratio1;//where pfraction=py/p0,x = pIII/p0,x

+//at x2, Max2 is found as

+Max2=2.14;

+//for this value of Max2, the following are found

+prat2=5.2;

+prat22=0.66;//where prat22=p0,y/p0,x

+May=0.56;

+//for this valur of May, Ay/Acritical is found as

+Aratio=1.24;

+Arat=(Acritical+(x1^2))/(Acritical+(x2^2));//where Aratio=A2/Ay

+Afraction=Aratio*Arat;//where Afraction=A2/Acritical

+A2=Acritical+(x1^2);//m^2

+Acritical1=A2/Afraction;//where Acritical1 critical area for the isentropic flow downstream of the shock

+//with the value of Afraction, the following are found

+Ma2=0.26;

+pfraction=0.95;//where pfraction=p2/p0,y

+//hence,

+pfrac=pfraction*prat22;//where pfrac=p2/p0,x

+disp(pfraction1,"The ratio of back pressure to inlet stagnation pressure that will result in a normal shock at the exit of the duct=")

+disp(pfrac,"The value of back pressure to inlet stagnation pressure required to position the shock at (x=0.3 m)=")

diff --git a/51/CH11/EX11.2/11_2.sce b/51/CH11/EX11.2/11_2.sce
new file mode 100755
index 000000000..bf0f47b03
--- /dev/null
+++ b/51/CH11/EX11.2/11_2.sce
@@ -0,0 +1,11 @@
+clc;

+clear;

+D=4;//in

+T1=540;//degree R

+p1=100;//psia

+T2=453;//degree R

+p2=18.4;//psia

+

+dratio=(p1/T1)*(T2/p2);

+sdif=(cv*(log(T2/T1)))+(R*(log(dratio)));//ft*lb/lbm*(degree R); change in entropy

+disp("ft*lb/lbm*(degree R)",sdif,"The change in entropy between (1) and (2)=")

diff --git a/51/CH11/EX11.3/11_3.sce b/51/CH11/EX11.3/11_3.sce
new file mode 100755
index 000000000..74fbd02d6
--- /dev/null
+++ b/51/CH11/EX11.3/11_3.sce
@@ -0,0 +1,7 @@
+clc;

+clear;

+T=0;//degree C

+R=286.9;//J/(kg*K)

+k=1.401;

+c=(R*(T+273.15)*k)^0.5;//m/s

+disp("m/sec",c,"The speed of sound for air at 0 degree C =")

diff --git a/51/CH11/EX11.4/11_4.sce b/51/CH11/EX11.4/11_4.sce
new file mode 100755
index 000000000..96c7ce83d
--- /dev/null
+++ b/51/CH11/EX11.4/11_4.sce
@@ -0,0 +1,19 @@
+clc;

+clear;

+z=1000;//m

+Ma=1.5;

+T=20;//degree C

+//alpha=atan(z/x), x=V*t,and Ma=(1/sin(alpha)); where alpha is the angle of the Mach cone

+//V=Ma*c

+c=343.3;//m/s found from the value of temperature

+V=Ma*c;//m/sec

+t=z/(Ma*c*tan(asin(1/Ma)));//sec

+disp("sec",t,"The number of seconds to wait after the plane passes over-head before it is heard=")

+Mach=0.01:0.01:4;

+count=1;

+for i=0.01:0.01:4

+    time(count)=z/(i*c*tan(asin(1/i)));

+    count=count+1;

+end

+plot2d(Mach,time,rect=[0,0,4,3])

+xtitle("t vs Ma","Ma","t, sec")

diff --git a/51/CH11/EX11.4/11_4graph.jpg b/51/CH11/EX11.4/11_4graph.jpg
new file mode 100755
index 000000000..86254ad74
--- /dev/null
+++ b/51/CH11/EX11.4/11_4graph.jpg
diff --git a/51/CH11/EX11.5/11_5.sce b/51/CH11/EX11.5/11_5.sce
new file mode 100755
index 000000000..7e55bc0d7
--- /dev/null
+++ b/51/CH11/EX11.5/11_5.sce
@@ -0,0 +1,23 @@
+clc;

+clear;

+A=1*(10^(-4));//m^2

+p1=80;//kPa(abs)

+p2=40;//kPa(abs)

+p0=101;//kPa(abs)

+pcritical=0.528*p0;//kPa(abs)

+k=1.4;

+//for (a) pth=p1>pcritical

+Math1=((((p0/p1)^((k-1)/k))-1)/((k-1)/2))^0.5;//Math=Mach number at throat

+//dth/d0=p1/p0; dth=density at throat

+dth1=(1.23)*(1/(1+(((k-1)/2)*(Math1^2))))^(1/(k-1));//kg/(m^3); density at throat

+Tth1=(288)*(1/(1+(((k-1)/2)*(Math1^2))));//K; temperature at throat

+Vth1=Math1*(286.9*Tth1*k)^0.5;//m/sec

+m1=dth1*A*Vth1;//kg/sec

+disp("kg/sec",m1,"a) The mass flowrate through the duct=")

+//for (b) pth=p2<pcritical, hence

+Math2=1;

+dth2=1.23*(1/(1+(((k-1)/2)*(Math2^2))))^(1/(k-1));//kg/(m^3); density at throat

+Tth2=(288)*(1/(1+(((k-1)/2)*(Math2^2))));//K; temperature at throat

+Vth2=Math2*(286.9*Tth2*k)^0.5;//m/sec

+m2=dth2*A*Vth2;//kg/sec

+disp("kg/sec",m2,"b) The mass flowrate through the duct=")
\ No newline at end of file
diff --git a/51/CH11/EX11.6/11_6.sce b/51/CH11/EX11.6/11_6.sce
new file mode 100755
index 000000000..8e2f90af9
--- /dev/null
+++ b/51/CH11/EX11.6/11_6.sce
@@ -0,0 +1,28 @@
+clc;

+clear;

+A=1*(10^(-4));//m^2

+p1=80;//kPa(abs)

+p2=40;//kPa(abs)

+

+p0=101;//kPa(abs)

+k=1.4;

+//for (a)

+pratio1=p1/p0;

+//for this value of p1/p0, 

+Math1=0.59;

+Tratio1=0.94;//=Tth/T0

+dratio1=0.85;//=dth/d0

+Tth1=Tratio1*(288);//K

+dth1=dratio1*(1.23);//kg/(m^3)

+Vth1=Math1*(286.9*Tth1*k)^0.5;//m/sec

+m1=(dth1*A*Vth1);//kg/sec

+disp("kg/sec",m1,"a)The mass flowrate=")

+//for (b)

+Math2=1;

+Tratio2=0.83;//=Tth/T0

+dratio2=0.64;//=dth/d0

+Tth2=Tratio2*(288);//K

+dth2=dratio2*(1.23);//kg/(m^3)

+Vth2=Math2*(286.9*Tth2*k)^0.5;//m/sec

+m2=(dth2*A*Vth2);//kg/sec

+disp("kg/sec",m2,"b)The mass flowrate=")
\ No newline at end of file
diff --git a/51/CH11/EX11.7/11_7.sce b/51/CH11/EX11.7/11_7.sce
new file mode 100755
index 000000000..c0efb19e8
--- /dev/null
+++ b/51/CH11/EX11.7/11_7.sce
@@ -0,0 +1,19 @@
+clc;

+clear;

+pratio=0.82;//ratio of static to stagnation pressure

+T=68;//degree F

+//for (a)

+//for the value of pratio given Ma is calculated as

+Ma1=0.54;

+k1=1.4;

+Tratio1=0.94;//T/T0

+T1=Tratio1*(T+460);// degree R

+V1=(Ma1*(53.3*T1*k1)^0.5)*(32.2^0.5);//ft/sec

+//for (b)

+k2=1.66;

+Ma2=((((1/pratio)^((k2-1)/k2))-1)/((k2-1)/2))^0.5;

+Tratio2=1/(1+(((k2-1)/2)*(Ma2^2)));//T/T0

+T2=Tratio2*(T+460);//degree R

+V2=(Ma2*(386*T2*k2)^0.5)*(32.2^0.5);//ft/sec

+disp("ft/sec",V1,"The flow velocity if fluid is air=")

+disp("ft/sec",V2,"The flow velocity if fluid is helium=")
\ No newline at end of file