2 files changed, 37 insertions, 0 deletions
diff --git a/2223/CH18/EX18.20/Ex18_20.sav b/2223/CH18/EX18.20/Ex18_20.sav
new file mode 100755
index 000000000..1783b434b
--- /dev/null
+++ b/2223/CH18/EX18.20/Ex18_20.sav
diff --git a/2223/CH18/EX18.20/Ex18_20.sce b/2223/CH18/EX18.20/Ex18_20.sce
new file mode 100755
index 000000000..fc5844b8b
--- /dev/null
+++ b/2223/CH18/EX18.20/Ex18_20.sce
@@ -0,0 +1,37 @@
+// scilab Code Exa 18.20 Calculation on an axial turbine cascade
+
+beta1=35; //  blade angle at entry
+beta2=55; // blade angle at exit
+i(1)=5; // incidence
+i(2)=10;
+i(3)=15;
+i(4)=20;
+delta=2.5; // deviation
+alpha2=beta2-delta; // air angle at exit
+a_r=2.5; // aspect ratio(h/l)
+
+n=4;
+for m=1:n
+//part(a)
+printf("\nfor incidence=%d\n",i(m))
+alpha1=beta1+i(m); // air angle at entry
+ep=alpha1+alpha2; // deflection angle
+disp("degree",ep,"(a)flow deflection is")
+p_c=0.505; //(s/l)
+
+//part(b) loss coefficient from Hawthorne relations
+
+z_p=0.025*(1+((ep/90)^2)); // Hawthorne's relation
+disp (z_p,"(b)the profile loss coefficient from Hawthorne relation is")
+z=(1+(3.2/a_r))*z_p; // the total cascade loss coefficient
+disp (z,"and the total loss coefficient is")
+Y=z; 
+
+// part(c)drag and lift coefficients
+alpham=atand((0.5*(tand(alpha2)-tand(alpha1))));
+C_D=p_c*Y*((cosd(alpham)^3)/(cosd(alpha2)^2));
+disp (C_D,"(c)the drag coefficient is")
+
+C_L=(2*p_c*(tand(alpha1)+tand(alpha2))*cosd(alpham))+(C_D*tand(alpham));
+disp (C_L,"and the Lift coefficient is")
+end