diff options
Diffstat (limited to '3511/CH11')
| -rw-r--r-- | 3511/CH11/EX11.1/Ex11_1.sce | 19 | ||||
| -rw-r--r-- | 3511/CH11/EX11.2/Ex11_2.sce | 33 | ||||
| -rw-r--r-- | 3511/CH11/EX11.3/Ex11_3.sce | 22 | ||||
| -rw-r--r-- | 3511/CH11/EX11.4/Ex11_4.sce | 26 | ||||
| -rw-r--r-- | 3511/CH11/EX11.5/Ex11_5.sce | 41 |
5 files changed, 141 insertions, 0 deletions
diff --git a/3511/CH11/EX11.1/Ex11_1.sce b/3511/CH11/EX11.1/Ex11_1.sce new file mode 100644 index 000000000..ad4d88b9c --- /dev/null +++ b/3511/CH11/EX11.1/Ex11_1.sce @@ -0,0 +1,19 @@ +clc;
+p02=6; // Inlet pressure in bar
+T02=900; // Inlet temperature in kelvin
+p0fs=1; // Outlet pressure in bar
+eff_isenT=0.85; // insentropic efficiency of turbine
+alpha_2=75; // Nozzle outlet angle in degree
+u=250; // Mean blade velocity in m/s
+Cp=1.15*10^3; // Specific heat in J/ kg K
+r=1.333; // Specific heat ratio
+
+T0fs=T02/(p02/p0fs)^((r-1)/r); // Isentropic temperature at the exit of the final stage
+Del_Toverall=eff_isenT*(T02-T0fs); // Actual overall temperature drop
+c2=2*u/sind (alpha_2); // absolute velocity
+c3= c2*cosd (alpha_2);// absolute velocity
+c1=c3; // From velocity triangles
+Del_Tstage=(c2^2-c1^2)/(2*Cp); // Stage temperature drop
+n=Del_Toverall/Del_Tstage; // Number of stages
+
+disp (round (n),"Number of stages n =");
diff --git a/3511/CH11/EX11.2/Ex11_2.sce b/3511/CH11/EX11.2/Ex11_2.sce new file mode 100644 index 000000000..f2db7249d --- /dev/null +++ b/3511/CH11/EX11.2/Ex11_2.sce @@ -0,0 +1,33 @@ +clc;
+N=10000; // Speed of gas turbine in rpm
+T01=700+273.15; // Total head temperature at nozzle entry in kelvin
+P01=4.5; //Total head pressure at nozzle entry in bar
+P02=2.6; // Outlet pressure from nozzle in bar
+p3=1.5;// Pressure at trbine outlet annulus in bar
+M=0.5; // Mach number at outlet
+alpha_2=70; // outlet nozzle angle in degrees
+D=64; // Blade mean diameter in cm
+m=22.5; // Mass flow rate in kg/s
+eff_T=0.99; // turbine mechanical efficiency
+Cp=1.147; // Specific heat in kJ/kg K
+r=1.33; // Specific heat ratio
+fl=0.03; // frictional loss
+R=284.6; // characteristic gas constant in J/kg K
+
+eff_N=1-fl; // Nozzle efficiency
+T_02=(P02/P01)^((r-1)/r)*T01; // Isentropic temperature after expansion
+T02=T01-eff_N*(T01-T_02); // Actual temperature after expansion
+c2=sqrt (2*Cp*10^3*(T01-T02)); // Absolute velocity
+u=(3.14*D*10^-2*N)/60; // Mean blade velocity
+// From velocity triangles
+wt2=c2*sind (alpha_2)-u;
+ca=c2*cosd (alpha_2);
+beta_2=atand((wt2)/ca);
+T3=T02/(P02/p3)^((r-1)/r); // Assuming rotor losses are negligible
+c3=M*sqrt (r*R*T3); // Absolute velocity
+beta_3=atand(u/c3);
+ct2=c2*sind(alpha_2);
+P=eff_T*m*(ct2)*u/1000; // Power developed
+
+disp ("degree",beta_3,"Gas angle at exit = ","degree",beta_2,"Gas angle at entry","(i).");
+disp ("kW (roundoff error)",P,"Power developed = ","(ii).");
diff --git a/3511/CH11/EX11.3/Ex11_3.sce b/3511/CH11/EX11.3/Ex11_3.sce new file mode 100644 index 000000000..a35b84b23 --- /dev/null +++ b/3511/CH11/EX11.3/Ex11_3.sce @@ -0,0 +1,22 @@ +clc;
+alpha_2=65; // Nozzle discharge angle in degree
+c3=300; // Absolute velocity in m/s
+alpha_3=30; // in degrees
+
+ca2=c3*cosd (alpha_3); // Axial velocity
+c2=ca2/cosd(alpha_2); // Absolute velocity
+// ca3=ca2=ca and equal blade angles then
+ca=ca2;
+beta_2=atand((c2*sind(alpha_2)+c3*sind(alpha_3))/(2*ca)); // Blade angle
+beta_3=beta_2; // equal blade angles
+u=c2*sind(alpha_2)-ca2*tand(beta_2); // Mean blade velocity
+// From velocity triangles
+ct2=c2*sind(alpha_2);
+ct3=c3*sind(alpha_3);
+WT=u*(ct2+ct3)/1000; // Work done
+sigma=u/c2; // optimum speed ratio
+eff_B=4*(sigma*sind(alpha_2)-sigma^2);
+
+disp ("degree",beta_2,"Blade angle = beta_2= beta_3 = ");
+disp ("kJ/kg (roundoff error)",WT,"Power Produced = ");
+disp ("%",eff_B*100,"Blade efficiency = ");
diff --git a/3511/CH11/EX11.4/Ex11_4.sce b/3511/CH11/EX11.4/Ex11_4.sce new file mode 100644 index 000000000..8aef2d3bd --- /dev/null +++ b/3511/CH11/EX11.4/Ex11_4.sce @@ -0,0 +1,26 @@ +clc;
+P01=7; // Pressure at inlet in bar
+T01=300+273.15; // Temperature at inlet in kelvin
+P02=3; // Pressure at outlet in bar
+alpha_2=70; // Nozzle angle in degree
+eff_N=0.9; // Isentropic efficiency of nozzle
+WT=75; // Power Produced in kW
+Cp=1.15; // Specific heat in kJ/kg K
+r=1.33; // Specific heat ratio
+
+T_02=T01*(P02/P01)^((r-1)/r); // Isentropic temperature after expansion
+T02=T01-eff_N*(T01-T_02); // Actual temperature after expansion
+c2=sqrt (2*Cp*10^3*(T01-T02)); // Absolute velocity
+// For optimum blade speed ratio
+u=(c2*sind (alpha_2)/2); // Mean blade velocity
+beta_2=atand((c2*sind(alpha_2)-u)/(c2*cosd(alpha_2))); // Blade angle
+// From velocity triangles
+ct2=c2*sind(alpha_2);
+w2=c2*cosd(alpha_2)/cosd(beta_2);
+w3=w2; // Equal inlet and outlet angles
+beta_3=54; // in degrees
+ct3=w3*sind(beta_3)-u;
+m=(WT*10^3)/(u*(ct2+ct3)); // Gas mass flow rate
+
+disp ("degree",beta_2,"Blade angle = ");
+disp ("kg/s",m,"Gas Mass Flow Rate = ");
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");
|