Численное исследование пространственных двухфазных течений и горения в пылеугольной топке с учетом шлакоулавливания тема диссертации и автореферата по ВАК РФ 05.14.14, кандидат технических наук Красильников, Сергей Витальевич

• Красильников, Сергей Витальевич
• кандидат технических наук
• 2003

Оглавление диссертации кандидат технических наук Красильников, Сергей Витальевич

Список сокращений и обозначений.

Введение.

Глава 1. Современное представление и уровень моделирования топочных процессов.

1.1. Обзор работ по моделированию трехмерных течений и теплообмена в топках.

1.2. Свойства минеральных компонентов пылевидного топлива как фактор шлакования.

1.3. Основные задачи моделирования.

Глава 2. Математическая модель пространственного двухфазного турбулентного потока в камере сгорания с учетом шлакования.

2.1. Описание аэродинамики, теплообмена и горения в газовой фазе.

2.2. Модель турбулентности.

2.3. Представление движения и теплообмена полидисперсных частиц в топке.

2.4. Учет влияния дисперсной фазы на несущую.

2.5. Формулировка граничных условий (в том числе процесса шлакоулавливания).

2.6. Осаждение дисперсной фазы на стенки. Модифицированная модель Маршака.

2.7. Выводы.

Глава 3. Численный метод расчета двухфазных турбулентных потоков в областях сложной геометрии.

3.1. Уравнение переноса для обобщенной переменной.

3.1.1. Выбор конечно-разностной сетки.

3.1.2. Получение дискретных аналогов дифференциальных уравнений.

3.1.3. Аппроксимация конвективных членов уравнений.

3.1.4. Расчет поля давления. Выбор итерационного метода для решения сеточного уравнения Пуассона.

3.2. Метод расчета движения и тепломассообмена топливных частиц и их воздействия на движение несущей среды.

3.3. Тестирование численной модели и метода расчета.

3.4. Выводы.

Глава 4. Исследование влияния условий подачи топливо-воздушной смеси, дисперсного состава топлива и геометрии топки котла БКЗ-220-100ЖШ на аэродинамику,

• теплообмен, горение и шлакоулавливание.

4.1. Геометрические характеристики и условия ввода топлива и воздуха.

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

4.2.1. Аэродинамика топки и осаждение шлака в топочной камере при изменении распределения нагрузки топлива на горелочные устройства.

4.2.2. Аэродинамика и осаждение шлака в топочной камере при разном угле наклона горелок по отношению к горизонту.

4.3. Влияние дисперсного состава топлива на интенсивность

• осаждения золовых частиц на стенах топки.

4.4. Влияние изменения геометрии топки (установки дополнительного пережима) на шлакоулавливание.

4.5. Выводы

Глава 5. Пакет прикладных программ Fire 3D для численного моделирования и визуализации результатов расчетов пространственной аэродинамики в пылеугольных топках.

5.1. Общая характеристика пакета прикладных программ Fire 3D.

5.2. Режимы организации вычислений.

5.3. Графическое представление результатов.

5.4. Выводы.

Введение диссертации (часть автореферата) на тему «Численное исследование пространственных двухфазных течений и горения в пылеугольной топке с учетом шлакоулавливания»

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

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

Работа выполнялась в соответствии с основными направлениями НИР Томского политехнического университета в рамках темы 6.71.02 ЕЗН Минобразования РФ и в продолжение исследований, выполненных ранее в ТПУ в соответствии с научно-техническими программами «Исследование и освоение сжигания канско-ачинских углей на электростанциях КАТЭКа на 1981-1985 годы» ГКНТ при СМ СССР, «Сибирь» СО АН СССР, государственной НТП «Экологически чистая энергетика».

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

Для достижения данной цели сформулированы следующие основные задачи исследования:

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

• применение новейших вычислительных схем для повышения точности результатов и сокращения времени счета;

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

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

