blob: 6adefee16add08832cb4f3aee4cc4d3cf4d7d512 (
plain)
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
52
53
54
55
56
|
clc;
KE=1/2*(800*1000/3600)^2/1000;
T0=-50+273;
cp=1.005;
T0_=T0+[24.7/cp];
Ieff=0.9;
T0_s=Ieff*(T0_-T0)+T0;
y=1.4;
pa=0.24;
p0_=[(T0_s/T0)^[y/(y-1)]]*pa;
p0_2!p0_=10;
T0_2s=T0_*[p0_2!p0_^([y-1]/y)];
T0_2=T0_+(T0_2s-T0_)/Ieff;
p0_2=10*p0_;
p0_3=p0_2-(0.14);
T0_3=820+273;
meff=0.98;
cp2=1.15;
T0_4=T0_3-[cp*(T0_2-T0_)/(cp2*meff)];
T0_4s=T0_3-[cp*(T0_2-T0_)/(cp2*meff)]/0.92;
y2=1.333;
p0_4=3.24/[(T0_3/T0_4s)^{(y2/(y2-1))}]
pc=p0_4*([2/(y2+1)]^{y2/(y2-1)});
T0_5=[2/(y2+1)]*T0_4;
T0_5s=T0_4-{(T0_4-T0_5)/0.92};
p5=p0_4/[(T0_4/T0_5s)^(y2/{y2-1})];
R=cp2*(y2-1)/y2;
v5=R*T0_5*1000/10^5/p5;
T5=741.3//K
Cj=(y2*R*1000*T5)^0.5;
A=0.08;
m=A*Cj/v5;
Cg=222.2;
mt=m*(Cj-Cg)
pt=(p5-pa)*A*10^5;
Tt=pt+mt;
Q=m*cp2*(T0_3-T0_2)
C=43300;
mf=Q/meff/C;
SFC=mf*10^3/6453
disp("kg/kNs",SFC,"specific fuel consumption is")
|
