initial commit / add all books

12 files changed, 151 insertions, 0 deletions

diff --git a/51/CH3/EX3.10/3_10.sce b/51/CH3/EX3.10/3_10.sce
new file mode 100755
index 000000000..a8c5e298e
--- /dev/null
+++ b/51/CH3/EX3.10/3_10.sce
@@ -0,0 +1,16 @@
+clc;

+clear;

+T=60;//degree farenheit

+z1=5;//ft

+atmp=14.7;//psia

+//applying bernoulli equation at points 1,2 and 3

+z3=-5;//ft

+v1=0;//large tank

+p1=0;//open tank

+p3=0;//open jet

+//applying continuity equation A2*v2=A3*v3; A2=A3; so v2=v3

+v3=(2*32.2*(z1-z3))^0.5;

+//vapor pressure of water at 60 degree farenheit = p2=0.256 psia

+p2=0.256;

+z2=z1-((((p2-atmp)*144)+(0.5*1.94*v3^2))/62.4);

+disp("ft",z2,"The maximum height over which the water can be siphoned without cavitation occuring=")
\ No newline at end of file
diff --git a/51/CH3/EX3.11/3_11.sce b/51/CH3/EX3.11/3_11.sce
new file mode 100755
index 000000000..c36f7b19d
--- /dev/null
+++ b/51/CH3/EX3.11/3_11.sce
@@ -0,0 +1,15 @@
+clc;

+clear;

+sg=0.85;

+Q1=0.005;//m^3/s

+Q2=0.05;//m^3/s

+dia1=0.1;//m

+dia2=0.06;//m

+

+//A2/A1=dia2/dia1

+d=sg*1000;

+Arat=(dia2/dia1)^2;

+A2=%pi/4*(dia2^2);

+pdiffs=(Q1^2)*d*(1-(Arat^2))/(2*1000*(A2^2));

+pdiffl=(Q2^2)*d*(1-(Arat^2))/(2*1000*(A2^2));

+disp("kPa",pdiffl,"to","kPa",pdiffs,"kPa","The pressure difference ranges from =")
\ No newline at end of file
diff --git a/51/CH3/EX3.12/3_12.sce b/51/CH3/EX3.12/3_12.sce
new file mode 100755
index 000000000..da4255811
--- /dev/null
+++ b/51/CH3/EX3.12/3_12.sce
@@ -0,0 +1,21 @@
+clc;

+clear;

+z1=5;//m

+a=0.8;//m

+b=6;//m

+Cc=0.61;//since a/z1=ratio=0.16<0.2; Cc= contracction coefficient

+z2=Cc*a;

+//Q/b=flowrate

+flowrate=z2*((2*9.81*(z1-z2))/(1-((z2/z1)^2)))^0.5;

+//considering z1>>z2 and neglecting kinetic energy of the upstream fluid

+flowrate1=z2*(2*9.81*z1)^0.5;

+disp("m^2/s",flowrate,"The flowrate per unit width=")

+disp("m^2/s",flowrate1,"The flowrate per unit width when we consider z1>>z2=")

+count=1;

+j=5:15;

+for i=5:15

+    fr(count)=z2*((2*9.81*(i-z2))/(1-((z2/i)^2)))^0.5;

+    count=count+1;

+end

+plot2d(j,fr,rect=[0,0,15,9])

+xtitle("Q/b vs z1","z1,m","Q/b, m^2/s")
\ No newline at end of file
diff --git a/51/CH3/EX3.12/3_12graph.jpg b/51/CH3/EX3.12/3_12graph.jpg
new file mode 100755
index 000000000..900df5d77
--- /dev/null
+++ b/51/CH3/EX3.12/3_12graph.jpg
diff --git a/51/CH3/EX3.13/3_13.sce b/51/CH3/EX3.13/3_13.sce
new file mode 100755
index 000000000..f007c3f95
--- /dev/null
+++ b/51/CH3/EX3.13/3_13.sce
@@ -0,0 +1,7 @@
+clc;

+clear;

+//Q=A*V=(H^2)*tan(theta/2)*(C2*(2*g*H)^0.5)

+//Q3H0/QH0=(3H0)^2.5/(H0)^2.5=3^2.5

+Qrat=3^2.5;

+disp("The flowrate is proportional to H^2.5")

+disp("times.",Qrat,"When depth is increased from H0 to 3H0 Q increases ")

diff --git a/51/CH3/EX3.15/3_15.sce b/51/CH3/EX3.15/3_15.sce
new file mode 100755
index 000000000..c2b13d7b1
--- /dev/null
+++ b/51/CH3/EX3.15/3_15.sce
@@ -0,0 +1,16 @@
+clc;

+clear;

+h=10;//Km

+//air is in a standard atmosphere

+p1=26.5;//kPa

+T1=-49.9;//degree celcius

+d=0.414;//Kg/m^3

+k=1.4;

+Ma1=0.82;//Mach

+//for incompressible flow,

+pdiff=(k*Ma1^2)/2*p1;

+//for compressible isentropic flow, 

+pdiff1=((1+((k-1)/2)*(Ma1^2))^(k/(k-1))-1)*p1;

+disp("Stagnation pressure on leading edge on the wing of the Boeing:")

+disp("kPa",pdiff,"flow is imcompressible =")