• выявление возможностей управляющего воздействия режимно-конструктивных факторов на распределение минеральной части пылеугольного факела между шлаком и уносом на примере экспериментально [45] изученной топки;

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

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

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

Разработанный программный комплекс FIRE 3D передан для использования в ООО «Сибтерм», в специализированную научно-исследовательскую организацию по обследованию объектов котлонадзора ООО «Теплоуниверсал». Результаты вариативных расчетов для котла БКЗ-220-100ЖШ используются ТЭЦ ОАО «Юргинский машиностроительный завод» для анализа эксплуатационных режимов и выработки наладочных мероприятий. Материалы выполненных исследований включены в лекционный курс и лабораторный практикум дисциплины учебного плана для специальности 101300 (котло- и реакторо-строение) в Томском политехническом университете. Пакет программ опробован в НИРС кафедры ПГС и ПГУ Томского политехнического университета и на занятиях по вычислительным методам студентами и преподавателями кафедры.

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

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

На защиту выносится:

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

• результаты моделирования поведения газодисперсной смеси в топках с жидким шлакоудалением на примере котла БКЗ 220-100ЖШ.

Основные положения диссертационной работы докладывались и обсуждались на международной конференции «Сопряженные задачи механики и экологии» (Томск, 1998 г.), IV научно-техническом семинаре и V, VI Всероссийской научно-технической конференции «Энергетика: экология, надежность, безопасность» (Томск, 1998, 1999, 2000 г.г.), V областной научно-практической конференции студентов, аспирантов и молодых ученых «Современные техника и технологии» (Томск, 1999 г.), Всероссийской научно-практической конференции «Проблемы использования канско-ачинских углей на электростанциях» (Красноярск, 2000 г.), II семинаре вузов Сибири и Дальнего Востока по теплофизике и теплоэнергетике (Томск, 2001 г.), 10 научной конференции по двухфазным течениям «10th Workshop on Two-Phase Flow Predictions» (Merseburg, Германия, 2002 г.), XXVI Сибирском теплофизическом семинаре (Новосибирск, 2002 г), на научных семинарах кафедры парогенераторостроения и паро-генераторных установок Томского политехнического университета (1998 -2002 г.г.).

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

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

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

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

В четвертой главе приведены результаты исследования влияния условий подачи топлива, полученные на основе численного моделирования топки котла БКЗ-220-1ООЖШ при различных ее геометрических характеристиках, аэродинамики вдува и конструктивных особенностях этого котлоагрегата.

В пятой главе описан собственно пакет прикладных программ, реализующих представленный в этой работе метод расчета и способ визуализации результирующих данных на основе COM (Component Object Model) DirectX Windows. Описан графический интерфейс и режимы организации вычислений.

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

Автором выполнены работы по созданию программного комплекса для расчета объемных топочных процессов на основе пакета прикладного программирования DELPHI. Разработаны методики проведения вычислений и обработки полученных данных, проведены вычислительные эксперименты и анализ их результатов. Разработка вычислительной задачи, постановка задач исследований, обсуждение методики вычислительных экспериментов и полученных результатов выполнены с участием научных руководителей к. т. н. Заворина А.С. и д. ф.- м. н. Старченко А.В.

Заключение диссертации по теме «Тепловые электрические станции, их энергетические системы и агрегаты», Красильников, Сергей Витальевич

5.4. Выводы

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

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

ЗАКЛЮЧЕНИЕ

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

Численное интегрирование системы дифференциальных уравнений сформулированной математической модели аэродинамики, горения и теплообмена в пылеугольной топке выполнено с использованием хорошо зарекомендовавших себя в практических расчетах метода контрольного объема и алгоритма SIMPLE Патанкара и Сполдинга. При аппроксимации конвективных членов уравнений применены схемы MLU и QUICK с минимальной схемной вязкостью, которые дают близкую скорость сходимости глобальной итерационной процедуры и неплохое воспроизведение экспериментальных данных. На основе проведенных сравнительных расчетов было установлено, что для решения сеточных уравнений для поправки давления лучше использовать итерационные методы сопряженных невязок MRES и метод Ван дер Ворста с предобуславли-вающей матрицей BICGSTAB, которые по скорости сходимости превосходят обычный явный метод Булеева на 30-50 процентов.

