183498 (629852), страница 3

Просмтор этого файла доступен только зарегистрированным пользователям. Но у нас супер быстрая регистрация: достаточно только электронной почты!

Текст из файла (страница 3)

{Плотность газов в расчётном сечении факела, кг/м3}

ROF:=ROV*273/TK;

GasPost; Gasheat; HeatLine;

until ABS (1-ZT/TF)<0.0001; {конец итераций}

{Радиус фронта пламени, м}

FR:=XX/7.9*sqrt (exp(ln (LX/XX)/SC2) – 1);

{Вывод результатов на печать}

Result;

end;

{–}

{Параметры факела в начальном сечении зоны догорания}

GTB:=GT; GB:=GB0/G0; HHB:=HHF; RF:=DW/2;

XB:=X; XXB:=XB+X0; LXB:=LX;

MG:=(1-exp (ln(XXB/LXB)/SC2))/NB/X0/GTB;

While GT>0.01 do begin {переход к очередному сечению}

{Расстояние от горелки до расчётного сечения, D0*м}

X:=X+1;

{Параметры факела в предыдущем сечении}

TP:=TF; QP:=QF;

repeat {начало итераций}

{Расход несгоревшего топлива через расчетное сечение, кг/с}

GT:=GTB*exp (MG*(XB-X));

{Средняя энтальпия газов в расчетном сечении факела, кДж/кг}

HHF:=HHB+QT*(GTB-GT)/(1+GB); {без теплообмена}

HF:=HHF-QF/(1+GB)/G0; {c теплообменом}

{Средняя температура в расчетном сечении факела, °C}

ZT:=TF; TF:=(HF-fHDis)/CF; TK:=TF+273;

{Плотность газов в расчетном сечении факела, кг/м3}

ROF:=ROV*273/TK;

Gasheat; HeatLine;

until ABS (1-ZT/TF)<0.0001; {конец итераций}

{Вывод результатов на печать}

Result;

end;

Nomina;

End;

{============================================================================}

Procedure tFurnace. GasPost;

{============================================================================}

{Расчет параметров диффузионного факела}

var RO1, RO2, RO3, ROG, VT, VB, VG, VF: real;

Begin

P1:=V1/VG0; P2:=V2/VG0; P3:=V3/VG0;

RO1:=44/22.41;

RO2:=18/22.41;

RO3:=28/22.41;

{Плотность стехиометрических продуктов горения, кг/м3 (н.у.)}

ROG:=RO1*P1+RO2*P2+RO3*P3;

{Объемный расход несгоревшего топлива, м3/с (н.у.)}

VT:=GT*G0/RO;

{Объемный расход воздуха, м3/с (н.у.)}

VB:=(GB – (1-GT)*NB)*G0/ROB;

{Объемный расход продуктов горения, м3/с (н.у.)}

VG:=(1-GT)*(1+NB)*G0/ROG;

{Объемный расход газообразной среды, м3/с (н.у.)}

VF:=VT+VB+VG;

{Объемные доли компонентов в газообразной среде}

PT:=VT/VF; {Топливо}

PB:=VB/VF; {Воздух}

PG:=VG/VF; {Продукты горение}

{Плотность газообразной среды, м3/кг (н.у.)}

ROV:=RO*PT+ROB*PB+ROG*PG;

End;

{============================================================================}

Procedure tFurnace. GasHeat;

{============================================================================}

{Расчет теплоемкости газообразной среды}

var CB, CT, C1, C2, C3: real;

Begin

{Объёмная теплоёмкость природного газа, кДж/м3*К (н.у.)}

IF TF=0 THEN CT:=2.522+0.0005815*TF ELSE

CT:=(-317.9+2.522*TF+0.0005815*sqr(TF)+86840/(TF+273))/TF;

{Объемная теплоемкость воздуха, кДж/м3*К (н.у.)}

CB:=1.287+0.0001201*TF;

{Объёмная теплоёмкость углекислого газа, кДж/(м3*К) (н.у.)}

IF TF=0 THEN C1:=2.081+0.0002017*TF ELSE

C1:=(-139.5+2.081*TF+0.0002017*sqr(TF)+38110/(TF+273))/TF;

{Объёмная теплоёмкость водяного пара, кДж/(м3*К) (н.у.)}

C2:=1.49+0.0002303*TF;

{Объёмная теплоёмкость азота, кДж/(м3*К) (н.у.)}

C3:=1.28+0.0001103*TF;

{Объёмная теплоёмкость продуктов горения, кДж/м3*К (н.у.)}

CG:=C1*P1+C2*P2+C3*P3;

{Объёмная теплоёмкость газов в факеле, кДж/(м3*К) (н.у.)}

