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

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

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

а оптимальное число М, соответствующее максимальному значению функции F\, имеет вид

Мопт

(4-78)

Уравнения (4-77) и (4-78)

дают значения числа М для

любой точки канала генератора с заданными параметрами торможения ро и То, в которой имеет место максимум функции F\, а следовательно, и максимум локальной мощ­ ности, отнесенной на единицу длины (а при постоянной единичной ширине канала и удельной объемной мощ­ ности). На рис. 4-22 даются зависимости оптимального

числа М от Т0 для слабоионизированиой плазмы с добав­ ками цезия и калия при у=1,66. Приводится также кривая для присадки калия и у=1,3.

Из приведенных на рис. 4-22 графиков можно сделать вывод, что максимальная локальная мощность достигается во всех практически интересных случаях в дозвуковой области движения плазмы от чисел М, равных 0,4 до 0,8, причем с увеличением температуры То и уменьшением показателя адиабаты у оптимальное число М увеличи­ вается.

Здесь уместно отметить, что при использовании уравне­ ния для электропроводности плазмы с присадками допус­

122

калось, что между компонентами плазмы равновесие в каждой точке канала достигается мгновенно, т. е. время, затрачиваемое на рекомбинацию при столкновении ионов с электронами на несколько порядков меньше времени пребывания частиц в генераторе.

Приведенный в данном параграфе анализ носит при­ ближенный характер. Кроме того, мы рассмотрели только одну прямую задачу: в канале заданной длины с произ­ вольным заданным изменением сечения найти оптимальные величины чисел М для заданных известных параметров торможения 7’о и роПоказатель адиабаты также считается известным. Однако, при проектировании МГД-преобразо- вателя часто необходимо найти оптимальную форму канала преобразователя, обеспечивающую при заданных начальных параметрах газа наибольшую выработку мощ­ ности. Решение задачи выполняется, как правило, вариацио!':Iыми методами.

3. Изотермическое движение газа в электромагнитном поле

Часто в практике встречаются случаи движения иони­ зированного газа (жидкости) в магнитном поле при посто­

янной температуре. Такой случай имеет

место,

например,

в жидкометаллических МГД-генераторах,

где

электриче­

ская энергия генерируется за счет изменения кинетической энергии жидкости.

Прежде всего запишем уравнения сохранения для изо­ термического одномерного установившегося потока в маг­ нитном поле

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

ри — — 1— — = — ]'В — уравнение количества дви-

dx

dx

 

жения;

(ill

■n

— уравнение энергии;

ри*---- = —

dx

р — р R Т — уравнение состояния;

j = о иВ{ 1 — т(э) —закон Ома.

 

Уравнение энергии можно преобразовать к виду

 

_______Р_______ JLd = — 1.

(4-79)

c £zi], (1-щ ) dx

 

123

Так как температура газа остается постоянной, то прово­ димость в соответствии с уравнением (4-42) зависит толь­ ко от давления

Ф 1 = (PiPi)~* ■

Для слабоиопизпрованного газа 2—0,5. т. е. проводимость равна

о = CTj

(4-80)

Уравнение состояния для Т — const дает

 

Р/Pi = PiPi-

(4-81)

Подставляя проводимость из уравнения (4-80) и плотность

газа из (4-81) в уравнение

энергии

(4-79),

получаем

 

рпп

1_____

/ р \3 2 Л

/ а \

 

(4-82)

°iВг

71,(1— г1й)

[ р, /

dx

(. /а )

 

 

 

Рассмотрим

изотермическое

движение

газа в

канале

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

Как и прежде

считаем, что током

Холла

и трением на стенке генератора можно пренебречь. Кроме того, считаем, что коэффициент нагрузки генератора т)п

постоянен.

В этом случае уравнение

(4-S2) интегрируется

в следующем виде:

 

 

 

.VА',

Q

(4-83)

 

11,(1 —11,)

 

 

 

где Q=.f

и

 

(4-84)

l f T d

 

 

 

и,

и,

 

 

Р Ш 1

(4-85)

aiB"

 

Уравнение (4-85) дает зависимость между длиной генера­ тора х, на которой скорость изменяется от и\ до г/2, и ре­ жимными параметрами. В качестве базовой длины выбрано расстояние х\, полностью определяемое из входных усло­ вий (p\U\ — приведенный массовый расход газа через еди­ ницу площади, щ задано начальным давлением и темпера­ турой, магнитная индукция считается заданной). Тогда длина генератора л: определяется при помощи базовой ве­

124

личины

х\,

помноженной на

два коэффициента:

первый

1

 

-

о

п

-------------

задается нагрузкой

генератора, другой

Ц — па-

( 1

Т|э)

 

 

 

дением давления по длине канала (обычно близок у к 1). На рис. 4-23 представлено изменение базовой длины х\ от проводимости газа ai для различных величин магнит­

ной индукции. Из графика видно, что с ростом проводи­ мости газа и возрастанием магнитной индукции длина Х \ уменьшается.

Коэффициент Q определяется следующим образом. Исключим из уравнений энергии и импульса плотность тока /. Тогда получаем выражение

 

Р - ^ О - Ъ )

