Научная статья на тему 'Метод построения конформного отображения единичного круга на риманову поверхность'

Метод построения конформного отображения единичного круга на риманову поверхность Текст научной статьи по специальности «Математика»

CC BY
194
26
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
КОНФОРМНОЕ ОТОБРАЖЕНИЕ / РИМАНОВА ПОВЕРХНОСТЬ / АНАЛИТИЧЕСКАЯ ФУНКЦИЯ / УРАВНЕНИЕ ФРЕДГОЛЬМА / CONFORMAL MAPPING / RIEMANN SURFACE / ANALYTIC FUNCTION / FREDHOLM EQUATION

Аннотация научной статьи по математике, автор научной работы — Иваньшин Петр Николаевич

В статье приведен метод построения конформного отображения единичного круга на риманову поверхность (отображение с неоднолистным образом). В случае области, лежащей на римановой поверхности, построение конформного отображения единичного круга на нее сведено к решению интегрального уравнения. Дан вывод необходимых соотношений из формул Сохотского. Построен пример, иллюстрирующий метод для двулистно накрывающей плоскость римановой поверхности. Необходимое и достаточное условие того, что заданная на замкнутой кривой функция является граничным значением некоторой функции, аналитической в находящейся на римановой поверхности области, ограниченной данной кривой, примененное для отображения единичного круга на односвязную и однолистную область и обеспечивающее появление интегрального уравнения, следует несколько изменить. Для функции ϕ( z ) = ln (ζ( z )/ z ) можно выписать уравнения, аналогичные уравнениям для однолистной области, но на участках контура, ограничивающих область двулистности, необходимо поделить правую часть на три.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A Construction Method for Conformal Mapping of the Unit Disk onto a Riemann Surface

A construction method for approximate conformal mapping of the unit disk onto a Riemann surface (a map with a self-overlapping image) has been described. An example has been provided to illustrate the applicability of the method to conformal mapping of the unit disk onto a two-sheeted covering of the domain by a Riemann surface. The function construction is based on the approximate solution of the second kind Fredholm integral equation by reducing it to the finite linear equation system, so the construction is easily programmable. The necessary and sufficient condition for the function given on the closed curve to be the boundary value of some function analytic in the region on the Riemann surface bounded by the given curve is naturally somewhat different from that for one-connected and one-sheeted domains. We have applied this condition for a multiply-sheeted region. Let z (ζ) be the function that maps the unit disk onto a multiply-sheeted region conformally. For the function ϕ( z ) = ln (ζ( z )/ z ), we write the equations similar to that for one-connected and one-sheeted domains. Note that for our example with two-sheeted domain it is necessary to divide the right-hand side of our relations by 3 for the points on the contour sections bounding the nonunivalent region. The solution then repeats the steps of the one-sheeted domain situation for our case.

Текст научной работы на тему «Метод построения конформного отображения единичного круга на риманову поверхность»

2018, Т. 160, кн. 4 С. 738-749

УЧЕНЫЕ ЗАПИСКИ КАЗАНСКОГО УНИВЕРСИТЕТА. СЕРИЯ ФИЗИКО-МАТЕМАТИЧЕСКИЕ НАУКИ

ISSN 2541-7746 (Print) ISSN 2500-2198 (Online)

УДК 517.54

МЕТОД ПОСТРОЕНИЯ КОНФОРМНОГО ОТОБРАЖЕНИЯ ЕДИНИЧНОГО КРУГА НА РИМАНОВУ ПОВЕРХНОСТЬ

П.Н. Иваньшин

Казанский (Приволжский) федеральный университет, г. Казань, 420008, Россия

Аннотация

В статье приведен метод построения конформного отображения единичного круга на риманову поверхность (отображение с неоднолистным образом). В случае области, лежащей на римановой поверхности, построение конформного отображения единичного круга на нее сведено к решению интегрального уравнения. Дан вывод необходимых соотношений из формул Сохотского. Построен пример, иллюстрирующий метод для двулистно накрывающей плоскость римановой поверхности.

