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

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

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

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

Добавлен: 19.10.2024

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

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

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

где

 

 

 

д ( -

дТ' \ ,

д [ - . „ д Т '

М

/т = а2 sin d sin d д Х

д Х

öd (РіР8Ш^ öd

)J '

Ha втором шаге решаем задачу

 

 

 

öpu

ctg Фру = О,

 

~дГ

 

 

 

 

дрѵ

ctg Фри == 0,

(16.4)

 

dt

ана третьем шаге — систему уравнении адаптации

дри

2(0 cosd Of

 

----- т—5г 4?-»

 

a t

 

r

 

a sin О

 

öpy

-2(0 cos Фри =

J_

дф

 

 

öj

 

 

 

 

 

a

ö d

 

 

ö(p

 

g

Г ,

 

 

 

 

 

ÖZ

Д y 2

 

 

 

 

 

 

 

 

 

 

 

div pu = О,

 

 

 

 

0 р Г

,

ч -

 

а

-

ат'

(16.5)

a t

(уа- у ) р и ; = —

ѵ,р

dz

Системы уравнений (16.3), (16.5) решаются методами, рассмотрен­ ными выше. Новым является уравнение (16.4). Решение системы уравнений (16.4) на интервале tf sg t + 1 найдем в аналитиче­ ской форме:

 

u

=

v! sin b (t tj) + и' cos b (t tj),

 

 

v

=

vi cos b (t tj) — u' sin b (t tj) .

 

Полагая t =

t j + 1 и

tj + 1 t} = т, получим:

 

 

 

 

ui*1= ui cos Ьт+ vi sin bx,

 

 

 

 

v i + 1 =

v i c o

s b x u i

sin b x .

(16.6)

Ввиду того, что |sin b x \

«gl,

|cos b x \

«gl, решение

(16.6) будет

ограниченным

и корректность

решения задачи (16.6)

обеспечена.

 

 

 

7.17.

ОСНОВНЫ Е И С О П РЯ Ж Е Н Н Ы Е У РА В Н Е Н И Я

 

 

 

 

 

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

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

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

19 Заказ 674

289


Для атмосферы

öpu ' - div puu — ^2(ocos Ф

ctg

pu

fn

1

t?cp

Fu,

а sin 0

dX

 

 

dpи

, div puf -r (2(ö cos d H— ctg ft ) pu =

i_

dtp

 

1FV 1

 

d t

 

&

J

а

д&

 

 

 

 

<?ф

дуг

Г ,

 

 

 

 

 

 

dz

 

 

 

 

 

 

div pu = 0

 

 

 

 

 

 

f divриГ + {ya- y ) p w = - ^ + FT;

 

С17-1)

для

океана

 

 

 

 

 

 

 

дрѵ

div puu — ^2(ocos d + -^- ctgft^ рг; = —

1___

др'

Fu,

 

dt

sin Ф

дХ

 

дрѵ

div риг; + (2® cos ft + ^-ctg ft) pu =

 

+ Fv,

 

dt

 

 

 

 

 

 

 

 

 

 

 

d p '

-oT’,

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divpu = О,

 

 

 

 

 

 

- 4 div puT' + Гргн = FT;

 

 

(17.2)

для

континента

 

 

 

 

 

 

 

 

i Pd tl ^

- . F

r .

 

 

 

(17.3)

Следует иметь в виду, что коэффициенты ѵ, ѵх, р и px в выраже­

ниях Fu, Fu и для атмосферы,

океана и континентов берутся

в зависимости от своей области определения и, кроме того, рг для почвы принимается равным нулю.

К системам уравнений (17.1)—(17.3) необходимо присоединить граничные условия: при z = Нт

ди

дѵ __~

дТ'

= 0 , w = 0;

(17.4)

dz

' dz

dz

 

 

при z = —h*

du

~

du

d l" = 0, tv = 0;

~dz ~

U ’

~ d z ~ V ’

~ d z ~

при z О на поверхности океана

= F,

(17.5)

iu = 0 ,

(17.6)

290



где символом [аі обозначен разрыв функции а, т. е.

[а] = lim |г=+6 — а |г=_Е};

£->-0

при г = 0 на поверхности континента

и = О, V= 0, w' = О,

-

д Т '

(17.7)

 

F .

Сопряженная задача, соответствующая задачам (17.1)—(17.7), имеет вид.

Для атмосферы

-^ Г + div P™* - (2« cos в + 1 ctg fl) pv* = -

öpü • + div рик* + ^2ö) cos ft + -j

ctg ft^ pu* = — -j

dt

 

d(f*

8

 

dz

I— 71*,

 

 

 

div pu* = 0,

.d p T*

div риГ* -f (Уа— Y) Pw* = Ft*-

dt

 

 

F u*.

----Fa*,

(17.8)

где Fu*, Fv* и ^ * берутся в форме (16.2), где вместо и, ѵш Т следует подставить и*, ѵ* и Т*. Что касается коэффициентов р и plt то они по-прежнему вычисляются по формуле (16.2').

Для океана

Т

+

div pun* - (2® cos ft + 1

ctg ft) pv* = -

-

F„.,

—f p

+ div pun* -f ^2(0 cos ft +

ctg ft) pu* = - j - %

F°*'

 

 

 

dp*

oT*,

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

div pu* = О,

 

 

 

 

dpT*

+ div риГ* + Грю* =■ — F t * ’,

 

(17.9)

 

 

dt

 

 

 

 

для

континента

 

 

 

 

 

 

 

dpT*

 

 

(17.10)

 

 

 

dt

 

 

 

 

 

 

 

 

К системе сопряженных уравнений (17.7)—(17.9) присоединим следующие граничные условия:

при z — Нт

ди*

Q дѵ*

■0, - ^ = 0, w* = 0;

(17.11)

dz

dz

 

 

19*

291


при

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ди*

п

ди*

_

п

дТ *

= 0 ,

іо* =

0;

(17.12)

 

 

 

dz ~

U’

dz

~

U’

dz

 

 

 

 

 

 

 

при z

О на поверхности океана

 

 

 

 

ди*

=

,

 

 

 

 

 

дТ* ~

= F*,

w* = 0;

(17.13)

[ѴР ^ Г

 

на

поверхности

 

 

dz

 

 

 

при

z =

0

континентов

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

и * =

0,

ѵ* =

0 ,

ш* =

0, £ѵхр

] =

F * -

(17.14)

Начальными

данными являются

 

 

 

 

 

 

 

и = и°,

V= ѵ°,

Т' — Т'° при t = О,

 

 

 

 

 

и* —О,

V* = О, Т* 0

при t — t0

 

(17.15)

в атмосфере и

океане и

 

 

 

 

 

 

 

 

 

 

 

 

Т' —То

при

t = О,

 

 

 

 

 

 

 

 

 

Т*

= 0 при t t0

 

 

(17.16)

в почве.

Итак, задача поставлена. Следует лишь дать определение вели­

чине F*.

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

f«4;T / 2

 

 

^2

 

 

 

\ dt

J F*T‘dS = J

dt

j FT*dS.

(17.17)

—co

S+ C

- с о

 

S +

C

 

Если выбрать

 

 

 

 

 

 

F* =

.7 ,—

при T

г

t

T -j- xtг

 

 

ДСгТ

 

 

 

 

 

и F* = 0 вне этого интервала, где АG — область континента, для которой дается прогноз, тогда будем иметь

*0~\~ХІ г

\ d t j FT*dS.

- c a

S+ C

В отклонениях от климата имеем

/ 2

бTj.%, — j dt

J ÖFT* dS.

(17.18)

-00

s+c

 

Аналогичные формулы можно получить для аномалий вертикаль­ ных токов (см. гл. 6).

292