Численная оценка улавливания шлака в топке открытого типа с жидким шлакоудалением

СПИСОК ЛИТЕРАТУРЫ

1. Гладышев В.П., Левицкая С.А., Филиппов Л.М. Аналитическая химия ртути. — М.: Наука, 1974. — 228 с.

2. Meyer S., Scholz F., Trittler R. Determination of inorganic ionic mercury down to 510—4 mol—1 by differential-pulse anodic stripping voltammetry // Fres. J. Anal. Chem. —1996. —V. 356. — № 2. — P. 247-252.

3. Захарчук Н.Ф., Илларионова И.С., Юделевич И.Г. Некоторые закономерности формирования и разрушения слоя с аномальными электрохимическими свойствами в системе C-Hg-Hg(II), HCl // Электрохимия. — 1982. — Т. 18. — № 3. — С. 331—338.

4. Выдра Ф., Штулик К., Юлакова Э. Инверсионная вольтампе-рометрия. — М.: Мир, 1980. — 278 с.

5. Брайнина Х.З., Нейман Е.Я. Твердофазные реакции в электро-аналитической химии. — М.: Химия, 1982. — 264 с.

6. Витер И.П., Каменев А.И. Определение компонентов системы Cd-Hg-Te методами инверсионной вольтамперометрии и хро-нопотенциометрии // Журн. аналит. химии. — 1998. — Т. 53. — № 11. —С. 1199-1204.

7. Надежина Л.С., Грихилес М.С., Демин В.А., Солодухина Е.В. Инверсионное вольтамперометрическое определение следов ртути в хлоридных растворах // Журн. аналит. химии. —1994. — Т 49. — № 9. — С. 974-980.

8. Захарова Э.А., Пичугина В.М., Толмачева Т.П. Определение ртути в водах и алкогольных напитках методом инверсионной вольтамперометрии // Журн. аналит. химии. — 1996. — Т. 51. — № 9. —С. 1000-1005.

9. Zen J-M., Chung M.J. Square-wave voltammetric stripping analysis of mercury (II) at a poly(4-vinylpyridine)/Gold film electrode // Anal. Qem. — 1995. —V. 67. — № 19. — Р. 3571—3577.

10. Salinas D.R., Cobo E.O., Garcia S.G., Bessone J.B. Early stages of mercury electrodeposition on HOPG // J. Electroanal. Chem. — 1999. —V. 470. — № 2. —P. 120-125.

11. Sahlin E., Jagner D., Ratana-Ohpas R. Mercury nucleation on glassy carbon electrodes // Anal. Chim. Acta. — 1997. — V. 346. — № 2. —P. 157-164.

12. Serruya A., Mostany J., Scharifker B.R. The kinetics of mercury nucleation from Hgf and Hg2+ solutions on vitreous carbon electrodes // J. Electroanal. Chem. —1999. —V. 464. — № 1. — P. 39 -47.

13. Гамбург Ю.Д. Электрохимическая кристаллизация металлов и сплавов. — М.: Янус-К, 1997. — 384 с.

14. Вол А.Е. Строение и свойства двойных металлических систем. — М.: Физматгиз, 1979. — 576 с.

15. Свидетельство на полезную модель № 12862 / Б.Ф. Назаров,

B.И. Чернов, Ю.А. Иванов. — Б. И. № 4. —2000.

16. Назаров Б.Ф., Виштакалюк Л.Д. Изучение механизма образования новой жидкой фазы на твердых электродах // Известия Томского политехнического института. —1976. —Т 302. —С. 18—19.

17. Мальков Е.М. Исследование по теории ртутно-графитного электрода в методе амальгамной полярографии с накоплением и его применение к анализу природных вод: Автореферат. — Дис. ... канд. хим. наук. —Томск: Изд-во ТГУ, 1970. — 24 с.

