автореферат диссертации по энергетике, 05.14.04, диссертация на тему:Исследование теплообмена в плотном слое влажного угольного шлама и разработка методики теплового расчета бункера сырого угля

кандидата технических наук
Гордобаева, Татьяна Владимировна
город
Череповец
год
2005
специальность ВАК РФ
05.14.04
цена
450 рублей
Диссертация по энергетике на тему «Исследование теплообмена в плотном слое влажного угольного шлама и разработка методики теплового расчета бункера сырого угля»

Автореферат диссертации по теме "Исследование теплообмена в плотном слое влажного угольного шлама и разработка методики теплового расчета бункера сырого угля"

На правах рукописи

Гордобаева Татьяна Владимировна

Ч/

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

]

Специальность 05.14.04 - Промышленная теплоэнергетика

I

Автореферат

диссертации на соискание ученой степени кандидата технических наук

Череповец - 2005

Работа выполнена в Череповецком государственном университете

Научный руководитель: - доктор технических наук, профессор

Синицын Николай Николаевич

Официальные оппоненты: - доктор технических наук, профессор

Русаков Павел Григорьевич

- кандидат технических наук, доцент Запатрина Наталия Владимировна

Ведущая организация: ОАО « Череповецкая ГРЭС»

Защита состоится 27 января 2006 диссертационного совета Д 212.297.01 университете по адресу 162600, г. пр. Луначарского, д.5

года в 14.00 на заседании в Череповецком государственном Череповец Вологодской обл.,

С диссертацией можно ознакомиться в библиотеке Череповецкого государственного университета

Автореферат разослан 23 декабря 2005 г.

Ученый секретарь диссертационного совета

Никонова Е. Л.

20о£А

Ъ Ъ 2) ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность работы.

Запасы шлама коксохимического производства, которые скопились на ОАО «Северсталь», достигают громадных размеров. Этих запасов хватит более чем на 30 лет, если использовать их в качестве энергетического топлива ТЭЦ-ПВС и ТЭЦ-ЭВС ОАО «Северсталь».

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

Цель работы.

Исследование теплообмена влажного угольного шлама в бункере сырого угля и разработка методики теплового расчета бункера.

Задачи работы.

1) Разработать математическую модель теплообмена во влажном слое угольного шлама.

2) Разработать математическую модель прогрева отдельных частиц в слое шлама с учетом выхода влаги.

3) Экспериментально изучить процессы высыпания, обрушения и зависание шлама в бункерах.

4) Провести расчетное исследование процессов нагрева отдельных частиц в плотном слое угольного шлама.

5) Разработать инженерную методику теплового расчета поверхностей внутренних стенок бункера сырого угля.

6) Определить характеристики бункера, обеспечивающие движение материала без зависания.

Научная новизна.

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

2. На основе экспериментальных исследований установлены зависимости: *

• Времени высыпания угольного шлама из бункера от влажности, размера и формы выпускного отверстия и угла наклона боковой стенки бункера;

• Угла обрушения шлама от влажности и зольности.

РОС. НАЦИОНАЛЬНАЯ БИБЛИОТЕКА

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

Практическая значимость.

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

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

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

Апробация работы и публикации.

Результаты исследований и основные положения обсуждались на XXXIII конференции молодых специалистов (Череповецкий военный инженерный институт радиоэлектроники, 2005 г.), III Международной научно- технической конференции «Проблемы экологии на пути к устойчивому развитию регионов» ( Вологодский государственный технический университет, 2005), V Международной научно- технической конференции «Прогрессивные процессы и оборудование металлургического производства» (Череповец, 2005), Международной научно- технической конференции «Автоматизированная подготовка машиностроительного производства, технология и надежность машин, приборов и оборудования» (Вологда, 2005), IV Всероссийской научно-технической конференции «Современные промышленные технологии» (Нижний Новгород, 2005). По теме диссертации опубликовано 6 работ.

Структура и объем диссертации.

Диссертация состоит из введения, пяти глав, заключения, списка литературы из 110 наименований и приложений. Объем диссертации 144 страницы текста, включая 4 страницы приложений, 59 рисунков.

СОДЕРЖАНИЕ РАБОТЫ

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

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

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

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

Уравнение теплового баланса для частицы, находящейся в потоке горячих газов при отсутствии градиента температуры по сечению частицы (В1 <0,1), имеет вид:

с1Т ах

где т,- масса частицы, кг;

кс1г ^

сч- удельная теплоемкость материала частицы, кДж/(кг-К); Т., -температура частицы, К; т-время нахождения частицы в потоке горячих газов, с; а/. -коэффициент теплоотдачи , Вт/(м 2 К);

а =

!

N11

Тт - температура горячих газов, К; Р,юв - площадь поверхности частицы, м2 (/•"„„в = пкг), с/- диаметр частицы, м.

Дифференциальные уравнения прогрева частицы:

с1Т„ 6-Ш-1

£-(ТГ-Тч), при Тч <100°С, ах сц • рч •

где ¿/-диаметр частицы, м, рч- плотность материала частицы, кг/м3, Хг -коэффициент теплопроводности газов, Вт/(м К), Ый- критерий Нуссельта;

•га =аР -(Тг-Тч)^пов, при Тч =100°С,

йтч

■гп =а/, .(Гг-Гч).^ов, при Гч > 100°С,

где гп - теплота парообразования, Дж/кг.

Прогрев частицы при В1 > 0,1 описывается нестационарным уравнением теплопроводности с переменными коэффициентами температуропроводности, зависящими от температуры, и переменными граничными условиями.

В начальный момент времени (при т = 0) все точки частицы угля имеют одинаковую температуру Го. Частица имеет сферическую форму с радиусом г 0. После подачи горячего газа в слой угольного шлама в бункер частицы начинают прогреваться.

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

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

ох дг дг г дг

дТ2(г,х)_ 8 ЭГ2(г,х) 1к2 дТ2(г,х)

2Р2 Г"-~ -^ ) +--1- '

дх дг дг г дг

дТ^,х)_ р db дт2(г, X) nÇ,r) = TUL„ = const = 100° С;