+disp("kPa",pdiff1,"flow is compressible and isentropic =")
\ No newline at end of file
diff --git a/51/CH3/EX3.17/3_17.sce b/51/CH3/EX3.17/3_17.sce
new file mode 100755
index 000000000..c60937030
--- /dev/null
+++ b/51/CH3/EX3.17/3_17.sce
@@ -0,0 +1,9 @@
+clc;

+clear;

+V=5;//m/s

+sg=1.03;

+h=50;//m

+//since static pressure is greater than stagnation pressure, Bernoulli's equation is incorrect

+//p2=(d*(V1^2)/2)+(d*g*h) ; V1=V

+p2=(((sg*1000)*(V^2)/2) + (sg*1000*9.81*h))/1000;//kPa

+disp("kPa",p2,"The pressure at stagnation point 2 =")
\ No newline at end of file
diff --git a/51/CH3/EX3.6/3_6.sce b/51/CH3/EX3.6/3_6.sce
new file mode 100755
index 000000000..7fcef1d6b
--- /dev/null
+++ b/51/CH3/EX3.6/3_6.sce
@@ -0,0 +1,26 @@
+clc;

+clear;

+v1=100;//mi/hr

+ht=10000;//ft

+//from standard table for static pressure at an altitude

+p1=1456//lb/ft^2(abs)

+P1=1456*0.006947;//psi

+d=0.001756;//slugs/ft^3

+//1 mi/hr = 1.467 ft/s 

+p2=p1+(d*(v1*1.467)^2/2);//lb/ft^3

+//in terms of gage pressure p2g

+p2g=p2-p1;//lb/ft^2

+//1lb/ft^2 = 0.006947 psi

+P2=p2*0.006947;//psi

+P2g=p2g*0.006947;//psi

+//pressure difference indicated by the pitot tube = pdiff

+pdiff=P2-P1;//psi

+disp("psi",P1,"Pressure at point 1 =")

+disp("psi",P2g,"Pressure at point 2 in terms of gage pressure=")

+disp("psi",pdiff,"pressure difference indicated by the pitot static tube=")

+v1=0:1:600;

+for i=0:600

+   prat(i+1)=p1/(p1+(d*(i*1.467)^2/2));

+end

+plot2d(v1,prat,rect=[0,0,600,1]);

+xtitle("v1 vs p1/p2","v1, mph","p1/p2")
\ No newline at end of file
diff --git a/51/CH3/EX3.6/3_6graph.jpg b/51/CH3/EX3.6/3_6graph.jpg
new file mode 100755
index 000000000..9dfb8ebcd
--- /dev/null
+++ b/51/CH3/EX3.6/3_6graph.jpg
diff --git a/51/CH3/EX3.7/3_7.sce b/51/CH3/EX3.7/3_7.sce
new file mode 100755
index 000000000..a07ee461c
--- /dev/null
+++ b/51/CH3/EX3.7/3_7.sce
@@ -0,0 +1,26 @@
+clc;

+clear;

+dia=0.1;//m

+dia1=1.0;//m

+h=2.0;//m

+//bernoulli's equation: p1+(0.5*d*V1^2)+(sw*z1)= p2+(0.5*d*V2^2)+(sw*z2)

+//assuming p1=p2=0, and z1=h and z2=0

+//(0.5*d*V1^2)+(g*h)= (0.5*d*V2^2)

+//assuming steady flow Q1=Q2, Q=A*V. hence, A1*V1=A2*V2

+//V1=((dia/dia1)^2)*V2

+//hence V2=((2*g*h)/(1-(dia/dia1)^4))^0.5

+V2=((2*9.81*h)/(1-(dia/dia1)^4))^0.5;

+Q=(%pi/4*(dia)^2)*V2;

+disp("m^3/sec",Q,"The flow rate needed is=")

+//let Q0 be the flow rate when v1=0, i.e. dia>>dia

+//Q0=(2*g*h)^0.5 and Qrat=Q/Q0

+count=1;

+i=0:0.05:0.8;

+

+for k=0.00:0.05:0.80

+    Qrat(count)=1/((1-(k^4))^0.5);

+    count=count+1;

+end

+

+plot2d(i,Qrat,rect=[0,1,0.8,1.1])

+xtitle("d/D vs Q/Q0","d/D","Q/Q0")

diff --git a/51/CH3/EX3.7/3_7graph.jpg b/51/CH3/EX3.7/3_7graph.jpg
new file mode 100755
index 000000000..481b10e8f
--- /dev/null
+++ b/51/CH3/EX3.7/3_7graph.jpg
diff --git a/51/CH3/EX3.8/3_8.sce b/51/CH3/EX3.8/3_8.sce
new file mode 100755
index 000000000..ea14437a5
--- /dev/null
+++ b/51/CH3/EX3.8/3_8.sce
@@ -0,0 +1,15 @@
+clc;

+clear;

+dia=0.03;//m

+dia1=0.01;//m

+p=3;//kPa(gage)

+//density of air d is found using standard temp and pressure conditions

+d=(p+101)*1000/((286.9)*(15+273));

+//applying Bernoulli's equation at points 1,2 and 3; p=p1

+v3=((2*p*1000)/d)^0.5;

+Q=%pi/4*(dia1^2)*v3;

+//by continuity equation, A2*v2=A3*v3

+v2=((dia1/dia)^2)*v3;

+p2=(p*1000)-(0.5*d*(v2^2));

+disp("m^3/s",Q,"Flowrate =")

+disp("N/m^2",p2,"Pressure in the hose=")
\ No newline at end of file