Investigation of the problem of representation of a linear functional in the form of a scoal product
- Authors: Ushakov A.L.1
-
Affiliations:
- South Ural State University
- Issue: Vol 18, No 3 (2022)
- Pages: 152-162
- Section: 1.2/1.2.2 Mathematical Modeling and Information Technology (Technical Sciences)
- Published: 08.10.2022
- URL: https://vestnikugrasu.org/byusu/article/view/111788
- DOI: https://doi.org/10.18822/byusu202203152-162
- ID: 111788
Cite item
Full Text
Abstract
Subject of research: the problem of representing a linear functional in the form of a scalar product in the Hilbert space. For example, this may be the Dirichlet problem for a second-order elliptic equation in variational form.
Purpose of research: to present a general scheme of the method of iterative extensions for solving elliptic boundary value problems with Dirichlet conditions.
Methods and objects of research: the object of study in such a problem can be the deformation of the membrane. A continuation of the considered problem in the Hilbert space is given. The extended problem is considered in a subspace of the Hilbert space. The extended problem is studied by the method of iterative extensions in the Euclidean space. An algorithm for implementing the method of iterative extensions is given.
Main results of research: as a result, it turns out that the general scheme of the method of iterative extensions as applied to the solution of a boundary value problem with Dirichlet conditions for an elliptic equation does not depend on the order of this equation. An example is given to illustrate the conclusion about the efficiency of the method of iterative extensions.
Full Text
Введение
Рассмотрим задачу представления линейного функционала в форме скалярного произведения в гильбертовом пространстве. Примером такой задачи является краевая задача с условиями Дирихле для эллиптического уравнения в плоской ограниченной области. Основные трудности при решении этой задачи возникают от сложности геометрии области, высоты порядка уравнения и наличия краевых условий Дирихле [1–5]. Метод решения такой задачи должен быть асимптотически оптимальны, являться вполне универсальным и иметь достаточно простую реализацию. Будем решать рассматриваемую задачу методом итерационных расширений, являющимся обобщением метода фиктивных компонент [4–7]. Решение исходной задачи рассматривается как обобщение решений методом итерационных расширений краевых задач с условиями Дирихле для эллиптических уравнений. Предполагается, что после сведения решения задачи для эллиптического уравнения к решению соответствующих задач в прямоугольной области можно будет применять оптимальные маршевые методы [8–10].
Результаты и обсуждение
1. Задача представления линейного функционала в форме скалярного произведения в гильбертовом пространстве
Пусть задана первая ограниченная область и выберем вторую ограниченную область
Области не пересекаются, а объединение их замыканий является замыканием прямоугольной области
Предполагаем, что границы этих областей имеют пересечение
На первой области рассмотрим задачу представления линейного функционала в форме скалярного произведения в гильбертовом пространстве, на второй области введем фиктивную задачу представления нулевого функционала в форме скалярного произведения в гильбертовом пространстве
(1)
Предполагается, что рассматриваемые скалярные произведения задаются с помощью симметричных операторов и удовлетворяют соответствующим свойствам скалярного произведения, являются скалярными произведениями, рассматриваемых гильбертовых пространств и обладают свойствами скалярных произведений, задаваемых посредством скалярных произведений функций суммируемых в квадрате
Примерами рассматриваемых пространств являются Соболевские пространства функций. Примерами операторов будут симметричные операторы, возникающие в полигармонических краевых задачах. Примером правой части у задач будет скалярное произведение функций, которые суммируемы в квадрате
У каждой из задач существует и единственное решение при этом, у фиктивной задачи это решение нулевое [1].
- Продолженная задача представления линейного функционала в форме скалярного произведения в гильбертовом пространстве
На прямоугольной области будем рассматривать пространство Гильберта из функций На этом пространстве зададим скалярное произведение как сумму скалярных произведений
Определим следующие подпространства
Отметим, что во введенном скалярном произведении
Обычно полагать, что для функций выполняются предположения о продолжении функций в следующей форме
Заведем оператор проектирования
Будем совместно рассматривать решаемую задачу и фиктивную задачу как продолженную задачу представления линейного функционала в форме скалярного произведения
(2)
Предложение 1. Имеют место равенства
Утверждение 1. Задача (2) имеет единственное решение , которое на . будет решением задачи (1) при , а на . будет нулевым решением задачи (1) при .
Для решений исходной и продолженной задач можно использовать одно обозначение как для функции и продолжения функции
3. Продолженная задача представления линейного функционала в форме скалярного произведения в конечномерном подпространстве
Рассматриваем введенную прямоугольную область и части ее границы в прямоугольных координатах
Полагаем, что введена прямоугольная сетка с узлами, постоянными и положительными шагами в направлении соответствующих осей координат
Введем сеточные функции, определенные в узлах сетки
При восполнениях сеточных функций используем в качестве базисных функций функции в узлах сетки, имеющие локальные носители
Полагаем, что эти базисные функции вне прямоугольной области будут равны нулю
Рассматриваем линейные комбинации из базисных функций, которые являются конечномерным подпространством
Важным примером аппроксимации пространств Соболева является кусочно-полиномиальная аппроксимация в виде конечных элементов при решении эллиптических краевых задач [1].
Приведем продолженную задачу в конечномерном пространстве.
Введем следующие подпространства
Отметим, что
Полагаем, что для функций и на конечномерном подпространстве выполняются предположения о продолжении функций в следующей форме
Будем использовать такой оператор проектирования под действием, которого обнуляются коэффициенты при базисных функциях, если носители этих функций не лежат полностью на первой области
Приведем продолженную задачу в матричном виде, определив матрицу и правую часть системы следующим образом
Так аппроксимируя продолженную задачу с помощью конечномерного подпространства, получим систему уравнений
(4)
Занумеруем первыми базисные функции с носителями в первой области. Вторыми занумеруем базисные функции с носителями, пересекающими границы первой области и второй области сразу. Третьими занумеруем базисные функции с носителями во второй области. При этой нумерации векторы имеют такой вид.
Отметим, что решается продолженная задача в следующем матричном виде
Так решается исходная задача в матричном виде и фиктивная задача в матричном виде.
Введем подпространства векторов
Определим матрицы
Для конечномерных подпространств можно записать прежние предположения, касающиеся продолжения функций в матричной форме
4. Метод итерационных расширений
Введем матрицу
Считаем, что выполняются предположения, записываемые в матричном виде
Рассмотрим для решения задачи (4) метод итерационных расширений как обобщение известного метода фиктивных компонент, используя введение дополнительного параметра. Метод фиктивных компонент при единичном значении этого параметра получается из метода итерационных расширений, если не учитывать выбор итерационных параметров
(5)
где при вычислениях итерационных параметров надо находить невязки, поправки, эквивалентные невязки
Введем норму
Лемма 1. Для метода итерационных расширений (5) имеет место оценка
Доказательство. Обозначим ошибки для итерационного процесса (5)
При начальном применении итерационного процесса выполняются следующие равенства
Можно заметить, что выполняется неравенство
Используя прежнее неравенство, получим следующие неравенства
При сокращении получим такие неравенства
Теорема 1. Для метода итерационных расширений из (5) имеют место такие оценки для сходимости приближенных решений к точному решению задачи из (4)
где относительные ошибки оцениваются сверху в более сильной норме, чем энергетическая норма бесконечно убывающей геометрической прогрессией.
Доказательство. В итерационном процессе получаем соотношения для ошибок, невязок
Будем минимизировать невязки
Определяем итерационные параметры при условии минимизации невязок
Получим равенства
Введем обозначения
Замечаем положительность итерационных параметров
Для выбранных итерационных параметров получается, что
Выписываем отношение скалярных произведений невязок для соседних итераций
Вводим обозначения
учитывая, что
Так устанавливаем неравенства
Учитывая, что
получаем неравенство для оценки сходимости у метода итерационных расширений
- Алгоритм метода итерационных расширений
Для выбора итерационных параметров используем метод минимальных невязок.
I. Берем начальное приближение, итерационный параметр
II.Вычисляем невязку
III. Вычисляем квадрат нормы абсолютной ошибки
IV. Находим поправку
V.Вычисляем эквивалентную невязку
VI. Вычисляем итерационный параметр
VII. Вычисляем очередное приближение
VIII. Устанавливаем критерий завершения работы по заданной оценке ошибки
6. Пример использования метода итерационных расширений
Рассмотрим следующие области
Пусть у областей следующие границы.
Рассмотрим задачу
Выберем правую часть, коэффициент в уравнении
Известно решение задачи
При дискретизации задаем шаги сетки
В прямоугольной области задаем сетку с узлами
Применяем метод итерационных расширений с параметром при нулевом начальном приближении. Процесс завершает работу всегда на восьмой итерации при заданной оценке ошибки На восьмой итерации при имеют место неравенства
Заключение и выводы
Разработанный асимптотически оптимальный метод итерационных расширений можно использовать при численном решении краевых задач Дирихле для эллиптических уравнений.
About the authors
Andrey L. Ushakov
South Ural State University
Author for correspondence.
Email: ushakoval@susu.ru
Candidate of Physical and Mathematical Sciences, Docent, Senior Research Officer at the Department
of Mathematical and Computer Modelling
References
- Aubin, J.-P. Approximation of Elliptic Boundary-Value Problems / J.-P. Aubin. – New York : Wiley-Interscience, 1972. – 360 p.
- Sorokin, S. B. An Economical Algorithm for Numerical Solution of the Problem of Identifying the Right-Hand Side of the Poisson Equation / S. B. Sorokin // Journal of Applied and Industrial Mathematics. – 2018. – Vol. 12, № 2. – P. 362–368.
- Sorokin, S. B. An Efficient Direct Method for the Numerical Solution to the Cauchy Problem for the Laplace Equation / S. B. Sorokin // Numerical Analysis and Applications. – 2019. – Vol. 12, №12. – P. 87–103.
- Ushakov, A. L. Investigation of a Mixed Boundary Value Problem for the Poisson Equation / A. L. Ushakov // 2020 International. Russian Automation Conference (RusAutoCon), Sochi, Russia. – 2020. – P. 273–278.
- Ushakov, A. L. Research of the boundary value problem for the Sophie GermainEquationinin in a cyberphysical system / А. L. Ushakov // Studies in Systems, Decision and Control. Springer. – 2021. – Vol. 338. – P. 51–63.
- Мацокин, А. М. Метод фиктивного пространства и явные операторы продолжения / А. М. Мацокин, С. В. Непомнящих. – Текст: непосредственный // Журнал вычислительной математики и математической физики. – 1993. – Т. 33, №1. – С. 52–68.
- Marchuk, G. I. Fictitious Domain and Domain Decomposion Methods / G. I. Marchuk, Yu. A. Kuznetsov, A. M. Matsokin // Russian Journal of Numerical. Analysis and Mathematical Modelling. – 1986. – Vol. 1, Iss. 1. – P. 3–35.
- Bank, R. E. Marching Algorithms for Elliptic Boundary Value Problems / R. E. Bank, D. J. Rose // SIAM Journal on Numerical Analysis. – 1977. – Vol. 14, № 5. – P. 792–829.
- Manteuffel, T. An Incomlete Factorization Technigue for Positive Definite Linear Systems / T. Manteuffel // Mathematics of Computation. – 1980. – Vol. 38, № 1. – P. 114–123.
- Swarztrauber, P. N. The Method of Cyclic Reduction, Fourier Analysis and FACR Algorithms for the Discrete Solution of Poisson’s Equations on a Rectangle / P. N. Swarztrauber // SIAM Review. – 1977. – Vol. 19, № 3. – P. 490–501.