Файл: Марчук, Г. И. Численное решение задач динамики атмосферы и океана.pdf

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

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

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

Добавлен: 19.10.2024

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

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

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

О

при

z' = Ä,

 

dz

 

 

 

йс = йсо

при

t = 0,

(9.20)

где h — глубина почвы, на которой годовой ход температуры отсут­ ствует.

Сопряженную задачу по отношению к (9.20) определим следу­ ющим образом:

 

 

дЪ*с ѵ іс

дЩ*с

о,

 

 

 

 

dt

 

 

dz'2

 

 

 

 

йС* =

 

при

z' = 0,

 

 

 

= °

"Р*

z' = k,

 

 

 

■&С= 0

при

 

t > T .

(9.21)

Обычным способом приходим к функциональному уравнению

 

ГТ

г кг

J

г -

д&*

 

 

J йсойсо dD.

(9.22)

\ d> \ T S ? L » * d s = S d l l ® i r ds

Ѵ1С DC

 

о

с

о

с

 

 

 

 

Далее, полученные выражения (9.19) и (9.22) используем для исклю­ чения соответствующих выражений из основного соотношения (9.11). Тогда получим:

На основе (9.23) приходим к формуле малых возмущений. Для этой цели предположим, что в начальный момент времени начальное состояние океана и почвы климатические. Тогда, предполагая, что

192


истинное состояние системы мало отличается от климатического, получим формулу для средней аномалии температуры

6( < 1 *_) = - f ( Ч ^ + «ВД

ЧА*) Р dD -

 

 

 

(9.24)

Здесь,

по-прежнему, бй — заданное отклонение относительной

температуры при z

= 0 от климатического значения.

Если,

наконец,

интервал 0 sc t

Т достаточно велик и началь­

ные данные в атмосфере уже не играют существенной роли в прогнозе аномалий температуры, то формулу (9.24) можно заменить более простой

_ .

J

/ г - дй*

г - дй*

\

іЙ Ц '"*

\dtU^i^ds + l m -a/is W

\

2 '

-о о

\ S

С

/

+ 7

^ Т

J *

