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

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

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

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

Добавлен: 19.10.2024

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

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

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

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

qJгdt J6as (ft' — $') $* dS

-с о S

описывает эффекты, связанные с отклонением теплопередачи атмо­ сфера — океан за счет штормов, нестандартной динамики ледового покрова и т. д.

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

Наконец, если пренебречь штормами и считать динамику льдов известной, то формула для прогноза аномалии температуры при­

нимает простейший вид

г

 

__

 

8 (рФу) = q J dt j as6M* dS.

(6.14)

 

-со S

 

В заключение следует отметить, что в настоящем параграфе мы рассмотрели только случай, когда источники в сопряженных уравне­ ниях задаются начальными данными (3.11). Все рассуждения остаются в силе для случая задания источников в граничных усло­ виях в форме (3.13). В этом случае следует лишь G заменить на Go и в качестве сопряженных функций взять решение задачи (3.7), (ЗЛЗ), (3.14).

6.7. ПЕРЕНОС ТЕП В БА РО КЛИ Н Н О М О КЕА Н Е

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

целесообразным поставить задачу о формировании теплового режима

в океане под действием процессов, протекающих в атмосфере. С этой

целью будем предполагать, что нам известны в любые моменты вре­

мени, предшествующие данному, температура

(или поток тепла)

на поверхности океана и напряжение ветра.

Выбирая в качестве

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

Итак, рассмотрим систему уравнений динамики океана в линей­ ном приближении:

Р

—(vsu* .)*'-^ Au = 0,

 

- J - + І и 1 - у

- (V s tv )z ' - VS]A V = ° .

12*

179


1

dp

go # s = 0,

 

 

 

p

dz'

 

 

 

P

 

 

 

 

 

du

.

du .

dw _л

 

 

 

dx

'

dy '

dz'

 

 

 

+ V W ~

(V iS®Sz’)z'—

HiS

=

0,

(7.1)

где z’ = —z, üs — относительная

температура

воды,

связанная

со стандартной температурой атмосферы соотношением

 

Т — уже ранее введенный параметр атмосферы. Такое определение удобно по той причине, что оно дает возможность сохранить непре­ рывность относительных температур атмосферы и океана на поверх­ ности океана: а — заданная константа, связывающая температуру воды и ее плотность по формуле р = cr&s- Влиянием солености мы пренебрегаем, хотя ее учет был бы тривиален.

Для удобства в системе (7.1) мы использовали те же обозначения для коэффициентов турбулентного обмена, что и для атмосферы, снабдив их только индексами S. Граничные условия для системы (7.1) выберем в виде

VSUz' = —

vszv = - ^ ß - t

w = 0,

 

 

Р

 

 

Р

 

 

 

^ p r = ßs (^s — Ä)

при

z' = 0,

 

 

u = 0, y = 0,

w — 0,

'

dzf = 0 при

z* = # s,

(7.2)

где й — заданная температура

поверхностного

слоя океана, совпа­

дающая с температурой воздуха,

Hs — глубина

океана,

ßs —ко­

эффициент теплопередачи океана.

 

 

 

 

 

На береговой цилиндрической поверхности будем иметь

 

и = 0,

у--0,

- ^ = 0

на

 

 

(7.3)

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

теристики

v = vQ,

®s = ®so

 

t = 0.

 

и = и0,

при

(7.4)

Далее, введем в рассмотрение систему сопряженных уравнений:

ди*

lv* -

1

dp*

