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

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

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

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

Добавлен: 19.10.2024

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

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

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

Предположим, что нас интересует прогноз среднего поля темпе­ ратуры по области G {х, у £ а, Os ^ z s ^ k ) . Тогда «начальные» для сопряженных уравнений условия выберем в виде

iij* = 0,

= О,

 

——— 19г =~7^> если

X £ G и Фу = 0 вне G,

(3.11)

Ye - Y 1 G

 

ѵ

где X — совокупность координат (х , г/, z).

Введем в рассмотрение обозначение для средней по области G аномалии температуры в момент времени t = Т:

J рb TdD =рйу.

 

G

 

 

 

Тогда формулу (4.10) перепишем в виде

 

 

_

_

 

г

 

Р$г = j (и 0и*0 + ѵ0ѵ*0 уаеТу

tt0K ) p d D +

q j dt

J asM * dS. (3.12)

D

a

 

o s

 

Выражение

означает,

что средняя

аномалия температуры

рассчитывается по данным на интервале 0 ^

t ^

Т .

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

задача прогноза средней аномалии температуры

свелась к решению сопряженной задачи (3.7), (3.8) при условии

(3.11).

В настоящем параграфе была построена теория возмущений при специальном задании начальных условий для системы сопряженных уравнений и при однородных граничных условиях. Покажем, что возможны и другие постановки задач для системы сопряженных уравнений, которые приводят к формулам теории возмущений, удобным для практического использования. В самом деле, рассмо­ трим систему уравнений (3.1) вместе с граничными условиями (3.2) и начальными данными (3.3). Рассмотрим далее сопряженную систему уравнений (3.7).

Граничные условия для системы (3.7) определим следующим

образом:

 

 

- dz = agd*

