Файл: Марчук, Г. И. Численное решение задач динамики атмосферы и океана.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 19.10.2024
Просмотров: 138
Скачиваний: 1
возмущений уже оказывается недостаточной. Поэтому мы приходим к необходимости создания более полной теории, не связанной с пред положением о малых возмущениях. С этой целью в качестве основ ного «невозмущенного» состояния атмосферы будем рассматривать климатическое состояние. Реальное состояние атмосферы будем называть возмущением. Пусть невозмущенное состояние атмосферы описывается основной задачей (3.4)
B l n + A(P==f' |
|
|
Вф = Вф0 |
при t = 0 |
(5.1) |
и сопряженной (4.1) |
|
|
- В ^ ~ |
М*<р* = °, |
|
Вф*=Вфу |
при t = T. |
(5.2) |
В фазовом пространстве D X Т введем в рассмотрение скалярное произведение
5 |
|
т |
|
|
(g, h)DXT = 2) |
J |
dt J gihi dD. |
(5.3) |
|
‘_1 0 |
|
D |
|
Умножим далее скалярно уравнение (5.1) на ф*, а уравнение (5.2) на ф и результаты вычтем друг из друга. Тогда, используя начальные данные и условия, связывающие компоненты ф и ф* на границах
области D*, |
приходим к соотношению |
|
|
|
|
|
г |
|
|
(Вфг , |
Фт)д — (Яф0, фоЬ f { *[(ф*, ^ф)о — (ф, ^*ф*)о1 = |
|||
|
|
О |
|
|
|
|
= Jг (/, <p*)Ddt. |
|
(5.4) |
|
|
о |
|
|
Учитывая (3.9), получим |
|
|
||
|
т |
|
|
|
|
} [(ф*, Ац>)в — (ф, П*ф*)д]Л = 0 |
(5.5) |
||
|
О |
|
|
|
Тогда соотношение (5.4) несколько упростится: |
|
|||
|
(Вфт , |
Фг)н —(Вф0, фИ)л = (/, |
Ф*Ь- |
(5.6) |
Итак, операторы |
решения и входные |
данные для |
основных |
и сопряженных уравнений в невозмущенном состоянии связаны зависимостями (5.4)—(5.6).
Наряду с невозмущенным состоянием теперь рассмотрим воз
мущенное состояние атмосферы. |
Пусть оно описывается |
задачей |
в % + л у = г , |
|
|
Вц>‘ = Вфо |
при t = 0. |
(5.7) |
174
К задаче (5.7) присоединим сопряженную задачу, соответствующую невозмущенному состоянию атмосферы
~ В ^ - + А*ф* = 0, |
|
7?Ф* = By? при t = Т. |
(5.8) |
Как и прежде, скалярно умножим уравнение (5.7) на ф*, уравнение (5.8) на <р' и результаты вычтем. Тогда приходим к соотношению, аналогичному (5.4)
т |
|
(Яфт, Фг)о-(5фо, <PÜb+ j dt[(ф*, Л'ф'Ь —(ф', Л*ф*)л] = |
|
О |
|
т |
|
= } (f,<p*)Ddt + R, |
(5.9) |
о |
|
где В — некоторый функционал, связанный с тем, что граничные условия для компонентов решения ф могут оказаться неоднородными.
Вид этого функционала для |
конкретных |
случаев будет |
приведен |
в дальнейшем. |
|
|
|
Введем теперь обозначения |
|
|
|
А‘ = А + 8Л, |
ф' = ф+ 6ф, |
f —f + б/, |
(5.10) |
где А, ф и / — оператор и векторы, соответствующие невозмущен
ному состоянию. |
воспользуемся соотношением |
(5.4) |
Подставим (5.10) в (5.9) и |
||
и тем очевидным фактом, что |
т |
|
т |
|
|
J * ( ф*, А<р')п= |
j dt(cp\A*<p*)D + R. |
(5.11) |
о |
о |
|
Тогда приходим к формуле возмущений в виде |
|
|
(В5фт, фг)в — |
{ (б/, Ф *Ь *+ -бя. |
|
|
|
(5.12) |
Формула (5.12) будет основной для получения различных про гностических выражений для искомых функционалов задач.
6.6. ТЕО РИ Я ВОЗМ УЩ ЕН Д Л Я ЗА ДА Ч ПРОГНО ЗА ПОГОДЫ
Переходим к покомпонентной записи формул теории возмущений. С этой целью рассмотрим возмущенную систему уравнений атмосфер ных движений
д ри' |
ЛѴ — Ipv’— р |
|
— pp Au' = 0, |
dt |
|
||
|
|
|
|
дрѵ' |
ЛѴ — Zpu' — р |
дф' |
pp Au' = 0, |
dt |
ay |
175
|
g p ^ - ~ P ~ - o , |
|
|
|||
|
dpи' |
ору |
. |
opw |
|
|
|
дх |
dy |
|
dz |
|
|
öpö |
. Л ' д ' + |
^ р ^ - ^ - р ѵ , - ^ — рціА <К = 0 |
(6. 1) |
|||
щ |
||||||
с граничными условиями |
|
|
|
|
|
|
— |
= ag (# '— 6''), . |
pip' = 0 |
при 2 = 0, |
|
||
|
db' |
pw' = 0 |
при |
-Н |
(6.2) |
|
|
— = 0, |
|||||
|
dz |
|
|
|
|
|
и условиями периодичности решения по (х, у). Здесь ф' = О -f бй,
Ь' = б1 б$, O' и ф — климатические температуры соответственно воздуха на уровне будки и поверхностного слоя трения океана,
а 6# и бі} — отклонения от климатических значений. В качестве начальных данных примем
u' = uq, ѵ' = ѵ'о, ■&'=$'() при t —0. |
(6.3) |
Далее, рассмотрим сопряженную систему, соответствующую не возмущенному состоянию атмосферы
|
-----Au* + l9 v * ~ P ~ - W Ли* = 0, |
|
|||
|
Лѵ* — Іри* — р |
---- рр Аѵ* = 0, |
|
||
|
— |
р |
= О, |
|
|
|
дри* . âpv* . <9pu>* |
|
|
||
|
dx ' ây |
' |
dz |
|
|
öp#* |
— Лб'* — Уд- |
|
д$* |
рр, Дй* = 0 |
(6.4) |
dt |
Т |
|
dz |
|
|
при граничных условиях
â§* |
ОдФ*, |
pw* = 0 при z — О, |
|
dz |
|
||
о, |
|
|
|
д&* |
рір* = 0 при z= H , |
(6.5) |
|
dz |
а также в предположении периодичности решения по (х , у) и началь ных данных
и* = и*т, ѵ* — ѵ*т, •ö'*='ö'p при t = Т, |
(6.6) |
176
где u£, Vj и 'öy — функции, которые будут определены в даль нейшем.
Прежде чем переходить к конструкции формул теории возмуще ний, введем в рассмотрение следующие обозначения:
Л = А -j- 6А , ctg = ctg -]- öctg.
В соответствии с рассмотренной формальной процедурой основное
уравнение |
|
системы (6.1) |
последовательно |
умножим на и*, |
ѵ*, |
w*, |
|
— |
gT |
Ф*, затем |
полученные выражения сложим и |
вычтем |
|||
RTq>* и ■ |
|
||||||
аналогичное выражение, |
полученное умножением уравнений сопря- |
||||||
женной системы (6.4) соответственно на и |
ѵ', w', R T cp' и |
gT |
^ Ф' |
||||
_ |
с последующим сложением. Результат проинтегрируем по всему фазовому объему D X Т . Тогда с помощью простых преобразований приходим к выражению
j* |
~Т V'j'Vj’ —I-------- Ф р dD — j" ^HqHq—j—иqVQ■ |
||||
gT |
P dD-\- j |
Л j |
(и*8Аи' -f"v*8Av’ |
gT |
AHSAO' ] dD |
Уа~У |
0 |
D |
|
У a |
Y |
|
|
|
|
||
|
|
T |
|
|
(6.7) |
|
- q I dt j |
(otgfP — 6авФ') 0* dS —0. |
оs
В(6.7) предположим, что возмущения операторов и решений отсутствуют. Тогда получим формулу
J (и тііт-f- ѵтѵт+ |
“ |
T^r'^r'ö'r) PdD — j [ u 0u l -f п0Уо + |
|
||
+ -E l— ®0® l) p d D - g |
\dt f a sM *dS = 0. |
(6.8) |
|||
Уа У |
/ |
|
оJ |
sJ |
|
Учитывая обозначения |
|
|
|
|
|
и' — и -j~ 6м, |
vr = v-\ -6v, |
Ф' = Ф-: 8Ф |
|
||
и вычитая из (6.7) соотношение (6.8), приходим к формуле |
|
||||
j ^ 6ит• и*т-А ÖHy • ѵ*т ^ |
|
|
pdD — ^ ^öu0uj -f- 6k0• |
-r |
|
|
|
|
|
D |
|
+ y^—y |
) P ^ '' j ^ J ( и*ЬАи' + v*8Av’ |
|
|||
|
|
0 |
D |
|
|
7~ ГЙ*6ЛЙ') pdD — q j dt J [otgöfl-j |
8as (tf' — &')]&* dS = 0. |
(6.9) |
|||
|
0 |
s |
|
|
|
12 Заказ 674 |
177 |
Формулу теории возмущений перепишем в окончательной форме:
I |
^6ит• uj -'г 8ѵт• V*’-j- уд—у |
‘ |
Р dD = I ^6mq• и* -f- 8v0-н£-f- |
D |
“ |
|
D |
|
|
|
T |
|
p dD — I dt I ^и*8Аи' -f- v*8Av’ + |
||
|
__ |
0 |
D |
|
т |
|
-j— й*бЛГ) рсШ + д f dt f [asöä + öasflF' —#')]#* dS. (6.10)
Ya-Y |
1 |
J J |
Если теперь в качестве начальных условий для сопряженных уравнений возьмем условие (3.11), то приходим к формуле для ано малий средней температуры
8 (рф£) = j |
(8и0-11,1 + 8 0 0 ^ 1 + YagI ?-S V ft* ) РdD — |
|
||
Т |
D |
— |
|
|
— Jdt ^ (и*8Аиг+ ѵ*8Аѵ‘-|--- ^ — б^бАй'^ pdD ^ |
|
|||
0 |
D |
|
т |
|
|
|
|
|
|
|
+ |
q j |
dt j [as66 -f 8as (O' — O')] 0* dS. |
(6.11) |
оs
Проведем некоторые упрощения формулы (6.11) в том случае, когда срок прогноза достаточно велик (например, месяц). В этом случае, как было отмечено раньше, влияние начальных данных будет мало, поскольку интенсивность величин и*, ѵ* и О* со временем за счет диссипации будет уменьшаться. Тогда формула (6.11) при
нимает вид |
т |
_ |
|
__ |
|
||
б (рО?) = - |
f d t\ (и*8Аи' f v*8Av’+ —Ц— 0*6Л0') dD + |
||
|
-со D |
|
|
|
T |
|
|
|
~-q\ A |
J ]a560 + 6as (0' — O')]0*^5. |
(6-12) |
|
—co |
S |
|
Если предположить, что сопряженная задача решается при
фактических значениях и, |
ниш, то 6Л = 0 и мы приходим к формуле |
||
__ |
т |
|
|
8 (рб£) = q \ d t [ [as66 + 6as (б' - б')] б* dS. |
(6.13) |
||
|
- с о |
S |
|
Смысл этой формулы весьма прозрачен, а именно первый член в пра вой части (6.13) описывает вклад в аномалию температуры за счет взаимодействия атмосферы и океана, при этом член
г
q j dt j as6M* dS
-c o |
£> |
178