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

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

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

du

ft

 

dt ~

Ho

dx

dv

11

d -

~dt

 

(6.6.19)

Q0 ‘ d y

dw

n

d -

dt

Ho

d z

 

 

Введем понятие о невозмущенных полях метеорологических элементов, т. е. полях, которые имели бы место при отсутствии мезометеорологчнеских возмущений связанных с термической неоднородностью подстилающей поверхности. Обозначим невоз­ мущенные значения метеоэлементов через U, V, Ѳ, II, W. Так как невозмущенное поле можно отождествить с крупномасштаб­ ной циркуляцией, для которой вертикальная скорость мала по сравнению с мезомасштабными процессами, то будем считать, что \Ѵ= 0. Запишем систему уравнений для определения невоз­ мущенных значений метеоэлементов:

DV

Ѳ

dll

 

 

Dt

H0

dx

 

 

D V

ft

d\\

 

 

Dt

H0

dy - 2 шг и + У а- и ,

 

ft

 

 

 

0 == - W0 ‘

dz

g;

 

о ѳ

-5 j _ L £ ll4 . y

ft

Dt

- 'Ф Т 'Л I -1

a

Ё Е

-j-

± ^ o

 

dx

'

dv

 

 

R - T ’

 

 

 

— «Ф ‘

a

 

где

( 6.6. 20)

(6. 6. 21)

Q _ JL j_ TJ j?_ !

dy'

 

D t ~ dt

'

дх ~

 

 

d

\r d

,

— -

dx

-r-

N l -X---

± N

dy

1 dy

1 dz

dz

149

IV/, N — коэффициенты горизонтальной и вертикальной турбу­ лентности при невозмущенном состоянии метеорологических полей. В уравнении неразрывности отброшено малое слагаемое.

Подставим в систему (6.6.15), (6.6.16), (6.6.17) и (6.6.19) вы­ ражения для значений метеоэлементов в виде суммы невозму­ щенной величины и мезометеорологического возмущения

Т — Т+Т', р ^ Р + р', 0 = Ѳ+ Г , я —П + л/,

 

р = Р + р', q — Q+ q'

(6.6.22)

к отброаим малые члены на основании неравенств

Г '<<'/",

р '« Я , 0 '< < Ѳ , я 'С С П , )Ѳ — Ѳ0 |«Ѳо.

Вычитая из полученной при этом системы уравнения для невозмущенного процесса, приходим к искомой системе уравнений мезометеорологии (при записи ее отбросим штрихи)

ди

,

дѵ

,

dw

.(6.6.24)

дх

1

Оу

1

dz

 

 

t

 

 

(

P

A R

(6.6.27)

P

150

г д е

g - f l Т — const,

ий

- J r (Та t )> T=

c/z ’

 

cr

плавучести,

(6.6.28)

—const—параметр

-h - f(Ä -,V )

д .

dz'

dz

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

соты

процесса пренебрегаем

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

(а =

0) в уравнении (6.6.24);

вместо третьего уравнения движе­

ния используем уравнение статики; не будем учитывать боковой

турбулентный обмен; предположим, что в

невозмущенном сос­

тоянии атмосфера находилась в состоянии

покоя (П=Ѵ’ = 0 );

не будем учитывать силу Кориолиса; введем понятие среднего коэффициента турбулентности k = const. Если учесть сделанные

допущения и направить ось у

вдоль берега, а ось .ѵ по нормали

к нему (начало координат на поверхности

земли),

тогда сие»

теіѵіа уравнений примет следующий вид:

 

 

ди

Ои .

ди

ди ,

д2и

(6.6.29)

Tt

U — -J-'ZC

dz

T x ' ^ k ~dé'

дх

 

6 Я

W

ГІЙ

 

 

(6.6.30)

dt

-----SU' — ko.t -^-T.

 

dz

 

dz-'

 

 

dz.

 

 

(6.6.31)

 

dz

 

 

 

 

 

 

 

du M d-E.

- - 0.

 

(6.6.32)

 

öx 1

dz

 

 

 

Приступим к формулировке пограничных и начальных усло­

вий: при г 0

 

 

 

 

 

 

H = ffii = 0 .

 

(6.6.33)

151

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

Z—О $= (а0 + хаі) sin о»/,

(6.6.34)

С>о

где а0 и cti — заданные константы. За ал бугдем принимать мак­ симальную разность температуры суши и моря, деленную на ха­ рактерную длину явления. Очевидно, что условие (6.6.34) прав­ доподобно только в достаточно малой области вблизи берега (вероятно, в нескольких километрах на обе стороны от берега).

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

ы = 0 = я = О.

(6.6.35)

Вследствие локального характера явления при х = ± с о

■и = Ь= я = 0.

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

С учетом (6.6.34) будем искать решение системы (6.6.32) в виде разложений в ряд:

ti = u(zi /),

6 '= 'Ö'o(2'i t)+ xßi(z, t ),

