Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,686

ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ТЕПЛООТДАЧИ В КАНАЛАХ С ЭЛЕМЕНТАМИ ИНТЕНСИФИКАЦИИ

Лаптев А.Г. 1 Башаров М.М. 2 Рунов Д.М. 3
1 ФГБОУ ВПО «Казанский Государственный энергетический университет»
2 ОАО «ТАНЕКО»
3 Татарский научно-исследовательский и проектный институт нефти ОАО «Татнефть»
На основе использования модели турбулентного пограничного слоя, с функцией турбулентной вязкости Дайслера, получено выражение для расчета среднего коэффициента теплоотдачи при обтекании плоской поверхности. Показано согласование с экспериментальными данными для пластины и трубы. Выполнена корректировка параметров полученного уравнения для расчета коэффициентов теплоотдачи от поверхностей с элементами интенсификации (шероховатость, закрутка потока и т.п.). Корректировка заключается в введении относительных коэффициентов сопротивления возмущенных и невозмущенных потоков в степени 0,5 и динамической скорости. Получены выражения для чисел Нуссельта и Стантона для каналов с шероховатой поверхностью и закруткой потока. Показано согласование с известными экспериментальными данными различных авторов.
пограничный слой
коэффициент теплоотдачи
интенсификация теплообмена
закрученный поток
шероховатость
гидравлическое сопротивление
1. Дзюбенко Б.В., Кузма-Кичта Ю.А., Кутепов А.М. и др. Интенсификация тепло- и массообмена в энергетике. – М.: ФГУП «ЦНИИАТОМ-ИНФОРМ», 2003. – 232 с.
2. Дзюбенко Б.В., Мякочин А.С., Щербакова Н.У. Тепломассообмен в каналах с закруткой потока // Сборник научных статей «Современная наука». – 2011. – № 2(7). – С. 17–22.
3. Дьяконов С.Г., Елизаров В.И., Лаптев А.Г. Теоретические основы и моделирование процессов разделения веществ. – Казань: Изд-во Казанского университета, 1993. – 438 с.
4. Жукаускас А.А. Конвективный перенос в теплообменниках. – М.: Наука, 1982. – 472 с.
5. Кутателадзе С.С. Анализ подобия в теплофизике. – Новосибирск: Наука, 1982. – 280 с.
6. Кутутеладзе С.С., Леонтьев А.И. Тепломассообмен и трение в турбулентном пограничном слое. 2-е изд., перераб. и доп. – М.: Энергоатомиздат, 1985. – 320 с.
7. Кутателадзе С.С. Теплопередача и гидромеханическое сопротивление: Справочное пособие. – М.: Энергоиздат, 1990. – 367 с.
8. Лаптев А.Г. Модели пограничного слоя и расчет тепломассообменных процессов. – Казань: Изд-во Казан. ун-та, 2007. – 500 с.
9. Лаптев А.Г., Николаев Н.А., Башаров М.М. Методы интенсификации и моделирования тепломассообменных процессов: учебно-справочное пособие. – М.: «Теплотехник», 2011. – 335 с.
10. Разаев А.М., Филатов Л.Л. Теплопередача и гидравлическое сопротивление воды в трубах со спиральными канавками // Теплоэнергетика, 1986. – № 1. – С. 44–46.
11. Шлихтинг Г. Теория пограничного слоя / Пер. Вольперта Г.А. – М.: Главная редакция физико-математической литературы издательства «Наука», 1974. – 711 с.

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

Приближенное математическое описание процессов переноса в пограничном слое связано с моделями Прандтля, Кармана, Ландау и Левича, и др., а также с развитием гидродинамической аналогии Рейнольдсом и Чилтоном – Кольборном. Причем наиболее теоретически обоснованной и перспективной является модель диффузионного пограничного слоя Ландау – Левича [3 ,8, 9]. Известно, что турбулентный пограничный слой, как и всякая устойчивая статистическая система, имеет некоторые консервативные свойства [5, 6, 11]. На важную особенность пристенной турбулентности – весьма слабую зависимость некоторых характеристик осредненного течения по отношению к внешним возмущениям – особое внимание обратил С.С. Кутателадзе совместно с А.И. Леонтьевым. На основе предельных относительных законов теплообмена и трения были созданы расчетные методы [5–7].

Теоретическая основа рассмотренного ниже подхода заключается в использовании известных свойств консервативности законов трения к продольному градиенту давления в пограничном слое, т.е. структура математического описания элементарных актов переноса инвариантна к различным возмущениям и масштабу аппарата. Влияние этих факторов не изменяет структуру математического описания пограничного слоя, а учитывается параметрически [3, 8, 9].

Определение коэффициентов теплоотдачи

Допущения к математической модели теплоотдачи:

– входной участок значительно меньше длины обтекаемого тела (канала) и поэтому им можно пренебречь. При небольшой длине канала влияние входного участка можно учесть известными поправками;

– теплофизические свойства среды поперек пограничного слоя изменяются незначительно. При больших градиентах температур влияние можно учесть известными поправками;

– турбулентное число Прандтля близко к единице;

