From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 2223/CH18/EX18.6/Ex18_6.sce | 30 ++++++++++++++++++++++++++++++ 1 file changed, 30 insertions(+) create mode 100755 2223/CH18/EX18.6/Ex18_6.sce (limited to '2223/CH18/EX18.6/Ex18_6.sce') diff --git a/2223/CH18/EX18.6/Ex18_6.sce b/2223/CH18/EX18.6/Ex18_6.sce new file mode 100755 index 000000000..fae4d033b --- /dev/null +++ b/2223/CH18/EX18.6/Ex18_6.sce @@ -0,0 +1,30 @@ +// scilab Code Exa 18.6 Calculations on a Gas Turbine + +m=472; // flow rate of hot gases in kg/s +T01=1335; // Turbine inlet temp in Kelvin +p01=10; // Turbine Inlet Pressure in bar +c2=150; // exit velocity in m/s +pr0=10; // Turbine pressure ratio +gamma_g=1.67; +T2=560; // Temperature of gases at exit in Kelvin +cp_g=1.157; // Specific Heat of gas at Constant Pressure in kJ/(kgK) + +// part(a) Determining total to total efficiency +T02=T2+(0.5*(c2^2)/(cp_g*1e3)); +T02s=T01/(pr0^((gamma_g-1)/gamma_g)); +n_tt=(T01-T02)/(T01-T02s); +disp("%",n_tt*100,"(a)total to total efficiency is") + + +// part(b) Determining total to static efficiency +T2s=T02s-(0.5*(c2^2)/(cp_g*1e3)); +n_ts=(T01-T02)/(T01-T2s); +disp("%",n_ts*100,"(b)total to static efficiency is") + +// part(c) Determining the polytropic efficiency +n_p=((gamma_g)/(gamma_g-1))*((log(T01/T02))/(log(pr0))); +disp("%",n_p*100,"(c)polytropic efficiency is") + +// part(d) Determining power developed by the turbine +P=m*cp_g*(T01-T02); +disp("MW",P/1e3,"(d)Power developed by the turbine is") -- cgit