начальные условия: T(r,Q)~T0-, f~0 (Ç-текущая координата фронта испарения).

dm dW dr dr '

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

дТ(г0,т)

- M-;- r=r„ = \Tmr-T ( Pq , т)] ,

дг

где Г (г, х) - текущая температура, К; г - текущая координата частицы; Т0 - начальная температура частицы, К; а Р - коэффициент теплоотдачи,

Вт/(м2 К); Г^ - температура потока, обтекающего частицу; р - плотность влажной частицы, кг/м3; с/=с/(Т/, IV р) и с2 = с2 (Т^ - теплоемкости влажного и сухого угля, кДж/(кг-К); Х1 = Я., (Ти Ш р) и Х2 = -

теплопроводности влажного и сухого угля, Вт/(м-К), р1= р] (Ти 1УР) и р2=р2(72)- -плотности влажного и сухого угля, кг/м3 ; гп -теплота парообразования, Дж/кг, IV р- влажность на рабочую массу, %, IV- текущее значение влажности, т - текущее значение массы.

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

Рис. 1. Схема бункера в системе координат (1 .боковая стенка бункера; 2: гипербола; 3: плоскость выпускного отверстия)

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

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

Уравнение данной линии было получено для конкретного бункера, имеющего такие размеры: диаметр выпускного отверстия 1м, длина боковой стенки 11,4 м, угол ее наклона 82

При допущении, что гипербола расположена так, как показано на рис. 1, совершив поворот системы координат, определили вид уравнения гиперболы.

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

Для этих условий решения Шумана представляются в виде:

1'г ' м / п

V г гм) О

где 9 м, 0 г - безразмерные температуры материала и газа по высоте слоя, t м - текущее значение температуры материала по высоте слоя, гг -

текущее значение температуры газа по высоте слоя, ^ и ^ - температуры

куИ

материала и газа в начальный момент времени, / =

сгм>г

- безразмерная

высота слоя, ку - объемный коэффициент теплопередачи, Вт/(м К), Сг -

_ ку х

удельная теплоемкость, кДж/(кг-К), IV г -скорость газа, м/с; £ =-

- безразмерное время, т - время (переменная величина), с, См -удельная теплоемкость материала частицы, кДж/(кг-К),/ -порозность слоя, /0 (х) -функция Бесселя первого рода от мнимого аргумента.

Для расчета ку необходимо определить поверхностный и объемный коэффициенты теплоотдачи.

Применяя к решениям Шумана метод интегрирования по частям и используя известное разложение функции Бесселя в ряд, получили:

6м —

(^г-С)

« -О

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

Ыи=0,61 -Яе0,67 (для Яе>200) N4=0,106Яе (для Яе<200). Здесь Яе = и>г • с//Уг; Ыи= а ^ • с1 Дг ,

где А.г - коэффициент теплопроводности газов, Вт/(м-К), Уг- коэффициент кинематической вязкости газа, м2/с, ¿/-диаметр частиц сферической формы, м. Свойства газа рекомендуется выбирать для средней температуры.

Отсюда можно получить поверхностный коэффициент теплоотдачи:

Для перехода к объемному значению коэффициента теплоотдачи использовали формулу

а, = а, -6-(1-/)/</,

М * 1ч

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

=

_

Вт

К

.а у л{\-г)хы)

где А =60 для частиц сферической формы.

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

у = ^ ' Н , м, сг -м>т

где Н - высота слоя, в котором осуществляется прогрев шлама, м. После этого нашли время прогрева шлама т по формуле:

= 2 См •(! - /) к,

Решение задачи дает возможность получить расчет температур материала и газа по высоте слоя.

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

