![](/user_photo/_userpic.png)
книги из ГПНТБ / Эккер, Г. Теория полностью ионизированной плазмы
.pdf§ 1. Вывод кинет, уравнений из цепочки уравн. Б БГК И 213
Используя разложение (1.22) и пренебрегая, как всегда, корреляционными функциями третьего порядка, можно первые два уравнения цепочки (1.20) и (1.21) пере писать в виде х)
dt дт1 V оt /столки
|
J ^ S t a ( r t, . . . , v 2;t)d r2dv2 |
(1.65) |
|||
и |
|
|
|
|
|
д |
Vz- i;)иг2 gi2' |
|
|
|
|
\ dt ^ V i - 3ri |
|
|
|
|
|
|
д |
v*; |
|
v2; O + ftsH - |
|
|
г /9/, 9дфФ 1Яз . |
|
|||
+ -N |
, dfi |
дф2з |
£13) dr3dv3. |
(1.66) |
|
|
аГ1 ё23 |
1 av2 |
аг5 |
|
|
Здесь по-прежнему не учитываются внешние силы. Кроме того, в этих уравнениях пренебрегается влиянием само согласованного поля на поток в пространстве скоростей, что дает возможность отбросить третий член в левой части уравнения (1.20) и четвертый и шестой члены в ле вой части уравнения (1.21). Очевидно, что для однородной системы такая процедура выполняется точно. Неоднород ность системы можно считать допустимой, если соответ
ствующее самосогласованное |
поле мало по сравнению |
со средним корреляционным |
полем. |
Следует подчеркнуть, что малую пространственную зависимость ft в (1.65), (1.66) и последующих формулах следует понимать в указанном выше смысле.
Как и в случае уравнения Больцмана, мы стремимся исключить член gi2 из (1.65) с помощью уравнения (1.66). Теперь вспомним, что в безразмерном представлении (1.4)
первый член в правой части |
(1.66) |
пропорционален |
Псв, |
|
а второй — произведению |
ПсвПпл. |
Следовательно, |
вто |
|
рым членом в (1.66) можно пренебречь, |
как членом более |
|||
высокого порядка малости [см. (1.64)]. |
|
|
||
1) Хотя используются частные функции |
распределения, но |
в силу основного предположения, данного на стр. 126, эффект от всех полевых частиц, представленный правой частью уравнения, записывается в виде произведения N на эффект только от одной Нолевой частицы с координатами (г2, v2).
214 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций
Физически это означает, что пренебрегается взаимо действием полевой частицы с третьей частицей, находя щейся вблизи пробной частицы, что приводит к ослабле нию механизма экранирования.
Учитывая (1.23), пренебрежем также членом gl2 в пра вой части (1.66). Тогда мы получим уравнение в частных производных относительно функции gl2:
Это уравнение можно переписать в виде
При интегрировании вдоль прямолинейных траекторий
v1= v2 = 0, г12 (<') = r12 (t) + \i2(t'— t) |
(1.69) |
второй член в правой части уравнения (1.68) дает нуль. Однако, используя приближенное интегрирование вдоль прямолинейных траекторий, следует помнить, что пере менные гг и уг связаны законом Ньютона. Из оценки вели чины
можно показать, что относительная ошибка при таком приближенном интегрировании вдоль прямолинейных траекторий порядка Псв. Обозначая оставшиеся члены
вправой части уравнения (1.68) через
А(гь г2, vt, v2;' t) =
после интегрирования получаем
|
t |
°°)* (1.72) |
g}2 |
t) — J -'1 ( • • • i t ) dt + gj2 (• • M |
§ 1. Вывод кинет, уравнений из цепочки уравн. |
Б Б Г К И 215 |
|||||||
Коррелированное |
|
состояние, соответствующее |
време |
|||||
ни |
t, |
получается |
|
из некоррелированного |
состояния |
|||
#12 (■ |
• •> |
—оо) = |
0 |
при t = —оо. |
на |
длине |
||
Предполагая, |
что |
изменением |
функции f t |
|||||
взаимодействия rc (L |
гс) можно |
пренебречь, |
получаем, |
что подынтегральное выражение зависит только от отно сительного расстояния между частицами. Вводя новую
переменную |
х = |
t — t', |
находим |
|
|
|
||||
gi2(|ri — г2|, vt, v2; t) = |
|
|
|
|
|
|||||
|
оо |
|
|
|
|
|
|
|
|
|
|
= j |
^ (ri(0 —ViT, r2(<) — v2T, v1? v2; t |
t) dx. |
(1.73) |
||||||
Определяя поток J |
в пространстве скоростей как |
|
||||||||
|
|
|
|
(Ч г) |
|
= — |
|
|
(1.74) |
|
|
|
|
|
\ |
o t /столки |
|
|
|
|
|
из формул |
(1.65), |
(1.71) |
и (1.73) находим |
|
|
|||||
|
|
|
|
|
о |
|
|
|
|
|
X R |
r ( ^ |
I - ^ |
) |
/l(ri’ Vi; 0 /i(rb V 2; |
|
dx. |
(1.75) |
|||
L |
1 |
' |
1 |
2 ' |
|
|
|
J ri2—V12T |
|
|
Используя |
разложение |
|
|
|
|
|
||||
l/i (Гь V t ; |
0/i(*i, |
V 2 ; / ) ] t _ T = |
|
|
|
|
||||
|
|
|
|
= |
( /i/i) f - ( i- x ) [ - ^ - ( /1/i)]t + . . . |
(1-76) |
и подставляя его в выражение (1.75), получаем уравнение Ландау с запаздыванием.
Если функция распределения заметно не изменяется в течение времени взаимодействия {Т^> tc), можно огра ничиться первым членом разложения (1.76), т. е. пре небречь запаздыванием в системе. В результате мы полу чим
•>= - N Ш 2 J <Е><Е>.. {/• ('*• ™ <)a/i(x i; 0 -
- / , ( r . . v , ; l ) l),l(y |
' ,>}<iv2. (1.77) |
216 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций
Здесь коэффициент ( Е)(Е )Г2 определяется |
выраже |
нием |
|
<Е)(Е)Г2= j dr2 j Е [г12])(Е [г12 —v12t] dx, |
(1.78) |
о |
|
аЕ является кулоновским полем.
Тензор ( Е)(Е )Г2 представляет собой автокорреля
ционный тензор второго порядка, зависящий от vt и v2 и вычисленный в приближении прямолинейных траекто рий. Индекс г2 указывает, что усреднение проводится только по г2.
Уравнения (1.74) и (1.77) вместе представляют собой уравнение Ландау, хотя в них не приводятся еще точные значения автокорреляционных функций поля. Уравнение Ландау в виде, приведенном здесь, соответствует урав нению Больцмана для слабых взаимодействий [см. (1.8)]. Поэтому интересно отметить, что при его выводе мы вынуждены были сделать те же самые основные допу щения. Действительно, условие gl2 —оо) = 0 соответ ствует предположению Больцмана о начальном моле кулярном беспорядке. Условия (1.53) для пространствен ной и временной зависимости функции распределения отражаются здесь в требованиях L гс и tc.
Вычисление корреляционных функций поля
Это вычисление представляет собой геометрическую задачу. Для 6е решения удобно выбрать цилиндрическую систему координат с осью z, направленной вдоль вектора относительной скорости g = v2 — vt. Радиальную коор динату мы вновь обозначим через Ь (см. стр. 207).
В этой системе координат поле выражается следующим образом:
|
Ъcos ф' |
|
Б(г» —г* + gr) |
е |
(1.79) |
[b2+ (Z- gT)2]3/2 b sin ф |
,2 —ГР
§ 1. Вывод кинет, |
уравнений из |
цепочки уравн. |
Б Б Г К И |
217 |
При этом корреляционный тензор записывается как |
|
|||
Е |0)(Е |, = е2 {[Ь2 + z2] [Ь2 + (z- gx)2]}- 3/2 х |
|
|
||
Ъ2 cos2 ф |
Ь2 sin ф cos ф |
(z— gx) Ъcos ф\ |
|
|
Ь ^ тф совф |
Ь2з т 2ф |
(z — §х)Ь8т ф | , |
(1.80) |
|
(bz cos ф |
bz sin ф |
z (z — gx) |
/ |
|
а элемент объема как |
|
|
|
|
|
dr2 = dzbdbdq. |
(1.81) |
Что касается пределов интегрирования по прицельному параметру Ъ, вспомним обсуждение области применимо сти метода Ландау — Фоккера — Планка, которое при вело к значениям Ьмин = rw и Ьмакс = X,D. (Точное дока зательство пределов применимости этого метода будет дано в п. 1.4 в связи с выводом уравнения Боголюбова — Ленарда — Балеску.) Интегрирование по углу ф в (1.80) тривиально и дает для недиагональных членов значение, равное нулю. Интегрирование по z дает для члена A zz тоже нулевое значение вследствие соотношения симметрии
j |
Azz ( |
z, gx)dx = |
J |
Azz (z, gx) dx. |
(1.82) |
||
о |
|
|
о . |
|
|
|
|
Члены A xx |
и A jn, интегрируются |
элементарно, |
так что |
||||
корреляционные |
функции поля |
принимают вид |
|
||||
|
|
|
|
/1 |
0 |
0\ |
|
|
<Е)(Е>Г2 = 2яе2- ^ Ц о |
1 |
0 , |
(1.83) |
|||
или |
|
|
|
\0 |
0 |
0 / |
|
|
|
|
|
|
|
|
|
|
(Е)(Е>Г2= 2яе2 ( ^ |
А - ) |
(йг1—g ) ( g ) , |
(1.84) |
где через I обозначается единичный тензор. Подставляя
этот результат и выражение (1.77) в (1.74), получаем окончательный вид уравнения Ландау [16]
(М±.\ \ б* )с
х |/ i (ть
О ЛТ I «2 \ 2 1
:2я г ( - )
v2; t) 9fi (ft, vi; d\i
А |
9 |
f |
g l ^^ |
-^ х |
|
t) |
7i (4, |
Vj! t) |
9fi (r«, v2; t) |
} d \2. |
|
|
|
|
|
d\2 |
(1.85) |
|
|
|
|
|
218 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций
Уравнение Ландау является кинетическим уравне нием для случая слабого взаимодействия, не учитываю щего эффекты экранирования. Интересно из цепочки уравнений ББГКИ получить эквивалентное уравнение Фоккера — Планка. Это можно сделать, исходя из (1.77) и связывая корреляционные функции поля с коэффи циентами переноса в пространстве скоростей.
Связь корреляционных функций поля с коэффициентами переноса в пространстве скоростей
"В последующем изложении мы используем два первых момента изменения скорости в единицу времени: коэффи циент динамического трения
|
о |
(1-8б> |
|
|
|
и коэффициент диффузии в пространстве скоростей |
||
|
Тс |
Тс |
(Av1)(Avi) = ± |
\ dt' |
J dt"E |4<)(Ег") . (1.87) |
|
Ь |
о |
Здесь Е — полное поле, действующее на рассматривае мую частицу 1; скобки () означают усреднение по фазо вым координатам всех остальных частиц; тс — среднее время свободного пробега, которое пропорционально вели
чине |
TpA/lnA [см. (3.1.7)]. |
|
|
|
в |
Удобно рассмотреть сначала коэффициенты диффузии |
|||
пространстве скоростей. |
Простая замена переменной |
|||
t" |
= |
t' + т в (1.87) дает |
|
|
|
|
Тс |
т с - t ' |
|
|
(Av1)(Av1) = ^ - ( ^ - ) 2 j dt' |
j <E|r )(E|r+ t>dt, |
(1-88) |
где подынтегральное выражение представляет собой авто" корреляционную функцию поля. Оценим (1.88), исполь" зуя два приближения, которые уже применялись ранее-
§ 1. Вывод кинет, уравнений из цепочки уравн. БВГКИ 219
а) |
Взаимодействия |
на |
расстояниях, |
больших |
гс, не |
|||
дают вклада в интеграл, |
так что существует максимальное |
|||||||
время |
Д = r j v c, |
в |
течение |
которого |
могут взаимодей |
|||
ствовать частицы, |
|
причем |
|
|
|
|
||
|
( Е |Г)(Е \t.+x ) = |
0 для |
т > Д. |
(1.89) |
||||
б) |
Изменения |
Д |
малы |
в |
смысле, |
указанном |
выше, |
|
т. е. Д |
(t — х) ж Д |
|
(t), |
когда использовалось «приближе |
ние в отсутствие запаздывания». Точнее говоря, если Т есть характерное время изменения Д , то условие Т^> tc
приводит |
к |
тому, что ( Е |
| г)(Е |
| t'+x)\ не зависит |
от t' |
||
для t' < |
Т. |
Отсюда |
следует, |
что |
|
|
|
|
|
|
То |
|
|
|
|
(Av1)(Av1) = |
( ^ ) 2^ |
\d t' |
j |
(Е |0)(Е |t) dx = |
|
||
|
|
с |
о |
|
1 Ц |
< ( С |
|
|
|
|
|
|
|
оо |
|
|
|
|
= 2 ( ^ ) 2 |
j< E |0)(E|t)dr. |
(1.90) |
||
|
|
|
|
|
|
о |
|
В этом выражении поле Е представляет собой линейную суперпозицию кулоновских полей всех частиц. Так как в методе Ландау предполагается, что тождественные поле вые частицы не коррелируют между собой, коэффициент
диффузии в |
пространстве |
скоростей |
(1.90) может |
быть |
||
записан как сумма N средних полей, возникающих при |
||||||
парных взаимодействиях, т. е. |
|
|
|
|||
<Av1)(Av1>= 2N ( - ^ ) 2 j |
(Е)(Е>Г2Д (г15 v 2; t)dvz. |
(1.91) |
||||
Здесь мы использовали выражение |
(1.78). |
|
||||
Рассмотрим теперь коэффициент динамического трения |
||||||
|
|
Тс |
|
|
|
|
|
( A v ^ - ^ - i - |
[ <E[r12(f')]>d*'. |
(1-92) |
|||
|
п г Т с |
J0 |
|
|
|
|
в котором подынтегральное |
выражение |
можно записать |
||||
в виде |
|
|
v |
г |
|
|
|
|
|
|
|
||
Е [г12 (*')] = |
Е [г12 |о + g |о |
- |
2 - М |
j |
Е [г12 (Г")] dt”] . |
|
|
L |
|
0 |
0 |
|
(1.93) |
|
|
|
|
|
|
220 Гл. 4. Неравное, состояния кулон, сист. с учетом корреляций
Поскольку мы ограничиваемся случаем слабых взаи модействий, приводящих только к малым отклонениям от прямолинейного пути, то имеем
Е [г12 (£')] Е [г12 -f- gt'] |
(1.94) |
или, интегрируя по частям, |
|
E[r)2(0 ] = E [ri2 + g r ] - ^ | ( / 2 -£_ j (*'_ **)Е ||.Л \ |
(1.95) |
о
Первый член, стоящий в правой части этого равенства (возникающий в приближении прямолинейных траекто рий), не дает вклада в среднее значение коэффициента динамического трения (1.92). Поэтому получаем
<Лт*> ■= - (= •)1i |
} л ' |
1 |
(»' ■-п < Ц | , ■:Е |-> ■**• |
(1.9!)) |
Вводя новую переменную x = t' — t", находим |
|
|||
|
То |
Г |
|
|
- ( - S - ) 4 |
J |
\ |
|
(1.97) |
|
о |
о |
|
|
или, усредняя по значениям т — t',
<ЛЧ> = 2 ( £ ) 2-dvik - тv, ] |
*' 1 <Е|.)(Е|»>*= |
|
|
|
0 |
0 |
|
е \2 |
д |
tc |
|
|
(1.98) . |
||
= 2 Ш |
^ Г - 1 ( 1 - ^ ) < Е |« ) ( В |.) Л . |
||
Так как, согласно |
общему предположению, отношение |
||
xc/tc л; Л/In Л 1, окончательно получаем |
|
||
|
оо |
|
|
<Ду,) = 2 ( ^ . ) г ^ - . |
|< Е |,)(Е |,> Л = |
|
|
|
о |
|
|
2 N ( ~k) |
|
j <Е)(Е)Гя/1 (rj, v„; t)dv9, |
(1.99) |
§ 1 . Вывод кинет, уравнений из цепочки уравн. В В Г К И 221
что приводит к важному соотношению между коэффи циентом динамического трения и коэффициентом диффу зии в пространстве скоростей:
<Avi>=±' ^ r ‘<Avi^ Avi^ |
(1.100) |
Уравнение Фоккера — Планка для слабого взаимодействия
С помощью соотношений (1.91) и (1.100) легко преобра зовать уравнение Ландау для слабого взаимодействия к уравнению Фоккера — Планка. С этой целью в урав нения (1.74) и (1.77) подставим выражение (1.91) и в ре зультате получим
( i r |
L |
[ т <Av.)<Av, >- ! £ - |
|
- j v ( - i -)2 J |
<E)(E)r,/i(r„ v,; t). |
?>riv,] . (1.101) |
Интегрирование по частям и использование выражения
(1.91) дает
|
6 t |
|
= •/- •f4 (AviXAvj)' ' |
• d \ i — |
{-гг) |
|
|
||
\ |
|
/столки d \ i \ 2 > |
|
|
|
|
|
|
■ ( 1- 102) |
Отсюда, учитывая (1.100), находим |
||||
( ж - ) ™ |
„ „ н - - ^ г '« Дт‘>« + |
|
Это уравнение является уравнением Фоккера— Планка
для слабого взаимодействия. Впервые полученное Фоккером и Планком для случая броуновского движения больших молекул [17, 18], оно в общем случае описывает движение частиц, испытывающих большое число слабых отклонений.
Приведенный выше вывод уравнения Фоккера — Планка не совпадает с выводом, данным самими автора ми. Следуя общей линии нашего изложения, мы получили
222 Гл. 4. Неравное, состояния кулон, сист. с учетом Корреляций
это уравнение из цепочки уравнений ББГКИ для част ных значений параметров Псв = О (е) и Ппл = О (е). Кроме того, мы использовали приближение парных столк новений, интегрирование по прямолинейным траекто риям и приближение в отсутствие запаздывания.
Для сравнения ниже приведем другой вывод урав нения, данный Чандрасекаром [19], который является более-коротким и не требует приближения парных столк новений, интегрирования по прямолинейным траекто риям и отсутствия запаздывания. Но при этом, как мы увидим из дальнейшего, уравнение не будет содержать всей той информации, которая была получена выше. Будем исходить из основного уравнения Чепмена — Колмого рова для марковских процессов
/ (v; £-f-A<)= j /(v — Av; t) PAt (v — Av, Av)dAv. (1.104)
Здесь функция / есть одночастичная функция распределе ния однородной системы в отсутствие внешних сил. Через PAt (v, Av) обозначена вероятность того, что частица, которая в момент времени t имела скорость v, за интервал времени At изменит свою скорость на величину Av. Рас кладывая подынтегральное выражение (1.104) в следую щий ряд:
/ (v — Av; t) PAt (v — Av, Av) =
|
|
= / (v; t) PAt (v, Av) - |
Av |
{fPAt) + |
|
|
||
|
|
+ 4 ( Av- | t)2( / ^ ) - |
+ |
--- |
, |
(1-Ю5) |
||
из (1.104) |
получаем |
|
|
|
|
|
|
|
/(v; |
■+ A«)->(v; . ) = _ ^ , , [(Av)/1 + |
|
|
|
|
|
||
|
|
+ T ^ - { - i r - 1<AvXAv>fl } - |
+ |
•■• |
• |
(‘ -106) |
||
При |
этом |
были использованы |
следующие |
соотношения: |
||||
|
|
P a, cZAv = |
| |
At (Ах ) |
\ . |
|
(1.107) |
|
|
|
|
[ Ai(Av)(Av) J |
|
|