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

книги из ГПНТБ / Бекнев В.С. Газовая динамика газотурбинных и комбинированных установок учеб. пособие

.pdf
Скачиваний:
49
Добавлен:
25.10.2023
Размер:
14.39 Mб
Скачать

§ 7. ПОТЕНЦИАЛЬНОЕ ДВИЖЕНИЕ

ИЕГО СВОЙСТВА

Выше было введено понятие о потенциальном движении и

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

ср.

Если известен ср, то путем дифференцирования можно найти

про­

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

Потенциал скорости есть кинематическая функция, характери­ зующая собой потенциальное движение жидкости. Каждому потенциальному потоку соответствует собственный потенциал скорости ф (.V, у, z, і), причем жидкость может быть как несжи­ маемой, так и сжимаемой.

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

факта, что потенциал скорости ср для

несжимаемой

жидкости

является

решением

уравнения Лапласа.

 

 

 

 

В самом деле, для несжимаемой жидкости уравнение нераз­

рывности

имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

дѵх

I

дѵц_

I

_dtfe

_

ç.

 

 

 

 

дх

~ і

ду

 

dz

~

 

'

 

 

HO

 

 

 

 

 

 

 

 

 

 

 

 

Èk.

JL

( Ü!L\

 

г

П

 

 

dx

~

dx

\ dx J

~

dx°-

11

д '

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

$ - ' - ^ + ^ = ^ = » -

 

Полученное уравнение называется уравнением Лапласа, а

функция

ф, удовлетворяющая

этому

 

уравнению,

гармони­

ческой функцией.

Докажем, что потенциал скорости ф для несжимаемой жидкости не может иметь экстремума внутри жидкости. Для этого покажем, что проекция скорости на любое направление равна производной от Ф по этому направлению, т. е. ѵп = дуідп.

Возьмем в жидкости произвольное направление я и вектор скорости V и составим скалярное произведение

va = ѵп = ѵх cos (il, x) -f- vy cos (л, y) - j - vz cos (п., z).

Затем продифференцируем функцию ф по п, тогда

_оф_

дер

дх

.

ду

I

_оф_

_ô_£ _

дп

дх

дп

ду

дп

 

дг

дп

= vx

cos (п,

х) - j -

Vy cos

(n,

y)

- j - vz

cos (n, z).

40

Сравнивая полученные выражения, будем иметь

Воспользуемся формулой Остроградского — Гаусса и заменим в ней ѵп = dcp/d/7, учитывая, что div ѵ = 0, тогда

J vnds = J - ^ - ds = J divörfV = О,

s

s

1'

т. е. для несжимаемой жидкости

s

Итак, предположим, что в некоторой точке M значение ср имеет максимум. Окружим эту точку достаточно малой сферой так, чтобы вся сфера была внутри жидкости. Тогда для каждой точки сферы будем иметь

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

s

что противоречит условию (26). Аналогичные рассуждения можно провести и для случая минимума.

Следовательно, функция ср не может достигать экстремума внутри жидкости. Однако она может получить наибольшее или

наименьшее значение

на

границе

области, занятой

жидкостью,

но это значение не

обязательно

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

аналитическому

экстремуму.

 

 

 

 

 

 

 

 

Рассмотрим свойства

потенциального движения

несжимаемой

жидкости.

 

 

 

 

 

 

 

 

1. Величина скорости

внутри жидкости не может иметь макси­

мума.

 

 

 

 

 

 

скорости ѵх, ѵу

Предварительно

покажем,

что

компоненты

и ѵг при потенциальном течении удовлетворяют

уравнению Лап­

ласа, т. е. являются гармоническими

функциями.

 

В самом деле, продифференцировав уравнение Лапласа по х,

получим

 

 

 

 

 

 

 

 

£ < А » > - А ( - £ ) = А . , = 0.

 

 

Поступая аналогично,

будем

иметь

 

 

 

 

Аѵу

= 0;

 

àvz

= 0.

 

 

Следовательно, компоненты скорости при потенциальном тече­ нии не могут иметь экстремума внутри жидкости.

41

Посмотрим, может ли скорость иметь экстремум внутри потока жидкости.

Пусть в некоторой точке А скорость имеет максимальное зна­ чение. Помещая начало осей координат в точке А и направив

ось X по вектору ѵА, получим, что ѵхА = ѵА = i>r a a x .

Однако в силу доказанного ранее в жидкости найдется точка В, где ѵхВ > ѵхА. Следовательно,

ѵв = У ѵхВ + v2yB + ѵ\в

^ ѵхВ > ѵхА = ѵА,

 

 

т. е. при потенциальном течении

скорость не может

иметь

макси­

мума внутри потока жидкости.

 

 

 

Относительно минимума таких рассуждений провести

нельзя.

В самом деле, допустив минимум скорости в точке А , получим

для соседней точки В ѵхВ •< ѵхА

= ѵА.

 

 

Но вследствие существования

компонентов ѵуВ

и ѵгВ

можем

получить ѵв > ѵА, т. е. минимум возможен.

Поскольку минимум величины скорости это нуль, то, следова­

тельно, внутри

потока скорость

может иметь нулевое

значение.

2. Свойство

единственности

потенциала скорости.

Докажем,

что значение функции ф внутри какой-то замкнутой области опре­

деляется:

1) единственным

