Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Прикладная гидрогазодинамика..pdf
Скачиваний:
87
Добавлен:
15.11.2022
Размер:
19.64 Mб
Скачать

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

Докажем, что циркуляция скорости по любому замкнутому кон­ туру, не охватывающему точечный вихрь, например, по элементар­ ному контуру 1, 2, 3, 4, равна нулю

d Г1,2,3,4 ==(Wa+ dWu) (г4- dr) cte - Wurdb=

=dbd (Wmr)=d6d (const)=0 .

Контур 1, 2, 3, 4 выбран произвольно, поэтому на основании перво­ го следствия теоремы Стокса заключаем, что вся область течении, за исключением точечного вихря, потенциальна. В этой области все жидкие частицы движутся поступательно по криволинейным траек­ ториям, деформируются, но не вращаются около собственных осей. Если мысленно провести на поверхности элемента линию, напри­ мер 13, то во время движения эта линия будет параллельна своему начальному положению, как стрелка компаса, вращаемого по окружности.

Определим потенциал

скорости

и функцию тока

с точностью

до постоянной, интегрируя уравнения (3.57) и (3.48):

 

ф= - ^ -

б = -^ — arctg — ;

Inг= ———In

у?. (3.59)

т 2я

лг т 2я

*

В этих формулах циркуляция скорости Г характеризует интенсив­ ность точечного вихря.

В реальных случаях потенциального вращения жидкости вмес­ то точечного вихря имеет место ядро вихря с конечным радиусом г0. В ядре вихря жидкость вращается по закону вращения твердо­ го тела © = H7u/r=const и максимальная скорость имеет конечное значение U^uo=a>r0 (пунктир на рис. 3.8).

3.9. СИНТЕЗИРОВАНИЕ БОЛЕЕ СЛОЖНЫХ ТЕЧЕНИЙ

ИЗ ПРОСТЕЙШИХ

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

алгебраическая сумма расходов источников и стоков должна равняться нулю;

одна из линий тока (нулевая) должна совпадать с контуром тела;

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

Пример 1. Диполь. На рис. 3.9 даны линии тока течения, полученного нало­ жением источника А и стока В одинаковой мощности Q. Вследствие смещения

источника и стока от начала координат, их потенциалы скорости и функции то­ ка в соответствии с формулами (3.56) примут вид

п = (Q/2л) In V ( X + е)2

= (С?/2я) arctg [у/(х + в)];

