Моделирование радиационного теплообмена с учетом неравновесности и селективности излучения в топках энергетических котлов тепловых электрических станций. П. Радиационные характеристики

Москаленко Николай Иванович; Загидуллин Рустем Абдулхакович; Паржин Сергей Николаевич; Зарипов Алмаз Вилевич; Локтев Николай Федорович

УДК 629.125:551.521

МОДЕЛИРОВАНИЕ РАДИАЦИОННОГО ТЕПЛООБМЕНА С УЧЕТОМ НЕРАВНОВЕСНОСТИ И СЕЛЕКТИВНОСТИ ИЗЛУЧЕНИЯ В ТОПКАХ ЭНЕРГЕТИЧЕСКИХ КОТЛОВ ТЕПЛОВЫХ ЭЛЕКТРИЧЕСКИХ СТАНЦИЙ II. РАДИАЦИОННЫЕ ХАРАКТЕРИСТИКИ

Н.И. МОСКАЛЕНКО, Р.А. ЗАГИДУЛЛИН, С.Н. ПАРЖИН, А.В. ЗАРИПОВ,

Н.Ф. ЛОКТЕВ

Казанский государственный энергетический университет

Рассматривается моделирование радиационных характеристик газовых компонентов продуктов сгорания и дисперсной фазы, представляющей собой сажевый аэрозоль, в математической модели замкнутого его радиационного теплообмена, предусматривающего предварительный расчет распределения температуры по объему топочной камеры. Острая селекция спектров поглощения и излучения учитывается применением двухпараметрического метода эквивалентной массы при вычислении функции спектрального пропускания оптически активных газовых компонентов продуктов сгорания топлива. Моделирование полидисперсного состава сажевого аэрозоля и его оптических характеристик выполняется с учетом различных механизмов его генерации.

Ключевые слова: радиационный теплообмен, перенос излучения, функция спектрального пропускания, продукты сгорания, дисперсная фаза.

Введение

Расчет радиационного теплообмена в топках энергетических котлов в настоящее время является актуальной задачей современной теплофизики в связи с созданием малогабаритных котлов с высокой производительностью, высокой теплонапряженностью топочного пространства, а также повышенным КПД, что требует оптимальных конструкторских решений при разработке или модернизации энергетических установок. Для сооружения и эксплуатации котельных агрегатов с высокими технико-экономическими показателями необходимо, чтобы камера сгорания (топка) удовлетворяла комплексному решению, отвечающему условиям оптимального способа сжигания топлива, интенсификации процесса теплопередачи к тепловоспринимающим поверхностям топки, рациональному распределению тепла по трубам экранов. В статье [1] рассмотрена математическая модель радиационного теплообмена в топках с учетом неравновесности и острой селекции спектров излучения в пламенных топках с применением замкнутой схемы моделирования, предусматривающей расчет распределения температуры по объему топочной камеры и химического состава продуктов сгорания топлива. Настоящая работа является продолжением работы [1] и рассматривает моделирование радиационных характеристик, входящих в алгоритмы радиационного теплообмена в объеме топки и на ее тепловоспринимающих поверхностях. Разработанное моделирование радиационного теплообмена базируется на дифференциально-интегральных соотношениях, записанных в алгоритмах функции спектрального пропускания (ФСП), что позволяет учитывать острую селекцию спектров излучения газовой фазы продуктов сгорания, неравновесность излучения зоны химических реакций при горении топлива и температурную неравновесность дисперсной фазы продуктов сгорания [2,4]. Основной

радиационной характеристикой газовой фазы продуктов сгорания является ФСП, а для дисперсной фазы - объемные спектральные коэффициенты поглощения и рассеяния, индикатрисы рассеяния по всему спектру теплового излучения, на моделировании которых остановимся ниже.

Функции спектрального пропускания (ФСП) паров Н20, С02 и малых компонентов продуктов сгорания

Разработанная методика применима для вычисления ФСП (основной радиационной характеристики) многокомпонентной неоднородной по температуре и эффективному давлению атмосферы дымовых газов продуктов сгорания в камере сгорания и газоходах котлов. Рабочий диапазон эффективных давлений 0,01 < Рэ < 100атм, что обеспечивает использование методики при решении задач

