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

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

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

Из-под знака производной -щ- выйдет масштабный множитель

-^г-,

а из-под знака

операции rot

{v X В) множитель

.

' о

 

 

 

 

' о

Тогда

будем иметь

в

безразмерном

виде следующее уравнение:

 

^ • ^ - = ^ r o t ( w x î ) + ^ # - А Ѣ

 

 

' о

dt

' о

/ 5

 

или

= r o t ( ' 0 x ß ) + ^ - A 5 ,

где Rem =

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

 

-ПГ-^

= rot (ихсо) -4- p^-Acû.

 

 

' о « о

dt

 

R e o

 

Числовое соотношение между членами, ответственными за

процессы конвекции

и

диффузии,

определяется масштабными

множителями, стоящими

перед ними: единицей и дробью

или

Если процесс конвекции преобладает над процессом диф­

фузии,

то

 

 

 

 

 

1

»

Re7' Т - е -

^ » 1 "

 

Если же процесс диффузии преобладает над процессом конвек­ ции, то -p-î— ^> 1, т. е. Rem <^ 1. В данном случае характерный безразмерный комплекс

Rem = Т~ = ѵ1^та^

называется магнитным числом Рейнольдса (по аналогии с гидро­ динамическим числом Рейнольдса Re = vl/v).

Рассмотрим отдельно процессы конвекции и диффузии. Для

определенности будем рассматривать вектор В.

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

ПО

В рассматриваемом случае (когда число Rem велико) говорят об эффекте «вмороженности» магнитного поля в среду, поэтому

размерное уравнение индукции упрощается

и имеет вид

f - = rot(äx5) .

(97)

Для конвекции гидродинамического вихря со необходимо, чтобы было велико гидродинамическое число Re, это будет, в част­ ности, при вязкости, равной нулю.

Рассмотрим на примере вектора В магнитной индукции тео­ рему, имеющую место при больших числах Rem .

Рис. 50. Эффект «вмороженности» магнитного поля в проводящую среду

При движении жидкости с бесконечно большим магнитным числом Рейнольдса поток вектора В магнитной индукции через любую «жидкую» поверхность s не изменяется с течением времени.

В частности, это означает следующее: пусть в начальный мо­ мент времени tx взята какая-либо «жидкая» поверхность s, це­ ликом состоящая из магнитных силовых линий вектора индук­ ции В (рис. 50). Поток J Вп ds в момент времени tx равен нулю

s

в силу того, что в каждой точке поверхности векторы В и п пер­ пендикулярны. Согласно сформулированной теореме этот поток

сохранит

постоянную величину,

если число Rem велико; т. е.

в момент

t2 он будет также равен

нулю, а это возможно, если век­

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

Если условие Re,„ > сю не соблюдается, то жидкость при дви­ жении проскальзывает через магнитные силовые линии, или же можно сказать, что магнитное поле диффундирует внутрь объема вещества, выделенного для рассмотрения.

Подобная теорема для гидродинамического вихря со скорости известна как вторая теорема Гельмгольца о вихрях.

111

В случае, если жидкость идеальна, массовые силы имеют по­

тенциал

и жидкость

баротропна

[р = / (/?)], то вихревая

трубка

состоит в течение движения из одних и тех же частиц.

 

В этой формулировке требование идеальности жидкости экви­

валентно

условию, что число Re очень велико,

а два других тре­

бования

следуют как ограничения

при доказательстве формулы

(95)

Гельмгольца.

Утверждение,

что вихревая

трубка

состоит

во время

движения

из одних и тех же частиц,

эквивалентно ра­

зобранному

эффекту «вмороженности».

Обычно

применительно

 

 

 

 

к вихревым

трубкам

говорят об эф­

 

 

 

 

фекте их сохраняемости,

имея в виду,

 

 

 

 

что

при сделанных

ограничениях

 

 

 

 

вихри перемещаются вместе с жид­

 

 

 

 

костью.

 

 

 

 

 

 

 

 

 

 

Можно

также

показать, что при

 

 

 

 

тех же предположениях

при движе­

 

 

 

 

нии

сохраняются

не только

вихре­

 

 

 

 

вые

трубки

(векторов В и со), но и

 

 

 

 

их интенсивности.

 

 

 

 

 

 

 

 

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

теормы о «вморо­

 

 

 

 

женности» вектора

В начнем с вычи­

 

 

 

 

сления потока J Вп ds через

поверх-

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

ность цилиндра (рис. 51), образован­

 

 

 

 

ного движением его нижнего осно-

Рие. 51. Поток

вектора

магнит-

с о

с к о р

о с т ь ю

„ среды. В силу

нон индукции через поверх-

 

^

 

'

 

>-

J

 

ность

цилиндра

теоремы

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

 

 

 

 

J" Вп ds = j

div В dV = 0 этот

поток

 

 

 

 

s _

 

v

 

 

 

 

 

равен

нулю, так

как div 5 =

0,

где V — объем

рассматри­

ваемого цилиндра,

ограниченного

поверхностью s=sH -|-s6oK -f-sB .

Для момента времени

t2, когда

 

с момента времени tx нижнее

основание пройдет

путь

v At, этот

поток

— { В (t2) nads + J Б {Qna

dt+ J В (t2) пбок ds =0.

S 1I

5 B

 

s6OK

(Знак минус у первого члена означает, что для нижней поверх­ ности угол между нормалью л н и вектором В тупой).

Интеграл по боковой поверхности s60 K с учетом того, что пбок ds = dl X (v At), преобразуем так:

\ В (к) пб0К

ds = —At ф В (4) (о X dl) = Atj) X В (t2)]~dî,

s6OK

1

1

где / — контур

основания

цилиндра.

112

Итак,

— J B(t2)7ids + \B(ti)nBdt+At§

[vxB(t2)]Il

= 0.

Добавим к этому равенству члены ± | В (tx) па ds, в которых

величина В взята в начальный момент времени t1. Перегруп­ пируем и разделим полученное уравнение на At, тогда

J В (t2)

п в ds — j В

(tj) nHds

J В (tt)

nu

d s - ^ B

(f{) nu ds

S B

f ] I

 

^11

 

f l l

I

 

M

 

 

 

M

"г"

 

 

-fcf) \vxB{t2)]dl

