книги из ГПНТБ / Мельникова И.И. Динамическая метеорология учеб. пособие для океанологов
.pdf
|
Эта сила называется си |
|||||
к |
лой |
барического |
градиен |
|||
та |
(вернее |
ее проекцией |
||||
|
||||||
|
на ось х). Аналогично 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