CV:=CT*PT+CB*PB+CG*PG;

{Массовая теплоёмкость газов в факеле, кДж/кг*К}

CF:=CV/ROV;

End;

{============================================================================}

Function tFurnace.fHdis:real;

{============================================================================}

{Расчет теплоты диссоциации факела}

var PCO2, PH2O, KCO2, KH2O, ACO2, AH2O, PS, LR, K1, K2: real;

Begin

fHdis:=0;

if TF>1500 then begin

{Объёмные доли трехатомных газов}

PCO2:=P1*PG; {…углекислого газа}

PH2O:=P2*PG; {…водяного пара}

{Константа равновесия для углекислого газа}

KCO2:=exp (ln(10)*(4.47–14700/TF));

{Константа равновесия для водяного пара}

KH2O:=exp (ln(10)*(3.05–13160/TF));

{Степень диссоциации углекислого газа}

ACO2:=exp (ln(2*sqr(KCO2)/PCO2)/3);

{Степень диссоциации водяного пара}

AH2O:=exp (ln(2*sqr(KH2O)/PH2O)/3);

{Теплота диссоциации трехатомных газов, кДж/кг}

fHdis:=(12630*ACO2*PCO2+10800*AH2O*PH2O)/ROV;

end;

End;

{============================================================================}

Procedure tFurnace. HeatLine;

{============================================================================}

{Расчет результирующего излучения факела}

var TG, FF, FW, PC, PH, PS, LR, K1, K2: real;

Begin

{Средняя температура факела на малом участке,°C}

TG:=(TP+TF) /2;

{Расчетная поверхность малого участка факела, м2}

FF:=2*PI*RF*D0;

{Внутренняя поверхность стенки печи, м2}

FW:=PI*DW*D0;

{Объёмные доли излучающих газов}

PC:=P1*PG+PT; {углекислый газ + топливо}

PH:=P2*PG; {водяной пар}

PS:=PC+PH; {трёхатомные газы}

{Эффективная толщина излучающего слоя, м}

LR:=1.8*RF;

{Коэффициент поглощения при температуре газов, 1/м}

K1:=PS*((0.78+1.6*PH)/sqrt (PS*LR) – 0.1)*(1–0.37*(TG+273)/1000);

{Степень черноты факела}

EG:=1-exp (-K1*LR);

{Коэффициент поглощения при температуре стенки, 1/м}

K2:=PS*((0.78+1.6*PH)/sqrt (PS*LR) – 0.1)*(1–0.37*(TM+273)/1000);

{Поглощательная способность факела}

AG:=1-exp (-K2*LR);

{Приведенная степень черноты}

ES:=1/(1/AG+(1/EW1–1)*FF/FW);

{Поток результирующего излучения факела, кВт}

QF:=QP+ES*SI*(EG/AG*sqr (sqr(TG+273)) – sqr (sqr(TM+273)))*FF/1000;

End;

{============================================================================}

Procedure tFurnace. Names;

{============================================================================}

Begin

Assign (File1,'Furnace.txt'); Rewrite(File1);

Writeln(File1); Writeln (file1,' ИСХОДНЫЕ ДАННЫЕ:');

Writeln (file1,' Внутренний диаметр печи DW:=', DW:5:2,' м');

Writeln (file1,' Toлщина футеровки HW:=', HW*1000:4:0,' мм');

Writeln (file1,' Температура обжигаемого материала TM:=', TM:5:0,'°C');

Writeln (file1,' Низшая теплота сгорания топлива QT:=', QT:5:0,' кДж/кг');

Writeln (file1,' Плотность топлива R0:=', R0:3:3,' кг/м3 (н.у)');

Writeln (file1,' Расход топлива через горелку G0:=', G0:5:3,' кг/с');

Writeln (file1,' Скорость истечения топлива из горелки U0:=', U0:5:1,' м/с');

Writeln (file1,' Диаметр горелки D0:=', D0*1000:3:0,' мм');

Writeln (file1,' Коэффициент избытка воздуха ALB:=', ALB:3:2);

Writeln (file1,' Температура горячего воздуха TB:=', TB:4:0,'°C');

Writeln(File1);

Writeln (file1,' X LF GT GB EG AG ES HF TF QF ');

End;

{============================================================================}

Procedure tFurnace. Result;

{============================================================================}

Begin

Append(File1);

if X mod 2 = 0 then

Writeln (file1, x:3,' ', lf:5:1,' ', gt:6:4,' ', gb:6:3,' '

eg:5:3,' ', ag:5:3,' ', es:5:3,' ', hf:6:0,' ', tf:6:0,' ', qf:5:0);

End;

{============================================================================}