=

ü.

 

 

 

i

 

 

 

 

Первые две дроби представляют собой полную и частную про­ изводные, если величина Ar стремится к нулю

 

•^- \~Bnds—

-~-^Bnds

+ j> (üxß)uf/

= 0,

(98)

 

 

s

s

'

 

 

 

где s

— поверхность замкнутого

контура

I.

 

 

Обратимся

теперь

к уравнению (97) индукции. Левая его

часть

является

вектором, поток

которого

через

поверхность

s:

 

 

J Щ-пds

— J [rot (ÜX В)}

nds^O.

 

 

 

 

s

s

 

 

 

 

Так как подынтегральные функции непрерывны и дифферен­ цируемы, в первом интеграле можно поменять порядок интегри­ рования и дифференцирования, а во втором выполнить преобра­ зование по теореме Стокса. Тогда

-^- JBnds

<j) (vxB)di = 0.

(99)

s

'

 

 

Сложив уравнения (98)

и (99),

получим

 

£\BKds

= 0.

 

 

s

 

 

Это и является математической формулировкой нашей теоремы. Для формирования производных было использовано условие,

что жидкость состоит из одних и тех же частиц жидкости.

Диффузия вектора индукции В магнитного поля наиболее четко проявляется в электропроводной среде, когда мала величина

магнитного

числа Рейнольдса Re,„ =

ѵІ\ьтаа.

Это возможно,

в частности,

в неподвижной жидкости

или при

небольшой вели-