Тестирование численного метода на результатах холодных продувок показало, что выбранная вычислительная процедура удовлетворительно предсказывает сложные пространственные течения и показывает неплохую сходимость расчетов с данными визуализации, полученными в ИТФ СО РАН и ТПУ. Применение разработанной численной модели к расчету реальных процессов в топках с жидким (БКЗ 220-100ЖШ) и твердым (БКЗ-210-140Ф) шлакоудалением показало хорошее соответствие предсказанных и измеренных температур топочной среды по высоте камеры сгорания.

Для котла БКЗ-220-100ЖШ исследовано влияние условий подачи топли-вовоздушной смеси, дисперсного состава пылеугольного топлива и геометрии топочной камеры на аэродинамику, теплообмен, горение и интенсивность шла-коулавливания в топке. Оценка роли перераспределения топлива по горелкам показала, что этот фактор в значительной степени связан с длиной факела. Чем длиннее реагирующая газодисперсная горелочная струя, тем лучше в ней прогорает топливо и меньше инерционных частиц осаждается на стенки. Соответственно уменьшается недожог и увеличивается количество прореагировавшего кислорода. Наиболее подходящими в этом смысле условиями подачи топлива являются ситуации, когда работают крайние горелки топки.

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

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

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

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

Программный комплекс FIRE 3D используется в конструкторской деятельности в ООО «Сибтерм», специализированной научно-исследовательской организации по обследованию объектов котлонадзора ООО «Теплоуниверсал», а также в дисциплинах учебного плана студентов ТПУ по специальности «Кот-ло-и реакторостроение».

На основании полученных результатов можно сделать следующие основные выводы:

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

2. Тестовые расчеты, выполненные с использованием разработанной численной модели для гидромоделей и натурных пылеугольных топок (БКЗ-220-100ЖШ, БКЗ-210-140Ф), показали хорошее воспроизведение результатов экспериментальных исследований.

3. Проведенные вычислительные эксперименты по расчету топочных процессов в топке парового котла с жидким шлакоудалением БКЗ-220-100ЖШ подтвердили правильность принятого конструкторского решения по уменьшению угла наклона горелочных струй с 30° до 15°, следствием которого стало существенное снижение сепарации коксовых частиц в шлаковую пленку с последующей металлизацией пода топки.

4. Анализ влияния условий подачи топливо-воздушной смеси, дисперсного состава и геометрии топочной камеры котла БКЗ-220- 100-ЖШ показал возможность изменения коэффициента шлакоулавливания от 12% до 66%. Определено, что чем меньше дисперсная фаза встречает на своем пути полета препятствий и вихрей, и тоньше помол топлива, тем меньше коэффициент шлакоулавливания, больше коэффициент уноса и, соответственно, больше загрязнений поверхностей нагрева за пределами зоны футерованных экранов.

5. На основании вычислительных экспериментов можно сделать вывод о возможностях управляющего воздействия на поведение топливовоздушной смеси и распределение минеральной части между шлаком и уносом при помощи конструктивно-наладочных факторов. Так, при перенаправлении топлива в разные горелки шлакоулавливание изменяется от 22% до 66%, при изменении угла наклона горелок - от 25%о до 29%, тонкости помола от 23% до 31% и установке дополнительного пережима от 12% до 26%.

6. На основе использованной вычислительной модели создан пакет прикладных программ FIRE 3D, ориентированный на персональные компьютеры, с дружественным пользовательским интерфейсом, удобной визуализацией пространственных полей параметров топочных процессов.

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

1. Сполдинг Д.Б. Горение и массообмен. М.: Машиностроение, 1985. —240 с.

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