18. Хустенко Л. А., Ларина Л.Н., Назаров Б. Ф. Экспресс-определение ртути в природных водах методом инверсионной воль-тамперометрии на углеродном электроде, модифицированном золотом // Журн. аналит. химии. — 2003. — Т. 58. — № 3. —

C. 297—302.

19. Полинг Л., Полинг П. Химия. — М.: Мир, 1978. —683 с.

20. Энергии разрыва химических связей. Потенциалы ионизации и сродство к электрону / Под ред. акад. В.Н. Кондратьева. — М.: Наука, 1974. —351 с.

УДК 662.933.1:004.42

ЧИСЛЕННАЯ ОЦЕНКА УЛАВЛИВАНИЯ ШЛАКА В ТОПКЕ ОТКРЫТОГО ТИПА С ЖИДКИМ ШЛАКОУДАЛЕНИЕМ

А.В. Старченко, А.С. Заворин, С.В. Красильников

Томский политехнический университет E-mail: ghost@tpu.ru

Приводится описание математической модели и результаты численного расчета доли минеральной части угля, улавливаемой в топке, применительно к технологии факельного сжигания с удалением шлака в расплавленном состоянии. Рассмотрено влияние некоторых конструктивных и режимных факторов на эффективность шлакоулавливания.

Введение

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

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

ние условий шлакообразования и компьютерный расчет коэффициента шлакоулавливания являются ключевыми, поскольку от результата исследования зависит надежность моделирования взаимодействия минеральной составляющей топлива с поверхностями нагрева в топке и газовом тракте. Учитывая определяющую роль аэродинамики для организации пылеугольного факельного сжигания в энергетических котлах, необходимо исследовать шлакообразование в комплексе с моделированием двухфазных течений, тепло- и массообмена и горения. Для этого применен пакет прикладных программ FIRE 3D [1], развитый на основе програм-

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

Математическая модель

В основу модели положен метод, совмещающий эйлеров и лагранжев подходы для описания движения и теплообмена реагирующих газов и взвешенных частиц в топке парового котла [2]. При этом общие уравнения движения, теплообмена и горения в газовой фазе описаны на основе эйлерова способа представления, т.е. используются пространственные уравнения баланса массы, импульса, концентраций газовых компонентов и энергии для газовой смеси. Лагранжев подход применяется для описания движения и тепломассообмена одиночных частиц топлива вдоль их траекторий с учетом обратного влияния дисперсной фазы на несущую среду. Турбулентные характеристики газа рассчитываются с использованием двухпараметрической "k-s" модели турбулентности, также учитывающей влияние движущихся частиц в несущем потоке. Радиационный теплообмен представляется в рамках P1-приближения метода сферических гармоник, который показывает хорошие результаты применения к пылеугольным топкам [2].

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

Перемещаясь по топочному объему, угольные частицы нагреваются за счет радиационно-конвективного теплообмена, начинается выход летучих компонентов, их воспламенение и горение, а также догорание коксового остатка. Предполагается, что газовая среда в топке может состоять из химически инертных молекулярного азота N2, двуокиси углерода СО2, паров воды H2O, а также реагирующих О2, СО и летучих.

В этом случае математическая модель включает следующие блоки [4].

1) Газовая фаза

Уравнение неразрывности

_j + j

dx i vap char'

Уравнения баланса массы газовых компонентов

дрСщи, _ д f ^ dCN

