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

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

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

ние геоетрофического ветра совпадает с осью х. Считая для простоты, что и si = v ^2 = 0, получим

tg,3 =

Vo

1 -\-m а*

(6.3.11)

 

Uq.

 

Из анализа (6.3.10) и (6.3.11)

видно, что в северном полуша­

рии (шг>0) дрейф льда происходит вправо от направления геоетрофического ветра, т. е. в сторону более высокого давления. Интересно, что угол ß зависит только от массы льда, плотности воды и коэффициента турбулентности в воде. Оценим угол ß

при некоторых

характерных

значениях

величин,

входящих в

(6.3.11).

 

 

г/см3,

k2= 100 см2/сек,

 

 

Допустим,

что

р2 = 1

ф = 90°,

о>, = 7,29• 1 (Г5

1/сек,

р, = 0,9 г/слі3,

/г = 3

м,

5=1

км2.

Тогда

т — 2,7 • ІО12 г/см2, а2 = 8,5 - І О 4, tg ß= —0,187

и ß= ll°.

Наблю­

дения показывают, что ß = 5

— 15°,

 

 

 

 

 

Соотношение (6.3.11) можно использовать и для обратной задачи: по измеренному углу ß и известным tp, т, р2 найти сред­ ний в слое трения коэффициент турбулентности

H -sin 23

(6.3.12)

sin2ß

 

На основании (6.3.10) и (6.3.11) установим связь между ско­ ростью дрейфа льда и углом отклонения дрейфа от геострофи-

ческого ветра. Обозначим

т— 2=Ь

и выразим входящие в

(6.3.10) величины b и (1+6)

р2

через tg ф

 

 

 

tgß

14-b--

1 -)- tg В•

 

 

 

1 + tg ß ’

1

В таком случае, подставляя эти соотношения в формулы для

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

-о •

 

А і

У7U + tg ß)2~Hg2ft U~Hg ft)2

 

k2 '

 

l+ tg 23

 

некоторых преооразовании

 

или, после

 

° і - ѵ

 

г

(6.3.13)

 

 

: Cg2 +

 

 

 

где с „2 — модуль скорости геоетрофического*течения.

119

Для малых величин cg2 соотношение (6.3.13) можно пред­

ставить в виде

 

 

 

c0 = AG,

(6.3.14)

где А «= I

— .^ — коэффициент

дрейфа. З а ­

висимость Л от целого ряда

параметров (плотности

воды и воз­

духа, коэффициентов турбулентности, угла ß) объясняет замет­ ные расхождения в величинах А, определенных по наблюдениям за дрейфом льда в различных физико-географических условиях. Оценим величину .4 для некоторых характерных значений вхо­

дящих в нее параметров: ß = 15°,

рі = 1,3 • ІО“3 г/см3, р2=1 г/см3,

&і = 4- ІО4 см2/сек, &2 =Ю 2 см2/сек.

Для этих условий А =3.2- ІО“2,

т. е. скорость дрейфа льда составляет около 3% от скорости геострофического ветра, что неплохо согласуется со средними оценками, полученными из наблюдений.

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

о. ио

В заключение этого параграфа остановимся еще на одном интересном факте, вытекающем из теории дрейфа льда. Если в пограничном слое над неподвижной поверхностью наблюда­ ется правое вращение ветра с высотой, то над дрейфующим льдом, в слое толщиной несколько метров наблюдается левое вращение ветра, которое затем переходит в обычное правое вра­ щение. Это явление связано с тем, что при zi = 0 воздух движет­

120

ся со скоростью и в направлении дрейфа льда, т. е. влево от геострофического ветра. В таком случае годограф скорости вет­ ра будет иметь вид, изображенный па рис. 45.

§ 4. Трансформация воздушных масс под влиянием подстилаю- * щей поверхности

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

суша

море

Рис. 46. Качественная картина трансформации воздушного потока при переходе с суши на море

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

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

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

влажности и т. д.

С учетом этого сложного взаимодействия

между изменением

различных свойств воздушной массы, дадим вначале наиболее полную постановку задачи о трансформации.