3. Лебедев А.Н., Шейнин Г.А., Шлакование пылеугольных топок и борьба с ним. M.-J1.: Энергия, 1966. -113 с.

4. Шагалова С.Л., Шницер И.Н. Сжигание твердого топлива в топках парогенераторов. ~Л.: Энергия, 1976.-176 с.

5. Померанцев В.В., Шагалова С.Л., Арефьев К.М. Приближенная методика расчета выгорания пылеугольного факела. // Теплоэнергетика. 1958. -№11. - С. 33-41.

6. Довжик Б. В., Мигай В. К. Численное моделирование трехмерных вязких течений в топочных камерах //ИФЖ. 1988. - Т. 55, № 1. - С. 42-50.

7. Егоров Ю. Э., Колешко С. Б. Применение метода дробных шагов для численного решения уравнений несэ/симаемой вязкой жидкости в естественных переменных В кн.: Динамика неоднородных и сжимаемых сред. - Л.: Изд-во ЛГУ, 1984. -С. 80- 92.

8. Дектерев А. А., Ковалевский А. М. Программный комплекс МАТОК по расчету аэродинамики топочных камер энергетических котлов // Сибирский физико-технический журнал. 1992. - №6. - С. 146-152.

9. Патанкар С. Численные методы решения задач теплообмена и динамики жидкости. —М.: Энергоатомиздат, 1984. 152 с.

10. Leonard B.P.A. Stable and accurate convective modeling procedure based on quadratic upstream interpolation // Computational Methods in Applied Mechanics and Engineering. 1979. - v. 19, 4. - P. 59-98

11. Hirt C. W., Cook J. L. Calculating three-dimensional flows around structures and over rough terrain //J. of Comput. Phys.-1972.-Vol.10.

12. Белоцерковский O.M., Давыдов Ю.М. Метод крупных частиц в газовой динамике. М.: Наука. — 1982. — 392 с.

13. Amsden А.А., Harlow F.H. The SMAC method: a numerical technique for calculating incompressible fluid flows. — Los Alamos Scientific Lab., Kept. N LA-4370. -Los Alamos: 1970.

14. Закрученные потоки: Пер. с англ./Гупта А., Лилли Д., Сайред Н. — М.: Мир, 1987.- 588 с.

15. Старченко А. В., Федецкий И. И. Численное моделирование трехмерных течений в топках паровых котлов // Сибирский физико-технический журнал. 1993. -№ 5.-С. 118- 124.

16. Launder В.Е., Spalding D.B. The numerical computation of turbulent flows // Computational Methods in Applied Mechanics and Engineering. 1974. - Vol.3, N.2.-P.269- 289.

17. Ильин В.П., Юдин А.Н. Решение трехмерных разностных уравнений методом Булеева с сопряженными градиентами //Технология моделирования задач математической физики. 1989. - Новосибирск: ВЦ СО АНСССР. - С. 152-165.

18. Заворин А.С., Старченко А.В., Красилъников С.В., Оценка степени шлакоулавливания пылеуголъных топок //Сопряженные задачи механики и экологии: Материалы международной конференции. — Томск: Изд. Том. ун-та., 1998. С. 109-110.

19. Бубенчиков A.M., Старченко А.В., Ушаков В.М. Численное моделирование горения пылеугольного топлива в топках паровых котлов //Физика горения и взрыва. 1995. - Т. 31, № 2. - С. 23-31.

20. Старченко А.В. Математическое моделирование образования оксидов азота при горении пылеугольного топлива //Физика горения и взрыва. -1998. Т.34, № 6. С. 3-13.

21. Fischer К., Leithner R., Muller H. Three-dimensional of the gas-solid flow in coal-dust fired furnace //Proceeding of the First Int. Symposium on Two-Phase Flow Modelling and Experimentation. Rome. Italy. 1995. - V. 1. — P. 1387 - 1393.

