Материал: Морозов Введение в теорию горячей плазмы Ч.1 2011

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

повторяющимся индексам. В магнитном поле тензор вязкости имеет гораздо более сложный вид

[13].

Легко видеть, что вязкость связана с

 

неоднородностью потока скорости. В теории

 

горячей плазмы не менее важным может быть вклад

 

от вязкости, связанный с неоднородностью потока

 

тепла, который в обычной гидродинамике

 

пренебрежимо мал.

 

 

 

 

Теперь остаётся

вычислить

столкновительный

 

член

Stvdp.

Умножим

подынтегнальное

 

выражение в (V.1.5)

на v, проинтегрируем по р и во

 

втором слагаемом переименуем немые переменные

 

pp,

p′′ p, p p′′. В результате получаем:

 

 

 

1

 

 

 

Stvdp

m

 

(V.2.15)

dpdpdp′′f p f pw p, p′′, p" p p" .

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

131

Stvdp mn V V .

(V.2.16)

Величина здесь представляет эффективную

частоту столкновений данных частиц с частицами сорта . Сечение кулоновских столкновений быстро падает с ростом относительной скорости частиц. В связи с этим в плазме интеграл столкновений даёт ещё один член, так называемую термосилу. Природа её такова. Рассмотрим незамагниченную плазму. Пусть в плазме находится пробный ион. С ним сталкиваются электроны, приходящие из области более высоких температур с той тепловой скоростью, которую они приобрели при последнем столкновении. Эта скорость соответствует температуре в точке, отстоящей на длину свободного пробега от места положения пробного иона. Они стремятся сдвинуть ион в область более низких температур. Из области более низких температур также приходят частицы, стремящиеся сдвинуть ион в противоположную сторону. Пусть потоки тех и других частиц равны. Тогда, в силу того, что сечение рассеяния более холодных частиц больше, на пробный ион действует результирующая сила RT, стремящаяся сдвинуть его в сторону более высокой температуры. Ясно, что эта сила должна быть пропорциональна температурному градиенту. По третьему закону Ньютона точно такая же сила, но с противоположным знаком, должна действовать на электроны. Точный расчёт показывает, что в

132

направлении, параллельном магнитному полю, эта сила для иона с Z 1 равна

RT 0. 71n T.

(V.2.17)

В замагниченной плазме выражение (V.2.1) дает параллельную составляющую термосилы, действующей на ионы. Перпендикулярная составляющая имеет более сложгый вид и приведена в работе [13].

Термосила является чисто столкновительным эффектом, но частота столкновений в неё не входит.

Cобирая вместе члены (V.2.2), (V.2.7), (V.2.8), (V.2.12), (V.2.13), и (V.2.16), в результате получаем уравнение баланса импульса:

 

mn

 

Vj

 

 

Vi

(V.2.18)

 

t

xj

 

 

 

 

 

 

 

 

 

 

 

 

P

 

Zen

Ei

 

1

V, H i

xi

 

c

 

 

 

 

 

 

 

 

 

mn Vi V i ij Ri

xj T

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

133

V.3. Второй момент (уравнение для температуры)

Второй момент получается аналогично нулевому и первому. Для этого надо умножить кинетическое уравнение на mv2/2 и проинтегрировать по импульсу. Снова разбивем скорость частицы на среднюю и хаотическую. В результате получаем:

 

3

n dT

PdivV divq Q.

(V.3.1)

2

 

dt

 

 

Здесь q

тепловой поток, Q

член,

описывающий источники тепла, потери тепла за счёт излучения, а также обмен энергией между разными компонентами плазмы, dtd t V . В

уравнение снова входит следующий момент от функции распределения q nm u2u . Однако цепочку уравнений можно оборвать, если предположить диффузионный характер теплопереноса. В этом случае поток тепла выражается через температуру

q T T

5

 

cnT

h, T .

(V.3.2)

2

ZeH

 

 

 

 

В это уравнение обычно не включают вязкое тепловыделение. В термоядерной плазме оно обычно мало по сравнению с другими источниками.

В гидродинамике обычно ограничиваются тремя моментами. Эта система уравнений в полном

134

виде приведена в работе [13]. Но довольно широкое распространение получил 13-моментный метод Грэда. Однако он годится лишь для численных расчётов и совершенно неприменим для качественного анализа процессов в плазме.

Задача

 

 

 

 

 

 

Доказать, что v v, H

f

dp n V, H .

p

 

 

 

 

 

 

Решение

 

 

 

 

 

 

В декартовых координатах

 

v, H

f

 

f

vyHz vzHy

p

 

 

 

px

 

pfy vxHz vzHx pfz vxHy vyHx .

Вычислим интеграл, соответствующий

x-компоненте уравнения (V.1.1).

dpxdpydpzpx

f

vyHz vzHy

 

 

px

pfy vxHz vzHx pfz vxHy vyHx .

Второй и третий члены при интегрировании дают ноль аналогично интегралам от похожих выражений в предыдущей задаче. Интеграл от

135

Источник: https://studfile.net/preview/16708796/