dpCcoUi _ d ( U cC,

dx , dxi I Sc, dx,

dpCO U. d ( U dCO , O

У O2 1 O O2 _ n J _ в J _ J O2 .

dx. dx. I Sc, dx. J PvolJvo‘ PcOJ CO Jchar ;

dPCvoiUi _ d ( U dCv

dx.. dx. 1 Sc, dx. ¡''aoapJoap Joo‘

dx. dx. 1 Sc, dx,

dPCCO2Ui _ d ( U dC

dx,.

dx,. Sc, dx,

+ Jchar + (1 + Pco) JC

E C _ 1>

j_N ,O2 ,vol,CO,CO 2,H 2O

Уравнение движения

dpUUj

dx,

d P , d

dx, dx,

dU. dUj

dx, + dx,

+ Pg, +F, ;

j _ 1, 2, 3. Уравнение энергии

dpU ,cT d

U + U

Pr Pr

dx, dx,

Уравнение состояния p=

C dT ) + Q'o,J'o‘ + QcOJco dx,

+ 0.

Co

Cn

Co.

Cc

Cc

Ch

Mo Mn

Moo, M CO M CO2 M H2O

Здесь р-плотность газовой смеси; U (/=1,2,3) — компоненты скорости газа; x (/=1,2,3) — декартовы координаты; Jvap, Jchar — массовые скорости выхода летучих при пиролизе угольного топлива и догорания образовавшегося коксового остатка (гетерогенные реакции); по повторяющимся индексам производится суммирование от 1 до 3 (соглашение Эйнштейна); Cj — концентрации газовых компонентов; л, Sc,—турбулентная вязкость газа и турбулентное число Шмидта (Sc, =0,9); Jvoi — массовая скорость горения летучих, зависящая от скорости химической реакции горения углеводородов CMHN—Jvolick и интенсивности турбулентного смешения реагирующих компонентов Jvoiidif, [3wl — количество O2, необходимое для сгорания 1 кг летучих; к, s — энергия турбулентности газа и скорость ее диссипации; /со — массовая скорость реакции горения CO+1/2O2=CO2; q[ad — компоненты вектора радиационного теплового потока; 0 — интенсивность теплообмена между несущей средой и частицами. В уравнении энергии моделируется перенос тепла за счет конвекции, термического излучения, поступления тепла вследствие горения летучих и CO, а также обмен теплом с дисперсной фазой; M — молекулярные веса компонентов газовой смеси; R — универсальная газовая постоянная.

2) Дисперсная фаза

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

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

Уравнения движения частицы-маркера /-ой группы [4]

йП1 и -и‘

— + gl; ] = 1,2,3,

йг

3свр\и‘-и\

\ ^(1 + 0,15 Ие,0’687 ); I 0,44;

Яе, < 1000;

Яе, > 1000;

Ие.

Е и  -и  )2 й,

=1

\и‘-и\

булентности.

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

Уравнения тепломассообмена и горения частиц

^ = -м "ар -мс1" йг ' ' ’

~асоп, (Т - Т ) + "

+8р

йТ.

тср = ОсЛат4,4

-пй2

-ОрМ

Здесь - масса одиночной частицы с диаметром ф; ср - теплоемкость угля; <2скар - теплота сгорания кокса; <2шр -тепло, необходимое для выхода летучих; бр - степень черноты частицы; асою - коэффициент конвективной теплоотдачи одиночной частицы (асжуф)/А=2+0,459Ке/!,'55Рг0'33); Я=^(с/Рг) - коэффициент теплопроводности газа; Т - температура /-ой частицы угля.

При определении массовой скорости выхода летучих ИГ используется соотношение (закон Аррениуса),

М7 = Кр ехр I -

Ер

КТ,

Лот

где ты - начальная масса сухой частицы с диаметром ф; - массовая доля летучих в сухом угле;

куар - предэкспонент; Бтр - энергия активации.

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

М сНаг = пй2.рС° КкКв

КК = ксНаг еХР| -

^КХ^В

°2 КК + К

ЕсИаг

В

КТ,

где Щ - компоненты скорости частиц; Ц=Ц+и) -компоненты актуальной скорости газа; рр - плотность угля; ф - размер частиц /-ой группы; св. - коэффициент сопротивления частиц; т - время их релаксации; турбулентные пульсации скорости газа и' выбираются случайным образом на основе принятого нормального распределения Гаусса со среднеквадратичным отклонением, равным 2/3к,

т.е. и' = ^^|к,|^2|=1 к - энергия турбулентности.

Значение ^ изменяется с помощью датчика случайных чисел в процессе движения частицы в исследуемом пространстве через отрезок времени