(vsuZ')z' —Ps Лм* = 0,

dt

p

dx

 

 

 

 

 

 

du*

 

1

dp*

(vsy*')2' — Hs Ау* = 0,

dt

 

!>

dy

 

 

180



 

1

dp*

go

 

 

 

p

dz'

 

 

 

 

du*

du*

âw*

= o,

 

dil*

äx

ây

d z’

 

 

■W* - ' (ViS^*Sz’)z' -

f^is A'ö’s = 0-

(7.5)

dt

Граничные условия для сопряженных уравнений сформулируем следующим образом:

du* = o,

Ms

dz' Ps®*s

d$*s _

dz' ßsOs,

и* = 0, ѵ* — 0,

dv*

II

если t

u>* = 0,

o'

o' II* â

 

 

если

tm^ t

z’ = 0,

 

trn

 

 

 

дЬ*

при z’ = H s.

(7.6)

-^т- = 0

Здесь a sd* — заданная функция, которая находится из решения сопряженной задачи (3.7), (3.13), (3.14) для атмосферы. На береговой цилиндрической поверхности поставим условия

и* = 0, к* = 0, - ^ - = 0 на 2 -

(7.7)

Начальными данными для системы (7.5) возьмем

ц*-=Ку, ѵ* = ѵт,

(7.8)

где uj , и? и ö’sj. — пока произвольные функции координат. Естественно, что мы всюду предполагаем выполнение условий

гладкости решений основной и сопряженной задач, аналогично тому, как мы это имели при рассмотрении уравнений динамики атмосферы.

Умножим далее уравнения системы (7.1) соответственно на и*, ѵ*, w*, L р* и ■dg, затем их сложим. После этого умножим

Ррг

Уравнения сопряженной системы (7.5) соответственно на и, v, w, L p и ■Йо и также сложим. Затем результаты вычтем один из

РрГ

Другого и окончательное выражение проинтегрируем по всему фазо­ вому пространству Do X Т, где Do — область, занятая океаном. Тогда с учетом граничных условий получим

go2 A A *

( (u0u*0-\

-

V S ^ S T ) d D —

pr

 

 

рг

T-</z' „ * ___ P±m*s) d s - L

 

r

]

r

p

181


где

 

 

 

 

 

 

рГѵ IS

 

при

z' = 0.

 

 

 

 

 

 

 

 

 

go"

 

 

 

Отсюда следует

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J dt j as®®* dS = r

f I Itj^UjT+ vTv*T+ —— 'fl's'ösy J dD

tm

S

 

 

 

 

 

 

 

 

 

 

 

 

рГ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J (u0u*0 + v0v*0+

 

Osodsoj dD j +

 

 

 

0

 

s

 

 

 

 

 

 

>*)

ßsO^sJ dS.

(7.10)

Пусть

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(7.11)

 

 

 

 

 

 

 

 

 

 

 

 

® s t 0.

Тогда формула (7.10) несколько упростится

 

 

г

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

J dt ^as®®*dS— J dt j

 

 

(txz'u* -j- Xyz'v*) + ßsM i

dS

*m

S

 

 

 

о

s

*- p

 

 

 

 

 

 

 

 

 

r

j (u au*0 +

v 0v%+

jr^ - ®so® so j d D .

(7.12)

Введем теперь в

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

возмущенную систему

основных

уравнений

du'

 

 

 

 

 

dp'

 

 

 

 

 

 

 

 

 

l v f

i

1

(У s u z')z'’

 

 

 

dt

 

 

 

 

p

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ді/

t

о /

;

1

dp'

■(Vs^')z'

1

 

 

 

dt

 

 

 

 

P

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

dp'

 

 

 

— 0,

 

 

 

 

 

 

 

 

P

dz'

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

du'

dv’ -1... dw'

- o

 

 

 

 

 

 

 

 

dx

'

dy

 

'

d z’

 

 

 

 

 

8

,

Г

 

in' — (v1s'6's2')z—

Pis

— О

(7.13)

при условии

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ѵ 8и'г, =

— >

 

VSVz’

 

^yz’

 

in' =

0

 

 

 

 

1

 

 

дй’

 

р

 

 

 

 

 

 

 

P

 

 

 

при z' = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

's

P s ^ s -# * ),

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

u* = 0,

v‘ = 0,

 

117'=0,

 

-

^

= 0

при z’ = H8, (7.14)

182