Необходимое и достаточное условие того, что заданная на замкнутой кривой функция является граничным значением некоторой функции, аналитической в находящейся на римановой поверхности области, ограниченной данной кривой, примененное для отображения единичного круга на односвязную и однолистную область и обеспечивающее появление интегрального уравнения, следует несколько изменить.

Для функции ф(г) = !а(^(г)/г) можно выписать уравнения, аналогичные уравнениям для однолистной области, но на участках контура, ограничивающих область двулистно-сти, необходимо поделить правую часть на три.

Ключевые слова: конформное отображение, риманова поверхность, аналитическая функция, уравнение Фредгольма

Введение

Известно много методов построения конформного отображения односвязной области на единичный круг и, наоборот, единичного круга на область комплексной плоскости. Рассмотрим сначала отображение области на круг. В случае метода оскуляции [1] происходит отображение области Б на подобласть единичного круга при помощи функций или 1пг, затем на большую подобласть и т. д., до тех пор, пока не получится отображение области на единичный круг. Развитием метода оскуляции является метод Грассмана [2], когда сначала ищется возможность отображения области на круг без сектора, потом применяется преобразование Мебиуса. Рассмотрим теперь метод [3]. Этот метод основан на том, что в логарифми-

ческом методе Кебе на определенном этапе на каждой итерации логарифмический образ области обычно не только лежит в левой полуплоскости, но и в некоторой горизонтальной полосе, то есть мнимая часть всех ее точек ограничена сверху и снизу. Соответствующее вещественно-аффинное преобразование отображает эту половину полосы на нормированную полуполосу {1ш[г] € (-п/2,п/2)}. В свою очередь, функция вшИ отображает эту полуполосу на всю левую полуплоскость. В результате гораздо большая часть полуплоскости, чем в логарифмическом методе Кебе, занята образом области до того, как полуплоскость отображается на единичный круг на последнем шаге итеративной процедуры. Методы Шварца -Кристоффеля для многоугольных областей [4, 5] работают для областей, которые

можно эффективно приблизить многоугольниками. Сложность этого метода заключается в поиске вспомогательных параметров для удачной первой итерации. Г.Т. Симм предложил решать задачу построения конформного отображения сведением к интегральному уравнению [6]. Им было предложено несколько линейных уравнений, но не был предложен алгоритм решения этих уравнений.

Существуют также методы построения конформного приближенного отображения единичного круга на односвязную область. Метод Теодорсена [7] применим только для звездных областей. Этот метод использует гармонически сопряженные функции для нахождения перепараметризации границы контура и последовательные итерации. Метод Форнберга [8] - итерационный метод построения полиномиального отображения - требует хорошего первого приближения. Все эти методы требуют хорошо подобранного начального приближения отображения. Итерационный метод Вегмана [9, 10] предлагает для численного построения конформного отображения решение вспомогательных краевых задач для построения соответствующих итераций. Алгоритм zipper [11] отображает единичный круг на многоугольную область. Этот алгоритм строит отображение как композицию конформных отображений на полуплоскости с разрезами, то есть опять представляет собой последовательность итераций.

Аналитический метод перепараметризации граничной кривой, предложенный впервые в [12], строго обоснованный в [13] и примененный в настоящей работе, имеет следующие преимущества перед перечисленными методами: решение имеет полиномиальный вид, что облегчает приложения решения к другим задачам [14, 15], решение не использует вспомогательные построения (см. [8, 16-18]) и итерации (см. [10, 19, 20]). Аналитический метод основан на приближенном решении интегрального уравнения Фредгольма 2-го рода сведением к системе линейных уравнений и легко программируется. Интегральное уравнение Фредгольма - следствие необходимого и достаточного условия того, чтобы заданная на контуре функция была граничным значением функции, аналитической в соответствующей области [21]. Это соотношение дано в [21] для произвольной однолистной области и основано на формуле Сохоцкого и интегральной формуле Коши. Здесь для случая самопересекающейся границы и многолистной области задача тоже сводится к интегральному уравнению Фредгольма, для вывода которого применяется измененная интегральная формула Коши для разных подобластей. Метод демонстрируется на примере области, расположенной на двулистной римановой поверхности.