где ь = - масштаб тур-

Б

Для Хв используется зависимость

КВ = ^(2,0 + 0,16 Яе,0’667);

В° = 16-Ю-6^-^)1'9.

°2 ч273^

3) Влияние дисперсной фазы на газовую

В приведенной выше математической модели необходимо определить члены в уравнениях переноса, которые представляют влияние частиц на характеристики несущей среды: 1Лаг, ¥ &, 8к, Бн.

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

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

Бк, Бн по следующим формулам [4].

, Ат”р

¥,,

¥,

Jcн

=-V Е

и‘-и,

Ат,

то,и‘

Ат

= 1,2,3;

= -V Е[&*» Атснаг +тгср АТ О АтГР ] т -

- VЕ сТ + 1Еи-и, )2 т„,+ VЕт-Еи - и,)

-БР

т„г,;

Р

Т

mt

t, =

0,41k

/=1 i

N

Sh = f M(4*T4 -H)ms

Здесь SH - поглощение излучения топливными частицами, Sk - дополнительная диссипация турбулентной энергии, обусловленная влиянием дисперсной фазы, N - число частиц-маркеров, прошедших через рассматриваемый контрольный объем, имеющий размер V; Am=Amvap+Amf'ar - изменение массы /-ой частицы в процессе выхода летучих и горения коксового остатка за время прохождения ею данного контрольного объема; md - начальная масса частицы-маркера; moi и m - начальный и текущий массовые расходы частиц данной фракции; ti - время прохода i-ой частицы через контрольный объем V; - - средняя масса i-ой частицы в контрольном объеме; AT - изменение температуры частицы в контрольном объеме V.

4) Шлаковая пленка

Учет образования и движения шлаковой пленки осуществляется с использованием решения Маршака [5], распространенного на случай пространственных топочных процессов. Согласно предложенному подходу температура поверхности шлаковой пленки определяется из уравнения теплового потока на поверхности

qn =оет(Тф -К) + а(Тф -Тп),

где а - постоянная Стефана-Больцмана, sT - степень черноты топки, T0 - локальная температура в топке, а - коэффициент теплоотдачи, ТП - температура поверхности движущейся шлаковой пленки.

Решение задачи об определении средней по сечению пленки скорости имеет следующий вид [5]

рш8^п ____________ш

Нп

Но

ln Н

Нп

РшЕ51

Нп

Нп

Нп

+21 -Î2L-11-2ln^ -

Нп

ln Е

Нп

(Рш^хЗп )Dow„ - Но | lnН

Рш

Рш81

Нп

+2g I Е -11-2g ln Е - g f ln Е

Нп

Е

X

здесь рШ - плотность шлаковой пленки, - относи-

тельная вязкость шлака на границе перехода в пластичное состояние, ¡лП - вязкость шлака при

температуре Гд: и-п = ^(Тп) = №> ; ^(Т)= Tg -ка-

сательное напряжение на границе шлаковой пленки.

С другой стороны, уравнение изменения массы шлаковой пленки на участке " Up-Down" (верх-низ ячейки) в интегральной форме записывается:

где SW- интенсивность осаждения шлаковых частиц на стенку или в пленку из топочного пространства на рассматриваемом участке поверхности топки. Из (1) и (2) получим:

Таким образом, порядок расчета температуры шлаковой пленки следующий. Начиная с верхних сечений топки, где не образуется шлаковая пленка, из (2) получаем ¿П=0 (т.к. SW=0). Если же осаждение частиц есть, то расчет производится по вышеприведенным формулам. Заметим, что формулу (3) легко можно распространить на случай наклонной поверхности, для чего следует заменить g на gsinа, где а - острый угол между горизонтом и поверхностью.

Результаты вычислительных экспериментов

