Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Методы математической физики и задачи гидроаэродинамики

..pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
14.42 Mб
Скачать

Эта теорема указывает, в частности, что не существует ана­ литической функции, осуществляющей конформное отображе­ ние конечной области на кривую или точку.

Принцип сохранения границ. Если функция w=w {z) кон­ формно отображает область D на область Д, то эта функция отображает взаимно однозначно границу То области D на гра­ ницу Гд области Д с сохранением ориентации.

Теорема существования и единственности. Если границы од­ носвязных областей D u b . состоят более чем из одной точки, то существует и притом одна функция w—w(z), конформно отображающая область D на область Д так, что w(zo)=Wot

w (zrD )= w гд.

где z0,

w0,

z го» ^гд—заданные числа (точки),

причем Z o ^ D ,

W 0œ A ,

Z T D

е Г д, шгдс:Гд.

З а м е ч а н и я . Условия

w (z0) — W О, W (Zrz) ) = w c д ,

задающие соответствие пары внутренних точек z0, w0 и пары

граничных точек zjb, ^гд» обеспечивают единственность отобра­ жения. Без учета этих условий существует бесконечно много функций w —w(z), конформно отображающих односвязную область D с границей, состоящей более чем из одной точки, на однос^язную область Д с аналогичной границей. В качестве ус­ ловий, обеспечивающих единственность конформного отображе­ ния, можно использовать также следующие

w (zQ)= w 0, aigw' (zQ)= aQ

или

w (z-pui)=WrM, i = 1 »2,3.

В первом случае задается соответствие пары внутренних точек Zo, wo и угол поворота направлений (векторов) в точке ZQ. Во

втором — соответствие трех пар граничных точек.

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

Из условия теоремы следует также, что она относится к односвязным областям. В случае неодносвязных областей воп­ рос существования конформного отображения является более сложным. Даже для простейших двусвязных областей не всегда существует конформное отображение.

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

 

Конформное отображение

с помощью

степенной функции.

Степенная функция

 

 

 

w —z nt

(6.117)

где

2 — натуральное число,

принадлежит к неоднолистным

(многолистным) или не взаимно однозначным функциям.

 

Найдем область однолистности функций (6.117), однознач­

ность которой при любых z очевидна (как

очевидна и много-

 

 

 

Л г т—*

листность, поскольку/обратная функция z —y w имеет п раз-

а)

5)

В)

Рис. 6.21

личных значений в любой фиксированной точке w^O). Одно­ листность нарушается там, где одному значению функции w соответствуют, по крайней мере, два значения аргумента, т. е. там, где

z4=z%, z xф ^ г.

(6.118)

Отсюда, полагая Zi=rie**', на основании правила вычисления корня находим, что

£ = 0 , 1, п — 1.

 

При Æ=0 и k— \ соответственно получаем

 

/—

(6.119)

z 2= z lt z 2= z xe п,

т. е. однолистность нарушается в точке 2г, имеющей

тот же

модуль, что и z u но аргумент, отличающийся на /л (точки А, В на рис. 6.21, а). Следовательно, однолистность имеет место для

всех точек z таких, что

argzi< argz< argzi + 2л/п, т. е. для

всех точек z, лежащих

внутри угла с центром в О раствора

2njn (рис. 6.21, а).

 

(6.117)

Если положить z = r e /4> о>=ре‘°, то из соотношения

имеем

 

'

 

р = г » , б=Я<р,

(6. 120)

откуда видно,

что внутренности

угла 2я/я

соответствует вся

плоскость w с разрезом по действительной

полуоси

Imtü=0,

Rete^O (при

движении точки

z вдоль дуги АВ радиуса г

угол ф изменяется от 0 до 2я/я

и точка w совершает

полный

обход вокруг начала координат вдоль окружности радиуса р = = г п, так как угол 0 изменяется от 0 до 2я).

На сторонах угла АОВ однолистность нарушается и точки

А, В с одинаковым модулем

перейдут в

одну и ту же точку т=гп, лежащую

на действительной оси

(точки А', В' на рис. 6.21, б). Если же в плоскости w провести