образом,

если на границах

области

заданы значения

функции

ср, 2) с точностью до произвольной по­

стоянной,

если

на границах области

заданы значения

где

п — внешняя нормаль

к контуру.

 

 

Для доказательства

этого свойства

воспользуемся интегралом

 

J [ ( * ) ' + ( f ) " + ( * ) > = J * * * .

m

где

V

s

 

фгармоническая

функция.

 

Рассмотрим первую

часть свойства.

 

Пусть на контуре

поверхности s задано значение

функции

Ф =

ф5 . Предположим, что внутри области имеются два решения q>1

и ф2 , но на поверхности они совпадают, т. е. ф 1 5 = cp2s. Рассмотрим новую функцию ср3 = ср1 — ф2 , причем на поверх­

ности s значение ф 3 з = 0, и применим к ней равенство (28). Правая часть равенства равна нулю, так как ф3 5 = 0, а значит, равна нулю и левая часть. Но подынтегральное выражение представляет собой сумму квадратов, следовательно, каждое слагаемое должно

быть

равно

нулю:

 

 

 

 

 

дфз

дфз _

дфз _ _

g

 

 

дх

ду

dz

'

т. е. фз = const по всей области.

 

Но на поверхности ср3

= 0, следовательно, const = 0 и ф х =

= ф 2

по всей

области течения жидкости.

доказана.

Таким образом, первая

часть

свойства

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

 

 

Фз = Фі — Фа-

 

 

На этот раз на границе задано значение производной

Тогда

дфі

дф2

получим, что на поверхности

 

 

 

Un

 

 

 

 

 

 

 

 

 

 

 

 

 

дфз

_

Q

 

 

 

 

 

 

дп

s

 

 

 

 

Правая

часть равенства

(28) опять обращается в нуль, но уже

вследствие того, что <5сра/<3/г = 0. Тогда

будет равна

нулю и левая

часть:

 

 

 

 

 

q

 

 

 

 

дфз _

дфз _

дфз

 

 

 

 

дх

ду

 

dz

'

 

 

т. е. ф 3 = const

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

ц>1

= Ф2 + const.

 

Таким образом, доказана и вторая

часть свойства.

 

3. Кинетическая энергия массы несжимаемой

жидкости при

потенциальном

движении

минимальна. Докажем,

что кинетиче­

ская энергия Т массы жидкости при потенциальном течении меньше кинетической энергии Т' при любом непотенциальном течении, т. е. Г' — Т >> 0 при условии, что нормальные составляющие ско­

рости к поверхности выделенного объема в обоих случаях

равны:

^ - = vn = v'n = v'xcos[n, X) - j - v'y cos(ii, y) - j - v'zcos(n,

z).

Кинетическая энергия массы при потенциальном потоке

г - J ' 4 " - - И [ ( £ ) * + ( £ ) ' + ( £ ) > .

 

V

V

 

Для

вихревого

потока

 

 

Г

= ±\(ѵ? + ѵ* +

ѵ?)аѴ,

 

 

V

 

откуда

их разность

запишется в виде

)'+*-(-£)']dV-

г-т = цp-(^)4^-(f

 

 

 

(29)

Как известно, уравнение неразрывности для несжимаемой жидкости при любом течении имеет вид div ѵ = 0, т. е. для вихре­ вого потока

дѵ'

дѵ'„

дѵ,

—— -\

£- -1

= 0

дх '

ду *

dz

43

а для потенциального потока

 

 

 

 

Дер = О

Сделаем

 

следующие преобразования

-

Щ

2

V

х

 

дх

)

