1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
|
clear;
clc;
disp('Example 11.2')
// aim : To determine the
// (a) the diameter of the cylinder
// (b) piston stroke
// (c) actual steam consumption and indicated thermal efficiency
// Given values
P = 900;// inlet pressure, [kN/m^2]
Pb = 140;// exhaust pressure, [kN/m^2]
cut_off =.4;// [stroke]
k = .8;// diagram factor
rs = 1.2;// stroke to bore ratio
N = 4;// engine speed, [rev/s]
ip = 22.5;// power output from the engine, [kW]
// solution
// (a)
r = 1/cut_off;// expansion ratio
Pm = P/r*(1+log(r))-Pb;// mean effective pressure, [kN/m^2]
Pm = Pm*k;// actual mean effective pressure, [kN/m^2]
// using ip=Pm*L*A*N
// and L=r*d; where L is stroke and d is bore
d = (ip/(Pm*rs*%pi/4*2*N))^(1/3);// diameter of the cylinder, [m]
mprintf('\n (a) The diameter of the cylinder is = %f mm\n',d*1000);
// (b)
L = rs*d;// stroke, [m]
mprintf('\n (b) The piston stroke is = %f mm\n',L*1000);
// (c)
SV = %pi/4*d^2*L;// stroke volume, [m^3]
V = SV*cut_off*2*240*60;// volume of steam consumed per hour, [m^3]
v = .2148;// specific volume at 900 kN/m^2, [m^3/kg]
SC = V/v;// steam consumed/h, [kg]
ASC = 1.5*SC;// actual steam consumption/h, [kg]
mprintf('\n (c) The actual steam consumption/h is = %f kg\n',ASC);
m_dot = ASC/3600;// steam consumption,[kg/s]
// from steam table
hg = 2772.1;// specific enthalpy of inlet steam, [kJ/kg]
hfe = 458.4;// specific liquid enthalpy at exhaust pressure, [kJ/kg]
ITE = ip/(m_dot*(hg-hfe));// indicated thermal efficiency
mprintf('\n The indicated thermal efficiency is = %f percent\n',ITE*100);
// End
|