Файл: Динамика и управление ядерным ракетным двигателем [Текст] 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