повторяющимся индексам. В магнитном поле тензор вязкости имеет гораздо более сложный вид
[13].
Легко видеть, что вязкость связана с |
|
||||
неоднородностью потока скорости. В теории |
|
||||
горячей плазмы не менее важным может быть вклад |
|
||||
от вязкости, связанный с неоднородностью потока |
|
||||
тепла, который в обычной гидродинамике |
|
||||
пренебрежимо мал. |
|
|
|
|
|
Теперь остаётся |
вычислить |
столкновительный |
|
||
член |
Stvdp. |
Умножим |
подынтегнальное |
|
|
выражение в (V.1.5) |
на v, проинтегрируем по р и во |
|
|||
втором слагаемом переименуем немые переменные |
|
||||
p′ → p, |
p′′ → p′, p → p′′. В результате получаем: |
|
|||
|
|
1 |
|
|
|
|
Stvdp − |
m |
|
(V.2.15) |
|
∑ dpdp′dp′′f p f p′ w 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