8 B . C . Бекнев

113

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

g - = v m A ß .

Известно, что такое дифференциальное уравнение описыв-ает закон распространения (диффузии) рассматриваемого параметра (в данном случае магнитного поля В) в неподвижной среде. Ис­ ходя из такой трактовки процесса диффузии, величину ѵ„, следует

 

 

рассматривать как коэффициент

диф­

 

 

фузии

магнитного поля.

!

 

 

 

Таким

образом,

распространение

 

 

магнитного поля в пространстве не­

 

 

может

происходить

мгновенно.

Д л я

 

 

примера можно отметить, что в мед­

 

 

ном шаре

диаметром

1 м

затухание

 

 

магнитного

поля длится

около

10 с,

 

 

а в том же объеме

морской

воды:

 

 

около

10~7 с.

При

диффузии;

т'. е.

Рис. 52.

Перемещение элемента

движении магнитного поля по элек­

тропроводящей

жидкости,

наводится

«жидкой» линии

электрический

ток,

который вызы­

 

 

вает нагрев среды джоулевым теплом. Таким образом,

энергия

магнитного поля переходит в теплоту, нагревающую среду.

Диффузия гидродинамического вихря со описывается анало­

гичным

уравнением

 

 

 

 

 

 

 

 

_ = ѵ Дсо,

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

Докажем важную теорему о сохранении циркуляции скорости по замкнутому «жидкому» контуру.

Теорема Томсона. Если массовые силы имеют потенциал, жидкость идеальна и течение баротропно, то циркуляция скорости по любому замкнутому жидкому контуру остается постоянной при движении жидкости.

Предположим, что среда и поле скоростей непрерывны. Для доказательства теоремы выделим в движущейся жидкости линию, состоящую из одних и тех же частиц, т. е. «жидкую» линию AB' (рис. 52). В силу непрерывности среды и поля скоростей через

промежуток

времени

А? линия

AB

может переместиться в поло­

жение А1В1,

причем

ее длина

может

измениться.

 

 

 

в

 

 

Рассмотрим интеграл / =

J vor,

являющийся циркуляцией

 

_

 

А

_

 

вектора скорости ѵ вдоль AB, где ôr берется вдоль «жидкой» линии.

114

Составим производную по времени от этого интеграла:

 

IL

в

' d v

fsl

I 7 d

(ôr)

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

г д е

 

 

 

 

 

 

 

 

 

= lim ^ -

6 7

^ l

i m

К - 7