разрез вдоль положительной части действительной

оси

(рис.

6.21, в),

то, поскольку на верхнем берегу разреза 0 = 0 ,

а

на

нижнем

0 = 2 я , преобразование

(6.117) станет однолистным

и

на сторонах угла АОВ (точки

Af, В' перестанут

совпадать).

Итак, преобразование (6.117) устанавливает взаимно однознач­ ное соответствие между внутренностью угла раствора 2л/л (включая стороны) и всей плоскостью ш, разрезанной вдоль положительной части действительной оси.

Но функция w = z n определена и вне угла АОВ, причем сле­ дующему углу ВОС с тем же раствором 2я/я вновь соответст­ вует вся плоскость w [см. формулы (6.120)], которая, однако, «занята» образами точек угла АОВ. Чтобы преодолеть это за­ труднение, поступим следующим образом. Возьмем я экземпля­ ров плоскостей w, пронумеруем их и разрежем каждый эк­ земпляр вдоль положительной части действительной оси. Да­ лее склеим нижний берег разреза на плоскости 1 с верхним берегом разреза на плоскости 2, нижний берег разреза на пло­ скости 2 с верхним берегом разреза на плоскости 3 и т. д. и, наконец, нижний берег разреза на плоскости я склеим с верх­ ним берегом разреза плоскости 1. Теперь каждому из я углов

раствора 2я/я

(АОВ,

ВОС,...) соответствует свой лист плоско­

сти w, разрезанной

вдоль полуоси м>0, о = 0 . В результате

преобразование

w—zn при любом z станет однолистным на

построенной я-листной поверхности, называемой римановой по­ верхностью функции w —zn.

Обратим теперь внимание, что функция (6.117) отображает

внутренность

угла 2я/я на всю плоскость с разрезом и > 0,

i)= 0

(на угол 2я), т. е. увеличивает углы с вершиной в нуле в

п раз

[см. также формулу (6.120)]. Следовательно, преобразо­

вание

(6.117)

не конформно в нуле. Это можно было ожидать,

поскольку производная - ^ - = n z n~l —о при z = 0 и я ^ 2 .

dz

Преобразование (6.117) позволяет отображать углы разно­ го раствора и различным образом расположенные на плоскости друг на друга. Проиллюстрируем это на примере, в котором

требуется конформно отобразить угол раствора п[п с вершиной в точке ЬфО (рис. 6.22) на верхнюю полуплоскость. Последо­ вательность решения этой задачи показана на рис. 6.22, откуда следует, что функция, реализующая • требуемое отображение, имеет вид

гv= z% = (zl e~/a)n= ( z — b)ne-/rte.

Функция Жуковского. Рассмотрим преобразование

® = Ÿ ( z + 7 ~ ) -

(6Л21)

Функцию (6.121) называют функцией Жуковского. Однознач­ ность функции (6.121) очевидна, как очевидна и ее многолистность, поскольку обратная функция

z = w - \ - y rw l —1

двузначна (квадратный корень имеет два различных значения для каждого фиксированного w) и, следовательно, риманова поверхность функции Жуковского состоит из двух листов. Най­ дем область однолистности функции Жуковского. Однолист­ ность нарушается в точках z\^=z2i в которых

Отсюда

 

 

 

 

(

1----- )= 0 ,

 

\

* 1*2

)

и так как

то ziz2= l , т. е.

 

UiI |*2l=l» arg«!+argz2=0.

Таким образом, однолистность имеет

место для точек |z |< l ,

лежащих внутри единичного

круга,

или для точек |z |> l , ле­

жащих вне единичного круга

(в точках единичной окружности

\z\ = 1, симметричных относительно оси Ох,

однолистность на­

рушается). Будем рассматривать внешность

единичного

круга