<РС = —(Q/2JI) In V(x е)2 -ьУ2; Фс = —(0/2я) arctg [*//(* — е)] •

Для результирующего течения ф= <ри+ фс и ф =,фи + фс. Если источниками сток сближать, сохраняя постоянной мощность, то при их совпадении (е *0 ) тече­ ние прекратится — сток поглотит источник.

Рис. 3.9. Взаимодействующие источник А, сток В и ди­ поль

• Диполем (см. рис. 3 .9 ) называется течение, возникающее при одновремен­ ном стягивании источника и стока в начало координат (е— »-0 ) и увеличении их

мощности (Q-*оо), но так, что момент диполя iVf=lim2 Q. е сохраняет посто-

е->0

Q—►оо

янное значение:

мIn V (X + е)2 + #2 — 1п 1Л* — £)2 + У2

«Рдип = Пт (<fu+ <РС) = о7 Iim

 

 

 

8—>-0

^

2 Л

е->0

 

 

 

Q—Vоо

 

 

 

 

 

 

 

 

Фдии = Иш (фв 4- Фс) =

м

lim ■

arctg х

4- е

— arctg JC— е

 

 

 

6-+-0

 

и

2 л:

е-*0

 

 

 

Рассматривая в этих уравнениях 2е как

приращение

аргумента,

а числитель

как приращение соответствующих функций, получим

 

 

М

д

,

______

 

М

д (

У \

* » - s r

571(|п ^

+ ^

-

аГ

37 Г * 8 т ) ■

Выполняя дифференцирование, получим потенциал скорости и функцию тока ре­ зультирующего течения — диполя

 

 

М

 

х

М

cos 6

#

 

 

 

 

<frtHU

х°- + y i

г

 

 

 

(з. 60)

 

Фдип =

М

и

 

уИ

sin О

 

 

 

х^ А- У2

 

2л;

 

г

 

 

 

 

 

 

 

 

где 0 — угол между радиусом-вектором точки и осью х.

Семейство . ^

пляшмр

ф = С, х2 + у2 = Су

представляют

окружности

с центрами

Ha2^

2_ ^ ^ P _

окруЖ.

через начало координат; эквипотенциальные линии

Ф =

с >

х

т}РПр Пгт п ш »

«ости с центрами

на оси

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

начало

координат, вследствие

равенства^ мощности источника и стока результирующий поток жидкости через> замкнутый контур, охватывающий диполь, равен нулю. Это свойство обеспечи­

вает широкое применение диполя для синтезирования еще более сложных тече­ ний жидкости около твердых тел.

Задача 3.22. Докажите, что фДИц и фдИП (3.60) удовлетворяют уравнению

Лапласа.

Пример 2. Поперечное обтекание бесконечно длинного кругового цилиндра, радиус которого г0. Рассмотрим это течение как плоское и покажем, что она

Рис. 3.10. Поперечное обтекание ци-

Рис. 3.11. Циркуляционное обтекание

линдра идеальной жидкостью

иилинппя

(^i; Ф,; ф,) и диполя

Л*—о

Р! $ ЛЬтат наложения плоскопараллельного поток»

функцию тока ф синтрэип,,о,,„

0

W{' ^2’ Фг. Определим потенциал скорости ® и

и значения М;

Руемого

течения с использованием формул (3.54), (3.60),

* = Ь + 42 = Г , + rlx/(x2 + у2)] r= Wx (г + гЦг) cos в,

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

<3' 61>

Нулевая л н н и я^ка^реясгя .л ^

3

Ч’= С

нУлевУю

У'= 0 : У11—л«2/(*2+г/2)]= 0 .

изводящей контуп O6 TP**AJ

Яет

совокУпность

окружности радиуса г0, воспро-

(рис. 3.10),

т. е улойлртппп™ цилинДРа, и прямой у=0,

совпадающей с осью х

окружную Wu скоппрти „®?ряет граничным условиям задачи. Радиальную

1РГ н

 

 

 

 

уисги найдем по

определению

 

 

 

 

 

 

 

 

 

 

 

 

 

^

= df/dr = r j

(1 гУг2) cos 0;

 

 

 

 

г

>()

 

 

Г “ =

^ /<?/ =

(ft/rdb) =

- Г

х ( 1 +

'о/'"2)

sin 0;

 

 

 

положительно^уРняВппя»поВ„ СТорону

Увеличения л

и

JFu> 0, если соответствует

ном удалении от цилий™.

враи*ения

(против часовой стрелки). На бесконеч-

оси y(Q = n/2 и

Д

ндра

*-оо) как вдоль оси х (0 = 0

и 0 = Jt), так и вдоль

б

2

тг|

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л}и при любом 0 имеет место плоскопараллельное те-

бесконечност/^тНтеля

г^Т0Ка

со скоРостью

 

 

 

Поэтому

«скорость

на

синонимами

Ня

 

 

то>>

и <<ск°РОсть

невозмущенного

потока

являются

 

 

 

па стенке Дилиндра при г=г0; И7Г~ 0

и

 

 

 

 

 

 

Формул

36

 

 

Г =

= -

Sin в.

 

 

 

 

(3. 62)

В

точках ^

и2 В °ппи bfl—Т

РаспРеделение

скоростей

по поверхности цилиндра.

передней критической

 

 

о

0

скоРость равна

нулю.

Точка А называется

торой потоки вновь спрпни

В этои точке поток раздваивается. Точка В,

в ко-

 

Задача

3 23

Лл

единяются» называется задней критической

точкой.

3.10)

требуется-

 

* ’

МЛЯ

случая обтекания бесконечного

цилиндра

(см. рис.

1.

Доказать,

что при г=г0 и углах 0= 90 и 270°

окружная

скорость Wu =

— =F2WU при 0=30°

Wu= — Wi и радиальная Wr= 0.

 

 

2.

Объяснить причину увеличения скорости при 0=90 и 270°

 

3.

Доказать,

что

уравнение линии тока, проходящей

через точки х = 0 , у =

= 2г0

имеет вид

у—r02yl(x2+yz) =Зг0/2, а расход жидкости между

этой линией

тока и нулевой линией тока Q= 3WVo/2.

4 . Указать положения критических точек А и Б в случае, когда поток будет

•направлен под углом 30° к оси х.

Пример 3. Циркуляционное обтекание цилиндра (рис. 3.11),. Циркуляционное обтекание бесконечного цилиндра получим наложением полей соответствующего

бесциркуляционного обтекания цилиндра (q>i,

и присоединенного потенциаль­

ного вихря

с циркуляцией Г. В

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

формулами (3.61) и (3.59), по-

.лучим

 

 

 

 

 

 

 

 

 

‘Р = 91+Ч>2=^Гоо (г + гЦг) со& b +

|

(о. Ои)

 

 

Ф= Ф1 + Ф2=

 

 

/

 

 

— Гр/г) sin 0 4- Г In г/2л. J

 

‘Уравнения эквипотенциальных линий и линий

тока получим,

положив в

(3.63)

ф= const

и

ф = const.

 

определяют

радиальную

и ок­

Частные

производные потенциала скорости

ружную

составляющие скорости

 

 

 

 

 

W г =

ду/дг = W <a(l — гЦг2) cos 0;

 

(3. 64)

 

 

W a =

dy/rdf) ——

