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

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

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

где у(х) —высота канала.

Рассмотрим для простоты пря­

моугольный канал постоянной ширины b— const.

основных

Анализ будем производить

при следующих

допущениях:

т. е.

и = и{х), р= р(.т),

р = р ( х) .

1. Поток одномерный,

Рис. 4-11. АУГД-генератор

2. Ток Холла пренебрежимо мал, т. е. индуцируемый в направлении движения ток jx является малым по сравне­ нию с /у.

3. Газ совершенный, т. е. для плазмы можно использо­ вать уравнение Клапейрона p = pRT.

4.Проводимость плазмы а постоянна по каналу.

5.Показатель адиабаты к остается величиной постоян­

ной.

6.Трение и теплообмен на стенках МГД-канала отсут­ ствуют.

7.Магнитное и электрическое поля стационарны.

8.Магнитное число Рейнольдса мало (Re„<$C 1), т. е. можно пренебречь влиянием на поток индуцированного магнитного поля.

Рассмотрим более подробно допущение 8. Уравнение

Максвелла (1-13) дает

rotB =/.

Так как магнитное поле имеет только одну составляющую,

направленн}чо

вдоль

оси

z и

изменяющуюся только по

I

о

---- =

п

то

 

оси л '----

— 0

О

 

и л и

dB

j = — (Ё — uB) a.

(4-18)

 

dx

Приведем уравнение (4-18) к безразмерному виду

dB*

n

...

—r

= - R e ,

■]*,

dX*

 

 

где

j 3“цВц

Безразмерная величина плотности тока /* по абсолют­ ной величине меньше единицы, так как

7

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

1

з£ 0 GWo^o

= 1 ^ - 1 | < 1-

зи0В0

ouoBq

 

dB'

= Re*/*

 

dx*

 

 

dB _ : 0 B = B0=rCOnst.

dx

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

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

ри Е =

 

const —уравнение неразрывности;

da

 

 

 

 

 

 

(4-19)

_

<<P

)В— уравнение

сохранения ко-

ри----

:

dx

 

dx

 

 

 

личества движения;

(4-20)

 

 

 

 

 

S . m [crT + J ^

) = -

 

jE

—уравнение

энергии;

(4-21)

 

1

 

 

 

p — p R T —

к -

 

 

 

к

 

-f'СрТ —уравнение состояния.

(4-22)

 

 

 

 

 

 

 

93

кроме указанных

уравнений

следует

написать

закон

Ома для плотности тока /

 

 

 

 

/ =: о (иВ —Е) = о (иВ----- — )

а [иВ------ -—■') ,

(4-23)

V

I

2y(x)j

\

к (л-) /

 

где h(x) —высота канала МГД-геператора.

Уравнения (4-19) —(4-23) представляют собой систему

уравнений с пятью неизвестными функциями: и, р, р, Т и /. Перепишем систему уравнений, использовав уравнения неразрывности, состояния для идеального газа и закон Ома. Для капала постоянной ширины, равной b= 1 м уравнения магнитной газовой динамики принимают вид

т = puli,

du

Р“ — dx

. dp

+ ~~r 4- a (uB ---- —) = 0, dx

d

m

/

к

p .

u-

+ я | „ в - Д ) Д = о. (4-24)

 

---------— 4-----

dx

 

\H- — 1

p

2

 

 

 

 

 

 

=

RT.

 

 

 

 

 

P

)

Исключая плотность p из уравнений импульса и энергии, при помощи уравнения неразрывности, находим

 

т

du

а ( и В ----—^ В — аВ3и

1—

 

п

dx

 

dx

I

к )

 

киВ

 

 

 

 

= — аВ2и (1 ----—

(4-25)

 

 

 

 

 

V

Цк

 

где

11

=

иВ

 

 

 

 

v

 

 

 

 

 

 

 

 

 

Интегрируем последнее уравнение,

находим

 

ft - Л =°в' f(дг ” 1)

i Т

(4-26>

о

о

 

где а —общая длина генератора.

94

Если задан закон

изменения скорости и{х), можно

найти распределение

давления по длине генератора, так

как при постоянном магнитном поле В и постоянном э. д. с.

V есть

функция только и(х). Условие

т] =

const соблю­

дается

точно

только

для

активного генератора,

для

кото-

рого скорость и по длине канала постоянна,

т. е.

clu

п

-----=

U.

Для

активного МГД-генератора

давление

 

с/х

 

находится как

Р ~

Ра =“в‘ I (■ Т|/(

 

11udx =

з Вги

с/х

 

 

оВ’иа. (4-27)

 

1]

 

h

 

 

Преобраз)'ем уравнение энергии

 

 

 

 

 

 

 

 

к

d

, _ , ч ,

d

 

 

 

 

 

 

 

 

 

 

dx "

dx \ 2 j

 

 

(

 

 

ц/г j t]

 

Последнее уравнение легко интегрируется

 

 

 

 

 

 

 

J i -

puh +

 

= а В2 f

± -

- \

\

j

-

L d

x -

 

 

 

 

 

 

T\h

 

 

 

Ti

 

 

 

 

 

 

-Рои»К +

 

 

 

 

 

 

 

(4-28)

где po,

uq —давление

и

скорость

на

входе

в

МГД-гене-

 

ратор;

 

 

 

 

 

 

 

 

 

 

ho — высота входного сечения в МГД-генератор.

Из уравнения (4-28) находим давление

к — 1

- ')■ г 1л + ДЬ

+

кик сВ'

+ Y

(“I - “■)].