J - (

' ' - ' )

=

= lim —

-irr1

=

lim - r

lim 1 .

= и Ü = OÜ.

Используя уравнение движения в форме Эйлера в векторной записи

do

-,

I ,

dT =

/ « - - p - g r a d p ,

преобразуем выражение (100) и получим, что

По предположению массовые силы имеют потенциал, т. е.

Im = grad U.

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

fmôT = grad c/Sr = % ôx + | - by - f ^ ôz = W;

grad p-ôr = ôp; u&v = -|- .

Подставляя эти выражения в формулу (100), получим

J i

или

W - U B — UA j — +

Таким образом, изменяется значение интеграла, взятого вдоль «жидкой» линии AB.

Если точки А и В совпадают, т. е. «жидкая» линия AB замы­ кается, то интеграл / дает циркуляцию Г вектора скорости по замкнутому контуру.

8*

115

В этом случае UA = UB;

vA =

ѵв

и тогда

dt -

y

p

(101)

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

В этом случае - ^ - = /(/?) и интеграл по замкнутому контуру

так как в одной и той же точке не может быть двух разных давле­ ний.

Следовательно, получаем - ^ - = 0, т. е. Г = const. Теорема

доказана.

Отсюда следует, что если в идеальной жидкости при баротроп- • ном и непрерывном течении и при наличии потенциала массовых сил в начальный момент времени не было вихрей, то их не будет

ив дальнейшем.

§14. СКОРОСТЬ РАСПРОСТРАНЕНИЯ МАЛЫХ ВОЗМУЩЕНИЙ

В сжимаемой жидкой среде могут возникать от тех или иных

причин возмущения параметров, которые характеризуют

эту

среду. Возникнув, эти возмущения будут распространяться

по

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

Важное

значение для

понимания процессов, происходящих

в жидкой

среде, имеют

малые возмущения параметров среды.

Под малым возмущением будем понимать такое изменение началь­ ного параметра среды (например, давления р), при котором абсо­ лютная величина изменения параметра (возмущение давления р') неизмеримо мала по сравнению с исходным его значением в данной

точке в рассматриваемый момент времени, т. е. р' <^

р.

В случае, когда абсолютная величина возмущения

соизмерима

с рассматриваемым параметром, возмущение называется конеч­ ным.

Малые возмущения распространяются по среде в виде волн. Например, распространение волны малых возмущений давления воспринимается нами как явление прохождения звука. Физически явление это происходит так. В состоянии покоя (t = /„) частицы в газе (рис. 53) распределены равномерно. Выделим в газе на­ правление и рассмотрим поведение частиц, например на оси х, когда в точке 0 возникнет повышение давления. Из-за сжимаемости частицы газа, расположенные вблизи точки 0, уплотнятся, а частицы, расположенные на некотором удалении, пока никак не

116

отреагируют. Возникшее местное уплотнение начнет

двигаться

в обе стороны с какой-то конечной скоростью, передавая

состояние

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

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

О

О

О

t "tj —О-О-О

о — о — ö *

с — с — — о

ооо*-

 

 

 

X

Рис. 53. Схема

распространения

звуковой волны во времени

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

Рассмотрим свойства малых возмущений с помощью общих уравнений движе­ ния электрически проводящего газа, находящегося в магнитном поле, характери­ зуемом вектором магнитной индукции В. Если не изучать вопрос, каким образом затухают малые возмущения, то можно отбросить все диссипативные явления. Будем полагать, что вязкость жидкости равна нулю (ѵ = 0), ее электрическая проводимость бесконечна (сг -> оо), газ нетеплопроводящий (Я, = 0) и все тепловые потоки Q отсутствуют; так как проводимость бесконечна, отсутствует и выделение джоулевой теплоты Qy. Эти ограничения позволяют записать необходимые основ­ ные уравнения:

уравнение движения с учетом выражений (62) и (51)

 

ой-

- —

1

1 1

_ _

 

 

_ _ j _ ( ö v ) ü = _

g radp + p ^ ( r o t ß ) x S ;

(102)

уравнение

(24) неразрывности;

уравнение

энергии

dsldt=0\

уравнение (97)

индукции;

уравнение

(4) Максвелла.

 

 

 

І17

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

ѵ=о

+

ѵ', р =

р0 + р',

р =

р0 +

р', (

(

ftna

=

(Рта) О +

*W

В = В0

+

В'.

 

Для простоты будем рассматривать неподвижную жидкость п = 0), а поле В0 •будем считать лежащим в плоскости дг, у, так что В0 ох, В0у, 0). Заранее никаких ограничений на возмущения накладывать нельзя, так как их следует считать функциями координат х, у, г и времени t, поэтому векторами возмуще­ ния скорости и магнитного поля будут v' [у'х, v', v'^j и В' (В'х, В', B'J

Система пяти записанных выше уравнений позволяет определить:

1)какие типы малых возмущений могут существовать и каковы их свойства;

2)каковы скорости распространения малых возмущений.

Для этого перейдем от уравнений, где неизвестными являются

параметры

среды,

к уравнениям, где неизвестными

будут малые

возмущения

этих пара­

метров.

 

 

 

 

 

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

что исходя

из уравнения

состояния,

записанного

в форме

р = р (р, s), можно

получить дифференциал

давления в

виде

