Journal of Siberian Federal University. Engineering & Technologies, 2016, 9(4), 603-610
УДК 697.9
To Calculation Heat- and Water Emission at Design of Ventilation in Conveyor Galleries for Transportation of Wet Heated Materials
Ilya V. Kalinich and Gennady V. Smolnikov
Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041, Russia
Received 27.01.2016, received in revised form 13.03.2016, accepted 21.05.2016
Article is devoted to questions of a heat and mass transfer of wet heated materials at conveyor transportation. Charge material and sinter return of are investigated. The data necessary at design of systems of ventilation in conveyor galleries for transportation of these materials are obtained.
Keywords: emission of water vapour, heat transfer, coefficient of heat transfer, coefficient of heat conductivity, mass transfer coefficient.
Citation: Kalinich I.V., Smolnikov G.V. To calculation heat- and water emission at design of ventilation in conveyor galleries for transportation of wet heated materials, J. Sib. Fed. Univ. Eng. technol., 2016, 9(4), 603-610. DOI: 10.17516/1999-494X-2016-9-4-603-610.
К расчёту тепло- и влаговыделений при проектировании вентиляции в конвейерных галереях для транспортировки влажных нагретых материалов
И.В. Калинин, Г. В. Смольников
Сибирский федеральный университет Россия, 660041, Красноярск, пр. Свободный, 79
Статья посвящена вопросам тепломассообмена влажных нагретых материалов при конвейерной транспортировке. Исследованы шихта и возврат агломерата. Получены данные, необходимые при проектировании систем вентиляции в конвейерных галереях для транспортировки данных материалов.
© Siberian Federal University. All rights reserved
Corresponding author E-mail address: il-k@mail.ru, genad123@mail.ru
*
Ключевые слова: выделение водяных паров, теплоотдача, коэффициент теплоотдачи, коэффициент теплопроводности, коэффициент массообмена.
На горнообогатительных комбинатах широко распространена конвейерная транспортировка влажных нагретых материалов. При проектировании систем вентиляции конвейерных галерей для транспортировки влажных нагретых материалов возникает необходимость определения количества теплоты и влаги, выделяющихся от поверхности материала в объём помещения. Для этой цели необходимы данные об изменении температуры поверхности материала при его транспортировании.
Суммарные поступления явного тепла от единицы поверхности материала в объем помещения определяются по формуле
q* = q + qT > (!)
где q - общая теплоотдача материала конвекцией и излучением, Вт; дт - количество явного тепла, поступающего в воздух за счёт остывания водяных паров, Вт.
Общая теплоотдача материала конвекцией и излучением рассчитывается как
q = а • (tM -tB), (2)
где а - суммарный коэффициент теплообмена, Вт/(м2*0С); /м - температура материала, 0С; 4 - температура окружающего воздуха, 0С.
Количество явного тепла, поступающего в воздух за счёт остывания водяных паров,
qT = с • m • (tM - tB), (3)
где c - теплоёмкость водяного пара, Дж/(кг*0С); m - интенсивность выделения водяных паров, (кг*м2)/с.
Интенсивность выделения водяных паров при внешнем массообмене подчиняется закону Дальтона:
m = Р • Р -ФВ • Рв )' (4)
где в - коэффициент массообмена, зависящий от крупности и влажности материала, а также от скорости движения конвейерной ленты [2], c/м; рм - парциальное давление насыщенных водяных паров при температуре поверхности материала, Па; рв - парциальное давление насыщенных водяных паров при температуре воздуха в помещении, Па; фв - относительная влажность воздуха в помещении, %.
По характеру тепло- и массообмена участки транспортировки влажных нагретых материалов в конвейерных галереях на горнообогатительных комбинатах можно разделить на две группы.
К первой группе относятся конвейеры с двухслойной укладкой материала. Нижним слоем является влажная шихта, имеющая температуру около 20 0С и влажность до 16 %. Верхний слой состоит из сухого возврата агломерата с температурой 400 - 800 0С.
Ко второй группе относятся конвейеры с равномерно прогретой до температуры 60 - 80 0С шихтой, влажность которой составляет 6 - 10 %.
Исследованиями авторов в производственных условиях установлено, что при транспортировании материалов второй группы изменение температуры и влажности материала по длине конвейера в пределах одной галереи незначительно. В связи с этим процесс поступления тепла и водяных паров от поверхности материала в объём галереи можно рассматривать как стационарный процесс при постоянной температуре и влажности материала.
При транспортировке материалов первой группы теплоотдача от поверхности в окружающую среду и в слой шихты имеет существенно нестационарный характер, что связано с остыванием поверхности верхнего слоя частиц.
Оценка интенсивности охлаждения возврата агломерата была произведена авторами с использованием закономерностей нестационарной теплопроводности, руководствуясь данными работы [3]. В расчётах приняты следующие допущения:
1) материал состоит из шаров различного диаметра соответственно гранулометрическому составу;
2) охлаждение шаров водяными парами является несущественным.
Процесс нестационарной теплопроводности определяется критериями Био и Фурье:
ш= ^; (5)
Ро=^м_1, (6)
l2
где а - коэффициент теплоотдачи от поверхности шара в окружающую среду, Вт/(м2*0С); l - определяющий размер, за который в рассматриваемом случае принят радиус шара [3], м; Хм - коэффициент теплопроводности материала, принятый по данным работы [4] для возврата агломерата, равный 0,65 Вт/(м2*0С); aм - коэффициент температуропроводности возврата агломерата, равный, по данным работы [3] 93*10-9, м2/с; т - время транспортировки, с.
При малых значениях числа Био (Ы < 1) решение дифференциального уравнения тепло -проводности для охлаждения шара приводится к виду
0Т =00 • ехр • (-3 • В{ ■ Бо), (7)
где 9т - избыточная температура поверхности шара в любой момент времени, 0С; 90 - избыточная температура поверхности шара в начальный времени, 0С.
Избыточная температура поверхности шара в любой момент времени определяется по формуле
0х = К - К, (8)
где ^ - температура поверхности шара в любой момент времени, 0С.
Избыточную температуру поверхности шара в начальный момент времени находим по формуле
0 0 = 'о - 'в, (9)
где 4 - температура поверхности шара в начальный момент времени, 0С.
- 605 -
В соответствии с уравнениями (5) - (7) интенсивность остывания шаров зависит от их диаметра, поэтому следует рассмотреть остывание отдельных шаров различного диаметра. Расчёт производили в следующем порядке. Определяли коэффициент теплоотдачи а как
а=ал +ак, (10)
где ал - коэффициент теплоотдачи излучением, Вт/(м2*0С); ак - коэффициент теплоотдачи конвекцией, Вт/(м2*0С).
Коэффициент теплоотдачи излучением, согласно данным работы [5], равен:
а„ =■
(tx + 273 ^4 Ч + 273 V1
1 --1 +1 --1
^ 100 ) .100 J J
(11)
t -1
м в
где епр - приведённая степень черноты обменивающихся теплом поверхностей, равная для материалов ограждающих конструкций и возврата агломерата 0,95; с0 - коэффициент излучения абсолютно чёрного тела, равный 5,67 Вт/(м2* 0С4).
При определении коэффициента конвективного теплообмена было принято допущение, что при высоких начальных температурах теплоотдающей поверхности в условиях конвейерной транспортировки преобладает теплоотдача за счёт естественной конвекции, поэтому величины ак находили, согласно данным работы [6], из критериального уравнения
№=0,54 • (вг • Рг)0Д5 при 500 < Ог • Рг<2 -107, (12)
где Ми - число Нуссельта; вг - число Грасгофа; Рг - число Прандтля; Число Нуссельта устанавливали как
а • I
^ = —, (13)
где I - определяющий размер, принятый, согласно данным работы [6]: ¡=ё, м; X - коэффициент теплопроводности воздуха, равный 0,026 Вт/(м*0С). Число Грасгофа определялось как
Gr = JXCP2 " , (14)
• I Ч^, - ^)
Т ■ V2
где g - ускорение свободного падения, g = 9,81м/с2; 4ср - средняя температура материала на участке, 0С; Тв - температура окружающего воздуха, К; V - кинематическая вязкость воздуха, V = 15*10-6 м2/с.
Число Био рассчитывали при а = ал + ак для каждого диаметра.
Расчёт выполнялся по участкам, разбитым на длины, соответствующие промежуткам времени транспортировки т, равным 5, 10 и 20 с. За начальную температуру материала на последнем участке принималась температура в конце предыдущего участка.
Расчёты проводили для возврата агломерата с начальной температурой t0 = 800 0С и t0 = 400 0С при температуре окружающего воздуха tB = 15 0С и tB = 20 0С.
Как и следовало ожидать, частицы меньшего диаметра остывают более интенсивно. Для них влияние начальной температуры на конечную незначительно при одинаковом времени транспортировки. Температура окружающего воздуха в пределах 4 = 15 - 20 0С не оказывает влияния на конечную температуру поверхности частиц.
Для дальнейших расчётов необходимо произвести усреднение температуры по поверхности слоя возврата агломерата. Допустим, что распределение частиц по поверхности слоя возврата такое же, как и в массиве, и для усреднения температуры используем формулу
X^ • F =WP -XF, (15)
i=1 i=1
где Fi - площадь поверхности частиц размером di , м2; - температура поверхности частиц размером di , 0С; 4>ср - средняя по поверхности температура частиц, 0С.
Используя соотношения
F~ di2 • n (16)
M N. M
Р' d
получим
N ■ М, (18)
1 Л р
где щ - количество частиц размером di ; Mi - масса частиц размером di , кг; M - масса всего возврата агломерата, кг; N - массовая доля частиц размером di ; р - плотность возврата агломерата, кг/м3.
Подставим (18) в (15):
(X • Л
'«.ср = ^=4- • (19)
(X N1)/Л,
1=1
Усреднение температуры проводили для следующего гранулометрического состава возврата [2], приведённого в табл. 1
Изменение средней температуры поверхности возврата агломерата по длине конвейера приведено на рис. 1. Наибольшее различие в темпах остывания характерно лишь для начальных участков транспортировки, поэтому на этих участках будут более интенсивные тепловыделения.
Полученная зависимость температуры поверхности возврата агломерата от времени транспортировки позволяет определить по формуле (2) количество тепла, поступающего в помеще-
и
n
Таблица 1. Гранулометрический состав возврата агломерата
Класс крупности частиц, мм 20-10 10-5 5-2,5 2,5-1,25 1,25-0,63 0,63-0
Средняя крупность частиц ( , мм 15 7,5 3,76 1,83 0,94 0,32
Массовое содержание частиц N , % 2 31 25 24 8 10
Рис. 1. Изменение температуры поверхности возврата агломерата при конвейерной транспортировке: 1 - при начальной температуре возврата агломерата 800 0С; 2 - при начальной температуре возврата агломерата 400 0С
ние. При этом тракт транспортировки материала делится на участки, на которых изменение температуры материала сравнительно невелико, определяется средняя температура материала на участке и количество тепла, поступающего в помещение на данном участке.
Согласно зависимости (4) для определения интенсивности влаговыделений от конвейера необходимо знать парциальное давление насыщенных водяных паров - рм при температуре на границе слоёв шихты и возврата агломерата - /гр .
Оценим температуру на границе слоёв исходя из условия равенства подводимого и отводимого тепла:
? - ? ? - ?
м,ср гр гр ш
—• <20)
где Яш - термическое сопротивление слоя шихты, Вт/м2; Яа - термическое сопротивление слоя возврата агломерата, Вт/м2; /ш - температура шихты на ленте конвейера, 0С. Термическое сопротивление слоя шихты определяется по формуле
^ = ^, (21)
где 5ш - толщина слоя шихты, м; 5ш = 0,16 м; - коэффициент теплопроводности шихты, Вт/м*0С; = 1,54 Вт/м*0С.
Рис. 2. Изменение температуры поверхности слоя шихты при конвейерной транспортировке: 1 - при начальной температуре верхнего слоя 800 0С; 2 - при начальной температуре верхнего слоя 400 0С
Термическое сопротивление слоя возврата агломерата можно вычислить по формуле
R = ^, (22)
где 5а - толщина слоя возврата агломерата, м; 5а = 0,04 м; Хм - коэффициент теплопроводности возврата агломерата, Вт/м*0С; Хм = 0,65 Вт/м* 0С.
Исходные данные для 5ш ; 5а ; ; Хм приняты согласно работе [4].
Решение уравнения (20) с использованием данных о распределении температуры возврата агломерата по длине конвейера (рис. 1) позволило получить изменение температуры поверхности слоя шихты по длине конвейера (рис. 2).
По температуре поверхности шихты /гр можно определить давление рм и количество водяных паров, поступающих в помещение. При этом на участках конвейера, где температура шихты выше температуры кипения при данном барометрическом давлении, величину парциального давления насыщенных водяных паров рм в формуле (4) следует принимать равной величине барометрического давления ра . На участках, где температура шихты ниже температуры кипения, парциальное давление насыщенных водяных паров рм определяется, как обычно, в зависимости от температуры.
Поученные данные об изменении температур поверхности возврата агломерата и шихты по длине конвейера могут быть использованы для расчёта тепло- и влаговыделений от поверхности транспортируемого материала при проектировании вентиляции конвейерных галерей для транспортировки влажных нагретых материалов.
Список литературы
[1] Шумилов Р. Н. Расчёт влаговыделений при транспортировке зернистых материалов конвейерами. Водоснабжение и санитарная техника. 1968, 3(8), 30-31 [Shumilov R. N. Calculation
- 609 -
of Moisture Allocations When Transporting Granular Materials by Conveyors, Water Supply and Sanitary Equipment, 1968, 3(8), 30-31 (in Russian)]
[2] Афанасьев И. И., Ващенко В. С., Генералов Г. С. Обеспыливание воздуха на фабриках горнообогатительных комбинатов. М.: Недра, 1972. 184 с. [Afanas'ev I. I., Vashhenko V. S., Generalov G. S. Obespylivanie vozduha na fabrikah gornoobogatitel'nyh kombinatov. Moscow, Nedra, 1972, 184 p. (in Russian)]
[3] Шорин С. Н. Теплопередача. М.: Высшая школа, 1964. 490 c. [Shorin S. N. Teploperedacha. Moscow, Vysshaja shkola, 1954, 490 p. (in Russian)]
[4] Марииченко В. М. Электромоделирование нестационарного теплообмена между шихтой и возвратом агломерата при конвейерном транспортировании. В кн.: Охрана труда и техника безопасности в горнорудной промышленности. М.: Недра, 1978. 69 c. [Mariichenko V. M. Jelektromodelirovanie nestacionarnogo teploobmena mezhdu shihtoj i vozvratom aglomerata pri konvejernom transportirovanii. V kn.: Ohrana truda i tehnika bezopasnosti v gornorudnoj promyshlennosti. Moscow, Nedra, 1978, 69 p. (in Russian)]
[5] Исаченко В. П., Осипова В. А., Сукомел А. С. Теплопередача. М.: Энергия, 1981. 416 c. [Isachenko V. P., Osipova V. A., Sukomel A. S. Teploperedacha. Moscow, Jenergija, 1981, 416 p. (in Russian)]
[6] Лыков А. В. Теоретические основы строительной теплофизики. Минск: Издательство АН БССР, 1961. 520 с. [Lykov A. V. Teoreticheskie osnovy stroitel'noj teplofiziki. Minsk: Izdatel'stvo AN BSSR, 1961, 520 p. (in Russian)]