Файл: Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 05.08.2024
Просмотров: 68
Скачиваний: 0
не. Результаты расчета для случая канала постоянного сечения представлены на рис. 4-13. Расчет был выполнен
для |
паров цезия |
при |
2’оо = 2000°К |
и 5000° К |
при |
дозвуко |
вом |
течении М= |
0,3. |
Статическое |
давление |
на |
входе в |
генератор равно 0,1 атм.
На рис. 4-13 представлены зависимости внутреннего относительного к. п. д. r|0i и коэффициента полезного пре-
Рпс. 4-13. Дозвуковые рабочие характеристики МГДпреобразователя с каналом постоянного сечения, но проводимостью, изменяющейся по длине канала
образования энергии T)lY о т т)э Д л я различных относитель ных длин генератора
m
На графике пунктиром для сравнения даны аналогич ные характеристики для постоянной проводимости. Для М= 0,3 скорость газа возрастает вдоль длины канала гене ратора, при этом его температура падает. Уменьшение тем
102
пературы плазмы является причиной уменьшения местной проводимости плазмы вдоль канала и ухудшения рабочих характеристик генератора. Так же как в случае постоянной проводимости к. п. д. преобразования Г|* ограничен усло вием запирания (кривая 1). Величина т]а, при которой т]*
0 1 Z 3 8 5 6 7 8 3
Безразмерная длина канала
Рис. 4-14. Зависимость плотности мощности от безразмерной длины канала о
достигает максимума, больше для переменной а, чем для случая а= const.
Для канала постоянного сечения средняя плотность мощности N может быть определена в соответствии с
уравнением (4-35) в виде |
|
|
-г, |
т |
г |
|
к —1 |
или заменяя расход т через безразмерную длину б
B q GO III
получаем
ii ~ т Д т в1 °о « г , |
. |
(М6-) |
На рис. 4-14 показана зависимость средней плотности мощности от безразмерной длины канала 6. Видно, что
103
для кривой с постоянной проводимостью а = const макси мальная плотность мощности выбрана для максимальной величины г|Л' для каждой б из кривых рис. 4-13. Для случая переменной проводимости характеристики существенно ухудшаются.
Следует особенно остановиться на анализе зависимости плотности мощности от входного числа М. Для малых температур торможения (7о~2000°К) с ростом Мо резко падает статическая температура потока и уменьшение про водимости оказывается настолько большой, что для около звуковой и сверхзвуковой областей течения плазмы вели
чины N близки к нулю. Максимальная величина N имеет место для Т0 = 2000оК при М= 0,4 (рис. 4-15). Возрастание
Рис. 4-15. Максимальная плотность мощно сти при 7'о=2000°К
удельной мощности до точки максимума объясняется тем, что к. п. д. преобразования цлувеличивается быстрее, чем уменьшается По-
Однако, если входное число М будет возрастать и далее, то удельная мощность достигнет максимума и начнет
уменьшаться из-за резкого уменьшения |
температуры, |
с падением которой по экспоненциальному |
закону (урав^ |
нение 2-5) падает а. |
|
104
Сравнение расчетных кривых для постоянной и пере менной проводимости показывает, что максимальная плот ность мощ ности при с г = const на 404-50 % превыш ает соот ветствующую величину для переменной проводимости.
Оптимальное входное число М возрастает с увеличе нием начальной температуры То. На рис. 4-16 показаны зависимости плотности мощности для постоянной и пере менной проводимости при То= 5000°К. Видно, во-первых, что резко вырос уровень удельной мощности из-за увели
чения |
а (N |
для |
То=5000°К |
составляет 1004-120 |
мвт/м3, |
в то |
время |
как |
W для Го |
=2000°К равно 0,06 |
мвт/м3). |
Во-вторых, оптимальное число Мо, при котором достигается
максимальная величина N, возрастает. При Го= 5000°К МОПт= 0,9. При сверхзвуковых режимах (Мо<11,5) еще имеет место достаточно высокий уровень мощности.
О 0,5 1,0 /,5 2,0 2,5 3,0
Рис. 4-16. Максимальная плотность мощности при Го=5000°К
В-третьих, сравнение расчетных кривых с постоянной и переменной проводимостями показывает, что при высоких температурах допущение о постоянстве проводимости не вносит серьезной ошибки в характеристики преобразо вателя.
На основании проведенного анализа можно сделать вы вод, что для представляющей интерес области температур Го= 20004-5000° К оптимальный МГД-генератор должен быть дозвуковым. Расчет такого преобразователя в первом приближении можно выполнять по уравнениям (4-304-4-34) для постоянной проводимости, причем с ростом темпера
105
туры торможения на входе в преобразователь неточность расчетных характеристик будет уменьшаться.
4-6. Частные случаи движения газа
вэлектромагнитном поле
1.Движение вязкого газа в магнитном поле
Если отсутствует внешнее электрическое поле (£ = 0), то полная энтальпия газа сохраняется постоянной вдоль
канала, так как |
уравнение энергии |
принимает вид |
|||||||||||
|
и |
dx |
|
|
dx |
- |
0 |
или |
i |
— const. |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||
Запишем уравнение количества движения |
|
|
|
||||||||||
|
и JlL + |
±JJL + Ж. + i J fl = |
о. |
|
|||||||||
|
|
dx |
|
р |
dx |
р |
|
2D |
|
|
|
||
Здесь D —гидравлический |
диаметр |
канала, |
где |
р — пери |
|||||||||
|
метр, |
a F — сечение канала; |
|
|
|
|
|||||||
£ — коэффициент сопротивления трения. |
|
||||||||||||
Закон Ома для Е — 0 имеет вид |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
j=ouB. |
|
|
|
|
|
|
Тогда уравнение количества движения |
|
|
|
|
|||||||||
|
|
udu |
, |
1 |
dp |
|
аиВг . |
, |
и1 |
— 0 |
|
||
|
|
dx |
|
р |
dx |
|
р |
|
2D |
|
|
|
|
Из уравнения неразрывности получим |
|
|
|
|
|||||||||
|
|
_1_ dp_ ___/_1_ du_ |
|
d F \ |
|
|
|||||||
|
|
р |
|
dx |
|
|
V u |
dx |
F |
dx / |
|
|
|
J _ _dp_ _ J _ dp_ dp_ _ al _dp_ |
a2 f_\_ _du_ |
J _ dF \ |
|||||||||||
p dx |
p |
dp |
dx |
|
|
dx |
|
\ и |
dx |
|
F dx j |
||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
da |
|
/ |
1 |
da |
j |
1 |
dF \ . |
auB2 |
. t |
a* |
|
||
и --------a' |
V a |
dx |
|
F |
dx ) |
p |
|
2D |
|
||||
dx |
|
|
|
|
Так как to=const, то можно использовать уравнение энер гии в форме
U I |
а* к + 1 |
h — —— h |
= con st, |
к - 1' Т к — 1 |
106
где а. —критическая скорость |
|
|
|
|
|
|
|
|
|||||
|
|
а* |
0 ± 1 |
а2 |
|
к —1 |
|
|
|
|
|||
|
|
|
|
|
2 |
‘ |
|
2 |
|
|
|
|
|
du. |
л : f |
1 |
о |
|
к — |
1 |
Л |
1 |
r fw |
|
£ ± 1 а2 — |
|
|
U-------- |
—— |
а2, ---------- и * ------------ |
|
||||||||||
d.v |
2 |
|
|
|
2 |
/ |
i |
i |
dx |
|
2 |
' |
|
|
к - 1 |
л 1 |
dF |
|
аиВ2 |
Е и’ |
|
|
|
||||
|
2 |
|
) |
F |
<1х |
|
|
р |
|
2D |
|
|
|
Преобразуем левую часть |
1 |
da |
|
н+ 1(и* |
|
|
|
||||||
|
|
|
|
|
|
|
du |
||||||
|
|
|
|
|
U |
dx |
|
|
2 |
|
|
dx |
|
Преобразуем правую часть |
|
|
|
|
|
|
|
|
|||||
/л г + 1 |
2 |
к — 1 |
« \1 |
dF |
|
аиВ1 |
|
^ и2 |
|
|
|||
[ 2 |
а ' |
2 U j ~ F ~ d I |
|
|
р |
|
6 2D ~ |
|
|||||
__ к + I а2 Л |
|
/с — 1 (г3 N 1 dF |
|
оиВ2 |
^ иа |
|
|||||||
2 |
( |
|
л: + 1 а2 ] F dx |
|
р |
|
2D |
|
|||||
Я, = —-----коэффициент скорости |
разделим |
на |
£ ± i |
а2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
(A,3 - |
= |
(1 _ |
Л и 1 А'Л -L i f --------- — |
сцР3 |
, |
||||||||
ра2 |
|
||||||||||||
X dx |
|
у |
|
/с + 1 I F dx |
к + 1 |
|
|||||||
|
|
|
|
+ |
А* |
|
|
|
|
|
|
(4-38) |
|
|
|
|
|
2D |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
Анализ уравнения (4-38)
На основании уравнения (4-38) можно определить, если известно F(x) и принять £ = const.
Два последних члена отрицательны. Следовательно, влияние магнитного поля при Е — 0 качественно такое же, как и влияние сил трения.
Для канала постоянного сечения |
= 0.j уравнение |
(4-38) имеет вид
dX |
ейВ2 |
Хг |
dx |
к + 1 ра2 |
+ |
~2D |
107