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 --- 32/CH9/EX9.11/9_11.sce | 68 ++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 68 insertions(+) create mode 100755 32/CH9/EX9.11/9_11.sce (limited to '32/CH9/EX9.11/9_11.sce') diff --git a/32/CH9/EX9.11/9_11.sce b/32/CH9/EX9.11/9_11.sce new file mode 100755 index 000000000..09abfcb56 --- /dev/null +++ b/32/CH9/EX9.11/9_11.sce @@ -0,0 +1,68 @@ +//pathname=get_absolute_file_path('9.11.sce') +//filename=pathname+filesep()+'9.11-data.sci' +//exec(filename) +//Pressure at which air is supplied(in bar): +p1=1 +p8=p1 +//Temperature at which air is supplied(in K): +T1=27+273 +//Maximum temperature in the cycle(in K): +T5=1000 +//Pressure at state 6(in bar): +p6=3 +p7=p6 +p4=10 +p5=p4 +p3=3 +//Temperature at state 7(in K): +T7=995 +//Calorific value of fuel(in kJ/kg): +c=42000 +//Value of Cp(in kJ/kg): +Cp=1.0032 +//Air flow in compressor(in kg/s): +m=30 +//Isentropic efficiency of compression: +nc=0.85 +//Isebtropic efficiency of expansion: +ne=0.90 +//Adiabatic index of compression: +r=1.4 +//Pressure ratio for perfect intercooling: +rp=sqrt(10) +//Temperature at state 2'(in K): +T21=T1*rp^((r-1)/r) +//Temperature at state 2(in K): +T2=(T21-T1)/nc+T1 +//For perfect intercooling: +T3=T1 +//Temperature at state 4'(in K): +T41=T3*(p4/p3)^((r-1)/r) +//Temperature at state 4(in K): +T4=(T41-T3)/nc+T3 +//Total compressor work(in kJ/kg): +Wc=2*Cp*(T2-T1) +//Temperature at state 6'(in K): +T61=T5*(p6/p5)^((r-1)/r) +//Temperature at state 6(in K): +T6=T5-(T5-T61)*ne +//Temperature at state 8'(in K): +T81=T7*(p8/p7)^((r-1)/r) +//Temperature at state 8(in K): +T8=T7-(T7-T81)*ne +//Expansion work output per kg air(in kJ/kg): +Wt=Cp*(T5-T6+T7-T8) +//Heat added per kg air(in kJ/kg): +qa=Cp*(T5-T4+T7-T6) +//Fuel required per kg of air: +mf=qa/c +//Air-fuel ratio: +afr=1/mf +//Net output(in kW): +Wnet=(Wt-Wc)*m +//Thermal efficiency: +nth=(Wt-Wc)/qa +printf("\nRESULT") +printf("\nThermal efficiency = %f percent",nth*100) +printf("\nNet output = %f kW",Wnet) +printf("\nA/F ratio = %f",afr) \ No newline at end of file -- cgit