Курсовая работа: Математическое моделирование тепловой работы вращающейся печи

GB:=GP+0.322*sqrt (ROD/R0);

{Средняя плотность факела на расчетном участке, кг/м3}

Ro:=R0*sqr (GB/X/0.322);

RS:=sqrt (Ro/R0); {Соотношение плотностей}

{Расстояние до полюса струи, D0*м}

X0:=2.37/RS;

{Расчетная длина факела, D0*м}

LF:=X0*((NB+1)*(2*SC+1) – 1);

{Расстояние от полюса струи до расчетного сечения и конца факела, D0*м}

XX:=X+X0; LX:=LF+X0;

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

RF:=0.211*XX*D0;

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

SC2:=2*SC; SC1:=SC2+1;

GT:=XX/X0/NB*(LX/XX/SC1+SC2/SC1*exp (ln(XX/LX)/SC2) – 1);

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

HHF:=(H0+HB*GB+QT*(1-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;

GasPost; Gasheat; HeatLine;

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

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

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

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

Result;

end;