(4-23)

Решаем совместно уравнение импульсов (4-27) и урав­ нение энергии (4-29) и находим искомое уравнение для

95

определения профиля МГД-генератора, обеспечивающего заданное распределение и { х ) .

к -

1

{) dx + 7 ~ i Pu"o/'“ ь Т ("» - Mi)

ки

I (тГ -

h[x) — -----

О

 

 

h

а

 

 

dx_

 

 

dtt

 

 

itdx -|- гп I

h

 

0

6 dx

(4-;ю)

Трансцендентное уравнение (4-30) представляет собой общее выражение для профиля канала МГДГ при произ­

вольном

изменении

и(х) и решается

методом

последова­

тельных

приближений.

 

 

Для

активного

МГД-генератора

(z^ = const,

r]= consl)

уравнение (4-30) решается в квадратурах. Действительно,

уравнения импульса и энергии

для и — const

имеют вид

 

 

 

«P- = i L L - auB\

(4-31)

 

 

 

dx

h

 

 

и

к — 1

-±- (ph) = ----а-~ (uBh - V).

(4-32)

 

 

dx

 

п

 

Решая уравнения (4-31)

и

(4-32) совместно, находим

 

а =

крп

 

h — lla

 

 

VBs

 

ц1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к - 1

 

 

 

 

 

(1

—-r/0 ( I — —

ц/;0 )

1—

*

 

цА.) In ------------ *----------!----- . (4-33)

 

К—1

(1

- v g ( \

11/!)

 

 

 

 

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

крп

кр0В _

nit]

_ А

VBs

VB

аВ* Afjj

Ml

96

Где А —величина постоянная для заданных а, В, г| и равна

пщ/оВ2.

После нахождения профиля канала h(x) находим изме­ нение плотности, температуры и числа М по длине канала МГДГ следующим образом:

плотность

Р (*) =

ГП

и (х) h (х)

температура

 

Т (х) _ Р(х) .

число М

Rp (х) '

 

и (X)

М( х ) = -

VKRT (х)

Впараграфе 4-3 было показано, что плотность генери­

руемой мощности (мощность,

снимаемая с единицы

объе­

ма генератора) равна

 

 

 

 

N =

аи2В2

“ Лэ) Лз,

 

И Л И

 

 

 

 

N —— | ou*B2h (х) 1 (1 — щ) dx.

(4-34)

Представляет интерес также средняя плотность мощно­

сти, снимаемой со всего объема генератора

 

— _

yV

тИдг

 

1

V —

 

V

 

Полная энтальпия на входе в канал равна

 

7(1 +

М' \ =

- ^ — R

М*}.

i o - i + — - c p

 

 

 

 

Так как

 

 

 

 

TJ T =

1 + - — - М \

 

то

к -

177’o'rljv

 

N =

(4-35)

 

 

 

7—599

97

Подставляя уравнение (4-30) для профиля МГД-генератора в последнее соотношение находим мощность генератора при заданной его длине или длину а при заданной мощ­

ности

N.

 

(w=const) электрическая

энергия

В

активном канале

генерируется

только за

счет

температуры

потока, так как

в выражении для энтальпии торможения

t'o= ^срТ + — 'j

кинетическая

энергия

газа

не меняется

по всей

длине

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

N = гпср -Г°- С Ъ.. .

(4-36)

То

 

Мощность может быть выражена также через давление на входе и выходе из МГД-генератора

N = тСрТо

(4-37)

После определения параметров потока и мощности, находим показатели эффективности МГД-генёраторов: внутренний относительный коэффициент полезного дей­ ствия rjoi и коэффициент полезного преобразования мощ­ ности rpv по уравнениям (4-38), (4-39).

Эффективность генератора можно наглядно проследить на примере МГД-преобразователя постоянного сечения ( h = const) при постоянном магнитном поле B = 5o= const и постоянной проводимости плазмы о= const. Для общно­ сти численный расчет выполнялся для безразмерных вели­ чин. В качестве безразмерных переменных использовались:

длина

5qЛ 1а

о — ----------

 

m

число

М М — Аи ;

 

V kRT

электрический к. п. д.

внутренний

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

к. п. д.

 

коэффициент полезного преобразования энергии

V .

uBh

Too Т\

Ло; — То.~Т\

N

mi„о

98