Выпишем замкнутую систему уравнений для планетарного пограничного слоя атмосферы, сохраняя в левых частях полную производную от свойства

du

д

,

ди

,

п

 

1

дп

-77- =

-г-

k

т р + 2 < V > --------

дх

dt

dz

 

dz

 

 

г

 

[j

dv

d

,

dv

 

 

 

-

1

"P.

Tt

-n—k

---- 2wji

 

dz

 

dz

 

 

 

 

 

dy'

 

du

 

dv

 

.

dw

 

 

 

 

dx

'

d v '

 

dz

 

 

ß?0 _

ö

 

dtd_

 

 

1

dR

 

jn_ '

dt ~

dz

T dz

 

'

pcp

dz

 

pc0 '

d q

 

d

 

dq .

m _

 

 

dt

dz 9

dz ~t~

p ’

 

 

 

g

,

 

 

f?0

 

b'“ , d , db

 

~ T ' kj’~dz ~ C~k ■*" ~dz kb 7z’

 

 

 

M

 

 

1 ' b\

 

 

 

 

 

 

 

 

'j

 

 

 

 

 

 

 

 

 

d'j d:

 

 

(6.4.1)

(6.4.2)

(6.4.3)

(6.4.4)*

(6.4.5)

(6.4.6)

(6.4.7)

(6.4.8)

 

 

du\

■dv\*

g

dt)

1

d

db

(6.4.9)

 

 

.dz 1

.dz!

J~ ~T *X*~dz

Т

~dz

b~dz'

 

 

 

dS

dS ,

dS

.

dS .

 

dS

(S — температура,

влаж-

dt

dt

dx

ду

 

dz

 

 

 

 

 

ность, скорость ветра и т. д.).

 

 

 

 

 

 

* Нетрудно показать, что первое начало термодинамики (2.1.І0) можио

записать

в

виде

dQ=cp • -q~dt) (так как Ср~^~ dB — cp-dT

—dp),

откуда

при

Г /6

следует

(6.4.4).

 

 

 

 

 

122

Если считать известными градиенты давления, радиацион­ ный и фазовый притоки тепла, а также фазовый приток влаги, то для определения 9 неизвестных (и, v, w, Ѳ, q, k, b, /, ф) име­ ется 9 уравнении н система замкнута.

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

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

[~щ—oj • Предположим, что начало координат находится на границе раздела между двумя поверхностями, ось у направлена

вдоль границы р а зд е л а (^ = о ) а ось х направлена вдоль вет­ ра (рис. 46).

Задача о трансформации поля температуры и влажности

Предположим, что в-атмосфере не происходит фазовых ,пе­

реходов

воды (

L— t=0 , - -

= о)и что можно пренебречь

радиа-

 

\

р

р

 

/

 

ционным

притоком

і

1

dR

верти­

тепла)—

• - ^ = 01 и упорядоченным

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

нулю и для не очень больших высот

можно считать'

малыми

di

dq

 

и ранее сделанных

допущений

w т и г ф . С учетом этих

dz

dz

(6.4.5) можно записать в виде

 

 

уравнения (6.4.4),

 

 

 

 

ДЯ

 

д%

 

(6.4.10)

 

 

а • дх

 

dz’

 

где ■&— потенциальная температура

или удельная

влажность;

а» — соответственно

