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/CH7/EX7.4/Ex7_4.sce | 35 +++++++++++++++++++++++++++++++++++ 1 file changed, 35 insertions(+) create mode 100755 2223/CH7/EX7.4/Ex7_4.sce (limited to '2223/CH7/EX7.4/Ex7_4.sce') diff --git a/2223/CH7/EX7.4/Ex7_4.sce b/2223/CH7/EX7.4/Ex7_4.sce new file mode 100755 index 000000000..e654fb602 --- /dev/null +++ b/2223/CH7/EX7.4/Ex7_4.sce @@ -0,0 +1,35 @@ +// scilab Code Exa 7.4 Calculation on a single stage gas turbine + +gamma_g=1.33; +gamma=1.4 +R_g=284.1; +R=287; +P=1e3; // Power Output in kW +N1=3e3; // Speed in RPM +n_t=0.87; // efficiency +cp_g=1.145; // Specific Heat of gas at Constant Pressure in kJ/(kgK) +cp_a=1.0045; // Specific Heat of air at Constant Pressure in kJ/(kgK) + +// part (a)mass flow rate of the gas through the turbine +T01=1000; // Entry Temperature in Kelvin +p01=2.5; // Entry Pressure in bar +T01a=500; // Entry Temperature of air in Kelvin +p01a=2; // Entry Pressure of air in bar +p02=1; // Exit Pressure in bar +pr0=p01/p02; +T02=T01*(pr0^(-((gamma_g-1)/gamma_g))); +delh_s1=cp_g*(T01-T02)*n_t; +m_g=P/delh_s1; +disp("kg/s",m_g,"(a)mass flow rate of the gas through the turbine is") + +// part (b)speed, mass flow rate, pressure ratio and power required +N2=sqrt(1/2)*5*N1; +disp("rpm",N2,"(b)(i)speed of a geometrically similar compressor is") +a=0.2; // a=D2/D1; +m2=(a^2)*sqrt(R_g/R)*sqrt(T01/T01a)*(p01a/p01)*m_g; +disp("kg/s",m2,"(b)(ii)mass flow rate of a geometrically similar turbine is") +delh_s2=0.5*delh_s1; +P2=m2*delh_s2; +disp("kW",P2,"(b)(iii)Power developed is") +pr=(1-(delh_s2/(cp_a*T01a*n_t)))^(-1/((gamma-1)/gamma)); +disp(pr,"(b)(iv)pressure ratio of a geometrically similar turbine is") -- cgit