Procedure tFurnace. Nomina;

{============================================================================}

Begin

Append(File1); Writeln(File1); Writeln (File1,' ОБОЗНАЧЕНИЯ:');

Writeln (File1,' X – относительное расстояние x/d0 (номер расчетного сечения)');

Writeln (File1,' LF – относительная длина факела LФ/d0');

Writeln (File1,' GT – доля несгоревшего топлива от поступившего в горелку Gт/G0');

Writeln (File1,' GB – относительное количество воздуха, вовлеченного в факел, Gb/G0');

Writeln (File1,' EG – степень черноты газообразной среды в факеле');

Writeln (File1,' AG – поглощательная способность факела');

Writeln (File1,' Es – приведенная степень черноты факела');

Writeln (File1,' HF – энтальпия газообразной среды в факеле');

Writeln (File1,' Tf – средняя температура в поперечном сечении факела,°C');

Writeln (File1,' Qf – теплоотдача факела, кВт');

Writeln(File1);

Close(file1);

End;

{============================================================================}

BEGIN

Furnace. Fakel;

END.

Результаты выполнения программы №1:

X LF GT GB EG AG ES HF TF QF

2 148.6 0.9708 0.475 0.254 0.123 0.122 1195 306 -11

4 171.4 0.9472 0.824 0.286 0.146 0.146 1737 390 -31

6 184.8 0.9259 1.147 0.308 0.164 0.163 2073 454 -62

8 193.9 0.9059 1.457 0.325 0.178 0.176 2308 512 -102

10 200.5 0.8866 1.761 0.337 0.190 0.188 2483 565 -151

12 205.6 0.8679 2.062 0.347 0.201 0.199 2619 615 -211

14 209.6 0.8498 2.359 0.354 0.211 0.208 2727 663 -279

16 212.9 0.8320 2.655 0.360 0.220 0.216 2814 710 -356

18 215.6 0.8146 2.949 0.364 0.228 0.223 2886 755 -442

20 217.8 0.7976 3.243 0.367 0.235 0.230 2946 798 -534

22 219.8 0.7809 3.535 0.369 0.242 0.236 2996 840 -634

24 221.5 0.7645 3.827 0.370 0.249 0.242 3038 881 -739

26 223.0 0.7484 4.118 0.370 0.254 0.247 3073 921 -849

28 224.3 0.7325 4.409 0.370 0.260 0.251 3103 960 -962

30 225.5 0.7169 4.700 0.369 0.265 0.255 3128 997 -1078

32 226.5 0.7016 4.990 0.367 0.270 0.259 3149 1034–1195

34 227.5 0.6866 5.279 0.366 0.274 0.263 3166 1069–1312

36 228.3 0.6717 5.569 0.363 0.278 0.266 3180 1103–1428

38 229.1 0.6571 5.858 0.363 0.282 0.269 3191 1136–1544

40 229.8 0.6428 6.147 0.358 0.286 0.272 3199 1168–1654

42 230.5 0.6287 6.435 0.355 0.289 0.274 3205 1199–1759

44 231.1 0.6147 6.724 0.352 0.293 0.277 3209 1229–1859

46 231.7 0.6010 7.012 0.348 0.296 0.279 3210 1258–1953

48 232.2 0.5875 7.300 0.345 0.299 0.281 3210 1286–2038

50 232.7 0.5743 7.588 0.341 0.301 0.282 3208 1313–2116

52 233.2 0.5612 7.876 0.337 0.304 0.284 3204 1339–2184

54 233.6 0.5483 8.163 0.334 0.306 0.285 3199 1364–2242

56 234.0 0.5356 8.451 0.330 0.308 0.286 3192 1388–2289

58 234.4 0.5231 8.739 0.326 0.310 0.288 3184 1411–2325

60 234.8 0.5108 9.026 0.322 0.312 0.289 3174 1433–2350

62 235.1 0.4986 9.314 0.319 0.314 0.289 3164 1453–2363

64 235.4 0.4867 9.601 0.315 0.316 0.290 3152 1474–2364

66 235.7 0.4749 9.889 0.311 0.317 0.291 3140 1493–2353

68 236.0 0.4633 10.177 0.308 0.319 0.291 3127 1510–2329

70 236.3 0.4518 10.465 0.304 0.320 0.292 3113 1528–2293

72 236.5 0.4406 10.753 0.301 0.321 0.292 3098 1544–2244

74 236.7 0.4295 11.041 0.297 0.323 0.292 3082 1560–2184

76 236.9 0.4185 11.330 0.294 0.324 0.293 3066 1574–2111

78 237.1 0.4078 11.618 0.291 0.325 0.293 3049 1588–2026

