Динамика вязкого газа, турбулентность и струи
.pdf12.9. Универсальные законы распределения скоростей для больших чисел Рейнольдса 141
При увеличении n экспериментальные данные лучше согласуются с теоретическими для больших чисел Рейнольдса, хуже – для малых (см. рис. 12.4, кривая 5, которая получена для n = 10).
Из (12.65) выведем выражение для v* :
|
v = 0,15u |
7/8 |
ж |
|
ц1/8 |
|
|
|||||
|
|
|
|
|
|
|
. |
|
(12.66) |
|||
|
|
з |
|
ч |
|
|||||||
|
* |
|
|
|
|
|
|
|
||||
|
|
|
|
и y |
ш |
|
|
|
|
|||
Следовательно, касательное напряжение на стенке |
|
|
||||||||||
|
2 |
|
|
|
|
7/ 4 |
ж |
ц1/ 4 |
|
|||
0 = |
v* |
= 0, 0225 |
u |
|
|
|
з |
|
ч |
, |
||
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
и y |
ш |
|
|
или |
|
|
|
|
|
|
|
ц1/ 4 |
|
|
||
|
0 |
= 0, 0225 |
ж |
|
|
. |
(12.67) |
|||||
|
з |
|
|
|
ч |
|
|
|||||
|
V 2 |
|
|
|
|
|
||||||
|
|
|
иVR ш |
|
|
|
|
При увеличении n в показателе степенной зависимости распределения скоростей должны стремиться к логарифмической зависимости.
12.9.Универсальные законы распределения скоростей для больших чисел Рейнольдса
При увеличении Re показатель степени 1/n уменьшается. Поэтому реальным кажется предположение о том, что в этом случае должен быть справедлив логарифмический закон, как предельный для очень малой степени. С физической точки зрения такие законы характеризуются тем, что учитывается только турбулентное трение, ламинарное же становится несущественным.
Вернемся к рассмотрению логарифмического закона в виде (12.50):
Aln B .
Сравнение с опытами Никурадзе для гладких труб дает A = 2,5, B = 5,5. Логарифмический закон распределения скоростей применим также к шероховатым трубам. В зависимости от величины шероховатости мы получим разные значения B. Значение A остается неизменным.
142 |
12. ТУРБУЛЕНТНЫЕ ТЕЧЕНИЯ |
В непосредственной близости от стенки в ламинарном подслое выполняются равенства
0 = dudy , или 0 = uy = v*2 .
Из этих равенств получим
u |
= |
yv* |
, или = . |
(12.68) |
|
|
|||
v* |
|
На рис. 38 формулам (12.68) соответствует кривая 1. Измерения показывают, что:
при < 5 – ламинарное течение;
при 5 < < 70 – ламинарно-турбулентное течение; при > 70 – турбулентное течение.
12.10.Универсальный закон сопротивления для гладких труб при больших числах Рейнольдса
Выведем из логарифмического закона для распределения скоростей закон сопротивления. Пусть
= Aln ( ) + B , |
(12.69) |
где |
= |
u |
, = |
yv* |
. |
|
|
||||
|
|
v* |
Вычислим среднюю по сечению трубы скорость, приняв во внимание, что y в (12.69) отсчитывается от стенки, а в интеграл (12.58) входит радиус, т.е. r = R – y.
|
|
ж Aln |
Rv* |
+ B - |
3 |
Aц . |
|
|
u = v |
(12.70) |
|||||||
|
|
|||||||
|
|
* з |
|
ч |
|
|||
|
|
и |
2 ш |
|
12.11. Турбулентный пограничный слой на плоской пластине |
143 |
Из выражений 0 = v*2 и |
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
||
0 = |
|
|
u |
8 получим |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
|
= |
u |
|
|
|
, |
или |
u |
= |
8 |
. |
(12.71) |
||||||
* |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
8 |
|
|
|
|
|
v* |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
Преобразуем выражение Rv* / , входящее в (12.70), к следующему виду:
Rv* |
|
|
|
|
|
|
|
|
|
|
|
= |
Du |
|
|
|
. |
(11.72) |
|||||
|
|
|
|
|
|
|
|
||||
|
4 |
2 |
|||||||||
|
|
|
|
|
Подставив выражение (11.71) и (11.72) в (11.70), получим соотношение
|
|
|
|
|
|
|
|
|
|
|
8 |
|
= A ln (Re |
|
) - Aln (4 |
|
) + B - |
3 |
|
|
|
|
|
2 |
A , |
(11.73) |
||||||
|
2 |
|||||||||
|
|
|
|
|
|
|
|
|
которое можно записать как
1 |
|
= l ln (Re |
|
) - b . |
(11.74) |
|
|
|
|||||
|
|
|
||||
|
Выражение (11.74) известно как универсальный закон сопротивления Прандтля для гладких труб. Для A = 2,5 и B = 5,5 расчет дает l = 2,035, b = 0,91. Из экспериментов Никурадзе получены достаточно близкие к ним значения l = = 2,0, b = 0,8, справедливые до Re = 3,2 106. На рис. 2.3 закону сопротивления Прандтля соответствует кривая 3.
12.11.Турбулентный пограничный слой на плоской пластине
Расчет турбулентного пограничного слоя на плоской пластине проведем на основе приближенного интегрального метода. Сложность расчета заключается в том, что в интегральное уравнение импульсов входит напряжение трения на стенке, которое будет определяться течением в ламинарном подслое, аппроксимация же профиля скорости зависимостью (12.57) или (12.50) справедлива в развитом турбулентном течении. Необходимо смыкание ламинарного и турбулентного течений в некоторой правильно выбранной точке течения.
144 |
12. ТУРБУЛЕНТНЫЕ ТЕЧЕНИЯ |
V
л uл
Рис. 12.5
Рассмотрим течение, показанное на рис. 12.5. Обозначения на рисунке следующие: – толщина пограничного слоя; т – толщина ламинарного подслоя; V – скорость внешнего невозмущенного течения; uл – скорость на границе ламинарного подслоя.
Введем число Рейнольдса, рассчитанное по параметрам на границе ламинарного подслоя:
Re2 = |
u |
/ = u |
л |
/ . (12.75) |
л |
л л |
л |
|
|
В ламинарном подслое скорость линейно зависит от координаты |
|
|||
u / uл = y / л . |
|
|
|
(12.76) |
В остальной турбулентной части пограничного слоя будем считать справедливым степенной закон
|
u |
|
ж y ц1/ n |
|
|
|
|||
|
|
= з |
|
|
ч . |
|
|
(12.77) |
|
|
|
|
|
|
|
||||
V |
|
и |
ш |
|
|
|
|||
Напряжение трения на стенке |
|
|
|
|
|
|
|
|
|
w = |
ж du ц |
= |
uл |
. |
(12.78) |
||||
з |
|
ч |
|
||||||
|
|
и dy ш |
w |
л |
|
|
|||
|
|
|
|
|
|
|
|
|
Выразим uл / л через параметры невозмущенного течения. На границе ламинарного подслоя выполняются обе формулы (12.75) и (12.77). Тогда
|
uл |
ж |
л ц1/ n |
|
|
|
|
= з |
|
ч . |
(12.79) |
|
|
|
|||
|
V |
и |
|
ш |
|
Система двух уравнений (11.75) и (11.79) служит для определения |
л и uл через |
параметры невозмущенного течения, которые после подстановки в (12.78) дают
|
|
|
|
|
|
1-n |
|
|
|
|
|
Re2л |
|
|
|
|
|
w |
= |
1 ж |
ц1+n |
, |
(12.80) |
|||
V 2 |
|
з |
|
ч |
||||
|
Re и |
Re |
ш |
|
|
12.11. Турбулентный пограничный слой на плоской пластине |
145 |
где Re = V / . Из опытных данных следует, что переход ламинарного течения в турбулентное происходит при Reл = 12,5. Если принять n = 7, то формулу (12.80) можно преобразовать к выражению
w |
= 0, 0225 |
ж |
ц1/ 4 |
|
||
з |
|
ч |
, |
(12.81) |
||
V 2 |
|
|||||
|
иV |
ш |
|
|
что с точностью до четвертого знака после запятой совпадает с (12.67). Уравнение импульсов пригодно и для турбулентных течений. Для пло-
ской пластины оно будет иметь вид
( |
** )ў = |
w |
. |
(12.82) |
2 |
||||
|
|
V |
|
Для вычисления * и ** можно воспользоваться степенной аппроксимацией (12.77) во всем пограничном слое. Ошибка из-за пренебрежения ламинарным подслоем в интегралах не будет слишком велика. После вычислений получим
* = /8, |
** = 7 /72. |
(12.83) |
Подстановка (12.81) и (12.83) в (12.82) приводит к уравнению |
|
|
7 /72 = 0,0226 ( /V )1/4. |
(12.84) |
Будем считать, что турбулентное течение начинается непосредственно от передней кромки пластины ( = 0 при x = 0). Решение уравнения (12.84) с учетом принятого предложения дает
( |
|
|
) |
|
|
|
( |
|
) |
-1/5 , |
|
||
|
x |
|
|
= |
0,37x Vx |
|
|
||||||
* |
( |
|
|
|
) |
|
|
( |
|
) |
-1/5 , |
(12.85) |
|
|
|
x |
|
= 0, 046x Vx |
|
|
|||||||
** |
( |
|
|
) |
|
( |
|
|
) |
-1/5 . |
|
||
|
|
x |
|
= 0, 036x Vx |
|
|
|
||||||
Полное сопротивление пластины |
|
|
|
|
|
|
|
||||||
L |
|
|
V 2 |
** (L) |
= 0, 036 V 2bL (VL )-1/5 . |
|
|||||||
W = b т wdx = |
|
|
(12.86) |
||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
146 |
12. ТУРБУЛЕНТНЫЕ ТЕЧЕНИЯ |
Коэффициент сопротивления |
|
C f = 0,072(VL )-1/5 . |
(12.87) |
Опыт дает значение постоянной в (12.87), равное 0,074.
В том случае, если переход ламинарного течения в турбулентное происходит на некотором расстоянии xn от передней кромки, сопротивление пластины по Прандтлю
|
Wпл = W – |
W, |
(12.88) |
|||
где |
|
|
|
|
|
|
W = |
1 |
V 2bxn (C ft - C fl ), |
|
|||
|
|
|||||
2 |
|
|
|
|
|
|
C ft |
= 0, 074 / Re1/5n , |
|
||||
Cfl = 1,328 / Re1/n 2 , |
|
|||||
|
|
Ren = |
Vxn |
. |
|
|
|
|
|
|
|||
Коэффициент сопротивления в этом случае дают формулы |
||||||
C f = 0, 074 / Re1/5L |
- C1 / ReL |
, |
||||
C1 = Ren (C ft - C fl ). |
(12.89) |
|||||
|
В реальных течениях число Рейнольдса перехода может изменяться в широких пределах. Его значение зависит от уровня и спектрального состава турбулентности набегающего потока и состояния поверхности пластины.
Полученные формулы справедливы до тех пор, пока справедлив закон Блазиуса, т.е. для Re< 105, что дает ReL 107.
Задача 12.1
Вывести коэффициент сопротивления гладкой трубы для степенного распределения скорости. Показать, что при n = 1/7 получим закон Блазиуса.
|
12.11. Турбулентный пограничный слой на плоской пластине |
147 |
|
|
|
|
|
|
|
|
|
|
|
|
|
13. СВОБОДНАЯ ТУРБУЛЕНТНОСТЬ
Термин «свободная турбулентность» обычно используют для обозначения потоков, которые непосредственно не соприкасаются с твердыми границами. Турбулентность в таких течениях может распростра-
няться в поперечном направлении в окружающую среду. Наиболее характерные примеры таких потоков: свободная граница струи, свободная струя, след за телом. Такие течения обладают свойствами, характерными для пограничных слоев; протяженность течения в поперечном направлении мала по сравнению с продольным. Следовательно, для их расчета можно использовать уравнения пограничного слоя.
Рассмотрим плоскую двумерную задачу. Уравнения Прандтля для нее будут иметь следующий вид:
¶¶ut + u ¶¶ux + v ¶¶uy = 1 ¶¶y ,
(13.1)
¶¶ux + ¶¶yv = 0.
Здесь под мы будем подразумевать только турбулентное трение. В отсутствие твердых границ оно намного превосходит ламинарное. В уравнениях (13.1) опущен член ¶P¶x , так как в первом приближении в задачах о свободной
турбулентности давление можно считать постоянным.
Для того чтобы проинтегрировать уравнения (13.1), необходимо знать связь со скоростью течения. В данном разделе рассмотрим три предложенные Прандтлем формулы:
= |
l |
2 |
¶u |
¶u |
, |
(13.2) |
|
¶y |
¶y |
||||
|
|
|
|
|
148 |
13.СВОБОДНАЯ ТУРБУЛЕНТНОСТЬ |
= l2 ¶u
¶y
=¶¶uy =
ж |
¶u ц2 |
2 |
ж |
¶2u ц2 |
||
з |
ч |
+ l1 |
з |
¶y2 |
ч |
, |
и ¶y ш |
|
и |
ш |
|
1b (umax - umin ) ¶¶uy ,
(13.3)
(13.4)
здесь l и l1 представляют собой функции только точки, выбор которых должен производиться исходя из условий задачи; b означает ширину зоны перемешивания; 1 – эмпирическая постоянная; – коэффициент кажущейся кинематической вязкости, который зависит только от x.
13.1. Развитие во времени слоя раздела
Пусть в момент времени t = 0 имеется два потока со скоростями V1 и V2, показанных на рис. 13.1, а, которые приходят в соприкосновение. Вследствие турбулентного перемешивания течение будет превращаться в слой с непрерывным изменением скорости (рис. 40, б), ширина которого b с течением времени растет. Возникает плоскопараллельное неустановившееся слоистое течение u = u(y, t), v = 0.
Уравнение (13.1) преобразуется к виду
¶u |
= |
1 ¶ |
|
|
|
(13.5) |
|||
¶t |
|
|
|
¶y |
|
|
|
||
|
|
|
|
|
|
|
|
||
с граничными условиями |
|
|
|
|
|
|
|
|
|
u (b) = V1, u (-b) |
= V2, |
|
¶u |
(±b) = 0. |
(13.6) |
||||
|
|
||||||||
|
|
|
|
|
|
|
¶y |
|
Примем для напряжения сдвига формулу (13.2), тогда уравнение (13.5) можно записать в виде
¶u |
= l |
2 |
|
¶u |
|
¶2u |
. |
(13.7) |
|
|
|||||||
¶t |
|
|
¶y |
|
¶y2 |
|||
|
|
|
|
|
|
13.1. Развитие во времени слоя раздела |
149 |
Ширина зоны перемешивания растет со временем: b = b(t). Примем, что длина пути смешения пропорциональна ширине слоя, l = b. Будем искать автомо-
дельные решения такие, что u f( ), где
= y /b, а b t p.
Подстановка этих выражений в уравнение (13.7) приводит к выводу: p = 1.
Таким образом, ищем решение задачи при следующих допущениях:
|
|
|
V1 |
|
|
|
V1 |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y = 0 |
|
|
|
|
2b |
|||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
a |
b |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
V2 |
V2 |
|
|
|
|
|
|
|||||
|
|
|
|
|
Рис. 13.1 |
|
|
|
|
|
|
b = Bt, |
= y /b = y /Bt. |
(13.8) |
Представим решение задачи в виде |
|
|
u = vm + Af( ), |
(13.9) |
где vm = (V1 + V2) /2, A = (V1 – V2) /2.
Такое представление решения удобно потому, что граничные условия за-
дачи (13.6) будут иметь следующий вид: |
|
|
||||
f(1) = 1, f(–1) = –1, |
f ( 1) = 0. |
(13.10) |
||||
Для определения f( ) получим уравнение |
|
|
||||
f ў + |
|
2 A |
f ў f |
ўў = 0 . |
(13.11) |
|
|
|
|||||
|
B |
|||||
|
|
|
|
|
|
|
Проинтегрировав (13.11), при условии, что |
f ў № 0 , получим решение задачи |
|||||
f = C |
3 + C |
+ C , |
|
|||
0 |
|
1 |
2 |
|
||
C0 = - |
B |
. |
|
(13.12) |
||
|
|
|||||
|
|
|
||||
|
|
6 2 A |
|
|
Условие f ў = 0 дает тривиальное решение f = const, которое нас не устраивает. Из граничных условий (13.10) находим значение констант C0, C1 и C2:
C0 = –0,5, C1 = 1,5, C2 = 0.
150 |
|
|
|
|
|
|
|
13.СВОБОДНАЯ ТУРБУЛЕНТНОСТЬ |
|||||||||||
Теперь можно выписать решение задачи |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
1 |
|
|
1 |
|
й |
3 |
ж y ц |
|
1 |
ж y ц |
3 |
щ |
|
|
||||
u ( y, t) = |
(V1 |
+V2 ) + |
(V1 |
-V2 ) к |
- |
|
ъ |
, |
(13.13) |
||||||||||
|
|
|
з |
|
ч |
|
з |
|
ч |
|
|||||||||
|
|
|
|
|
|
|
|||||||||||||
|
2 |
|
|
2 |
|
к |
2 |
и |
|
ш |
|
2 |
и |
|
ш |
|
ъ |
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
|
|
|
|
ы |
|
|
причем = 32 2 (V1 -V2 )t .
В решение задачи входит эмпирическая постоянная = l/b, которую можно определить только из опыта. Скорость в зоне смещения переходит в скорость невозмущенного течения не асимптотически, в точках y = b произ-
водная ¶2u¶y2 терпит разрыв. Такой недостаток присущ всем решениям, по-
лученным на основе формулы Прандтля (13.2).
13.2. Плоский след
Рассмотрим стационарное течение на больших расстояниях x позади тела толщиной d и высотой h, где скорость спутного течения u1 = V – u << V (рис. 13.2). Вспомним решение для ламинарного следа. В разделе 9, не принимая каких-то специальных предположений о том, что течение ламинарное, мы получили из теоремы импульсов некоторый закон сохранения, справедливый и для турбулентного течения в следе:
Ґ |
1 |
|
|
|
т u1dy = |
CwdV = const . |
(13.14) |
||
2 |
||||
-Ґ |
|
|
||
|
|
|
Из оценки интеграла, входящего в (13.14), следует, что хорошим приближением будет равенство
u1maxb = const. |
(13.15) |