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

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

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

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

относительно магнитного поля. В данном случае пондеромоторная сила F, дей­ ствующая на этот заряд,

F = q(vXB),

что примем как опытный факт.

Магнитное поле в материальной среде характеризуемое вектором В магнит­ ной индукции, в отсутствие тока смещения слагается из двух полей (причем, и в том и в другом случае в их основе лежат движения каких-либо электрических зарядов):

1. Магнитного поля, приложенного извне, а также магнитного поля, создан­ ного движущимися по объему, занятому веществом, свободными электрическими зарядами (электрическим током);

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

 

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

 

 

 

В = цтН

+ 7,

где

H — вектор напряженности магнитного

поля (характеризует первое поле).

 

Если магнитное поле рассматривается в вакууме, то намагничение отсутствует

и вектор магнитной индукции В =

u m o # ,

где ц т о — коэффициент магнитной

проницаемости вакуума.

 

 

 

 

 

В системе СИ

 

 

 

 

 

 

 

ц.,„0

= 4я-10"7

=

1,256637-10-« Гн/м.

 

Для ряда сред между

векторами

/

и H существует связь:

где kM — коэффициент

магнитной восприимчивости.

 

Если kM^> 0, то

намагничение усиливает внешнее поле (парамагнетизм),

если

ku < 0, то ослабляет

(диамагнетизм).

Наконец, для ферромагнитных тел

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

 

Если для материальной

среды

ввести

соотношение

 

 

 

 

В = Ѵ-таН,

 

 

то

получим, что

 

 

 

 

 

где

Дта коэффициент абсолютной

магнитной

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

данной среды.

 

Удобно

считать, что j . i m a

= |xm fxm o , где

ц т

— коэффициент

относительной

магнитной

проницаемости среды.

 

 

 

 

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

Электрический ток является направленным движением (переносом) электри­ ческих зарядов. Если носители электрических зарядов qk движутся с индивидуаль-

20

ными

скоростями

Vk, то для объема

ДѴ за вектор плотности электрического

тока

принимается

величина

 

 

 

 

-

,.

к

 

 

1 =

1 1 1 1 1

" - Л І / • •

Иначе говоря, вектор у плотности электрического тока — это количество электричества, переносимого через единицу поверхности за единицу времени в на­ правлении, перпендикулярном этой поверхности.

Если для простоты считать, что все скорости частиц <?/г одинаковы, то вели­

чину -Ok можно вынести за знак предела. Тогда оставшаяся величина

S 4k

pe = I im —^j— Д1'->0 д і '

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

/ = РеѴк*

Отметим, что для рабочих тел КУ носителями электрических зарядов обычно являются отрицательно заряженные электроны и однократно ионизованные поло­ жительные ионы.

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

перемещения для одного

заряда, равного

будет величиной ^F/zV^,

где

 

 

 

 

к

сила F^

суммирование происходит по всем зарядам объема ДУ. Пондеромоторная

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

действующая

на

рассматриваемый заряд

 

 

 

 

?fc = <fe£ +

<7fc(ufcXß).

 

 

Величины Е и В в выбранной системе координат полагаем постоянными.

Тогда ^ F^v/t =

^q/cEvk,

так как смешанное произведение % (ѵ^ ХВ)

Vk=

О-

k

к

 

 

 

 

 

Работа электромагнитного поля над электрическими зарядами, заключенными

в единице объема, за единицу времени

 

 

 

 

lim

k . . . — = £ l i m

-^-тт}— = ~Ej.

 

 

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

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

rot7? = 7 +

 

(i)

r o t £ = _

f ;

(2)

div D =

pe;

(3)

div 1 =

0.

(4)

21

Эту систему скалярных уравнении дополняют соотношения:

 

 

 

 

 

 

 

(5)

 

 

 

В =

ц,І І а //.

 

(6)

Связь электрического тока с электрическим и магнитным полем выражается

обобщенным

законом Ома:

 

 

 

 

 

 

]=аЕ

+ о&хВ)

— -Ц-

(jxB)

+ aÇjxB)B

+ p,u

(7)

 

 

 

 

 

 

 

где каждый

член

отражает

определенное

физическое

явление, происходящее

в среде, движущейся со скоростью ѵ и обладающей удельной проводимостью о. Первый член в правой части выражения (7) является током проводимости, вызванным перемещением электрических зарядов под действием приложенного электрического поля Е\ второй член — индукционным током и отражает эффект взаимодействия движущейся электропроводной среды как целого с магнитным по­ лем В; третий член — током Холла, где е — заряд электрона, п — число заряжен­