Расчеты выполнены для топки с жидким шлакоудал ением котла БКЗ-220-100ЖШ, по которой имеется обширный банк данных, полученных в ходе натурных экспериментов и физического моделирования аэродинамики [6]. Топка открытого типа с прямым вдуванием пыли имеет четыре плоскоструйные горелки, установленные на фронтовой стене с наклоном к горизонту и разворотом к продольной оси котла (на 15° для крайних и на 7° для средних горелок). За счет взаимной ориентации горелок и стен в нижней части топки создается своеобразное взаимодействие плоских струй, приводящее к возникновению в зоне активного горения горизонтальной вихревой структуры факела, названной разработчиками аэрошатровой .

Исходя из вариативности воспроизведения условий работы топочной камеры, вычислительная оценка осаждения шлака выполнена несколькими этапами: 1) исследовано влияние перераспределения угольной пыли по горелкам, а также изменение вертикального угла оси горелок; 2) проведен анализ роли дисперсного состава топлива; 3) исследовано влияние изменения геометрии топки путем установки дополнительного пережима разной глубины напротив верхнего среза горелок. При этом базовым для сравнения вариантов принят режим несения полной нагрузки при работе всех четырех горелок (сочетание 1-2-3-4), наклоне горелок на 15° и при одноплоскостном исполнении противоположной горелкам стены топки. Оценка коэффициента шлакоулавливания произведена по осаждению частиц в пределах футерованных экранов, ограничивающих объем камеры сгорания.

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

п

ры. Трехмерная визуализация распределения шлако-осаждения по стенам при различных вариантах сочетания работающих горелок показана на рис. 1.

Можно видеть, что наибольшая локализация осаждения шлака в камере сгорания достигается при рассредоточенном вводе угольной пыли в топку (в ра->0,0042

боте 3 или 4 горелки). В базовом режиме (4 горелки) расчетный коэффициент шлакоулавливания составляет 0,258. При обеспечении нагрузки котла тремя горелками при варианте с двумя крайними и одной средней (сочетание 1-2-4) коэффициент шлакоулав-ливания уменьшается до 0,245 за счет увеличения вы-

0,004

0,0033

0,0038

0,0034

0,0032

0,003

0,0023

0,0026

0.0024

0,0022

0,002

0,0013

0,001В

0,0014

0.0012

0,001

0,0003

0,0005

0.0004

0,0002

- 0

Горелки 1

Горелки 1-2

Г орелки 1 -2-3

Режим показа:

Т олщина пленки, м

Горелки 1-2-3-4

Горелки 1-2-4

Горелки 1-3

-Л

Горелки 1-4

Горелки 2

Горелки 2-3

Рис. 1. Визуализация распределения шлакоосаждения в топке котла БКЗ-220-100 ЖШ при разных сочетаниях задействованных горелок

носа частиц из камеры сгорания через пространство вблизи неработающей горелки. Вариант с большей асимметрией (сочетание 1-2-3) увеличивает коэффициент шлакоулавливания до 0,346 за счет концентрирования области аэродинамического взаимодействия факела с тыльным экраном с соответствующим увеличением толщины шлаковой пленки.

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

Наименьший коэффициент шлакоулавливания имеет режим работы с крайними горелками (сочетание 1-4), однако за счет удлинения траектории факела обеспечивается более полное выгорание топлива. По сравнению с этим при прочих равных условиях

коэффициент шлакоулавливания возрастает более чем в два раза при работе с одиночной горелкой или с двумя центральными горелками (сочетание 2-3).

Влияние угла наклона горелок вниз от горизонтали в визуализированном представлении результатов сказывается на сокращении зоны осаждения и смещении ее в подовую часть топки (рис. 2).

В целом это сопровождается тенденцией увеличения коэффициента шлакоулавливания (рис. 3, а). Однако в области угла наклона 20° и более происходит повышение сепарации недогоревших коксовых частиц на под топки с внедрением их в шлаковую пленку. На практике это приводит к развитию восстановления железа в расплаве шлака и "металлизации" пода с угрозой взрывоопасности системы механизированного удаления шлака. Уменьшение угла наклона до 10° также несколько увеличивает коэффициент шлакоулавливания за счет появления эффекта "наброса" на противоположный экран при той же дальнобойности горелочных струй.