^ -

f = 0

(4-86)

или

dx

 

 

и

dx

 

1

dp _

1 — г\э

dit

 

 

(4-87)

 

р

dx

 

Т

dx

 

 

 

Уравнение (4-87) легко интегрируется в виде

 

^ ■

= Н

Ф ^

М

1 2 “, - т Ф ) -

<4-88>

Тогда коэффициент Q находится

 

 

1

ехр |з, 4 (1 -

ты

_

 

 

Q = |

(4-89)

"1 Hi

 

ть

 

k RT

К1

'

 

 

 

 

 

 

125

Введем обозначения

 

а = 3 4 1-1Ь к М 2 и л- = — .

«1

Перепишем затем уравнение (4-89) в виде

 

1

 

Q = [ exp а (.V* — 1) dx.

(4-90)

■г,

 

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

генератора при и = и\.

Случай %2 = 0 практически

не реали­

зуем, так как при

этом площадь поперечного

сечения

канала становится бесконечно большой. Однако вычислен­

ный на основании

Х2 = 0

интеграл дает

максимальное зна­

чение Q, не очень

сильно

отличающееся

от истинного.

Выражение (4-90) может быть представлено

в виде

I

 

Q<Qn,0x= [ exp а (я— 1) (x+l)dx.

(4-91)

о

 

Интеграл (4-91) не представляется в виде элементарных

функций, поэтому рассмотрим предельные значения,

меж­

ду которыми лежит функция Qwax.

 

 

 

1

 

 

1

 

 

 

j" с . х р [а (л- 1) ]

d x >

Qmax >

j*

e x p f 2 u ( . v —

1 )] t/ л ;

 

о

 

 

о

 

 

 

exp a (x — 1)

>

QmiU->

exp 2a (к -

1] ;

(4-92)

1 — exp (— a)

^

1 — exp (— 2a)

 

 

b

^

^

 

2a

 

 

Для малых величин a последнее неравенство может быть заменено на неравенство вида

1 2a>Qmax> \ — а.

126

Такой вывод может быть получен при прямом интегриро­ вании первого уравнения системы (4-92).

Для больших величин a Qmax равно

_

у~

Q m ax = exp ( - a)lYa

j

exp л-1 dx.

о

При а-> QmaX -»■ 0, так что соответствующие коэффи­ циенты нагрузки г) также стремятся к нулю. Остается открытым вопрос о том, к чему стремится произведение

коэффициентов при х\, равное Q

------------- . Коэффициент

 

Ъ (1 —4 а )

нагрузки т)л- молсет быть выражен через величину коэффи­ циента а

3к A4j

Тогда коэффициент перед длиной х\ лежит в пределах

1

Зк Л4j

 

exp (— f t ) >Q— >

4/3/с М\

М

 

 

1 - 4 а

> -----Ц -

fl +

4a

{1 — exp ( - 2 a ) }

З/2/с M\

\

I

Минимальная величинa Q 1—4a имеет место при больших

значениях a (a—»^) и равна —

На первый взгляд

3 7 М г

 

удивительным является результат, что минимальная длина характеризуется нулевым значением, т. е. нулевой напря­ женностью электрического поля Е и нулевой мощностью, генерируемой в единице объема. Однако при стремлении т)э к нулю стремится к бесконечности объем, так что мощ­ ность на единицу объема остается конечной не нулевой величиной, равной в точности изменению кинетической энергии. Здесь уместно отметить, что при малых значениях очень быстро возрастает площадь проходного сечения канала, и одномерное приближение становится несправед­ ливым. Однако минимальной длиной при т] -*-0 удобно пользоваться как первым достаточно грубым приближе-

127

нием, дающим нижнюю границу для длины канала. Следо­

вательно, минимальная длина для г|

 

О равна

 

2

Р

|»Ч

2

(4-93)

 

 

Зт Лф

ЗТ -'И/

 

