Файл: Дейч М.Е. Элементы магнитной газодинамики конспект лекций учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 05.08.2024
Просмотров: 72
Скачиваний: 0
Тогда уравнения импульсов (3-8), (3-9), (3-10) можно за писать в такой форме:
|
|
|
|
dCr |
(-Ъ |
R --------- |
dP |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
dt |
|
|
г |
dr |
|
|
|
|
||||
|
|
|
|
|
|
|
Р |
|
|
|
|
|||||
|
|
|
|
dC$ |
|
, |
Сг Сц |
н |
1 |
|
dP |
|
|
|
|
|
|
|
|
|
dt |
|
1 |
г |
и |
|
r d (1 |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
dw _ |
- ± ^ |
- |
|
|
|
|
|
|
||
|
|
|
|
|
|
dt |
Z |
|
’ |
|
|
|
|
|||
|
|
|
|
|
|
|
р |
dz |
|
|
|
|
|
|||
где R, 0 |
и Z — проекции |
единичной |
массовой |
силы на |
оси |
|||||||||||
Подставив |
координат /', 0 и z. |
|
|
полных |
производных |
|||||||||||
в |
|
(3-12) |
значения |
|
||||||||||||
dCr |
dCa |
|
dw |
через |
частные, |
окончательно находим |
|
|||||||||
----, |
---- |
и ---- |
|
|||||||||||||
tit |
dt |
|
ett |
|
|
|
|
|
|
|
|
|
|
|
|
|
dSz. |
|
C |
-1- |
|
r |
<!£l- 4- Щ, ?£l- _ |
|
— f t ___L — |
; |
|||||||
|
dt |
|
r dr |
|
|
dti |
|
dz |
|
|
r |
|
p |
dr |
|
|
dCg |
^ |
dCq |
| |
Cq dC{j |
| ^ |
dCQ |
^ |
Cr C0 __л |
|
1 |
dP |
|
||||
~~df T L r~dT |
' |
7~ И Г |
+ ^ |
И Г |
|
|
~r — W |
|
f |
dOr |
’ |
|||||
|
|
da' |
dw |
|
Сп |
dw |
|
dw |
|
1 |
dP |
(3-13) |
||||
|
|
|
|
|
|
|
||||||||||
|
|
dt |
^ r dr |
|
|
r |
oU |
W |
dz |
|
^ |
p |
dz |
|
|
3-3. Уравнение энергии
Уравнение энергии представляет собой первый закон термодинамики, примененный к конечному движущемуся объему жидкости. Этот закон можно сформулировать следующим образом: изменение полной энергии данного объема газа (внутренней и кинетической) происходит под действием подводимого (отводимого) к веществу тепла и работы внешних объемных и поверхностных сил. Измене ние полной энергии равно
v |
v |
где U —внутренняя энергия единицы массы газа.
52
Из уравнения неразрывности следует
d Р dV ___ у
dt
и соотношение (3-14) принимает вид
|'|> (у + - f ) ^ ' = ] ' р Т ±{и + -у -) ■ (3-5)
V
Работа внешних поверхностных сил в общем случае описывается некоторым тензором Е, который для вязкого газа имеет вид
Е = 2р5—(Р+ — pclivcj, 3
где 5 —тензор скоростей деформации, равный
|
|
|
ди. |
|
-2( I— + |
д'о |
ди |
|
да |
|
|
|
дх |
|
дг |
дг |
|
дх |
|
S = |
1 |
,/ |
дv |
|
dv |
|
dv |
-Е |
да |
2 |
|
дх |
д у ) |
ду |
|
— — |
ду |
||
|
|
|
д~г |
|
|||||
|
2 |
I |
+ |
ди |
да + |
dv |
да |
|
|
|
дх |
дг |
ду |
дг |
дг |
|
|
р —динамическая вязкость. |
|
|
|
||
Используя теорему |
Гаусса —Остроградского, |
получаем |
|||
|
J С ЕdS = |
j div ЕС dV. |
(3-16) |
||
|
s |
|
г |
|
|
Работа |
внешних объемных |
сил состоит из работы: |
|||
а) гравитационных |
сил |
|
|
|
|
|
^ E ;Pa„CdV- |
(3-17) |
|||
|
v |
|
|
|
|
б) работа электромагнитного поля |
|
||||
|
|
\E ~ jd V . |
(3-18) |
||
Тепло, подводимое к газу теплопроводностью с пото- |
|||||
ком <7 пт |
, равно по |
теореме |
Гаусса — Остроградского |
||
|
•— j q dS г=— j |
div q dV, |
(3-19) |
||
|
S |
|
V |
|
|
53
В плазме важную роль играет тепловое излучение, однако для простоты анализа мы будем им пренебрегать.
Так как изменение полной энергии равно сумме работ внешних (поверхностных и объемных) сил и количеству подведенного извне тепла, то уравнение энергии в инте гральной форме имеет вид
j |
+ - f - ) d ^ = |
|
\ [ d i v E C t W — |
к |
~ |
i' |
\ |
|
— j* div qdV -|- |
j F,PaoC dV. |
\v
Учитывая произвольность выбора объема V, ' ы можем записать уравнение энергии в дифференциалы mil форме
P~dt{U + ~2 ~) = |
+ div ЕС - dive/ 4- |
С. (3-20) |
Заметим, что
DivEC = Div[2pS — — pdivC] —divPC. 3
Следовательно, работа, совершаемая поверхностными си лами, состоит из работы давления плюс работы сил трения. Для краткости вязкий член в уравнении энергии обозна чают буквой Ф, называемой диссипативной функцией Рэлея. Полностью вид диссипативной функции подробно рассматривается в курсах теоретической гидромеханики.
Напряжения давления совершают над газом работу, равную
diуp c = P d ivС+ CgradP. |
(3-21) |
Но из уравнения неразрывности следует
— -I- divp С:—0;
д(
— + р div С -\~Сgrad р = — + р div С — 0
или
div С —------
(IP
о at
54
Подставляя последнее соотношение в (3-21), получаем
divPC = - — - ^ - + CgradP = P p — f— ) + — |
- |
|
— |
|||||||||||
|
|
р |
dt |
|
|
dt |
\ |
(j I |
dt |
|
|
dt |
||
— О P |
d |
1 |
1 |
dP |
dP |
V± ( J L ) _ J U L . |
(3,22) |
|||||||
dt |
||||||||||||||
|
'• |
dt |
p |
p |
dt |
|
dt |
\ p j |
|
dt |
|
' |
||
Уравнение энергии в этом случае приобретает вид |
|
|
|
|||||||||||
d j ,, |
Р |
^ |
С- \ |
dP |
. т? |
|
,. |
— |
, |
- ■ —'■ |
||||
P Y l [ U + |
Т |
|
~ |
Ь 7• £ — div (/ -г И ’ -г РгРл,-С. |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
р |
(3-23) |
||
Как |
известно |
из |
термодинамики |
сумма |
|
|
||||||||
UН-----равна |
||||||||||||||
|
|
|
• |
г, |
|
|
|
|
|
|
|
р |
тормо |
|
энтальпии газами Вводя для краткости энтальпию |
|
|||||||||||||
жения |
/о = Н—— получаем |
уравнение |
энергии для энталь |
|||||||||||
пии торможения |
|
|
|
|
|
|
|
|
|
|
||||
|
Р |
-77- |
= |
В - ] + |
ц Ф |
+ В грав - С — |
d |
i v<7 - f |
F t f - . |
|
3(- 24) |
|||
|
|
dt |
|
|
|
|
|
|
|
ot |
|
|
|
Для того, чтобы уяснить принцип, приводящий к рас сеянию энергии, необходимо найти, какие причины приво дят к рост}' энтропии системы. Для этой цели домножим скалярно уравнение количества движения в виде:
р — = — grad Р 4 Grad 2 р 5 ---- р div С |
+ F,гРйв i |
7 X В |
||
dt |
|
|
||
на скорость С, получим |
|
|
|
|
Р — |
= — С-grad Р -f С Grad 2 p S ------ u div С |
+ |
||
dt 2 |
|
|
3 |
|
|
-j- F-C -|- C (j X B) |
|
|
|
и вычтем из уравнения энергии, получим |
|
|
||
Р at = / в — С - и х В) -р Div 2p S ----— р div С |
|
|||
—С ‘Grad |
2pS - — pdivC |
dive/ — div PC -|- Cgrad P; |
3
так как
/ • p — — ;2-i- с (7 x в)
65
и |
|
|
|
|
|
|
|
|
|
|
- div PC -Ь Cgrad P = — P div C = — |
rfp- |
— Pp — — ; |
||||||
то |
имеем |
|
|
|
|
rj |
dt |
cli |
p |
|
i |
P |
|
|
О |
—> |
|
||
|
~dU + |
P — |
d'v q |
|
Div C\ |
||||
P |
2pS ----— и di v C |
||||||||
dt |
dt |
p |
a |
|
|
3 |
|
|
|
|
P T |
dS |
-i------div <7 -j- |
2|.iS |
p div C |
Div C. |
(3-25) |
||
|
di |
В уравнении (3-25) вязкая диссипация так же, как и омический нагрев, всегда положительны. Уравнение опи сывает влияние теплопроводности и диссипации на энтро пию движущегося элемента газа.
Если тепловой поток q равен нулю, то очевидно, что
— >0, т. е. уравнение (3-25) выражает тот факт, что dT
омическая проводимость и вязкость приводят к необрати мому увеличению энтропии.
3-4. Уравнение состояния ионизированного газа
Уравнение состояния в общем случае имеет вид
F(P, р, Т) = 0 .
Однако в большинстве случаев достаточно точно плазму можно описать уравнением Клапейрона для совершенного газа
P = RpT. |
(3-26) |
В некоторых случаях, когда вязкий тепловой и магнитный пограничные слои малы так же, как и омический нагрев плазмы, то для областей, удаленных от пограничных слоев, течение можно считать изоэнтропическим. В этом случае давление Р и плотность р непосредственно без уравнения энергии связываются между собой. В этом случае в основу решения кладется не уравнение состояния, а уравнение изоэнтропического процесса в виде
— = const,
Рк
где k —показатель изоэнтропы, равный для совершенного газа отношению теплоемкости при постоянном давлении СР к теплоемкости при постоянном объеме Cv-
56
3-5. Система уравнений магнитной газодинамики
Неизвестными величинами в магнитной газодина
мике являются: Р, р, С, Т, Е, Н, рэ и /. Для нахождения этих восьми неизвестных рассматривается система газо динамических уравнений ((3-5) —(3-26)) совместно с урав нениями Максвелла и законом Ома. Система уравнений магнитной газовой динамики имеет вид
|
-j- divpC = О, |
|
|
—- = ------ grad Р |
|
Ра |
|
dt |
р |
|
|
р |
= Е - Т + \ 1Ф + F |
ав-С — div q f |
, |
dt |
|
н |
dt |
|
Р = RpT, |
(3-27) |
|
|
|
|
|
|
+с |
дВ |
|
|
r o t £ = — — , |
|
|
|
rot // ="7"-1 |
> |
|
|
|
dt |
|
dive0£ = р5,
1 — о (Е + С X В).
Эта система состоит из нелинейных уравнений, поэтому в общем случае решения этой системы возможно только численным путем. Однако, чтобы найти существенные свойства течения плазмы в магнитном поле, можно не прибегать к численному решению, а определить основные безразмерные величины, характеризующие это течение.
Существуют два метода получения безразмерных пере менных: обобщенный анализ в теории подобия и анализ размерностей. В теории подобия основная система урав нений приводится к безразмерной форме и выявляются безразмерные переменные —комплексы, характерные для рассматриваемых процессов. Переход от размерных вели чин к комплексам, называемым критериями подобия, при водит, во-первых, к уменьшению числа переменных и, во-вторых, отражает внутренние связи для подобной груп
57