Файл: Марчук, Г. И. Численное решение задач динамики атмосферы и океана.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