(l + rl/r?)

sin 0 + Г/2лг.

 

 

 

 

 

На поверхности цилиндра r— r0, Wr= 0 и распределение скорости имеет вид:

 

 

W =

W u = —2W oo sin 0 +

Г/2яг0.

(3. 65)

Задача 3.24. Для условий рис.

3.11 и циркуляции

Г*——2nroWoo доказать,

что: 1) скорость жидкости на поверхности цилиндра (г=/*о)

при различных

уг­

лах 0 имеет следующие значения:

 

 

 

 

 

 

 

 

 

0

90

 

180

210

270

330

360

 

Wu

—W„ - 3W m

 

 

0

+ W „ 0

- W "

 

'2) вектор скорости

в

точке

г= 2 г0

и

0= 30°

равен W «1,3

и составляет с

осью х угол —‘30°.

 

 

 

 

 

 

 

 

 

 

 

Вопрос 3.25. При

каком

значении

циркуляции

скорости

Г, критические

точ­

ки Л и В совпадут

в точках х = 0 , у = —г0? Ответ: Г= —4яг0№оо.

 

3.10. О МЕТОДЕ КОНФОРМНЫХ ОТОБРАЖЕНИЙ

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

Ф о р м у л и р о в к а з а д а ч и .

Пусть в физической плоскости z комплекс­

ного переменного z = x + iy задан

произвольный крыловой профиль I, обтекае­

мый плоским потенциальным потоком идеальной несжимаемой жидкости, ско­ рость которого на бесконечности Woo составляет угол 0 с осью х (рис. 3.12).

Требуется определить поле скоростей во всей области течения, внешней по от­ ношению к контуру /.

Пусть уже решена сложная, но чисто геометрическая задача конформного отображения —найдена однозначная аналитическая функция комплексного переменного:

С= F{z) или z = f { С),

(3.66)

конформно отображающая область внешнюю, относительно иссле­ дуемого контура I на плоскости комплексного переменного z= x + Л-iy на область, внешнюю относительно круга /*,. во вспомогатель­

ной плоскости комплексного переменного £=|-Ить а также осу­ ществляющая обратное отображение. Для того, чтобы конформное отображение было единственным, функции (3.66), в соответствии с теоремой Римана, выбираются так, чтобы точкам z = oo в физи­ ческой плоскости соответствовали точки £=оо во вспомогательной плоскости и чтобы в этих точках производная dz/d£ была положи­ тельна, т. е.

(dz/d(.)oo = nioo > 0.

(3.67)

В этом случае, каждой точке z = x + iy крылового профиля I в плос­ кости z будет соответствовать одна определенная точка £ = £ + гп окружности /* во вспомогательной плоскости и наоборот. При кон­ формном отображении контуры изменяются так, что их бесконечно малые элементы остаются геометрически подобными и углы между касательными в точке пересечения двух кривых не изменяются и в данном случаеравны я. Исключение составляет особая точка В, для которой конформность отображения нарушается и острый угол отображается в угол равный я (точка В *).

Метод конформных отображений используется для решения гидродинамических задач потому, что вместе с контуром тела, при использовании той же самой функции (3.66), отображается и поле скоростей течения около него на поле скоростей циркуляционного

обтекания цилиндра во вспомогательной плоскости £ ,и

наоборот.

П р и м е н е н и е фу н к ц и й к о м п л е к с н о г о п е р е м е н н о ­

го при к о н ф о р м н о м о т о б р а ж е н и и течений.

Форму­

лы (3.50) в теории функций комплексного переменного называют­ ся условиями Коши — Римана, которые необходимы и достаточны,

чтобы выражения фН- г-ф и ф*-Н'ф* являлись аналитическими функ­ циями x()z и х* (?) комплексных переменных z= x + iy и £=|-И'т1,

Х (* )= Х (* 4 -^ )= ? + % Х*(5)=Х*(5+«Т1):=?* + ^*, (3.68)

где ф, ф —потенциал скорости и функция тока течения около про­ филя в плоскости z, соответственно; ф*, ф* — известные потенциал скорости и функция тока циркуляционного обтекания цилиндра (3.63), соответственно; %(z) — комплексный потенциал течения около профиля в физической плоскости z.

Комплексный потенциал циркуляционного обтекания кругового цилиндра во вспомогательной плоскости х*(?) получим, подставив

в (3.68) значения ф* и ф* из

(3.63) и производя преобразования:

 

Х*(С)= ^ : с +

Г :^ + - £ 1 - 1 п С ,

(3.69)

где W„ —сопряженная скорость на бесконечности

в плоскости £;

Г* — циркуляция скорости вокруг контура /*.

 

Комплексный

потенциал

течения около исследуемого профиля

X ( z ) неизвестен

и задача

сводится к его отысканию с последую­

щим определением ф и ф.

Соотношение между комплексными потенциалами течений при конформном отображении установим, используя отображающую функцию (3.66):

х(*)= х[/(0]= х*Й )= т+ *'!»= ?*+ /ф * -

(3.70)

Символ комплексного потенциала % сам по себе не подразуме­ вает какой-либо определенной функции, а лишь указывает, что он относится к определенному течению, рассматриваемому в данной

комплексной плоскости.

точках течений z и £ по­

Равенство ф=ф* в соответственных

казывает, что на обоих контурах ф и

имеют одно и то же по­

стоянное значение, т. е. контуры I и ■/* являются нулевыми линия­ ми тока.

С о о т н о ш е н и е м е ж д у с к о р о с т я м и н е в о з м у щ е н - н ы х п о т о к о в при конформном отображении найдем, продиф­ ференцировав (3.70) по £, учтя (3.67), и что сопряженная ско­ рость на бесконечности равна производной от комплексного потен­ циала по комплексной переменной:

Wl> = [dx*(l)/dQ'0 = [dx(z)/dZ]'o = [dx(z)/dz

dzld^\oo = W ооТПоо.

 

 

(3.71)

При отображении направление комплексных

и сопряженных ско­

ростей невозмущенных потоков не изменяется

(/Поо>0 ), а их вели­

чины изменяются в т » раз.

 

 

Циркуляция скорости вокруг контуров I и I* при конформном

отображении сохраняет неизменное значение.

Действительно, ис­

пользуя (3.71), (3.66), (3.67) и

(3.70) и учтя, что_циркуляция рав­

на действительной части (д. ч.)

(j) W dz и д. ч. | W*d%, получим

Г = д .

d z

 

3

d z

^ . л _

j

 

dtl

=

Д.

 

<*x* (0 £/C =

r*

(3.72)

 

 

/*

diC

 

 

 

 

 

 

 

Итак, при конформном отображении сохраняются неизменными на­

правление невозмущенного потока и циркуляция скорости и в rriooidzldt^oo раз изменяется абсолютная величина скорости.

П о р я д о к о п р е д е л е н и я к о м п л е к с н о г о п о т е н ц и ­ а л а о б т е к а н и я з а д а н н о г о п р о ф и л я :

1) заданный профиль конформно отображается на вспомога­ тельную плоскость £:

а) определяется радиус окружности г0 и расположение окруж­ ности на плоскости £• Для этого используются геометрическое опи­

сание профиля в плоскости z = x + iy и функция (3.66);

скорости

б) определяется масштаб конформного отображения

в соответствии с (3.66) и (3.71) W*00IW 00= {dzld(,)00= m0o

и учи­

тывается (3.72). Все эти данные подставляются в (3.69) и находит­ ся явный вид комплексного потенциала циркуляционного обтека­ ния круга

Х*(С)= /га-^Д + от.Г.г?/С + (Г/2щ)/1п<;;

(3.73)

2) определяется искомый комплексный

потенциал течения око­

ло профиля I в физической плоскости z.

Для этого учитывается,,

что x(z)=X*(t)> и в (3.73) подставляется функция £,=F(z):

 

X{z)= m„W»F (z) + m „ W y 0/F (z) +

(Г/2яг) InF (z).

(3.74)

Определение ф и ф сводится к выделению действительной и мни­ мой частей (3.74).

Как видим, гидродинамическая задача отыскания комплексного потенциала обтекания заданного профиля заданным потоком не представляет труда, если известны отображающая функция (3.66) и циркуляция скорости Г.

Функция z —0,5 (£“Ьг02/£), отображающая круг на профиль кры­ ла, была найдена Н. Е. Жуковским в 1910 г. и названа его именем. Применяя эту функцию Н. Е. Жуковский и С. А. Чаплыгин получи­ ли серию теоретических крыловых профилей. Профили, отличаю­ щиеся от теоретических, при отображении дают искаженный круг и метод конформного отображения применим лишь для прибли­ женного исследования их обтекания. Циркуляция Г для круга мо­ жет иметь произвольное значение и поэтому должна быть задана такой, какая действительно возникает при обтекании профиля I.

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