i-/*, рw* = 0 при z = 0,

 

ЯА#

(3.13)

—^ -= ^0,

рц;* = 0 при z = H ,

где /* (X, у, t) зададим в форме

/* = б(t — Т), если (х, у) G0,

/* = 0, вне указанной области.

Здесь G0 — некоторая область на поверхности Земли, где требуется дать прогноз средней аномалии температуры.

170



В качестве начальных данных для сопряженной задачи возьмем

и* = 0, у* = 0, й* = 0 при любом t^>T.

(3.14)

С помощью основной и сопряженной задач и уже рассмотренной техники получим функционал

Р^т° = j

( и0и*0-і- vüv l Ч-

«

) p d D -f

 

 

q j" dt

j a sM * dS,

 

(3.15)

 

0

s

 

 

 

 

рйг° =

j P®TdS.

 

(3.16)

 

 

Go

 

 

 

Предположим, что мы имеем дело с малыми возмущениями и' =

= и -f би, у' = у +

6у,

б$\ й1' = й + бй. Тогда мето­

дами, изложенными выше, приходим к формуле

 

б ( і а д = } ( б и 0^ + 6 у0у0* +

 

бЙ 0Й0* ) р

+

D

 

 

ü

 

 

 

q JГdt

Ja s6M* dS.

(3.17)

 

о

s

 

 

 

Сравнивая формулу (3.15) с (3.12), видим, что они совпадают, если выбрать в (3.12) G = Go. В этом случае, как нетрудно убе­ диться, будут совпадать и решения сопряженных уравнений, по­ скольку безразлично, задать ли мгновенный «источник» при t = Т на Go в граничных условиях или задать его на Go в момент времени t — Т как условие Коши.

До сих пор мы предполагали, что задача состоит в рейіении долго­ срочного прогноза погоды на момент времени t = Т по отношению к начальному t = 0. Естественно, что при рассмотрении месячных или сезонных прогнозов нет смысла так точно определять момент предсказания аномалий температуры или других элементов. Мето­ дически было бы более правильно находить прогноз средней анома­ лии температуры за некоторый интервал времени. Например, при прогнозе погоды на ближайший месяц было бы целесообразно дать

его осредненным по декадам:

средний прогноз за первую, вторую

и третью декады. Если нас

интересует прогноз погоды на сезон,

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

171


С этой целью в качестве граничных условий для системы сопря­ женных уравнений выберем следующие:

 

ДА*

 

+

 

 

 

=

 

pw* = 0 при z = 0,

 

 

д

 

 

 

(3.18)

 

-^— = 0,

рш* = 0 при z = # ,

где /* (х, у,

t) задается в форме

 

 

 

/*(*, г/.

-7É*(*»

y)n*(t),

если (х, y)eG 0, t £ \T —т,

Т],

0 =

 

 

 

 

 

 

0

 

 

вне этой

области.

 

Относительно функций

I* (х , у) и rj* (t)

предположим только,

что они нормированы, т.

е.

 

т

 

 

 

 

 

 

 

 

 

§ l* (x ,y )d S = 1,

J* r]*(t)dt = l.

 

 

Go

 

 

T - т

 

 

Вчастном случае, естественно, они могут быть постоянными.

Вкачестве начальных данных для сопряженных уравнений примем

и*т= 0, [4 = 0, $4 = 0 при t> T .

(3.19)

При граничных условиях (3.18) и начальных данных (3.19) обыч­

ным методом получим функционал

 

 

р^ г -— ~

1 ( иоцо + ѵоѵ*«+

у * W ) РdB -f

 

2

D

Т

 

 

 

 

 

j «SM * d s ,

(3 .20)

 

+

Ч j j d T

где

 

0

s

 

 

 

T

 

 

 

_

 

 

 

 

P^r-2 L =

I ^ { ^ d t

Jp^*(a:, y)dS.

(3.21)

 

2

T—T

 

Go

 

Сравнивая формулу (3.20) с (3.15) и (3.12), видим, что формально они совпадают друг с другом, изменяется лишь смысл сопряженного решения. На практике, конечно, удобно выбирать функции и р* гладкими и положительными. В этом случае соотношение (3.21)

будет иметь смысл осредненного с весом значения

температуры

в заданной области Go и в интервале времени t — т ^

t ^ Т.

6.4.

М ЕТЕОРОЛОГИЧЕСКОЙ ИНФ ОРМ АЦИИ ПО ОТНОШ ЕНИЮ К АНОМ АЛИЯМ СРЕДНЕЙ ТЕМ П ЕРА ТУ РЫ И ДО ЛГОСРО ЧН Ы Й ПРОГНОЗ

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

172


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

- B f - M

V = o ,

(41)

Вц>* =

при t = T,

 

где компонентами вектора В ф£, например, выбраны условия (3.11). Нетрудно представить себе, что при t — Т, согласно условиям (3.11),

имеем Uj = 0, ѵ? = 0, а

отлична от нуля только в области G,

где она является константой.

Предполагая функции и,

v,

w задан­

ными, решаем систему сопряженных уравнений для t <;

Т.

При этом

за счет переноса субстанций область ненулевых начальных данных переместится на запад примерно на расстояние и At, на север (или юг) на V At и по z на w At, где At = Т t. Конечно, сразу же при этом возникнут гравитационные волны — волны Россби, которые «размажут» эту картину, расширив область возмущений и т. д. В результате турбулентной диффузии интенсивность полей компо­ нентов сопряженных функций постоянно будет уменьшаться, в пре­ деле при t -V — оо стремясь к нулю. Это значит, что при достаточно далеких прошлых моментах времени по отношению к моменту t = Т за счет диссипативных процессов информация о начальных полях уже не будет полезной, так как превращается в метеорологические шумы. В математической модели это и выражено тем, что первый интеграл в формуле (3.12) при t ->■ —°° будет стремиться к нулю, а вариации температурной аномалии определятся только потоком тепла из океана, т. е.

__ т

8 (р'б'г) = q J dt I a s6M* dS.

(4.2)

-с о S

 

Это предельное соотношение позволяет сделать весьма важное за­ ключение о роли океана в формировании аномалий температуры в долгосрочном прогнозе погоды. Более того, структура формулы (4.2) указывает на тот факт, что сопряженная функция ■&* — решение сопряженной задачи — является функцией влияния в отношении

кпрогнозируемым вариациям.

Всвязи с тем что сопряженное решение является в конечном

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

6.5. О БЩ АЯ ТЕО РИ Я ВОЗМ УЩ ЕНИЙ Д Л Я ЭВОЛЮ ЦИОННЫ Х ЗАДАЧ

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

173