Для испытания использовались установки типа воронки. Исследования проводились на двух образцах (обозначим их №1 и №2). Лабораторными методами были определены влажность (IV"), зольность (А*), выход летучих веществ (V), выход веществ на условную горючую массу.

Результаты эксперимента представлены на рис. 2- 5.

Рис. 2. Зависимость времени высыпания шлама от угла наклона боковой стенки бункера.

Рис. 4. Зависимость времени высыпания от влажности шлама.

Рис. 3. Зависимость времени высыпания шлама от диаметра выпускного отверстия.

Рис. 5. Зависимость угла обрушения от зольности шлама.

т,с

0,025 0,02 0,015 0,01 0,005

3 3,4 3,8 4,2 4,6 5 \л/ш,%

Зависимости, представленные на рис. 2 - рис. 5, записали в виде функций:

т=2,845 10"4- р -6,288-Ю'3 (р-угол наклона боковой стенки бункера, т=2,233-Ю32- ¿„л, 16>972 е'463>327'ают (</внп - диаметр выпускного отверстия, м); т=3,519-Ю"3- IV* +1,066-10"3 ( ^"-влажность шлама, %,т - время, с); аоб=-9,061 • А* +365,993 (А* - зольность шлама, %, Ооб - угол обрушения,

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

В 4 главе приведено расчетное исследование сыпучести шлама с учетом его прогрева в бункере.

Определили продолжительность прогрева шлама, характеризуемого насыпной массой Рте- 0,997 т/м3. Влажность шлама 9%. Порозность слоя определяли по формуле

где У0 -первоначальный объем шлама, м3, V- объем шлама, м3, т - его масса, кг, р - плотность слоя, т/м3.

Исходные значения шлама:

Ас = 33%, Ср= 53,6%, Яр = 3%, ИГ" = 7%, Ар = 30,7%. Сг= 86%, НТ = 4,8%, рорГ = 1,43 т/м3.

С учетом этих данных р ист = 1,75 т/м3, /= 43%. В исходном состоянии температура влажного шлама 0°С. Температурный интервал прогрева находится в пределах 0 - 120 0 С. Прогрев шлама осуществляется в слое высотой 80 мм при просасывании продуктов сгорания со скоростью 0,04 м/с с температурой газа на входе в слой 120° С. Эта температура соответствует верхнему пределу прогрева.

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

Затем получили зависимость температуры газа от времени прогрева на выходе из слоя (рис. 6).

Рис 6. Зависимость температуры газа от времени прогрева на выходе из слоя угольного шлама.

Аналогично установили зависимость температуры материала от этого же параметра. Результат представлен на рис. 7.

Кроме того, установлены зависимости температур газа и материала от времени прогрева в различных участках слоя высотой 0,08 м при скорости газа на входе в слой 0,04 м/с.

Для этого рассматривались участки слоя, имеющие высоту 0,01 м, 0,07 м, 0,08 м и на каждом из них проводились вычисления 6Г, 0„, /г и /м. Графики получившихся зависимостей представлены на рис. 8 и рис. 9.

Рис. 7. Зависимость температуры материала от времени прогрева.

Рис. 8. Температура газа в разных Рис. 9. Температура материала участках слоя высотой #=0,08 м. в разных участках слоя 0,08 м. ( 1:А=0,01 м; 2 :А=0,02 м; 3 : /И),07 м ).

Получены также зависимости критерия Ро от Ш при различных значениях высоты слоя (рис. 10). Для вычислений использованы формулы

Ро =

4 X г т

а

В1 = —— 21,

с„ рч</2 ¿Л.м

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

Рис. 10. Зависимости Ро от В1 для слоев высотой 0,01 м, 0,07 м, 0,08 м. ( 1:Я=0,01 м; 2 :Я=0,07 м; 3: #=0,08 м )

Функциональные зависимости Го от В1, полученные методом наименьших квадратов, представлены в таблице 1.

Таблица 1. Функциональные зависимости Ро=/(Вр.

Высота слоя, м Ро=/(В0 е, %

0,01 Ро=4,381-ВГ1-018-е45Л87В' 0,7

0,07 Ро=25,624-ВГ°л"Ч-'м"в' 9

0,08 Ро=30,813-ВГ0'966-е-35'778В' 7

Изучение полученных функциональных зависимостей позволило записать обобщенную формулу Ро=/ (В!) в виде Ро=/*В1яе'"в\ Для этого были найдены формулы связи А, В, С с высотой слоя Н. Обобщенная формула приобрела вид

Ро=( 369,349Я+0,571)В10'697И" |017-е (е=7„/о)

На рис. 11 показаны зависимости Ро от М), найденные для слоев высотой О м, 0,01 м, 0,07 м, 0,08 м при выходе 50%, 75%, 100% влаги из частицы.

5000 -4500 -4000 3500 -

Ро 3000 " 2500 -

2000 -

1500

1000 -

500

0 -

0

Рис.11. Зависимости Ро от ЯЬ для слоев высотой 0м, 0,01 м, 0,07 м, 0,08 м (I: Я=0 м; 2: Я=0,01 м; 3: Я=0,07 м; 4: Я=0,08 м).

Для нахождения значения ЛЬ применялась формула

яь = "»•Д,АГ, Ди - г,

где ДГ =100 °С; Л/и - убыль массы частицы после выхода влаги; га -теплота парообразования; т0 - начальная масса частицы; с- теплоемкость влажного шлама.

Функциональные зависимости Ро от Шэ представлены в таблице 2.

Таблица 2. Функциональные зависимости Ро=/(Щ>).

Высота слоя, м Ро=/(ЯЬ) е,%

0 Ро=-271,4176ЯЬ+728,116 5

0,01 Ро=-209,4483 -Шз+1689,143 0,8

0,07 Ро=-372,7957-ЯЬ+4770,144 0,5

0,08 Ро=-319,2977ЯЬ+4862,290 0,5

Полученные зависимости позволили записать обобщенную формулу Ро=/(Я, ЯЬ), где //-высота профеваемого слоя: Ро=(3-106Я5-370159-Я2+9595,1 Я-271,42)ЯЬ-620629Я-2+10Ш0Я+732,83. Аналогичным образом установлена и связь критериев Ро и Яе с учетом высоты профеваемого слоя:

Ро=(677795,2-Я+6637,594)-Яе"0'0258124|Н'0'9951233 (е=0,4 %).

На основании проведенных исследований при помощи метода наименьших квадратов установлен и характер зависимости относительного времени 2 от относительной высоты слоя К:

2= 1,2849- У+4,19 (е=0,33%). График этой зависимости представлен на рис. 12.

Рис. 12. График зависимости относительного времени 2 от относительной высоты слоя К

Пятая глава посвящена разработке опытно- промышленной установки бункера сырого угля.

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

у

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

Известно, что в гиперболических выпускных воронках достигается

постоянное значение коэффициента сужения ССУЖ = (Р — F')/F, где Г и

/<" - площади двух соседних поперечных сечений бункера, м2. Зная площадь нижнего основания бункера и используя то, что коэффициент сужения Ссуж не меняется, методом подбора установили вид гипербол, соответствующих тем или иным значениям коэффициента. Выяснили, что оптимальный вид гиперболы получается при значении коэффициента сужения, равного 0,15 (рис.13). При этом длина участка стенки, подвергающейся прогреву, равна 6 м, а максимальная высота прогрева шлама Н = 0,08 м.

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

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

Рис. 13. Изображение труб отсоса дымовых газов.

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

1) Определяем Аов и угол обрушения аОБ (рис. 14), зная влажность и зольность шлама, используя зависимости Аоь =/ (№*) и Аов =/ (А?)

Рис. 14. Изображение угла обрушения (ЛОБ-размер отверстия, образующегося при обрушении, а0Б - угол обрушения, И0Б - высота свода обрушения).

2) Определяем величину Иоь для данной модели бункера.

3) При известном значении Иов определяем длину £ части боковой стенки бункера, подвергающейся прогреву.

4) Находим высоту Я прогреваемого слоя шлама по формуле

2 н

где р-угол наклона боковой стенки бункера, ¿1 и ¿2 -длины сторон квадратных сечений ^ и проходящих через вершину гиперболы.

5) Определяем значение ЯЬ.

6) Зная Я и ЯЬ, находим величину Ро.

7) Определяем время т, необходимое для прогрева слоя высотой Я.

8) Используя зависимость Ро=/(Яе, Я), зная Я и Ро, определяем значение Яе

и находим скорость газа .

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

Основное содержание диссертации опубликовано в следующих работах:

1. Гордобаева Т.В., Синицын H.H. Нагрев неподвижного слоя кусковых материалов потоком газа с учетом изменения высоты слоя по линейному закону // Материалы 33 военно- научной конференции молодых специалистов института.- Череповец: ЧВИИРЭ, 2005,- С. 10

2. Гордобаева Т.В., Синицын H.H. Расчет температуры неподвижного слоя кусковых материалов потоком газа при изменении высоты слоя по гиперболическому закону // Материалы 33 военно- научной конференции молодых специалистов института,- Череповец: ЧВИИРЭ, 2005.- С. 11

3. Гордобаева Т.В., Синицын H.H. Расчет температуры неподвижного слоя кусковых материалов потоком газа при изменении высоты слоя по линейному закону // Материалы 33 военно- научной конференции молодых специалистов института.- Череповец: ЧВИИРЭ, 2005.-С.12

4. Гордобаева Т.В., Синицын H.H. Расчет подогрева слоя кусковых материалов потоком газа с учетом изменения высоты слоя по линейному закону // Проблемы экологии на пути к устойчивому развитию регионов: Материалы III Международной научно- технической конференции,-Вологда: ВоГТУ, 2005,- С. 169-171

5. Гордобаева Т.В., Синицын H.H. Расчет прогрева слоя кусковых материалов потоком газа с учетом изменения высоты слоя по гиперболическому закону // Проблемы экологии на пути к устойчивому развитию регионов: Материалы III Международной научно- технической конференции.- Вологда: ВоГТУ, 2005,- С. 172-175

6. Гордобаева Т.В. Исследование зависания угольных шламов в бункерах // Автоматизированная подготовка машиностроительного производства, технология и надежность машин, приборов и оборудования: Материалы международной научно-технической конференции,- Вологда: ВоГТУ, 2005.-С. 47-48

i

(

Щ--333

Подписано в печать 21.12.05 г. Формат 60x84/16. Гарнитура «Times». Ксерокопия. Уч.-изд. Л. 1,00. Тираж 100 экз. Заказ №

Отпечатано в РИО ГОУ ВПО ЧГУ 162600, г. Череповец, пр. Луначарского, 5

Оглавление автор диссертации — кандидата технических наук Гордобаева, Татьяна Владимировна

Введение.

Глава 1. Состояние вопроса.

1.1. Объект исследования.

1.2. Расчет бункеров.

1.3. Побудители сыпучих материалов в бункерах.

1.4. Питатели сыпучих материалов.

1.5. Основные положения механики насыпных грузов.

1.6. Расчет прогрева плотного слоя материала и газа.

1.7. Нагрев тел сферической формы с фазовым переходом.

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

1.9. Выводы по главе и постановка задачи исследования.

Глава 2. Математическое моделирование прогрева одиночной частицы угля из шламов в плотном слое.

2.1. Прогрев частицы угольного шлама при Bi < 0,1.

2.2. Прогрев частицы угольного шлама при Bi > 0,1.

2.3. Методика расчета температурного поля слоя угольных шламов.

2.4. Расчет гиперболической поверхности.

2.5. Выводы по главе.

Глава 3. Экспериментальное исследование сыпучести угольного шлама.

Ф 3.1. Методика проведения эксперимента и описание экспериментальной установки.

3.2. Анализ полученных экспериментальных данных.

3.3. Выводы по главе.

Глава 4. Расчетное исследование сыпучести шлама с учетом его прогрева в бункере.

4.1. Исследование температурного поля прогреваемого слоя угольного шлама.

4.2. Исследование температурного поля материала в прогреваемом слое угольного шлама.

4.3. Закономерности, полученные при исследовании прогрева слоя угольного шлама.

4.4. Выводы по главе.

Глава 5. Инженерная методика расчета бункеров для угольных шламов котла БКЗ 210-140 ТЭЦ-ПВС ОАО «Северсталь».

5.1. Расчет конструктивных параметров бункера.

5.2. Расчет параметров прогрева прогреваемого слоя материала.

5.3. Пример расчета бункера для котла БКЗ 210-140 ТЭЦ-ПВС ОАО

Северсталь».

5.4. Выводы по главе.

Введение 2005 год, диссертация по энергетике, Гордобаева, Татьяна Владимировна

Актуальность темы.

Запасы шлама коксохимического производства, которые скопились на ОАО «Северсталь», достигают громадных размеров. Этих запасов хватит более чем на 30 лет, если использовать их в качестве энергетического топлива ТЭЦ-ПВС и ТЭЦ-ЭВС ОАО «Северсталь».

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

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

Научная новизна.

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

2. На основе экспериментальных исследований установлены зависимости:

- времени высыпания угольного шлама из бункера от влажности, размера и формы выпускного отверстия и угла наклона боковой стенки бункера;

- угла обрушения шлама от влажности и зольности.

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

Практическая значимость.

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

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

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

Апробация работы и публикации. Результаты исследований и основные положения обсуждались на XXXIII конференции молодых специалистов (Череповецкий военный инженерный институт радиоэлектроники, 2005 г.); III Международной научно- технической конференции «Проблемы экологии на пути к устойчивому развитию регионов» ( Вологодский государственный технический университет, 2005); V Международной научно- технической конференции «Прогрессивные процессы и оборудование металлургического производства» (Череповец, 2005); Международной научно- технической конференции «Автоматизированная подготовка машиностроительного производства, технология и надежность машин, приборов и оборудования» (Вологда, 2005), IV Всероссийской научно-технической конференции «Современные промышленные технологии» (Нижний Новгород, 2005).

Структура и объем диссертации.

Диссертация состоит из введения, пяти глав, заключения, списка литературы из 110 наименований и приложений. Объем диссертации 144 страницы текста, включая 4 страницы приложений, 58 рисунков. По теме диссертации опубликовано 6 статей.

Заключение диссертация на тему "Исследование теплообмена в плотном слое влажного угольного шлама и разработка методики теплового расчета бункера сырого угля"

5.4. Выводы по главе.

Предложена инженерная методика расчета бункера сырого угля с учетом прогрева плотного слоя шлама потоком дымовых газов. V г

Представлены результаты расчета устройства согласно предложенной "3 | методике.

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

Библиография Гордобаева, Татьяна Владимировна, диссертация по теме Промышленная теплоэнергетика

1. А.С. № 1495241А1 СССР В 65 D 88/64. Способ выгрузки сыпучих материалов из бункера / В.У. Волощук, М.М. Вещезеров, В.П. Егоров.

2. А.С. СССР № 1512869, кл. В 65 D 88/64. Способ выгрузки трудносыпучих материалов из бункеров / В.И. Казаков, Ж.Э. Афлятунов,1988.

3. А.С. СССР № 1399220, кл. В 65 D 88/16. Бункер для слеживающихся сыпучих материалов / Г.М. Безкровный, С.В. Ковалевский, 1988.

4. А.С. СССР №979234, кл. В 65 D 88/66. Устройство для обрушения сводов в бункерах / В.Д. Пушкарных, Ю.В. Горошников, 1981.

5. А.С. СССР № 1472370, кл. В 65 D 88/64. Бункер для влажных сыпучих материалов / Ж.З. Афлятунов, 1987.

6. А.С. СССР № 1027099, кл. В 65 D 88/68. Устройство для сводообрушения трудносыпучего материала в бункере с выпускным отверстием / В.Н.Гутик, Л.И. Смелов,1983.

7. А.С. СССР № 1433869, кл. В 65 D 88/64. Бункерное устройство / Вл.В. Кривчиков, Вяч.В. Кривчиков,1968.

8. А.С. СССР № 491545, кл. В 65 D 88/64. Устройство для разгрузки силосов / Н.П. Тихомиров, В.В. Щеблыкин, В.Е. Трофимов, 1982.

9. А.С. СССР №979231, кл. В 65 D 88/64. Устройство для сводообрушения плохосыпучих материалов в бункерах, 1982.

10. А.С. СССР № 785119, кл. В 65 D 88/64. Устройство для обрушения сводов сыпучих материалов в бункерах, 1986.

11. Агроскин, А.А. Физика угля / А.А. Агроскин-М.: Недра, 1965.-365 с.

12. Алферов, К.В. Бункеры, затворы и питатели / К.В. Алферов M.-JL, Машгиз.-178 с.

13. Алферов, К.В. Бункерные установки / К.В. Алферов, P.J1. Зенков М.: Машгиз,1955.-308 с.

14. Андреев, Ю.В. Режимы горения, обеспечивающие снижение выбросов оксидов азота в дымовых газах паровых котлов, работающих на твердом топливе/Ю.В. Андреев- дисс.канд. техн. наук.- Череповец, ЧГУ, 2001.

15. Арутюнов, В.А. Металлургическая теплотехника. T.l / В.А. Арутюнов, В.И. Миткалинный, С.Б Старк-М.:Металлургия, 1974.-672 с.

16. Бабушкин, Н.М. Охлаждение агломерата и окатышей / Н.М. Бабушкин и др.-М.:Металлургия, 1975.-208 с.

17. Байков, В.Н. Проектирование железобетонных тонкостенных пространственных конструкций: Учеб. Пособие для вузов / В.Н. Байков, Э. Хампе, Э. Рауэ; Под. ред. В.Н. Байкова- М.: Стройиздат, 1990.-231 с.

18. Бакластов, A.M. Проектирование, монтаж и эксплуатация тепломассообменных установок: Учеб. пособие по специальности «Промышленная теплоэнергетика» / A.M. Бакластов-М.: Энергоиздат,1981.-336 с.

19. Бакластов, A.M. Промышленная теплоэнергетика и теплотехника / А. М. Бакластов, В. М. Бродянский, Б. П. Голубев и др.-М.: Энергоатомиздат, 1983.- 552 с.

20. Баскаков, А.П. Котлы и топки с кипящим слоем / А.П. Баскаков, В.В. Мацнев, И.В. Распопов- М.: Энергоатомиздат, 1996.-352 с.

21. Белеия, Е.И. Предварительно-напряженные несущие металлические конструкции / Е.И. Беленя- М.: Стройиздат, 1975.-415 с.

22. Вельский, М.Р. Усиление металлических конструкций под нагрузкой / М.Р. Вельский- Киев: Буд1вельник, 1975.-119 с.

23. Будак, Б. М. Численные методы в газовой динамике / Б. М. Будак, Ф. П. ;• • -- • S• rrj.

24. Васильев, А. Т. Егорова- М.: МГУ, 1965.- 158 с. |

25. Галин, Н.М. Тепло- массообмен (в ядерной энергетике) / Н.М. Галин, Л.П. Кириллов- М.: Энергоатомиздат, 1987.-376 с.

26. Глинков, М.А. Общая теория печей / М.А. Глинков, Г.М. Глинков-М.: Металлургия, 1978.-266 с.tZr

27. Гольдфарб, Э.М. Теплотехника металлургических процессов / Э.М. гф Гольдфарб- М.: Металлургия, 1967.-439 с. •т:.

28. Гордобаева, Т.В. Расчет температуры неподвижного слоя кусковых материалов потоком газа при изменении высоты слоя по гиперболическому £ закону/ Т.В. Гордобаева, Н.Н. Синицын // Материалы 33 военно- научнойг.;

29. Гордобаева, Н.Н. Синицын // Проблемы экологии на пути к устойчивому развитию регионов: Материалы III Международной научно- технической f1. Г" гконференции.- Вологда: ВоГТУ, 2005.-С. 172-175 vf

30. Гутер, Р.С. Элементы численного анализа и математической обработки результатов опыта/Р.С. Гутер, Б.В. Овчинский- М.: Наука, 1970.-432 с.

31. Днепров, Ю.В. Монтаж котельных установок малой и средней мощности / Ю.В. Днепров, Д.Н. Смирнов, М.С. Файнштейн- М.: Высшая школа, 1985.272 с.

32. Задачник по технической термодинамике и теории тепломассообмена: Учеб. пособие для вузов / В.Н. Афанасьев, С.И. Исаев, И.А. Кожинов; Под ред. В.И. Крутова , Г.Б. Петражицкого-М.: Высшая школа, 1986.-383 с.

33. Зенков, P.JI. Бункерные устройства / P.JI. Зенков, Г.П. Гриневич, B.C. Исаев- М.: Машиностроение, 1977.- 223 с.

34. Зенков, P.JI. Механика насыпных грузов / P.JI. Зенков-М.: Машиностроение, 1964.-251 с.

35. Зобнин Б.Ф. Теплотехнические расчеты металлургических печей. Учебное пособие для студентов вузов / Б.Ф. Зобнин, М.Д. Казяев, Б.И. Китаев, В.Г. Лисиенко, А.С. Телегин, Ю.Г. Ярошенко-М.:Металлургия, 1982.-360 с.

36. Исаков, В.Н. Элементы численных методов / В.Н. Исаков-М.: издательский центр «Академия», 2003.-192 с.

37. Исаченко, В.П. Теплопередача. Учебник для вузов / В.П. Исаченко, В.А. Осипова, А.С. Сукомел-М.: Энергия, 1975.-488 с.

38. Ицкович, A.M. Основы теплотехники / A.M. Ицкович-М.: Высшая Школа, 1975.-343 с.

39. Казанцев, Е.И. Промышленные печи / Е.И. Казанцев-М.: Металлургия, 1975.-366 с.

40. Калиткин, Н.Н. Численные методы: Учеб. пособие для вузов / Н.Н. Калиткин; Под ред. А.А. Самарского- М.: Наука, 1978.-512 с.

41. Карпин, Е.Б. Расчет и конструирование весоизмерительных механизмов и дозаторов / Е.Б. Карпин-М.: Машгиз, 1963.-524 с.

42. Карслоу, Г. Теплопросводность твердых тел / Г. Карслоу, Д. Егер- М.: Наука, 1964.-487 с.

43. Кацевич, J1.C. Теория теплопередачи и тепловые расчеты электрических печей / J1.C. Кацевнч-М.: Энергия, 1977.-304 с.

44. Квапил, Р. Движение сыпучих материалов в бункерах / Р. Квапил- М., 1961.

45. Кирпичев, М.В. Теплопередача / М.В. Кирпичев, М.А. Михеев, Л.С. Эйгенсон- M.-JL: Госэнергоиздат, 1940.-242 с.

46. Китаев, Б.И. Тепло- и массообмен в плотном слое / Б.И. Китаев, В.Н. Тимофеев, Б.А.Боковиков, В.М. Малкин, B.C. Швыдкий, Ф.Р. Шкляр, Ю.Г. Ярошенко-М.: Металлургия, 1972.-446 с.

47. Китаев, Б.И. Теплообмен в плотном слое / Б. И. Китаев, В. II. Тимофеев, Б. П. Боковиков и др.-М.: Металлургия, 1972.-431 с.

48. Китаев, Б.И. Теплотехнические расчеты металлургических печей / Б.И. Китаев, Б.Ф Зобнин, В.Ф. Ратников и др.-М.: Металлургия, 1970.-528 с.

49. Китаев, Б.И. Теплотехника доменного процесса / Б.И. Китаев, Ю.Г.

50. Ярошенко, E.J1. Суханов-М.: Металлургия, 1978.-242 с.

51. Китаев, Б.И. Теплообмен в доменной печи / Б.И. Китаев, Ю.Г. Ярошенко, Б.Л. Лазарев- М.: Металлургия, 1966.- 355 с.

52. Клипов, А.В. Исследование теплообмена в плотном и взвешенном слоях твердого топлива и разработка системы пылеприготовления, обеспечивающей снижение выбросов оксидов азота / А.В. Клинов-Дисс.канд.техн.наук.- Череповец, ЧГУ, 2004.-135 с.

53. Ключников, А.Д. Теплотехническая оптимизация топливных печей / А.Д. Ключников-М.: Энергия, 1974.-343 с.

54. Кноффе, Г.Ф. Теория топочных процессов / Г.Ф. Кноффе, К.М. Арефьев, А.Г. Блох; Под ред. Г.Ф. Кноффе- Л.: Энергия, 1966.-151 с.

55. Краснощеков, Е.А. Задачник по теплопередаче: Учеб. пособие для вузов / Е.А. Краснощеков, А.С. Сукомел- М.: Энергия, 1980.- 288 с.

56. Кривандин, В.А. Металлургические печи / В.А. Кривандин, Б.Л. Марков-М.: Металлургия, 1977.-463 с.

57. Кривандин, В.А. Металлургическая теплотехника.Т.2 / В.А. Кривандин, И.Н. Неведомская, В.В. Кобахидзе-М.: Металлургия, 1986.-592 с.

58. Кривандин, В.А. Тепловая работа и конструкция печей черной металлургии: Учебник для вузов / В.А. Кривандин, А.В. Егоров- М.: Металлургия, 1989.-462 с.

59. Кришер О. Научные основы техники сушки / О. Кришер-М.: Издатинлит, 1961 .-539с.

60. Кузнецов, А.Ф. Строительные конструкции из сталей повышенной и высокой прочности/А.Ф. Кузнецов- М.: Стройиздат, 1975.-81 с.

61. Кузовлев, В.А. Техническая термодинамика и основы теплопередачи / В.А. Кузовлев; Под ред. JI.P. Стоцкого- М.: Высшая школа, 1983.-335 с.

62. Кутателадзе, С.С. Основы теории теплообмена / С.С. Кутателадзе-Новосибирск: Наука, Сибирское отд., 1970.- 660 с.

63. Кушнырев, В.И. Техническая термодинамика и теплопередача: Учеб. Для вузов по спец. «Теплогазоснабжение и вентиляция» / В.И. Кушнырев, В.И. Лебедев, В.Л. Павленко- М.: Стройиздат, 1986.-461 с.

64. Лашутипа, Н.Г. Техническая термодинамика с основами теплопередачи и гидравлики / Н.Г. Лашутина, О.В. Макашова, P.M. Медведев- Л.: Машиностроение, 1988.-335 с.

65. Лыков, А.В. Тепломассообмен: Справочник / А.В. Лыков- М.: Энергия, 1078.-480 с.

66. Мак- Кракен, Д. Численные методы и программирование на фортране / Д. Мак- Кракен, У. Дорн- М.: Мир, 1977.- 584 с.

67. Малис, А.Я. Пневматический транспорт сыпучих материалов при высоких концентрациях / А.Я. Малис- М.: Машиностроение, 1969.-183 с.

68. Мастрюков, Б.С. Теория, конструкция и расчеты металлургических печей.Т.2 / Б.С. Мастрюков-М.: Металлургия, 1978.-270 с.

69. Мастрюков, B.C. Теплотехнические расчеты промышленных печей / B.C. Мастрюков-М.: Металлургия, 1972.-368 с.

70. Матвеев, Н.М. Обыкновенные дифференциальные уравнения: Учеб. Пособие для студентов пед. ин-тов по физ.-мат. спец. / Н.М. Матвеев- С.Петербург: специальная литература, 1996.-372 с.

71. Математическое и физическое моделирование и оптимизация Щи-'• *—

72. Межвузовский сборник научных трудов // Ивановский энегет. ин-т. им. В.И. ^ Ленина- Иваново: Ивановский гос. ун-т., 1981118 с.

73. Металлургические печи / Д.В. Бурдин, М.А. Глинков, М.В. Канторов; Под ред. М.А. Глинкова-М.: Металлургиздат, 1963.-687 с.

74. Михеев, М.А. Основы теплопередачи / М.А. Михеев, И.М. Михеева- М.: ТЖti t.1. Энергия, 1973.-319 с. ;|

75. Невский. А.С. Лучистый теплообмен в печах и топках / А.С. Невский- М.: | Металлургия, 1971.-433 с.

76. Новиков, А.Н. Методы борьбы со сводообразованием сыпучих материалов / А.Н. Новиков // Сборник НИИинформстройдоркоммунмаш..::1. М., 1966.-68 с.т

77. Оборудование для механизации перегрузки большегрузных контейнеровз

78. М.: НИИинформтяжмаш, 1973.-105 с. |v

79. Основы практической теории горения: Учеб. пособие для вузов / В.В. |

80. Померанцев, К.М. Арефьев, Д.Б. Ахмедов; Под ред. В.В. Померанцева-Энергоатомиздат, 1986.-312 с.

81. Павлов, И.И. Котельные установки и тепловые сети / И.И. Павлов, М.Н.

82. Федоров- М.: Стройиздат, 1986.- 232 с.

83. Петухов, Б.С. Теплообмен в ядерных энергетических установках / Б.С. Петухов, Л.Г. Гении, С.А. Ковалев- М.: Атомиздат, 1974. -408 с.

84. Практикум по теплопередаче: Учеб. пособие для вузов / А.П. Солодов, Ф.Ф. Цветков, А.В. Елисеев, В.А. Осипова; Под ред. А.П. Солодова-М.: Энергоатомиздат, 1986.- 296 с.

85. Промышленные тепломассообменные процессы и установки: Учеб. для вузов / A.M. Бакластов, В.Л. Горбенко, О.Л. Данилов- М.: Энергоатомиздат, 1986.-326 с.

86. Равич, М.Б. Газ и его применение в народном хозяйстве / М.Б. Равич- М.: Наука, 1974.-367 с.

87. Свенчанский, А.Д. Электрические промышленные печи. 4.1. / А.Д. Свенчанский-М.: Госэнергоиздат, 1958.-288 с.

88. Синицын, Н.Н. Исследование теплообмена и разработка методики теплового расчета плотного слоя скрапа / Н.Н. Синицын-Дисс.канд.техн.наук.- Череповец, ЧГУ, 2004.-133 с.

89. Синицын, Н.Н. Теплофизические процессы при движении одиночных частиц в газовом потоке/Н.Н. Синицын- Череповец, ЧГУ, 2001.-154 с.-~~-Ч-

90. Синицын, Н.Н. Теплофизические процессы при движении одно- и Ng многокомпонентных одиночных частиц из различных материалов в газовом Ш потоке и защита окружающей среды / Н.Н. Синицын- Дисс. уч. степ, д.т.н.- v Череповец, ЧГУ, 2001.- 477 с.

91. Спиваковский, А.О. Вибрационные конвейеры, питатели и вспомогательные устройства / А.О. Спиваковский, И.Ф. Гончаревич- М.:1. Щ.

92. Машиностроение,1972.-327 с.•1г

93. Тарасов, Ю.Д. Очистка подконвейерного пространства на предприятиях ;f-rr.нерудной промышленности / Ю.Д, Тарасов- JL: Стройиздат, 1983.-193 с. 1• —:т.

94. Телегин, А.С. Термодинамика, тепломассопередача / А.С. Телегин, B.C. Швыдкий, Ю.Г. Ярошенко-М.: Металлургия, 1980.-264 с.

95. Теория тепломассообмена: Учебник для вузов / С. И. Исаев, И. А.- ~

96. Кожинов, В. И. Кофанов и др.-М.: Высшая школа, 1979.- 495 с. -3

97. Тепловые процессы в технологических системах / Сборник научных zS. трудов- Череповец: ЧГИИ, 1996.-76 с. |1. ТГ1

98. Теплотехника. Учебник для студентов втузов / A.M. Архаров, С.И. Исаев, Lri И.А. Кожинов; Под общ. ред. В.И. Крутова- М.: Машиностроение, 1986.432 с.

99. Теплотехника и теплоэнергетика металлургического производства: ;•:1. У

100. Сборник задач: Учеб. Пособие для вузов / В.В. Белоусов, А.Г. Клевцов, И.А. 2

101. Прибытков, Г.С. Сборщиков- М.: Металлургия, 1993.-335 с.

102. Теплотехнический справочник / Под общ. ред. В.Н. Юренева и П.Д.1Й: а- r-—zz

103. Лебедева // В 2-х т. T.l.-М.: Энергия, 1975.-744 с.

104. Тихонов, А.Н. Уравнения математической физики: учебное пособие / ^1. Ч

105. A.Н. Тихонов, А.А. Самарский-М.: Наука, 1972.-736 с. *-■ТЛ.гг£

106. Успехи теплопередачи / Сборник статей // Пер. с англ.- М.: Мир,1971.- f 576 с.

107. Филимонов, Ю.П. Топливо и печи: Учебник / Ю.П. Филимонов, Н.С. Громова- М.: Металлургия, 1987.-320 с.- »:-г

108. Хзмалян, Д.М. Теория горения и топочные устройства: Учеб. пособие /1. Z-rr

109. Д.М. Хзмалян, Я.А. Каган; под ред. Д.М. Хзмаляна-М.: Энергия, 1976.-487 с. S• *тт

110. Холоднов, В.А. Математическое моделирование и оптимизация химико- ' технологических процессов. Практическое руководство / В.А. Холоднов, 1;

111. B.П. Дьяконов, Е.Н. Иванова, Л.С. Кирьянова СПб.: АНО НПО «Профессионал», 2003.-480 с.

112. Чернобыльский, И.И. Сушильные установки химическойл•-'спромышленности / И.И. Чернобыльский, Ю.М. Танайко- Киев: Технжа, г?1969.-280 с. Ч• -t

113. Шираке, З.Э. Теплоснабжение: Пер. с латыш. / З.Э. Шираке- М.:Энергия,т: ■г;1979.-255 С.

114. Щукин, А.А. Промышленные печи и газовое хозяйство заводов / А.А. Щукин- М.: Энергия, 1973.-223 с.

115. Эккерт, Э.Р. Теория тепло- и массообмена / Э.Р. Эккерт, P.M. Дрейк- М.:г -iir

116. Госэнергоиздат, 1961.-680 с.•3139 :1•• *f4