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

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

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

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

Добавлен: 19.10.2024

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

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

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

где

Txz' — Т*хг' :г’ > lyz’ - ' y z бтУ* '

■в1' ='& -}-бй.

Здесь для простоты предположим, что ßs = ßs< т- е- Sßs = 0. Учет вариации ößs новых трудностей в рассмотрение не вносит. Услб-1 виями на боковой границе примем

 

«* = 0, z/ = 0,

- ^ - = 0 на 2

(7.15)

и начальными данными — климатические значения —

 

 

и' = и0, V‘ =V0, tfs^ S o -

(7.16)

Умножим уравнение (7.13) на сопряженные функции и*,

v*, w*,

4- р* и

й'д, соответствующие

невозмущенной (климатической)

Рр г

задаче, и почленно сложим, а сопряженные уравнения (7.5) умножим

1

ог(т2

на u', ѵ’, w', — p" и

~— йй и также сложим, затем результаты

РРг

этих операций вычтем друг из друга и проинтегрируем поflo X Dt. Тогда аналогично предыдущему получим выражение

т

т

 

 

 

I dt I аф'Ь* dS — I dt j ~

{l'xz'U * + l'yz'V * ) - f ß s l l ^ S dS —

 

 

. p

 

 

 

r

u0ul + v0v*0

go2

S0u S0

(7.17)

рг

 

 

 

 

Используя выражение (7.12)

и распространяя интегрирование

по t до —оо, соотношение (7.17) запишем для отклонений:

 

г

т

 

 

 

I dt I a s6M* dS

I dt j[-^(ÖT^^+ÖT^y^+ßsSMsdS.

(7.18)

- C O S L P

Учитывая (7.18), прогностическое выражение (6.14) можно пере­ писать в виде

б (рд£)

m

Jdt

г

J dt jа ^ Ш * dS +

jI-■=- (öXxz'U* + ötx2>v*)+

 

+ ßs6M s* dS

(7.19)

Формулу (7.19) можно несколько упростить. В самом деле, поскольку отклонения температуры поверхностного слоя трения океана от

183


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

бй = бй. Пренебрегая далее отклонением динамики ветровых тече­ ний в океане будем иметь

Распространение пределов интегрирования от нуля до Т на интер­ вал (—оо, Т) не вносит существенных погрешностей, если исходный интервал так велик, что за рассматриваемый отрезок времени началь­ ные данные «затухают» и ими можно пренебречь.

Поскольку здесь рассматриваются вопросы долгосрочного про­ гноза погоды, то такая замена является вполне оправданной. Если величины й* и й-s заранее рассчитаны по климатическим данным, то первый интеграл в правой части (7.20) вычисляется просто. Что касается второго интеграла, то для него величина бй оказывается неизвестной. Ее нужно найти на основе совместного решения уравне­ ний динамики атмосферы и океана.

6.8. ТЕРМ ИЧЕСКОЕ В Л И Я Н И Е КОН ТИН ЕН ТОВ

НА ДО ЛГОСРО ЧН Ы Й ПРОГНОЗ АНОМ АЛИЙ ТЕМ П ЕРА ТУ РЫ

При решении задач прогноза погоды обычно на поверхности Земли ставят условие тепловой изоляции

= 0

при z = 0.

~дг

 

 

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

Итак, изучим процесс формирования температурного режима почвы под влиянием заданного потока тепла на поверхности. Пусть й с — относительная температура почвы в точке (х, у , z') и в момент

времени t, отнесенная к средней температуре атмосферы Т, т. е.

Т — Т

Ас

т

184


Рассмотрим уравнение притока тепла в почве

 

 

 

д®с

дЧс

(8. 1)

 

~ді

 

О

 

VlC b z '2

 

при условии

 

 

 

 

 

 

=

ßc (бс — 'ö')

при z' = О,

 

дЬр

= 0 при

л' = h.

( 8. 2)

-

dz

,

 

 

г

 

 

Ось координат z' здесь направлена от поверхности z = 0 вер­ тикально вниз. В уравнении (8.1) ѵіс — коэффициент теплопровод­

ности почвы, ßc в (8.2) — коэффициент теплопередачи, й — темпе­ ратура поверхности Земли.

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

температуры

'&с = 'б'со

при t — 0.

(8.3)

 

Введем далее в рассмотрение сопряженное уравнение

 

 

dt

VlC

92Ч

= 0

(8.4)

при условии

 

 

d z '2

 

 

 

 

 

 

 

 

dz

№ с ~ Г с

при

z = 0,

(8.5)

 

дК

Л

при

= h,

(8. 6)

где

- з Г = 0

 

 

 

 

 

у) 6 С, tm

 

ас$*, если (х,

(8.7)

Гс = 0,

вне

этой

области.

Здесь й* — сопряженная функция, являющаяся решением задачи динамики атмосферы, а С — область, занятая континентами.

Начальные данные для сопряженной задачи возьмем в виде

 

 

 

йс = 0 при t = T.

(8.8)

Умножим уравнение (8.1) скалярно на йс , а уравнение (8.4) —

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

 

г

 

 

(8.9)

! dt

I Ж

dD

 

О

Г>с

 

 

Здесь Dc — область D, занятая континентами. Используя гранич­ ные условия и начальные данные функции (8.9)^преобразуем к виду

 

т

J (ссг йй£ - /*й) dS = 0.

 

- } (йсй£)0<Ш -

vlC J dt

(8.10)

De

О С

 

 

185


С помощью соотношения (8.10) получим

 

I

 

j

('ö'c'ö'c)odD -j- j

I

J

dt I ас$с$* dS =

dt J ßcM c dS. (8.11)

n

C

1C D C

0

C

Если временной интервал 0 ^

t ^ Т сделать достаточно большим,

то роль начальных данных будет несущественной и вместо соотнощения (8.11) получим

 

т

т

 

 

\ d t \ ас®с$* dS =

\ d t \ ßcM £ dS.

(8.12)

tm

C

-o o

C

 

На основе этого соотношения можно получить формулу теории

дозмущений (отклонений от климата)

 

 

т

т

 

 

j

dt / acöftcü* dS =

J dt j ßcöMZ dS.

(8.13)

t „

C

-oo

C

 

В дальнейшем будем считать, что б^с —

бй1.

погоды мы

При рассмотрении динамики атмосферы и прогноза

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

где — температура поверхности почвы, то, очевидно, в построен­ ных формулах для прогноза средних аномалий температуры конти­ нент будет отличаться от океана только коэффициентом теплопере­ дачи: на поверхности океана его следует взять равным ßs , а на поверхности континентов ßc. В результате основная прогности­ ческая формула (7.20) может быть обобщена и примет вид

 

'пг

f6ft(<xs$*+ßs®s)dS+

m

+ßc^c)dS

iS

J dt

J di j6^(ac&*

 

 

 

т

 

 

 

(8.14)

 

~

j

dt f f

ßs M*s dS-h f ßcöW cds)

 

 

tm

[ s

c

J

 

 

 

Выражение, стоящее

в квадратных

скобках в (8.14),

известно

и рассчитывается по фактическим отклонениям температуры от климатической 60 как на поверхности океана S , так и на поверхности континента С. Что касается второго слагаемого в правой части (8.14), то оно неизвестно и должно быть найдено в результате совме­ стного решения задачи динамики атмосферы и океана и теплопро­ водности почвы. Этот вопрос будет обсужден в п. 6.9 и 6.10.

Если предположить, что роль последнего слагаемого в правой части (8.14) мала, то, так же как при выводе (7.21), приходим к про­

186


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

б (рйj°)= q J dt fJö'ö' (asü* -fßs^g)

dS-\- J бй (асй* +ßc^c) dS (8.15)

-oo U

c

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

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

Врамках теории малых возмущений в данном случае всюду в фор­

мулах для средних аномалий температуры величины ßgöH и ßcöft нужно заменить на б (ай) и б (асй) соответственно. Например, фор­ мула (8.15) примет вид:

иш

°) =q J dt { J [б (а5й) й* + 6 (ßs0)d£]ÄS+J [6 (асй)й*

-оо

S

_

С

 

 

бфсй)й£] dS Г

(8.16)

 

 

6.9. Д РУ ГА Я Ф О РМ УЛИ РОВКА ЗА ДА ЧИ

 

 

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

При рассмотрении

задачи долгосрочного

прогноза погоды мы

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

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

Итак, рассмотрим систему основных уравнений:

J Au — гру + р Ц - — цр Ан = 0,

^ - + A v + lp u + P — РР Дн = 0,

187