1 files changed, 29 insertions, 0 deletions
diff --git a/2223/CH18/EX18.22/Ex18_22.sce b/2223/CH18/EX18.22/Ex18_22.sce
new file mode 100755
index 000000000..e42240371
--- /dev/null
+++ b/2223/CH18/EX18.22/Ex18_22.sce
@@ -0,0 +1,29 @@
+// scilab Code Exa 18.22 Isentropic or Stage Terminal Velocity for Turbines 
+
+T01=1273;  // in Kelvin
+funcprot(0);
+p01=5; //  Initial Pressure in bar
+p02=3.5; //  exit gas Pressure in bar
+cp=1.005; // Specific Heat at Constant Pressure in kJ/(kgK)
+gamma=1.4; // Specific Heat Ratio
+m=28; // mass flow rate of the gas in kg/s
+n_tt=0.84; // stage efficiency
+shi=1.7; // stage loading coefficient
+pr_0=p01/p02;
+delh01_03ss=cp*T01*(1-(pr_0^((1-gamma)/gamma)));
+
+//part(a)stage terminal velocity
+c0=sqrt(2*delh01_03ss*1e3);
+disp("m/s",c0,"(a)stage terminal velocity is")
+
+// part(b)isentropic blade to gas speed ratio
+sigma_s=sqrt(0.5*n_tt/shi);
+disp(sigma_s,"(b)the isentropic blade to gas speed ratio is")
+
+//part(c) peripheral speed of the rotor
+u=sigma_s*c0;
+disp("m/s",u,"(c)peripheral speed of the rotor is")
+
+//part(d) the power developed
+P=m*n_tt*delh01_03ss;
+disp("MW",P*1e-3,"(d) the power developed is")