книги из ГПНТБ / Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие
.pdfГ л а в а 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 (u3— ut) и, |
|
|
||||||
|
|
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 |