Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Батыгин&co

.pdf
Скачиваний:
374
Добавлен:
10.02.2015
Размер:
7.66 Mб
Скачать

610

Глава 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)

612

Глава XIV

При таком столкновении электрон подходит к иону на расстояние, на котором его кинетическая энергия — порядка потенциальной:

mv2 ^

е2

/о\

— ~т-

(3)

Это приближенное равенство позволяет оценить сечение столкновения

 

 

4

(4)

 

ТТГ1Г

 

и среднее время между столкновениями

 

 

T r a J^ M i

 

( 5 )

nav

4тгпе4

 

Подставляя г в (2), находим F те 4тгпе v

или, учитывая тормозящий ха-

рактер силы,

mv

 

 

 

 

Р и

( 6 )

что отличается от (1) отсутствием кулонова логарифма А. Но это естественно: при оценке по формулам (2)-(5) мы не учитывали далеких столкновений с малыми передачами импульса, вклад которых и дается кулоновым логарифмом.

Усредним теперь формулу (1) по возможным скоростям электрона. Для этого положим

v = u + vT,

(7)

где vT — тепловая скорость, и — скорость, приобретаемая под действием электрического поля Б. При и < « т можем положить в знаменателе выражения (1) v3 RitiJ.B числителе же нельзя пренебречь и по сравнению с vT, так как при усреднении по направлениям тепловой скорости получим vT = 0. В итоге будем иметь

^ щ

( )

mv*

 

где под vT теперь нужно понимать величину порядка средней тепловой скорости. В случае распределения Максвелла v2 = ЗкТ/т. Таким образом, при и <g.vT получаем F ~ и.

При и » vT полагаем v и и и получаем

ти

§ 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 В действительности из-за коллективных эффектов электронный газ как целое не ускоряется, его сопротивление может даже возрастать.

614

ГлаваXIV

§2. Коллективные движения в плазме

860.Естественно предположить, что скорость движения жидкости направлена вдоль оси z и зависит только от поперечной координаты х. Поскольку движущаяся проводящая жидкость увлекает за собой силовые линии магнитного поля, то при движении должна возникнуть про-

дольная составляющая магнитного поля Hz(x). Таким образом, неизвестные функции v и Н ищем в виде v(0,0,v(x)); H(HQ,0,HZ(X)); при этом уравнения (XTV.9), (XIV. 10) удовлетворяются тождественно. Уравне-

ния (XTV.7), (XTV.8) принимают следующий вид:

dv

с2

d?Hz

 

dx +

4тп7Я

d 2 ~ '

( )

d?v__Ho_dHz__l±(n,H?_\

dx* ~ *Щ dx ~ г, dz \ Р + 8тг Г

и2

Из последнего равенства следует, что р + ^ - зависит только от z. Но

-й- p + f -

= ^ = const,

(4)

dz V 87Г /

dz

 

так как Я2/8тг от z не зависит. Поэтому равенства (1) и (2) представляют собой систему обыкновенных линейных уравнений для определения неиз-

вестных функций v(x), Hz(x). Исключая из них —г-5-, получим уравнение

относительно и = ^-: dx

d?u

1

п

т

с

Рп

/гч

dx

XQ

 

 

tiQ

у

а

 

из которого следует

 

 

 

 

 

 

 

/

Г

— Л

 

Г _

П

\

 

 

 

 

 

 

 

-С.

(6)

Граничные условия имеют вид v(±a) = 0, так как вязкая жидкость у стенки неподвижна. Кроме того, из соображений симметрии следует v(x) = v(—x).

§ 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Я(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.