diff options
| author | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
|---|---|---|
| committer | prashantsinalkar | 2017-10-10 12:27:19 +0530 |
| commit | 7f60ea012dd2524dae921a2a35adbf7ef21f2bb6 (patch) | |
| tree | dbb9e3ddb5fc829e7c5c7e6be99b2c4ba356132c /3733/CH25/EX25.4 | |
| parent | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (diff) | |
| download | Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.gz Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.tar.bz2 Scilab-TBC-Uploads-7f60ea012dd2524dae921a2a35adbf7ef21f2bb6.zip | |
initial commit / add all books
Diffstat (limited to '3733/CH25/EX25.4')
| -rw-r--r-- | 3733/CH25/EX25.4/Ex25_4.sce | 62 |
1 files changed, 62 insertions, 0 deletions
diff --git a/3733/CH25/EX25.4/Ex25_4.sce b/3733/CH25/EX25.4/Ex25_4.sce new file mode 100644 index 000000000..1f3f16aac --- /dev/null +++ b/3733/CH25/EX25.4/Ex25_4.sce @@ -0,0 +1,62 @@ +// Example 25_4
+clc;funcprot(0);
+//Given data
+W_t=100*10^3;// kW
+W_g=(60/100)*W_t;// kW
+T_1=300;// K
+p_1=1;// bar
+T_3=1000+273;// K
+p_r=8;// Pressure ratio
+n_c=0.85;// Isentropic efficiency of compressor
+n_t=0.90;// Isentropic efficiency of both turbines
+n_com=0.95;// Combustion efficiency
+Gc=2500;// Rs./ton
+T_7=600+273;// K
+T_6=200+273;// K
+p_9=0.05;// bar
+C_pa=1;// kJ/kg.K
+C_pg=1.1;// kJ/kg.K
+r_a=1.4;// Specific heat ratio for air
+r_g=1.33;// Specific heat ratio for gases
+CV=40*10^3;// kJ/kg
+dT=10;// °C
+C_pw=4.2;// kJ/kg°C
+
+// Calculation
+// Considering compressor
+p_2=p_1*p_r;// bar
+T_2a=T_1*(p_r)^((r_a-1)/r_a);// K
+T_2=((T_2a-T_1)/n_c)+T_1;// K
+// Considering turbine
+T_4a=T_3/(p_r)^((r_g-1)/r_g);// K
+T_4=T_3-((T_3-T_4a)*n_t);// K
+T_5=T_3;// K
+// Considering heat balance in CC-I
+function[X]=mass(y);
+ X(1)=(y(1)*CV*n_com)-((y(2)+y(1))*C_pg*(T_3-T_2));
+ X(2)=(W_g)-(((y(2)+y(1))*C_pg*(T_3-T_4))-(y(2)*C_pa*(T_2-T_1)));
+ X(3)=(y(3)*CV*n_com)-((y(2)+y(1)+y(3))*C_pg*(T_5-T_4));
+endfunction
+y=[1 100 1];
+z=fsolve(y,mass);
+m_a1=z(2);// kg/sec
+m_f1=z(1);// kg/sec
+m_f2=z(3);// kg/sec
+AF_1=m_a1/m_f1;
+m_f=m_f1+m_f2;// kg/sec
+Q_s=(m_f*CV);// kW
+n=(W_t/Q_s)*100;// Efficiency of the plant in %
+// From h-s chart:
+h_7=3610;// kJ/kg
+// From steam table
+h_9=32.6;// kJ/kg
+m_s=((m_a1+m_f1+m_f2)*C_pg*(T_5-T_6))/(h_7-h_9);// kg/sec
+Afsf=m_a1/m_s;// Air flow to steam flow ratio
+Cf=((m_f*3600)/1000)*Gc;// Cost of fuel per hour in rupees
+E_g=W_t;// Energy generated per hour kWh
+Cg=Cf/E_g;// Cost of generation in rupees/kWh
+// From h-s chart:
+h_8=2220;// kJ/kg
+m_w=(m_s*3600*(h_8-h_9))/(C_pw*dT*1000);// Quantity of cooling water required in tons/hr
+printf('\n(i)Overall efficiency of the plant=%0.1f percentage \n(ii)Air flow to steam flow ratio=%0.2f \n(iii)Cost of generation=%0.2f rupees/kWh \n(iv)Quantity of cooling water required=%0.0f tons/hr',n,Afsf,Cg,m_w);
+// The answer provided in the textbook is wrong
|