Откуда

Аналогично

V-(f)! 4°-f)2 + 2»»f-2(f)';

Оф

~~дг

Уравнение (29) с учетом преобразований можно записать в виде:

а (г- _ г ) = j [(„;, _

 

+ ( » і -

f

у + (»; -

tvi

dV

оф

,

'

Оф ,

'

оф

dV~

 

 

I F

~ г

y

l)y~*

 

2~дТ

 

 

 

 

 

 

 

 

 

 

 

dz

dV.

 

(30)

 

 

 

 

 

 

 

 

Преобразуем в правой части последние два интеграла. Послед­ ний член в равенстве (30) можно преобразовать так:

! I[(ï)'+(îJ+(î)V-! I'îf c

Далее, так как

 

 

 

 

UV

то

£_

 

 

 

с»Ф

( 0 я ф )

ср

дх

öx

дх

 

 

Аналогично

 

 

 

 

Оф

 

 

 

 

~~

ду

 

 

 

дф

а

/ • \

 

дѵ

ôz

öz

 

 

 

44

Второй член правой части равенства (30) можно написать в виде

дц>

 

 

дх

Ѵ у ду + V z

dz

 

2 1 [ - е - ( ^ ) + і - К ф ) + - е - ^ ф )

dV —

 

 

 

- 2

 

dvr

 

 

дѵп

дѵ^

 

dV.

(31

 

 

 

— - Ч

 

-

+ • dz

 

 

 

 

 

дх

'

 

ду

>

 

 

 

 

Но так

как

по формуле

Остроградского — Гаусса

 

 

 

 

 

 

 

 

 

 

 

 

 

dV =

 

 

v

 

 

 

 

 

 

 

 

 

 

 

 

 

=

j div (ц>ѵ ) dV =

J ери,, ds =

j

ф

ds

 

и, кроме того,

div

v

0,

то

равенство (31)

перепишется в

виде

 

J L

 

 

 

 

 

ду

dV

 

J<p - |r d s -

 

2

~дх

~

J ~dy

^

 

z

dz

 

 

 

 

 

 

После

всех

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

равенство

(31)

будет иметь

вид

т-т-Ш''-%)'Н"'-Ъ)'Н*-Ъ)'

dip \2" dV

 

 

ил и

 

 

 

( Г -

 

Т)

>

0,

 

 

 

 

 

 

 

 

 

 

 

 

 

так как подынтегральное выражение представляет сооои сумму квадратов, т. е. всегда положительная величина.

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

Третьему свойству можно дать наглядное механическое толко­ вание.

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

ческая энергия диска Тл = т-^-, где ѵс — скорость центра диска.

Вихревое движение (диск при движениии вращается вокруг своего центра) будет соответствовать качению диска с трением вдоль той же направляющей. Кинетическая энергия диска в этом слу-

чае Гд = m + I-пг, т. е. будет больше на величину / - О - .

Следовательно, для движения" жидкости с потенциалом ско­ рости надо затратить меньше энергии, чем для движения с теми же скоростями, но без потенциала скорости, т. е. для вихревого дви­ жения. Вихри всегда аккумулируют в себе часть энергии потока, поэтому течение жидкости в турбомашине целесообразно по воз­ можности приближать к потенциальному.

§ 8. ПОТЕНЦИАЛЬНОЕ ДВИЖЕНИЕ ЖИДКОСТИ ПРИ НАЛИЧИИ ДИСКРЕТНЫХ ВИХРЕВЫХ ШНУРОВ. ФОРМУЛА БИО—САВАРА

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

Понятие о циркуляции скорости. Теорема Стокса

В векторном анализе имеется понятие о циркуляции Г век­ тора а по замкнутому контуру с, которое формально напоминает понятие работы некоторой силы на пути dl

с

Применяя это понятие к полю скоростей, получим

с

или, раскрывая скалярное произведение векторов ѵ и dl, можем записать

с

Для потенциального потока

 

ду .

ду_

~~ ду

~ dz '

тогда

 

dtp

 

дх

с

с

46

В случае однозначного потенциала скорости <р, когда интеграл по замкнутому контуру с от ср равен нулю, циркуляция по этому контуру тоже равна нулю.

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

У

•s.

т.

 

 

^

v-r=const

 

\

X

п, J

Рис.24. Поле скоро­ стей в окрестности вихревого шнура

ѵг = const (рис. 24). При этом скорость жидкости постоянна на окружности радиуса г и циркуляция скорости по контуру с

Г = (j) V aï = (j) vrda — 2nrv = const =j= 0.

с0

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

величину 2пгѵ.

Этот результат сохраняется до тех пор, пока vr

~

= const. При

выборе контура с внутри вихревого шнура, где,

как в

твердом

теле, ѵ = cor, циркуляция будет зависеть от

ра­

диуса

г.

 

 

Соотношение между циркуляцией скорости и завихренностью внутри контура с устанавливается теоремой Стокса, которая гла­

сит, что циркуляция

скорости

по

замкнутому

контуру

равна

сумме удвоенных интенсивностей

вихрей, охваченных этим

контуром.

Рассмотрим элементарный контур в плоскости хоу (рис. 25).

Пусть в точке А скорость имеет проекции

ѵх и ѵу.

Составим выра­

жение для циркуляции по контуру

ABCD

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

против

часовой стрелки. Так

как стороны

контура равны dx

и

dy, то,

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

для AB

ѵ2 ад: '

47

для

ВС

д ля

CD

для

DA

Составляя

интеграл

по контуру

ABCD,

получим

 

 

 

 

 

 

 

 

 

 

 

А ВСDA

 

Л

Г )

 

с

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'

 

• дѵх

,

I 1 дѵх

, \ ,

А

ал

 

8

X

 

^ . v - r - ^ d y - r ^ - - - ^ d x ) d x -

 

 

 

 

 

 

 

 