я = яо(г, t) + хяі(г, t).

(6.6.36)

решение

(6.6.29) —

(6.6.37)

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

На основании (6.6.32), (6.6.37) и (6.6,33)

w 5= 0.

(6.6.38)

Подставив (6.6.37) — (6.6.38) в (6.6.29) — (6.6.32), получим систему, в которой переменные не зависят от х

ди

_

 

д2и

 

Hi

~ ~

1

к l)zv

 

^ 4 . « » 1 ==Ая.£!в?

 

dt

1

1

1 dz*

 

£»,

 

d%

(6.6.39)

dt

 

dz2’

 

 

 

_____j,

(,

 

 

152

и которую нужно решать при следующих граничных условиях:

2 = 0, и = 0,

0 , = öi sin о/, Оо = йо sin с>/;

(6.6.40)

2 =оо

и = О'о —■0 1= ло—л 1= 0.

(6.6.41)

Перейдем к безразмерным переменным (обозначим их индек­ сом п) на основании следующих соотношений:

,

/ 2k

2 „, * =

А , a ,= «1 V

1

/ —

1

(•)

 

0)

и

X-rtj

 

f T k

 

to

\

/

— ■«„,

 

 

(0

 

'

l a *

 

/ 2

*

0„>

II ОоС

м2~

 

/

w

 

 

 

 

 

а0 =

«l2 >«, f~2k

&Qn'

)2

1I

0

 

О

1/

)

 

/~оь

V ~ ‘■іяі

(6.6.42)

2k

(2)'

В таком случае уравнения (6.6.33) примут вид (опустим значки п):

 

Ли

 

1

Л2 Я

 

 

 

 

0»,

1

d2 8,

 

 

(6.6.43)

 

dt

 

о

Z2

 

 

 

 

 

 

 

 

dz

-lM>

 

 

(6 .6.44)

 

 

 

 

 

 

 

du

 

'-l~r

1

d2ii'

 

(6.6.45)

 

~ët ~

 

2

dzv

 

 

 

 

 

 

dü0 ,

 

2

dz2'

 

( 6 6.46)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dz

'V

 

 

(6.6.47)

 

 

 

 

 

 

 

Граничные условия запишутся так:

 

 

 

з—О,

! = sin

t,

i40= f l 0-sin t,

u — 0;

(6.6.48)

 

z — oo,

ы = А0 = 6 і = ло= Яі = О.

(6.6.41)

Система уравнений

(6.6.43) — (6.6.47)

отражает,

хотя и

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

153

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

Решим систему (6.6.43) — (6.6.47) при граничных условиях (6.6.41) и (6.6.48). Для решения (6.6.43) перепишем граничное условие при z = 0 для От

= і (cos t i sin t) = /cos / + sin t.

(6.6.49)

Вследствие линейности уравнения действительная и мнимая часть полученного комплексного решения должны каждая в от­ дельности удовлетворять уравнению. Если это решение удовлет­ воряет .условию (6.6.49), то действительная часть удовлетворяет

(6.6.48), а мнимая fh = cos/.

Будем искать решение в виде

 

=

(6.6.50)

где а и ß — неизвестные

постоянные. Подставив

(6.6.50) в

(6.6.43) и (6.6.49), определим а и ß

 

ct2 = 2 ß;

ß = —г,

 

откуда

 

 

а= ± У —2 г =±(1 — г).

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

= i£l(1^

2 •e~'lt — ie 2 .е~ 1^

г'г ^ =

/e± r[cos(± 2 -j-£)~-

— «sin ( ±

z y t) \ ~ e :-2[i cos( ±

z -4■0 +

sin ( ± 2 -f£)|.

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

 

 

§i = e

sin (^ — г).

(6.6.51)

Теперь из (6.6.44) нетрудно определить .щ

, "г — Je

^(si,n t - cosz — cos t - sin z)~dz-\-c— sin

&— Z

 

 

£— Z

-)- COS 2 ) -|- C—

X —2 ~ (—

COS 2 -f sin 2 )-f- COS t— (sin 2

 

V 2

cos

t

+ c >

 

 

154

так как при г — со

лі —0 , то с = 0 и

 

 

 

 

 

_

__ г

ѵ

соз

t —- z

4

і 4

]•

(6.6.52)

 

 

 

^

 

-

 

 

Подставим (6.6.52) в (6.6.45); тогда

 

 

 

 

ди

1

д-и

I о

-

2

 

 

 

(6.6.53)

dt ~

9

б>2 2

Ч

' е

• cos

 

 

z ~ 4 )■

 

 

 

Будем искать решение в виде

 

 

 

 

 

 

 

 

и — М (г) cos t + N(z) sin t.

 

(6.6.54)

Если подставить (6.6.54) в (6.6.53)' и приравнять коэффициенты при cos t и sin/1, то для определения M(z) и N (г) получим сле­ дующую систему уравнений:

— М ---- -- y V "= + L ? - е " г sin

- г);

 

