АНАЛИЗ СМЕШАННОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ПУАССОНА Текст научной статьи по специальности « Математика»
Аннотация научной статьи по математике, автор научной работы — Ушаков А.Л.
Смешанная краевая задача для уравнения Пуассона рассматривается в ограниченной плоской области. Проводится продолжение этой задачи в вариационном виде через границу с условием Дирихле до прямоугольной области. Для решения продолженной задачи формулируется модифицированный метод фиктивных компонент в вариационном виде. Продолженная задача в вариационном виде рассматривается на конечномерном пространстве. Для решения предыдущей задачи формулируется модифицированный метод фиктивных компонент на конечномерном пространстве. Для решения продолженной задачи в матричном виде рассматривается известный метод фиктивных компонент . Показывается, что в методе фиктивных компонент абсолютная ошибка в энергетической норме сходится со скоростью геометрической прогрессии. В качестве обобщения метода фиктивных компонент предлагается новый вариант метода итерационных расширений . Продолженная задача в матричном виде решается методом итерационных расширений . Показывается, что в предложенном варианте метода итерационных расширений относительная ошибка сходится в норме более сильной, чем энергетическая норма задачи со скоростью геометрической прогрессии. Итерационные параметры в указанном методе выбираются с помощью метода минимальных невязок. Указываются условия достаточные для сходимости применяемого итерационного процесса. Выписан алгоритм, реализующий предложенный вариант метода итерационных расширений . В данном алгоритме производится автоматический выбор итерационных параметров и указывается критерий остановки при достижении оценки требуемой точности. Приводится пример применения метода итерационных расширений для решения частной задачи. В расчетах ставится условие достижения оценки относительной ошибки в норме более сильной, чем энергетическая норма задачи. Но приводятся относительные ошибки полученного численного решения примера исходной задачи и другими способами. Например, вычисляется поточено относительная ошибка в узлах сетки. Для достижения относительной ошибки не более нескольких процентов требуются всего несколько итераций. Вычислительные эксперименты подтверждают асимптотическую оптимальность метода, полученную в теории.
Похожие темы научных работ по математике , автор научной работы — Ушаков А.Л.
ANALYSIS OF THE MIXED BOUNDARY VALUE PROBLEM FOR THE POISSON’S EQUATION
The mixed boundary value problem for the Poisson’s equation is examined in a bounded flat domain. The problem is continued in a variational form through the boundary with the Dirichlet condition to a rectangular domain. To solve the continued problem, a modified method of fictitious components in a variational form is formulated. The continued problem in a variational form is considered on a finite-dimensional space. To solve the previous problem, a modified method of fictitious components on a finite-dimensional space is formulated. To solve the continued problem in matrix form, the known method of fictitious components is considered. It is shown that in the method of fictitious components the absolute error in the energy norm converges with the speed of a geometric progression. To generalize the method of fictitious components , a new version of the method of iterative extensions is proposed. The continued problem in matrix form is solved using the method of iterative extensions . It is shown that in the proposed version of the method of iterative extensions , the relative error converges in a norm that is stronger than the energy norm of the problem with a geometric progression rate. The iterative parameters in the specified method are selected using the minimum residual method. The conditions which are sufficient for the convergence of the applied iterative process are indicated. An algorithm which implements the proposed version of the method of iterative extensions is written. In this algorithm, an automated selection of iterative parameters is conducted, and the stopping criterion is established when achieving an estimate of the required accuracy. An example of the application of the method of iterative extensions for solving a particular problem is given. In the calculations, the condition for achieving an estimate of the relative error in the norm that is stronger than the energy norm of the problem is set. However, the relative errors of the obtained numerical solution of the example of the original problem are shown in other ways. For example, the relative error in grid nodes is calculated pointwise. To achieve a relative error of no more than a few percent, just a few iterations are required. Computational experiments confirm the asymptotic optimality of the method obtained in theory.
Текст научной работы на тему «АНАЛИЗ СМЕШАННОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ПУАССОНА»
АНАЛИЗ СМЕШАННОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ УРАВНЕНИЯ ПУАССОНА
Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: ushakoval@susu.ru
Смешанная краевая задача для уравнения Пуассона рассматривается в ограниченной плоской области. Проводится продолжение этой задачи в вариационном виде через границу с условием Дирихле до прямоугольной области. Для решения продолженной задачи формулируется модифицированный метод фиктивных компонент в вариационном виде. Продолженная задача в вариационном виде рассматривается на конечномерном пространстве. Для решения предыдущей задачи формулируется модифицированный метод фиктивных компонент на конечномерном пространстве. Для решения продолженной задачи в матричном виде рассматривается известный метод фиктивных компонент. Показывается, что в методе фиктивных компонент абсолютная ошибка в энергетической норме сходится со скоростью геометрической прогрессии. В качестве обобщения метода фиктивных компонент предлагается новый вариант метода итерационных расширений. Продолженная задача в матричном виде решается методом итерационных расширений. Показывается, что в предложенном варианте метода итерационных расширений относительная ошибка сходится в норме более сильной, чем энергетическая норма задачи со скоростью геометрической прогрессии. Итерационные параметры в указанном методе выбираются с помощью метода минимальных невязок. Указываются условия достаточные для сходимости применяемого итерационного процесса. Выписан алгоритм, реализующий предложенный вариант метода итерационных расширений. В данном алгоритме производится автоматический выбор итерационных параметров и указывается критерий остановки при достижении оценки требуемой точности. Приводится пример применения метода итерационных расширений для решения частной задачи. В расчетах ставится условие достижения оценки относительной ошибки в норме более сильной, чем энергетическая норма задачи. Но приводятся относительные ошибки полученного численного решения примера исходной задачи и другими способами. Например, вычисляется поточено относительная ошибка в узлах сетки. Для достижения относительной ошибки не более нескольких процентов требуются всего несколько итераций. Вычислительные эксперименты подтверждают асимптотическую оптимальность метода, полученную в теории.
Ключевые слова: уравнение Пуассона; метод фиктивных компонент; метод итерационных расширений.
Рассмотрим смешанную краевую задачу с однородными краевыми условиями Дирихле и Неймана для эллиптического уравнения Пуассона в ограниченной области. Основные проблемы при решении эллиптических задач обычно связаны со сложностью геометрии области, высотой порядка дифференциального уравнения и наличием краевого условия Дирихле 3. Будем исходить из желаемых положений, что предлагаемые численные методы должны являться устойчивыми к вычислительным ошибкам округления, быть асимптотически оптимальными по вычислительным затратам в арифметических операциях, являться достаточно универсальными и при этом иметь не сложную реализацию при автоматизации вычислительного процесса. Для возможного преодоления указанных трудностей и реализации указанных положений при решении рассматриваемой задачи используем в разных вариантах методы типа фиктивных компонент с применением модернизации и обобщения указанных методов, что приведет к их усложнению, но не кардинального характера 4. Отметим, что для решения задачи в прямоугольной области, к решению которой будет сводиться решение рассматриваемой задачи, можно использовать, например, известные маршевые методы оптимальные по вычислительным затратам 8.
Пусть задана первая ограниченная плоская область и выбирается вторая ограниченная плоская область.
Требуется, что бы пересечение этих областей было пусто, а объединение их замыканий было замыканием прямоугольной области.
w1 n Wii =0, W1 и WII = п.
У каждой из этих трех областей граница состоит из замыкания объединения двух открытых непересекающихся частей.
эп=s, s = Г1 и г2, Г1 n г2 =0.
= sa, sa = Га,1 U Га,2, Га,1 П Га,2 = 0.
Полагаем, что пересечение границ первой и второй области является замыканием непустого пересечения первой части границы первой области со второй частью границы второй области.
ЭП1 П ЭЙП = S, S = Г1Д П Г112 Ф 0.
Все рассматриваемые части границ у всех областей являются объединением конечного числа непересекающихся открытых дуг достаточно гладких кривых. Рассматриваются области границы, которых не имеют самопересечений и самокасаний.
2. Краевые задачи в операторном и вариационном видах
В первой области рассматриваем смешанную краевую задачу для уравнения Пуассона. Во второй области вводим смешанную краевую задачу для однородного экранированного уравнения Пуассона. На первых частях границ областей задаем однородное условие Дирихле. На вторых частях границ областей рассматриваем однородное условие Неймана. Задача на первой области является решаемой задачей. Задачу на второй области рассматриваем в качестве нулевого фиктивного продолжения решаемой задачи. Приведем решаемую и фиктивную задачу в операторном виде:
DUa + kaUa = fa, k = 0, kn > 0, /„ = 0,
Рассмотрим решаемую и фиктивную задачи в вариационном виде как задачи представления линейных функционалов в виде скалярных произведений в функциональных пространствах.
иае На : Аа(йаУа) = Ра(уа) «уае На. (2)
Пространства решений для таких задач будут следующие пространства функций Соболева
Правые части этих задач являются линейными функционалами
В левых частях этих задач стоят билинейные формы
Полагаем, что билинейные формы задают в пространствах решений рассматриваемых задач нормировки эквивалентные нормировкам соответствующих пространств Соболева
$С1,С2 > 0: С1 ||Уа| 12(«а) £ Аа(Уа,Уа) £ с2 |Ы»Уае На.
Эти предположения обеспечивают существование и единственность решения у каждой из рассматриваемых задач, физических систем [1]. Заметим, что решение фиктивной задачи нулевое.
3. Продолженная задача в вариационном виде
Возможно совместное рассмотрение решаемой и фиктивной задач в вариационном виде. Такую задачу будем называть продолженной задачей
UeV: АДы, I1v) + ЛП (U,v) = F1(I1v) «v eV. (3)
Расширенное пространство решений для такой задачи будет следующее пространство функций Соболева.
В расширенном пространстве решений находится подпространство являющееся пространством решений продолженной задачи. Это пространство решений исходной задачи в первой области, проложенное нулем на остальную часть прямоугольной области
В формулировке продолженной задачи используем оператор проектирования из расширенного пространства решений продолженной задачи на пространство ее решений
I1: V ^ V1, V1 = imI1, I1 = I12.
Дополнительно введем подпространства в расширенном пространстве решений
Мы использовали билинейную форму являющуюся суммой билинейных форм
А(ы,v) = А1(ы,v) + Ajj(и,v) «и,v e V. Полагаем, что билинейная форма задает в расширенном пространстве решений нормировку эквивалентную нормировке соответствующего пространства Соболева.
2 2 3c c0 > ü: c v 1 „ £ A(v v) £ Со v 1 „ «v e V
^ » «1 ir Ilw2 (П) v ‘ ‘ 2\г Ilw2 (П) Предполагаем, что используемые пространства Соболева таковы, что в них возможно продолжение функций с сохранением нормы. Обычное в таких случаях это положение будем использовать в указываемом виде
$Д e (ü;1], b e [Д;Ч: AA^,v2) £ An(v2,v2) £ bA(y2,v2) «v2 e V2. В таком случае у продолженной задачи существует и единственное решение. Решение продолженной задачи — это решение исходной задачи в первой области, продолженное нулем на остальную часть прямоугольной области.
4. Модифицированный метод фиктивных компонент в вариационном виде
Рассмотрим итерационный процесс, на каждом шаге которого решаем расширенную задачу с билинейной формой из продолженной задачи, но без оператора проектирования. Решение этой задачи ищем в расширенном пространстве решений продолженной задачи, в пространстве решений расширенной задачи
uk eV: A(uk — uk -1, v) = -tk _1 (A1 (uk-1, I1v) + An (uk-1, v) — F1 (I1v)), k e N «v e V,
«uüe V, t = 1.1_1 = 2/(b + Д), ke N\<1>. (4)
Этот итерационный процесс имеет оценки сходимости
V, e = cqk-1, ke N, c =>/ЙП,ü £q = (fa _Д)/(Д + b) Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
+; = V’, 1, Ф«(г-1)+; = Ф 1 (х» ^ У] ^ » = 1,2. га 1 = 1,2. «
у=ЕЕ у», 1Ф» 1 (х; у)=Е V ф I.
Но нумеруем сначала базисные функции, носители которых полностью содержатся в первой области. Далее продолжаем нумеровать базисные функции, носители которых пересекают границу первой области и второй области одновременно. Завершаем нумерацию на базисных функциях носители, которых полностью содержатся во второй области. При такой нумерации возникающие векторы имеют следующий вид
V = у, , Уз’)’, й=(«1′ ,0′ о’)’, у=(у ,0′ о’)’.
Вычисляем элементы матрицы, компоненты вектора правой части, приведенной системы Ьу = ¿Г^ЧАДФ,, 11Ф1) + Ап(Ф,,Ф1)), у = /^^ВДФ,),», 1 = 1,2. N.
9. Метод фиктивных компонент в матричном виде
Аппроксимируя модифицированный метод фиктивных компонент в вариационном виде с помощью конечномерного подпространства при указанном ранее операторе проектирования, получим известный метод фиктивных компонент в матричном виде
йк е МN : А(йк — йк-1) = -тк-1(Вйк-1 — у), к е М,
«й0 еУх, То = 1. Тк-1 = 2/(Д + Д2), к е М\ <1>. (8)
На каждом шаге этого итерационного процесса получаем расширенную задачу в матричном виде с расширенной матрицей
(Ай,V) = А(й,V) «й,Vе I7. Вычисляем элементы этой матрицы
а,! = АТ1Л2-1А(Ф1., Ф1), 1,1 = 1,2. N. Возникающие матрицы имеют известную структуру
А11 А12 0 А11 А12 0
А21 А22 А23 ,В= 0 А02 А23
0 Аз2 Азз _ 0 Аз2 Азз
Используем подпространство векторов
Дополнительно введем подпространства векторов, как и предыдущие соответствующие конечномерным подпространствам введенным ранее
W,v’)’e MN : Allvl + A12v2 = 0, A23v2 + A33v3 = o>, MN = V © V © V.
Отметим, что в методе фиктивных компонент решается продолженная задача в матричном
«An A12 0 » u1 fx
Bu = f, 0 A02 A23 0 = 0
0 A32 A33 _ 0 0
По решению продолженной задачи в матричном виде получаются решения соответственно исходной задачи в матричном виде и нулевое решение фиктивной задачи в матричном виде
A02 A23 A32 A33
Введем нормы, порождаемые единичной матрицей, расширенной матрицей и ее квадратом
Лемма 1. В итерационном процессе метода фиктивных компонент (8) выполняется оценка
0: (А у1,у1) Не можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
фи еУ ((0; Ьх) х (0; Ъ2)):-Л^ = Л,Д,, фи * 0, Л,, = 0,25р2 ((2/ — 1)Ь-2 + (2; — 1)Ь-2), /, , е N.
Мы рассматривали случай, когда константа экранирования в фиктивной задаче равна нулю, полагая, что иначе рассматриваемая в лемме 2 оценка все равно не должна зависеть от параметров дискретизации.
Теорема 1. В итерационном процессе метода фиктивных компонент (8) выполняются оценки сходимости
, $ae (0; ■+»): ^Ajv2, AjV2)
Далее используем только последнее неравенство (CVj,CV2) = (AjV2, AjV2) + 2gAjV2, Ajjv^ + g2 (Anv2, Anv^ a, 4 = 1, 4—1 =(rk—1,nk—1)/
Для вычисления итерационных параметров необходимо вычислять невязки, поправки и эквивалентные невязки соответственно
гк-1 = Бйк-1 — 7, *>к-1 = С
1гк-1, Г!к-1 = В^к-1, к е N. Лемма 3. в методе итерационных расширений (9) выполняется оценка
Приводим скалярные произведения невязок при выбранных итерационных параметрах
тгк -к\ -к -1 —к -1х \ 11
Выписываем отношение квадратов норм невязок на соседних итерациях
2 (гк, гк) _ (Л11с-1гк-1,гк-и2
Гк -1, Гк-1) (Л11С»Ук-1, Л11С «Ук-1) <гк -1, Гк-1
(лп як-1,л11 як-1)(с^к-1,сяк-1)- (лп як-1,сяк-1 )2 (Ъ,Ъ)(а + уь,а+уь)- (Ъ,а+уь)2
(л11 Як -1, Л11 як-1> (сяк -1, сяк-1) (ь, ь)(а + уЬ, а+уь) ‘
(а, а) = а, (ь,ь ) = ь, (о,ь ) = г,
2 аь — г2 ^ 2 Г -а ^ а а2 2
ь(а + у2ь + 2уг) И^ у ) у^ь у2
2^п/ ^ч-2у(г + а/у)(г + уь) , а + у2ь ¡—г а г—т
4к ^ 0, (4к (г)) =-1: \ ч2 , — уь Надоели баннеры? Вы всегда можете отключить рекламу.