книги / Методы математической физики и задачи гидроаэродинамики
..pdfЭта теорема указывает, в частности, что не существует ана литической функции, осуществляющей конформное отображе ние конечной области на кривую или точку.
Принцип сохранения границ. Если функция 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 но аргумент, отличающийся на 2л/л (точки А, В на рис. 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) |
точки (|, |
TÎ) œ 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)] полностью известен.