Батыгин&co
.pdf610 |
Глава XIV |
и те частицы, которые пролетают от нее сколь угодно далеко. Фактически в плазме любой заряд экранируется зарядами противоположного знака, поэтому с любой частицей взаимодействуют только те частицы, которые пролетают от нее на расстояниях, не превышающих радиуса экранировки. Для статистически равновесной плазмы радиус экранировки был вычислен в задаче 308 (радиус Дебая-Хюккеля):
где е и е' — заряды электронов и ионов, п и N — их концентрации. Величина Л называется кулоновым логарифмом. Пренебрегая слабой
зависимостью А от v, можно считать А = const, где const — число порядка 10.
855. |
|
|
|
|
|
F(v) =-fe*e*\J ^-Z^_/(v')(dv'), |
(1) |
где fj. = |
ттп |
— приведенная масса. |
|
|
m-\-m |
|
|
Полезно иметь в виду следующую электростатическую аналогию: выражение (1) можно записать в виде электрической силы F = qE, с которой
действует на заряд q = — ^е 2 е' 2 А «электростатическое поле»
E(v) = -gfadytp(v), |
(2) |
где
— «электростатический потенциал», удовлетворяющий уравнению Пуассона
856. Энергия пробной частицы не меняется при столкновениях с неподвижными бесконечно тяжелыми частицами. Изменение среднего импульса описывается уравнением
где F — средняя сила. Ее удобно вычислять с помощью электростатической аналогии, указанной в решении предыдущей задачи. Распределение
§ 1. Движение отдельных частиц в плазме |
611 |
по скоростям частиц среды описывается функцией /(v) = n5{\). Поэтому y(v) = n/v, E(v) = nv/v3,
F = -|£eV2 nAj, |
(2) |
F имеет характер «силы трения», стремящейся уменьшить направленную скорость частицы. Но это трение тем меньше, чем больше скорость частицы (F ~ 1/v2, «падающее трение»).
Интегрируя уравнение (1), найдем
v(*)=voexp[-|r],
где т = |
— - т ^ — ~ |
характерное |
время потери частицей направленной |
||||||||
скорости. |
4тге е п\ |
|
|
|
|
|
|
|
|
|
|
857. |
|
{ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
О |
|
|
|
при |
v < vn, |
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
/1 |
1 |
\ |
Щ- |
при |
v > VQ. |
|
|
|
4тге2е'2А( — + —) |
|
|||||||
858. |
|
|
|
|
|
\m |
m |
) v |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( |
|
|
|
|
/ 1 |
1 \ V |
|
|
||
|
ъ — ) |
-4тге |
2 |
е'2 |
Ап(— + -^)-% |
при |
v • v0 > vh |
||||
|
|
|
|
|
Vm |
m / v |
|
|
|||
|
|
|
|
|
|
1 |
7 ) ^ |
|
При |
V - VO<UQ - |
mm / v
859.На электрон, движущийся со скоростью v в среде неподвижных однозарядных ионов, действует сила трения
р _ 4тге4пА v |
Q4 |
m v3
(см. решение задачи 856). Отметим, что зависимость силы F от скорости можно получить и из следующих полукачественных соображений. Сила трения есть потеря импульса частицей в единицу времени из-за столкновений. Если среднее время между столкновениями т, а при каждом столкновении теряется импульс порядка полиого импульса частицы mv (это означает, что в результате столкновения электрон отклоняется на большой угол), то
F « ^ . |
(2) |
§ 1. Движение |
отдельных частиц в плазме |
613 |
т. е. F ~ 1/и2. Максимуму |
F, очевидно, будет соответствовать |
значе- |
ние и ~ г>т; при этом обе формулы (8) и (9) дадут одинаковое значение
(10)
mvr
Примерный ход функции F(u) представлен на рис. 134.
F
Рис. 134
Если поле в плазме Е < Fmax/e = Ей} то сила торможения при некотором значении и, удовлетворяющем равенству F(u) = еЕ, превысит ускоряющую электрическую силу еЕ, и электроны не смогут больше ускоряться. Это — область значений поля Е, при которых имеет место обычный закон Ома. В случае Е > ED ускоряющая сила превышает торможение, и электроны получают возможность ускоряться неограниченно2. Это явление получило название «убегающих электронов».
Подставляя v2 = ЪкТ/т в формулу (10), получим
i. ел |
гуг |
кТ |
(И) |
|
Ажпе2 |
||||
|
|
Точный расчет для этого же случая дает ([28], в. 1)
ED = 0,214^ . |
(12) |
Наша порядковая оценка дала результат, близкий к точному значению.
1 Критическое значение поля Е = ED называется драйсеровским.
2 В действительности из-за коллективных эффектов электронный газ как целое не ускоряется, его сопротивление может даже возрастать.
§ 2. Коллективные движения в плазме |
615 |
Из граничных условий и (6) находим
где «о — новая постоянная, имеющая смысл скорости в средней плоскости х = 0. vo можно выразить через градиент давления:
ах0 ch^~1dp
Магнитное поле определяется из (2), (7) и граничных условий Hz(±a) = 0:
|
|
(х\ |
а |
х |
|
|
Hz(x) = |
Awn |
\а) |
хо |
хо |
. |
(9) |
^-y/^rjvo |
|
|
||||
Отношение а/хо = М называется числом Гартмана. При М с |
1 имеем |
|||||
a2 |
dp |
^ " " " l 1 |
|
|
|
|
~ Vdz' |
a?)' |
|
|
|
как в обычной гидродинамике. Магнитное поле Hz = 0 в первом порядке по числу Гартмана. Продольное поле Hz появляется только в следующих приближениях.
В противоположном предельном случае М > 1 получаем
Сравнение (10) с (11) показывает, что средняя скорость движения уменьшается с ростом Но, а профиль скоростей становится более плоским в средней части потока, но резко меняется в слое толщиной х0 у стенок. Продольное магнитное поле в этом пределе имеет вид
= ^Р-ЛШ |
-expf-^1 Л 1 ) . |
(12) |
сМ dz\a |
L xoi XQ) |
V |
Как видно из формулы, оноубывает с ростом числа Гартмана. Наибольшую величину Hz имеет при М и 1.
616 |
Глава XIV |
Плотность тока в движущейся жидкости вычисляется из уравнения Максвелла j = -£- rot Н. Отлична от нуля только у-компонента тока:
Создаваемое им магнитное поле Hz равно нулю всюду вне области, занятой жидкостью. Там остается только поперечное поле HQ.
|
h — |
|
|
|
h — |
|
861. |
v(x) = vo |
-. |
Плотность тока j(x) = |
° |
-. Этот ток |
|
|
s h ^ |
|
|
Н |
Х Ъ |
% |
создает магнитное поле |
|
|
|
|
|
|
|
„ |
, ч |
Awqvoc h |
Ш ~ c h Ш |
|
|
|
Л |
' |
Нохо |
s h ^ . |
|
|
обращающееся в нуль при \х\ ^ а. |
|
|
|
|||
862. |
Магнитное поле имеет одну проекцию |
|
|
о
Интегрируя уравнение движения (XTV.7)с граничным условиемр|г > а = = 0, находим
где Н = Щ:J rj(r) dr при г < а, Н = 2J/cr при г > а.
о
Чтобы связать силу тока J с Т и N, полагаем р = 2п(г)кТ, где к —по- стоянная Больцмана, и интегрируем обе части (1) по площади поперечного сечения столба плазмы. Получим
J = 2cVNkT. |
(2) |
При Т и 108 °К и N и Ю1 5 частица/см3 (значения, характерные для термоядерных исследований) имеем
§ 2. Коллективные движения в плазме |
617 |
863.Ток должен течь потонкому поверхностному слою. Тогда внутри столба будет постоянное давление
864.Беря проекцию уравнения (XIV. 12) наосьг иподставляя v = и£,
v = const, получим уравнение для определения Нг:
_ = __Яг -„—. |
(1) |
Решение этого уравнения выражается через произвольную функцию F от аргументов г —vt, fl и а:
|
Hr(r,d,a) = ±F(r-vt,d,a). |
|
|
(2) |
||
Граничное условие имеет вид |
|
|
|
|
||
Я г | г = а |
= НОг{-д, а +ГО)= ^F(a |
- |
vt, -д, а) |
|
(3) |
|
(аргумент a—uty |
Hor |
написан в связи с переходом в неподвижную систему |
||||
координат). Таким образом, |
|
|
|
|
||
|
F(a - vt, ti, a) = о2Я0г(1?, а + ГО). |
|
|
|||
Следовательно, (2) запишется в виде |
|
|
|
|
||
Нг(г, |
0, a, t) = ( f ) 2 Я 0 г (*, а - {Г |
~уа)П |
+ |
ГО). |
(4) |
|
Таким же путем находим |
|
|
|
|
,а - ( Г ~„°) П + |
ГО). |
(5) |
Из уравнения divH = 0 вытекает следующая связь между проекциями вектора Но:
618 |
ГлаваXIV |
При Нов = 0 находим Н# =О,
если положить /(#) = 0, то будем иметь
Я„(г, 0, a, t) = ^ Я О г (Ч а - ( Г ~ / ) П + ta) sin tf. (6)
Паркер использовал рассмотренную модель для описания межпланетного магнитного поля, создаваемого потоками солнечной плазмы (солнечным ветром). В модели межпланетного магнитного поля Паркера Н# = О,а Нг и На даются формулами (4), (6). Измерения межпланетного магнитного поля, произведенные на спутниках и ракетах, показывают, что усредненное магнитное поле вблизи орбиты Земли удовлетворительно описывается моделью Паркера.
865. Силовые линии имеют вид спиралей Архимеда:
г = ^г(а- а0), |
а0 = const, |
в = arctg ^ф- « 56°; |
Я а 4,5 • 1(Г5 э. |
866. е± = 1 Н Y~,где р — плотность плазмы. Найденное значение е±_ получается из результатов задачи 321 в предельном случае из —> 0.
867. из = изр = у га6 • г д е m ~~м а с с а электрона.
868.При из < изр, R = 1,
Е= - — ^ Е о exp[-qz - iwt],
где q = ^ \ / - | ~ 1 > ' г = %•> Ео — амплитуда падающей волны. Глубина проникновения Ь = \ = — с ; 5 Й ^ - при из <Сизо. Затухание поля
вызвано недиссипацией энергии, а возникновением токов в плазме, которые создают поле противоположного знака.
§ 2. Коллективные движения в плазме |
619 |
При UJ > и)р
где q = лЛ^2 —ш%/с, волна распространяется в плазме без затухания.
869. Представим радиус-вектор частицы в виде
(1)
где vo — скорость частицы в отсутствие поля (тепловая скорость); Ro — радиус-вектор при t = 0; Ri (t) —добавка, обусловленная действием электрического поля плоской волны (магнитным полем пренебрегаем, считая частицы нерелятивистскими).
Величина Ri удовлетворяет уравнению
mRi = еЕ0 ехр[г(к • Ro + к • vot + к • Ri - ut)}. |
(2) |
В показателе экспоненты можно пренебречь слагаемым к • R b считая выполненным неравенство kR\ « 1. В этом приближении, линейном по Ео, решение (2), соответствующее вынужденным колебаниям, имеет вид
|
|
R o - i ( w - k . v 0 ) t ] |
|
Скорость частицы выражается в виде1 |
|
||
v(t) = vo + И |
_ ^ |
_ exp^ . Ro - Г(Ш - k • vo)t]. |
(4) |
Ток, создаваемый одной |
частицей, начальные значения радиуса-вектора |
||
и скорости которой равны соответственно Ro и vo, запишем в виде |
|
||
h(r,t) |
= ev(t)(r-R(t)), |
(5) |
'Расходимость выражений (3) и (4) при kvo = ш связана с некорректнымрассмотрением «резонансных» частиц, т.е. частиц, скорость которых удовлетворяет условию kvo = ш. Чтобы избежать этой трудности, предположим, что в плазме отсутствуют частицы с такими скоростями, т.е. исключим из рассмотрения интервал скоростей, удовлетворяющих неравенству |v —vo| < vo.