Файл: Динамика и управление ядерным ракетным двигателем [Текст] 1974. - 253 с.pdf

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

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

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

Добавлен: 23.06.2024

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

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

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

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

Наиболее общая форма закона Ома, приемлемая для исследования динамики течения в реакторе, следующая

[18]:

j = а [Е +

Етэ + (1/с) V X В] + р е V - (1/с) RH(j X В).

(5.10)

Здесь а

— проводимость; R h — постоянная

Холла;

Етэ — термо-э. д. с.

Для замыкания общей системы уравнений в том случае, когда электрические токи, протекающие в газе, толь ве­ лики, что наведенное ими магнитное поле существенно иска­ жает магнитное поле, создаваемое внешней магнитной си­ стемой, необходимо в эту систему включить уравнения Мак­ свелла. В магнитогидродинамическом приближении без учета токов смещения эта система имеет вид [17—20, 23]

 

С

ли)

 

 

rot Е =

( — 1/с) (дВ/д/);

(5.12)

diVD = 4яре;

(5.13)

 

div В = 0;

(5.14)

 

D = еЕ;

(5.15)

 

В = pH.

(5.16)

Обычно для газов в

системе единиц

Гаусса е = 1 и

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

187


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

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

Qß = E . j .

(5-17)

Этот член должен быть добавлен к уравнению баланса энер­ гии.

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

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

(др/ді) + div (pV) = 0.

(5.18)

Во-вторых, тепловое излучение обладает собственным тензором напряжений [22] и в изотропном и равновесном случае оптически толстой среды создает дополнительное давление излучения:

PR = {\ß )a RT \

(5.19)

где aR = 7,67 • ІО-15 эргІ(см3 ■град) — постоянная

Сте­

фана — Больцмана. Это давление (а в общем случае тензор напряжений излучения) должно быть добавлено к урав­ нению сохранения количества движения. Однако эти допол­ нительные члены существенны только при очень высоких температурах и в условиях газового реактора существенной роли не играют.

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

новесном случае оптически толстой среды

определяется

по формуле

 

Er/p = aR Tl!p,

(5.20)

188


где р — плотность газа; второй член связан с потоками из­ лучения. Опять же в случае равновесного излучения в оп­ тически толстой среде

 

 

Чк = о Т Ѵ р,

(5.21)

где

а = caRI4

= 5,75 • ІО" 5

эрг!(см1 ■сек град);

о

так­

же

называется

постоянной

Стефана — Больцмана

(с —

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

L * v = l/(p /Q ,

(5.22)

возможно то или иное описание процесса переноса энергии излучения. Здесь Кѵ — приведенный коэффициент погло­ щения для частоты излучения ѵ; р — плотность газа [2 2]. В случае оптически толстой среды для наиболее энерго­ несущей части спектра излучения возможно диффузионное описание лучистого переноса энергии. Когда средняя по спектру длина свободного пробега излучения мала по сравнению с размерами области, в которой оно происходит:

L * = l/(p K * )« L ,

(5.23)

можно ввести коэффициент диффузии

излучения [22]

DR = LRc/3 = c/(3KRp)

(5.24)

и представить радиационный тепловой поток в виде

qr ^ - D r VEr = - X r VT.

(5.25)

Здесь XR — эффективный коэффициент лучистой теплопро­ водности

XR = 4DRaRT*.

(5.26)

В этом случае член, описывающий в уравнении баланса энергии лучистый теплообмен, имеет такой же вид, как и член, описывающий обычную теплопроводность:

Q 7 = d i v ( ^ g r a d T ) ;

QR = div(XRg r a d Т).

189


В данном случае можно ввести понятие эффективного теп­ лового потока и эффективного коэффициента теплопровод­ ности:

<2Эф = Qt -f Qr = diV (V + KR) grad T = div (Хэф grad T). (5.27)

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

=

І Я я у Ь Ч ? г>

(5.28)

где

qRy —J /mdco; qRz = ^Indm.

 

qRx--=\ I Ійщ

(5.29)

Здесь I — интегральная интенсивность излучения по всем частотам спектра; I, ш, п — направляющие косинусы ли­ нии луча, а интегрирование ведется по всевозможным те­ лесным углам для всех направлений излучения. Интеграль­ ная интенсивность излучения получается интегрированием по частотам спектральной интенсивности / ѵ:

оо

 

/ = j' / ѵ dv.

(5.30)

о

 

Последняя же подчиняется специальному уравнению балан­ са излучения:

(1/с) (dlv/dt) -j~ (dly/ds) =

p/v — pkv / ѵ,

(5.31)

где / ѵ — коэффициент излучения;

kv — коэффициент по­

глощения; ds — расстояние между двумя точками по лучу. В большинстве газодинамических задач не рассматриваются

явления, протекающие с

чрезмерно

большой частотой.

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

(1 /c)(dlv/dt)

можно пренебречь.

Тогда уравнение переноса излучения в направлении еди­ ничного радиуса-вектора имеет вид

п .Ѵ/ѵ =

р/ѵ —р/гѵ/ѵ-

(5.32)

Коэффициент излучения /ѵ

в общем случае

определяется

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

190