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

книги из ГПНТБ / Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие

.pdf
Скачиваний:
14
Добавлен:
23.10.2023
Размер:
4.12 Mб
Скачать

Г л а в а 4

ОДНОМЕРНОЕ ТЕЧЕНИЕ ГАЗА В ЭЛЕКТРОМАГНИТНОМ ПОЛЕ

4-1. Основные уравнения

Основные уравнения выводятся при следующих допу­ щениях:

а) движение газа считается установившимся (d/dt = 0); б) течение одномерное, хотя стро­ го говоря одномерность следует рас­ сматривать только к средним по се­

чению параметрам; в) существует только одно внеш­

нее поле — электромагнитное (грави­ тационным полем (Ггр= 0) пренебре­ гаем, пренебрегаем также силой вяз­ кости и теплопроводностью.

Исходными уравнениями являются

уравнения

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

количества

движения, сохранения энергии, урав­

нение состояния:

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

а)

уравнение

Fpu = const.

 

J _ dF_ ,

1 dp

+

da

=

0.

( 4 1 )

, 1

 

n

F dx +

dx

и

dx

 

 

 

Рис. 4-1. Струйка то­ ка в магнитном поле

Скорость с имеет составляющую, на­ правленную вдоль оси х или с (и,0,0);

вектор магнитного поля В направлен

вдоль оси z, т. е. В (0, 0, В ) и вектор электрической напря­

женности Е направлен вдоль оси у, т. е. Е { 0, Е, 0) (см. рис. 4-1);

72

6) уравнение импульсов для одномерного течения в электромагнитном поле

“ -7L + - - £ -

+ * . = о,

ах

р

dx

 

ХУ1—проекция Ры на

ось

ж. Пондермоторная сила равна

FM=.~JxB = ° 0 +~cXB\XB,

векторное произведение равно

 

 

 

с Х В

=

i

/

k

 

 

 

и 0 = -

и Bj,

и

0

0

= ( - 1 П

 

0

0

В

 

 

 

0 В

 

 

 

 

 

 

 

а плотность тока /

выражается как

 

 

Тогда

 

 

 

j = a[E uB]j.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

/

k

 

~РЛ = а [В — и В] j X В

 

0; а — мВ]; 0 — аВ[Е иВ] i,

 

 

 

 

 

0

 

0

0

 

т. е. магнитная сила действует

в

направлении

оси х. Учи-

cl

д

 

 

 

 

 

 

 

 

тывая, что —

■— = 0, можно записать

 

ду

dz

 

 

dp

 

 

 

 

 

du

 

 

+ а[Е — иВ]В;

(4-2)

 

ри ---- = -

 

dx

 

dx

 

 

 

 

 

 

в) уравнение энергии

 

 

 

 

 

PU~dv[i +

 

 

 

 

В^ Е '

 

Так как i— cpT, при c„=

const получаем

 

Со ои dx

+

ри2 dx =

а [Е — и В\ Е;

(4-3)

г) уравнение состояния. Считая газ совершенным, будем

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

Клапейрона p — pRT или в диффе­

ренциальной форме

 

 

 

 

 

t

dp

1

dp

1

dT

(4-4)

р

dx

р

dx

Т

dx

 

73

В уравнении энергии (4-3) заменим теплоемкость с„ через газовую постоянную и показатель адиабаты. Изве­ стно, что

сР

 

си—i?,

или

 

tC

 

____

Т~)

- •

Ср

R

 

 

 

к —1

С учетом последнего выражения уравнение (4^3) прини­ мает вид

ри

R

