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

книги из ГПНТБ / Мельникова И.И. Динамическая метеорология учеб. пособие для океанологов

.pdf
Скачиваний:
36
Добавлен:
23.10.2023
Размер:
5.62 Mб
Скачать

 

Эта сила называется си­

к

лой

барического

градиен­

та

(вернее

ее проекцией

 

 

на ось х). Аналогично F

 

можно получить и другие

 

составляющие силы бари­

 

ческого

градиента

 

F,ру~

L

ÈP

р _

 

!>

dy'

I'z

 

 

 

1

Op

 

 

Рис. 5. Действие давления на эле­

 

Р

dz

 

 

ментарный объем

 

 

 

 

Главной причиной воз­ никновения силы бариче­ ского градиента в атмосфере является неравномерность нагре­

вания подстилающей поверхности.

*

Сила трения

Рассмотрим поток жидкости, в котором скорость растет с высотой (рис. 6). За счет теплового движения молекулы с уров­ ня 1 могут попадать на уровень 2 и переносить сюда некоторое

дополнительное

количество

движения.

И

наборот,

молекулы с

\ ровня 2 будут

переносить

на уровень

/

меньшее

количество

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

ний ясно,

что поток будет

тел? больше, чем больше градиент

z

 

скорости. Поток количества движе­

 

 

ния можно рассматривать как отне­

 

 

сенную к единице поверхности каса­

 

 

тельную силу, называемую каса­

 

 

тельным напряжением т. Итак, за

 

 

счет молекулярного перемешивания

 

 

в потоке жидкости с вертикальным'

 

 

градиентом скорости возникает ка­

 

 

сательное напряжение тм

 

Рис. 6.

Молекулярное

de

(2.2.8)

dz’

перемешивание

 

20

где коэффициент пропорциональности р называется динамиче­

ским коэффициентом вязкости и зависит от свойств жидкости;