80 237.3 0.3971 11.907 0.288 0.326 0.293 3032 1602–1929

82 237.5 0.3867 12.196 0.285 0.327 0.293 3014 1614–1821

84 237.6 0.3764 12.485 0.282 0.327 0.293 2996 1626–1701

86 237.8 0.3662 12.775 0.280 0.328 0.293 2977 1636–1571

88 237.9 0.3562 13.065 0.277 0.329 0.293 2958 1647–1430

90 238.0 0.3464 13.355 0.275 0.329 0.292 2939 1656–1279

92 238.1 0.3367 13.645 0.273 0.330 0.292 2920 1665–1118

94 238.2 0.3271 13.936 0.270 0.330 0.292 2900 1673 -947

96 238.3 0.3177 14.227 0.268 0.331 0.291 2880 1681 -767

98 238.4 0.3085 14.519 0.267 0.331 0.291 2861 1688 -579

100 238.5 0.2994 14.810 0.265 0.331 0.291 2840 1694 -382

102 238.5 0.2904 15.103 0.263 0.332 0.290 2820 1700 -177

104 238.6 0.2816 15.395 0.262 0.332 0.290 2800 1705 36

106 238.6 0.2729 15.689 0.261 0.332 0.289 2780 1710 256

108 238.7 0.2643 15.982 0.259 0.332 0.288 2759 1714 483

110 238.7 0.2559 16.276 0.258 0.332 0.288 2739 1717 716

112 238.7 0.2477 16.571 0.257 0.332 0.287 2719 1721 956

114 238.7 0.2395 16.866 0.256 0.333 0.287 2698 1723 1201

116 238.7 0.2316 17.161 0.256 0.333 0.286 2678 1726 1452

118 238.7 0.2237 17.457 0.255 0.332 0.285 2658 1727 1708

120 238.7 0.2160 17.753 0.255 0.332 0.285 2638 1729 1969

122 238.7 0.2084 18.050 0.254 0.332 0.284 2618 1730 2234

124 238.7 0.2010 18.348 0.254 0.332 0.283 2598 1730 2503

126 238.7 0.1937 18.493 0.235 0.309 0.273 2593 1733 2745

128 238.7 0.1867 18.493 0.233 0.309 0.273 2602 1738 2959

130 238.7 0.1800 18.493 0.232 0.309 0.273 2611 1743 3176

132 238.7 0.1735 18.493 0.230 0.309 0.273 2619 1748 3395

134 238.7 0.1672 18.493 0.229 0.309 0.273 2627 1752 3617

136 238.7 0.1612 18.493 0.228 0.309 0.273 2634 1756 3842

138 238.7 0.1554 18.493 0.227 0.309 0.273 2640 1759 4068

140 238.7 0.1498 18.493 0.226 0.309 0.273 2645 1762 4296

142 238.7 0.1444 18.493 0.225 0.309 0.273 2651 1765 4526

144 238.7 0.1392 18.493 0.224 0.309 0.273 2655 1767 4758

146 238.7 0.1341 18.493 0.224 0.309 0.273 2659 1770 4990

148 238.7 0.1293 18.493 0.223 0.309 0.273 2663 1772 5224

150 238.7 0.1246 18.493 0.223 0.309 0.273 2666 1773 5458

152 238.7 0.1201 18.493 0.222 0.309 0.273 2669 1775 5694

154 238.7 0.1158 18.493 0.222 0.309 0.273 2671 1776 5930

156 238.7 0.1116 18.493 0.221 0.309 0.273 2673 1777 6166

158 238.7 0.1076 18.493 0.221 0.309 0.273 2674 1778 6404

160 238.7 0.1037 18.493 0.221 0.309 0.273 2675 1778 6641

162 238.7 0.1000 18.493 0.221 0.309 0.273 2676 1779 6879

164 238.7 0.0964 18.493 0.221 0.309 0.273 2676 1779 7117

166 238.7 0.0929 18.493 0.221 0.309 0.273 2676 1779 7355

168 238.7 0.0895 18.493 0.221 0.309 0.273 2676 1779 7593

170 238.7 0.0863 18.493 0.221 0.309 0.273 2676 1779 7831

172 238.7 0.0832 18.493 0.221 0.309 0.273 2675 1778 8068

174 238.7 0.0802 18.493 0.221 0.309 0.273 2674 1778 8306

176 238.7 0.0773 18.493 0.221 0.309 0.273 2673 1777 8543

178 238.7 0.0745 18.493 0.222 0.309 0.273 2671 1776 8779

180 238.7 0.0718 18.493 0.222 0.309 0.273 2669 1775 9016

Характеристики

Список файлов курсовой работы