2
Введение
В современной науке и технике требуется решать задачи математической физики в областях сложной геометрической формы. При этом требования к точности получаемого решения зачастую превосходят возможности современной вычислительной техники. Процесс численного решения многих проблем, содержащих дифференциальные уравнения, включает, по крайней мере, две трудоемкие задачи: построения вычислительной дискретной сетки в области решения задачи и построение и реализацию алгоритма, вычисляющего приближенное решение исходной задачи. Очень часто проблема построения вычислительной сетки сравнима, а иногда превосходит по затратам компьютерного и человеческого труда задачу получения и непосредственного решения дискретных уравнений, аппроксимирующих исходную задачу.
В то же время, современные алгоритмы вычислительной математики позволяют "быстро" решать только ограниченный класс проблем. В частности, к ним относятся задачи в областях простой формы, допускающие разложение решения в ряды Фурье. Другую возможность эффективного решения проблем математической физики в областях сложной формы дают многосеточные алгоритмы, но они требуют построения иерархической последовательности сеток в области решения задачи, что в свою очередь существенно ограничивает и усложняет их применимость.
Возможной альтернативой задаче построения сложных сеток является метод фиктивных областей. Общая схема применения метода фиктивных областей для решения краевых задач математической физики состоит в следующем. Пусть в некоторой области С Н’г ищется решение и(х), которое удовлетворяет уравнению с частными производными
Ьи = /(я), X = (яь £2.Хп) € Сь и граничным условиям
1и = д(х), х € дО\.
Основная идея метода фиктивных областей состоит в том, чтобы решать задачу не в исходной сложной области (?ь а в некоторой другой, более простой области С такой, что С С. Таким образом, мы будем иметь дело только с регулярными областями, что позволяет создавать программное обеспечение сразу для достаточно широкого класса задач с произвольными расчетными областями. В качестве (7 можно выбрать, например, п-мерный параллелепипед.
з
Задачу в расширенной области (7 для приближенного решения иш(х) можно записать в виде
їй\Jllu; X Є (?,
км* = 9и{х), х £ дЄ.
Здесь отмечено, что краевая задача в Є содержит большой параметр и, который определяет величину разрыва коэффициентов оператора Ьи на границе исходной области (?і и фиктивной С\С[. Необходимо так по-строить вспомогательную задачу для определения ии(х), чтобы имела место сходимость приближенного решения иш(х) к точному в области С1 при и —► сю. При этом выбор вспомогательной задачи не единственен.
Сдерживающим моментом для практического применения метода фиктивных областей является то, что задача для определения и„(х) - это задача с сильноразрывными коэффициентами, задача с большим параметром ш. Поэтому затрудняются вопросы построения соответствующих разностных схем и численного решения сеточных задач. В практическом плане все преимущества метода фиктивных областей могут быть сведены на нет, если скорость сходимости применяемых итерационных процессов для определения приближенного решения и^(х) будет слишком медленной.
Приведем краткий обзор литературы по вопросам, затронутым в диссертации. Подробный обзор и ссылки на работы, посвященные обоснованию и практическому применению метода фиктивных областей можно найти в книге Г1.Н. Вабищевича [18]. В работе Н С.Бахвалова [3] наиболее полно отражены современные итерационные методы для решения задач с сильноразрывыми коэффициентами.
Впервые метод фиктивных областей, как метод приближенного решения краевых задач в сложных областях с помощью ЭВМ формулируется в работах В.К.Саульева [31 ]. Им рассмотрена задача Дирихле для эллиптического уравнения второго порядка, предложен вариант с продолжением по старшим коэффициентам и получены оценки близости приближенного и точного решений непрерывных задач. В.И.Лебедевым [26] предложен вариант метода фиктивных областей с продолжением по младшим коэффициентам. В работе В.Я.Ривкинда рассмотрена третья краевая задача. Развитие общих направлений метода фиктивных областей для различных задач математической физики отражено в работах А.Н.Бугрова, А.Н.Коновалова, С.А.Войцеховского и других авторов.
Одним из крупных направлений практического использования метода фиктивных областей являются задачи гидродинамики. Теоретическому
4
исследованию и обоснованию метода фиктивных областей для уравнения Навье-Стокса посвящены работы А.Н.Бугрова [16], Ш.Смагулова [33], М.К.Оруханова [32] и др. Практическое применение и численный эксперимент отражены в работах Л.А.Руховца, П.Н.Вабищевича, 1Магйкатеп [29] и др. В работе А.Н. Бугрова и Ш.Смагулова [17] приведено обоснование метода фиктивных областей для уравнения Навье-Стокса с продолжением по старшим коэффициентам в односвязном случае и по младшим коэффициентам в неодносвязном случае.
Вопросы построения итерационных алгоритмов для решения задач метода фиктивных областей с продолжением по старшим коэффициентам наиболее полно рассмотрены в работе Н С.Бахвалова [3]. Впервые специальный метод для решения задачи метода фиктивных областей с продолжением по старшим коэффициентам для уравнения Лапласа предложен в работе Г.М.Кобелькова [22]. Задаче Дирихле для эллиптических операторов посвящены работы Г.М.Кобелькова [21], Н.С.Бахвалова [8], задача Неймана рассмотрена в [4]. Смешанной задаче для квазилинейного эллиптического уравнения посвящена работа К.Ю.Богачева [11]. Вопросы построения и обоснования итерационных алгоритмов для решения методом фиктивных областей задач упругости рассмотрены в [7, 2, 5, 6]. Нестационарная задача Стокса рассмотрена в [3].
Настоящая работа посвящена построению, исследованию и применению метода фиктивных областей с продолжением по старшим и младшим коэффициентам одновременно для решения задач математической физики в неодносвязных областях. Исследование проводилось в непрерывном случае. Особое внимание уделялось применению метода фиктивных областей к уравнению Навье-Стокса и эффективному решению получающихся задач с сильноразрывными коэффициентами.
Основные результаты работы заключаются в следующем:
1. В работе предложен новый вариант метода фиктивных областей для решения задачи Дирихле для различных типов уравнений математической физики в многосвязном случае. Достоинствами предложенного продолжения уравнений с частными производными являются возможность решать первую краевую задачу в неодносвязных областях и высокая скорость сходимости решения продолженной задачи и^х) к решению исходной и(х).
2. Обосновано применение построенного метода для решения уравнений эллиптического типа, линеаризованной задачи Стокса и уравнения Навье-Стокса, ког да граница области может быть неодносвязна. В случае уравнений эллиптического типа и задачи Стокса показана скорость
5
сходимости метода фиктивных областей не ниже ш 1.
3. Построены итерационные алгоритмы решения задач метода фиктивных областей с сильноразрывными коэффициентами. Показано, что скорость сходимости предложенных методов не зависит от разброса коэффициентов.
Диссертация состоит из введения, пяти глав, заключения и списка цитируемой литературы. Общий объем диссертации - 89 печатных страниц.
Результаты диссертационной работы докладывались и обсуждались на семинаре кафедры Вычислительной математики Механико-математического факультета МГУ под руководством академика РАН Н С. Бахвалова, научных семинарах Института Вычислительной Математики под руководством академика РАН Н.С. Бахвалова, В.И. Лебедева, Института Математического Моделирования под руководством A.B. Гулина, П.Н. Вабишевича, Вычислительного центра РАН под руководством Б.В.Пальцева. Основные результаты были доложены на Лобачевских чтениях КГУ (Казань, 2001), Ломоносовских чтениях МГУ (Москва, 2002), V Международном Конгрессе по Математическому Моделированию (Дубна, 2002).
По материалам работы имеется четыре публикации ([15, 14, 13, 12]).
Автор выражает искреннюю признательность научному руководителю академику РАН Бахвалову Николаю Серг еевичу за внимательное отношение к работе и ценные советы.
Глава 1
Обоснование метода фиктивных областей в неодносвязном случае
В этой главе предлагается вариант метода фиктивных областей с продолжением по старшим и младшим коэффициентам. Показано, что решение вспомогательной задачи в расширенной области приближает решение исходной при и —> ос для первой краевой задачи в случае эллиптических уравнений, стационарной и нестационарной задачи Стокса и нелинейного уравнения Навье-Стокса.
1.1 Эллиптические уравнения
Рассмотрим задачу Дирихле для эллиптического уравнения
-1г (ф)!^) = /(*). и 6 Я01(СО, /(*) € Ь2(С,), (1.1)
в двусвязной области й\, где а(х) € ЬХ(С\), 0 < ао < а(я) < а\. Пусть (7 - липшицева область такая, что С?1 С С и задача Дирихле для уравнения Лапласа в этой области допускает эффективное решение (например, С - параллелепипед).
Введем обозначения. Пусть (7 = (?2 и дИч и &1 и дО\ и (7з, где Счб\С\ = 0, С1 ПОз = 0, односвязная область (72 со всеми граничными точками лежит строго внутри 012 = <7г и <9(7г и бь а 6*12 = С\2 и дО^ строго внутри О = С$ и дО\2 и (712.
6
- Київ+380960830922