* - ( $ ) , * + ( £ ) . * •

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

/ Яг> \

(104)

= const

С учетом уравнения (92) и условий (103) запишем используемое в дальнейших преобразованиях соотношение

dp' = a1 dp'.

(105)

Система уравнении для определения возмущения получается из пяти исход­ ных уравнений при подстановке в них выражений (103):

-g-, = - - p g r a < h p ' +

'

(гЫД')ХД,;

(106)

Ol

Ро

 

 

Ро крта) о

 

 

^ = - р 0 п і ѵ й ' ;

(107)

 

- | ^ - =

rot(u'Xßo);

(108)

 

div

 

В' =

0,

(109)

Особенности данной подстановки заключаются в следующем: во-первых, при дифференцировании исчезают постоянные значения невозмущенных параметров, например в уравнении (102)

rot В = rot 0 + В') = rot В0 + rot В' = rot В',

при этом rot В0 = 0, так как В0=. const;

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

118

Например, в том же уравнении

(102) имеем

 

 

 

 

 

 

 

 

 

 

 

1

( r o t ß ' ) X ß .

 

 

 

 

 

 

 

 

( r o t ß ) X ß ' = • РН-п

 

 

 

 

Это произведение состоит из суммы членов структурно подобных следующему

члену:

 

 

 

 

 

 

 

 

 

 

 

 

 

1

дв'х

 

1

 

 

дВ^

 

 

 

 

 

РЦя

W B x ~

 

(Po+P')[(^ma)o +

l * m . ]

' S f (В<*

+

В*)-

 

Разложим в ряд

 

 

 

 

 

 

 

 

 

 

 

 

 

|(Р0 +

Р') [0Л„а)о +

V-'тЛ Г

=

[Ро 0*т«)о

+

 

 

 

+

РоИша +

Р' ( Л « а ) о + Р > т а Г =

 

1

 

 

 

 

 

Po (llma)o

 

 

 

 

 

 

 

 

 

1

В

 

 

 

 

 

 

 

 

 

 

 

 

Ро (Цта)о

Ро (Hffla)o

 

 

 

 

 

 

где е — сумма членов

первого и второго порядка

малости.

 

 

 

 

После отбрасывания

членов с порядком малости выше

первого:

 

 

1

дВх

 

 

1

 

дВх

Box

+

 

 

 

Pollma

ду

 

Po (Мліа)о

дУ

 

 

 

 

 

 

 

 

 

и окончательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(rot ß ) X ß

1

 

( r o t ß ' ) X ß 0 .

 

 

 

 

 

Робота

Ро

 

 

 

 

Проектирование уравнений для возмущений на координатные оси проведем

с учетом следующего упрощающего предположения. Будем считать,

что малое

возмущение

распространяется плоской волной.

 

 

 

 

 

 

Это означает, что в любой момент времени во всех точках любой

плоскости,

перпендикулярной

к направлению движения

волны (например,

к оси х), все

 

 

 

 

 

 

 

 

 

 

 

д

 

д

параметры сохраняют постоянное значение, т. е. производные —^- и - ^ - равны

нулю.

 

 

 

 

 

 

 

 

 

В

результате

получим такие уравнения:

 

из

уравнения

(106) движения

 

 

dBи

 

 

 

дѵх

1

dp'

В,

 

 

 

dt

Ро

дх

 

âx

PO(LI„

 

 

 

 

 

 

 

о у

 

 

 

 

 

 

дѵу

 

 

ß,0.r

 

dB и

(ПО)

 

 

 

dt

Po(llma)o

dx '

 

 

 

 

 

 

 

dv'z

 

 

ßp.t

_

dBz

 

 

 

 

dt

 

PoOma)o '• дх

 

из

уравнения

(107)

неразрывности

 

 

 

 

 

 

dt

 

 

дѵх

 

(111)

 

 

 

 

 

 

 

 

11&

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