А. Г. Лаптев, Т. М. Фарахов
МОДЕЛИ ТРЕНИЯ В ТУРБУЛЕНТНЫХ ПОТОКАХ
ПРИ ОБТЕКАНИИ РАЗЛИЧНЫХ ПОВЕРХНОСТЕЙ
Ключевые слова: касательное напряжение, перенос импульса, диссипация энергии, динамическая скорость, пограничный
слой.
Рассмотрены процессы переноса импульса в турбулентном пограничном слое. Получены оригинальные модели для расчета касательного напряжения, коэффициентов трения, толщины пограничного слоя и коэффициентов переноса импульса в пограничных слоях при различных условиях движения сред. Показаны примеры расчета данных характеристик и согласование с известными экспериментальными данными. Полученные уравнения и методы определения характеристик пограничного слоя могут использоваться при расчетах промышленных тепло- и массообменных аппаратов различных конструкций и другого оборудования.
Keywords: shear stress, momentum transfer, dissipation of energy, dynamic velocity, boundary layer.
The processes of transfer of momentum are considered in a turbulent boundary layer. Obtained the original model for the calculation of the shear stress, the coefficient offriction, the thickness of the boundary layer and the transport coefficients of the impulse (impulsootdachi) in the boundary layers at different conditions offluid flow. Shows the examples of calculation of these characteristics and agreement with the experimental data. The resulting equations and methods for determination of the boundary layer can be used in the calculation of industrial heat and mass transfer devices of various designs and other equipment
Введение
Теоретические методы моделирования и исследования массо- и теплообменных процессов условно подразделяются на точные, асимптотические, численные и приближенные. В связи с разнообразием конструкций контактных устройств и одновременно происходящих процессов обмена импульсом, массой и теплотой в большинстве задач химической технологии получить точные аналитические решения невозможно, поэтому наибольшее применение получили последние три метода. Так, например, среди различных асимптотических методов применяется метод функциональных параметров. Для этого строится разложение оператора относительно малой шкалы сравнения. Зависимость членов асимптотической последовательности от малого параметра осуществляется с помощью процедуры сращивания. Получаемые асимптотические ряды часто расходятся или очень медленно сходятся. Кроме этого удается вычислить только несколько первых членов разложения. Эти обстоятельства ограничивают использование асимптотических формул для инженерных расчетов.
Для моделирования и исследования процессов тепломассообмена в химической технологии используются чаще приближенные и численные методы. К приближенным методам относятся, например, однопараметрические интегральные методы в теории пограничного слоя, пленочная и пенетра-ционная модели, методы линеаризации уравнений и др. Приближенные методы позволяют получать необходимые формулы для выполнения конкретных инженерных расчетов.
В рамках приближенных методов находит применение подход, когда сложное явление заменяют совокупностью «элементарных процессов (актов)». Такими элементарными актами, прежде всего,
являются процессы переноса импульса, массы и тепла в пограничном слое.
Приближенное математическое описание процессов переноса в пограничном слое в первую очередь связано с моделями Прандтля и Кармана, гидродинамической аналогией Рейнольдса и Чилто-на-Кальборна и далее моделью диффузионного пограничного слоя Ландау-Левича. Кроме этого важное значение имеют консервативные свойства пограничного слоя - т.е. весьма слабые зависимости некоторых характеристик осредненного течения по отношению к внешним возмущениям, на которые одними из первых обратили внимание С.С. Кутате-ладзе и А.И. Леонтьев.
В данной статье рассматриваются различные модели пограничного слоя для вычисления напряжения трения и других характеристик на поверхности различных обтекаемых тел при турбулентном режиме. Касательное напряжение трения очень важно при определении коэффициентов гидравлического сопротивления обтекаемых тел, а также при решении тепломассообменных задач с использованием различных вариантов гидродинамической аналогии.
Рассматриваются пограничные слои при движении несжимаемой вязкой жидкости с постоянными физическими свойствами. В статье широко используется понятие коэффициента переноса импульса, которое позволяет получить в аналитическом виде расчетные формулы для коэффициента трения и толщины пограничного слоя. Даны результаты расчета этих характеристик и сравнение с экспериментальными данными.
Определение касательного напряжения трения
Для расчета коэффициентов переноса в тепло- и массообменных аппаратах различных конструкций необходимо определить среднее значение
касательного напряжения на стенке тст (или динамической скорости и* = Я ст / р, ), где р -
плотность среды, кг/м3.
Рассмотрим общий метод определения касательного напряжения трения (локального потока импульса) на стенке - поверхности обтекаемого тела. Пусть ось Оу направлена по нормали к стенке. В этом случае касательное напряжение в любой точке пограничного слоя будет определяться по известному выражению с учетом молекулярного и турбулентного переноса
йи . . йи . ,йи
Т = -Р~— Мт (У)— = -рУ + Ут
ау ау ау
(1)
где М,Мт - динамическая молекулярная и турбулентная вязкость, Па-с;У,Ут - коэффициенты кинематической молекулярной и турбулентной вязкости, м2/с; йи /йу - градиент скорости; у - поперечная координата, м.
Касательное напряжение на стенке, то есть при у = 0 , может быть выражено, используя уравнение импульсоотдачи (потока импульса) в виде:
т = -ру(иж - ист ) = руи^,
(2)
где у - коэффициент переноса импульса, м/с;
,ист - скорость в ядре потока и на стенке, м/с. В уравнении (2) учитывая условие прилипания, имеем ^ст = 0.
В литературе [1] у имеет размерность кг/(м2-с). Однако, удобнее в выражении (2) использовать размерность у - м/с.
Еще в середине XIX столетия изучением процесса трения в трубах занимались сначала Дарси, а потом Фаннинг, которые и ввели понятие о коэффициенте трения
= т = 2у, (3)
ры\2 и»
при иж = иср , где иср - средняя скорость потока, м/с.
Запишем значение динамической скорости, используя выражение (2)
и* = ■\[уиж • (4)
Рассмотрим подход определения коэффициента переноса импульса.
Введем относительный поток импульса
т* (у) = т(у)/тст . Тогда из выражения (1) запишем
“ПО
’1йи
• т * (у)ау
Т * •’V + V
со о о
(5)
где 5 - толщина гидродинамического пограничного слоя, м.
Сравнивая (5) и (2) приходим к следующему выражению для определения коэффициента переноса импульса [3]
у =
г Т*(у)йу
У + Ут (у)
-1
(6)
Выражение (6) является достаточно общим и позволяет рассчитать коэффициент переноса импульса на основе известного коэффициента молекулярной вязкости, а также известного характера изменения коэффициента турбулентной вязкости
Ут (у) и относительного потока импульса т (у) в пограничном слое. Применение выражения (6) рассмотрено ниже.
В тех случаях, когда значения С ^ неизвестны (например, в аппаратах с вводом внешней энергии в контактирующие фазы), находит применение подход вычисления динамической скорости (касательного напряжения), с использованием средней диссипируемой энергии в единице объема среды [2, 5-7].
Скорость диссипации энергии £ в жидкости, обусловленной вязкостью, при градиенте скорости йи / йу записывают в виде [8]
ёе &и
£ = — = т------, (7)
ёу ёу
где е - плотность потока энергии, е = ти; т - касательное напряжение, Па; и - скорость среды, м/с; у - поперечная координата, м.
Получим значение динамической скорости на основе уравнения (7) и известного характера переноса импульса в турбулентном пограничном слое.
В уравнении (7) значение производной найдем из выражения:
,йи йи т
т = -р(У + Ут )—, — = —---------------, (8)
йу йу р(у + ут)
Из (7) и (8) получим: .2
Т
Б = -
Ут + ут (у).
(9)
р(у + Ут )
Среднее значение скорости диссипируемой энергии в пограничном слое толщиной 8 запишем в виде т = тст :
2 8
Б =
Тст Г Гр Г і
ёу
5р 0у + ут (у)
(10)
где подынтегральная функция определяет сопротивление переносу импульса в турбулентном пограничном слое (6), и тогда уравнение (10) получит форму:
_ Т2
Б= Тст
4
и* р
(11)
8ру Гу Отсюда найдем значение динамической скорости:
и
u * =
\0
єбу
P
8—У
P
(12)
V ’ /
где Яд = и*8/у — безразмерная толщина пограничного слоя.
По данному выражению можно вычислять
и* (или С ^ ) для аппаратов с подводом внешней
энергии или при обтекании тел сложной геометрии, например, насадочных или зернистых слоев, пучков труб и т. п.
Коэффициенты переноса импульса
В общем случае в выражении (12) коэффициент переноса импульса определим путем интегрирования зависимости (6):
5 г* (y)dy
У 0К+ — (у)
(13)
Для наглядности получим формулу для у на основе применения относительно простой двухслойной модели пограничного слоя Прандтля. Запишем интеграл (13) в виде сопротивлений переносу импульса в вязком подслое (т = 0) и турбулентной области (т = 2и*у)
51 dy + 5 dy
У =1 - +1
О
51
Хп* y
получим
У =
п*
Ri +—ln (5/Ri) X
(14)
(15)
где Я§ = и*Г/у- безразмерная толщина турбулентного пограничного слоя;
Я1 = и*Г1 /у = 11,6 - безразмерная толщина вязкого подслоя на пластине; % = 0,4- константа Прандтля.
Далее с применением трехслойной модели Кармана выражение (13) получит вид:
Г
і 51 52 , 5 ,
I =f І. + 5 f А + 2,5 f А,
у О — 5u*y 5u*y
(16)
где 81 = 5у/и* ;82 = 30у/и*.
Первое слагаемое в выражении (16) определяет сопротивление переносу импульса в вязком
подслое толщиной 81, второе - в переходной (буферной) области толщиной 82 — 81, а третье - в
турбулентной области толщиной 8 — 82.
Для удобства далее используем безразмер-
ные числа
~ u *5. ~ u *5 2
R = * 1 • R = * 2
где согласно модели Кармана Я1 = 5; Я2 = 30 . После интегрирования (16) получим
или
1 R1 R^ 2,5,
— = —L + —— ln 6 + ^ ln
У = ■
(17)
(18)
2,792Я1 + 2,51п(8 /Я2)
Как известно, в моделях Прандтля и Кармана в вязком подслое Ут = 0.
Согласно теории Ландау и Левича, подтвержденной Дайслером, турбулентность в вязком подслое описывается функцией [9]
4
при У < 6,
(19)
тогда после интегрирования аналогичного выраже-
ния
(16) с —г (y) (19) получим
У =
2,288R1 + 2,5ln(R5 / R2)
(2О)
где R1 = 6, для плоского пограничного слоя.
Как следует из расчета по выражениям (18) и (20) учет затухания турбулентных пульсаций в вязком подслое в модели Кармана практически слабо влияет на значение коэффициента переноса импульса. Наиболее значительное влияние наблюдается для двухфазных сред [3, 5, 6].
В работах [2,3] выполнено интегрирование выражения (6) с различными функциями турбулентной вязкости в пограничном слое. Например, с
функциями Ут (y) Owen P. получено:
_ u^ (21)
Г= 11,73 + 2,5 ln[(0,4u*S - 3v)/(5v) ^
Или с функциями Ут (у) Дайсслера и Хан-
ратги
у=-
u*
(22)
5,309 + 2,5іиЯГ
Результаты расчета у для пластины по (15), (18), (20), (21) и (22) согласуются с расхождением ± 7%.
Ниже показано, что использование данных выражений для коэффициента импульсоотдачи позволяют получить в явном виде ряд важных характеристик турбулентного пограничного слоя.
Толщина пограничного слоя
Локальная и средняя толщина турбулентного пограничного слоя на пластине имеет вид [9]: 0,37х
ReX’2
г 1 r0,37xdx 0,205L
5=L J-
(23)
LJ ReX’2
О ^сх
Re0’2
где L - длина пластины, м; х- продольная координата; числа Рейнольдса
Rex = x/v; Re[ = uюL/v.
0,2
У
u
u
u
u
С применением выражений для у, толщину пограничного слоя найдем из потокового соотношения (2):
2 _
тст = и *р = руи да , (24)
На основе (24) с учетом коэффициента переноса импульса (22) получено [3]:
’ (и У
0,4 — — 5,309
_ Vи * )_
Аналогичные выражения для 8 или R8 можно получить из (24) с коэффициентами у (15) и (18):
Яд = ехр
(25)
Яд = Я^ехр
Ґ
0,4
и г
- Я
Л
1
0,4
--Я,
(26)
(27)
где= 5 ^2 = 30 ^1 = 11,6 .
Выражения (25) - (27) дают практически одинаковые результаты и согласуются с результатами расчета по выражению (23).
Коэффициенты трения
Рассмотрим применение выражений для коэффициентов импульсоотдачи в форме коэффициентов трения.
Из (2) с у (15) получим динамическую скорость в виде:
и иг . (28)
и * —---------------
Я + -іп(Яг /Я1)
2
Выражение (28) связывает динамическую скорость трения и * с характеристиками турбулентного пограничного слоя в двухслойной модели Прандтля -Яд, Яь %. Расчет и* по (28) удовлетворительно согласуется с и* , полученному по выражению (3) и* = Д(С д).
С использованием выражения (28) можно записать коэффициент трения
-1-2
Я1 + !іп((г/Я1)
ж
(29)
Аналогично с выражением (18)
С = 2|2,792Я, + 2,5 1п(яг / Я
(30)
(31)
и с формулой (22)
С у = 2[5,309 + 2,51п Я8 ]—2.
На рис. 1,2 представлены результаты расчета по выражению (29) коэффициента трения на пластине и коэффициента сопротивления (X = 4С у) для труб. Расхождение с эксперимен-
тальными данными в пределах 10-15%. Аналогичные результаты дают и выражения (30) и (31).
|5 2 4 6 8-ю6 2 4 6 8 Ш7 2 4 6 8
Рис. 1 - Коэффициент сопротивления плоской пластины. Сплошная линия - расчет по (29); точки экспериментальные данные [12]
Рис. 2 - Коэффициент сопротивления для гладких круглых труб. Сплошная линия - расчет по (29) (Х = 4Су); точки экспериментальные данные [12]
При обтекании тел сложной геометрии можно воспользоваться выражением (12). Получим
с =-2-
С = и!
бЯ 5уу р
(32)
По данному выражению с Яд (26) и у (15) выполнены расчеты для различных случаев турбулентного обтекания поверхностей: шахматных и коридорных пучков труб; зернистых и насадочных слоев; в аппарате с механическими мешалками и в гидроциклоне. Средняя скорость диссипации энергии Б записывается через перепад давления (мощность N)
ХТ Л
(33)
Б=N=Ар£
V V
где Q - объемный расход, м3/с.
В результате многочисленных расчетов различного оборудования и обобщения полученных результатов получены выражения
С 2А
С Д =~Т
Г-Л0,5
бу
Р
и * = А
и X V Г J
бу
^-л0,25
\Р;
(34)
(35)
и
и
X
и
А
0,4
где коэффициент пропорциональности является слабо меняющийся функцией А = /(Яе)и с погрешностью ±12% можно принять: А = 2,0 — для шахматных и коридорных пучков труб; А = 1,8 — для насадочных и зернистых слоев; А = 2,0 — для аппаратов с мешалками;
А = 1,55 — для гидроциклонов.
Выражение (35) соответствуют известной зависимости [7,10,11], где коэффициент А находится экспериментально для каждого типа аппарата и установлен в пределах А = 1,7 + 2,4.
Основные результаты и выводы
Выполненные в данной работе аналитические исследования моделей пограничного слоя и сделанные сопоставления с экспериментальными данными и известными расчетными зависимостями позволяют сделать следующие обобщения и выводы:
1. Для турбулентного режима получено аналитическое выражение (12), где основными параметрами являются средняя скорость диссипации энергии, толщина пограничного слоя и коэффициент переноса импульса.
2. Путем использования моделей Прандтля, Кармана, Ландау, Левича и Дайсслера и интегрирования уравнения (6), получены выражения для коэффициентов переноса для турбулентного режима.
3. На основе баланса импульса в пограничном слое с использованием выражений для коэффициентов переноса импульса получены уравнения (25) - (27) для расчета средней толщины гидродинамического пограничного слоя.
4. Получено обобщенное выражение для коэффициента трения (32) и ряд частных случаев (34) при турбулентном обтекании различных тел.
Приведенные в данной работе уравнения и методы расчета гидромеханических характеристик пограничного слоя могут использоваться в расчетах промышленных тепло- и массообменных аппаратов и другого оборудования в зависимости от решаемых задач [4-7,13,14].
Литература
1. П. Бенедек, А. Ласло Научные основы химической технологии. Под ред. П.Г. Романкова и М.И Курочкиной. - Изд. «Химия», 1970.376 с.
2. А.Г. Лаптев Модели пограничного слоя и расчет теп-ломассобенных процессов. Казань: Изд-во Казанск. унта, 2007.500 с.
3. А.Г. Лаптев, Н.А. Николаев, М.М. Башаров Методы интенсификации и моделирования тепломассообменных процессов. М.: Теплотехник, 2011.287 с.
4. А.Г. Лаптев, Т.М. Фарахов Математические модели переноса импульса в пограничном слое. Инженернофизический журнал, 86, 3, (2013)
5. С.Г. Дьяконов, В.И. Елизаров, А.Г. Лаптев Теоретические основы и моделирование процессов разделения веществ. Казань.: Изд-во Казанск. ун-та, 1993.437 с.
6. Л.В. Брагинский, В.И. Бегачев, В.М. Барабаш Перемешивание в жидких средах: физические основы и инженерные методы расчета. Л.: Химия, 1984.389 с.
7. В.Н. Соколов, И.В. Доманский Газожидкостные реакторы. Л.: Машиностроение, 1976.216 с.
8. Л.Д. Ландау, В.М. Лившиц Теоретическая физика: учеб. пособие: в 10 т. Т.6. Гидродинамика. 5-е изд. М.: Физматлит, 2003.560 с.
9. Г. Шлихтинг Теория пограничного слоя. М.: Наука, 1974.714 с.
10. И.В.Доманский, В.Н.Соколов Обобщение различных случаев конвективного теплообмена с помощью полу-эмпирической теории турбулентного переноса. Теор. основы хим. технол., 2, 5, 761-767 (1968).
11. П.П. Ермаков Влияние воздействия акустических колебаний на процесс массопередачи .Теор. основы хим. технол., 25, 2, 198-203 (1991).
12. Г.Н. Абрамович Прикладная газовая динамика. Изд-е 4-е: Главная ред. физико-математической лит-ры изд-ва «Наука», М.,1976. 888с.
13. А.Г. Лаптев, М.М., Башаров, А.И. Фарахова Математическое моделирование разделения тонкодисперсной фазы в гидроциклоне. Вестник Казан. гос. технол. ун-та. 15, 24, 132-134 (2012)
14. А.Г. Лаптев, М.И. Фарахов Энергоэффективное оборудование разделения и очистки веществ в химической технологии. Вестник Казан. гос. технол. ун-та, 9, 152158, (2011)
© А. Г. Лаптев - д-р техн. наук, проф., зав. каф. КГЭУ, М_к£еи@таі1.га; Т. М. Фарахов - канд. техн. наук, вед. инж. ООО ИВЦ «Инжехим».