радиационного теплообмена как в современных котлах, так и в перспективных разработках энергетических и энерготехнологических агрегатов.

В общем случае при вычислении ФСП тду (где V - волновое число, Д -спектральное разрешение) необходимо выделить вклады в поглощения, обусловленные крыльями отдаленных спектральных линий атмосферных газов тДV , индуцированным

давлением поглощением тДV, селективным поглощением тДV спектральных линий, входящих в выбранный спектральный интервал. Тогда для заданного компонента

. (1)

Функция:

тД-у» • тДу = ехр [-[(Т) + руи(Т))шР], (2)

где Рук (Т) и Руи (Т)- коэффициенты континуального и индуцированного давлением поглощения(ИДП), зависящие от температуры Т; ш - содержание компонента; Р -парциальное давление. Для получения ФСП тДу предложено использовать общее соотношение:

1

1п ТДу

(

\

(

V

1 1

11п Т 'Ду ) 11п т''Ду]

М

(1П Т'Ду )(1П Т "Ду )

(3)

где

'Ду = ехр[-ку (Т )ш]

(4)

определяет ФСП в условиях слабого поглощения и при повышенных давлениях (Р > 10атм) в условиях смазанной вращательной структуры спектра поглощения;

,,С о ™ Шу,-» «V

Ду

ехр[-руС (Т )шШ Р«у ]

(5)

— ФСП при малых Рэ < 1атм в условиях сильного поглощения. Параметр М характеризует изменение скорости роста ФСП при переходе из области слабого поглощения в область сильного поглощения. Параметры ку, , , Рус определяются из экспериментальных данных, полученных с помощью измерительных комплексов. В соответствии с теорией модельного представления спектров поглощения ку = определяет отношение средней интенсивности к расстоянию между линиями, а

39

величина kvДv - характеризует интенсивность группы спектральных линий, расположенных в выбранном спектральном интервале Ду .

Было показано [2], что соотношение (3) описывает любую модельную структуру спектра, в том числе и закон Бугера для континуального спектра сильно перекрывающихся спектральных линий. Действительно, в этом случае m = 1, п = 0, РУС = ку, M = 1. Чем сильнее перекрывание спектральных линий, тем больше параметр m и тем меньше параметр п , и тем ближе параметр к единице. Для реальных спектров параметр M е{0,-1}. Континуальное поглощение крыльями линий и ИДП

описывается следующим набором параметров: m = 1, п = 1, к%, = РУ, M = -1.

Заметим, что спектры ИДП подчиняются другим правилам отбора, по сравнению с колебательно-вращательными спектрами, и полосы поглощения, запрещенные правилами отбора в колебательно-вращательных спектрах, становятся разрешенными в спектрах индуцированного давлением поглощения. В связи с этим учет ИДП может стать необходимым в радиационном теплообмене в силовых топках. В силовых топках более важен учет и континуального поглощения крыльями сильных линий, и полос поглощения.

С ростом температуры увеличивается плотность спектральных линий и,

следовательно, изменяются параметры mv, пл, ,Mv. В связи с этим при расчетах тДу в

условиях неоднородной по температуре и давлению среды используются средние значения этих параметров в определенном диапазоне температур.

Для расчета тДv в условиях неоднородной по температуре и давлению среды удобно ввести температурные функции:

Fl.iT) = ^ , F2c(T) = ^^7^ . (6)

Kv (Г)) Pvc(To)

Тогда:

- 1пт^ = ^с(То)т, - 1пТ"Дv=Pv(T)Wmlv , (7)

где

Пу 1

Wl = 1р(е)Ло[/(Т)]сИ, W2 = |р(е)т F¿У[I(Т)]Ш . (8)

e e V "0 у

Здесь эффективное давление:

N-1

Рэ = "N2 + В0 • Р02 + 2 ВгкРгк , (9)

2 г=1

где PN2 - давление N2, РО2 - давление О2, В^ - фактор уширения (отношение средней полуширины линий в выбранном интервале спектра для столкновений молекул г - к к средней полуширине спектральных линий в случае соударения молекул типа г с молекулами азота N2).

Аналогично для индуцированного и континуального поглощения:

Pvu (Т) = Р™ (70)^1 (Т), Р™ (Т) = Р™ (То№ (Т). (10)

40

Используемые для расчетов температурные функции: Еи (Т), Ек (Т), (Т), ^2с (Т) могут быть представлены в численной форме или в виде простых аналитических аппроксимаций, например в экспоненциально-степенной форме.

Экспериментально показано, что для многокомпонентной атмосферы полная ФСП определяется законом произведения функций по всем газовым компонентам:

т^у = Пт^ , (11)

г

где г - номер компонента.

Параметры Ку, РУ, шл,, пЛ!,М определены из измеренных спектров излучения и поглощения излучения высокотемпературными газовыми средами, моделированных с помощью подогревных кювет и пламенных измерительных комплексов, и численного моделирования тонкой структуры спектров [4]. На рис.1 иллюстрируется сохранение вращательной структуры спектра паров воды при высоких температурах [3].

Ми

Рис. 1. Восстановленный монохроматиеский спектр |1п| в области спектра 3064-3072см 1; 1 -полученный путем исключения аппаратурных искажений; 2 - измеренный спектр; 3 -индивидуальные линии. Температура пламени 2230K

Моделирование радиационных характеристик дисперсной фазы продуктов сгорания

Моделирование радиационных характеристик дисперсной фазы продуктов сгорания: спектральных коэффициентов поглощения ст^^, рассеяния а^х, ослабления стах - выполняется с учетом механизмов генерации золя в топочном объеме, процессов гетерогенной коагуляции первичного золя, седиментации золя в периферийных зонах топки [1,2,5]. Экспериментальные исследования микроструктуры дисперсной фазы продуктов сгорания [2,5] газового топлива в пламенной зоне оптическим методом дистанционной диагностики показали, что распределение числа частиц по размерам

^^ (где г - эффективный радиус частицы) сильно зависит от ингредиентного состава

дг

топлива и режима горения, которые определяют электрические свойства первичного золя и скорость его коагуляции. Первичный сажевый золь образуется в процессе горения углеводородного топлива в результате ионной нуклеации из газовой фазы и содержит около 10% растворимой в воде фракции углеводородно-кислородных ионов, тонкодисперсную нестабильную фракцию сажевого золя С2П, обладающую

41

свойствами газовой фазы и имеющую полосовую структуру электронного спектра. Частицы сажевого золя становятся устойчивыми, если их модальный радиус гт > 0,001мкм. Наиболее вероятный модальный радиус первичного золя составляет 0,003мкм [6,7].

Золь продуктов сгорания является полидисперсным и по выбранному нами моделированию представляется суперпозицией отдельных фракций, микроструктура N (г) которой определяется соотношением N

N (г) = 2 Щ (г), (12)

г=1

где г - номер фракции, N - число фракций. Как и в работах [2,3,5-7], спектральные коэффициенты ослабления, поглощения и рассеяния излучения золем будем определять с их нормировкой на длину волны ^ = 0,55мкм, а массовое содержание золя на оптическом пути будем выражать через оптическую толщину |- 1п т| , где т- ФСП на длине волны ^ = 0,55мкм. Математическое описание оптической плотности в пламенной зоне горения метана (дТ / д1) = 0,099мкм 1. При сжигании пропан-бутана в

воздухе наблюдается увеличение оптической плотности, так что = 0,2м 1.

Наиболее высокие значения I I = 0,4м наблюдаются в пламёнах при сжигании

ацетилена С2Н2, когда микроструктура сажевого золя наблюдается более грубодисперсной. При сжигании С2Н2 наблюдается и более высокое значение недогорания топлива. Микроструктура сажевого золя для различных компонентов газового топлива в пламенной зоне представлена в работах [2,6,7], а радиационные характеристики сажевого золя - в [5].

По разработанному нами моделированию предполагается, что спектральные радиационные характеристики для ансамбля частиц, генерируемых определенным механизмом образования золя, описываются эквивалентным ансамблем частиц сферической формы с комплексным показателем п = т - гх, где т их - его действительная и мнимая части, которые для многокомпонентного химического состава частиц описываются соотношениями [3]:

2 тгРг 2 ХгРг

т -; Х = -, (13)

2Рг 2Рг

г г

где Рг - массовая концентрация соединения г в частице. Суммирование распространяется на все соединения, входящие в состав частицы, в том числе и воздуха для рыхлых частиц, а также влаги, при учете воздействия конденсационных процессов на оптические характеристики дисперсной фазы в атмосферных выбросах продуктов сгорания. Спектральные радиационные характеристики сажевого золя для каждого компонента продукта сгорания вычислены по измеренной микроструктуре [5-7] и представлены в виде электронной базы данных. Далее предполагается, что процесс горения каждого компонента и образование сажевого золя происходит независимо, что позволяет использовать экспериментальные данные по оптической плотности для

42

нормировки радиационных характеристик. Вес каждой микроструктуры в оптической плотности сажевого золя определяется долей каждого компонента в массе сжимаемого топлива.

Микроструктура сажевого золя при движении продуктов сгорания изменяется в связи с выгоранием его тонкодисперсной фракции и захватом мелких частиц сажи крупными частицами в процессе коагуляции. Эти эффекты можно учесть путем изменения веса оптической плотности различных фракций сажевого золя по объему топки. Если В, (Ь) - оптическая плотность фракции г сажевого золя на длине волны

X = 0,55мкм в направлении Ь , то спектральные коэффициенты ослабления стах

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

А—

поглощения ста^

А—

, рассеяния ста^

(-

и индикатрисы рассеяния /,

V У ЧУ

^ - угол рассеяния, описываются соотношениями :

аБ сХа

(— N (— Ь

V у

= 1 В Ь с^ / X В,

г=1 V У г=1

N (— Ь

V У

Л

Ь, ^

V У

Ь

V У

где

(14)

( — \ N

А

ь, V

V У

= 1В

г=1

Ь

V У

Аг

Ь, V

N

/ X В,

г=1

(— Ь

V У

где г - номер компонента; N - число независимых механизмов генерации золя.

(15)

Особенности радиационного теплообмена в пламенных топках

В настоящее время разработанные различные методы моделирования переноса излучения в двухфазных средах, обзор которых дан в книге [8], не позволяют учитывать острую селекцию спектров излучения газовых компонентов продуктов сгорания, что приводит к значительному завышению, достигающему 100%, в потоках излучения на тепловоспринимающие поверхности. Расчеты температуры в топке, выполненные с использованием средних спектральных коэффициентов поглощения для оптически активных газовых компонентов с нормировкой результатов расчета на изменение энтальпии продуктов сгорания на выходе из топки, приводят к занижению температуры пламенной зоны. Только вычисление интенсивностей теплового излучения в алгоритмах ФСП двухпараметрическим методом эквивалентной массы позволяет правильно вычислить трехмерное поле температуры в топочном объеме и потоке теплового излучения на тепловоспринимающие поверхности. При этом необходимо учитывать неравновесность температуры дисперсной фазы в топочном объеме. Самой сложной в расчете радиационного теплообмена является задача вычисления стационарного трехмерного поля температуры в камере сгорания с использованием временной итерационной процедуры (соотношение (25) в [1]) и воздействий диффузии, турбулентной диффузии, естественной или принудительной рециркуляции продуктов сгорания в топочном объеме. Развитие и перемешивание струйных течений в топке может быть определено путем решения уравнения энергии, уравнения движения с учетом уравнений неразрывности и состояния газа [9,11] и привлечением теории подобия при инженерных расчетах [12] струйных течений.

Применение итерационных процедур в трехмерном варианте выполнения моделирования радиационного теплообмена в топочной камере требует задания

43

нулевого приближения поля температуры, химического состава продуктов сгорания и структурных характеристик дисперсной фазы, продуктов сгорания. Неравновесная температура дисперсной фазы вычисляется по радиационному балансу поглощенного в элементарном объеме излучения золем и собственным его излучением в пределах телесного угла 4п при неравновесной температуре Тн . Выполненный анализ влияния рассеяния излучения на радиационный теплообмен показал, что его влияние на потоки и притоки тепла не превосходит 2-3% при камерном сгорании газового топлива. Это позволяет при вычислении интенсивностей излучения использовать более простое соотношение (7) в [1], когда учитываются только излучающие свойства сажевого золя и оптически активных газовых компонентов.

Наиболее эффективна многокамерная конструкции топки или радиационного теплообменника при сжигании топлива в предтопках [2,5]. На рис.2 показана зависимость скорости радиационного теплообмена в пламенной зоне многокамерной топки с многорядным горелочным устройством матричного типа в зависимости от расстояния Ь между стенками топки. Из представленного рисунка видно, что скорость

дТ (Ь)

радиационного выхолаживания

дл

достигает максимального значения при малых

расстояниях Ь между стенками топки, когда ФСП т^ (Ь) стремится к единице. С

увеличением расстояния Ь, между стенками топки, скорость радиационного

дТ (Ь) Ч г й

уменьшается. Чем меньше величина Ь в многокамерной топке,

выхолаживания

дл

тем на меньшей высоте от пода топки достигается значение температуры Ту% на

выходе из топки. Детальные расчеты спектральных интенсивностей и пространственного распределения спектральных и интегральных потоков теплового излучения на тепловоспринимающих поверхностях ячеек многокамерной топки рассмотрены в [2,5].

эта-у гат(р\

Ш 7 V в* к

Рис.2. Относительная зависимость средней скорости радиационного выхолаживания дТ(Ь) К дТ (Ь) ^

—^^ / I —^^ I от расстояния Ь между стенками тепловоспринимающих поверхностей. дл / К дл /шах

Топливо - природный газ, подовое горелочное устройство матричного типа с многорядными горелками, коэффициент избытка воздуха а=1,03

44

В общем случае полная величина радиационного выхолаживания ДТ = Та - Тух зависит от расположения и функционирования горелок, формирующих

потоки пламени и продуктов сгорания и трехмерное распределение поля температуры в топке для стационарного режима работы. На рис.3 для примера иллюстрируются результаты расчетов вертикальных профилей скорости радиационного выхолаживания дТ(х)/дл, дТ(х)/дх и стационарное распределение температуры Т (г) от высоты топки

2 над срезами капилляров матричного горелочного устройства. Топливо - природный газ газопровода Шебаловка-Брянск-Москва, размер горизонтального сечения ячейки многокамерной топки 1,25х1,6 м2.

Скорость подачи продуктов сгорания на начальном участке топки составляет и>0 = 25м/с и ио = 20м/с при давлении в топке 1 • 105 Па. Высота пламенной зоны Дг = 0,7. В расчётах учтены равновесные и неравновесные процессы излучения по алгоритмам, рассмотренным выше. Предполагается, что процесс горения различных ингредиентов газового топлива происходит независимо при оптимальном значении коэффициента избытка воздуха а = 1,03 .

Микроструктура сажевого золя измерена при горении метана, пропан-бутана и ацетилена [2]. Оптические характеристики сажевого золя вычислены для измеренных микроструктур дисперсной фазы продуктов сгорания [13]. Объемные коэффициенты ослабления, поглощения и рассеяния нормированы на измеренные значения оптической плотности золя [2,5]. Ингредиентный состав продуктов сгорания вычислен по [13,14].

а) б)

Рис.3 Результаты расчёта радиационного теплообмена в многокамерной топке с размером горизонтального сечения ячейки 1,25х1,6 м2 для начальной средней скорости течения продуктов

сгорания 25 м/с (а) и 20 м/с (б). дТ(г)/дл, дТ(х)/дх - скорости радиационного выхолаживания, Т(х) -

температурный профиль средней по сечению температуры в зависимости от высоты ъ над срезами

капилляров многорядной горелки. 1 - дТ(х)/дл ; 2 - дТ(х)/дх для начальной средней скорости

течения продуктов сгорания 25 м/с; 3 - Т(ъ) для начальной средней скорости течения продуктов

сгорания 25 м/с; 4 - дТ(х)/дх для начальной средней скорости течения продуктов сгорания 20 м/с; 5

- Т(х) для начальной средней скорости течения продуктов сгорания 20 м/с

В расчетах учтено выгорание тонкодисперсной фракции золя и гетерогенная коагуляция тонкодисперсной фракции сажевого золя при движении продуктов сгорания.

45

Заключение

1. Оперативное моделирование радиационного теплообмена в топочном объеме и тепловосприятие тепловоспринимающими поверхностями требует одновременного расчета трехмерного поля температуры газовой и дисперсной фаз продуктов сгорания. Неравновесное радиационное выхолаживание в пламенной зоне горения топлива имеет противоположный характер для газовой и дисперсной фаз продуктов сгорания. Если неравновесность излучения газовой фазы повышает радиационное выхолаживание, то неравноверность температуры дисперсной фазы снижает радиационное выхолаживание пламенной зоны. В пламенной зоне температура дисперсной фазы Тд всегда ниже

температуры газовой фазы Тр, а величина ДТ = Тр - Тд может достигать 200К .

2. Интенсивности и потоки теплового излучения с достаточной точностью могут быть вычислены путем применения спектральных моделей расчета, записанных в алгоритмах ФСП с разделением на части за счет поглощения излучения аэрозолем и рассеяния излучения. ФСП тгх для оптически активного компонента г вычисляется двухпараметрическим методом эквивалентной массы в структурно неоднородных средах. Дисперсная фаза продуктов сгорания моделируется в виде суперпозиции нескольких фракций, каждая из которых формируется соответствующим механизмом его генерации. Неравновесность теплового излучения может быть учтена введением измеренных спектральных коэффициентов неравновесности или введением колебательной неравновесной температуры для функции источника излучения.

3. Наиболее эффективны конструкции многокамерных топок с восходящим движением продуктов сгорания в топке и вертикальным развитием пламени подовых многорядных горелок с общим газовым коллектором для подачи газового топлива и коллектора для подачи окислителя (воздуха или кислорода). Конструкция многокамерной топки позволяет повысить КПД на 2-3 % и увеличить паропроизводительность в 2-3 раза при сохранении параметров пара и габаритов котла.

4. Рассмотрение оптических свойств газовой и дисперсной фаз продуктов сгорания углеводородных топлив и предложенных алгоритмов численного моделирования позволяет сделать следующие выводы:

- неравновесное излучение снижает концентрацию NОх;

- неравновесное излучение приводит к нагреванию части топлива и ускоряет их воспламенение тем интенсивнее, чем мельче частицы и чем больше их сечение поглощения излучения;

- в случае сжигания газового топлива неравновесное излучение практически не трансформируется продуктами сгорания и без ослабления поглощается стенками камеры сгорания (экранами);

- в случае наличия дисперсной фазы неравновесное излучение поглощается аэрозолями (частицами сажи и топлива) и его роль в процессах радиационного теплообмена ослабляется, так как энергия неравновесного излучения в результате поглощения излучения переходит в тепловую энергию частиц. Основным компонентом, определяющим неравновесное излучение в пламёнах, является гидроксил ОН. Квантово-механическое рассмотрение формирования неравновесного излучения показывает, что неравновесное излучение проявляется как в электронных так и колебательно-вращательных спектрах основного и возбуждённого электронных состояний: полосы vl,2vl,3vl (где V! - частота нормального колебания). Неравновесное излучение ОН обнаружено в окрестности длин волн © Проблемы энергетики, 2012, № 3-4

46

1;1,43;2,1;2,7;4,1мкм водородно-кислородного пламени и при сжигании других органических топлив.

Summary

Modeling of radiate characteristics of gaze component of combustion products and disperse phase, representing sooty ashes, in a mathematic model of radiate heat exchange, reserved preliminary calculation of temperature distribution on volume of a fire chambers is considered. Sharp selection of radiate and absorption spectrum of gases is accounted by application of two-parametrical method of equalent mass on calculation of spectral transmission functions of optical active gaze components of combustion products of fuel. Modeling of polydispersion of composition sooty ashes with accounting of different mechanisms its generation.

Keywords: radiating heat exchange, radiation caring over, function spectral transmission, combustion products, a disperse phase.

Литература

1. Москаленко Н.И., Загидуллин Р.А., Зарипов А.В., Локтев Н.Ф., Моделирование радиационного теплообмена с учетом неравновесности и селективности излучения в топках энергетических котлов тепловых электрических станций. I. Математическая модель // Проблемы энергетики. 2012. №1-2. С. 23-33.

2. Moskalenko N.I., Zaripov A.V., Loktev N.F., Parzhin S.N., Zagidullin R.A. Transfer over nonequilibrium radiation in flames and high temperature mediums. - Optoelectronics-devices and applications. Croatia. 2009. P. 469-526.

3. Кондратьев К.Я., Москаленко Н.И., Поздняков Д.В. Атмосферный аэрозоль. Л.: Гидрометеоиздат, 1983. 224 с.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4.Москаленко Н.И., Локтев Н.Ф. Моделирование переноса селективного излучения в структурно-неоднордных средах // Тепловые процессы в технике. 2009. V. 1. №10. С. 432-435.

5.Москаленко Н.И., Зарипов А.В., Загидуллин Р.А. Спектры эмиссии и радиационный теплообмен в высокотемпературных средах, пламенах и топочных камерах. // - Техника и технологии в 21 веке: Современные состояния и перспективы развития. Новосибирск, ЦРНС. Т4. М. 48-87.

6.Москаленко Н.И., Локтев Н.Ф. Диагностика зольной фракции продуктов сгорания. / -Труды V-ой Рос. нац. конф. по теплообмену. М.:Изд.дом МЭИ, Т.8.

7.Москаленко Н.И., Локтев Н.Ф., Сафиуллина Я.С., Садыкова М.С. Идентификация ингредиентов и определение ингредиентного состава атмосферных выбросов и продуктов сгорания методом тонкоструйной спектрометрии // Междунар. научн. журнал «Альтернативная энергетика и экология». 2010. №2. С.43-54.

8. Сушкевич Т.А. Математические модели переноса излучения. М.: БИНОМ. Лаборатория знания. 2006. 664 с.

9.Сурженков С.Т. Тепловое излучение газов и плазмы. М.: Изд. МГТУ им. Н.Э. Баумана, 2004.

544с.

10. Вафин Д.Б. Дифференциальный метод теплового расчета топок: научное издание. Казань, 2008. 115с.

11. Ахмедов Р.Б. Основы регулирования топочных процессов. М.: Энергия,1977. 280с.

12. Абрамович Г.Н. Прикладная газовая динамика. М.: Наука, 1969. 824с.

13. Алемасов В.Е., Дрегалин А.Ф., Тишин А.П. Термодинамические и теплофизические свойства продуктов сгорания. М. ВИНИТИ, 1972. 490 с.

14. Алемасов В.Е., Дрегалин А.Ф., Крюков В.Г., Тишин А.П. Математическое моделирование высокотемпературных процессов в энергосиловых установках. М.: Наука, 1989. 254 с.

Поступила в редакцию 21 марта 2012 г.

Николай Иванович Москаленко - д-р физ.- мат. наук, профессор кафедры «Котельные установки и парогенераторы» (КУПГ) Казанского государственного энергетического университета (КГЭУ). Тел.: 8 (843) 234-64-06; 8 (843) 519-43-16.

47

Рустем Абдулхакович Загидуллин - аспирант кафедры «Котельные установки и парогенераторы» (КУПГ) Казанского государственного энергетического университета (КГЭУ). Тел.: 8 (843) 562-59-83; 8 (843) 519-43-16; 8 (917) 9195171.Б-шай: rz.13_.-@mail.ru.

Паржин Сергей Николаевич - канд. техн. наук, доцент кафедры «Котельные установки и парогенераторы» (КУПГ) Казанского государственного энергетического университета (КГЭУ). Тел.: 8 (843) 234-58-06; 8 (843) 519-43-16; 8 (917) 2457376.

Николай Федорович Локтев - аспирант кафедры «Котельные установки и парогенераторы» (КУПГ) Казанского государственного энергетического университета (КГЭУ). Тел.: 8 (843) 255-44-74; 8 (843) 519-43-16; 8(917) 9146768.Б-шаП: niklok@yandex.ru.

Зарипов Алмаз Вилевич - канд. техн. наук кафедры «Котельные установки и парогенераторы» (КУПГ) Казанского государственного энергетического университета (КГЭУ).

48

Modeling of radiate heat exchange with accounting of nonequilibrium and selection of radiation in fire chambers of energetic electric stations