– основное сопротивление процессу теплоотдачи сосредоточено в пограничном слое, где происходит молекулярный перенос тепла в сочетании с затухающей турбулентной диффузией (модель Ландау-Левича).

Сопротивление переносу тепла в турбулентном пограничном слое записано в виде [8, 9]:

lap01.wmf (1)

где α – средний коэффициент теплоотдачи, Вт/ (м2·К); q* – относительная плотность теплового потока; ρ – плотность среды, кг/м3; cp – теплоемкость среды, Дж/(кг К); a, aт – коэффициенты молекулярной и турбулентной температуропроводности, м2/с; δ – толщина пограничного слоя, м; y – поперечная координата, м.

Коэффициент турбулентного переноса aт = vт/Prт, принят в форме функции Дайслера [11]

lap02.wmf, (2)

где lap03.wmf lap04.wmf lap05.wmf u* – динамическая скорость, м/с; v – коэффициент кинематической вязкости, м2/с, /Prт = vт/aт – турбулентное число Прандтля; vт – коэффициент турбулентной вязкости, м2/с, (Prт ≈ 1).

После интегрирования (1) с функцией (2) получена формула для коэффициента теплоотдачи при осевом движении потока в канале в виде:

lap07.wmf, (3)

где безразмерные величины связаны с областью интегрирования и установлены в виде:

k = 13,91; b = 2,5; R’ = 30; с = 0,14; Pr – число Прандтля.

Показатель степени при числе Прандтля Pr0,66 следует из закона затухания турбулентных пульсаций в пограничном слое.

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

Теплоотдача от пластин в трубе

Первоначально выполним проверку уравнения (3) для пластины и трубы.

При движении среды в турбулентном режиме вдоль плоской поверхности параметры уравнения (3) имеют вид [7, 11]:

динамическая скорость:

lap09.wmf (4)

коэффициент трения:

lap10.wmf (5)

число Рейнольдса:

lap11.wmf (6)

где u∞ – средняя скорость среды, м/с; L – длина пластины, м.

Значение средней безразмерной толщины пограничного слоя Rδ можно определить по следующим формулам [11], [8]:

Rδ = exp[0,4(u∞/u*–5,0)], (7)

или lap14.wmf, (8)

lap15.wmf (9)

Выражение (7) следует из логарифмического профиля скорости при u = u∞, у = δ.

На основе (3) запишем число Нуссельта NuL для пластины:

lap16.wmf (10)

Для проверки адекватности полученных значений чисел Нуссельта произведем их сравнение со значениями чисел Нуссельта по следующему известному уравнению:

lap17.wmf (11)

Результаты расчета чисел Нуссельта по полученным значениям Rδ, которые были рассчитаны по (7) и (9) представлены в табл. 1 и табл. 2 соответственно. Длина пластины 1 м.

Таблица 1

Значения чисел Нуссельта в зависимости от числа ReL

Re

Nu по (10)

Nu по (11)

Расхождение, %

105

347,8

370

6

106

1957,4

2334,5

16,2

107

14791,5

14730

0,4

Таблица 2

Значения чисел Нуссельта в зависимости от числа ReL

Re

Nu по (10)

Nu по (11)

Расхождение, %

105

341,8

370

7,6

106

1970,2

2334,5

15,6

107

15138,2

14730

2,7

На основании результатов расчета чисел Нуссельта по (10) и (11) можно сделать вывод, что их значения имеют удовлетворительную сходимость от 0,4 % до 16,2 % при значениях средней безразмерной толщины пограничного слоя Rδ рассчитанной по (7) и от 2,7 % до 15,6 % при Rδ рассчитанной по (9).

Для круглой трубы с гладкими стенками при осесиметричном движении среды динамическая скорость и Rδ в уравнении (3) имеют вид:

lap18.wmf (12)

lap19.wmf, (13)

где ξ – коэффициент гидравлического сопротивления; ucp – средняя скорость, м/с.

Значение средней безразмерной толщины пограничного слоя Rδ можно вычислить по формуле (7) при u∞ = ucp. В качестве примера определим Rδ используя безразмерный профиль скорости в виде:

lap20.wmf, (14)

где c = 8,74, n = 7.

На границе вязкого подслоя при у = δ1 функция, (14) имеет значение:

lap21.wmf. (15)

Отсюда запишем:

lap22.wmf. (16)

При у = δ из (15) имеем:

lap23.wmf (17)

Тогда

lap24.wmf (18)

В результате из (15) и (17) получим:

lap25.wmf (19)

Для пластины lap26.wmf и тогда:

lap28.wmf (20)

или lap29.wmf (21)

При R1 =11,6 получим (n = 7)

lap30.wmf (22)

Аналогично получим для трубы:

lap31.wmf. (23)

Расчеты по формулам (13) и (23) дают близкие значения Rδ. Уравнение (3) для трубы в безразмерной форме запишется в виде

lap32.wmf. (24)

Для сравнения результатов расчетов числа Нуссельта (Nud = αd/λ) по (24) использовалась известная формула для трубы:

lap33.wmf, (25)

а также уравнение Петухова:

