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

Здесь у (X, у ) — заданный на поверхности океана поток плотности. Кроме того, мы предполагаем, что океан имеет постоянную глубину Н, совпадающую со средней толщиной слоя главного термоклина,

ао — цилиндрическая береговая поверхность.

Вкачестве одного из граничных условий в (2.5) при ъ = 0 мы выбрали условие заданного на поверхности океана потока плот­ ности. Такое условие выбрано только ради простоты и определен­ ности. Все последующие результаты легко обобщаются на тот слу­

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

р = р при z = О,

либо более сложное

= а (р р) при z = О,

где а и р — заданные функции координат и времени.

Докажем теперь единственность решения задачи (2.4), (2.5). С этой целью, как обычно, предположим существование по крайней мере двух различных решений системы уравнений (2.4), удовлетво­ ряющих одним и тем же граничным условиям. Составим разность таких решений. Тогда разность будет удовлетворять системе (2.4) и однородным граничным условиям (2.5). Построим теперь квадра­ тичный функционал. Для этого уравнения системы (2.4) умножим

С помощью формулы Гаусса—Остроградского и при учете гра­ ничных условий

врdz — 0, w — 0 при z = 0,

D

Отсюда, а также принимая во внимание граничные условия, следует, что для разностей компонентов решения имеют место следующие тождественные для всей области D равенства:

101

Что касается давления, то оно, как и плотность, определяется с точ­

Таким образом, единственность решения задачи (2.4), (2.5) в ука­ занной системе доказана.

Переходим к формулировке модели динамики с явным выделением баротропной составляющей. Для этой цели решение системы урав­

Соотношения (2.12) подставим в (2.4) и результат проинтегри­ руем по всей глубине океана, используя граничные условия (2.5). Тогда для баротропной составляющей решения получим систему

102

Система уравнений (2.14) незамкнутая, поскольку в нее входят величины р ' .

Для получения уравнений бароклинной составляющей решения выражения (2.12) снова подставим в систему (2.4) и для исключения баротропных составляющих течения и р0 воспользуемся получен­ ными уравнениями (2.14). Тогда приходим к следующей системе

больше, чем требуется. Однако следует помнить, что дополнитель­ ное условие для р' является естественным условием для данной постановки задачи.

Системе уравнений (2.16) придадим более удобный для решения вид. С этой целью с помощью уравнений статики и неразрывности

Ог

Спомощью соотношений (2.18) систему уравнений (2.16) преобразуем

к виду:

103

при условии

у при z = 0,

(2.20)

Кстати, необходимо отметить, что на основе специального пред­ ставления р в форме (2.12) следует, что

Этот вывод может иметь важное значение для рассмотрения других постановок задач динамики океана.

Далее рассмотрим задачу (2.14), (2.15). Она имеет только три­ виальное решение при любых jo'. Покажем это, предположив до­ статочную гладкость решения задачи и входных данных. С этой

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

где использовано обозначение

н

о

Второй интеграл в (2.22) преобразуем, записав в форме

а

Первый интеграл в правой части (2.23) равен нулю вследствие фор­ мулы Грина

а второй интеграл в (2.23) равен нулю вследствие уравнения нераз­ рывности

104

Из соотношения (2.24), а также из граничных условий на а сле­ дует, что

и = 0, V = 0,

а из уравнений (2.14)

р = р0 + р" = const.

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

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

Докажем теперь единственность решения задачи (2.16), (2.17) пли эквивалентной ей задачи (2.19), (2.20). Для этого необходимо показать, что при у = 0, например, задача (2.16), (2.17) имеет только тривиальное решение. В самом деле, уравнения системы (2.16)

и проинтегрируем по всей области определения решения D. Полу­ ченные интегралы преобразуем с помощью теоремы Гаусса—Остро­ градского и используем однородные граничные условия. Тогда получим

и используем тот факт, что компоненты u', v' и w' удовлетворяют уравнению неразрывности. В результате подстановки (2.28) в (2.26)

105