CC BY
62
УДК 517.3
М. Ю. Медведик
СУБИЕРАРХИЧЕСКИЙ МЕТОД РЕШЕНИЯ ИНТЕГРАЛЬНОГО УРАВНЕНИЯ ЛИППМАНА - ШВИНГЕРА
Аннотация. Рассмотрено решение интегрального уравнения Липпмана -Швингера. Представлен численный метод Галеркина. Получены численные результаты для расчета акустического поля внутри тела с использованием су-биерархического алгоритма.
Ключевые слова: субиерархический алгоритм, интегральное уравнение, численный метод, краевая задача.
Abstract. This paper is considered solution integral equation of Lippmann -Schwinger. The integral equation is solved by Galerkin method. Numerical results of solving of are obtained by using subhierarchical algorithm by body free form. Keywords: subhierarchical algorithm integral equations, numerical method, boundary value problem.
Введение
Определение рассеянного поля в различных материалах и средах является актуальной задачей акустики и электродинамики. Так как точные решения задач дифракции могут быть получены лишь для ограниченного числа тел правильной геометрии, то большое значение для практических приложений представляет развитие различных приближенных и численных подходов, справедливых для тел произвольной формы. Таким образом, возникает необходимость разработки новых методов решения подобных задач. Одним из перспективных методов является метод объемных сингулярных интегральных уравнений [1]. С помощью него краевая задача сводится к решению объемного сингулярного интегрального уравнения. Решение получающегося интегрального уравнения в общем случае возможно лишь численными методами, но благодаря снижению размерности задачи за счет сведения к интегралу по поверхности происходит значительное упрощение численных расчетов. Решение таких задач с приемлемой для практики точностью требует очень большого объема вычислений. Представленный в статье метод позволяет решать подобные задачи на телах сложной геометрической формы, опираясь на результаты, полученные при решении задачи на теле базовой формы [2-9].
Постановка задачи
Рассмотрим задачу дифракции акустической волны на теле Q, распо-
3
ложенном в свободном пространстве R (рис. 1).
Пусть дано неоднородное уравнение Гельмгольца
Ап + k2 (x^ju = f (x), (1)
где f (x) - известная функция с компактным носителем.
Будем предполагать, что на границе раздела двух сред выполняются условия сопряжения
[и ]эе = °.
и условия излучения Зоммерфельда
Эu
Эп
= °
ЭQ
Эи (1
— = іки + o\ — і, г := |х| Эг ^ г ,
(2)
(3)
E
III
(4)
Представим данное уравнение в виде (1)
Дu + к(и = (( -k2 (х))и + /(х),
где kQ - волновое число в свободном пространстве; k(х) - функция, характеризующая волновое число внутри тела Q .
Обозначим через F(х) = (( - k2 (x))и + /(x) правую часть уравнения (4). Тогда, используя вторую формулу Грина, получаем |( + kQ )и(х)О(х,у)-( + kQ )(х,у)и (х)с1:
ях =
V
= /( ^ - (х., )^ и (х)
где 51 - сфера, а V - ее объем (рис. 1). Учитывая, что F(х) = ( + к(). а ( + к( )(х,у) = §(х), приходим к следующей формуле:
и (х )= Г F (х )G (х, у )Лх - Г —(х, у)-( , У) и (х )ds.
* * Яп Яп
Интеграл по поверхности S в правой части уравнения стремится к нулю при r . Таким образом, задача свелась к следующему уравнению:
u (x) = J F (x )G (x, y )dx.
V
3
Устремим радиус S к бесконечности и перейдем от V к R . Учитывая, что F(x) = (o -k2 (x))u + f (x), получаем
u(x) = JG(x,y)( -k2 (y))u(y)dy + J f (y)G(x,)dy .
Q R3
Обозначим в правой части уравнения интеграл через
f0 (x )= Jf (y )G (x, y )dy.
R3
В результате приходим к уравнению, известному в литературе как интегральное уравнение Липпмана - Швингера:
u(x) = f0 (x) + JG(x,y)( -k2 (y))u(y)dy. (5)
Q
Это уравнение играет чрезвычайно важную роль не только в акустических задачах дифракции, но также в электродинамике, квантовой механике и во многих других областях физики.
Для однородного тела, т.е. такого, для которого k (x)- const, уравнение модифицируется и принимает вид
u (x ) = f0 (x) + k2 J G (x, y )u (y )dy. (6)
Q
Уравнение (5) является уравнением Фредгольма второго рода. К нему применима известная теория Фредгольма. Вопросы единственности и существования решения данного уравнения рассмотрены в [10].
Метод Галеркина
Рассмотрим «-мерное пространство Vn . Проведем аппроксимацию элементов ^ элементами ^n е Vn . Методом Галеркина находим ^n из системы уравнений
(Ly n,v) = (f,v). (7)
Эти уравнения определяют конечномерный оператор Ln : Vn ^ Vn, где Vn есть антидуальное пространство к Vn .
В качестве базисных функций выберем функции v\ . Будем считать, что Q - параллелепипед: Q = {x: a1 < x1 < a2, b1 < x2 < b2, C1 < x3 < C2} . Разобьем Q на элементарные параллелепипеды рис. 2:
Пкіт {х : х1.к ^ х1 ^ х1.к+1. х2,ї ^ х1 ^ х2,ї+1. х3.т ^ х3 ^ х3.т+1} .
хі,к = а + °2 а1 к. х2, = *1 + *2 Ь і.
с2 _ С х3. т = С1 +—-1 т.
п п п
где к.і.т = 0,...,п -1. Объем любого параллелепипеда Пкіт равен уоі .
Q
Базисные функции укіт . і = 1,2,3 . определяются следующим образом:
і I Г, х Є Пкіт ,
укіт
0, х ^П
кіт.
Построенное множество базисных функций удовлетворяет необходи-
3
мому условию аппроксимации [11] в Ь2 = 1,2 XІ XІУ2.
Описанный выше метод приводит к решению матричного уравнения. Произведем перебор всех пар элементарных параллелепипедов. Каждый элемент матрицы получается путем вычисления шестикратного интеграла
іі у = уоі -5г- у - | О(х, у)уі (х) • Уу (у^, имеющего слабую особенность в об-й
ласти интегрирования. Процедура избавления от особенности представлена в [4, 7]. Правая часть матричного уравнения задается формулой
Л\х-у\
/і = І /¥я . Здесь х = (хьх2,х3), у = (У1,У2,У3), а О(х,у) = ---------г
й Iх - у
вестная функция. Решение СЛАУ производится методом сопряженных градиентов.
- из-
Субиерархический алгоритм
Алгоритм расчета акустического поля внутри фигуры в форме параллелепипеда Q описан выше. Рассмотрим алгоритм расчета акустического поля для тела сложной геометрической формы. Будем предполагать, что решение задачи тела Q получено, и в нашем распоряжении находится матрица, составленная методом Галеркина. Для решения задачи дифракции акустической волны на теле сложной формы необходимо, чтобы тело целиком вмещалось в параллелепипед Q и состояло из элементов сетки [3-10]. Субиерархиче-ский метод позволяет составить матрицу для определения акустических полей внутри тела сложной конфигурации, используя матрицу, составленную для параллелепипеда Q. В построенной фигуре введем новую нумерацию элементарных параллелепипедов. Произведя полный перебор всех элементарных параллелепипедов, принадлежащих фигуре сложной конфигурации, получаем новую сетку. Эту сетку будем использовать для расчета поля на теле сложной формы. Решая СЛАУ для матрицы, составленной с использованием новой сетки, находим значения поля внутри фигуры сложной формы. Скорость построения новой матрицы будет напрямую зависеть от размера фигуры и размера сетки. Субиерархический подход позволяет избежать длительных расчетов, связанных с вычислением матричных элементов.
Численные результаты
На рис. 3-5 приведен расчет значений акустического поля внутри тела сложной геометрической конфигурации. Параллелепипед Q покрыли равномерной сеткой 11^11x11 элементов, используя субиерархический алгоритм, произвели расчет поля на представленной ниже фигуре. Первые четыре слоя по направлению оси 02 фигуры имеют вид, представленный на рис. 3, следующие три слоя фигуры представлены на рис. 4, последние четыре слоя изображены на рис. 5.
Рис. 3. Значение модуля акустического поля во втором сечение фигуры перпендикулярно оси 02
Преимущества субиерархического метода особенно хорошо заметны при расчете серии задач на телах различной конфигурации.
Рис. 4. Значение модуля акустического поля в пятом сечение фигуры перпендикулярно оси 02
Рис. З.Значение модуля акустического поля в восьмом сечение фигуры перпендикулярно оси OZ
Список литературы
1. Самохин, A. Б. Интегральные уравнения и итерационные методы в электромагнитном рассеянии I A. Б. Самохин. - М. : Радио и Связь, i99S.
2. Медведик, М. Ю. Субиерархический метод решения интегрального уравнения на плоских экранах произвольной формы I М. Ю. Медведик II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. -2009. - № 4. - С. 49-ЗЗ.
3. Медведик, М. Ю. Субиерархический параллельный вычислительный алгоритм для решения задач дифракции электромагнитных волн на плоских экранах I М. Ю. Медведик, Ю. Г. Смирнов II Радиотехника и электроника. - 200S. - Т. ЗЗ. -№ 4. - С. 441-446.
4. Медведик, М. Ю. Применение ГРИД-технологий для решения объемного сингулярного интегрального уравнения для задачи дифракции на диэлектрическом теле субиерархическим методом I М. Ю. Медведик, Ю. Г. Смирнов II Извес-
тия высших учебных заведений. Поволжский регион. Физико-математические науки. - 200S. - № 2. - С. 2-14.
З. Медведик, М. Ю. Субиерархический параллельный вычислительный алгоритм и сходимость метода Галеркина в задачах дифракции электромагнитного поля на плоском экране I М. Ю. Медведик, Ю. Г. Смирнов II Известия высших учебных заведений. Поволжский регион. - 2004. - № З. - С. 3-19. - (Естественные науки).
6. Медведик, М. Ю. Субиерархический подход для решения объемного сингулярного интегрального уравнения задачи дифракции на диэлектрическом теле в волноводе методом коллокации I М. Ю. Медведик, Д. А. Миронов, Ю. Г. Смирнов II Известия высших учебных заведений. Поволжский регион. Физикоматематические науки. - 20i0. - № 2. - С. З2-4З.
7. Медведик, М. Ю. Субиерархический метод для решения псевдодифференци-ального уравнения в задаче дифракции в слоях, связанных через отверстие I М. Ю. Медведик, И. А. Родионова II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2009. - № 3. - С. З9-71.
S. Антонов, А. В. Разработка Web-ориентированного вычислительного комплекса для решения трехмерных векторных задач дифракции электромагнитных волн на основе субиерархических параллельных алгоритмов и ГРИД-технологий I А. В. Антонов, М. Ю. Медведик, Ю. Г Смирнов. II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2007. - № 4. -С. З1-б0.
9. Медведик, М. Ю. Параллельный алгоритм расчета поверхностных токов в электромагнитной задаче дифракции на экране I М. Ю. Медведик, Ю. Г. Смирнов, С. И. Соболев II Вычислительные методы и программирование. - 200З. - Т. б. -С. 99-i0S.
10. Колтон, Д. Методы интегральных уравнений в теории рассеяния I Д. Колтон, Р. Кресс. - М. : Мир, i9S7.
11. Kress, R. Linear Integral Equations I R. Kress II Applied Mathematical Sciences. Springer-Verlag. - New-York, 19S9. - V. 82.
Медведик Михаил Юрьевич
кандидат физико-математических наук, доцент, кафедра математики и суперкомпьютерного моделирования, Пензенский государственный университет
Medvedik Mikhail Yuryevich Candidate of physical and mathematical sciences, associate professor, sub-department of mathematics and supercomputer modeling,
Penza State University
E-mail: [email protected]
УДК 517.3 Медведик, М. Ю.
Субиерархический метод решения интегрального уравнения Липпмана - Швингера / М. Ю. Медведик // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2010. -№ 4 (16). - С. 82-88.