lap34.wmf, (26)

где k1 = 1+ 3,4ξ; k2 = 11,7 + 1,8Pr–1/3, ξ – по формуле Блазиуса:

lap35.wmf, (104 < Re < 105). (27)

Значения чисел Нуссельта представлены в табл. 3.

Таблица 3

Значения чисел Нуссельта от чисел Рейнольдса (Pr = 5)

Re

Nu по (24)

Nu по (25)

Nu по (26)

Расхождение, %

104

72,5

72,8

74,6

2,8

2·104

128,4

126,8

132,7

4,45

4·104

227,9

220,8

235,7

6,3

5·104

274,1

263,9

283,4

6,9

6·104

319

305,4

329,4

7,3

8·104

405

384,4

417,4

7,9

В выражении (24) можно ввести известный множитель (Pr/Prст)0,25, который учитывает зависимость физических свойств жидкости от температуры и влияние направления теплового потока.

Теплоотдача для поверхностей с элементами интенсификации

Следует отметить, что выражения для чисел Нуссельта, аналогичные (10) и (24) были получены различными авторами еще 1950–70 гг., поэтому определение коэффициентов теплоотдачи для пластины и трубы с гладкой поверхностью не имеет особой новизны. Значительно более сложной задачей является определение теоретическим путем коэффициентов теплоотдачи для поверхностей с элементами интенсификации (шероховатость, выступы, кольцевые накатки, закрутка и т.д.). В настоящее время для этого используются в основном различные полуэмпирические подходы.

Рассмотрим применение выражения (3) для закрученного потока при стационарном режиме.

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

lap37.wmf, (28)

где ucp – средняя скорость в канале, м/с; θ – угол закрутки потока; ξ3 – коэффициент сопротивления потока с закруткой.

В пограничном слое с возмущениями (интенсификаций) параметры уравнения (3) имеют вид [8, 9]:

lap38.wmf lap39.wmf

lap40.wmf. (29)

где ξ – коэффициент сопротивления для осесимметричного потока, вычисляется по известному выражению, например (27). Запишем выражение (3) в более удобном для расчетов виде. Используя значения (28), (29), получим число Нуссельта:

lap41.wmf (30)

Аналогично запишем число Стантона lap42.wmf для канала с закруткой:

lap43.wmf. (31)

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

lap44.wmf,

где ξш, сfш – коэффициенты гидравлического сопротивления и трения с учетом шероховатости. Известно, что в режиме максимального проявления шероховатости наступает автомодельный режим и ξш ≈ 0,08.

Число Стантона для канала с шероховатой стенкой запишется в виде:

lap45.wmf. (32)

Результаты расчетов

Для расчета коэффициента гидравлического сопротивления закрученных потоков в работах [1, 2, 7, 11] представлены различные выражения и графики. На рис. 1 даны результаты расчета числа Нуссельта (30) и сравнение с обобщенными опытными данными для закрученного потока [11] Расхождение в пределах 10 %. Также удовлетворительное согласование (± 15 %) получено с экспериментальными данными, приведенными в работе [2].

lapt1a.wmf

а)

lapt1b.tif

б)

Рис. 1. а) зависимость комплекса Nu/Pr0,43 от числа Re в канале с ленточным завихрителем: 1 – расчет по уравнению (30); 2 – экспериментальные данные; 3 – для осевого потока; б) результаты расчета и опытные данные по теплоотдаче в канале с шероховатыми стенками 1,3 – расчет по уравнению (32), 2,4 –экспериментальные данные [4]

lapt2.tif

Рис. 2. Эффективность интенсификации теплообмена в круглом канале [2], [9]: 1 – область, занимаемая экспериментальными точками, 2 – при Nu/Nu0 = ξ/ξ0, 3 – расчет по формуле (30)

На рис. 2 представлена зависимость Nu/Nu0 от ξ/ξ0 для различных способов интенсификации теплообмена в каналах и результаты расчета с использованием формулы (30), где Nu0, ξ0 для каналов без интенсификации. Из рис. 2 следует, что рекомендуемое многими авторами отношение Nu/Nu0 = ξ/ξ0 справедливо до значения ξ/ξ0 ≤ (1,5÷2). При ξ/ξ0 > 2 начинается опережающий рост гидравлического сопротивления по сравнению с теплоотдачей.

Выводы

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

Полученные выражения для чисел Нуссельта и Стантона рекомендуется для практического применения при расчетах теплообменных аппаратов.

Работа выполнена в рамках проектной части государственного задания в сфере научной деятельности (Задание № 13.405.2014/К).


Библиографическая ссылка

Лаптев А.Г., Башаров М.М., Рунов Д.М. ОПРЕДЕЛЕНИЕ КОЭФФИЦИЕНТОВ ТЕПЛООТДАЧИ В КАНАЛАХ С ЭЛЕМЕНТАМИ ИНТЕНСИФИКАЦИИ // Международный журнал прикладных и фундаментальных исследований. – 2015. – № 3-3. – С. 311-315;
URL: https://applied-research.ru/ru/article/view?id=6534 (дата обращения: 20.06.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.252