ных частиц в единице объема (концентрация заряженных частиц).

Эффект

Холла

проявляется

в том, что электрический заряд, движущийся

в магнитном

поле,

отклоняется

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

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

Четвертый член (а — некоторый комплекс параметров) отражает эффект проскальзывания ионов относительно нейтральных атомов. В некоторых условиях такая разница в скоростях ионной и нейтральной компоненты для рассматривае­ мых КУ может иметь место.

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

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

впорядке их убывания.

Врассматриваемых І\У возможно упрощенное рассмотрение некоторых свойств электромагнитного поля без ущерба инженерной точности расчета.

Во-первых, характерное отношение

(у/с)2 (где с«=< 3-108 м/с — скорость

света) даже при скоростях газа порядка 103

м/с всегда чрезвычайно мало. Во-вто­

рых, в первом дифференциальном уравнении Максвелла можно пренебречь чле­ ном, называемым током смещения dD/dt, если сравнить его величину с величиной /.

В-третьих, в уравнении (7) можно пренебречь конвективной составляющей

электрического тока реи

в сравнении с током проводимости вЕ. Для большинства

задач, пользуясь

обобщенным законом

Ома, можно пренебречь токами

Холла

и эффектом проскальзывания ионов относительно нейтральных атомов.

 

Следовательно, в большинстве технических задач можно использовать обоб­

щенную форму

записи

закона

Ома

 

 

 

 

j =

a[E +

(vXB)].

(8)

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

Г Л А В А II

ОСНОВЫ КИНЕМАТИКИ ЖИДКОЙ СРЕДЫ

§ 4. ДВИЖЕНИЕ БЕСКОНЕЧНО МАЛОЙ ЧАСТИЦЫ ЖИДКОСТИ. ТЕОРЕМА ГЕЛЬМГОЛЬЦА

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

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

dl \M(X,lj,Z)

0(хо,Уо,го)

Рис. 13. Бесконечно малая: «жидкая» частица

Пусть

в точке О проекции скорости

равны

vx0,

vIJ0 и

 

ѵгог

тогда в соседней точке M с координатами

х = х0

- f

dx, у

=

у0

+

+ dy\ z

— z0 + dz проекции скорости будут равны

ѵх,

иу

и

vz.

Разлагая проекции скорости в точке M в ряд Тейлора с точностью

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

dx = £,

dy

=

ï],

dz = £, получим

 

 

 

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

" . - - + (%). s+ Ф ) . " + (Зг).ь

 

 

 

 

 

 

 

 

 

 

 

 

 

(9)

 

* = ' - + (тЁ).е + (тг).і + ( т ) Л

 

 

 

 

причем

ду

дц

II т. п.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

Эти равенства в силу малости частицы (£, т], Z, — бесконечно

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

 

Чем

меньше

т], £,

тем

точнее

выполнятся

равен­

ства (9).

 

 

 

 

 

Допуская известное неравенство правой и левой частей выра­

жений

(9), например,

равное

1%, 0,1% и т. п., тем самым опре­

деляем

размеры частицы.

 

 

 

 

Произвольная деформация «жидкой» частицы в рассматривае­

мом случае сводится

к растяжению

(сжатию)

«жидких»

отрезков

и к изменению углов между двумя «жидкими» отрезками. Искрив­ ление «жидкого» отрезка исключено в силу принятой линейной зависимости проекций скоростей от координат £, т), £. Поворот плоскости «жидкого» угла будет учитываться при оценке враща­ тельного движения частицы как твердого тела.

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

Деформация

«жидкого»

отрезка

 

 

 

Рассмотрим движение элементарного отрезка длиной dx, выде­

ленного мысленно в бесконечно малой жидкой частице.

 

Для простоты анализа оси координат

расположим так,

чтобы

ось X проходила

по отрезку, а начало отрезка совпадало бы с на-

у\

 

 

чалом координат (рис. 14).

О на

 

 

 

Если проекция скорости в точке

 

 

ось X равна

ѵх0,

то для точки А по

фор-

ог

А

А г

мулам (9)

получим

 

Рис. 14. Деформация (dy я dz в данном случае равны нулю), «жидкого» отрезка Через бесконечно малый промежуток времени dt отрезок OA переместится в но­ вое положение Ох Л х , а его проекция на ось х в положение 0 2 Л 2 .

причем в соответствии с формулами (9) он останется прямоли­ нейным, изменившим лишь свою длину и положение в про­ странстве.

За время dt точка О в направлении оси х пройдет расстояние vx0 dt, точка А отрезка OA за то же время переместится на рас­ стояние ^ ѵх0 + - ^ г dx^j dt.