|а |> 1 . Полагая r = r e i(P, w = u + iv ,

из выражения (6.121)

полу­

чаем

 

 

 

u-\-i'o=-^-\r

-j——

 

 

и, следовательно, ,

Рис. 6.23

Отсюда видно, что окружностям r= co n st> l соответствуют эл­ липсы с полуосями

При

увеличении

г эллипсы «округляются» (при

больших г

полуоси

а » 6 » г /2 ) ,

при уменьшении — «вытягиваются» и при

г = 1

вырождаются в отрезок [— 1, 1], так как при г = 1

из фор­

мул

(6.122) следует, что

 

 

 

 

 

 

 

 

a=cos<p>

-п=0.

 

 

 

Очевидно, внешности

единичного

круга

при преобразовании

(6.121)

соответствует вся плоскость w с выброшенным

отрез­

ком [— 1, 1], а самой

окружности

12г| = 1 соответствует

отрезок

[ - 1 , 1).

полном

обходе окружности г= 1

(рис. 6.23,

а), когда

При

угол <р изменяется от 0 до 2п, отрезок [— 1, 1] проходится дваж­ ды (рис. 6.23, б), т. е. однолистность нарушается, что отмеча­ лось выше. Если же сделать разрез вдоль отрезка [— 1, 1] и от­ нести (рис. 6.23, б) верхнему берегу разреза верхнюю полу­ окружность АВС, а нижнему берегу — нижнюю полуокружность CDA, то преобразование Жуковского станет однолистным в об­

ласти

| z | ^ l . Таким образом, преобразование (6.121) устанав­

ливает

взаимно однозначное

соответствие

между внешностью

единичного круга (включая

окружность

г = 1 ) и всей пло-

скостью

w, разрезанной

вдоль

отрезка [— 1, 1]. Функция

Жу­

ковского

определена и

внутри

единичного

круга,

причем

[см.

формулы

(6.122)] внутренности

единичного

круга

вновь соот­

ветствует вся плоскость w, которая, однако, «занята» образами точек внешности единичной окружности. Чтобы преодолеть это

затруднение, как и в случае функции

w = z n, возьмем два эк­

земпляра плоскости w, перенумеруем

их и проведем

на

каж­

дом из них разрез вдоль отрезка [— 1,

1]. Затем склеим нижний

берег разреза на плоскости 1 с верхним берегом разреза

на

плоскости 2, а нижний берег разреза

на плоскости 2 — с верх­

ним берегом разреза на плоскости 1.

Теперь

внешности

или

внутренности единичной окружности

соответствует

свой

лист

плоскости w, разрезанный вдоль отрезка [—1,

1] и преобразо­

вание Жуковского стало однолистным при любом z на постро­ енной двулистной поверхности, которая называется римановой поверхностью функции Жуковского. Применение функции Жу­

ковского в задачах гидроаэродинамики рассмотрено ниже. Рассмотрим теперь гидромеханический смысл производной

аналитической функции. Пусть течение идеальной несжимае­ мой жидкости в области D плоскости XOY стационарное и без­ вихревое и в этой области нет стоков и источников. Обозначим через V==uï+v] поле скоростей частиц жидкости. Тогда rot V— = 0 , div F = 0 . Следовательно, поле V потенциально и для его потенциала ф(дг, у) имеем

Vr=grad«>, div V = d iv grad <р=Дср=0.

Рассмотрим функцию ty(x, у), связанную с <р(х, у) условия­ ми Коши — Римана

 

д?

d<|j

д<р

dty

 

дх

ду *

ду

дх

Дифференцируя

первое равенство по у, второе — по х

 

д2<р

 

д2у

 

 

дхду

ду2

Ьудх

дх2

и учитывая, что

- ^

дудх

, получаем

 

дхду

 

 

 

 

дЦ> .

дЦ> _

п

 

 

дх2

ду2

 

т. е. ф(лг, у) — гармоническая функция.

Гармоническая функция ф, связанная с гармонической функцией <р условиями Коши — Римана,, называется сопряжен­ ной функции <р.

Заметим, что согласно

условиям Коши — Римана

grad<p и

grad-ф ортогональны:

 

 

 

 

