ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ
Кузнецов В. А., д-р техн. наук, проф., Трулев А. В., аспирант
Белгородский государственный технологический университет им. В. Г. Шухова
НЕСТАЦИОНАРНАЯ ТЕМПЕРАТУРА СТЕН ЦЕМЕНТНОЙ ВРАЩАЮЩЕЙСЯ
ПЕЧИ
kouzn@intbel.ru
Разработан аналитический подход к расчету нестационарной температуры стен вращающейся печи при численном моделировании в ней процессов горения и теплообмена. Адекватность подхода проверена численным решением дифференциального уравнения нестационарной теплопроводности стен в условиях вычислительного эксперимента.
Ключевые слова: цементная вращающаяся печь, нестационарная температура, математическая модель, аналитический подход, вычислительный эксперимент.
Постановка задачи.
При обжиге цементного клинкера внутренняя поверхность стен вращающийся печи подвергается тепловому воздействию на нее поочередно газообразной среды и слоя технологического материала, что вызывает цикличное изменение ее температуры, способное приводить в итоге к разрушению футеровки. Определение нестационарной температуры вращающихся стен является актуальной задачей, решение которой необходимо для создания условий, обеспечивающих более длительную безостановочную работу печей.
Цель достигается путем численного решения системы дифференциальных уравнений процессов движения газов, горения и теплопереноса в печи. При этом нестационарная температура Тст внутренней поверхности стен печи является одной из искомых функций.
Чтобы найти эту температуру в процессе выполнения итераций, составляется уравнение баланса тепла, подводимого к внутренней поверхности стен печи от газов и слоя клинкера и отводимого от нее через стены печи:
Чк + ^ + < + Чнст = Чотв, (!)
где чк , длг , длм - плотность тепловых потоков подвода теплоты соответственно конвекцией, излучением газов и излучением слоя материала, днст - плотность нестационарного потока теплоты, приносимой в участок I на поверхности стены из соседнего /+1 при вращении печи.
Как показано в [1], расчетные выражения для плотности тепловых потоков содержат искомую термодинамическую температуру Тст внутренней поверхности стен, так что ее можно выразить в явном виде и вычислить:
Чк =л/рстр| cpС (Тг " Тст VAw '
яЛ = 2ест° (т4 - «пТСТ/а)/(3аЛи),
= (1 -
а,
l/a) 8пра(у/,кл ткл
( т4 7^4 ) Wl кл ст /'
Здесь рст, рг - плотность газов на поверхности стен и вблизи нее, ср - удельная массовая теплоемкость газов при постоянном давлении, ^ -коэффициент в граничных условиях турбулентного движения газов, с - постоянная Стефана-Больцмана, вст , впр - степень черноты внутренней поверхности стен печи и приведенная, ап, а -планковский и локальный коэффициенты поглощения, Тг , Ткл - температура газов вблизи стенки и слоя обжигаемого клинкера, Т - определяющая температура радиационного переноса энергии, Ап - расстояние от ближайшего узла сетки в газовой среде до стенки печи, у1кл - угловой коэффициент излучения участка I стены на слой клинкера.
Плотность нестационарного потока теплоты днст в поверхностном слое Аг стены определяется с помощью дифференциального уравнения нестационарной теплопроводности:
dT д ( . dT
рс— =-1 rX—
дг rdr I дг
(2)
где р, с, X - плотность, удельная массовая теплоемкость и коэффициент теплопроводности материала поверхности стены, т - время, г - координата по радиусу печи, Т - температура.
Для численного решения задачи стены печи по толщине разбиваются сеткой на слои, причем толщина поверхностного слоя Аг не превышает 1 мм. Интегрирование дифференциального уравнения (2) по поверхностному слою Аг дает следующее дискретное выражение:
ЛТ(т,1 - TW+1 )Лг .Л (TU - TW ^
где Лг1 - интервал времени перемещения элемента внутренней поверхности футеровки из узла сетки l+1 в узел l при вращении печи, To l - температура в узле сетки j=0, лежащем на внутренней по-
dT
верхности участка / стены, Тц - температура на том же участке / в узле }=1, ближайшем по оси г к поверхности стены,.
Последнее слагаемое дискретного уравнения заменяется суммой тепловых потоков в соответствие с граничными условиями (1):
Чк + ^л + ?лм + (т0, /+1 - Т0, /) Лг = "Г (Т0,/ - ТИ ) •
Лт/ Лг
Из сопоставления уравнения (1) с полученным следуют выражения для плотности нестационарного и отводимого тепловых потоков:
Чнст
£С_
Лт/
(То, /+1 - То/ )Лг, (3)
Чотв - ~Г(т0,/ - Т1,/)•
Аг
(4)
При компьютерном моделировании приходится решать сопряженную задачу радиационно-конвективного теплообмена в пламенном пространстве вращающейся печи и нестационарной теплопередачи через ее стены при неизвестной температуре внутренней поверхности стен. При этом разность двух неизвестных температур в формуле (4) вносит в расчет итерационную погрешность, которая многократно возрастает при ее делении на малую величину Лг. В таких условиях трудно ожидать сходимости итерационного решения поставленной задачи. Поэтому для учета нестационарного изменения температуры стен при численном моделировании теплообмена в цементной вращающейся печи применен упрощенный аналитический подход.
Аналитическое выражение нестационарного теплового потока.
Пусть в момент входа вращающейся стены под слой клинкера температура ее внутренней поверхности скачкообразно изменяется от Т0 до более высокой температуры слоя клинкера Ткл. Применим к этому случаю аналитическое выражение плотности нестационарного потока аккумулируемой теплоты на поверхности полуограниченного твердого массива при внезапном изменении температуры его поверхности [2]:
Чак
Хрс
пт
(Тн - То ),
(5)
где т - время, отсчитываемое от момента температурного скачка, То , Тн - соответственно первоначальная и новая температура поверхности нагреваемого полуограниченного массива.
Теплоотдача от слоя клинкера к футеровке отличается от задачи [2] тем, что в печи непрерывно входят под слой клинкера и испытывают температурный скачок все новые участки стен, время т нагрева которых клинкером при прохождении угла Лф возрастает обратно пропорцио-
нально угловой скорости их вращения ю в соответствие с соотношением т = Лф/ю (рис. 1).
Поскольку температура стен зависит от времени их пребывания под слоем клинкера, то количество теплоты Qaк, аккумулированное стеной за один оборот на единице длины печи, находится интегрированием плотности теплового потока дак , определяемого формулой (5), по поверхности теплообмена Е, пропорциональной углу ф, с учетом зависимости времени т от угла ф:
<2ак - I Чак^ - Л / — юЯ('кл - 'о ) I Ф ^Ф
ЛФк
-12.
^ак I Чак'о >
0 * Ж 0
где Я - внутренний радиус печи, Лфкл - центральный угол, занимаемый слоем клинкера в печи. Физические свойства материала стен полагаются постоянными.
Рис. 1. Схема нестационарного теплообмена в печи: 4л - температура клинкера, , '0 - температуры поверхности стенки печи, Лфкл - центральный угол клинкерного слоя, / - номер расчетного участка
После математических преобразований получаем расчетную формулу:
0ак - 2 ^Дфи-го)-V п ,
- 2Ьл Дфклкл - г0)
(6)
где Лткл - максимальное время пребывания футеровки под слоем клинкера, 6кл - коэффициент, аналогичный коэффициенту теплоусвоения [2],
ькл - 2
Хрс
пДт„
Примем, что стены печи в течение одного оборота отдают теплоотдачей столько же теплоты, сколько они получили от клинкера. В этом случае можно приравнять друг другу правые части формул (3) и (6), учитывая при этом различие центральных углов стены, контактирующей с клинкером (Лфкл) и газовой средой (Лфгаз):
(г кл - г 0 )ЛГэк Лфгаз -
Лт
- 2
Яре
жЛт„
(гкл - г0 )Лфк
Отсюда следует оценка эквивалентной толщины Дгэк поверхностного слоя футеровки, пригодная для расчета нестационарности теплообмена:
= 2,
ДФкл
пРсДт кл Дфгаз
Дт Газ = 2
ХДт,,
.(7)
прс
Примем приближенно X =1 Вт/(мК), р = 2000 кг/м3, с = 1500 Дж/(кг К), Дт = 15 с. По формуле (7) получим эквивалентную толщину поверхностного слоя Дгэк « 2,5 мм, составляющую примерно сотую долю от толщины стен печи.
Поскольку толщина слоя с нестационарной температурой сравнительно мала, можно полагать, что теплота, аккумулируемая им в пределах одного оборота печи, слабо влияет на распределение температуры по толщине стен печи в целом. Подставив выражение (7) для эквивалентной толщины слоя Дгэк в равенство (3), получим приближенную расчетную фор формулу для плотности дополнительного теплового потока через поверхность расчетного участка стены, возникающего из-за нестационарности процесса:
Ян
2
Лре Дт кл
пДт кл Дт
= *ЬШ Ф (7)+1 _ т ) ДФ
(8)
Для оценки величины отводимого теплового потока использовано приближенное равенство
Яотв ~ к(т0,/ _ Тв ), (9)
где Тв - температура атмосферного воздуха, к -коэффициент теплопередачи через стены:
к = V 12(8,7 X)+1 «в ].
8, - толщина ,-го слоя футеровки, ав - коэффициент теплоотдачи от корпуса печи к воздуху.
Основной недостаток примененного аналитического подхода состоит в упрощенной оценке величины отводимого теплового потока (9).
Проверка адекватности упрощенного аналитического подхода.
Для разработанной математической модели процессов горения и теплообмена в цементной вращающейся печи составлена компьютерная программа и выполнен численный эксперимент [3] с использованием формулы (8) при расчете температуры стен. Вычисления выполнены для печи мокрого способа производства с диаметром корпуса 5 м, толщиной хромомагнезитовой футеровки 230 мм и минимальной расчетной толщиной слоя клинкерной обмазки 20 мм, предохраняющей футеровку от разрушения. Полученные результаты представлены на рис. 2 сплошными линиями.
Температурные кривые построены для тех поперечных сечений печи, в которых температура
клинкера приближенно равна 1300 и 1450 °С. Так как слой клинкера, перемещаясь вдоль печи, сначала нагревается в зоне спекания, а затем, охлаждается, то выбранные значения его температуры соответствуют четырем поперечным сечениям. В зоне спекания внутренняя поверхность стен печи при выходе из-под слоя клинкера нагревается газами до температуры большей, чем температура клинкера, а в зоне его охлаждения - наоборот, температура поверхности стен становится ниже температуры слоя клинкера.
1600
1400
и
о
я —
Р
а
ю С
£ 1200
1000:
4
3
... 1
0 60 120 180 240 300 360 Угол поворота печи, град
Рис. 2. Изменение температуры внутренней поверхности клинкерной обмазки в поперечных сечениях печи в зонах охлаждения (1, 2) и спекания (3, 4) клинкера, линии - по формулам (8), (9), пунктир - численное решение уравнения (2)
На рис. 2 видно, что внутренняя поверхность стен печи при попадании под слой материала испытывает значительные температурные скачки, превышающие 200 К в менее горячих частях печи. По мере приближения к поперечному сечению печи с наиболее горячим клинкером величина температурных скачков уменьшается.
Достоверность результатов, полученных по упрощенному аналитическому подходу, проверена путем численного решения дифференциального уравнения нестационарной теплопроводности (2) по толщине стены при сохранении неизменным температурного поля, определенного в численном эксперименте для газообразной среды. При этом в граничные условия задачи вводилось известное приближенное значение температуры Ту внутренней поверхности слоя клинкерной обмазки, что обеспечило устойчивую сходимость итераций и позволило найти с высокой точностью температуру Т1,1 в ближайшем к поверхности узле сетки.
В результате появилась возможность рассчитать поток отводимого тепла по более точному выражению (4), ввести его значение в уравнение (1) баланса тепловых потоков и вычислить более точное значение температуры поверхности Т0,г . Такой алгоритм расчета повторялся многократно
до тех пор, пока отличие вычисленной температуры от ее значения в предыдущей итерации не становилось менее 0,00001 определяемой величины. Достоверность расчетной температуры внутренней поверхности стен печи в данном случае не вызывает сомнений.
Результаты расчета температуры на поверхности слоя клинкерной обмазки по уточненному численному методу, изображенные на рис. 2 пунктирными линиями, близки к кривым, полученным на основе упрощенного аналитического подхода. Отклонение сплошных линий от пунктирных составляет, как правило, менее 10 °С, что подтверждает адекватность упрощенного аналитического подхода в пределах погрешности инженерных расчетов.
Результаты расчета по уточненному алгоритму показывают особенности неизотермического распределения температуры по толщине стен. На рис. 3 представлены кривые изменения температуры по толщине слоя клинкерной обмазки, охлаждаемой (рис. 3, а) и нагреваемой газами (рис. 3, б). В зоне охлаждения поверхностные слои стен, попадая под слой клинкера, нагреваются, а в зоне спекания клинкера, наоборот, охлаждаются (кривые 4... 6).
После выхода из-под слоя клинкера (кривые 1.3) температурная кривая 1 приобретает сложную волнообразную форму на фоне теплоотдачи к газам тепла, аккумулированного стенами. При всем многообразии форм температурных кривых, представленных на рис. 3, влияние нестационарности, вызывающей колебание температуры внутри стен вращающейся печи, распространяется на относительно небольшую толщину клинкерной обмазки, практически не превышающую 12 мм, оставляя неизменной температуру по толщине футеровки.
Таким образом, при численном моделировании процессов горения и теплообмена в цементных вращающихся печах допустим упрощенный аналитический подход к учету нестационарности теплообмена. Для определения распределения температуры по толщине стен печи возможно численное решение дифференциального уравнения нестационарной теплопроводности при сохранении неизменным температурного поля, определенного в численном эксперименте для газообразной среды.
Выводы
1. В математической модели теплообмена разработан упрощенный аналитический подход к расчету нестационарной температуры стен вращающейся печи, обеспечивающий достоверность результатов в пределах погрешности инженерных расчетов.
2. Стены цементной вращающейся печи испытывают значительные температурные напря-
жения при скачкообразном изменении температуры, превышающем 200 °С при входе их под слой клинкера..
3. Нестационарное изменение температуры распространяется на сравнительно тонкий поверхностный слой, находящийся, как правило, в пределах защитной клинкерной обмазки.
а)
1300
U 1200
cd"
а
£
& поо
с %
(D
н
1000
Расстояние по толщине стен, мм
б) 1600
1500
и
о
си
Р 1400 &
п> с
а
<о
Е- 1300
Расстояние по толщине стен, мм
Рис. 3. Температура в слое клинкерной обмазки в зоне охлаждения (а) и спекания (б) клинкера: 1 - 3 - при контакте с газами, 4 - 6 - под слоем клинкера
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Рязанцев О.А., Кузнецов В.А., Трубаев П.А. Математическая модель теплообмена во вращающейся печи с учетом переизлучения // Математические методы в технике и технологиях ММТТ-23. Т. 8. - Саратов, 2010. - С. 87 - 89.
2. Шорин С. Н. Теплопередача. - М.: Высшая школа, 1964. - 487 с.
3. Кузнецов В.А., Рязанцев О.А., Трулёв А. В. Численное моделирование горения и теплообмена в цементной вращающейся печи // Вестник БГТУ им. В.Г. Шухова. 2011. № 4. - С. 161 - 164.
\2
X
1 5/ ...