Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Динамика вязкого газа, турбулентность и струи

.pdf
Скачиваний:
112
Добавлен:
14.06.2019
Размер:
4.94 Mб
Скачать

12.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) опущен член Px , так как в первом приближении в задачах о свободной

турбулентности давление можно считать постоянным.

Для того чтобы проинтегрировать уравнения (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 произ-

водная 2uy2 терпит разрыв. Такой недостаток присущ всем решениям, по-

лученным на основе формулы Прандтля (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)

Соседние файлы в предмете Механика жидкостей и газов