В случае области, лежащей на римановой поверхности, сведение построения конформного отображения к интегральному уравнению возможно аналогично случаю однолистной области. При этом необходимое и достаточное условие того, что заданная на замкнутой кривой функция является граничным значением некоторой функции, аналитической в находящейся на римановой поверхности области, ограниченной данной кривой, примененное в для отображения единичного круга на односвязную и однолистную область и обеспечивающее появление интегрального уравнения, следует несколько изменить. В работе построен пример подобного отображения единичного круга на двулистную область.

Пусть Z(z) - искомое отображение. Функция ф(г) является аналитической в области D тогда и только тогда, когда выполнено соотношение

1 1 г 1пГ7~^+ i(0(T) " 'AIgz(T))

1 -iai4-\ i+w 1 \z(t

1п

\ + чт - argzO) = П-J 1 z(T' \z(T ) .zW-z'(T > *+

1 1

1 г 1пГТ^ + i(6(T) - argz(r))

+ — \z(T) \ , ,-Г^-z'(T) dr. (1)

3ni J z(r) - z(t) W W

Здесь в(Ь) - функция, выражающая зависимость полярного угла в единичной окружности от параметра Ь контура Ь = Ь, ограничивающего область Б. Рассмотрим теперь мнимую часть последнего соотношения, обозначим разность в(Ь) — argz(t) через д(Ь) и получим следующее уравнение Фредгольма второго типа на д(Ь):

2п

II >

д(arg(z(т) — г(Ь)))

дт

¿т+

2п

д(ь) =

+ -^У 1п \г(т

,д[1п\г(т) — г(1)\]

дт

¿т, Ь € Ь\,

2п

>

д(а^(т) — г(г)))

дт

¿т+

2п

д[1п\г(т) — г(г)\] 1 дт

¿т, Ь € Ь2.

(2)

Схема решения после построения основной системы уравнений для данной мно-голистной области повторяет шаги, представленные в [13] для односвязной области Б С С. Если решение интегрального уравнения представляет собой ряд Фурье

д(Ь) = ^^ ап сов(пЬ) + /Зп в[п(пЬ),

(3)

то приближенное решение этого уравнения можно представить в виде

м

д(Ь) = ^^ ап сов(пЬ) + /Зп в[п(пЬ).

(4)

п=1

Здесь число М нужно брать настолько большим, чтобы функция в(Ь) = д(Ь) + + arg(z(t)) была монотонно возрастающей. Решаем полученную конечную систему линейных уравнений на ап, [Зп и восстанавливаем д(Ь).

Построение конформного отображения единичного круга на область, расположенную на римановой поверхности

Рассмотрим сначала решение интегрального уравнения Фредгольма второго рода

2п

х(ь) = ! с(т,г)х (т) ¿т + У(г). (5)

о

Здесь X(Ь) - неизвестная функция, 0(т,Ь) - 2п-периодическая функция по обеим переменным, У (Ь) - 2п-периодическая функция по Ь. Соответственно, решение тоже 2п-периодично, то есть представимо в виде ряда Фурье. Заметим, что интегральное уравнение вида (5) появляется, например, при решении задачи Дирихле для плоской области методом потенциалов [22].

Пусть интегральное уравнение однозначно разрешимо в линейном пространстве интегрируемых функций. Сведем решение уравнения к решению бесконечной системы линейных уравнений. Тогда приближенное решение интегрального уравнения будет иметь вид конечной частной суммы ряда Фурье, так как бесконечная

система заменяется усеченной конечной системой. Оказывается, что для усеченной системы линейных уравнений достаточно высокого порядка решение тоже будет существовать и сходиться к решению интегрального уравнения.

Напомним основное утверждение, обеспечивающее существование конформного приближенного отображения [13].

Лемма. Пусть существуют такие числа з, р > 1 и константа и > 0, что

д^+рС(т,г)

дР дт р

<и,

и функция У(£) обладает ограниченной второй производной: \У< Т. Тогда приближенное решение однозначно разрешимого интегрального уравнения Фред-гольма второго рода (5), может быть сведено к решению конечной линейной системы с ошибкой, которую можно оценить как 0(1/N2). Здесь N - порядок конечной системы.

В случае области, лежащей на римановой поверхности, сведение построения конформного отображения к интегральному уравнению также возможно. При этом необходимое и достаточное условие того, что заданная на замкнутой кривой функция является граничным значением некоторой функции, аналитической в находящейся на римановой поверхности области, ограниченной данной кривой, примененное в [13] и обеспечивающее появление интегрального уравнения, следует несколько изменить. Рассмотрим это на примере двулистной области.

Пусть Ь = дВ - граница двулистной области В, тогда для интеграла типа Коши

1 Г Ф(т)

йт

определены Ф+(г), г € В, Ф-(г), г € В [21]. Пусть функция ф аналитична в области В. Обозначим часть контура, не накладываемую на сам контур, через Ь. Участки контура, ограничивающие область неоднолистности, обозначим через Ь2 (рис. 1).

Для точек г в зависимости от их расположения относительно В получим

йт

2пг, г € В+,

0, г € В-,

п г, г € Ьъ

3пг, г € Ь2.

т — г

т — г

Следовательно,

Иш

Л гт ¿т—¿т

2пг ] т — Ь 2пг ] т — Ь

ь ь

Иш

х—Л-

[Фт! ¿т — ¿„т

2пг ,] т — Ь 2пг ,] т — Ь

ьь

Ф+(Ь) — ф(Ь),

= Ф-(Ь),

Ы^т — м/-—¿т = т — т2 Ь€ Ь1,

ьь

^ ! Ф—1 ¿т — И ! ¿т = ™ - ЫШ ' € Ь2.

ьь

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Для Ь € Ь2 получим Ф+ (Ь) — ф(Ь) = Ф-(Ь) = Ф(Ь) — 3ф(Ь)/2. Тогда в терминологии [21]

ф+(ь)=— 1 ¿т, ь € Ь2,

ь

ф-(ь)=— 3м+Ъ!^ ¿т, Ь€ Ь2.

ь

Следовательно,

,, , 2 1 Г ф(т) 1 1 Г ф(т) ,

Ф(Ь) = ЗъТг .1 т—г ¿т =ьпЛ—г ¿т' Ь€ Ь2.

ьь

С (г)

1п- - аналитическая функция в Б в том и только том случае [21], когда

z

1 1 1п рт1-^ + г(в(т) — агёг(т)) 1п Ш + т) — аГёФ)) = Пг! \(г) — г(Ь)-z,(т) ¿т+