а г или а (/

(в дальнейшем будем

писать а

без значка 6).

Сформулируем теперь граничные условия. Бѵдем считать,

что над исходной подстилающей поверхностью (х<0)

воздуш­

ная масса находилась достаточно долго и Ь на любой

высоте

равно величине О на поверхности — іѴ В таком случае

 

Д| =»>!.

(6,4.11)

х = 0

 

2»0

123

С другой стороны, на большой высоте над новой подстилающей поверхностью воздушная масса также будет иметь A = öi, так как влияние подстилающей поверхности ограничено слоем ко­ нечной толщины, т. е.

» I

(6.4.12)

Z-- ОО

х>0

.Предположим, что нам известно значение flна новой под­ стилающей поверхности (л;>0) — Ао и что fl0 не зависит от х

» I

= А 0.

(6.4.13)

Z

О

 

X. О

 

Граничное условие (6.4.13)

в действительности следовало бы

заменить уравнением баланса свойства fl на подстилающей по­ верхности, так как в процессе трансформации меняется не толь­

ко fl, в воздушной массе, но и flo .

Если предположить,

что в на­

чальный момент ? = /0

(перед приходом

воздушной

массы

на

-л'>0) подстилающая

поверхность

имела

свойство fl0n = const,

то

в дальнейшем за счет обмена свойством между воздушной мас­

сой и поверхностью (различного

на

разных

расстояниях от

х = 0) установится flo = flo(*). С

другой

стороны,

пока недоста­

точно надежно определяется температура поверхности * Ѳо и со­ ответствующая ей удельная влажность q0 (для морских условий насыщающая удельная влажность). С этой точки зрения нельзя считать flo известной величиной, и при 2 = 0 нужно задавать ус­ ловия теплового (для температуры) и водного (для влажности) баланса

 

R = P + LE + B,

\

(6.4.14)

 

М = Е,

I

 

 

компоненты которых

(R — радиационный

баланс, В — теплооб­

мен с нижележащими

слоями, М — осадки, В — скорость испа­

рения и так далее) зависят от .ѵ. Однако

использование (6.4.14)

заметно усложняет решение задачи, к тому же, если воздушная масса поступает на водную поверхность, есть основания считать

fl°o = Ao = const (независящим от х).

Это

упрощение связано с

тем, что удельная теплоемкость воды

почти в

4 раза

больше,

удельной теплоемкости воздуха (т. е.

большим

изменениям тем­

пературы воздуха будут соответствовать

сравнительно

слабое

* Из-за больших градиентов температуры в самом верхнем слое почвы температуру поверхности суши стандартными приборами измерить практиче­ ски невозможно. В море за поверхностную температуру обычно принимают температуру, измеренную в слое от 0,5 до 1,5 м глубиной. При слабом ветре ома может почти на Г С отличаться от температуры поверхности.

124

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

й0 = const.

первого приближения^ что скорость

Предположим, в качестве

ветра и коэффициент турбулентности

не зависят

от высоты.

В таком случае (6.4.10) примет вид

 

 

и іт

 

52Я

\

■OL-k-

(6.4.15)

Ö X

 

dzv

 

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

*2

 

' = 4 г -

(6.4.16)

Перепишем частные производные по х и г через полные произ­ водные по в

 

С/11

db

öl

 

z2

 

db

 

 

 

 

дх

dl

дх

 

 

X1 '

Ir:

 

 

 

 

_

db

д'

_

2z

'

db

 

 

 

 

dz '

 

dl

dz

 

X

dl’

 

 

 

д

_

 

д / 2z

_ d»\

 

2_ cß t

_4z_2

 

dz \ dz) ~

dz\ X

dl)

 

X

dl

1

x l

 

и подставим полученные

 

выражения

в

(6.4.11),

которое после

некоторых преобразований примет вид

 

 

 

 

 

 

dl2

 

-•

+

- L i *

= о.

 

 

(6.4.17)

 

 

A'j.k '

211dl

 

 

 

 

Для новой переменной

 

 

граничные

условия

(6.4.1

(6.4.13)

перепишутся как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

&('): ! о=»о;

 

 

 

 

(6.4.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

Ч=* 00

 

 

 

 

 

 

(6.4.19)

 

 

 

 

 

 

 

 

 

 

 

Проинтегрируем (6.4.17)

один раз

 

 

 

 

 

 

 

4-J-k

е

125

Полученное выражение проинтегрируем от 0 до £ й исполь­ зуем (6.4.18)

 

Ч -) — V = c ' f

е

d' 1

 

(6.4.20)

 

 

 

 

Г Т Т

 

 

 

 

о

 

 

 

 

 

 

 

 

 

Введем новую функцию

 

 

 

 

 

 

-2 .

и-,.

 

 

 

 

 

 

4aJfe ’

 

 

 

тогда

 

г

I

I

 

a

8-а-А

/

d : , = ----------

ч

_,/■’«

с

V

4 ^

4

11

 

 

 

У n

При

:=о,

 

- 0

,

 

 

 

 

 

 

 

 

 

5_

1 /

j E l

 

 

 

 

 

у

4а£

 

и (6.4.20) примет вид

] ür

A a k

(6.4.21)

Ч')=&о+*Ѵ f ' е ’do’

4a&

где cx—2 ‘C- jy

и

Интеграл, входящий в (6.4.21), известен в математике как инте­ грал вероятности и затабулирован. Известно, что

00

_,.0

 

1'

■*“

(6.4.22)

.(«

*

= - £ ■

 

О

 

 

для определения С\

 

Используем (6,4.22) и

(6.4.19)

 

 

л _.ft

 

У it

 

 

с . - __

 

 

wi

ио — Іу1

2

 

или

„ .

(» i-» o )2

 

 

 

 

 

 

cl —

 

— = .

 

V к

В таком случае имеем окончательное выражение для опре­ деления 0- в любой точке над новой подстилающей поверхностью

126

Aakx

 

„з

( 6 .4 .2 3 )

1

e

-d*.

 

Из физических соображений ясно, что заметное влияние под­ стилающей поверхности распространяется не до бесконечности, а до какой-то конечной высоты — высоты внутреннего погранич­ ного слоя, которую можно определить как уровень, где t f s f h _.

Известно, что интеграл вероятности становится близким к - L -

\же при верхнем пределе около 2—3. Таким образом, обозначив в (6.4.23) верхний предел, обеспечивающий заданную точность равенства іТі, через N, получим условие для определения высоты внутреннего пограничного слоя

 

 

 

(6.4.24)

Если N = 2,0; м= 1 0

м/сек-, а=1,0;

k — 5 м2/сек, то на расстоя­

нии 1 км от

границы

раздела

90 м, а на расстоянии

10 км — 284 м.

 

 

В (6.4.24)

прямая зависимость от J х и обратная от ) и свя­

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

Используем (6.4.23) для определения турбулентного потока тепла и затрат тепла на испарение. Ограничимся для примера выводом формулы для расчета скорости испарения или затрат тепла на испарение

2=0

Определим производную из (6.4.23), подставив вместо & значение q и использовав правило дифференцирования интегра­ ла, пределы которого зависят от параметра

д ь Ь)

b

dj {x, a)

Ж j

J

а (а )

a (»)

 

В нашем случае

 

 

127

d g

и

 

 

Яо .

d z

\ x q k x

 

1

 

z= -0

 

 

 

 

LE — - u L - k - x q

 

 

 

 

=(?o—</i)-p-|

kll OLq

^

(6.4.25)

K-X

 

 

 

 

 

Увеличение скорости испарения при росте скорости ветра объясняется ускорением замены увлажненного воздуха новыми порциями сухого (при q0>qi), тогда как обратная зависимость

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

Влиянием X на скорость испаре­ ния объясняется зависимость пос­ леднего от направления ветра и формы испаряющего бассейна (рис. 47). Например, при прямоугольной

 

 

 

форме бассейна скорость испарения

 

 

 

со всего бассейна для ветра, дую­

 

 

 

щего вдоль меньшей стороны Е\,

 

 

 

будет больше, чем для ветра, дую­

 

 

 

щего

вдоль

большей

стороны

£?,

Рис. 47. Влияние кон­

так

как

во

втором

случае

воздух

проходит

большее

расстояние

над

фигурации

бассейна

на

испарение

 

водой и успевает сильнее увлажни-

Получим

 

 

ться.

 

 

испарения с

полосы

выражение для определения

шириной Ау и длиной L, расположенной вдоль оси х. Если обоз­

начить испарение

с

элементарной

площадки

dx-dy

через

dEI = Е • ах dy, то

 

 

 

 

 

 

 

 

 

 

Еі= Ы Ч о -Я іУ \ /

~

?- -Аѵ-1 Т .

 

(6.4.26)

По аналогии с (6.4.25) и (6.4.26) можно получить выражения для турбулентного потока тепла и потери тепла с полосы шири­ ной Ау и длиной L

128

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