summaryrefslogtreecommitdiff
path: root/689/CH9/EX8.3/3.sce
diff options
context:
space:
mode:
Diffstat (limited to '689/CH9/EX8.3/3.sce')
-rw-r--r--689/CH9/EX8.3/3.sce32
1 files changed, 32 insertions, 0 deletions
diff --git a/689/CH9/EX8.3/3.sce b/689/CH9/EX8.3/3.sce
new file mode 100644
index 000000000..a67630ce9
--- /dev/null
+++ b/689/CH9/EX8.3/3.sce
@@ -0,0 +1,32 @@
+clc; funcprot(0);
+//Example 9.3 Horseshoe Vortex
+
+// Initialisation of variables
+T = 250;
+
+// Calculations
+function[y] =Velocity(R, theta1, theta2)
+ y = (T/(4*%pi*R))*(cos(theta1)-cos(theta2))
+endfunction
+//Considering Leg AB
+ theta1 = 0; //From figure
+ theta2 = %pi - atan(4/5); //From figure
+ R = 4; //From figure
+ Vt1 = Velocity(R, theta1, theta2);
+
+//Considering Leg BC
+ theta1 = atan(5/4); //From figure
+ theta2 = %pi - atan(5/6); //From figure
+ R = 5; //From figure
+ Vt2 = Velocity(R, theta1, theta2);
+
+//Considering Leg CD
+ theta1 = atan(5/6); //From figure
+ theta2 = %pi; //From figure
+ R = 6; //From figure
+ Vt3 = Velocity(R, theta1, theta2);
+
+Vt = Vt1 + Vt2 + Vt3;
+
+//Results
+disp(Vt,"Velocity at point P (ft/sec):",Vt3,"Velocity due to leg CD (ft/sec):",Vt2,"Velocity due to leg BC (ft/sec):",Vt1,"Velocity due to leg AB (ft/sec):");