УДК 517.927.4 : 536.24 + 66.015.23
В. Н. Лаптинский, д-р физ.-мат. наук, проф., А. А. Романенко, канд. физ.-мат. наук, доц., В. Ю. Стеценко, канд. техн. наук, Р. В. Коновалов
РАСЧЕТ ДИНАМИЧЕСКОГО И ТЕПЛОВОГО ПОГРАНИЧНЫХ СЛОЕВ ДЛЯ ОЦЕНКИ ОХЛАЖДАЮЩЕЙ способности кристаллизаторов
Получено численно-аналитическое решение задачи о динамическом и тепловом пограничных слоях в автомодельном случае. Выполнены расчеты локального коэффициента теплоотдачи от рубашки кристаллизатора к основному потоку охладителя. Проведены экспериментальные исследования охлаждающей способности кристаллизатора для продольного и поперечного течений охладителя.
Введение
Исследования пограничных слоев, начатые Прандтлем в 1904 г., показали, каким образом проявляет себя вязкость при больших числах Рейнольдса и каким путем можно упростить дифференциальные уравнения Навье-Стокса для получения их приближенных решений в случае малой вязкости. Только идея Прандтля о пограничном слое и сделала доступным для теоретического изучения такого рода явления.
Благодаря усилиям многих исследователей к настоящему времени явление пограничного слоя достаточно хорошо изучено (см. [1-10] и приведенную в них библиографию). Разработаны аналитические и численные методы решения соответствующих задач. Однако следует отметить, что для ряда гидродинамических задач, например в задаче об оценке охлаждающей способности кристаллизаторов при непрерывном литье металлов и сплавов, данное явление мало изучено, хотя решение этих задач имеет большое значение для расчета затвердевания отливок. Сейчас такие оценки и расчеты проводятся по сложным численно-полуэмпирическим методикам.
Как показывают исследования, именно пограничный слой охлаждающей жидкости во многом определяет охлаждающую способность кристаллизатора, стабильность и качество литья. Но для изучения влияния этого слоя не-
обходимо разработать конструктивный метод анализа соответствующих математических моделей. Для повышения эффективности исследований следует создать физическую модель кристаллизатора из прозрачных материалов и соответствующую методику проведения модельных исследований. Для струйных кристаллизаторов подобных моделей и методик в мире не существует. Поэтому такие кристаллизаторы необходимо создавать, а также разрабатывать адекватные методики аналитического расчета гидродинамики охлаждения кристаллизатора. Это позволит изучить математическую модель процесса, понять сущность механизма струйного охлаждения и провести оптимизацию конструктивных и технологических параметров с целью создания кристаллизатора с максимальной охлаждающей способностью, что приведет к значительному повышению производительности процесса литья и улучшит качество получаемых заготовок.
Научная идея исследования состоит в развитии методов решения краевых задач [11] и подходов [12-14] применительно к краевым задачам типа Прандт-ля с последующим применением полученных результатов для расчета толщин гидродинамического и теплового пограничных слоев, а также локального коэффициента теплоотдачи для различных режимов охлаждения кристаллизатора и параметров потока охладителя.
Постановка задачи
Слой, в котором отчетливо проявляется действие вязкости и происходит наиболее существенное изменение скорости по нормали к обтекаемой поверхности, называется гидродинамическим пограничным слоем. Толщина динамического слоя - расстояние 5 от обтекаемой поверхности, на котором скорость потока изменяется от скорости движения обтекаемой поверхности до значения, весьма близкого к скорости
обтекающего потока [1-5].
Тепловым пограничным слоем называется пристенная область, в которой существенно проявляются тепловые возмущения. Толщина этого слоя - расстояние 5Т от стенки, на котором температура потока изменяется от Тст до значения, весьма близкого к температуре невозмущенного потока Тш [1-5].
Схема динамического и теплового пограничных слоев представлена на рис. 1.
Рис. 1. Схематическое изображение динамического и теплового пограничных слоев
Толщины динамического 5 и теплового 5Т пограничных слоев связаны соотношением [1, с. 266]
5Т
1
л/Рг
Рг =
V
а
где Рг - безразмерная величина, называемая числом Прандтля; а - коэффициент температуропроводности. Если Рг > 1, то 5Т < 5; если Рг < 1, то 5Т > 5 .
Известно, что для воды, как правило, Рг > 1. Отметим, что Рг ~ 7 для воды при 20 0С. При температурах от 130 до 310 0С значения числа Рг для воды изменяются незначительно и близки к единице. Характер зависимости Рг от температуры резко изменяется при давлениях, близких к критическим.
Поскольку число Прандтля выражает отношение способности жидкости переносить количество движения к ее
способности переносить тепло, то ввиду широких пределов изменения Рг для обычных веществ распределения температуры в конкретном течении могут быть весьма различными в зависимости от рабочего тела. Следует отметить, что в турбулентном потоке механизмы переноса тепла и количества движения приближенно аналогичны, поскольку «турбулентное число Прандтля» (см., например, [1, с. 631; 9, с. 195]) почти всегда близко к единице.
Хотя аналогия между процессами переноса тепла и количества движения обычно не является точной, она всегда достаточно верна, чтобы ее можно было использовать для качественных оценок даже в турбулентном течении, особенно при выборе методов решения уравнений теплообмена. Часто она настолько точна, что безразмерный параметр температуры можно приравнять соответствующему безразмерному параметру скорости, умноженному на «коэффициент аналогии», определяемый на основании эксперимен-
тальных данных [9, с. 23]. В связи с этим широкое применение находят полуэмпи-рические методы, которые состоят в том, что на основе некоторой модели явления теоретически устанавливается структура (общий вид) зависимости между искомыми параметрами, а входящие в нее константы или вспомогательные функциональные связи определяются экспериментально [8, с. 24].
Таким образом, для качественного анализа процесса теплообмена достаточно рассчитать толщину динамического пограничного слоя с учетом «коэффициента аналогии». Количественное определение толщины теплового пограничного слоя проводится на основе уравнения энергии с учетом решения уравнения движения. Толщина теплового пограничного слоя является оценочной характеристикой эффективности работы кристаллизатора [15].
Преобразованные уравнения переноса количества движения и энергии для плоского ламинарного течения в автомодельном случае имеют вид [5, 9]:
/' + ^ //' + т(1 - /'2) = 0; (1)
3" + Рг
2 т +1
/ 3-г/'(3-1)
= 0, (2)
где /(П)- безразмерная функция тока
(фактически / = /(г/, т) ), /'(П= 4(П/П П - безразмерная автомодельная пере-
менная,
П = У
и (х )
V xV У
(3)
3 - безразмерная температура, Т„ - Т
3
Т - Т '
ст да
(4)
Здесь и(х) = С хт ; Тст - заданная температура стенки; То — температура невозмущенного потока; Т - температура
в пограничном слое, при этом Т_ = Тда + ах7;
" ст да
(5)
краевые условия имеют вид:
I(0) = I'(0) = 0, /'() = 1; (6)
5(0)= 0, 5(оо) = 1. (7)
Заметим, что значения параметров С, т должны быть согласованы соответствующим образом.
Решение задачи
На основе подхода [11, гл. 1] разработан численно-аналитический метод решения задачи (1), (6). Согласно этому методу решение указанной задачи представимо в виде
Ч
/ (П = \
-} р(
йа
да -| р(а)йа
| е 0 йт
йт, 0 <^<да.
где р(П)- вспомогательная функция, определяемая как решение задачи
йр
ч
= р2 - аЛр\
(а \
| р(
0 4а
4т +
п №
+ (а - 2Ь)Д[ е 0 йт,
0
р(0)=Л,
при этом значение параметра Л находится из трансцендентного уравнения
т
да -| р(а)4а
Ле 0 йт-1 = 0;
т +1 7
где а =-------; Ь = т .
2
Заметим, что структура решения такова, что оно удовлетворяет краевым условиям (6). С помощью этого метода получены приближенные решения с различной погрешностью (от 2 до 0,1 %). Одно из приближенных решений (до 2 %) имеет вид:
f (п) = Ц jexP
О О
Л l! +UJ 2! +
( f І \3 І 3 \Л
2(a] +(a-2b)Л a
3!
da) dr, (В)
JJ
где значение параметра Л находится из уравнения
Aj exp
Г Г b a Г b 12 a2
О
Г Г b 1
+ 2
f Л J
A l!
+ (a - 2b )A
fA j
2!
V
, 3! y
J yJ
dr -1 = О.
На основе приближенного решения задачи (1), (6) в случае Тст = const
(у = 0 в (5)) соответствующее решение задачи (2), (7) получено в виде
& = -
e x 43 n m+1n(^ Pr • 2 G GO
e x 43 n m +1„/,.\ Pr • 2 G GO
(9)
ч
где g (n)=jf az )dz.
Функция О(п) получена численно и как решение задачи
О{1Г) + 0'О"' + т (1 - О"2) = 0;
О(0) = 0'(0) = 0"(0) = 0, 0"(да) = 1, с помощью разработанного численно-
аналитического метода, при этом в качестве функции / (г) можно взять ее
приближение (8).
Используя (9), из (4) находим распределение температуры в тепловом пограничном слое, при этом толщина теплового пограничного слоя и локальный коэффициент теплоотдачи от рубашки кристаллизатора к основному потоку охладителя даются соотношениями типа [5, с. 172, 175]
ST (x) = пТ (m)
XT І
vU (x h
1
/ 1-m
= пт(m)[T\ x 2;
(1О)
а
(x) = k (m) A
г U (x)l
:k (m )A (£ ]2 xt
(11)
где Пт (т) - значение безразмерной переменной п, определяемое на основе значений функции Э(П из (7) с учетом (3) в соответствии с принятой точностью (З(г/Т ) = 0,99); А0 - коэффициент теплопроводности охладителя (воды);
k (m) = •
1
exp
- Pr •
m +1
ъ
j f (z )dz
фактически а = а( х, т,А0,у, С).
Замечание - Параметр т связан с углом в потока охладителя соотношением [5, с. 177]
m =
в
2п - в
Применительно к параметру в
да
имеем
( л-в
8Т (х ) = П (в)[^ ]2 х-в;
а
(х) = к (в) А
где х > 0,
к (в) = к
V) ^ в
2, - в
[ в ' 2п-в,
1
в-п Ля-в
л
I ехр
Рг-п %
2п- в
Ниже приведена табл. 1 значений П, соответствующих различным углам
атаки Р = ~ [9, с. 108] и некоторым характерным значениям Тст.
В качестве примера приведем графики 5Т ( х) при характерных значениях
параметров в, С и Тст.
В табл. 2 приведены значения к , соответствующие различным углам атаки р и некоторым характерным значениям Т .
ст
В качестве примера приведем графики для а( х) при характерных значениях параметров в, С, Тст и А0.
Табл. 1. Значения ПТ (в)
вп ПТ (Т„ = 100 0С) ПТ (Тст = 120 0С) ПТ (Тст = 140 0С)
0,0 3,900 4,285 4,895
0,1 3,755 4,040 4,630
0,2 3,560 3,840 4,410
0,3 3,395 3,660 4,215
0,4 3,240 3,500 4,035
0,5 3,095 3,345 3,865
0,6 2,955 3,200 3,700
0,7 2,820 3,055 3,535
0,8 2,690 2,910 3,370
0,9 2,555 2,765 3,210
1,0 2,420 2,620 3,040
Табл. 2. Значения к (в)
в, к (Тст = 100 0С) к (Тст = 120 0С) к (Тст = 140 0С)
0,0 0,4033 0,3770 0,3332
0,1 0,4345 0,4056 0,3574
0,2 0,4624 0,4311 0,3791
0,3 0,4889 0,4555 0,3998
0,4 0,5152 0,4797 0,4205
0,5 0,5420 0,5044 0,4416
0,6 0,5700 0,5301 0,4638
0,7 0,5998 0,5576 0,4874
0,8 0,6320 0,5873 0,5129
0,9 0,6674 0,6199 0,5411
1,0 0,7070 0,6565 0,5726
А
х
Рис. 2. Зависимости 8Т (х) при Тст = 120 0С (Рг = 1,44), и = 0,246-10-6 м2-с-1 и различных значениях параметров в, С : 1 - в = 0, С = 1,6 с-1; 2 - в = ^, С = 4,7 с-1; 3 - в = я, С = 15,8 с'1
х ---------->-
Рис. 3. Зависимости а(х) при Тст = 120 0С (Рг = 1,44), и = 0,246-10-6 м2-с-1, Х0 = 0,677 Вт/(м-К) и различных значениях параметров в, С : 1 - в = 0, С = 1,6 с1; 2 - в = С = 4,7 с-1;3 - в = я, С = 15,8 с’1
Экспериментальные исследования коэффициента теплоотдачи
Для определения коэффициента теплоотдачи в кристаллизаторе при различных углах обтекания охладителя был изготовлен измерительный комплекс. Его общий вид представлен на рис. 4. Измерительный комплекс включает: экспериментальную отливку диаметром 70 мм высотой 160 мм из чисто-
го алюминия марки А7; в средней части заготовки установлены три термопары с открытым спаем диаметром 0,5 мм (схема расположения термопар представлена на рис. 5); устройство затоп-ленно-струйного охлаждения отливок (УОО), позволяющее подавать поток охладителя под углами 0 и 900 к охлаждаемой поверхности; блок марки ЬО регистрирования температуры и ноутбук.
Рис. 4. Общий вид измерительного комплекса
_ />~Л 3 2 1
Отливка Охлаждающая
жидкость
Рис. 5. Схема расположения термопар в экспериментальной отливке: 1\ = 5 мм, 12 = 5 мм, 13 = 5 мм
Эксперименты проводились на опытной установке литья закалочным затвердеванием. Заготовки нагревали до температуры 120 0С в печи сопротивления марки «8по1», после чего помещали в устройство охлаждения отливок (УОО). Расход охладителя составлял 15 м3/ч. Опыты проводились при двух различных углах подачи охладителя: 0 и 900. По результатам экспериментальных данных построены зависимости распределения температуры в отливке от времени. Установлено, что при охлаждении цилиндрической заготовки из алюминия в УОО распределение температуры вблизи наружной поверхности отливки можно принять линейным. На основании этого приближения из граничного условия 3-го рода на границе отливка-охлаждающая жидкость получена формула для расчета коэффициента теплоотдачи
а = -Л
Т - Т
11 2
1
Ті - у (Т2 - Т )-Тж
1-і
где X - теплопроводность алюминия, X = 213 Вт/(м-К); Т1, Т2 - показания термопар в отливке в точках 1 и 2; Тж -температура основного потока охладителя, 11 = 5 мм, 12 = 5 мм (см. рис. 5). Среднее значение результатов расчета по данной формуле коэффициента теплоотдачи для углов обтекания 0 и 900 составило 3000 Вт/(м2-К), что согласуется с результатами теоретических расчетов. Это свидетельствует о приемлемости выбранной математической модели, корректности полученных теоретических результатов.
Выводы
Таким образом, в результате проведенных исследований получено численно-аналитическое решение задачи о динамическом ламинарном пограничном слое в автомодельном случае. На основе этого решения выполнены рас-
четы толщины теплового пограничного слоя, а также значений локального коэффициента теплоотдачи от рубашки кристаллизатора к основному потоку охладителя при различных способах охлаждения кристаллизатора и параметрах потока охладителя. Выполнены экспериментальные исследования охлаж-
дающей способности кристаллизатора для продольного и поперечного течений охладителя, которые показали приемлемость выбранной математической модели и корректности полученных теоретических результатов. Некоторые результаты проведенных исследований опубликованы в [16-19].
СПИСОК ЛИТЕРАТУРЫ
1. Шлихтинг, Г. Теория пограничного слоя / Г. Шлихтинг. - М. : Наука, 1974.
2. Основы теплопередачи в авиационной и ракетно-космической технике / В. С. Ав-дуевский [и др.]. - М. : Машиностроение, 1975.
3. Лыков, А. В. Тепломассообмен : справочник / А. В. Лыков. - 2-е изд., перераб. и доп. / А. В. Лыков. - М. : Энергия, 1978.
4. Кутателадзе, С. С. Основы теории теплообмена / С. С. Кутателадзе. - М. : Атомиз-дат, 1979.
5. Теория тепломассообмена. Учебник для вузов / Под ред. А. И. Леонтьева. - М. : Высш. шк., 1979.
6. Войткунский, Я. И. Гидромеханика : учебник / Я. И. Войткунский, Ю. И. Фаддеев, К. К. Федяевский. - Л. : Судостроение, 1982.
7. Теория турбулентных струй / Под ред. Г. И. Абрамовича. - М. : Наука, 1984.
8. Емцев, Б. Т. Техническая гидромеханика / Б. Т. Емцев. - М. : Машиностроение, 1987.
9. Себиси, Т. Конвективный теплообмен / Т. Себиси, П. Брэдшоу. - М. : Мир, 1987.
10. Репик, Е. У. Турбулентный пограничный слой. Методика и результаты экспериментальных исследований / Е. У. Репик, Ю. П. Соседко. - М. : Физматлит, 2007.
11. Лаптинский, В. Н. Конструктивный анализ управляемых колебательных систем /
В. Н. Лаптинский. - Минск : ИМ НАН Беларуси, 1998.
12. Лаптинский, В. Н. Матрица типа Коши для нелинейных дифференциальных систем / В. Н. Лаптинский // VIII Белорус. мат. конф. : тез. докл. Междунар. конф. - Минск : ИМ НАН Беларуси, 2000. - Ч. 1. - С. 125.
13. Лаптинский, В. Н. Об одном представлении решений нелинейных дифференциальных систем / В. Н. Лаптинский // Дифференциальные уравнения и системы компьютерной алгебры : тез. докл. Междунар. матем. конф.- Брест : БрГУ, 2000. - С. 42-43.
14. Лаптинский, В. Н. К методу интегральных уравнений в теории краевых задач для дифференциальных уравнений / В. Н. Лаптинский // Математическое моделирование и дифференциальные уравнения : тез. докл. Первой Междунар. конф. - Минск : ИМ НАН Беларуси, 2007. - С. ВЗ.
1З. Марукович, Е. И. Повышение эффективности работы кристаллизатора при литье слитков /Е. И. Марукович, В. Ю. Стеценко // Литье и металлургия. - 200З. - № 2, Ч. 1. -С. 139-141.
1б. Лаптинский, В. Н. О методах анализа граничной задачи Блазиуса-Хартри / В. Н. Лаптинский, Е. И. Марукович, В. Ю. Стеценко // Математическое моделирование и дифференци-
альные уравнения : тез. докл. Первой Междунар. конф. - Минск : ИМ НАН Беларуси, 2007. -
С. 34-3б.
17. Стеценко, В. Ю. Методика расчета гидродинамического пограничного слоя при охлаждении кристаллизатора / В. Ю. Стеценко, В. Н. Лаптинский, А. А. Романенко // Материалы, оборудование и ресурсосберегающие технологии : материалы Междунар. науч.-техн. конф. - Могилев : Белорус.-Рос. ун-т, 2007. - Ч. 1. - С. 97.
1В. Лаптинский, В. Н. Расчет охлаждающей способности кристаллизатора при непрерывном литье / В. Н. Лаптинский, В. Ю. Стеценко, Р. В. Коновалов // Материалы, оборудование и ресурсосберегающие технологии : материалы Междунар. науч.-техн. конф. - Могилев : Белорус.-Рос. ун-т, 2009. - С. 7З.
19. К задаче о тепловом пограничном слое / В. Ю. Стеценко [и др.] // Математическое моделирование и дифференциальные уравнения : тез. докл. Второй Междунар. науч. конф. -Минск : ИМ НАН Беларуси, 2009. - Ч. I. - С. В0.
Институт технологии металлов НАН Беларуси Материал поступил 05.07.2010
V. N. Laptinsky, A. A. Romanenko,
V. J. Stetsenko, R. V. Konovalov Calculation of dynamic and thermal boundary layers for estimation of cooling ability of crystallizers
The numerically-analytical solution of the problem of dynamic and thermal boundary layers in a selfsimilar case has been obtained. Calculations of the local factor of heat-transfer between a ciystallizer jacket and main coolant stream have been performed. Experimental researches of cooling ability of the crystallizer for longitudinal and cross coolant flows have been made.