+ риа ——— а[Е иВ\ Е.

(4-5)

к — 1

dx

dx

 

Подставляя — из уравнения (4-4), находим

Тdx

 

ри

RT

1

dp

{

I d

/

7 .

1

du

 

 

 

 

'

!

“Г

z

;

Т

и

dx

 

 

 

 

к — I

 

р

dx

 

F

dx

 

 

 

 

 

+ pиг

= ст (E — и В) E.

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

Из уравнения движения (4-2)

выразим

 

и

подставим в

последнее выражение. Тогда имеем

 

 

 

 

 

 

 

 

к- D l l ----Р !LJ!lL + J - a(E UB ) B +-

dF

 

 

dx

 

 

к — 1

 

р

dx

 

р

 

 

 

 

F

 

 

I

1

du

 

 

du

= a[E -

и B\ E-

 

 

 

 

i

и

dx

+ ри* —

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

1

, .

к

 

\ du

 

к

 

1

 

dF

 

 

----------pu2

-|--------- p

----- = -----------pU—-----------

 

 

к —1

 

 

к —1

j dx

 

к — 1

F

 

dx

 

 

 

 

■a (E и В)

 

■uB — E ) .

 

 

 

 

 

 

 

 

 

 

к — 1

 

 

 

 

 

 

 

 

 

и2

. Л du

 

 

к

 

J _ dF_ _

 

 

к — 1

 

k R T +

) ~dx

 

 

к —■1

 

F

dx

 

 

 

a(E uB)

^

u B — E^.

 

 

 

 

; -----!---- [JL i L

+

a(* ~

(E — uB) l

u

B

— E

 

dx

уИа — 1

[ F dx

 

p к

 

 

 

Vк — 1

 

 

 

= — 1— Гл J L + ^

 

( * - - a) (a -

к -

\ E_

(4-6)

dx

 

к В

Мг — 1 |_ F dx

 

p \ В

 

J \

 

 

74

Приведем уравнение (4-6)

к

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

виду. Так

как и= Ма,

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

da

л л

da

,

 

dM

 

 

 

 

 

 

 

---- = М ------- b а ---- .

 

 

 

 

 

 

 

dx

dx

 

 

dx

 

 

 

 

В свою очередь скорость звука

 

 

 

 

 

 

 

 

 

 

 

 

da

 

 

 

V k R T

dT

 

 

 

a = V ' k R T и - j - =

 

2

T

 

dx

 

 

 

 

 

 

 

dx

~

 

 

 

 

da _

 

и

 

V k R T dT + V k R T dM

dx

V T r t

2T

 

dx

 

 

 

dx

 

 

dM

 

 

1

du

 

 

 

 

1

dT

 

 

 

dx

 

V k R T dx

 

 

2 V k R T T dx

 

 

d

' dT

 

 

 

 

 

 

 

 

 

.

находим

Выражая —- —— t из уравнения

энергии (4-5),

dM

 

 

 

da

 

 

 

 

к — 1 1 p ir

d'u

dx

V k R T dx + 2 V k R T

к pu

dx

-

и

 

к — 1 1

r „

 

„ , „

 

1

i +

.___ - —

-]—

a\E — u B \ E -

r____

2 V k R T

 

к

pa

 

 

 

1

 

V k R T

 

 

 

P«a \

da

 

1

 

 

ЛГ— 1 £

 

 

и

 

2p j

dx

 

 

 

 

к

В

p

 

 

 

 

 

 

 

 

 

 

1 +

h: — 1

аг

1

, / .

 

 

 

 

 

 

 

 

1

4“ (л*

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n r

 

k — 1

E

 

 

а, a

E

через

иг.

_

иоозначим — ---- — через

Тогда

имеем

 

1

 

К —1

 

 

 

 

 

 

 

dM

 

 

м /

 

 

a

dF

 

dx

V k R T (■\ +

 

/Иа-- 1

F dx

 

аВ2 и3) (и Ui)

 

 

 

ай2 («з — u)u1 =

 

 

 

 

 

 

 

2 V k R T

р

 

 

 

 

1 +

 

к —1

 

 

 

 

 

 

 

 

 

1

 

--------/VI2

 

a

dF

ай2

,

 

. ,

,

 

 

к

 

 

V k R T

 

М2 — 1

— —------------ (и — Из) (“ — Ui) —

 

 

F dx

р

 

 

 

 

 

 

 

М"-— 1

 

1

ofl2 (u3ut) и,

 

 

 

 

1+

* r i M.

2

 

f

 

 

 

 

 

75

об2

—-ы3)

и — ы.

/VI2 - 1

 

- к - 1М2

 

об2 (« — Ms)

1 +

 

и — Ыг 1+

Ма — 1

1

 

 

 

 

- — 1 Л'12 2

2

М2 —1

2 + /сМа - А12 + М2 — 1

1 + 2 + (к - 1) /И2

2 + - 1) М2

Обозначим

1 + к М2

2 + (/с — 1) Л'12

 

 

 

 

Мо —

1+ к/И2

и*

 

 

 

 

 

2 +

( к — О М 1

dM

_

1.+—2

М= Г Д-I

//б

 

об*

= (и — м3) (и — г/2) . (4-7)

d.v

_

Ма — 1

1

rf.v

 

p V k R T

 

 

 

 

б

 

Для простоты

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

сечения,

т. е.

 

dF

п

 

 

 

---- =

и

 

 

 

 

 

1

 

dx

 

 

 

 

dM _

 

 

1+

Л4!

об2

(и — и3) (и и.,). (4-7а)

dx ~

 

М* - 1

 

 

 

 

 

 

 

Vк R T p

Анализ течения газа в электромагнитном поле

Течение можно разделить на дозвуковое и сверхзвуко­ вое. Рассмотрим уравнения (4-6) и (4-7,а).

В дозвуковой области

В сверхзвуковой области

(М <1),

(М>1),

112<11\<U3,

Ui<ii2<ii3,

1/1

при М — 1; U2 — u\,

при М — 0; U2 — .

 

п р и М = оо и? —> — - — и,

 

к — 1

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

76

разбить

на 8 областей

(А \; А 2; 5i; В2;

С\; С2; D

Рассмотрим подробно эти области (рис. 4-2).

Область А\ (М >1)

Область А 2 (М <1)

 

и > щ М >1.

 

и > и 3,

Из ур-ния (4-7,а)

 

> 0.

 

dM/dx<0.

 

 

 

dx

Поток тормозится.

Поток ускоряется.

 

Зона В\

Зона В2

 

112<и<Щ.

U \ < U < l l 3 ,

 

Из ур-ния (4-6)

± ^ < 0

Ш < 0,

 

dn/dx>0.

dx

dx

Из ур-ния (4-7,а)

поток тормозится.

 

dM/dx>0,

 

 

поток ускоряется.

 

 

Переход через скорость звука в канале постоянного

сечения

в областях А\ и А 2 невозможен.

В областях В\ и

Рис. 4-2. Период через скорость звука при течении проводника в магнитном поле

В2 мы даже не можем приблизиться к скорости звука. Переход через скорость звука возможен только на стыке

областей В2 и А\ при

через особую точку. Имеем ана­

логично с соплом.

Область С2

Область С1

U \ < U < U 2.

u2< w< w1,

77

Из уравнения (4-6)

>

0,

du/dx> 0.

dx

 

Из уравнения (4-7,а)

dM/dx<0.

dMjdx СО.

 

 

Б этих областях поток

ускоряется при

уменьшении

числа М. Здесь скорость звужа растет быстрее скорости потока, так как выделение тепла Джоуля приводит к росту

Т, а так как a—Y kRT, т о

а растет быстрее, чем с.

Область D\

Область £>2

U < U \ ,

и < u%

— < 0,

du/dx>0,

dx

 

dM/dx<0.

^dx- > 0.

Переход через скорость звука на границе D\ —£>2 невоз­ можен. Можно перейти через скорость звука в особой точке из области £>2 в .Вь Возможен также переход из области D 1 в С%

Полученные данные можно пояснить, если вспомнить, что воздействие электромагнитного поля на течение выра­ жается в силовом воздействии пондермоторной силы (Ём) и в выделении тепла <2ДШ. Отношение работы механиче­ ских сил £м-с к полной подводимой энергии в данном случае равно параметру нагрузки ц, выражаемому как

_ F и _ а иВ) Ви _ иВ _ и

^

N

а иВ) Е

Е

«з

Если u>uz, то N<0; Fu<.0, т. е. механическая энергия при движении газа переходит в энергию электромагнитного поля. Если и<из, то энергия электромагнитного поля пере­ дается газу в виде тепловой и механической энергии. В последнем случае при и-*-из (т]-*-1) воздействие поля выражается в виде работы электромагнитной силы над газом, а при ц малом —в основном в виде подвода тепла.

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

78

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

Здесь уместно привлечь на помощь закон обращения воздействий газовой динамики, впервые сформулированный Л. А. Вулисом.

4-2. Закон обращения воздействий

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

силы совершают работу трения «dLTp».

привлечь:

Для анализа общего случая необходимо

1.

Уравнение

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

в

дифференциальной

форме

 

 

 

 

 

 

 

dF

1

dp

1

dc __ q

 

 

F dx

p

dx

2

dx

 

2.

Уравнение

сохранения

энергии

в форме первого

начала термодинамики, которое после

несложных преобра­

зований имеет вид

 

 

 

J_JP_ ц _

d Lmp +

d Lmex — 0.

 

рdx

3.Уравнение состояния идеального газа

п /г- 1

dp

1 dp

1 dT

=

n

p — p R T ------

dx----------------------------

p dx

T dx

0.

p

 

 

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

(

др

 

\

=

(

д Т \

находим

 

 

 

уравнении связи

\

----

 

 

----

j s

 

 

 

 

ds

 

) р

 

 

ov

 

 

 

 

 

1

du

_

1

dF

 

 

i

dm

1

d Lmex

+

и

dx

_

F

dx

 

 

m

dx

a2

dx

 

 

 

+

d Lmp

 

 

к — l dQ

 

 

(4-8)

 

 

 

 

dx

 

 

a2

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение (4-8) носит название уравнения воздействия.

79

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

компенсируются.

При

переходе через скорость

звука

(М = 1) одно и

то же

воздействие оказывает на

поток

обратный эффект (например, подвод тепла к газу в дозву­ ковом потоке приводит к возрастанию скорости, а в сверх­ звуковой — к убыванию).

Закон обращения воздействий: для перехода через ско­ рость звука необходимо изменить знак внешнего воздей­ ствия.

Для изменения давления уравнение воздействия имеет вид:

(М2-

1) J - J E - -

-

к

М2

dF

1 й LmeK

4

р

dx

 

 

F

dx

a2 dx

к

dZmp

к — 1

dQ

1

dm

(4-8a)

a2

dx

as

 

dx

m

dx

 

 

Для изменения плотности:

 

 

 

(Л42 -

1) -L-^P-

 

Г М 2

dF

1

d Lmex

 

р

dx

 

[ F

dx

a2

dx

 

 

к — 1

dQ

1 dm 1

(4-86)

 

 

a2

dx

m

dx ]

 

 

 

Под a^"‘ex можно подразумевать любую техническую

dx

работу (турбина, компрессор, работа магнитной силы). Частные случаи.

а) Движение газа при наличии одного геометрического воздействия (рис. 4-3).

 

( М 2 _

1) _ L

J 0 .

1

dF .

 

 

и

dx

F

dx

dti .

 

 

 

 

, . __ ,

----oo, но достаточно велика,

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

dx

 

 

 

 

 

 

 

= 0.

 

 

 

 

dx

 

 

 

6) Движение газа при наличии

теплообмена (рис. 4-4)

dF

_ dm

__ d L-mex __ ^ Lmp __ q

dx

dx

 

dx

 

dx

80

Существует только одно тепловое воздействие

(-'И*

1)

1

du

dx

 

 

 

 

и

dx

аz

 

В дозвуковом

потоке возрастание

скорости

связано с

подводом тепла

dQ

.

«

 

 

dQ

п

----

> 0, а в сверхзвуковом ----

<*. U.

 

dx

 

 

 

 

dx

 

 

М <о

 

 

 

d£>o

 

 

 

d x

 

 

 

d x

 

 

M‘l

W

/ У / У / , ' / / / / / / /

М>/

 

 

M <f

W

7

' \

\ \ \

0 > о

 

с& * 0

dec

 

-

Рис. 4-4. Тепловое сопло

Следует обратить внимание на множитель, стоящий перед тепловым воздействием: он очень мал (для воздуха

к = 1,4)

^ 0,4-10-3.

6 -5 9 9

81

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