Ь1

1 г 1пГ7"^ + г(в(т) — 'аг^(т))

+ Ы — ^-) ¿т, (6)

Ь2

где в( ) - зависимость полярного угла в на границе единичного круга от исходного параметра Ь контура Ь - границы Б.

Значит, для ф(г) = 1п(С(г)/г) получим уравнения, аналогичные уравнениям из [13], но на участках контура, ограничивающих область неоднолистности, необходимо поделить правую часть на 3 согласно формуле (2).

Далее, следуя рассуждениям, проведенным в [13] для односвязной области: 1) составляем бесконечную систему относительно коэффициентов Фурье искомой функции; 2) решаем конечную систему в соответствии с приведенной леммой; 3) найдя соответствие между параметрами Ь и в, приближенно восстанавливаем г(С) в виде полинома, перепараметризуя граничную кривую.

Затем строим приближенное решение интегрального уравнения по формулам

(3), (4).

Приближенное решение интегрального уравнения по лемме сводится к решению системы линейных уравнений

(ВВ АвБв) ® = (0) • (7)

где

/аЛ (¡Зг

а = I ... I , в = I •

\ам ) \вм /

Векторы

, = Р) о = (й

\1м/ \9м /

правой части системы согласно формуле (2) состоят из элементов

2п

11 = ( 1т ¿о + 1т \ + ~2 /1п\г(т )\ йт [ Ь(т,г)оов(Н) йг+

3 к ^ 0 Ьг

2п

+ -П ! 1п\г(т)\ йт ! Ь(т, г) соа(и) йг,

2п

91 = ( ^ Ке 3 + Ке 1 + 1 / !п \г(т)\ йт [ Ь(т,г)вт(Н) йг+

3 к ^ 0 Ьг

2п

+ П J 1п\г(т)\ йт I Ь(т,г)в!п(и) йг, I = 1,..., м.

з к 0

2п

1

0 Ь2

Блок-матрицы ВВ, ВВ, ВВ, ВВ порядка М согласно формуле (2) состоят из элементов

2п

1

ВВ = I 81п--2 У сов(пт) йт ! К(т,г)еов(О) йг—

2п \М

— ^ J сов(пт) йт ! К(т, г) еов(О) йг I .

0 Ь2 / 1,п=1

2п

1

ВВ = I —1 ! яш(пт) йт I К(т,г)еов(И) йЬ—

0 Ьг

2п

— ^ J эт(пт) йт ! К(т,г)еов( 1г) йг

Ь2 / 1,п=1

2п X М

2п

2п

1

В А = I —1 J соэ( пт) ¿^ К (т,1)вт(и) ¿Ь—

о ь1

2п

— 3—2 ^ сов(пт) ¿т ! К(т,Ь)5т(1Ь) ¿Ь

1,п=1

2п X м

2п

2п

В В = I $1п--вт(пт) ¿т ! К (т,Ь)5т(1Ь) ¿Ь—

2п

— 3—2 ^ вш(пт) ¿т ! К(т,Ь)вт(и) ¿Ь

1,п=1

Здесь §1п - дельта-функция Кронекера.

Пример. Построим приближенное отображение единичного круга на область, лежащую на двулистной римановой поверхности и ограниченную самопересекающейся кривой

г(Ь) = —0.05е4й — 0.55вш + еи — 0.05е-2й + 0.05в-Ш, Ь € [0, 2п].

Решим систему, соответствующую уравнениям (2). Получим для полинома степени 35, отображающего единичный круг на область, ограниченную заданной самопересекающейся кривой и лежащей на римановой поверхности, картину, представленную на рис. 2.

Полином, дающий отображение единичного круга на неоднолистную область с заданной границей, имеет вид

— 0.00164298С35 — 0.00136456С34 — 0.0011142С33 — 0.00120001С32 —

— 0.00138477С31 — 0.0011375С30 — 0.000504466С29 — 0.000313748С28 — — 0.000801963С27 — 0.00 1 26 1 93С26 — 0.000746909С25 + 0.000038404С24—

— 0.000517881С23 — 0.0020071С22 — 0.00287575С21 — 0.00138472С20 —

— 0.000646775С19 — 0.00278893С18 — 0.005697С17 — 0.00624325С16 —

— 0.000982626С15 — 0.00192618С14 — 0.00649747С13 — 0.0126356С12—

— 0.00905622С11 + 0.00769869С10 — 0.00675355С9 — 0.00999969С8 —

— 0.0398029С7 + 0.00415075С6 + 0.0478242С5 — 0.0585483С4+

+ 0.0140739С3 — 0.507787С2 + 1.01976С.

Наибольшее отклонение исходной граничной кривой от границы области, полученной построенным отображением единичного круга на двулистную область, равно 0.02 .

На рис. 2 видно, что образы окружностей радиусов 1 и 0.9 имеют самопересечения, а образ окружности радиуса 0.8 его не имеет.

Рис. 2. Аппроксимация двулистного отображения: слева - граница неоднолистной области при построенном полиномиальном отображении степени 35 и образ концентрических

окружностей, справа - исходная кривая

Литература

1. Conway J.B. Functions of One Complex Variable. - N. Y.; Berlin: Springer-Verlag, 1978. -XIV, 322 p. - doi: 10.1007/978-1-4612-6313-5.

2. Grassmann E. Numerical experiments with a method of successive approximation for conformal mapping // Z. Angew. Math. Phys. - 1979. - V. 30, No 6. - P. 873-884. - doi: 10.1007/BF01590486.

3. Porter R.M. An accelerated osculation method and its application to numerical conformal mapping // Complex Var. Theory Appl. - 2003. - V. 48, No 7. - P. 569-582. - doi: 10.1080/0278107031000110892.

4. Driscoll T.A., Trefethen L.N. Schwarz-Christoffel Mapping. - Cambridge: Cambridge Univ. Press, 2002. - XVI, 132 p. - doi: 10.1017/CB09780511546808.

5. Driscoll T.A., Vavasis S.A. Numerical conformal mapping using cross-ratios and Delaunay triangulation // SIAM J. Sci. Comput. - 1998. - V. 19, No 6. - P. 1783-1803. - doi: 10.1137/S1064827596298580.

6. Symm G.T. An integral equation method in conformal mapping // Numer. Math. - 1966. -V. 9, No 3. - P. 250-258. - doi: 10.1007/BF02162088.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

7. Henrici P. Applied and Computational Complex Analysis. V. 3: Discrete Fourier analysis, Cauchy integrals, construction of conformal maps, univalent functions - . N. Y.: John Wiley & Sons, 1986. - 637 p.

8. Fornberg B. A numerical method for conformal mappings // SIAM J. Sci. Stat. Comput. -1980. - V. 1, No 3. - P. 386-400. - doi: 10.1137/0901027.

9. Wegman R., Murid A.H.M., Nasser M.M.S. The Riemann-Hilbert problem and the generalized Neumann kernel // J. Comput. Appl. Math. - 2005. - V. 182, No 2. -P. 388-415. - doi: 10.1016/j.cam.2004.12.019.

10. Wegman R. An iterative methods for conformal mapping //J. Comput. Appl. Math. -1986. - V. 14, No 1-2. - P. 7-18. - doi: 10.1016/0377-0427(86)90128-7.

11. Marshall D.E., Rohde S. Convergence of a variant of the zipper algorithm for conformal mapping // SIAM J. Numer. Anal. - 2007. - V. 45, No 6. - P. 2577-2609.

12. Shirokova E.A. On approximate conformal mapping of the unit disk on an simply connected domain // Russ. Math. - 2014. - V. 58, No 3. - P. 47-56. - doi: 10.3103/S1066369X14030050.

13. Shirokova E.A., Ivanshin P.N. Approximate conformal mappings and elasticity theory // J. Complex Anal. - 2016. - V. 2016. - Art. 4367205, P. 1-8. - doi: 10.1155/2016/4367205.

14. Abzalilov D.F., Shirokova E.A. The approximate conformal mapping onto simply and doubly connected domains // Complex Var. Elliptic Equations. - 2017. - V. 62, No 4. -P. 554-565. - doi: 10.1080/17476933.2016.1227978.

15. Abzalilov D.F., Shirokova E. A. The approximate conformal mapping onto multiply connected domains // Probl. Anal. Issues Anal. 2019. - V. 8, No 1. - P. 3-16. - doi: 10.15393/j3.art.2019.5050.

16. DeLillo Th.K., Elcrat A.R. A Fornberg-like conformal mapping method for slender regions //J. Comput. Appl. Math. - 1993. -V. 46, No 1-2. - P. 49-64. - doi: 10.1016/0377-0427(93)90286-K.

17. Warschawski S.E. On conformal mapping of variable regions // Nat. Bur. Stand. Appl. Math. Ser. - 1952. - V. 18. - P. 175-187.

18. Li B.C., Syngellakis S. Numerical conformal mapping based on the generalised conjugation operator // Math. Comput. - 1998. - V. 67, No 222. - P. 619-639.

19. Hough D.M., Papamichael N. An integral equation method for the numerical conformal mapping of interior, exterior and doubly-connected domains // Numer. Math. - 1983. -V. 41, No 3. - P. 287-307. - doi: 10.1007/BF01418327.

20. Kress R. Inverse problem and conformal mapping // Complex Var. Elliptic Equations. -2012. - V. 57, No 2-4. - P. 301-316. - doi: 10.1080/17476933.2011.605446.

21. Gakhov F.D. Boundary Value Problems. - Pergamon, 1966. - 584 p.

22. Tricomi F.G. Integral equations. - N. Y.: Dover Publ. Inc., 1985. - 256 p.

Поступила в редакцию 13.09.17

Иваньшин Петр Николаевич, кандидат физико-математических наук, доцент кафедры геометрии

Казанский (Приволжский) федеральный университет

ул. Кремлевская, д. 18, г. Казань, 420008, Россия E-mail: [email protected]

ISSN 2541-7746 (Print)

ISSN 2500-2198 (Online)

UCHENYE ZAPISKI KAZANSKOGO UNIVERSITETA. SERIYA FIZIKO-MATEMATICHESKIE NAUKI

(Proceedings of Kazan University. Physics and Mathematics Series)

2018, vol. 160, no. 4, pp. 738-749

A Construction Method for Conformal Mapping of the Unit Disk onto a Riemann Surface

P.N. Ivanshin

Kazan Federal University, Kazan, 420008 Russia E-mail: [email protected]

Received September 13, 2017 Abstract

A construction method for approximate conformal mapping of the unit disk onto a Riemann surface (a map with a self-overlapping image)has been described. An example has been provided to illustrate the applicability of the method to conformal mapping of the unit disk onto a two-sheeted covering of the domain by a Riemann surface.

The function construction is based on the approximate solution of the second kind Fredholm integral equation by reducing it to the finite linear equation system, so the construction is easily programmable.

The necessary and sufficient condition for the function given on the closed curve to be the boundary value of some function analytic in the region on the Riemann surface bounded by the given curve is naturally somewhat different from that for one-connected and one-sheeted domains. We have applied this condition for a multiply-sheeted region.

Let z(Z) be the function that maps the unit disk onto a multiply-sheeted region conformally. For the function 0(z) = ln(Z(z)/z), we write the equations similar to that for one-connected and one-sheeted domains. Note that for our example with two-sheeted domain it is necessary to divide the right-hand side of our relations by 3 for the points on the contour sections bounding the nonunivalent region.

The solution then repeats the steps of the one-sheeted domain situation for our case.

Keywords: conformal mapping, Riemann surface, analytic function, Fredholm equation

Figure Captions

Fig. 1. The two-sheeted domain.

Fig. 2. Approximation of the two-sheeted mapping: on the left - the border of the self-overlapping domain in case of the constructed 35-degree polynomial mapping and the image of concentric circles, on the right - the source curve.

References

1. Conway J.B. Functions of One Complex Variable. New York, Berlin, Springer-Verlag, 1978. XIV, 322 p. doi: 10.1007/978-1-4612-6313-5.

2. Grassmann E. Numerical experiments with a method of successive approximation for conformal mapping. Z. Angew. Math. Phys., 1979, vol. 30, no. 6, pp. 873-884. doi: 10.1007/BF01590486.

3. Porter R.M. An accelerated osculation method and its application to numerical conformal mapping. Complex Var. Theory Appl., 2003, vol. 48, no. 7, pp. 569-582. doi: 10.1080/0278107031000110892.

4. Driscoll T.A., Trefethen L.N. Schwarz-Christoffel Mapping. Cambridge, Cambridge Univ. Press, 2002. XVI, 132 p. doi: 10.1017/CB09780511546808.

5. Driscoll T.A., Vavasis S.A. Numerical conformal mapping using cross-ratios and De-launay triangulation. SIAM J. Sci. Comput., 1998, vol. 19, no. 6, pp. 1783-1803. doi: 10.1137/S1064827596298580.

6. Symm G.T. An integral equation method in conformal mapping. Numer. Math., 1966, vol. 9, no. 3, pp. 250-258. doi: 10.1007/BF02162088.

7. Henrici P. Applied and Computational Complex Analysis. Vol. 3: Discrete Fourier analysis, Cauchy integrals, construction of conformal maps, univalent functions. New York, John Wiley & Sons, 1986. 637 p.

8. Fornberg B. A numerical method for conformal mappings. SIAM J. Sci. Stat. Comput., 1980, vol. 1, no. 3, pp. 386-400. doi: 10.1137/0901027.

9. Wegman R., Murid A.H.M., Nasser M.M.S. The Riemann-Hilbert problem and the generalized Neumann kernel. J. Comput. Appl. Math., 2005, vol. 182, no. 2, pp. 388-415. doi: 10.1016/j.cam.2004.12.019.

10. Wegman R. An iterative methods for conformal mapping. J. Comput. Appl. Math., 1986, vol. 14, nos. 1-2, pp. 7-18. doi: 10.1016/0377-0427(86)90128-7.

11. Marshall D.E., Rohde S. Convergence of a variant of the zipper algorithm for conformal mapping. SIAM J. Numer. Anal., 2007, vol. 45, no. 6, pp. 2577-2609.

12. Shirokova E.A. On approximate conformal mapping of the unit disk on an simply connected domain. Russ. Math., 2014, vol. 58, no. 3, pp. 47-56. doi: 10.3103/S1066369X14030050.

13. Shirokova E.A., Ivanshin P.N. Approximate conformal mappings and elasticity theory. J. Complex Anal., 2016, vol. 2016, art. 4367205, pp. 1-8. doi: 10.1155/2016/4367205.

14. Abzalilov D.F., Shirokova E.A. The approximate conformal mapping onto simply and doubly connected domains. Complex Var. Elliptic Equations, 2017, vol. 62, no. 4, pp. 554565. doi: 10.1080/17476933.2016.1227978.

15. Abzalilov D.F., Shirokova E. A. The approximate conformal mapping onto multiply connected domains. Probl. Anal. Issues Anal., 2019, vol. 8, no. 1, pp. 3-16. doi: 10.15393/j3.art.2019.5050.

16. DeLillo Th.K., Elcrat A.R. A Fornberg-like conformal mapping method for slender regions. J. Comput. Appl. Math., 1993, vol. 46, nos. 1-2, pp. 49-64. doi: 10.1016/0377-0427(93)90286-K.

17. Warschawski S.E. On conformal mapping of variable regions. Nat. Bur. Stand. Appl. Math. Ser., 1952, vol. 18, pp. 175-187.

18. Li B.C., Syngellakis S. Numerical conformal mapping based on the generalised conjugation operator. Math. Comput., 1998, vol. 67, no. 222, pp. 619-639.

19. Hough D.M., Papamichael N. An integral equation method for the numerical conformal mapping of interior, exterior and doubly-connected domains. Numer. Math., 1983, vol. 41, no. 3, pp. 287-307. doi: 10.1007/BF01418327.

20. Kress R. Inverse problem and conformal mapping. Complex Var. Elliptic Equations, 2012, vol. 57, nos. 2-4, pp. 301-316. doi: 10.1080/17476933.2011.605446.

21. Gakhov F.D. Boundary Value Problems. Pergamon, 1966. 584 p.

22. Tricomi F.G. Integral Equations. New York, Dover Publ. Inc., 1985. 256 p.

/ Для цитирования: Иваньшин П.Н. Метод построения конформного отображения ( единичного круга на риманову поверхность // Учен. зап. Казан. ун-та. Сер. Физ.-ма-\ тем. науки. - 2018. - Т. 160, кн. 4. - С. 738-749.

/ For citation: Ivanshin P.N. A construction method for conformal mapping of the unit ( disk onto a Riemann surface. Uchenye Zapiski Kazanskogo Universiteta. Seriya Fiziko-\ Matematicheskie Nauki, 2018, vol. 160, no. 4, pp. 738-749. (In Russian)

i Надоели баннеры? Вы всегда можете отключить рекламу.