22. Файвленд, Вассел. Численная трехмерная модель для расчета характеристик пылеугольной топки //Современное машиностроение. Сер. А. 1989. -Ml.-С. 158-169.

23. Emami M.D., Lazopulos G., Lockwood F.C. The computation of heat transfer in engineering combustion equipment //Proc. of the 37 EUROTHERM Seminar. -Sallugia: ENEA. 1994. - P. 83-105.

24. Broughton R.C., Oliver A.J., A numerical method for predicting three-dimensoinal flows and its application to boiler furnaces, 2nd UK Nat. Conf Heat Transfer, Glasgow, 14-16 Sept., 1988. Vol. 2. Sess 4A-6C. -London, 1988.-P. 12331242.

25. Leithner R., Doring M., Dampferzeugersimulation unter besonderer Berucksichtigung des Auskuhlverhaltens. IWBT-Tag, 30. November 1993 Tagungs-band (Teil II).

26. Lockwood F.C., Shah N.G. Evaluation of an efficient radiation flux model for furnace prediction procedures // Proceedings of Sixth Int. Heat Transfer Conference. 1978. -P. 33-41.

27. Бабий В.И., Куваев Ю.Ф., Горение угольной пыли и расчет пылеуголь-ного факела. М.: Энергоатомиздат, 1986. - 208 с.

28. Горбис З.Р. Теплообмен и гидродинамика дисперсных сквозных потоков. М: Энергия, 1970. - 397 с.

29. Chukhanow Z.F. Heat and Mass transfer between gas and granular material. Part III Intern. Jorurnal Heat and Mass Transfer, 1971. - V.14. - P. 337 -351.

30. Заворин А. С. Проявление свойств минеральной части углей в паровых котлах. Монография. — Новосибирск: Изд-во НГТУ, 1998. — 221 с.

31. Заворин А. С. Состав и термические свойства минеральной части бурых углей (теплотехнический аспект). — Новосибирск: ИТ СО РАН, 1997. 187 с.

32. Эпик И.П. Влияние минеральной части сланцев на условия работы котлоагрегата. Таллинн: ЭГИ, 1961. - 249 с.

33. Вдовенко М.И. Минеральная часть энергетических углей. Алма-Ата: Наука, 1973. - 256 с.

34. Залкинд И.Я., Вдовченко B.C., Дик Э.П. Зола и шлаки в котельных топках. М.: Энергоатомиздат, 1988. - 80 с.

35. Юровский A 3. Минеральные компоненты твердых горючих ископаемых. М.: Недра, 1968.-215 с.

36. Красилъникова Л.Г. Исследование минерального состава канско-ачинских углей как энергетического топлива (применительно к углеразрезам урюпскому и ачинскому). -Дисс. канд. техн. наук. Томск: 1982. - 207 с.

37. Карякин С.К. Исследование минерального состава канско-ачинских углей в связи с их энергетическим использованием. — Дисс. канд. техн. наук. -Томск: 1975. 156 с.

38. Лебедев И.К., Карякин С.К., Заворин А. С. Результаты исследования исходной минеральной части Березовского угля. — Известия Томского политехнического института. Т. 283. — 1974. - С. 37-40.

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

40. Lockwood F.C., Salooja А.P., Syed A. A. A prediction method for coal-fired furnaces //Combustion and Flame. — 1980 . Vol.38, N.l. — P. 1 - 15.

41. Fischer K., Leithner R., Muller H. Three-dimensional of the gas-solid flow in coal-dust fired furnace //Proceeding of the First Int. Symposium on Two-Phase Flow Modelling and Experimentation. Rome. Italy. 1995. - V.l. - P. 1387 - 1393.

42. Дектярев А.А., Каменщиков Л.П., Ковалевский A.M. Программа AEROCHEM для моделирования трехмерных турбулентных реагирующих течений излучающего газа при наличии распыленных частиц //Вычислительные технологии. 1994. - Т.4, №12. - С. 107- 111.

