Идентификация критериального уравнения теплоотдачи экспериментальным данным при охлаждении алюминиевых слитков Текст научной статьи по специальности «Строительство и архитектура»

УДК 621.78

ИДЕНТИФИКАЦИЯ КРИТЕРИАЛЬНОГО УРАВНЕНИЯ ТЕПЛООТДАЧИ ЭКСПЕРИМЕНТАЛЬНЫМ ДАННЫМ ПРИ ОХЛАЖДЕНИИ АЛЮМИНИЕВЫХ СЛИТКОВ

© 2013 А. С. Горшенин, В. Е. Кривошеев Самарский государственный технический университет

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

Критериальное уравнение, скорость охлаждения, аппроксимация, камера охлаждения.

Введение

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

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

Математическая модель теплообмена в камере охлаждения

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

ние многие режимные и конструктивные факторы. Для оценки влияния данных факторов на конвективный теплообмен в камере охлаждения алюминиевых слитков авторами предложена математическая модель конвективного теплообмена в системе «охлаждающий воздух - алюминиевые слитки» [1]:

tВN - + тТ~^(N-1))"

СРе В

(N-1) 'в( N-1);

-(То-1 В1) ехр(- 4 • Ы •

Тм - tв N. + (То -1В N. ) ехр(- 4 • Ы • М

— < 700С/ч в интервале 500-3000 С. г

Здесь tВN - температура охлаждающего

воздуха при омывании N - го слитка, 0 С ; tв - температура охлаждающего воздуха при омывании предыдущего N - му слитка, 0 С ; сРал - коэффициент теплоемкости алюминия, Дж /(кг С); шст - масса слитка, кг; сРв - коэффициент теплоемкости воздуха, Дж /(кг •0 С) ; GВ - расход воздуха, м3 / ч ; Т0 - начальная температура охлаждения слитка, С ; Тх - температура N - го слитка, 0 С ; Т'(N-1) - температура предыдущего N - му слитка, 0 С; Ы - критерий Био; ¥о -

критерий Фурье; Т - скорость охлажде-т

ния слитка, 0 С / ч .

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

Решение данной математической модели предполагает определение коэффициента теплоотдачи а . Движение воздуха происходит по каналам, ограниченным цилиндрическими слитками. Уравнения для расчёта коэффициента теплоотдачи при турбулентном режиме движения при омывании волновой поверхности даны в работах [2, 3]. В работе В.И. Гоме-лаури [4] предложено критериальное уравнение теплоотдачи при омывании волновой поверхности:

Шж^ = 0,021-Ке^ • РгГГ

0,25

(1)

— а • й,

где Ш жн =

1

среднее значение

критерия Нуссельта по длине канала; 1 -коэффициент теплопроводности,

Вт/ к ; йэкв - эквивалентный диаметр

Ж • й _

м

канала, м; Яе

жё э

среднее

значение критерия Рейнольдса по длине канала; Ж - скорость охлаждающего воздуха, м/с ; V - коэффициент кинематиче-

ской вязкости,

м

; Ргж - среднее зна-

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

= 1,04 • Рг^04 • ехр

0,85 •Г (% )

~ шаг волн, £ - расстояние между

соседними волнами, м; И - высота волновой поверхности, м .

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

Описание эксперимента

Экспериментальное исследование теплообмена при охлаждении слитков проводилось на смонтированной камере охлаждения. Камера охлаждения представляет собой конструкцию, состоящую из рабочего пространства, в котором располагается садка; раздающего и собирающего газоходов. Для подачи охлаждающего воздуха на раздающем газоходе установлены четыре осевых вентилятора. При проведении эксперимента определялись зависимость температуры слитков в течение времени охлаждения, скорости охлаждающего воздуха. Температура слитков замерялась термопарами градуировки ХА, зачеканенными в слитки. Вторичным прибором служил измеритель -регулятор температуры 2ТРМ1. Скорость охлаждающего воздуха определялась путём пересчёта динамического давления, измеренного напорной трубкой. При проведении экспериментальных исследований погрешность измерений температуры слитков составляла 1,5%, скорости охлаждающего воздуха ± 2,5%.

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

е

ш

ш

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

время охлаждения,ч

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

16 14 12 10

о а.

0,1 0,3 0,5 1 2 г 4 Б 6 7 3 9 10 11 12 13 14 15 16 17 13 19 20

время охлаждения, ч

Рис. 2. Изменение скорости охлаждающего воздуха в течение времени охлаждения

Идентификация критериального

уравнения теплоотдачи по методу наименьших квадратов

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

Для аппроксимации экспериментальных данных степени соответствия между данными эксперимента и результатами определения коэффициента теплоотдачи по формуле (1) воспользуемся методом наименьших квадратов [5].

Экспериментальные данные аппроксимируем к уравнению, выражающему зависимость между критериями Нуссельта Ыи , Рейнольдса Яе и Прандтля Рг :

Ыи = А ■ Яев ■ Ргс ,

где A, B, C - постоянные числа. Для упрощения расчётов опустим значение критерия Прандтля PrС, так как в интервале температур 0 - 6000 C его значение почти не меняется.

Функция, подлежащая минимизации для определения оценок параметров методом наименьших квадратов, имеет вид [6]

F(A, B) = ¿ (Nusг - A • ReB) , (2)

i=1

где Nudг - значение зависимой переменной в г -м эксперименте, A • ReB - значение независимой переменной, А, В - неизвестные коэффициенты, подлежащие определению.

Для определения неизвестных коэффициентов в (2) вычислим частные производные по А и В и приравняем их нулю:

2

— = — t(Nus - A • Re $ ) = 0, дА dAtl{ d i s/

— = ~1t (Nus - A • ReB ) = 0.

дБ дБ di

После проведения математических преобразований данных уравнений окончательно получим систему двух уравнений с двумя неизвестными:

¿(Nust • ReB )

А = —-,

Ê(Re2f )

г=1 (3) ±(Nusг • ReB-1 ) Б = -.

¿(ReB • Ref1 )

i=1

Решая уравнение (3) методом подстановки, получим значения коэффициентов: А = 0,352 и Б = 0,65 . Подставим данные коэффициенты в уравнение Nu s = 0,352 • Red °'65.

С учётом критерия Прандтля окончательно получим

Nus = 0,352 • Red 0,65 • Pr0,33. (4)

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

Заключение

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

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

Библиографический список

1. Горшенин, А.С. Исследование тепловых процессов при охлаждении алюминиевых слитков после их термической обработки [Текст]/ А.С. Горшенин,

B.Е. Кривошеев // Энергетика и энергоэффективные технологии: тр. V между-нар.науч.-практ. конф. — Липецк, 2012. —

C. 60-62.

2. Кружилин, Г. Н. Теория теплопередачи круглого цилиндра в поперечном потоке жидкости [Текст] / Г. Н. Кружилин // ЖТФ. — 1936. — Т. 6. — Вып. 5. — С. 858865.

3. Eckert, E. R. G. Die Berechnung des Warmeudergangen in der laminaren Grenzschicht umstromter Korper [Text] / E. R. G. Eckert // VDI — Forschungsh, 1942. —S. 1-26.

4. Гомелаури, В. И. Влияние искусственной шероховатости на конвективный теплообмен [Текст] / В. И. Гомелаури // Тр. института физики АН ГССР. — Тбилиси. — 1963. — Т. 9. — С. 111-145.

5. Горский, В.Г. Планирование промышленных экспериментов [Текст] I В.Г. Горский, Ю.П. Адлер, А.М. Талалай.— М.: Металлургия, 1978. — 279 с.

6. Соколов, А.К. Коэффициенты, упрощающие расчет теплопередачи через обмуровки печей [Текст] I А.К. Соколов. — Иваново: Ивановский государственный университет, 1980. — 11 с.

AGREEMENT BETWEEN CRITERIAL EQUATIONS OF HEAT TRANSFER IN THE PROCESS OF COOLING ALUMINIUM INGOTS AND EXPERIMENTAL DATA

© 2013 A. S. Gorshenin, V. Е. Krivosheev Samara State Technical University

The paper deals with the problem of agreement between the criterial equation for the determination of the heat transfer coefficient and the experimental data for the case of cooling aluminum ingots using the least squares method. An experimental unit for the cooling of aluminum ingots is described. Experimental data are presented.

Criterial equation, rate of cooling, approximation, cooling chamber.

Информация об авторах

Горшенин Андрей Сергеевич, заместитель заведующего кафедрой «Промышленная теплоэнергетика», Самарский государственный технический университет. Email: andersonag1@yandex.ru, andersonag 1@mail.ru. Область научных интересов: изучение и оптимизация процессов теплообмена в конвективных печах для термической обработки алюминия.

Кривошеев Владимир Евгеньевич, кандидат технических наук, доцент кафедры «Промышленная теплоэнергетика», Самарский государственный технический университет. E-mail: pt@samgtu.ru. Область научных интересов: изучение теплообмена в печах и котельных установках.

Gorshenin Andrey Sergeyevich, deputy head of the department of industrial heat power engineering, Samara State Technical University. E-mail: andersonag1@yandex.ru, andersonag1@mail.ru. Area of research: study of heat exchange in convective furnaces for aluminum thermal treatment and optimization of the process.

Krivosheev Vladimir Evgenievich, candidate of technical science, associate professor of the department of industrial heat power engineering, Samara State Technical University. Email: pt@samgtu.ru. Area of research: study of heat transfer in furnaces and boilers.