Файл: Мельникова И.И. Динамическая метеорология учеб. пособие для океанологов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 12.07.2024

Просмотров: 94

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.
д<ѵ дх

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

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

тицы вокруг вертикальной оси,

обозначается Й, и связана с

составляющими скорости ветра соотношением

 

()

дѵ

ди

(7.1.9)

дх

ду’

 

 

 

Q, имеет положительный знак при движении против часовой стрелки, т. е. в циклонических областях, и отрицательный —- в антициклонических.

Действительно, согласно рис. 56, в области низкого давления

ф

и

*і - ^8 > 0; 2г

 

 

ди

l i t - и4

 

 

 

 

 

 

 

д у ~

 

2г

< 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:довательно, Й,

 

 

 

 

 

 

 

В

свободной

атмосфере,

 

 

 

 

 

где выполняются

соотношения

 

 

 

 

 

(7.1.3),

(7.1.4), с

достаточной

 

 

 

 

 

точностью можно

принять

 

 

 

 

 

р

1 /д°-Ф д*Ф\

1

 

 

 

 

 

I

\дх'г

'

 

-1-

I V2 Ф.

 

 

 

 

 

 

дуд!

 

 

 

 

 

 

 

 

 

 

 

 

Рис. Ü6. Распределение ветра

в обла­

 

 

 

 

 

 

 

 

сти

низкого давления

 

 

Т-акнм образом, поле вихря однозначно связано с полем гео­

 

потенциала (или давления). -

 

 

 

 

ч

 

Уравнение

(7.1.2) для эволюции й г называют уравнением

вихря.

Впервые

оно

было

получено

и

проанализировано

 

 

А. А. Фридманом и в настоящее время является основным прог­

 

ностическим уравнением.

 

 

 

 

 

 

Согласно

равенства

(7.1.9),

уравнение

вихря может

быть

 

получено из уравнений движения, если первое из них дифферен­

 

цировать по у,

а второе — по х и вычесть

из

второго

первое.

 

Тогда после несложных

преобразований получается

 

160


dt "~г dx ‘ dy ~г " !!\дх ~ г'

ду)

р1 [дх ду

ду дх)

, I du

,

dv\

dl ■

dl

 

- И т -

+

7

—U ----- Ѵ ~ .

 

\дх

 

dvl

dx

