1 files changed, 30 insertions, 0 deletions
diff --git a/689/CH15/EX15.6/6.sce b/689/CH15/EX15.6/6.sce
new file mode 100644
index 000000000..7362b8662
--- /dev/null
+++ b/689/CH15/EX15.6/6.sce
@@ -0,0 +1,30 @@
+clc; funcprot(0);
+//Example 15.6 Thrust Horsepower of a constant speed propeller at Sea level
+// Initialisation of variables
+V = 125*1.467;      // Velocity in ft/sec
+n = 1800/60;        // Rotaion per second
+D = 9.5;              // Diameter in ft
+P = 600*550;
+rho = 0.002378;
+
+// Calculations
+//For Beta = 23 degree
+Cp = P/(rho*n^3*D^5);
+ThrustCoeff = rho*n^2*D^4;
+VeloCoeff = n*60*D/88;
+
+VBynD = [1.03 1 0.9 0.8 0.7 0.6];
+Ct = [0.055 0.056 0.062 0.068 0.076 0.086];              // From figure 15.4 for corrospondiong values of V/nD at Cp = 0.066
+T = Ct*ThrustCoeff;
+V = VBynD*VeloCoeff;
+
+Result = zeros(6,6);
+Result(:,1) = VBynD';
+Result(:,2) = Ct';
+Result(:,3) = Result(:,2)/Cp;
+Result(:,4) = diag(diag(Result(:,3))*diag(Result(:,1)));        // Using equation 15.5.5
+Result(:,5) = Result(:,4)*P/550;
+Result(:,6) = V';
+//Results 
+
+disp(Result,"!! V/nD   Ct        Ct/Cp        Efficiency   thp          V(mph)  !!") ;