Файл: Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 05.08.2024

Просмотров: 67

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Умножим все члены на 2D/X2

dX

1 ( 2auB*D

Xя dx

+

к + 1 ^ рп; <х"

где .v =

О ’

так как X

1

 

/ 20 а В~ +

 

(4-39)

к + 1

ри

 

 

 

>сг со

постоянна,

 

 

 

тате

получаем

 

 

X(A.j) — Х(Х\

2

 

( 2Da В2

ё -v,

(4-40)

к -i- 1

Р

 

 

 

где

Х{ Х ) = ± - \ п Х * - ,

?Ч —скорость на входе в канал.

Уравнение (4-40) дает возможность рассчитать распре­ деление скоростей вдоль канала постоянного течения в магнитном поле.

Итак, на поток в канале постоянного сечения действу­ ют силы трения и магнитное поле, причем магнитное поле воздействует на поток так же, как и силы трения.

В дозвуковом потоке (Х<1, X2—1<0) ---- >0, т. е. dx

внешнее магнитное поле в канале МГД (г|>1) ускоряет дозвуковое течение.

В сверхзвуковом потоке (А>1, X2—1>0) и ■----<0, т. е. dx

сверхзвуковой поток тормозится под действием магнитной силы и трения.

2. МГД-генератор с постоянным числом М

Рассмотрим МГД-преобразователь, состоящий из каме­ ры сгорания, сопла для ускорения газа и канала для гене­ рирования электрической энергии. Для исключения эффек­ та Холла электроды выполнены секционированными, т. е. каждая пара противоположных электродов имеет собствен­ ную нагрузку.

108


Исследуем оптимальные условия работы МГД-генера- тора, т. е. такие условия, при которых в заданном объеме генерируется максимальная мощность. Если рассматривать канал постоянной ширины, то это равносильно генерации максимальной мощности на заданной длине канала. Опти­ мальные условия, очевидно, можно сформулировать и другим образом: генератор называется оптимальным, если заданная мощность генерируется в минимальном объеме (а при постоянной ширине канала на минимальной его длине).

Оптимальным условиям соответствуют минимальные потери, так как больший объем канала требует больших затрат на магниты для создания магнитного поля. Тепло­ отвод от стенок МГД-канала тем больше, чем больше сум­ марная поверхность его стенок, окружающих данный объем. Кроме того, в оптимальном канале гидродинамиче­ ские потери в пограничном слое будут минимальны.

Задача об оптимальном МГД-генераторе может быть решена двумя путями: точным вариационным методом и приближенным.

В вариационном методе мощность, снимаемая с данного

объема, записывается в виде функционала

и нахо­

дится максимальное значение функционала методом мно­ жителей Лагранжа. Экстремум ищется как по числу М, так и по параметру нагрузки. Вариационная задача в настоящее время разрабатывается и не может считаться законченной.

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

Анализ проводится при следующих основных допу­ щениях:

1. Показатель адиабаты у по длине канала меняется незначительно, т. е. у const.

2.Квазиодномерное приближение магнитной газовой динамики считаем справедливым даже при сравнительно больших изменениях площади F/F\.

3.Принимаем степенной закон для проводимости а (ниже будет показано, что результаты существенно изме­

109

нятся, если проводимость а определять по уравнению Чэпмена —Каулинга).

Оптимальное число М не зависит от нагрузки генера­ тора, т. е. принимаем постоянным.

Оптимизацию можно проводить по дифференциалу электрической мощности dN0, генерируемой на элементар­ ной длине канала dx. В этом случае мы найдем оптималь­ ные условия в каждой точке канала.

Запишем электрическую мощность dN3 как

dN3 — ап2В2г) э (1 —цэ) Fdx,

или электрическая мощность, генерируемая в единице дли­ ны в любой точке канала равна

^

dx

=

(4-41)

 

 

 

Используя местные параметры торможения, выразим пра­ вую часть уравнения (4-41) через число М. Из газодина­ мики известно, что

^ = 1 + -I— - М2.

(4-42)

Т2

Электропроводимость термически ионизированного газа является функцией температуры и давления. В рассматри­ ваемой задаче удобно использовать простое приближение для проводимости в виде степенного закона

a = ( S T ) vp :,

где 5 —безразмерная постоянная.

У— 1d(\na)/d(\nT) )

z = ( d(lncr)/d(lnp)

Экспериментальные значения для аргона у, z и S пока­ заны в таблице.

Из приведенной таблицы видно, что показатели г, у в

пределах каждой области не изменяются, а

коэффициент

5 при изменении ст от 0,1 до

100 остается

постоянным.

Следовательно, проводимость

может быть

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

в виде

 

 

(4-43)

j- = ( £ r ( t ) " = ( £ r

110


 

 

 

 

 

 

 

 

 

 

 

Таблица 2

 

 

 

 

 

Проводи­

 

У

Z

Нй 2)

 

 

 

 

мость, мо/м

 

 

 

 

 

 

Слабоионизпрованиая

 

10

ОД

 

13

0,5

 

1/2000

плазма

 

 

 

 

 

Частично

ионизиро­

100

 

1

10

0,4

 

1/2000

ванная плазма

 

 

 

Тогда

мощность

на

единицу длины

может быть

записана

в виде

 

 

 

 

 

 

 

 

 

 

 

=

а0(ОСрТ 0)У2^-

Я Ч (1 -

Л.) (T/ToY (1 -

 

T!T0)V\

 

dx

 

И

ро

 

 

 

 

 

 

(4-44)

 

 

 

 

 

 

 

 

 

 

 

где

X=

(у + 1) -

(г +

1) —^

;

 

 

 

 

 

 

пг — расход газа.

(4-42)

в

(4-44),

получаем

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

 

 

СТ0 (2срТоу* f L В Ч (1

 

 

r - ^ i y %

 

=

-

1Ъ) -- ->

2,- ~

 

-----г •

rfx

 

 

Ро

 

 

 

(1+ х

^ 1 ^ + 4

П

'

 

 

 

 

 

 

 

(4-45)

Найдем оптимальное число М, при котором правая часть уравнения имеет максимум. Для этого проведем логариф­ мическое дифференцирование уравнения (4-45), в резуль­ тате которого находим

Monm = {(у + 2) (к -

1) -

(z + 1) Т}-1/2 •

(4-46)

Уравнение

(4-46) показывает,

что

оптимальное

число М

не зависит

от местных параметров торможения

потока,

а определяется, главным образом, физическими свойствами газа. Для постоянного заданного показателя адиабаты у число М0пт, очевидно, остается постоянным. Следовательно, канал генератора можно считать оптимальным, если в каж­ дой его точке выдерживать постоянное число М, опреде­ ляемое в соответствии с уравнением (4-46).

На рис. 4-17 представлена зависимость оптимального числа от показателя адиабаты у для типичных значений

111



t/= 13 и 2= 0,5. Анализ рис. 4-17 показывает, что для боль­ шинства используемых газов оптимальное число М оказы­ вается меньше единицы, т. е. поток дозвуковой (так для аргоновой плазмы у=1,25 и М Оп т = 0,73).

Итак, проведенная оптимизация работы МГД-генера- тора, показывает, что оптимальный генератор имеет посто-

Рис. 4-17. Зависимость /М0пт от показа­ теля адиабаты f

янное число М по длине канала. Исследуем такой гене­ ратор.

Основные уравнения одномерной магнитной газодина­ мики для этого случая имеют вид

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

du

, dp

 

(4-47)

— /В —уравнение сохранения ко­

Р«—

+ - г -

ах

dx

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

(4-48)

 

 

di0

d I

^

и2

ри

= ри

СрТ

4------ -

dx

d x \

 

2

^—]Е— уравнение сохранения энер­

гии;

(4-49)

112

гг, T — 1

 

p — pcp 1 --------- уравнение состояния;

 

/ = a \uB E] закон Ома.

(4-50)

(4-51)

Для случая слабоионизированной плазмы, когда прово­ димость ст зависит только от давления и температуры, най­ дем изменение параметров (р, р, Т, и) и зависимость для площади проходного сечения F в генераторе с постоянным числом М.

Исключим плотность тока j из уравнений (4-47) и (4-48)

- f

(СрТ) + (1 — Yl)'T' (“Т"')---- -- ~7~ — 0

(4'52)

dx

dx \ 2 /

р dx

 

Используя соотношение для числа М

V LT i c ' 7'

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

 

1

dT _

1

dp

 

(4-53)

 

Т

dx

Р/;

dx '

 

где

 

 

 

 

 

 

 

 

Р = (1 +

V (1 - ч){т- ' ) М1

T - - V

(4-53')

I

2

 

 

J

(Т—1)>Ь

 

Величина Р для постоянных М и у также постоянна. Тогда после интегрирования уравнения (4-53) для Р = const получаем

p /p ,= (r/Ti)3.

(4-54)

Для рассматриваемого случая | 1 —| представляет собой

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

dN,

= — р и Ad [с Т ,

-1- и2) j dx —

dx

 

\

"

2

) I

= - pu A d

срТ {l +

- j {к — 1) М1} j dx.

8 - 599

113