книги из ГПНТБ / Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие
.pdfгде у(х) —высота канала. |
Рассмотрим для простоты пря |
||
моугольный канал постоянной ширины 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 |
VB2а |
аВ* 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