Уравнение (4-93) для минимальной длины генератора может быть представлено в виде зависимости от началь­ ных параметров торможения (Т0, ро) и числа М|. При ана­ лизе считаем, что разгонное сопло генератора адиабатное, что дает возможность вычислить р\ и Т\ через ро и Т0. Так как проводимость газа может быть выражена в виде степенной функции от температуры и давления, то послед­ нее уравнение легко преобразуется к следзчощему виду:

xmtn^ 2 p l \ 3 S T y UiB ^

= 2/3 Ро32

S T 0B*

= 2/3

p f

 

(р,'р)32

( П / Г о У {ср (Т0 — П ))12

Зт

( T ilT о )

2 (т — 1)

S T ^ 12B2(2cp)12

f +>2\[ T ° ) ljl!

— А л 1 +

М г

T —I AH

= л'0ф(М), (4-94)

где

x0 = 2j3pl*ISTj'+'*B*(2cpy »;

а = у + 1/2-------- —

;

 

*

 

2 (к —1)

 

 

ф (Л4)

: м\' а

Г к 1

м ? 1

1,2

2

 

2

'J

 

Как видно из уравнения (4-94), минимальная длина генератора выражается в виде произведения некоторой функции .то, зависящей от параметров торможения перед генератором (р0, То), на функцию от числа Mi на входе в генератор.

128

На рис. 4-24 представлена зависимость функции хо от Го. для различных давлений торможения ро. С возрастанием начальной температуры длина %o, па которой срабаты­ вается заданная кинетическая энергия, превращаемая в

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

Для заданных начальных параметров ро и То длина

задается числом М|. Продифференцируем

уравнение

(4-94) по числу М|, находим,

что длина х тЫ имеет минимум

при

 

 

= (Т — I) - 1а (се---- * -

(4‘95)

Данные для гр(М) от числа

Mi представлены

на рис. 3

для у= 4/3. Из графика видно, что минимум имеет место

при Ali~0,65. Однако

можно

показать,

что при Mi = 0,65

к. п. д. цикла Карно, равный

 

 

-1

цт = 1 —

Т1

 

т -

- 1 - 1 +

 

М?

будет очень низким (для у= 4/3, г)Y=6,7%). Следователь­ но, для оптимального, с точки зрения генератора, мы полу­ чим очень низкую эффективность. Поэтому следует пере­ ходить на большие длины канала, приемлемые как с точки зрения сжатия мощности, так и к. п. д. генератора.

9 -5 9 9

129

При изотермическом движении газа в магнитном поле в канале постоянного сечения из уравнения количества движения и энергии можно легко вывести соотношение

du

jB

(4-96)

dx

 

1 - ТуИ2

 

Р

 

Выражая и -^ - при помощи уравнения энергии

^dx

 

 

 

 

т_

 

 

 

находим

 

 

Ри

 

 

 

 

 

 

1

 

 

 

 

Е

 

 

 

(4-97)

 

 

иВ

 

ууИа — 1

 

 

 

 

 

 

Подставляя

коэффициент

нагрузки

из

уравнения

(4-97)

в уравнение для изменения числа

М

(уравнение

27а),

получим

 

 

 

 

 

 

f аб2dx — 1 /2 / — ---- 1)

I

( ‘[М] -------— \ - 21п(Л4/Л4,).

L“l

J

\ М2

I

 

у

 

Р 1 « 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4-98)

В случае постоянства произведения сгВ2 имеем

а аВ*

Pi«i

-

f ^ - 1

T/Ms

21n f

(4-99)

Af,

V/и,

 

Для M[ = l число М изменяется по длине канала по сле­ дующему закону:

/VI = ] / — — а .

УТ

 

 

 

Л И ТЕР у\ ТУРА

 

 

 

1.

К у л и к о в с к и ii

А. Г.,

Л ю б и м о в Г.

А.,

Магнитная

гидро­

динамика, Фпзматгнз, 1962.

 

газодинамика

и

динамика

плазмы,

2.

Б a ii Ш и - И, Магнитная

Пзд-во «Мир»,

1964.

Курс

магнитной гидродинамики, Изд-ио

«Мир»,

3.

Дж. Ше рклиф,

1967.

Дж. Ш е р к л и ф,

«Инженерные вопросы магнитной гидродинамики»,

4.

Изд-во «Мир», 1965.

 

 

генерирование

электроэнергии»,

Изд-во

5.

«Магпитогндродииамнческое

«Мир», 1966.

Дж.,

Ше р ма н

А., Основы

технической магнитной

6 . Са т т о н

газодинамики, Изд-во «Мир», 1968.

 

 

 

 

9*

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