CC BY
22
ФИЗИКА
УДК 532.517.2
ТЕПЛООБМЕН ПРИ ЛАМИНАРНЫХ ТЕЧЕНИЯХ ВЯЗКОЙ ЖИДКОСТИ В КАНАЛАХ С ТРЕХЗАХОДНЫМИ ШНЕКОВЫМИ ВСТАВКАМИ ПРИ БОЛЬШИХ
ЧИСЛАХ ПЕКЛЕ
И.А. Конахина1, Э. М. Хуснутдинова1, А.И.Кадыйров2, Е.К.Вачагина2
1 Казанский федеральный университет, Казань, Россия 2 ФИЦ КазНЦ РАН, Казань, Россия
Резюме: В работе представлены математическая модель и результаты численных исследований теплообмена при ламинарных течениях вязких жидкостей в каналах с трехзаходными шнековыми вставками, которые широко распространены на практике в качестве интенсификаторов конвективного теплообмена.
Ключевые слова: гидродинамика, теплообмен, шнековые вставки, винтовая система координат.
Вв1:10.30724/1998-9903-2018-20-11-12-161-169.
Для цитирования: Конахина И.А., Хуснутдинова Э.М., Кадыйров А.И., Вачагина Е.К. Теплообмен при ламинарных течениях вязкой жидкости в каналах с трехзаходными шнековыми вставками при больших числах Пекле // Известия высших учебных заведений. ПРОБЛЕМЫ ЭНЕРГЕТИКИ. Т. 20. № 11-12. С. 161-169.DOI:10.30724/1998-9903-2018-20-11-12-161-169.
HEAT TRANSFER FOR LAMINAR VISCOUS FLOWS IN A TUBE FITTED WITH TRIPLE TWISTED FINS ON ROD UNDER LARGE PECLET NUMBERS
I.A. Konakhina',, E.M. Khusnutdinova1,, A.I. Kadyirov2, E.K. Vachagina2
1Kazan Federal University, Kazan, Russia
2FRC Kazan Scientific Center of RAS, Kazan, Russia
Abstract: The mathematical model and the numerical results of the heat transfer in laminar viscous flows in a tube fitted with triple twisted fins on a rod is obtained. The tubes fitted with triple twisted fins on a rod are widely used for heat transfer enhancement.
Keywords: hydrodynamics, heat exchange, screw inserts, helical coordinate system.
For citation: I.A. Konakhina, E.M. Khusnutdinova, A.I. Kadyirov, E.K. Vachagina. Heat transfer for laminar viscous flows in a tube fitted with triple twisted fins on rod under large Peclet
161
numbers. Proceedings of the higher educational institutions. ENERGY SECTOR PROBLEMS vol. 20. № 11-12. pp. 161-169. DOI: 10.30724/1998-9903-2018-20-11-12-161-169.
Введение
Как известно, винтовая закрутка потока жидкости по всей длине канала является одним из пассивных методов интенсификации теплоотдачи при течении вязких жидкостей в каналах теплообменного оборудования [1-6]. В случаях, когда не допускается сильное повышение гидравлического сопротивления канала для интенсификации теплообмена, используется метод разрушения пристенных слоев жидкости с помощью нанесения винтовой накатки [3] или использование коротких закручивающих устройств [6]. Однако для высоковязких жидкостей, для которых характерен ламинарный режим течения, эффект от действия винтовой накатки или коротких закручивающих устройств невысок. В таких случаях и когда допускается увеличение гидравлического сопротивления каналов используется закрутка всего потока жидкости. В качестве устройств, обеспечивающих закрутку всего потока жидкости, а не только ее пристенных слоев, наиболее часто используются винтовые ленточные и шнековые вставки. Исследованию применения ленточных вставок для интенсификации теплообмена посвящено большое количество работ [5]. Однако очевидно, что процессы переноса тепла в каналах с ленточными вставками менее эффективны, чем в каналах со шнековыми вставками. Площадь поверхности теплообмена, приходящаяся на единицу длины канала, в первом случае меньше, чем во втором. Работ, посвященных процессам переноса тепла в каналах со шнековыми вставками, несколько меньше, причем почти все они рассматривают однозаходные шнековые вставки [1, 4, 6]. При этом очевидно, что влияние закручивания потока на процессы переноса тепла в каналах с однозаходной шнековой вставкой меньше, чем в каналах с двузаходной или трехзаходной вставками: чем больше число ходов, тем интенсивнее происходят процессы переноса тепла. С другой стороны, при увеличении числа ходов возрастает гидравлическое сопротивление каналов. Таким образом, в случаях, когда запасы мощности на прокачку жидкости достаточны для обеспечения повышенного гидравлического сопротивления канала, использование трехзаходной шнековой вставки позволит достичь большей эффективности по сравнению с однозаходной шнековой вставкой. Этот факт диктует целесообразность исследования процессов переноса тепла в каналах с трехзаходными шнековыми вставками.
Постановка задачи
Рассматривается теплоперенос при стационарном ламинарном течении вязкой ньютоновской жидкости в канале с трехзаходной шнековой вставкой (рис. 1). Исследуемый канал имеет три равнозначные области, в связи с чем в качестве расчетной области рассмотрена одна из них. В оставшихся двух областях распределения скоростей и
температур будут идентичными.
1
Рис. 1. Продольное сечение канала с трехзаходной шнековой вставкой: 5 - шаг винтового канала (длина части канала, соответствующая повороту ребер шнека на 360 градусов); 1 - ребра трехзаходной шнековой вставки; 2 - внутренняя труба
В качестве основных допущений при постановке задачи принято: теплофизические характеристики (плотность р f, коэффициент теплопроводности X f, удельная
теплоемкость cf и динамическая вязкость ц f) жидкости в рассматриваемом диапазоне
температур меняются незначительно и могут быть приняты постоянными. Силы тяжести пренебрежимо малы, перенос тепла за счет теплопроводности в осевом направлении много меньше, чем в поперечных направлениях, и меньше конвективного переноса тепла как в осевом, так и в поперечных направлениях. На входе в канал задан сформировавшийся профиль скорости, на стенках канала принято условие прилипание жидкости. На границах канала поддерживается постоянное значение температуры.
Канал с трехзаходной шнековой вставкой обладает винтовой симметрией. Данное свойство позволяет ввести специальную винтовую систему координат [1-3], переменные ~ 12 3
которой r , r , r связаны с переменными декартовой системы координат соотношениями
У
1 /2 3 \ = Rr cos (r + k\ I
= Rrr1 sin (
Л2 ±к\31
(1)
г = Яг\
где к = 2пЯ / 5 , Я - больший радиус канала с трехзаходной шнековой вставкой; 8 - шаг винтового канала; г = Яг - ось канала. В (1) знак "+" выбирается для канала, винтовые ребра которого закручены по часовой стрелке, а знак "-" - для канала, винтовые ребра которого закручены против часовой стрелки.
Система уравнений, описывающая процессы гидродинамики и теплообмена при сделанных допущениях принимает вид
/ п? 1 /л'1
fa 1 а.
cv 1 cv 2 1 2 2 „ 123 2133 -v1 +-v - rvv + 2k\1v2v3 -K2r1v3v3
cr1
л
Ф SsV ds]2
-л 1 -л 1 -л 2 1
cr o\ o\ r
1 2 , K 3 11
— s2 ±"T S23 + -г S1 , r1 r1
(2)
ReR (r1) Rer (r1)
^я.2 я.2 0 A
Cv 1 Cv 2 212 ^v1 + — v2 + —v1v2
Cr Cr r
2K 1 3 +—vv+K
■n
Cp CS21 Cs2
^ dv' 1 Cv3
v+v
Crp Crp
//
12 _...13
Cr
(
- + —- +—Y + r[Sl +Krisi
(3)
Cr1 Cr
Re R
+к
M2
f dv2 1 Cv2 2.2 1 2 _ 2k 13 v H--—v H—-v v +—-V V
8r\ r| r|
лл
+ Re
R
1 + K
Cri1
(r1)
2\U .3
Cv3 1 Cv3
___ v1 + v2
Cr1 Cr2
dp C53 953 к i _ i
2 1 3
1 1
3+—j" + ——'s2 +Krl^i +к ri Sl +—53
Cr Cr Cr'
r
(4)
Рея
ду1 ду2 1 1 л -Г +д"Г + -Ту1 = 0 :
дл дл л
( 1 де 2 де 3 де л
у1 —Т+у2 —2 + у3 —з
дл дл дл
( (
д2 е
.(V)
д
2 +5^2
1 + К
(л)
2 ^
де
(л1 )2
дл2
1 де
+' (5) л дл
где Яея =р/УаЯ / М-/ - модифицированное число Рейнольдса, рассчитанное по среднерасходной скорости ¥а = Q / О и радиусу канала: Q - объемный расход жидкости;
О - площадь поперечного сечения; V = V1 / ¥а (1 = 1,2,3) - безразмерные контравариантные компоненты безразмерной скорости; V = V / Va ; V1 (1 = 1,2,3) - контравариантные компоненты скорости V в винтовой системе координат; р = РЯ / (ц/Уаа) - безразмерное давление; Р - давление; Рея =р/с^аЯ / X/ - модифицированное число Пекле; е = (Т - Т^) / (То - Т^) - безразмерная температура; Т - температура жидкости; -температура стенок канала; То - температура жидкости на входе в канал; = (¿га^)^ -
тензор.
В качестве условий однозначности взяты следующие условия:
2
V1 (б, л2 ) = V (1, -л2 ) = 0, (1 = 1,2,3) при -л/3 <л2 <л/3 , V1 (л1, +л /3) = 0 (г = 1,2,3) при 8<-1 < 1,
е(б,л2,л3) = е(1,л2,л3) = о,-л/3<л2 <л/3,о<-3 <»,
е(л\ +л/3,л3) = о,б<л1 < 1,о<л3 <■»,
е(л\л2,о) = 1,б<л1 < 1, -л/3 < л2 <л/3,
е(л!,л2,л3 ^<») = о,б<л1 < 1,-л/3<л2 <л/3, где б = г / Я ; г - меньший радиус канала с трехзаходной шнековой вставкой.
(6)
1 2
Компоненты скорости зависят только от двух переменных лил, а безразмерное
= / = сош!.
давление может быть записано как р = /о л3 + Уоо (л1, л2), где др / дл3 = /о
При сделанных допущениях система уравнений переноса количества движения и неразрывности решается отдельно от уравнения переноса энергии, затем полученные гидродинамические поля используются для решения уравнения переноса энергии. Для решения поставленной задачи адаптирован программный пакет "СОМБОЬ МиШрИу^ча'", использующий метод конечных элементов. Уравнения (2)-(5) записаны в винтовой системе координат. "СОМБОЬ МиШрИу^ча'" позволяет ввести в общую форму уравнений движения и неразрывности необходимые функции для перехода к винтовой системе координат (1).
Апробация пакета на примере решения задачи о течении вязкой жидкости в канале с однозаходной шнековой вставкой представлена в работе [2].
При обработке полученных результатов использованы традиционные числа Рейнольдса Яе^ =р/ ц^ и Пекле Ре^ =р/ X^ , где в качестве
определяющих размеров взят гидравлический диаметр В^ . Его значение для каналов с трехзаходной шнековой вставкой определено как В^ = 4У/2, где объем V и площадь смоченной поверхности 2 одного хода канала и одного витка канала определены следующим образом:
V = -521 Я2 Б/3,
2 = RS
2л(1 + 8) ( I о о arcsinh(к)-arcsinh(к8) —Ь-'-+ VK2 +1 -8л/к282 +1 +--
Безразмерные среднерасходная Qa и нормализованная »(л1,Л2,Л3) температуры
вычислены по формулам:
Qa = R2 Ц6(Л1' Л2 ) U3 (Л1, Л2 ^ Л2
(8)
□о
»(л1, л2, л3 ) = е(г Л2, Л3 )/Qa (л3 )
где »(л1, л2, Л3)
стабилизированному распределению (л1,л2) = lim »(л1,Л2,Л3) •
(9)
с увеличением расстояния от входа в канал стремится к
2\ п/Л„2
3 Л
Локальное Ыи^ и среднее Кий числа Нуссельта определены через среднерасходную температуру следующим образом:
Nuh = РеА
DK _L
4r e„
se.
бл
Nuh = PehDh |ln eal •
Результаты
Расчеты представлены для частного случая: геометрические размеры канала 8 = r / R = 0,5 и S / D = 4, где D = 2R - больший диаметр канала с трехзаходной шнековой вставкой; для двух значений чисел Рейнольдса Reh = 100; 500 и одного значения числа Пекле Peh = 100.
На рис. 2 показаны профили осевых физических составляющих вектора скорости. С ростом числа Рейнольдса картина распределения составляющих вектора скорости становится все более ассиметричной, при этом максимальные значения осевой составляющей смещаются по ходу закрутки потока.
Reh = 100 Reh = 500
Рис.2. Профили осевых составляющих скорости вдоль прямых л2 = ConSt:
, 2 4л , о 2л . , 2 2л , о 4л
1 - л2 =--; 2 - п2 =--; ^ - Л = 0 ; 4 - л2 = — ; J - л2 =---
15 '15 15 ' 15
На рис. 3 показаны профили поперечных радиальной и тангенциальной физических составляющих вектора скорости для различных значений чисел Рейнольдса. Видно, что величины поперечных составляющих вектора скорости соизмеримы с величинами осевой составляющей, следовательно, существенно влияют на процессы переноса тепла в рассматриваемом канале.
ЯеА = 100 ЯеА = 500
Рис. 3. Профили поперечных составляющих скорости вдоль прямых = сош!:
сплошная линия - радиальная составляющая ^и1 /д/о1"1 ^; пунктирная линия -
- тангенциальная составляющая I и /А а22 ; обозначения согласно рис. 2
На рис. 4 показано распределение безразмерной среднерасходной температуры в по
длине канала. Видно, что с увеличением числа Рейнольдса температура жидкости приближается к температуре стенки канала на меньшем расстоянии: для Яе^ = 500 это
расстояние равно 12, а для Яе^ = 100 равно 16.
Рис.4. Распределение безразмерной среднерасходной температуры по длине канала
Распределения локальных и средних чисел Нуссельта по длине канала со шнековыми вставками характеризуются резким падением вблизи области, непосредственно прилегающей к входу в канал, затем ростом и стабилизацией (рис. 5), что качественно согласуется с результатами работы [7]. При этом локальные числа Нуссельта принимают наименьшие значения вблизи входной области канала Nu h =5,51 ( Re h = 100 ) и Nu h =6,72 ( Re h = 500 ). Увеличение числа Рейнольдса в 5 раз приводит к увеличению
стабилизированного числа №^ в 1,44 раза. Несмотря на характерный спад интенсивности
теплообмена вблизи входной области, локальные и средние значения числа Нуссельта значительно превышают значения, соответствующие круглой трубе и плоской щели. Например, стабилизированное число Нуссельта в канале с трехзаходной шнековой вставкой при Яе^ = 500 на 154,1% больше соответствующего значения в круглой трубе и на 23,2%
больше соответствующего значения в плоской щели.
Нц, _
(N11,,)!
51_,_._,_,_,_Jn
0 10 20 30 Рис.5. Распределение локального (сплошная линия) и среднего (пунктирная линия) чисел Нуссельта по длине канала
На рис 6 показаны профили безразмерной нормализованной температуры для различных чисел Рейнольдса. С ростом числа Рейнольдса область с максимальными значениями нормализованной температуры смещается к ребру шнека по ходу закрутки потока, в отличие от малых чисел Рейнольдса, где данное распределение более равномерно распределено по сечению канала. Здесь кривая 1 соответствует распределению нормализованной температуры в области, наиболее близко расположенной к ребру шнека
по ходу закрутки (ц2 = - (4я/15)).
ЯеА = 100 ЯеА = 500
Рис.6. Профили безразмерной нормализованной температуры: обозначения согласно рис.1
Заключение
Разработана математическая модель теплопереноса при ламинарном течении вязкой жидкости в канале с трехзаходной шнековой вставкой, для записи которой использовалась винтовая система координат. Результаты численных исследований показали, что поперечные потоки существенно влияют на процессы теплопереноса в исследуемом канале. Стабилизированное число Нуссельта в канале с трехзаходной шнековой вставкой при Яе^ = 500 на 154,1% больше соответствующего значения в круглой трубе и на 23,2% больше соответствующего значения в плоской щели. Таким образом, использование трехзаходной шнековой вставки в каналах теплообменного оборудования, использующего в качестве рабочих сред вязкие жидкости, приводит к значительному увеличению интенсивности теплообмена.
Литература
1. Назмеев Ю.Г. Гидродинамика и теплообмен закрученных потоков реологически сложных жидкостей. М.:Энергоатомиздат, 1996. 304с. http://opac.mpei.ru/ notices/index/ IdNotice:1903 5/Source:default#.
2. Vachagina E.K., Kadyirov A.I. The use of helical coordinate systems // Quarterly Journal of Mechanics and Applied Mathematics, 2014. Vol. 67, № 4. P.553-566. DOI: 10.1093/qjmam/hbu017 https://elibrary.ru/item. asp?id=24064796.
3. Назмеев Ю.Г., Конахина И.А. Исследование процесса интенсификации теплообмена при ламинарном течении вязких жидкостей в трубах с винтовой накаткой // Теплоэнергетика. 1993. №11. С.59-62.
4. Кадыйров А.И., Халитова Г.Р., Вачагина Е.К. Гидродинамическая структура потока вязкой жидкости в каналах с интенсификаторами в виде шнековой вставки // Труды Академэнерго. 2016. №3. С. 7-16. http://elibrary.ru/item.asp?id=26743777.
5. Bergles A.E., Jensen M.K., Somerscales E.F.C. and Manglik R.M. Literature Review of Heat Transfer Enhancement Technology for Heat Exchanges in Gas-Fired Applications, Report GRI 91-0146, Gas Research Institute, Chicago. 1991.
6. Митрофанова О.В. Гидродинамика и теплообмен закрученных потоков в каналах с завихрителями (Аналитический обзор) // ТВТ. 2003. Т. 41, №4. С.587-633. https://elibrary.ru/item.asp?id=17294932.
7. Kadyirov A.I., Abaydullin B.R., Vachagina E.K. Hydrodynamic characteristics of viscous fluid flow in screw channels formed by two ribs // IOP Conf. Series: Journal of Physics: Conf. Series 980 (2018) 012011. DOI: 10.1088/1742-6596/980/1/012011.
8. Mulugeta T. Sisay, Shimelis A. Emire, Hosahalli S. Ramaswamy, Tilahun S. Workneh Residence time distribution and flow pattern of reduced-gluten wheat-based formulations in a twine-screw extruder // LWT - Food Science and Technology. 2017. №79. P. 213-222. DOI: 10.1016/j.lwt.2017.01.015.
9. Shyy Woei Chang, Wei Ling Cai, Ruo Sin Syu. Heat transfer and pressure drop measurements for tubes fitted with twin and four twisted fins on rod // Experimental Thermal and Fluid Science. 2016. № 74. P. 220-234. DOI: 10.1016/j.expthermflusci.2016.01.001.
10. Исаченко В.П., Осипова В.А., Сукомел А.С. Теплопередача: Учебник для вузов. Изд. 3-е, перераб. и доп. М.: «Энергия», 1975. 488 с.
Авторы публикации
Конахина Ирина Александровна - д-р техн. наук, заместитель директора комплексного центра обучения в сфере энергоэффективности Казанского (Приволжского) федерального университета. E-mail: [email protected]
Хуснутдинова Эльвира Мусавировна - аспирант Казанского (Приволжского) федерального университета. E-mail: [email protected].
Кадыйров Айдар Ильдусович - канд. техн. наук, старший научный сотрудник лаборатории «Теплофизика и волновые технологии» Казанского научного центра Российской академии наук. E-mail: [email protected].
Вачагина Екатерина Константиновна - д-р техн. наук, ведущий научный сотрудник лаборатории «Теплофизика и волновые технологии» Казанского научного центра Российской академии наук. E-mail: [email protected].
References
1. Nazmeev Yu.G. Hydrodynamics and heat exchange of swirled flows of Theologically complex media. M. Energoatomizdat, 1996. 304 p. http://opac.mpei.rU/notices/index/IdNotice:19035/Source:default#.
2. Vachagina E.K., Kadyirov A.I. The use of helical coordinate systems // Quarterly Journal of Mechanics and Applied Mathematics, 2014. Vol. 67, № 4. P. 553-566. DOI: 10.1093/qjmam/hbu017 https://elibrary.ru/item. asp?id=24064796.
3. Nazmeev Yu.G., Konakhina I.A. Investigation of heat transfer enhancement for laminar viscous flow in pipes with helical rolling // Teploenegretika. 1993. №11. P.59-62.
4. Kadyirov A.I., Khalitova G.R., Vachagina E.K. Hydrodynamic structures of viscous flow in a coaxial channel with screw inserts // Transactions of Academenergo. 2016. №3. P. 7-16. http://elibrary.ru/item.asp?id=26743777.
5. Bergles A.E., Jensen M.K., Somerscales E.F.C. and Manglik R.M. Literature Review of Heat Transfer Enhancement Technology for Heat Exchanges in Gas-Fired Applications, Report GRI 91-0146, Gas Research Institute, Chicago. 1991.
6. Mitrofanova O.V. Hydrodynamics and heat transfer in swirling flows in channels with swirlers (analytical review) // High temperature. 2003. Vol. 41, №4. P. 518-559. https: //elibrary.ru/item. asp?id= 13435839.
7. Kadyirov A.I., Abaydullin B.R., Vachagina E.K. Hydrodynamic characteristics of viscous fluid flow in screw channels formed by two ribs // IOP Conf. Series: Journal of Physics: Conf. Series 980 (2018) 012011. DOI: 10.1088/1742-6596/980/1/012011.
8. Mulugeta T. Sisay, Shimelis A. Emire, Hosahalli S. Ramaswamy, Tilahun S. Workneh Residence time distribution and flow pattern of reduced-gluten wheat-based formulations in a twine-screw extruder // LWT - Food Science and Technology. 2017. №79. P. 213-222. DOI: 10.1016/j.lwt.2017.01.015.
9. Shyy Woei Chang, Wei Ling Cai, Ruo Sin Syu. Heat transfer and pressure drop measurements for tubes fitted with twin and four twisted fins on rod // Experimental Thermal and Fluid Science. 2016. № 74. P. 220-234. DOI: 10.1016/j.expthermflusci.2016.01.001.
10. Isachenko V.P., Osipova V.A., Sukomel A.S. Heat transfer. Handbook for high schools, 3rd edition, publisher. M., «Enegriya», 1975. P. 488.
Authors of the publication
Irina A Konakhina - Grand PhD in Engineering sciences Kazan Federal University, Kazan, Russia. E-mail: [email protected].
Elvira M. Khusnutdinova - PhD student Kazan Federal University, Kazan, Russia. E-mail: [email protected].
Aidar I. Kadyirov - PhD in Engineering sciences Kazan Scientific Center of RAS, Kazan, Russia. E-mail: [email protected].
Ekaterina K. Vachagina - Grand PhD in Engineering sciences, leading researcher Kazan Scientific Center of RAS, Kazan, Russia. E-mail: [email protected].
Поступила в редакцию 11 октября 2018 г.