Вычитая из последнего выражения первое, получим линейную деформацию отрезка dx в направлении оси х в виде - ^ - dxdt.

Деля линейную деформацию на dx, получим относительную линейную деформацию, а деля ее на dt, получим скорость относи­ тельной линейной деформации отрезка в направлении оси х в виде гх = дѵхІдх (здесь и далее индекс О опущен).

24

Рассуждая аналогично для отрезков dy и dz, имеем скорости относительных линейных деформаций отрезков в направлении осей у и z:

 

дѴг

ду

dz

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

Деформация

«жидкого»

прямого

угла

 

 

 

 

 

 

 

 

 

 

Выделим

мысленно

в

жидкой

частице

прямой

угол

А OB

со

сторонами dx и dy и рассмотрим

его

деформацию за

 

время

dt.

 

Для упрощения анализа расположим оси координат так, чтобы

оси X и у совпали со сторонами угла,

а начало координат совпа­

дало

с вершиной угла (рис. 15).

 

 

О

 

 

ѵх0,

 

ѵу0

 

 

vz0,.

то

Если проекции скорости

в точке

равны

 

и

 

в точке А

 

в соответствии с формулами

(9),

получим

 

 

 

 

 

 

 

 

 

Зѵ V

 

 

 

 

 

dx\

 

 

 

 

 

 

 

 

 

ѵхА

= ѵх0

+ -g§- dx- V,

 

'У° 1

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'zO '

dx

dx.

 

 

 

 

 

 

 

 

 

 

 

а

в точке

В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

JxB

 

 

 

дѵх

 

 

 

дѵу

dy,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ду

 

 

 

 

 

 

 

 

 

 

 

 

 

V

zB

20 '

dvz

dy.

 

 

 

 

Рис.

15.

 

Деформация

 

 

 

 

= V.

ду

 

 

 

 

 

 

«жидкого» прямого

угла

 

 

 

 

 

 

 

 

 

 

 

 

 

За бесконечно малый промежуток времени dt прямой угол А OB

переместится

 

в положение J

4 1

0 1 ß 1

,

а его проекция на

 

плоскость-

хОу

в положение А202В2,

 

причем в общем случае угол

прямым

не

останется.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

хОу

 

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

 

за время dt,

которое складывается

из углов dyy и

dy2.

 

 

 

 

 

 

Из

рис.

15 следует,

что при малых

углах

 

 

 

 

 

 

 

 

 

 

 

 

 

ПіАа

 

 

 

dy2

tg (dy2)

= пВ

2

 

 

 

 

 

 

 

 

 

 

тОп

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее,

вычислив

длины

отрезков

тА2,

т02,

nBz

и

/г02 ,.

получим

 

 

дѵ,

 

 

 

 

 

 

dvx

 

 

 

 

 

 

 

 

 

 

 

 

dx

dt

 

 

 

 

dydt

 

 

 

 

 

 

 

 

 

dyx =

dx

 

 

 

dy.2

_=

~df

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx-

dv*-dxdt

 

 

 

dy-

dv ü-dydt

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

25-

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

dyx. Щи. dt;

d y ^ ^ d t .

дх

'

'*

ду

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

.формулой

 

dt

 

 

дх

^

ду

'

 

 

Обозначим половину

скорости

изменения

прямого угла

через

Ѳг, тогда

 

 

 

 

 

 

 

 

 

 

0I

_

_!_

(Еа±\

дѵ*

 

\

 

 

г

 

2 \ дх "т" ду )

 

 

Проведя аналогичные рассуждения для двух других плоско­

стей, получим

 

 

 

 

 

 

 

 

 

 

fl

JL

(ÈLl-lÊ!l!L\

)

 

 

х

~

2

\

ду

~f

dz

'

 

 

О

—-L(Èk.

 

- L * !

l V

І

(10)

fl — _і_

і(

àvu

_l_ d"

\

 

 

 

D z

~

2

\

дх

 

ду

)

 

 

 

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

Вращательное движение частицы без изменения формы

Как уже было сказано, движение жидкой частицы отличается от движения твердого тела наличием деформаций.

Рассмотренные выше скорости деформационного движения можно ввести в выражения для проекций скорости произвольной точки частицы жидкости.

Тогда прибавляя и вычитая члены ^ " ( " ^ j " ) 1 ! и ~іг(~~^")£ из первой строчки выражений (9), получим

Т) — 71

дѵх

t

, дѵх

m

I дѵх

y _ 1

дѵу ^

і

1 дѵ.

Т

х

Р - L

х

п

-1-

х

Г -и —

и

"

п

-4- —