. ( p . d T,gtad«—

аф

аф

дф

<?Ф_Q

ду

дх

дх

ду

 

следовательно, в каждой

точке области

D проходящие через

нее линии уровня функций ( риф также ортогональны. Но тог­ да линия уровня /: ф(я, #)==const (ортогональная gradгр) ка­

сается вектора

F = grad

ср (рис. 6.24) и, значит, линия / явля­

ется векторной линией поля V. Ее назы­

 

вают в гидродинамике линией тока. (При

 

установившемся течении траектории час­

 

тиц жидкости и линии тока совпадают.)

 

Функция ф_(л:, у)

называется

функцией

 

тока поля V.

 

 

 

 

 

 

 

 

Аналитическая функция

 

 

 

w №)=<?(*, y )+ iù (x ,у),

(6.123)

 

где z = x + iy , ф(лс,

у)

и

ф(х,

у ) — соот­

 

ветственно

потенциал

и

функция тока

 

поля скоростей

V,

называется комплексным потенциалом тече

ния.

 

производную комплексного потенциала

 

Рассмотрим

 

 

dz

= * L J L

_ IЛ.=а _ iv.

(6.124)

 

 

дх

 

дх

дх

ду

 

Если ввести

понятие

комплексной

скорости Кк движущейся

жидкости:

 

 

 

 

 

 

 

 

 

 

 

 

VK=V t^=u -\-iv,

(6.125)

где V— | V\ — модуль

вектора

скорости, а 0 — угол его

накло­

на к оси Ох, то, как следует из

(6.124)

и (6.125),

 

 

 

 

 

-£ ^ = К е -'8= К к,

(6.126)

 

 

 

 

dz

 

 

 

 

т. е. производная комплексного потенциала является сопряжен­ ной комплексной скоростью.

Рассматривая произвольную аналитическую функцию / (г) = как комплексный потенциал течения, устанавливаем гидромеханический смысл ее производной: /'(г) является со­

пряженной комплексной скоростью.

Заметим теперь, что при безотрывном обтекании твердой стенки идеальной жидкостью вектор скорости V направлен по касательной к стенке и, следовательно, она является линией тока. Поэтому в течении с каким-либо комплексным потенциа­ лом молено любую линию тока заменить твердой стенкой и

комплексный потенциал будет описывать обтекание этой твер­ дой стенки (принцип отвердевания). Поясним на примере.

Рассмотрим комплексный потенциал

 

 

W= k ( z + - Z L j .

 

(6.127)

Отделяя в нем действительные и мнимые части, получаем

w = k ( x + ly + J L ) = k ( x + R * —

) + i k ( y - & J-

).

I

l X+iyl

\

*

Х2 +

у2 f 1

\

x 2 + y lj

Этому комплексному

потенциалу соответствует

картина

тече­

ния, определяемая линиями тока

 

 

 

 

 

 

 

 

 

 

 

 

< 6

Л 2 8 >

где С — произвольная постоянная. Полагая

С = 0, получаем

 

 

y(x2+ y 2—R2) = О,

 

 

 

т. е. у =

0 и

 

 

 

 

 

 

 

 

 

x * + y 2= R 2.

 

(6.129)

Таким образом, линией тока рассматриваемого течения явля­ ется, в частности, окружность (6.129) радиуса R с центром в нуле. Следовательно (принцип отвердевания), комплексный потенциал (6.127) определяет обтекание бесконечного (в на­ правлении оси Oz) круглого цилиндра радиуса R. При таком обтекании поле скоростей определяется сопряженной комплекс­ ной скоростью

 

^ Г

= А ( 1 _ - ^

) = К е " "

(6Л30)

откуда при 2 = о о получаем Vm>é~lüae =

kt

т. е. 1

 

 

 

Vm= k , 6 „ = 0 ,

 

 

(6.131)

где Voo,

0оо — соответственно модуль

и направление

скорости

на бесконечности.

потенциал

(6.127)

определяет

обтека­

Итак,

комплексный

ние бесконечного круглого цилиндра плоскопараллельным по­ током жидкости, скорость которого на бесконечности равна Voo—k и параллельна действительнойоси. Формулу (6.127)

можно переписать в виде

. (6.132)

Из формулы (6.128) следует, что линии тока, соответствующие отличным от нуля значениям постоянной. С, располагаются, как изображено на рис. 6.25.

Найдем теперь распределение скорости Кц на поверхности

цилиндра. Очевидно [см. (6.130)], Кц= dw

Полагая в

d z

 

формуле (6.130) z = R e ‘°, k=Voo, получаем

 

dw

d z

и, следовательно,

d w

V 'u - d z

откуда

е- ш )=К «(1 —cos2O-f-/sin20)

=V( 1 - cos 20)2+ sin 220,

l/4= 2 K eo|sin0|.

(6.133)

Как видим,

в точках 0

= 0 , 0 = ji

цилиндра

скорость

жидкости

равна нулю

(критические точки),

а в точках

0= JI/2, 0= Зя/2 ско­

рость максимальна и равна 2КТО (удвоенной скорости на беско­ нечности). При произвольно за­ данной форме обтекаемого кон­ тура непосредственный подбор

соответствующего ему комплексного потенциала затруднителен. В этом случае поступают (рис. 6.26) следующим образом. Кон­ формно отображают внешность рассматриваемого профиля (контура), 5 на внешность какого-либо канонического профиля

Si (например, на внешность окружности), для которого

комп­

лексный потенциал обтекания известен:

 

«(С)=?(МЖ Ф(М)-

(6.134)

При этом профиль S конформно отобразится

на профиль S b

Так как Si — линия трка, то на ней, т. е. при

(£, TJ) œ S I,

Ÿ(?»4)=Q =const.

(6.135)

Пусть аналитическая функция £ = £ (z), или,

подробнее, £+

+ t il= I(л, У)+Щ(х* У)* что эквивалентно двум равенствам

£=£(*,*/), -п—ЧU .0 ).

(6.136)

реализует конформное отображение внешности профиля 5 на внешность профиля Si, причем бесконечно удаленные точки плоскостей z и £ переходят друг в друга. Рассмотрим аналити­ ческую функцию w(z) = w[t)(z)]i которая получается из ку(£) заменой £ на £(z), т. е. заменой переменных £, TJ функциями (6.136) . Выясним, комплексным потенциалом какого течения на плоскости z является функция w (z) =w[t,(z)]. Учитывая вы­ ражения (6.134) и (6.136), имеем

w (z) = w (z)]=