На рис. 4-12 показаны характеристики МГД-преобразо- вателя постоянного сечения, полученные в предположении постоянной проводимости плазмы, для дозвуковых и сверх­ звуковых потоков на входе в канал. Внутренний относи­ тельный к. п. д. г|о; и к. п. д. преобразования T]lV показаны как функции электрического к. п. д. Безразмерная длина канала задавалась: 6—0,3; 0,6; 1; 2; т)э изменялась от 0 (закороченная цепь) до 1 (разомкнутая цепь).

Анализ рис. 4-12,а для дозвукового режима на входе в канал М= 0,3 показывает, что при малых величинах 6 по мере уменьшения г)э от 1 до 0,5 т]0( падает, а к. п. д. преоб­ разования возрастает. При г|э, равном приблизительно 0,5 rjjV, достигает максимума, и при дальнейшем уменьше­

нии цэ

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

цо.- и цл-

вплоть до

нуля при

г)о —0.

В этом режиме канал

работает

в режиме

короткого

замыкания, и кинетическая энергия газа превращается в джоулево тепло. Для больших длин 6= 1 к. п. д. преобра­ зования ограничен сверху условием запирания, при кото­ ром поток становится сверхзвуковым на выходе из канала. В этом случае в канале или за его срезом возникает удар­ ная волна. Следовательно, для коротких каналов 6=1 электрический к. п. д. (параметр нагрузки) выбирается из условия минимальных джоулевых потерь, в то время как для более длинных каналов максимальный к. п. д. преобра­ зования определяется условием запирания.

Увеличение длины МГД-генератора при постоянной нагрузке (i]0 = const) приводит к увеличению к. п. д. пре­ образования i].v, причем преобразование энергии происхо­ дит с меньшей эффективностью (см. верхний рис. 4-12,а). Скорость потока в дозвуковом преобразователе возрастает

по длине

канала,

следовательно, возрастает

индуцируемая

э. д.с. иВ

и для

обеспечения постоянства

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

к. п. д. (постоянного параметра нагрузки) необходимо уве­ личить внутренний перепад напряжения.

В сверхзвуковом преобразователе (рис. 4-12,6) с увели­ чением длины генератора при Т]э= const возрастает как к. п. д. преобразования Г|л-, так и внутренний относитель­ ный коэффициент полезного действия МГД-генератора, так как в канале с увеличением его длины скорость потока падает, что вызывает уменьшение индуцируемой э. д. с. Еинд— 0- Следовательно, для соблюдения условия r)3=const необходимо меньшее падение внутреннего напряжения. Таким образом, в отличие от дозвукового генератора с

7*

99

Х.п.д. преоИраоВания, % Изотропичеагии клд., %

а)

too

80

60

40

20

О

Рис. 4-12. Дозвуковые (а) и сверхзвуковые (б) рабочие характеристики МГД-преобразователя с постоянным сечением канала и s=const

ростом длины генератора возрастает к. п. д. преобразова­ теля т]л-, причем преобразование энергии происходит с большей эффективностью. Условия запирания, т. е. усло­

вия, при которых число М на

выходе равно 1, показаны

на рис. (4-12) цифрой 1.

что максимум к. п. д. преоб­

Из анализа рис. 4-12 видно,

разования энергии как для дозвукового, так и сверхзвуко­ вого генератора имеет место при цэ приблизительно рав­ ном 0,70-т-О,5.

Для дозвукового канала (М = 0,3) газ можно принять несжимаемым (т. е. как твердый проводник). В этом слу­ чае максимум энергии имеет место при ц0 —0,5, т. е. когда энергия, подводимая к внешнему источнику, равна энергии, теряемой на омический нагрев газа в канале генератора. Для больших длин дозвукового генератора не удается достичь условия максимума г]л-, так как раньше наступает условие запирания. Кроме того, с возрастанием длины генератора возрастает число М на выходе из него и начи­ нает сказываться эффект сжимаемости газа, оптимальные условия при этом получаются более сложными при боль­ ших М.

В сверхзвуковом генераторе вследствие сжимаемости газов оптимальные условия для к. п. д. преобразования отличаются от условий для твердого проводника, и вели­ чины т}э0ПТ лежат в пределах 0,5-1-0,7.

4-5. Течение в канале МГД-генератора с переменной проводимостью газа

Электропроводность ионизированного газа зависит в первую очередь от его местных термодинамических пара­ метров (давления и температуры), которые в МГД-преоб- разователе достаточно быстро уменьшаются. Будем рас­ сматривать течение плазмы в канале МГД-генератора при тех же допущениях, что и в предыдущем параграфе за

исключением

допущения 4 о постоянной проводимости,

для которой

используем уравнение Чэпмена — Каулинга

(2-5).

Основные уравнения магнитной газовой динамики не изменяются и имеют такой же вид, что и уравнения (4-24). Однако в этих уравнениях величина а принималась зави­ сящей от местных величин давления и температуры по

формуле (2-5).

Эти

уравнения были запрограммированы

для числового

счета

на электронно-вычислительной маши­

101

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