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

задачи равно нулю, то для оптимизации следует брать наименьшее ненулевое.

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

т*

Что касается уравнения (7.21), то его решение было обсуждено выше. Некоторое отличие будет иметь место в алгоритме решения задачи (7.22). Оно состоит в том, что решение этих уравнений также следует теперь искать в виде рядов Фурье по собственным функциям задачи

Представив ит и ѵт в виде рядов

т* т*

мы приходим к системе уравнений для коэффициентов Фурье и",

н", которая решается на основе метода минимальных невязок в ком­ бинации с методом расщепления по схеме, близкой к рассмотренной выше.

Переходим теперь к теоретическому анализу сходимости итера­ ционного метода (7.1), (7.2 ) и (7.21), (7.22). Для того чтобы такой анализ провести до конца, мы вместо указанных задач рассмотрим другие, близкие к ним в дифференциальной форме, решения которых будем предполагать периодическими. Итак, сначала рассмотрим ите­ рационный процесс для задачи, соответствующей (7.1), (7.2), но записанной в форме, где р и w не исключены. Будем иметь

т

Будем считать, что любой из компонентов решения имеет вид

где ф — коэффициент Фурье, а а, ß и у — параметры волны

Подставим (7.26) в (7.25) и последовательно приведем систему <7.25) к соотношению для каждой из гармоник возмущения (7.26)

где Т — оператор шага, величина которого будет оценена на основе характерных масштабов возмущений.

Подставляя (7.26) в (7.25), получим

— (Ці7-2 + viY2) Pn+1 = Ггг" " ! * ,

где

г2= ct2 -j- ß2.

Будем теперь последовательно исключать все компоненты реше­ ния, кроме рп+1 и рп. Поскольку

YP

125

Для сходимости итерационного процесса необходимо выполне­

для всех гармоник, присутствующих в решении.

Если рассматривается более общая система уравнений, учиты­ вающая турбулентное трение в уравнениях Эйлера, то мы имеем более полную задачу:

Тогда аналогично предыдущему приходим к рекуррентному соот­ ношению (7.27), где оператор шага Т имеет вид

T = jL ll_________ (рг2+ѵѴа)________

р у 2 ( Р і '-2 + Ѵ1у 2 ) [ ( р г 2 - і- ѵ у 2 ) 2 - |- /2 ] • ( / . О О )

126

Дадим интерпретацию полученным критериям сходимости ите­ рационного процесса. С этой целью рассмотрим более общий слу­ чай (7.37). Прежде всего заметим, что для периодической задачи в разностной форме мы снова приходим к формуле (7.36), где пара­ метры г2 и тг связаны с h, Н, Ах = Ау и L, где L — максимальная длина волны в горизонтальной плоскости, формулами

Необходимо теперь с помощью формул (7.38) подсчитать границы изменения параметров а, ß и у и проверить выполнение условия (7.37). Конечно, условность в выборе L остается, так же как и пред­ положение о периодичности задачи, однако итерационный механизм,

•связанный с величиной Т в зависимости от cs, ß, у, и выводы о схо­ димости, а также достаточные условия применимости метода стано­ вятся более понятными и конструктивными.

4.8.

ДИ Н А М И КИ О КЕА НА С УЧЕТОМ Н ЕЛ И Н ЕЙ Н О ГО ТУ РБУ Л ЕН ТН О ГО ОБМ ЕНА

Следующим уточнением теории климата в океане является вве­ дение нелинейной турбулентности. Хорошо известно, что турбулент­

ный обмен в условиях устойчиво стратифицированной жидкости ~

0 может быть с достаточной точностью описан постоянным коэффициентом вертикального турбулентного обмена ѵ. Однако такое описание оказывается удовлетворительным только в этом слу-

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

где е фиксированная безразмерная константа феноменологического характера, которую можно выбирать порядка ІО-2 или 10_3. Будем

предполагать, что если условие ^ «S 0 реализуется на некоторой

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

127

подвергается весь слой 0 z ^ Н . Модель, близкая к такой модели, была рассмотрена Брайном.

Основные уравнения динамики при такой модели внешне не претерпевают каких-либо существенных изменений. В самом деле, возьмем в качестве исходной задачу (4.4) и определим метод после­

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

где р9, fg и I — сеточные функции к я I. Решение этой задачи про­ водится точно таким же методом, что и решение задачи (7.10). Важ­ но отметить, что принятое предположение о механизме турбулент­ ного обмена не нарушило возможности применения к расчетам метода разделения переменных по индексам к, I и т.

Аналогичным образом может быть рассмотрена более полная модель динамики (4.7).

4.9.

Н Е Л И Н Е Й Н Ы Е ПОСТАНОВКИ ЗА ДАЧ

Рассмотрим задачу динамики океана в более полной формули­ ровке без линеаризации, в квазистатическом приближении. Будем иметь систему уравнений:

цДи !- Ь, = -=” §£•+ Я И ,

128