5 files changed, 117 insertions, 0 deletions
diff --git a/1187/CH9/EX9.2/2.sce b/1187/CH9/EX9.2/2.sce
new file mode 100755
index 000000000..821fd2317
--- /dev/null
+++ b/1187/CH9/EX9.2/2.sce
@@ -0,0 +1,18 @@
+clc
+
+// p_a-p_b=-1/2*rho*C^2*(1/R_A^2-1/R_B^2)
+
+rho_w=1000; // kg/m^3
+g=9.81; // m/s^2
+h=0.0115; // m
+rho=1.22; // kg/m^3
+R_A=0.4; // m
+R_B=0.2; // m
+
+C=sqrt(rho_w*g*h*2/(rho*(1/R_B^2-1/R_A^2)));
+
+m=rho*C*R_B*integrate('1/R','R', R_B, R_A);
+
+disp("Mass flow rate =")
+disp(m)
+disp("kg/s")
diff --git a/1187/CH9/EX9.3/3.sce b/1187/CH9/EX9.3/3.sce
new file mode 100755
index 000000000..a22fcfd41
--- /dev/null
+++ b/1187/CH9/EX9.3/3.sce
@@ -0,0 +1,25 @@
+clc
+
+// p=1/2*rho*w^2*R^2 + C
+
+
+// At z=0
+rho=900; // kg/m^3
+g=9.81; // m/s^2
+h=0.6; // m
+
+C=rho*g*h;
+
+// p = -rho*K^2/(2*R^2)+D
+// From this we get, D = 9*w^2 + C
+
+// At z = 0
+// p = D - rho*K^2/2/R^2;
+p_max=150000; // Pa
+
+// From the above equation we obtain,
+w=135.6; // rad/s
+
+disp("The maximum speed at which the paddles may rotate about their vertical axis =")
+disp(w)
+disp("rad/s")
diff --git a/1187/CH9/EX9.4/4.sce b/1187/CH9/EX9.4/4.sce
new file mode 100755
index 000000000..9eef6a109
--- /dev/null
+++ b/1187/CH9/EX9.4/4.sce
@@ -0,0 +1,27 @@
+clc
+
+U=40; // m/s
+h=0.01; // m
+
+m=2*U*h;
+disp("(a) the strength of the line source =")
+disp(m)
+disp("m^2/s")
+
+s = m/(2*%pi*U);
+disp("(b) the distance s the line source is located behind the leading edge of the step =")
+disp(s*1000)
+disp("mm")
+
+x=0; // m
+y=0.005; // m
+
+u=U + m/(2*%pi)*(x/(x^2+y^2));
+v=m/(2*%pi)*(y/(x^2+y^2));
+disp("Horizontal component =")
+disp(u)
+disp("m/s")
+
+disp("Vertical Component =")
+disp(v)
+disp("m/s")
diff --git a/1187/CH9/EX9.5/5.sce b/1187/CH9/EX9.5/5.sce
new file mode 100755
index 000000000..8d46a31a9
--- /dev/null
+++ b/1187/CH9/EX9.5/5.sce
@@ -0,0 +1,13 @@
+clc
+
+b=0.0375; // m
+t=0.0625; // m
+U=5; // m/s
+
+m=2*%pi*U*t/atan(2*b*t/(t^2-b^2));
+
+L=2*b*(1+m/(%pi*U*b))^(1/2);
+
+disp("L =")
+disp(L)
+disp("m")
diff --git a/1187/CH9/EX9.7/7.sce b/1187/CH9/EX9.7/7.sce
new file mode 100755
index 000000000..1b3faec04
--- /dev/null
+++ b/1187/CH9/EX9.7/7.sce
@@ -0,0 +1,34 @@
+clc
+
+l1=10; // m
+r1=2; // m
+C_D1=0.0588;
+theta1=6.5; // degrees
+
+AR1=l1/r1; // Aspect ratio
+
+C_L=0.914;
+
+C_D2=C_L^2/(%pi*AR1);
+theta2=atand(C_L/(%pi*AR1))
+
+C_D3=C_D1-C_D2;
+theta3=theta1-theta2;
+
+AR2=8;
+
+C_Di=C_L^2/(%pi*AR2);
+C_D=C_Di+C_D3;
+
+theta4=atand(C_L/(%pi*AR2));
+theta=theta4+theta3;
+
+disp("Lift coefficient =")
+disp(C_L)
+
+disp("Drag coefficient =")
+disp(C_D)
+
+disp("Effective angle of attack =")
+disp(theta)
+disp("degrees")