43. Лебедев А.Н. Приготовление пылевидного топлива и сжигание его под котлами большой мощности. М: Госэнергоиздат, - 1948. - 178 с.

44. Повышение бесшлаковой мощности паровых котлов. М: Госэнергоиздат, 1958. - 224 с.

45. Устименко Б.П., Алияров Б.К., Абубакиров Е.К. Огневое моделирование пылеуголъных топок. — Алма-Ата: Наука, 1982. — 212 с.

46. Lockwood F.C., Shah N.G. A new radiation solution method for incorporation in general combustion prediction procedures // Proc. of the 18th Int. Symposium on Combustion. The Combustion Institute. Pittsburgh. - 1981. - P. 1405 - 1413.

47. Виленский Т.В., Хзмалян Д.М. Динамика горения пылевидного топлива.- М.: Энергия. 1978. - 248 с.

48. Основы практической теории горения / под. ред. В.В. Померанцева -Л.: Энергия, 1973. 264 с.

49. Горение натурального твердого топлива / под ред. А.Б. Резнякова -Алма-Ата: Наука, 1968. 410 с.

50. Горение твердого топлива //Материалы IVВсесоюзной конференции.- Новосибирск: Наука. 1974. - 205 с.

51. Блох А.Г. Теплообмен в топках паровых котлов. Л. : Энергоатомиз-дат, 1984. - 240 с.

52. Нигматулин Р.И. Динамика многофазных сред. 4.1. — М. : Наука, 1987.- 464 с.

53. Гришин A.M., Фомин В.М. Сопряженные и нестационарные задачи механики реагирующих сред. — Новосибирск: Наука, 1984. —319 с.

54. Алексеев Б.В., Гришин A.M. Введение в аэротермохимию. — Саратов: Изд-во СГУ, 1978. 419 с.

55. Зимонт В.Л., Трушин Ю.М. О суммарной кинетике горения углеводородных горючих // Физика горения и взрыва. 1969. - Т. 5, №4. - С. 567-573.

56. Ока С., Милошевич Д., Павлович П., Стефанович П. Математическое моделирование высокотемпературного турбулентного газового потока над поверхностью с интенсивной абляцией // Теплофизика и аэромеханика. — 1994. Т.1, №4. -С. 265-277.

57. Spalding D.B. Mathematical models of turbulent flames: a review // Combustion Science and Technology. — 1976. Vol. 13, N. 1-6. — P. 3-22.

58. Magnusen B.F. Hjertager B.H. On mathematical modeling of turbulent combustion with special emphasis on soot formation and combustion // Proceedings of 16th Int. Symposium on Combustion. Pittsburg. - 1976. - P. 719-727.

59. Волков Э.П., Зайчик Л.И., Першуков В.А. Моделирование горения твердого топлива. М.: Наука. 1994. - 320 с.

60. Оцисик М.Н. Сложный теплообмен. М.: Мир, 1976. - 616 с.

61. Pourahmadi F., Humpherey J.А.С. Modeling solid-fluid turbulent flows with application to predicting erosive wear // Physic-Chemical Hydrodynamics. — 1983. Vol. 4, N. 3. - P. 191-219.

62. Мостафа А.А., Монджиа Х.Ц., Макдонелл В.Г., Самуэлъсен Г.С. Распространение запыленных струйных течений. Теоретическое и экспериментальное исследование // Аэрокосмическая техника. 1990. - №3. - С. 65-82.

63. Асланян Г.С., Майков И.А. Численное исследование влияние турбулентности на процессы горения // Теплофизика высоких температур. — 1994. — Т.32, №6. -С. 892-901.

64. Кроу, Шарма, Сток. Численное исследование газокапельных потоков с помощью модели «капля-внутренний источник» // Теоретические основы инженерных расчетов. 1977. - Т. 99, №2. - С. 150-159.