0,0085

0,008

..0075

0,007

0,0005

0.006

0,0055

0,005

0,0045

0,004

0.0035

0,003

0,0025

0,002

0,0015

0.001

0,0005

0

Режим показа:

Т олщина пленки, м

10 градусов

15 градусов

20 градусов

25 градусов

Л

Рис. 2. Визуализация распределения шлакоосаждения в топке котла БКЗ-220-100 ЖШ при разных углах наклона горелок к горизонту

а)

б)

в)

Рис. 3. Изменение коэффициента шлакоулавливания (%) от конструктивных и режимных факторов топочного процесса: а) от угла наклона горелок; б) от тонкости помола угольной пыли (по остатку на сите R9G); в) от глубины пережима на противоположном горелкам экране

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

эффициента шлакоулавливания от 0,235 до 0,304 (рис. 3, б).

В связи с обсуждением целесообразности существенного изменения принципов организации аэродинамики топки путем реконструкции противоположного горелкам экрана с обустройством и нахождением оптимальной глубины аэродинамического пережима [7] показано (рис. 3, в), что увеличение его размера приводит к уменьшению коэффициента шлакоулавливания вдвое (от 0,258 до 0,125).

Лишь при перекрытии сечения топки пережимом на 30 % намечается последующее повышение коэффициента шлакоулавливания до 0,147. Выявленные эффекты связаны со сложной аэродинамикой, формирующейся вблизи пережима при различных его размерах. В частности, дополнительный пережим разбивает горелочные струи на два потока, один из которых уходит вверх над пережимом, другой - в нижнюю часть топки. При увеличении глубины пережима над ним вплоть до выходного окна топки отчетливо проявляется обширная рециркуляционная зона, осаждение частиц на стенки камеры сгорания уменьшается, зона повышенных температур (зона горения) смещается к низу топки, наблюдается рост скорости газодисперсного потока над подом топки. Установка дополнительного пережима уменьшает общий объем зоны реагирования горючей смеси, но увеличивает теплообменную поверхность, и при меняющейся аэродинамической картине это приводит лишь к незначительному уменьшению температуры на выходе из топки.

Заключение

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

СПИСОК ЛИТЕРАТУРЫ

1. Старченко A^., Заворин A.C., Красильников С.В. Применение пакета FIRE 3D к анализу процессов шлакоулавливания в пылеугольных топках // Известия Томского политехнического университета. -2002. -Т 305. - N 2. -С. 152-157.

2. Бубенчиков A.M., Старченко A^. Численные модели динамики и горения аэродисперсных смесей в каналах. - Томск: Изд-во ТГУ, 1998. -236 с.

3. Заворин A.C, Старченко A^., Красильников С.В. Оценка степени шлакоулавливания пылеугольных топок // Сопряженные задачи механики и экологии: Матер. Междунар. конф. — Томск: Изд-во ТГУ, 1998. - С. 109-110.

4. Starchenko A.V., Krasilnikov S.V. A Numerical Investigation of Particle Deposition in Pulverized Coal-Fired Furnaces // Two-Phase Flow Predictions. Proc. of the 10th Workshop. Merseburg, April 9-12, 2002. -P. 237-246.

5. Маршак Ю.Л., Рыжаков А.В. Шиповые экраны топок паровых котлов. - М.: Энергия, 1969. -240 с.

6. Федецкий И.И. Исследование процессов в минеральной части назаровского угля в топке с плоскими параллельными струями. -Дис. ... канд. техн. наук. - Томск: ТПИ, 1980. -255 с.

7. Федецкий И.И. Модернизация аэродинамической структуры факела и топочного процесса котла БКЗ-220-100 ЖШ // Исследование и конструирование котлов: Межвузовский научнотехнический сборник. -Томск: Изд-во ТПУ, 1993. -С. 39—45.