N — - I

 

е

cos (

т — г) ‘

Умножим

первое

на і и

сложим

со

вторым, обозначив

W= M + iN,

 

 

 

 

 

W" - f 2 ЛБ =

у 2 е

cos

-t sin

(6.6.55)

Решим сперва однородное уравнение W" + 2iW ==0, которому соответствует характеристическое уравнение s2 -f2 t = 0 с корнями s = ± V —2і = ± ( 1 — і). В таком случае решение однородного уравнения запишется следующим образом:

(1— i)z

— ''(1— І) г

.

(6.6.56)

Wü-==Cl-e

\-с,-е

Частное решение неоднородного уравнения (6.6.55) следует искать в виде

W\ — ze ~(А cos z+ B sin z ),

(6.6.57)

так как a + ib = — ( 1 і) является однократным корнем харак­ теристического уравнения; а — показатель степени у экспонен­ ты; b — коэффициент при аргументе sin и cos.

155

Подставив (6.6.57)

в (6.6.55) и приравняв коэффициенты при

z sin г, г cos г, sin г и

cos z, получим систему уравнений для оп­

ределения А и В

 

 

2Л+270 = 0; —2(Л + Д) = 1+І;

, ,

—2В + 2іЛ = 0; 2(В — Л ) = 1 — і.

 

Из этой системы независимыми являются только два уравнения

 

А А- В

 

 

и В -

А =

 

 

откуда находим

 

 

 

_1_

 

 

I

 

 

 

 

 

 

 

В:

 

 

 

 

 

 

2

;

т

 

 

 

 

 

 

 

 

 

и тогда

 

 

 

1

 

 

 

 

И',

z-e

 

 

 

г).

—fr ( cos z i sin

Общее решение уравнения

(6.6.55)

можно представить в ви­

де суммы Wo и W\

 

 

 

 

 

 

 

 

(1-/ )z ,

- ( 1 - 0 г

 

 

-i-(cos 2 ~f+ sin z).

W—cx-e

 

-f-c2e-

 

z-e

Для определения c\

и c2

используем

граничные условия

(6.6.41) и (6.6.48) и выражение

(6.6.54):

 

 

при

2 = 0

0 = М

(0 ) cos t + N ( 0 ) s in t:

при

2=

со

0 = A

f ( с о )

cos t - - А Д о э )

sin t,

откуда следует,

что /VI (0) —N(0) = А4(оо) = ЛД оо)=0 и

 

 

при 2 = 0

1Р = 0 , с1+ с2 = 0 ;

 

I ".

 

при 2 = ОО 1Г-~- 0, Г;

= 0 :

 

:>'А

-

 

 

с1= —Сг = 0 .

 

 

 

Таким образом,

 

 

 

 

 

 

 

 

\V = M --riN —

 

1

—г

 

. .

-к z-e

(cosz+г

sin г).

Разделив действительную и мнимую часть уравнения, най­ дем, что

1

М (z)= -^-2 -e^-cos 2 :

156

 

1

2

' s'n Z'

 

 

 

N (z) — ---- 2 Z‘e

 

 

С учетом выражений для

M(z)

и

N (z)

общее

решение

(6.6.54)

будет иметь следующий вид:

 

 

 

 

 

и ----IfZ-e "cos (t— z).

 

(6.6.58)

Аналогично предыдущему можно решить и

уравнения

(6.6.46)

и (6.6.47), однако, поскольку

решения

этих

уравнений

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

»о=а0-(? • sin (t —z)-f----— 6 • cos 2 (t — z)-j-

 

4 ]

2

 

cos 2 I t —z ----

) - Т - е~ ‘Ѵ 1 X

 

 

 

XCos2

 

(6.6.59)

a0‘ cos

-0 , 0 2 sin 2 1.

(6.6.60)

 

 

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

Если определять момент появления ветра у земли при воз­ никновении бриза из условия

ди

dz

= 0

2-0

 

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

157

VII. ОСНОВНЫЕ ПРИНЦИПЫ ЧИСЛЕННОГО ПРОГНОЗА

погоды

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

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

пзадача оказалась практически невыполнимой.

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

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

ваемость этих методов выше, чем обычных синоптических.

§ 1. Общая постановка задачи численного прогноза

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

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

du dt

 

Ti тз

 

II iti

 

-I

 

u

 

d r

*

dt

 

p

dx

1

■ d - ¥ - l u ,

 

p

ду

 

A R T dp

 

cpp

dt’

(7.1.1)

158

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