Файл: Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие.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