( \ 6F^ d D -

I

- (9.25)

 

-о о

\ D

S+ C

i

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

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

о,

то приближенно запишем

j ö F ^ d D - J ÖFhÜh dS = J ÖFfttdS

D

S+C

S+C

и, следовательно, формула (9.25) упростится:

б (рй^ Т

бй —

-

дЪ*

dS j »

с dS

2

 

 

dz'

 

 

 

8

Т

 

 

f dt [ бУ^й* dS.

(9.26)

Ср (уа—у)

»S+ C

Полезно отметить, что формула (9.26) могла бы быть получена строго в предположении, что лучевым притоком тепла в атмосфере31

13 Заказ 674

193


можно пренебречь, т. е. в уравнениях (9.1) с самого начала положено е = 0, а трансформация лучистой энергии в тепловую происходит только на поверхности океана и континентов, где выполняются гра­ ничные условия (9.8) и (9.9). Такое приближение может оказаться удовлетворительным.

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

ходятся выражения й*,

и

. Они могут быть затабулированы

заранее при фиксированных входных параметрах: Go времени года и интервалах осреднения Т — т ^ t sg Т . Затем рассчитываются величины

г -

дф*

 

Г -

a{t) = —q ] b b ~ - d S ,

b{t) = —q j 8ft

s

 

 

c

c (t) = q

г

f 6Ffi*dS.

1

4 MVa-Y)

J 0

 

 

s+c

Функции а (t), b (t) и c (t) рассчитываются по известным из наблю­

дений значениям бй и б /’о в интервале времени — < t < tm. Недостающие значения этих величин в интервале tm ^ t sg Т могут быть найдены, например, на основе теории экстраполяции случай­ ных функций, учитывая большую корреляционную связь метеоро­ логических полей во времени. Этот путь, по-видимому, является весьма рациональным. Однако простейшая формула прогноза полу­ чается в том случае, когда используются упрощенные формулы. В са­ мом деле, поскольку величины й^, й£ и й* со временем затухают (если, например, речь идет о прогнозе на месяц или сезон), то с не­ которой погрешностью вместо (9.26) можно воспользоваться следу­ ющим соотношением:

g

/»«•

 

J dt J бFfi*dS.

(9.27)

с р ( У а — У )

 

S4-C

 

 

6.10.

ДО ЛГОСРОЧНЫ

С РЕД Н И Х О Т К Л О Н Е Н И Й ОСАДКОВ ОТ НОРМ Ы

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

194


ную связь между вертикальной скоростью w и осадками. Такая связь указывает, что осадки, как правило, коррелируют в нижней тропо­

сфере с вертикальной скоростью w, обращаясь в нуль при

w <

0.

Если этот

факт принять за основу,

то можно построить теорию воз­

мущений для средней по региону Qh вертикальной

скорости

w,

осредненной также по интервалу времени

Т — т ^ К

f . Здесь

Qh — трехмерная область в тропосфере,

которая своей

проекцией

на плоскость z = 0 имеет Go.

 

задачу

 

 

 

 

 

 

 

 

Рассмотрим

сопряженную

 

 

 

 

 

 

 

 

 

 

 

-----Ли* 4- Ірѵ* р

----- pp Au* = 0,

 

 

 

 

 

 

~ Аѵ* * 1Ри* ~ Р

—W* Аѵ* = °*

 

 

 

 

 

 

 

 

gp#*

 

 

=

 

 

 

 

 

 

 

 

 

 

 

dpи*

 

дрѵ*

 

dpw*

 

 

 

 

 

 

 

 

 

 

 

dx

 

dy

 

dz

 

 

 

 

 

 

 

 

 

M ! . _ ЛФ* + J t

l

gw* -

(p\i»z)z-

P!p Лй* = 0

(10.1)

при условии

рц;*=0,

й* = 0

при

 

г — 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рw* = 0,

dz

0

при

Z — H

 

 

(10.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и начальных данных

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

u * = 0 ,

V* = 0,

й* = 0

при

t > T .

 

 

(10.3)

 

I* (х, у) и * (г) ц * (t)

при

(х, у, z) 6 Q a

и Т — x ^ t ^ T ,

 

 

(

 

f* — {

0 вне

этой

области,

 

 

 

 

 

 

 

 

 

 

1

I

 

 

 

 

 

 

 

 

 

 

a g* (X, у), я*

(z) и

г]* (t) — функции,

нормированные

в й ЛхД Т

на единицу (здесь АТ обозначает

интервал

Т

т

 

t ^

Т).

Уравнения системы (9.1)

умножим соответственно

по

и*,

ѵ*,

w* В Т ф* и еТ- й* и почленно сложим. Далее систему уравнений

т Уа - У

(10.1) умножим на и, ѵ, w, ВТ<р и

й и сложим. Результаты вы­

чтем друг из друга и проинтегрируем по всему фазовому простран­ ству, используя граничные условия (9.2), (10.2) и начальные данные (9.3), (10.3). Тогда приходим к функционалу

= [ ( „ о„; + ѵ Г Ь _ П - » р(» ;) р ід + , | *

J » % - d s . (іоА)

Т 2 D

0

S+C

13*

195


Для малых возмущений получим

6 (■ £ л ) = 1 ( 8иои*о+ 8ѵоѵо +

6W ) РdD +

т

+ q \ dt ^ ö ü - ^ d s . (10.5) 0 s+c

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

Если интервал времени 0 ^ t sg Т достаточно велик, то началь­ ными данными можно пренебречь и мы приходим к формуле

б

(10. 6)

На основе формулы (10.6) можно в рассматриваемом регионе в интервале времени установить превышение осадков над климатиче­

ским значением, если б / іѵйь \ > 0

и дефицит, если б / wah \ < 0.

\ т- т )

l r - f j

Формулу (10.6) упростим, учитывая, что й* затухает со временем, тогда имеем

w h т_) = ? \ d t \ w ^ d S .

т-

dz

S f C

6.11.

(10.7)

ПРОСТЕЙШ АЯ

ДОЛГОСРОЧНОГО ПРОГНО ЗА ПОГОДЫ

Переходим к расчету сопряженных функций. С этой целью пред­ положим, что основной нашей задачей является долгосрочный про­ гноз средней аномалии температуры в заданном регионе G0 и в ин­ тервале времени Т — т ^ t ^ Т. Поскольку пространственные и вре­ менные масштабы явлений, ответственные за такой прогноз, весьма велики, то естественно сосредоточить внимание только на наиболее крупных возмущениях фонового характера. Что касается масштабов возмущений циклонических размеров, то их будем описывать фено­ менологически с помощью системы макромасштабного турбулент­ ного обмена. В этом случае приходим к сопряженной задаче

dpи*

— дри*

f фу* — р-%

црДи* = 0,

dt

U дх

 

 

 

дри*

дрѵ*

Іри* р

<9ф*

др Аѵ* = 0,

dt

дх

W

196