{dy

 

Поскольку I — 2е> sin ср

можно считать величиной постоянной

во времени, это уравнение можно переписать в следующем виде:

dQ

dQ

 

dH

 

 

-1 st —I—//

 

-L'v -ІЛ :

 

d t +

дх ^

ду

 

 

 

 

 

1 (dp

dp

do dp\

 

 

 

\âx

ду

dy dx)’

Величина Qrt = Q ,,+ /, называемая абсолютным вихрем, представляет вращение частиц воздуха относительно неподвижной системы координат, не связанной с землей. Эмпирические дан­ ные показывают, что Q. в несколько десятков раз меньше /, по производные от о г, как правило, значительно больше производ­ ных от / (так, если ось х направлена по шпроте, то dl/dx = 0 и, как уже говорилось, dl/dt = 0). Поэтому последнее уравнение может быть записано

dila

.іди

.

d v

1 ,

1 [ d p

d p

d p

d p *

dt ~

'\дх ~ т d y

) + f > * ( d x d y

d y

d x ) '

Последнее слагаемое в равенстве (7.1.10) можно переписать,

используя геострофические соотношения и -----~

v = jj- X

dp следующим

образом:

 

 

 

 

 

J j ' Ф ЩР _

dy

2 E \ - J - (

È L 4- 7, A n

Ръ\'дх dy

д х ! ~ "

р \

d x ^

дуІ[

Как показывает выполненное преобразование, это слагаемое представляет собой (с точностью до коэффициента 1/р) адвек­ цию плотности. Из опытных данных следует, что оно значитель­ но меньше первого слагаемого правой части равенства (7.1.10).

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

11

16.1


Если подставить в уравнение (7.1.2) соотношение (7.1.8), го нетрудно получить

 

 

 

І|5= м , - №

 

 

 

 

dt

 

 

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

 

 

 

Г

дФ

І А * - І Ю .

(7.1.11)

 

 

 

dt''

Здесь А* —

dQ

 

 

— адвекция вихря.

а

дх

+

V

Если из уравнения

(7.1.11)

найти Ф,

то на основании (7.1.3)

п (7.1.4) нетрудно найти поле ветра, па

основании (7.1.6) — по-

• ле температуры,

а

из

(7.1.5)— вертикальную скорость.

Таким образом, уравнение (7.1.11) совместно с уравнениями (7.1.3) — (7.1.7) составляет систему прогностических уравнений, позволяющих находить будущие поля основных метеоэлементов, если известно их значение в данный момент.

§ 2. Прогноз давления на среднем уровне

Рассмотрим простейший случай решения поставленной зада­ чи— прогноз поля геопотенциала (что идентично полю давле­ ния) на среднем уровне.

Опытные данные показывают, что в атмосфере на высоте 5—6 км дивергенция скорости очень мала. Этот уровень, где можно пренебречь дивергенцией, называют средним уровнем и равенство (7.1.11) для него записывается

дФ

(7.2.1)

^ T t =1А:;,

что означает: изменение поля геоиотенциала на этом уровне оп­ ределяется только адвекцией вихря.

Если известно распределение Ф(л:, у) в момент t,

то величину

А :, нетрудно рассчитать и использовать уравнение

(7.2.1)

для

отыскания

<9Ф

л

интервала времени At

этѵ

.

В течение короткого

величину можно считать постоянной,

и определить

поле Ф (х,у)

в момент tA-At, считая, что

 

 

 

 

Ф , у, tA-At) =Ф (х, у,

,0 + ^ А t.

(7.2.2)

162


По вычисленному «юлю Ф (х,у)

можно вновь

рассчитать А<:,

дФ

 

 

— и Ф (х,у) в момент t -j- 2At.

 

 

Такой метод называется методом шагов по времени. Он дает

тем более точные результаты, чем меньше шаг

At (поскольку

равенство (7.2.2). представляет

разложение в ряд

по параметру

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

ляется интервал в один час.

 

уравнения

(7.2.1),

Остановимся теперь на интегрировании

которое необходимо выполнить, чтобы,

имея

дф

'

А ■>, получить-^ .

 

Для получения однозначного решения надо задать граннч-

ные условия, т. е. надо знать значения

дФ на границах области

 

 

д Ф

г

 

интегрирования. Допустим, что нам известію —^ и равно

Г ^ на

окружности радиуса R вокруг точки, для которой дается прог­ ноз.

Для интегрирования уравнения (7.2.1) его удобнее записать в полярных координатах. Поскольку известно, что

 

дг

1

1

+ 1

 

 

дг'і

г дг ' г2 д\}-’

 

где г (радиус-вектор) и #

(полярный угол) — координаты точки

в полярной системе.

 

 

 

 

 

 

Уравнение (7.2.1) записывается

 

 

 

ff-F_

J _

OF

J

_

cFF_

 

d r 2 '

r

i) r

^ r2 дЬ'-

 

дФ

 

 

 

 

 

 

где F = ~ді'

 

 

 

 

 

 

Полученное уравнение умножим на г/й, проинтегрируем по О

в пределах от 0 до 2л и разделим па 2я

 

1 2?d2F

öF,

1 2f

1

2-

2- ,

2~д' дг2

 

 

 

 

■sPrf,b

Поскольку г и ö — независимые

переменные, поменяем мес­

тами дифференцирование и интегрирование. Тогда

 

 

± ( ±

2-

 

2

д г \ 2 -

 

 

Fdb

-J lA9dti,

 

дг

2т:

 

 

163


так как

dp, 2-

= О (â — О и 0 = 2л соответствуют одной

 

о

точке).

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

1

2?

_

1

.

\ I'd' '

--г / і

2~

I l-A'jCf-ty LAs.

“ ’ о

 

 

Іі

Тогда уравнение принимает вид

</-/*’ , 1 оГ

,

дг* г

и может быть переписано в полных производных ,

d\F

,

1

rf/’"

ІА ,—-0.

(7.2.3)

rfr2

r

r

dr

 

 

Это линейное неоднородное уравнение второго порядка. Для его решения требуется два граничных условия.

Первое условие запишем на основании заданного значения

<ЗФ, на r = R, а именно

дГ

? (г ) I

2 -J

dl).

 

r — R

dt.

 

 

R

 

Нетрудно понять, что----:

|I

при разных значениях §

будет

<п г

R

 

 

иметь разные знаки и при достаточно большом R значение ин­

теграла будет близко к нулю, а потому можно принять

 

F(r)

 

г0.

(7.2.4)

 

г~rR

 

По физическому смыслу F(г) не может быть бесконечностью. Это свойство и используем в качестве второго граничного усло­ вия, а именно примем, что

F(r) I

='Д(0) Ф<х>.

(7.2.5)

г— 0

 

Таким образом, оба граничных условия сформулированы.

164