1 files changed, 41 insertions, 0 deletions
diff --git a/3511/CH11/EX11.5/Ex11_5.sce b/3511/CH11/EX11.5/Ex11_5.sce
new file mode 100644
index 000000000..893986aeb
--- /dev/null
+++ b/3511/CH11/EX11.5/Ex11_5.sce
@@ -0,0 +1,41 @@
+clc;
+P01=4.6; // Total head inlet pressure in bar
+T01=700+273.15; // Total head inlet temperature in kelvin
+P2=1.6; // Static head pressure at mean radius in bar
+Dm_h=10; // Mean blade diameter/blade height
+lc=0.1; // Nozzle losses coefficient
+alpha_2=60; // Nozzle outlet angle in degree
+Cp=1.147; // Specific heat in kJ/kg K
+r=1.33; // Specific heat ratio
+m=20; // Mass flow rate in kg/s
+R=284.6; // characteristic gas constant in J/kg K
+
+T_2=T01*(P2/P01)^((r-1)/r); // Isentropic temperature after expansion
+T2=(lc*T01+T_2)/(1+lc); // Actual temperature after expansion
+c2=sqrt(2*Cp*10^3*(T01-T2)); // Absolute velocity
+// From velocity triangles
+ca=c2*cosd(alpha_2);
+row=P2*10^5/(R*T2); // Density of gas
+A=m/(ca*row); // Area
+Dm=sqrt (A*Dm_h/3.14); // Mean Diameter
+h=Dm/10; // Blade height
+rm=Dm/2; // Mean radius
+// At root
+r_root=(Dm-h)/2;
+//At the tip
+r_tip=(Dm+h)/2;
+// Free vorte flow
+ct_mean=c2*sind (alpha_2);
+// At the root
+ct2_root=(ct_mean*rm)/r_root;
+alpha2_root=atand(ct2_root/ca);
+c2_root=ct2_root/sind (alpha2_root);
+T2_root=T01-c2_root^2/(2*Cp*10^3);
+// At the tip
+ct2_tip=ct_mean*rm/r_tip;
+alpha2_tip = atand (ct2_tip/ca);
+c2_tip=ct2_tip/sind(alpha2_tip);
+T2_tip=T01-c2_tip^2/(2*Cp*10^3);
+
+disp ("degree",alpha2_root,"Discharge angle at the root = ","m/s",c2_root,"Gas velocity at the root = ","K",T2_root,"Gas Temperature at the root = ","A the Root");
+disp ("degree",alpha2_tip,"Discharge angle at the tip = ","m/s",c2_tip,"Gas velocity at the tip = ","K",T2_tip,"Gas Temperature at the tip = ","A the tip");