Рис. 25.

Элементарный

контур

 

 

 

 

 

 

 

A BCD в

поле

течения

жидкости

 

 

 

 

 

 

 

Раскрывая

скобки

и приводя

 

подобные

члены, будем иметь

 

 

 

« г - ( - & - $ - ) * *

 

 

 

или, заменяя

! ,

j

 

и

г =

'

/ du,.

ou. \

П 0 Л У Ч И Ы

dxdy =

as и вводя

 

-тт^-^г

Л / / '

 

 

 

 

аТ =

2azds

 

 

 

 

 

Таким образом доказана теорема Стокса для элементарного прямоугольного контура. То же самое можно получить и для треугольного контура, например, ABD. Интеграл по BD следует вычислять как от суммы произведений проекций и и dl на осп координат.

Обобщение для произвольного контура проводится следующим образом. Разделим контур с на ряд элементарных контуров пря­ мыми, параллельными осям координат. Площадь s разделится на элементарные прямоугольники и треугольники, гипотенузы ко­ торых сколь угодно близко подходят к отрезкам контура с.

При суммировании циркуляции по элементарным контурам получим циркуляцию по внешнему контуру с, так как составляю­ щие циркуляции по внутренним контурам взаимно уничтожаются.

48

Следовательно, для произвольного контура в плоскости хОу

 

 

Г =

2 [ ю г

в Ь .

(32)

 

 

 

s"

 

 

Если контур

с

охватывает

я вихрей с интенсивностями

;.S,-,

то циркуляция

по

контуру с

 

 

 

 

 

 

п

 

 

 

 

Г =

2 £

tû,S j .

 

Из этой теоремы следует важный

вывод о том, что циркуляция

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

вращается

как твердое

тело.

 

Таким

образом, можно найти положение и размеры

вихревого

шнура в жидкости путем измерения скоростей в точках

выбранных

контуров и подсчета циркуляции по этим контурам.

 

Определение поля скоростей

 

в области

вихревой нити.

 

 

Формула

Био—Савара

 

 

 

В разлѵічных технических

задачах приходится определять поле скоростей

в области вихревых нитей.

К

этому приводит, например, задача

определения

так называемого индуктивного сопротивления крыла конечного размаха, впервые решенная С. А. Чаплыгиным в 1913 г.

Пусть имеется

бесконечно тонкая

вихревая нить с радиусом поперечного

сечения е; отрезок ab этой нити представлен на рис. 26.

Циркуляция,

взятая по

любому контуру, охватывающему нить, равна Г

и не зависит от s.

 

 

 

Скорость на поверхности

вихревой

нити, нормальная к оси нити,

 

 

 

Г

 

 

~~

2т '

4 B . C . Бекнев

49

Соседние файлы в папке книги из ГПНТБ