<р[g(я, у ), YJ(х, у)]

[g(х, у ), т] (л, у)] =

 

 

 

 

=?(•*:, у )+ $ ( х ,у ) .

 

 

 

 

Если точки

(х,

t/)e S , то

в силу

конформного

соответствия

(6.136)

точки (|,

) œ S I

и,

следовательно [см.

(6.135)],

ф[£(х,

1/)»

ÿ )]= ^ i

при

(х,

у) œ S или

ф(х,

г/) = C j

при

(х,

y)e=S,

т. е. контур S является линией тока течения с комплексным по­

тенциалом

до[£(г)].

Как

видим,

комплексный

потенциал

o>[£(z)] определяет обтекание заданного профиля S на плоско­ сти Oz. Из изложенного следует, что для получения потенциа­ ла w[£(z)], т. е. для отыскания потенциала обтекания контура произвольной формы, нужно знать функцию £(z), осуществляю­ щую конформное отображение внешности этого контура на внешность какого-либо канонического контура (например, круглого цилиндра) и комплексный потенциал о>(£) обтекания канонического контура. Предположим, что обе эти функции из­ вестны. Тогда комплексная скорость течения в плоскости z определится формулой

 

 

 

(6.137)

d z

d C d z

d C \ d C J

'

Отсюда, используя обозначения

 

 

W

* .

 

(6.138)

получаем следующие соотношения между абсолютными вели­ чинами и направлениями скоростей в плоскостях z и £:

V,= -!л£-. е2= ес+ а .

(6.139)

Параметры потока на бесконечности в плоскости £ определя­ ются формулами, вытекающими из равенств (6.139)

Vi<x>=A<>VZOB, 8Соо—

— а»,

(6.140)

после чего потенциал до(£) [см. (6.132)] полностью известен.