Файл: Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие.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