diff options
Diffstat (limited to '2006/CH5/EX5.16')
| -rwxr-xr-x | 2006/CH5/EX5.16/ex5_16.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/2006/CH5/EX5.16/ex5_16.sce b/2006/CH5/EX5.16/ex5_16.sce new file mode 100755 index 000000000..55d20839e --- /dev/null +++ b/2006/CH5/EX5.16/ex5_16.sce @@ -0,0 +1,23 @@ +clc;
+p1=3; // Pressre of air at state 1 in bar
+p2=p1; // constant pressure process
+T1=450; // Temperature of air at state 1 in kelvin
+T2=1250; // Temperature of air at state 2 in kelvin
+T3=1000; // Temperature of air at state 3 in kelvin
+V3=50; // Velocity of air at state 3 in m/s
+T4=800; // Temperature of air at state 4 in kelvin
+Cpo=1.0035; // Specific heat at constant pressure in kJ/kg K
+// (a).Combustion chamber
+q=Cpo*(T2-T1); // Heat added to air
+disp ("kJ/kg (round off error)",q,"Heat added to air = ","(a).Combustion chamber");
+// (b).Turbine
+k=1.4; // Index of adiabatic process
+w=Cpo*(T2-T3)-V3^2/2000; // Work done
+disp ("kJ/kg (round off error)",w,"Work done = ",("(b).Turbine)"));
+// (c).Nozzle
+V4=sqrt (2*Cpo*10^3*(T3-T4)+V3^2); // Velocity of air leaving the nozzle
+disp ("m/s (round off error)",V4,"Velocity of air leaving the nozzle = ","(c).Nozzle");
+// (d).Pressure drop
+p3=p2*(T3/T2)^(k/(k-1)); // Pressure of air leaving turbine
+p4=p3*(T4/T3)^(k/(k-1)); // Pressure of air leaving nozzle
+disp ("bar ",p4,"Pressure of air leaving nozzle = ","bar",p3,"Pressure of air leaving turbine = ","(d).Pressure drop");
|