65. Fischer К., Leithner R., Muller Н. Three-dimensional of the gas-solid flow in coal-dust fired furnace //Proceeding of the First Int. Symposium on Two-Phase Flow Modelling and Experimentation. Rome. Italy. 1995. - VI. - P. 1387 - 1393.

66. Тепловой расчет котельных агрегатов (нормативный метод) / Под ред. Н.В. Кузнецова. М.: Энергия, 1973. - 295 с.

67. Бубенчиков A.M., Старченко А.В. Пакет прикладных программ CHA1F для численного исследования аэродинамики и теплообмена внутренних течений. Томск: Изд. ТГУ. 1991. - 8с.

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

69. Finalayson В. A. The Method of Weighted Residuals and Variational Principles. -Academic, New York, 1972.

70. Harlow, F.N. and Welch, J.E. (1965). Numerical Calculation of Time-Depended Viscous Incompressible Flow of Fluid with Free Surface // Phys. Fluids. -vol. 8. P. 2182.

71. Caretto, L.S., Gosman, A.D., Patankar, S.V, and Spalding, D.B. (1972). Two Calculation Procedures for Steady, Three-Dimensional Flows with Recirculation. Proc. 3d Int. Conf. Num. Methods Fluid Dyn. - Paris. - Vol. II. - P. 60.

72. Patankar, S.V. and Spalding D.B. (1972a). A Calculation Procedure for Heat, Mass and Momentum Transfer in Three-Dimensional Parabolic Flows. Int. J. Heat Mass Transfer. - Vol. 15. - P 1787.

73. Красилъников С.В. Численное моделирование аэродинамики в топочных камерах парогенераторов // Энергетика: Экология, надежность, безопасность. Материалы докладов V Всероссийской научно-технической конференции. Томск: Изд-во ТПУ, 1999. - С. 104-105.

74. Noll В. Evaluation of a bounded high-resolution scheme for combustor flow computations //AIAA Journal. 1992. - V.30, №1. - P. 64-68.

75. Есаулов А.О., Старченко А.В. К выбору схемы для численного решения уравнений переноса //Вычислительная гидродинамика. Томск: Изд. Томск, унта. 1999. - С. 27-32.

76. Ильин В.П. Методы неполной факторизации для решения алгебраических систем. М.: Наука, 1995. - 288 с.

77. Самарский А.А., Николаев Е. С. Методы решения сеточных уравнений. М.: Наука, 1978. - 591 с.

78. Белов И.А., Коловандин Б.А., Кудрявцев Н.А. Развитие и перемещение вихревых структур у поверхности твердой стенки // Инженерно-физический журнал. 1989. Т. 56, №6. - С 900 - 909.

79. Fehlberg Е. Low order classical Runge-Kutta formulas with stepsize control. NASA TR R-315.

80. Shampine L.F., Watts H.A., Davenport S. Solving non-stiff ordinary differential equations The State of the Art, Sandia Laboratories Report SAND75-0182, Siam Review. - 18(1976), N3. - P. 376-411.

81. Красильников С.В. Численное моделирование аэродинамики, горения и степени шлакоулавливания в пылеугольных топках // Современные техника и технологии: Материалы V областной научно-практической конференции. — Томск: Изд-во ТПУ, 1999. С. 80-82.

82. Кеупо A.W., Krasinsky D.V., Salomatov V.V., Rychkov. Experimental and numerical modeling of the vortex furnace aerodynamics // Russ. J. Eng. Thermophys. 1996. - Vol. 6. -No. 1. -P. 47-62.

83. Маршак Ю.Л., Верзаков B.H. Исследование горения березовского угля в тангенциальной топочной камере с газовой сушкой топлива // Теплоэнергетика 1985. - №>1. - С. 4-9.

84. Томпсон Н. Секреты программирования трехмерной графики для Windows 95 /Перев. с англ. СПб.: Питер, 1997. - 352 с.