Iл]=ЛЬ/г‘-Г 2, [uj = ЛЬ/ / 1Г 1

Для большинства атмосфер­

 

ных движений главное значе­

г

ние имеет не молекулярное,

а

 

тVрбуле«гное перемешивание,

 

при котором свойства перено­

 

сятся в потоке не отдельными

 

молекулами, а частицами жид­

 

кости значительных размеров

 

(более подробно этот вопрос

 

будет рассмотрен в следующем

 

параграфе). В случае турбу­

 

лентного перемешивания

по

 

аналогии с (2.2.8) можно счи­ тать, что

Рис. 7. Действие касательных напря­ жений на элементарный объем

d e

(2.2.9)

"ch’

 

где А г — коэффициент турбулентной вязкости или турбулентно­ го обмена вдоль оси z. Обычно А г на много порядков величин

больше р || т. е. (р+ДД.

. Наряду с А г часто вво­

дят понятие коэффициента

турбулентности кг -

Д,

и тогда

 

■'j-k

de

(2.2. 10)

Tz

 

Физический смысл Аг и к г будет пояснен позже. Сначала полу­ чим выражение для проекции силы турбулентного трения. Для этого рассмотрим элементарный объем dx-dy-dz, расположен­ ный в потоке, скорость которого растет с высотой. Обозначим че­ рез ті и Т2 касательные напряжения (рис. 7), действующие на

нижнюю и верхнюю

грань

объема вдоль оси х: ті = тДг),

Т2 = т x(z + dz) и разложим

в ряд

•а “ Яг

—д (а )

21

В таком случае силу трения, действующую "на объем вдоль оси •V, за счет касательных напряжений т2 и п можно представить как

f-..x=(~2—’ i ) d x ’dy

или

 

 

f =

dz

dx-dy-dz,

J~x

 

(Здесь и в дальнейшем нижний индекс указывает проекцию си­

лы трения,

а

верхний — проекцию

градиента скорости).

Сила /5

,

отнесенная

к единице массы,

запишется

в виде

 

 

 

 

 

 

 

 

'1

1

д

/

ди

(2 .2 .1 I 1)

 

 

 

 

 

р

dz

 

1

1**

\ A*Tz

 

 

 

 

 

 

 

По аналогии

можно

найти

выражение

и для Р*

F2

 

 

р *

1

 

 

1

 

d

 

 

 

 

 

 

1

-у ---

p

dz

~

p

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

(2 .2 .1 1 ")

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

p z

 

= J _

±

(

 

 

 

 

 

 

~Z ---

p

dz

 

P

dz \

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим

теперь случай,

когда

поток имеет горизонталь­

ный градиент скорости. Тогда можно получить выражения для остальных составляющих

 

1

dzx

 

1

d

 

du

 

 

p

dx

 

p

dx

•А* dx'

 

 

1

dzy

 

1

â

 

dv

 

/ • - у -

p

 

 

fj

dx *

dx'

 

р-.*-=

1

dz.

 

1

d

Лл

dw

 

p

dx

~~ p

dx

dx'

(2.2.11 "О

 

1

Özv

 

1

â

•Ay du

 

'

 

 

p

ây

p

Jy

 

W

 

 

 

 

 

 

 

 

 

1

 

 

1

d

 

dv

 

 

p

 

 

p

dy Л

Ty'

 

F l =

1

 

 

1

d

Л

dw

 

 

p

dy

 

p

dy

"'y

dy

 

где Ах, А у — коэффициенты

турбулентного обмена вдоль осей

х и у.

 

 

 

 

 

 

 

 

22

Воспользовавшись полученными выражениями для сил, дей­ ствующих в атмосфере, (2.2.3), (2.2.6), (2.2.7), (2.2.11), перепи­ шем теперь уравнения движения (2.1.5) в развернутом виде

öu

 

du

{-V

öu ,

 

du

1

dn

,

1

dt +

«

dx

Ту

 

' Tz ~~

0

—— f-

 

 

 

dx

1

 

 

 

.i- d ,

du

,

d ,

Oll

d ,

du

 

 

 

1

dx

x dx

'

T ? y ~Ty -~r

dz

- dz’

 

dv

-f и

dv t-д' to ,

dv

1 A _ _

2wj

dt

dx

 

dy

dz ~

P

dy

 

 

 

i A b

dv , A u dv ± i k dv

 

 

 

dx

dx ^

dyky Ту

 

dz

; T z’

dw

4- и

dw

-f V dw ^ dw

1

dz

_£ -

~dt

1

dx

 

dy '

dz

0

 

О

 

 

 

д ,

dw

д , dw

,

д

 

dw

 

 

 

Т х х дх + д у у Ту " ‘ ö z ! Tz

( 2.2. 12)

*

I

Поскольку в (2.2.12) содержится новая неизвестная харак­ теристика — коэффициент турбулентности, то необходимо более подробно познакомиться с явлением турбулентности и найти спо­ соб замыкания системы уравнении.

§ 3. Понятие оо атмосферной туроулентности. Выражения для турбулентных потоков тепла, влаги и количества

движения

Если построить график зависимости скорости ветра, впрочем как и любой другой характеристики атмосферы, от времени ѵ поверхности земли (рис. 8), то будет видно, что скорость движеиия пульсирует, т. е. и резко изменяет свою ве­ личину и направление в течение коротких проме­ жутков времени. Режим движения, при котором отдельные частицы жид­ кости или газа имеют не­ правильные-, хаотические

траектории с

поперечны­

 

ми и даже

обратными

 

(по отношению к общему

Рис. 8.

23

движению) перемещениями, называется турбулентным. Из гид­ ромеханики известно, что характер движения жидкости или газа зависит от безразмерного числа Рейнолыса

Re = ^ ,

ц=ѵр,

(2.3.1)

Здесь и — характерное

значение

скорости

потока; L — ха­

рактерный размер потока;

р —- динамический

коэффициент вяз­

кости; р — плотность. При

малых Re движение имеет ламинар­

ный характер (частицы перемещаются параллельно друг другу по плавным траекториям), при больших Re движение имеет тур­ булентный характер.

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

Из рис. 8 видно, что мгновенную скорость ветра (скорость в

данный момент времени) можно представить как

 

и = и + и',

(2.3.2)

где и — средняя за некоторый.период времени скорость; и' — от­ клонение От средней. Так как средняя кинетическая энергия по­

тока пропорциональна и2, то с учетом того, что

(й:г)= /Г2 + 2тг'+ и4 = и2 у и

(2-3-3)

получаем: средняя кинетическая энергия потока

складывается

из кинетической энергии оередненного, сглаженного движения и кинетической энергии пульсаций, называемой энергией турбу­ лентности (Ь).

При получении этого вывода мы исподьзовали один из посту­

латов статистической теории турбулентности /•/' = () (считается, / = /+ /'). Напомним и другие постулаты:

fi+f2=f+l2; j=f,

fi */2 = /1 */2 + /7i ’ /V2

(2.3.4)

(если /1 и /2 независимы, то / 'і ■/^ = 0),

0 .

Все постулаты довольно легко доказываются. Например, справедливость последнего можно показать так:

//==/ — Т; ? = 7 —7= Г — F=o.

. ;

24

По аналогии с процессами молекулярного перемешивания, когда переносчиками свойства являются отдельные молекулы; в теории турбулентности вводится понятие турбулентного моля— отдельной частицы жидкости или газа, которая отрывается от общего потока в одной его точке и смешивается с ним в другой (перенося таким образом его свойства из одной точки в другую).

Хотя практически смешение идет непрерывно, вводится по­ нятие пути смешения (аналог пути свободного пробега в моле­ кулярной теории) как расстояния /, которое проходит турбулент­ ный моль от момента зарождения до полного смешения с окру­ жающей средой. Турбулентный моль представляет вихрь, воз­ никающий за счет гидродинамической неустойчивости основного потока. Прохождение через данную точку пространства вихрей разных размеров, несущих свойства из различных частей основ­ ного потока, и создает ту сложную картину изменения метеоро­ логического элемента, которая показана на рис. 8.

В дальнейшем мы будем рассматривать только вихри с го­ ризонтальной осью, определяющие перенос свойства в верти­ кальном направлении.

Получим общее выражение для турбулентного потока любой субстанции S. Если понимать под S количество субстанции в единице массы воздуха, то для потока тепла S = c BQ (Ѳ — потен­ циальная температура), для потока водяного пара S = g (удель­ ная влажность), для потока количества движения S = c (с—ско­ рость потока).

Рассмотрим на уровне г единичную горизонтальную площадку а, через которую снизу вверх проходят вихри со скоростью W. За единицу вре­ мени через площадку пройдет количество субстанций S, со­ держащееся в объеме aw, т. е. Spcrco (рис. 9).

/ '

7

ѵѵ

■-1

Рис. 9. Перенос свой­ ства через площадку

Так как под потоком субстанций Q обычно понимают ее ко­ личество, переносимое за единицу времени через единичную площадку в направлении нормали к ней, то'

ß =

\ ( р wSdt.

(2.3.5)

 

1 о

 

Выразим w и S через

средние величины

и флюктуации:

w — w + w'; S = S + S'. Тогда (2.3.5) примет вид

25

f p w Seit -f-

f p tc' S' di -j- f pw' S dt -f

 

 

о

0

 

- f \ pw' S' dt

—pw S ~j~ pw'S'.-

(2.3.6)

ö

 

 

 

Первый член этого соотношения определяет поток субстан­ ций, обусловленный средним движением, остальные три члена дают поправку за счет турбулентности. На основании ранее рас­ смотренных постулатов вторым и третьим членами можно пре­ небречь. Четвертый член обычно не равен нулю (так как а / и S' взаимосвязаны). Таким образом, поток свойства 5 складывается

из потока, обусловленного средней скоростью Qi = pSw, и тур-

Пульсацию свойства на уровне z можно представить как разность между свойством вихря и среды. Если вихрь пришел с уровня z — / и из пути сохранил свои начальные качества, то

он принес свойство S = S(z I) и тогда S' = S(z — / ) — S(z), при условии, что на исходном уровне вихрь обладает свойством этого уровня.

Для достаточно малого / —пути смешения S ( 2 I) можно разложить в ряд

S ( z - l ) = S ( z ) - ^ z i

Втаком случае

Вобщем случае следует считать, что / может быть различной для разных субстанций, ибо оно само зависит от градиента суб­ станции (действительно, вихрь, первоначально обладающий ко­ личеством субстанции S, должен пройти различное расстояние, чтобы смешаться с окружающей средой по содержанию количе­ ства движения, тепла или влаги, в зависимости от градиента этих субстанций). С учетом этого формулу для турбулентного потока субстанции S следует записать в виде

(2.3.8)

26

Из физических соображений можно считать, что w'-l, явля­ ется характеристикой турбулентного перемешивания, ее обоз­ начают через к, и называют коэффициентом турбулентности для потока субстанции S

Qi

и

д $

(2.3.9)

- A

- S -

 

В таком случае турбулентные потоки количества движения, тепла и влаги можно выразить так:

 

а ■k •

Г) С

 

 

 

 

dz:

с/Ѳ

 

р*ггр [ dz + ъ ];

dz

 

г-

,

 

Ö(1

£ = —pAf0

 

dz

 

1

4

 

В формуле (2.3.10) знак выбран из соображений

(2.3.10)

(2.3.11)

(2.3.12)

удобства

записи тѵ так как

в атмосфере

дбычно dc/dz>0. С учетом

(2.3.10—2.3.12)

притоки

тепла

и

влаги, входящие в уравнение

притока тепла

(2.1.18 или 2.1.21)

и притока влаги (2.1.23) мож­

но рассматривать

как

дивергенцию

соответствующих потоков

и для единицы массы представить в форме

 

 

 

 

dz

 

 

dz

 

(2.3.13)

 

 

 

 

 

 

 

 

 

 

- ? Т -

— k

4

dz

(2.3.14)

 

 

 

dz

 

 

§ 4. Факторы, определяющие интенсивность турбулентности. Уравнение баланса энергии турбулентности

Обозначим через Е\ энергию турбулентности, возникающую за счет кинетической энергии основного потока. Вихри, образо­ вавшиеся либо из-за обтекания неровностей поверхности, либо благодаря резкому изменению скорости с высотой, перемещают­ ся из слоя в слой. Если плотность окружающей среды отлича­ ется от плотности вихря, то на него начинает действовать сила Архимеда. Обусловленное ею ускорение либо совпадает с на­ правлением движения вихря, увеличивает скорость и запас его энергии, либо направлено в сторону, противоположную движе­

27

нию вихря, уменьшает скорость и запас энергии вихря. Обозна­ чим через Е2 приток (отток) энергии за счет действия силы Архимеда. Благодаря соприкосновению и поверхностному тре­ нию между вихрями часть их энергии может переходить в тепло­ вую. Такой переход называется диссипацией энергии турбулент­ ности— обозначим ее через Ез. Наконец, через Е4 обозначим приток (отток) энергии турбулентности за счет обмена вихревой энергией между соседними слоями (диффузия энергии турбу­ лентности).

С учетом введенных обозначений можно записать в общем виде уравнение баланса энергии турбулентности

f/- --/Т

: А,.

(2.4.1)

Уравнение (2.4.1) утверждает, что изменение энергии турбулент­ ности связано с притоком энергии турбулентности от среднего потока, действием силы Архимеда, диссипацией и диффузией. Получим выражения для Еи Е2, Еъ и Е4. Для определения Е\ рассмотрим взаимодействие элементарного объема dx-dy-dz

с окружающей средой (рис. 10). Обозначим через (z)c и

c{z + dz) — вектора скорости на уровне z и z + dz, а через т(-г)

Рис. 10. Действие сил турбулентного трения на элементарный объем

и т (z + dz) — касательные напряжения, действующие на ниж­ нюю и верхние грани объема. В таком случае работу, затрачи­ ваемую на преодоление сил турбулентного трения за единицу времени, можно выразить уравнением

где c = i-u+jv; %— ixx + j-x ,

 

 

—^

 

>

 

 

Разложим t(z-t-cte) и

c(z+dz) в ряд

 

1

 

 

T(z +

öfz) = T ( z ) +

dz

;

"c (z - M z )= 7 (z )+ - |^ -« rz -f . . . .

Воспользовавшись

выражением

для

"-?>

c(z+c?2 ) и x(z + fifz),

перепишем (2.4.2) в виде

W = -J—[c(z) -(z)] dz-dx-dy

или, подставив выражения для векторов скорости и касательно­ го напряжения, получим

W =

д ( у и )

(Ң-y-v)

dx • dy • dz.

(2.4.3)

dz

dz

 

 

 

При выводе (2.4.3) использовано свойство скалярного произве­

дения двух векторов: і-і = 1 и і • / = О и т. д. Для единицы массы (2.4.3) примет вид

\Ѵ =

_1_

д{*хи) I

d (ту г»)

(2.4.4)

Р

öz

c?z

 

 

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

Определим интересующую нас энергию турбулентности, воз­ никающую за счет кинетической энергии среднего потока, как разность между работой, совершаемой турбулентными напряже­ ниями, и той ее частью, которая идет на изменение кинетиче­ ской энергии среднего движения. Последнюю выразим из урав­ нения баланса кинетической энергии среднего движения, кото­ рое получается, если уравнения движения

d u _

dp

О I 1

ДІѴ

dt

dx

 

2шл» 4-----

dz

dt

p

dy

z

JLp dz

d v _

1

dp

0

 

 

 

- - 2ш и

 

29

Соседние файлы в папке книги из ГПНТБ