Файл: Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 05.08.2024
Просмотров: 70
Скачиваний: 0
нейтральных молекул (атомов) находим через статическое давление и температуру газа
— |
зр — |
р |
(2-15) |
' — |
т С т ~~ |
kT |
|
Используя численные значения заряда электрона, его мас сы и постоянной Больцмана, получаем следующее выраже ние для параметра Холла:
= 3,6-10-1G |
В Г12 |
(2-16) |
|
PQ |
|
Из уравнения (2-16) видно, что наибольшее влияние на р оказывает статическое давление Р и магнитная индукция В.
В случае, когда плазма содержит смесь рабочего веще ства и ионизирующие добавки, параметр Холла будем находить из уравнения
1
13,1-10° т 2 |
в |
(2-17) |
Р .= |
|
|
«.Qo + nsQs |
|
|
где индексы 0 и s относятся к молекулам рабочего газа и ионизирующих присадок.
При больших значениях параметра Холла р„>10 необ ходимо учитывать вклад ионной проводимости, вносимой в электропроводность газа. Уравнение для результирующей проводимости будет иметь вид
|
|
П: е* X; |
<* = |
+ ° е |
+ |
Плотность тока в направлении оси у будет равна
h |
— - --------h |
— ----1Е, |
|
или |
i + |
в? |
1+ (И-/ В)2 \ |
|
|
|
|
|
|
-L |
се |
7у = |
1 + |
Не |
|
В? I |
1+ (щ BY |
Подвижность ионов ц,- (существенно меньше подвижности электронов щ, т. е. щ/щ J> 1 и ( ) 2 1. Тогда прибли женное выражение для проводимости сг„ имеет вид
(2-18)
1+ Й
40
где
>Р« = М?,
■Pi=(Xf-S.
Уравнение (2-18) для сту отличается от соответствующего выражения (2-13) членом (1 + рс ■р.), который учитывает ионную проводимость при достаточно больших величинах параметра Холла (рс> 10).
Г л а в а 3
ОСНОВНЫЕ УРАВНЕНИЯ МАГНИТНОЙ ГАЗОДИНАМИКИ
Проводящие газы (плазма) представляет собой смесь положительно и отрицательно заряженных частиц и нейт
ральных атомов (молекул). В макроскопической теории плазмы для каждой t-ой компоненты необходимо знать температуру 7\, давление Pi, плотность р( и скорость Ci,
а также напряженности Е и Н. Если рассматривать плазму как смесь п компонент, то для определения параметров плазмы необходимо иметь (бтг+ б) уравнений. Этими соот ношениями являются уравнения сохранения: б уравнений неразрывности, выражающие закон сохранения массы каж дой компоненты, уравнения движения, уравнения энергии и состояния для каждой компоненты.
Газодинамические уравнения дополняются шестью урав нениями Максвелла, описывающими электромагнитное по ле. Следовательно, -система уравнений, описывающая дви жение плазмы в электромагнитном поле, является очень сложной. Даже в простейшем случае плазмы без присадок и, когда ионной проводимостью можно пренебречь ( п = 2), полная система состоит из 18 уравнений. Анализ такой системы, а тем более точные решения чрезвычайно трудо емки.
В магнитной газодинамике плазму рассматривают как единый газ неизменного фиксированного состава. Тогда вместо анализа смеси различных компонент можно рас сматривать газодинамические уравнения для суммарных величин .совместно с уравнениями электромагнитного поля. Связь между суммарными и парциальными величинами определяется следующим образом.
42
1. Плотность плазмы р. Определяется как сумма плотно стей отдельных компонент
р= 2 рс
2.Давление Р представляет собой сумму всех парциаль ных давлений Р,- компонент плазмы, т. е.
Р= 2 Р<-
3.Температура определяется как
T = i £ v , r „
где V,- и v — число частиц Той компоненты и сумма частиц
вединице объема.
4.Скорость течения для плазмы в целом определяют осреднением уравнения сохранения массы для всех компо нент, т. е.
С— 1 £ р , г ,
5.Плотность электрического заряда в соответствии с уравнением (1-22) определяется как
Р* :: |
Y, |
= <4 р«. |
|
л V |
|||
|
i |
6. Плотность электрического тока j находится по урав нению (1-23)
2 9/С,
/
Система основных уравнений магнитной газовой дина мики для суммарных величин, характеризующих плазму в целом, будет состоять из следующих уравнений, рассмот ренных ниже.
3-1. Уравнение неразрывности
Выделим в движущейся плазме элементарный объем в форме параллелепипеда (рис. 3-1) и запишем условие
.43
неизменяемости массы во времени для этого элемента. Это условие будет выражать закон сохранения массы; оно может быть представлено в таком виде:
rf(PAK) |
(3-1) |
|
dt |
||
|
где Д V — объем элемента;
р —средняя плотность элемента.
Дифференцируем, имея в виду, что р и Д У —перемен ные величины:
д У - ^ - + Р at
d а V
dt
Разделим это уравнение на рД V. Получим уравнение неразрывности в виде:
|
dp . |
1 |
d д V |
(3-2) |
|
Р |
dt |
д V |
dt |
||
|
|||||
Здесь производная |
d&V |
выражает скорость изменения |
|||
------ |
dt
объема или, следовательно, .скорость объемной деформации
„ |
тг |
1 d&V |
_ „ |
жидкой |
частицы. Член — — |
представляет собой ско |
рость относительной объемной деформации.
В частном случае несжимаемой жидкости, когда p=const, последнее уравнение принимает весьма простую форму:
d л V |
= О, |
dt |
|
выражающую условие постоянства объема элемента: ско рость объемной деформации элемента равна нулю. Отсюда следует, что частица несжимаемой жидкости деформируется в процессе движения так, что объем ее сохраняется неизменным.
Определим величину скорости относительной объемной деформации частицы, выразив ее через соответствующие проекции скорости и, v и до.
Подсчитаем линейную деформацию частицы в направ лении оси я: (рис. 3-1). Эта деформация будет происходить в связи с тем, что скорости граней ABCD и A'B'C'D' неоди наковы. Если скорость левой грани равна и, то скорость
правой будет и+ ——dx.
44
Предположим, что в пределах каждой из рассматривае мых граней параллелепипеда скорости постоянны. За эле мент времени левая грань ABCD переместится на расстоя ние udt вправо. За тот же отрезок времени грань A'B'C'D'
Рис. 3-1. К выводу уравнений Эйлера в декартовых координатах
переместится в том же направлении «а расстояние
(и-1----dx)dt. Следовательно, объем элемента изменится,
дх
так как скорости этих |
двух граней различны. Подсчитав |
|
абсолютное изменение |
объема частицы |
по направлению |
оси ж, получим |
|
|
^и + -rj— d.x'j dy dz dt — и dy dz dt = |
dx dy dt. |
Рассуждая аналогично, для других двух пар граней мож но получить приращения объема частицы по осям у и z в (следующем виде:
dx dy dz dt:
ду
dz dx-dy dz dt.
45
Полное изменение объема частицы определяется как сумма этих приращений.
Следовательно, скорость относительной объемной дефор
мации определяется весьма просто: |
|
||
1 |
d д V |
_ да , dv . dw |
n> |
д V |
7t |
~ ~7х + ~ду + ~дГ ' |
|
так как объем элемента Д У — dxdydz.
Подставив (3-3) в уравнение неразрывности (3-2),
получим |
|
+ f\-r- + |
|
+ -Т-) = °- |
<3-4> |
||
- |
|
- г - |
|||||
р |
dt |
|
од: |
ду |
dz |
J |
|
Заметим здесь, |
что входящие в уравнение (3-3) прямые |
||||||
частные производные |
да |
dv |
дш |
|
|
„ |
|
— |
dy |
— имеют |
конкретный меха- |
||||
|
|
дх |
дг |
|
|
|
«ический смысл. Из предыдущих рассуждений очевидно, что эти производные определяют величины скоростей отно сительных линейных деформаций граней параллелепипеда. Рассмотрим, например, линейную деформацию грани DCC'D'
в направлении оси х. Так как |
скорость ребра CD равна и, |
|||
а ребра CD' и + -—dx, то удлинение грани по х |
будет |
|||
(и |
dx] dt — tidt — |
dx dt. |
|
|
V |
дх ) |
|
дх |
|
Относительное |
удлинение |
|
да |
а скорость |
составляет — dt, |
ди
относительного удлинения — .
дх
Преобразуем теперь уравнение (3-4). Так как р= р(л:, у , z) , то полная производная плотности равна
dp _ dp dx , dp dtj , dp dz
|
dt |
dx |
dt |
dy |
dt |
|
dz |
dt |
|
Имея в виду, что |
dx |
|
dy |
|
dz |
=w, получим |
|||
— = щ |
— = v; |
— |
|||||||
|
|
|
dt |
|
dt |
|
dt |
|
|
|
|
|
, |
dp |
, |
dp |
dp |
|
|
|
dt |
|
+ |
v —t— -f |
w —t— -j-----’t- |
|
|||
|
dt |
ду |
|
dz |
dt |
|
|||
Подставив |
dp |
в уравнение |
(3-4) и |
преобразовав, |
будем |
||||
— |
|||||||||
иметь |
д? _1 <Нрц) _j_ d (ро) |
|
d (ри>) |
|
|
||||
|
+ |
= 0. |
(3-5) |
||||||
|
dt |
|
дх |
ду |
dz |
|
|
46