z

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ох = ѵх0

+ гЛ + Ѳ,С +

 

+ 4 - ( ^

 

-

 

 

) С -

 

•26

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

Слагаемые типа ъх\ - f + Ѳгг| представляют собой проекции на оси координат скоростей деформационного движения, состоя­ щего из растяжения (сжатия) отрезка длиной \ и изменений пря­

мых углов в плоскостях

xOz и

хОу.

 

 

Оставшиеся слагаемые

типа

 

 

 

l \ dz

дх I =

- 2 V дх

ду ) Л

(12>

 

могут выражать собой только скорость вращательного движения частицы жидкости как твердого тела, так как деформационное движение уже учтено.

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

полюс с

угловой

скоростью

со (см. рис.

 

13).

 

 

 

В этом случае скорость точки равна векторному

произведению

G) на dl,

т.

е.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

І

j

 

k

 

 

 

 

 

 

 

ü B p

=

(<DX<tt) =

t СО^ сог

 

 

 

 

 

 

 

 

 

 

 

 

I

T1

 

£

 

 

 

 

Проекции скорости

ивр

на оси координат

имеют вид

 

Таким

образом,

 

 

 

 

 

 

 

 

 

 

 

 

 

(

дѵхвР

дѵгпр

\

_ .

1

(' düyap

 

дѵхвр

\

_

со,

 

\

dz

дх

 

 

 

2

V дх

 

ду

j

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Vx =

Ухо +

Zxî +

%£, +

Q2il +

®yt, — «W. )

 

(13)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уг =

Уго +

е2С +

Ѳѵ.і-і +

Ѳу? - f

û^n — <fl

 

 

 

В то же

время

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

(iîr ~

i f )

=

 

~т (%+<»*—%+ау)

= % и

т- п->

т. е. выражения типа (12) представляют собой проекции скорости вращательного движения «жидкой» частицы, что и требовалось показать.

27

Поступая аналогично для осей і/ и z, получим

CD,.

1

(

дѵг

dz } '

2

\

ду

 

_

_ L ( дих

(14)

соУ~~

2 V dz

дх ) ' I

2 V дх

Следовательно, скорость произвольной точки частицы жидко­ сти можно записать в виде

Ѵх — Ѵх о " h Ѵх д е ф + <-'л пр-

Пользуясь круговой подстановкой, получим для проекций скорости на оси у и z:

или в векторной форме

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

§ 5. КЛАССИФИКАЦИЯ ДВИЖЕНИЙ ЖИДКОСТИ

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

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

ческим условием

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

является

тождество

 

 

 

со =

0.

 

 

 

 

 

Следовательно, условие потенциальности требует, чтобы поле

^скоростей в жидкости

подчинялось

равенствам,

 

 

 

 

x

= Loy =

coz =

0

 

 

 

-лли, согласно выражениям

(14),

 

 

 

 

 

 

дѵх

дѵу

т дѵх

дѵг

.

дѵу

dv

z

(15)

ду

дх

'

dz

дх

'

dz

 

ду

 

28

В тех областях потока, где со Ф 0, движение называется вихре­ вым или непотенциальным,

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

 

Зф

 

ôcp

 

оф

 

 

Ѵх-~Ш'

V»~W;

Vz~~te'

 

 

В самом деле, по определению

 

 

 

;

j

*>х = ЧѵО + гхх

+ Q,jZ -f- Qzy =

ѵв =

ѵуо+гуу

+

Вгх +

^

= ^ ;

\

(16)

=

Ü20 +

F«Z +

ѲѴГ/ +

0J,A;

=

.

 

Подчеркнем, что величины е и 0 в этих выражениях постоянны и относятся ж полюсу частицы. Координаты точки £, г), Ç заменены на А', (/, г.

Для построения функции ф проинтегрируем первую строку выражений (16)

по X, тогда

 

ф = ѵх0х + ~- 8v.ï3 + Qyzx + Qzyx + f{y,z).

(17)

Теперь продифференцируем полученное выражение для ф по у и z и, исполь­

зуя выражения (16) и (17), найдем

Откуда

jjjj = fi/O + Byy -f ѲЛг; -^- = ад + ег г -f Ѳ^.

или, интегрируя, получим

/ = ѴуОУ + \ - ЪуУ~ + Ѳд-zf/ + f1 (г) ;

f = адг + ~ sz z2 + Ѳ.ѵ(/г + ^ (у),

что дает

ѴуОУ + 4 ~ ey<r /а (У) = "гог + -^- ег г2 — /х (г).

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

Ѵуоу + -^- бу(/2 — /а (у) = — с,

29

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