CC BY
38
УДК 539.3+517.44
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ПОТЕНЦИАЛЬНЫХ ПОЛЕЙ МЕТОДОМ ГРАНИЧНЫХ ОПЕРАТОРОВ ПРЕОБРАЗОВАНИЯ
© 2012 Ю.А. Парфенова1
В статье предлагается метод граничных операторов преобразования, который дополняет классические методы и позволяет найти решение краевых задач в замкнутом виде, что открывает новые возможности в моделировании потенциальных полей. Применение разработанного нами метода позволило решить ряд задач теории потенциала и применить для исследования пакет прикладных программ MatLab.
Ключевые слова: моделирование, потенциальное поле, граничные операторы преобразования.
Введение
При исследовании различных проблем, связанных с излучением стационарных полей, например гравитационного или магнитного, с целью поиска полезных ископаемых, возникают задачи определения краевого условия. В работе предлагается метод векторных операторов преобразования, который дополняет классические методы и позволяет найти решение краевых задач в замкнутом виде, что открывает новые возможности в исследовании моделей потенциальных полей. Применение разработанного нами метода позволило найти в замкнутой форме решения прямых и обратных задач математической физики со сферической симметрией. Актуальность метода граничных операторов преобразования состоит в том, что он позволяет упростить вычислительные схемы при применении методов итерации и регуляризации.
1. Обобщение формулы Пуассона в круге как основа для конструкции граничных операторов преобразования
Ряд важных задач теории потенциала состоит в определении неизвестных граничных данных по другим известным данным. Например, по условиям Неймана
Парфенова Юлия Алексеевна ([email protected]), кафедра математического анализа Пензенского государственного педагогического университета им. В.Г. Белинского, 440026, Российская Федерация, г. Пенза, ул. Лермонтова, 37.
восстановить условие Дирихле. Будем обозначать соответствующую задачу символами i ^ j, где i,j = 1,2,3. Так, 2 ^ 3 — это обозначение задачи определения граничного условия третьего рода по условиям второго рода.
В самой общей постановке определение граничного оператора преобразования можно ввести следующим образом:
пусть мы имеем две краевые задачи с различными граничными условиями для уравнения Лапласа в единичном круге B:
Дм = 0, Дм = 0,
^s = f (ф), rju\s = f (ф),
где i,j = 1, 2, 3, i = j, тогда определен оператор Lji : м ^ м, который по решению j-й краевой задачи находит решение i-й краевой задачи. Если функции м и М допускают непрерывные продолжения в замкнутый круг, то возникают граничные операторы преобразования j : Г [м] \s ^ rj [м] \s, позволяющие по известным значениям на границе одного интегродифференциального выражения Г [м]^ вычислять значения другого интегродифференциального выражения rj [м]^.
В 70-х годах XX века И.И. Баврин [1] получил обобщение классической формулы Пуассона для круга:
оо / 2п \
м (r, Ф) = — e-£ \ 2веr cos(cos s— sin (ф — ф) — 11 м (r0, ф) 3,ф0,е. (1.1)
2п J \J —о /
о \о /
Формула (1.1) позволяет восстановить гармоническую функцию м в круге по ее значениям на внутренней окружности радиуса —о. На основе формулы (1.1) образованы явные конструкции для граничных операторов преобразования j : ri [м] \s ^ rj [м]\s :
о / 2п \
rj [м]^ = 2"/e-£Lij ¡J 2e£ cos(ф-'ф) cos (s sin(ф — ф)) — 11 Г, [м (1, ф)] ¿фйе. (1.2)
оо
В ядре формулы (1.2) действуют операторы преобразования Lj, определенные автором в статье [2].
2. Граничные операторы преобразования Г12 и Г21
Граничные операторы преобразования позволяют получить в замкнутом виде решения ряда граничных задач теории потенциала. Представим результаты.
Теорема 1. Для задачи 1 ^ 2 граничный оператор преобразования Г12 имеет вид:
Г12 : I ^ /,
1
оо 2п
f = ¡n7e"
L
12
ee cos(0-V>)
(2 cos s sin (ф — ф) — 1) f^^ds. (2.1)
оо
f предполагается такой, что сходится числовой ряд
J2k (\ak \ + \bk\) <
(2.2)
к=1
2п
ak = / COS к^ ^ f ф
о
2п
bk = / sin k^ • f ф dф
о
— коэффициенты ее ряда Фурье.
Доказательство. Пусть и (r, ф) — решение задачи Дирихле для уравнения Лапласа в круге:
Дм = 0,
U\s = f (ф) •
Из формулы (1.1) следует интегральное представление:
сю 2п
L0 [U] (r, ф) = J e-e J L12 e£r сов(ф-ф) (2cos(re sin^ - ф)) - 1) f (ф)dфd£.
оо
Выполним формальный предельный переход r ^ 1 в приведенной формуле. Учитывая, что L0 [U]\r=1 = f (ф) , находим:
с 2п
f (ф) = 2П / е-£ j L12 [e£ сов(ф-ф) (2 cos (е sin (ф - ф)) - 1)] f (ф^е.
00
Законность предельного перехода обосновывается требованием (2.2) и элементарно проверенным тождеством:
2п
у e£r с°*(ф-ф) (2 cos (е sin (ф - ф)) - 1) f (ф)dф =
0
с k k i k k
k k k k
ao , akе rk bkekrk
= y + ^ —k!—cos кф +--k! sin кф. (2.3)
k=i ' ' Замечание. Рассматриваемая задача некорректна, поэтому нами предложены ре-гуляризирующие алгоритмы, основанные на формуле (1.2). Представлена их реализация в вычислительной среде MatLab.
Пример 1. По значениям потенциала и = | sin ф - 4 sin 3ф на границе единичного круга нужно восстановить поток и'п на границе. Точное решение этой задачи имеет вид:
3 з
иП (г,ф)\г=1 = 4sin ф - 4sin3^ Опишем алгоритм решения. По значениям наблюдаемой функции
f (ф) = (j sin ф - 1 sin 3ф^ (1 + se),
где S = 10-4, e — случайное число, полученное с помощью генератора случайных чисел, находим коэффициенты ряда Фурье для функции f (ф): fcn,fsn. В результате расчетная формула примет вид:
No / 2п \
и'г (1,ф) « 2П / e-£L12 ¡J 2ee со<ф-ф) cos(e sin(ф - ф)) - 1 I и (1,ф) dфdе,
2No
u'n (Г,ф)\г=1 = fcnanncosиф + fsnannsinиф,
где
No
f e-xxndx
ап =-и\- _ регуляризирующии множитель.
На рис. 1 представлены графики наблюдаемой и расчетной зависимости от угла ф нормальной производной от потенциала и'п .
Рис. 1. и'п наблюдаемая и и'п расчетная зависимость от ф
Теорема 2. Для задачи 2 ^ 1 граничный оператор преобразования Г21 имеет вид:
Г21 : I ^ I,
2п
f=U'-
L
21
а{ф--ф)
(2 cos £ sin (ф - Ф) — 1) f (ф)0,ф0,£. (2.4)
Функция I предполагается такой, что сходится числовой ряд
I + \Ьк I) <
k = 1
(2.5)
2п
ak = 2П °OS к'Ф ^ f (Ф)
1
о
2п
bk = 2П J sin кФ • f (Ф) d-Ф
о
— коэффициенты ее ряда Фурье.
Пример 2. По значениям потока u'n = 3 sin ф — | sin3ф на границе требуется восстановить потенциал u.
С точностью до произвольной постоянной C решение будет иметь вид:
3 1
u (г,ф)\г=1 = 4 sinф — 4 sin3ф.
e
Опишем алгоритм решения. По значениям наблюдаемой функции
f (ф)= (jsin ф - 3 sin 3ф^ (1 + se),
где S = 10-4, e — случайное число, полученное с помощью генератора случайных чисел, находим коэффициенты ряда Фурье для функции f (ф): fcn,fsn. Расчетная формула примет вид:
No / 2п \
u (1,ф) « J e-eL2i (J 2e£ соз(ф-ф) cos (е sin(ф - ф)) - 11 u'r (1,ф) ¿ф0е. о \о )
Применив интегральную формулу (2.3), придем к разложению неизвестного потенциала в ряд Фурье:
1 1
u (т,ф)\г=1 = У] fcnan -cos пф + fsnan -sin пф. ' 1 ^ п п
n=l
В итоге расчетная формула может быть записана в виде:
2 No
u (т,ф)\г = л ~ У^ fcn^n 1Cos пф + fsnan -sin пф.
1 1 ¿—^ п п
n=l
Рис. 2. Графики наблюдаемой и расчетной зависимости от угла ф потока и
Указанная картина (рис. 2) подтверждает, что решение приведенной задачи определяется с точностью до произвольной постоянной О.
Заключение
Применение метода операторов преобразования позволило найти решение задач, в которых требовалось определить краевое условие, в замкнутом виде. Получено аналитическое решение, что позволило моделировать стационарные потенциальные поля, например, гравитационные или магнитные поля.
Литература
[1] Баврин И.И. Операторный метод в комплексном анализе. М.: Прометей, 1991. 200 с.
[2] Парфенова Ю.А. Векторные операторы преобразования функций, гармонических в шаре // Вестник СамГУ. Естественнонаучная серия. 2010. № 76. С. 48-56.
Поступила в редакцию 18/IX/2012; в окончательном варианте — 18/IX/2012.
MODELING OF POTENTIAL FIELDS BY METHOD OF BOUNDARY TRANSFORMATION OPERATORS
© 2012 J.A. Parfenova2
In the article the method of boundary transformation operators is offered. It supplements classical methods and allows to find the solution of regional tasks in the closed look that opens new opportunities in the simulation of potential fields. Application of the developed method allowed to solve a number of problems of potential theory and apply to the research a package of applied MatLab programs.
Key words: modeling, potential field, boundary transformation operators.
Paper received 18/IX/2012. Paper accepted 18/IX/2012.
2Parfenova Julia Alexeevna ([email protected]), the Dept. of Mathematical Analysis, Penza State Pedagogical University, Penza, 440026, Russian Federation.
