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