
книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений
.pdf3.2. Анализ Боголюбова ББКГИ-уравнений |
127 |
Кроме того, / 8 зависит от времени неявно, через функцию/!. Примером такой зависи мости может служить функция f s для s ста тистически независимых частиц:
|
|
|
|
|
|
/. = ,П /і (*)• |
|
|
(3.17) |
|
|
|
|||
|
|
|
|
|
|
|
h=l |
|
|
|
|
|
|
||
|
Наконец, на гидродинамической стадии |
|
|
|
|||||||||||
|
|
|
|
Т о все функции распределения f s |
|
|
|
||||||||
зависят только |
от |
п, и и Т, т. |
е. первых |
|
|
|
|||||||||
пяти моментов функции /4. |
изображена |
|
|
|
|||||||||||
|
На |
рис. |
3.3 схематически |
|
|
|
|||||||||
последовательность только что рассмотрен |
|
|
|
||||||||||||
ных случаев (когда при t = 0 только одна |
|
|
|
||||||||||||
частица из всего стационарного газа при |
|
|
|
||||||||||||
ходит в движение). |
|
|
|
|
|
|
|
|
|||||||
|
При изменении состояния системы от |
|
|
|
|||||||||||
начального |
хаотического до конечного хо |
|
|
|
|||||||||||
рошо организованного равновесного со |
|
|
|
||||||||||||
стояния нам требуются резко различаю |
|
|
|
||||||||||||
щиеся уровни описания системы. Вначале |
|
|
|
||||||||||||
необходимо |
знать |
не меньше |
чем N-ча- |
|
|
|
|||||||||
стичную функцию распределения. В конеч |
|
|
|
||||||||||||
ной, равновесной, стадии достаточно знать |
|
|
|
||||||||||||
термодинамические |
|
переменные, |
|
дающие |
|
|
|
||||||||
несравненно менее подробное описание си |
|
|
|
||||||||||||
стемы. |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Анализ Боголюбова иерархии ББКГИ- |
|
|
|
|||||||||||
уравнений касается |
кинетической стадии |
|
|
|
|||||||||||
и включает следующее ограничение: |
|
|
|
||||||||||||
|
|
|
lim |
( |
^ |
Ь |
Ч |
- г Ь |
т |
. |
(3.18) |
Р и с . |
3.3. а — началь |
||
|
|
|
N- |
|
|
||||||||||
|
|
|
Ѵ-*оо |
|
|
|
|
|
|
|
|
ная стадия: 0<1 < т0, /_ѵ |
|||
где і/ѵ является конечным малым пара |
определяет состояние га |
||||||||||||||
метром. |
(Соответствующим безразмерным |
за; б —кинетическая ста |
|||||||||||||
малым параметром будет г30/ѵ, где г0—ра |
дия: |
т0 < t < t 0. |
Части |
||||||||||||
цы приобретают |
некото |
||||||||||||||
диус взаимодействия.) В этом |
предель |
рый |
элемент сходства,. |
||||||||||||
ном случае (см. задачу 2.27) уравнение БИ г |
fs = fs [/il; 6 — гидроди |
||||||||||||||
имеет |
|
вид |
|
|
|
|
|
|
|
|
намическая стадия: t0< |
||||
|
|
|
|
|
|
|
|
|
<1 t <_ Т0. Состояние га |
||||||
dFt |
, |
|
f |
„ |
|
1 |
1 |
|
d(l + l)Qul+1Fi+ = 0. |
за близко к равновесно |
|||||
|
|
|
му и задается через функ |
||||||||||||
dt |
|
■LiF I |
|
|
■3 |
|
|
ции и, и, |
Т. |
||||||
|
|
|
|
|
|
|
i= l |
|
|
= |
|
(3.19) |
|
|
|
Напомним, |
|
что Fi = |
Vlfi и і |
( |
Xi, Pi). Функциональная зави |
||||||||||
|
= |
х |
|
|
|
||||||||||
симость между Fs и Fj имеет следующий вид: |
|
|
(3.20) |
||||||||||||
|
|
|
|
|
|
|
|
|
Fs (1, |
|
s: Fi). |
|
|
128 Гл. II I . Формальное развитие теории кинетических уравнений
Уравнения (3.19) для I = 1, . . s запишутся так:
dt |
— — - ! ^ + — ( d 2 ë 12F2, |
(3.21) |
|||
т ds.i |
V |
J |
и f |
|
|
9FS |
|
s |
|
|
|
— LSFs+ ~ |
2 |
|
®г. s + l - f 1S + I - |
(3.22) |
|
dt |
|
i=l
Используя (3.20), производную по времени от Fs можно предста вить следующим образом:
д І \ |
/ б F4 |
dF, |
(3.23) |
|
I t ~ |
XöFi ’ |
ИГ |
||
|
где скобками обозначено внутреннее произведение по всем функ циям, a 8FS/8F! — функциональная производная от Fs по F^ В простейшем случае, когда Fs является алгебраической функ цией от Flf эта производная запишется в виде
9FS _ |
-у |
dFg |
dF1 (I) |
dt |
Гл |
dFi (l) |
dt |
|
1=1 |
|
|
Исходными уравнениями анализа Боголюбова являются послед ние три уравнения (3.21) — (3.23). Формальное разложение, кото рое будет использовано для получения замкнутой системы уравне ний, задается равенством
Fs( 1, |
JP1) = ^ 0) + 4 ^ 1,+ . . . , |
(3.24) |
F! = Fi (Fi не разлагается, поскольку она является «независимой переменной»).
Подставляя (3.24) в (3.21), получим
dFj |
»L ~ Ч |
|
— |
( d2ëiz |
Г/1' |
. P<i> |
4- |
||||
dt |
1 |
||||||||||
т |
öxj |
|
V |
J |
л |
L 2 |
1 2 |
||||
|
Н(0, + — А ( |
|
|
|
|
|
|
|
(3.25) |
||
|
' |
V |
|
|
|
|
|
|
|
|
|
Это соотношение |
служит |
для определения коэффициентов ИФ; |
|||||||||
|
|
Ат |
|
|
|
Рі |
d F i |
|
■LiF1> |
|
|
|
|
|
|
|
т |
дхі |
|
|
|||
|
|
И(1) = |
+ |
j |
â2QaF ? \ |
|
(3.26) |
A(l)= + j d2ëlzF t i)-
3.2. Анализ Боголюбова ББНГИ-уравнений |
129 |
Если мы подставим разложение (3.25) для dFJdt в функциональ ное внутреннее произведение (3.23), то в результате будем иметь
9F7, |
/ bFs |
, |
А 1' т А " ’+ - - ) = |
|
dt |
\6 F i |
|||
/ бö[Fj<»Р |
+ (i/v) FSP+ . . . ] |
А™ +±-А™ + . . . у . (3.27) |
||
|
бFi |
|
, |
|
|
|
|
|
Путем перемножения членов, получим
* • > + * { < » ? . ~ > + ф . * • »
(3.28)
Это разложение можно записать в следующем виде:
dFs |
: F>(0>F '0>+ |
— [DWF ^ + |
F>(1)F '0>] |
(3.29) |
||
dt |
|
' |
V |
|
|
|
где операторы |
и |
Z)ll) |
определяются равенствами: |
|||
>> = |
/ |
|
5 |
f)(l) = |
/_б |
(3.30) |
Д" |
W |
i ’ |
/ ’ |
_ |
W i |
|
|
>■ |
Вернемся теперь к оставшимся уравнениям ББКГИ (3.22) и ис пользуем в них разложение (3.24). В результате получим
^ = - L , [ f i » + 4 - W + . . . ] +
2 |
j d ( s + i ) é {. e+1[/'i¥1+ 4 ^ i + - - - ] . |
(3.31.) |
г = |
1 |
|
Приравнивая коэффициенты при одинаковых степенях к-1 в раз ложениях (3.29) и (3.31), придем к замкнутой системе функци ональных уравнений для переменных F® [1, . . ., s; FJ:
b imF(am= — LsF™, |
(3.32а) |
|||
|
|
S |
|
|
D™F? + D™F? = - |
L,F?' + |
2 j |
d (s + 1) Ѳ*. S+1F'°>, |
|
|
I |
І = 1 |
|
|
n |
• |
s |
|
|
2 b mF ? - l) = - LsF™ + 2 \ d (s + 1) Ѳ,. s + t F ^ . |
(3.326) |
|||
J=0 |
i=l J |
|
|
|
Последовательный ход |
решения |
этих |
уравнений представлен |
|
на рис. 3.4. |
|
|
|
|
Резюмируя изложенное, отметим три этапа анализа Боголю |
||||
бова: 1) предполагается, |
что FS= FS [F^; 2) находится уравнение |
для Fs [F4]; 3) построение завершается выводом двух разложений для dFJdt. Ими являются уравнения (3.29) и (3.31).
9-01243
130 Гл. II I . Формальное развитие теории кинетических уравнений
В данном формализме мы вывели замкнутую систему уравне ний для функций Fs, но две задачи остались нерешенными. Это, во-первых, получение конкретных выражений для операторов
D (t) и, во-вторых, метод вычисления Fi.
Именно на этих двух проблемах мы сосредоточим сейчас наше внимание. В частности, дадим стандартный метод, с помощью
Р и с . 3.4. Метод решения системы уравнений (3.32).
которого получается уравнение для Fi, справедливое до членов порядка ѵ~г. Для этого используем замкнутую систему уравнений, получаемую из (3.25) и (3.32а),
9 F I |
I |
P i |
d F j |
1 |
j d2è12*T, |
(3.25') |
dt |
' |
т |
дХі |
V |
||
|
|
DmF ^ = |
|
|
(3.32a') |
Первый шаг состоит в том, чтобы решить второе из этих уравне ний для
F™ =F™( 1, 2; Fi)
при некоторых определенных «граничных условиях». Это дает явную функциональную зависимость F20>от Ft при соответствую щим образом выбранном начале отсчета времени. Последнее озна чает, что система частиц (со взаимодействиями, препятствующими существованию связанных состояний) становится некоррелиро ванной, когда протечет достаточно длительный интервал времени в обратном направлении. Уравнение (3.32а) решается именно при этом условии. Такое условие поставлено в связи с его дальней шим приложением к процессам столкновения и означает, что до столкновения частицы не коррелированъ!.
На втором шаге полученное таким образом решение для F20> подставляется в уравнение (3.25') и выводится соответствующее
выражение для Ѳ12ІД0>, что приводит к характерной форме боголюбовского обобщения уравнения Больцмана.
3.2. Анализ Боголюбова БВКГИ-уравнений |
131 |
Для решения этих задач вводится оператор сдвига |
|
= exp [tLi], |
(3.33) |
где Li — ^-частичный оператор Лиувилля (cp. с Â;, |
который |
задается уравнением (2.30)). Определение L t вытекает из уравне ния Лиувилля для изолированной системы из I частиц:
Ц г = \Ни Fü ^ - U F l |
(3.34) |
Решение этого уравнения имеет вид
Ft ( 1, |
I; 0 = А(Ѵ г(1, |
J; 0). |
(3.35) |
Оператор ДІ?Т действует на фазовые переменные Fi- Он «сдвигает» их обратно по времени на интервал т по отношению к значениям, определяемым Z-частичным гамильтонианом і / г. Для единичной свободной частицы мы получим
 - \ F 1 (х, р, 0 = ^ і ( х —-|;Т , р, tj . |
(3.36) |
В более общем случае (возвращаясь к обозначению фазовой переменной через z) для произвольной фазовой функции ф мы имеем
Â(t4>(zi > . . . , Z b г) = ф(Д(тг)2і, . . . , Д ^ 2 г; t). |
(3.37) |
Из (3.35) следует, что Ft постоянна вдоль траекторий системы в Г- пространстве І частиц. Ниже приведено несколько других свойств
оператора А:
|
До=1, |
(3.38) |
|
д 11 ^ т 2 = Ä ( T l + T 2)> |
(3.39) |
ад-т |
— LA_t = — Д-XL. |
(3.40) |
дх |
|
|
Последнее соотношение справедливо для операторов, не завися щих от времени. Решением неоднородного уравнения Лиувилля
дР (z, t) |
LF (z, t) = P (z, t) |
(3.41) |
|
dt |
|||
|
|
||
является функция |
t |
|
|
|
(3.42) |
||
F (z, t) = Â_tF (z, 0) + A_, j' dt'hfP (z, t') |
(здесь использовалось свойство (3.40)).
Упомянутое выше граничное условие для «-частичной функции распределения как функционала от Fi легко задать через опера
9*
132 Гл. III . Формальное развитие теории кинетических уравнений
торы Â. Как мы помним, это условие состоит в том, что система s взаимодействующих частиц в далеком прошлом была некор
релированной, т. е. |
, |
|
lim Â%FS[1, |
.. ., *; F, (1, t)] = lim Д« |
ft Fi (h <)• (3-43) |
Т-fc o |
T—У оо |
1 = 1 |
Поскольку данное условие справедливо для произвольной функ
ции Ft, оно будет справедливо и для F't = lim At^ ^ i (1, t).
T-»- со
Подставляя |
эту функцию в предшествующее равенство, получим |
|
|
|
S |
lim Д(1\F S[1, |
.. ., s; А(Л Л (1, t)) = lim |
П А(Л Ft (l, t). (3.44) |
Т-Юо |
T—>oo |
1 = 1 |
Такая форма граничных условий используется при решении функ ционального уравнения (3.32а). Явный вид правой части (ПЧ) равенства (3.44) легко получить, используя (3.36). Это дает
ПЧ (3.44) = lim Â(_?T ft F A ( x t + ^ - x ) , рЛ =
T-*oo |
г = 1 |
L ' |
H l |
I |
J |
— lim Ц Fi [ Â ^ (хг+ - ^ - т ) , А ^рЛ =
T-»oo |
L |
' |
i n |
! |
J |
= II Л [ х Л , р Л ]. (3 .45)
г=і
Переменная pls> — это постоянный импульс, который имеет 1-я частица, взаимодействующая с совокупностью s частиц, если про следить за этой частицей из состояния р г в далекое прошлое. Переменная x';S) — это координата, которую имела бы l-я частица, если бы она двигалась с импульсом p)S) при начальной координате
Ali'ooX;. Такая интерпретация обусловлена той операцией, которая определена последним равенством:
X;S) = lim Â(_s)x ( Щ+ ~ т ) = Н т ( А(Д хг + |
х ) . |
Как мы увидим, эта интерпретация весьма полезна при нахож дении явного выражения для в уравнении (3.25').
Перейдем теперь к интегрированию уравнения (3.32а):
DmF ^ = —L»F?\ |
|
|
|
Согласно (3.26) и (3.30), |
операция Dl )F f' |
определяется |
следую |
щим образом: |
|
|
|
DtmF™ = |
Л<°>(Fi) = |
( - LiFi). |
(3.46) |

3.2. Анализ Боголюбова Б Б К ГИ-уравнений |
133 |
Эту функциональную операцию можно привести к виду, удобному для дальнейшего исследования и в то же время позволяющему включить предшествующее граничное условие, если воспользо ваться следующим приемом. Поскольку (3.46) справедливо для
любой функции F1: оно будет справедливо также и для AliiFt. Таким образом,
DWF?' [1, . . s; А Э Д ]= - |
р — Г- U |
• (3.47) |
о |
[Л- т^і] L |
J |
Вспоминая свойство (3.40), мы можем записать это уравнение так:
/)«»/?'<». |
бFT dbll\Fi |
(3.48) |
|
' ÖIÂ-tFj] |
’ |
что дает нам желаемый вид функциональной операции:
Z)(0)F ‘0,[1, . . . , 8; ДЗД] - |
дх F T [ 1, |
ATXF, |
(3.49) |
|
Подставляя это выражение в (3.32а), получим |
|
|
||
F T [1, . . . , |
s; Д -’т-Fі] = |
- L SFT [1, . .., |
s; Д З Д . |
(3-50) |
Рассмотрим пробное решение |
|
|
|
|
F T [1 |
Д-т^і] = |
№ xF ? (1......... |
Л ). |
(3.51) |
Подстановка его в (3.50) дает
± № xF T (1, . .., в; Fd = - Ls№ xF T (1, - • •, s; F,). (3.51a)
Если мы вспомним (3.40) и заметим, что FT в таком виде, как она записана в уравнении (3.51а), не зависит от т, то увидим, что это уравнение является тождеством. Отсюда следует, что решением уравнения (3.50), а следовательно, и уравнения (3.32а) будет про извольный функционал от функционала, заданного соотношением
(3.51).
Другой формой решения может быть следующая:
F T (1, .. . , s; Fi) = Д( Ѵі® (1, . ■■, s; Дт’Еі). |
(3.52) |
Это решение пригодно для всех интервалов т, так что
F T ( 1, . .., s; ^ ) = И т A<Ѵі°>(1, ■■■, s; Д^ ) (3.53)
T->co
также является решением. Используя предшествующее граничное условие, а также (3.44) и (3.45), получим
F T (U . . . ,s ; E i ) = [I iM x(A Pis)). |
(3.54) |
i=i |
|
134 Гл. II I . Формальное развитие теории кинетических уравнений
Это решение удовлетворяет уравнению (3.32а) и граничному усло вию (3.44). Таким образом, первая часть задачи о получении одно го уравнения для Ft из уравнений (3.25') и (3.32а') завершена. Из (3.54) видно, что решение низшего порядка Р$т в анализе Бого любова обнаруживает s-частичную корреляцию между частицами (через переменные x(S), p(S)).
Следующий шаг состоит в том, чтобы получить в уравнении
(3.25') соответствующее выражение для 6>i2F'0). С этой |
целью |
||||
напомним сначала соотношение (2.155) и последующие: |
|
||||
Ь2 = К 2- & 12 ) |
|
||||
к 2 = Рі |
д |
Р2 |
д |
|
|
Зх-і |
т |
д %2 |
|
||
|
|
||||
Отсюда вытекает, что |
|
|
|
|
|
e 12F<°> = K2F‘°> -L 2F ‘0’. |
(3.55) |
||||
Чтобы получить надлежащее выражение для L2F'0), вернемся |
|||||
опять к уравнениям (3.32а') и (3.26): |
|
|
|||
- L 2F™ = D"'F™ = ( ^ |
- , Л«”) , |
(3.56) |
|||
/ісо) _____Рг_ |
dFi . |
ч |
|||
|
Ат ' дхі
Выполняя функциональное суммирование, соответствующее внут реннему произведению ( , ), и вспоминая решение F '0>, которое
задается посредством (3.54), получим явный вид для L 2Fl2m:
F2F r |
р'2> |
Р22)] * д а Л К 2), p H - |
|
|
(3.57) |
+ - |
Fi [х)2), р l^ — F, [x'2\ p H - |
Наконец, чтобы получить нормальную форму обобщенного
уравнения Больцмана, произведем в операторе К 2 замену перемен ных Х[ и х2 на Хі и г = х2 — х4, что дает
К = {Р2~ Рі ) .JL4_.P1__ |
К2, Б ■ к 2. |
(3.58) |
2 \ т ) д г ' m öxi |
|
|
Используя выражение (3.57) для L2F'0) и (3.58) для K 2F f \ под ставим (3.55) в (3.25') и придем к обобщенному уравнению Больц мана в форме Боголюбова:
dFj |
Pi |
djj |
V j d2Kz, bF ‘0) + І - j d2 (K2- L 2) F ‘0). (3.59) |
dt |
m |
Ö X j |
В случае, когда F ‘0> однородна по всей области взаимодействия, второй интеграл обращается в нуль (К2 и Ь2 оба содержат про
3.2. Анализ Боголюбова ББКГИ-уравнений |
135 |
странственные производные), в то время как К 2>Б переходит
вобычный больцмановский оператор столкновений. Чтобы это показать, в соотношении
J К 2. бF T = j dx2dp2V - ^ F {F1[р)2>] F, [р«’]} (3.60)
(где аргумент х в F\ выпадает, поскольку F^ однородна по всей области столкновения) выполним интегрирование по х2 в цилин дрических координатах с осью цилиндра г), параллельной отно сительной скорости
|
|
|
|
V = -?2~ Pi |
; т] = - ^ . |
|
(3.61) |
|||
|
|
|
|
|
|
т |
1 |
V |
|
N |
Радиус цилиндра Ъи азимут ф лежат в плоскости, |
нормальной к У. |
|||||||||
Отсюда следует, |
что |
|
|
|
|
|
|
|||
|
|
|
|
dx2= Ьсіф db dp, |
|
|
||||
|
|
|
|
|
|
= |
|
|
|
(3-62) |
Интегрируя сначала по т], получим |
|
|
|
|||||||
|
+°° |
|
|
|
|
|
|
|
|
|
|
|
j |
drj |
{F±[рП Ft [р™]} = Fi [pH Fi [pn + l - |
(3.63) |
|||||
|
—oo |
|
|
|
|
|
|
|
|
|
Области после столкновения соответствует ц = |
+оо (рис. 3.5). |
|||||||||
Если из этой области (ц = |
+ °°) проследить в бесконечно далекое |
|||||||||
прошлое, то каковы будут зна |
|
х2 |
|
|
||||||
чения рі и р2 |
в случае парного |
|
|
|
||||||
взаимодействия |
частиц? |
Ими |
|
|
|
|
||||
должны быть те значения им |
|
|
|
|
||||||
пульсов, которые перейдут |
в р4 |
|
|
|
|
|||||
и р2 после столкновения, т. е. |
|
|
|
|
||||||
импульсы р( |
и |
р' инверсивного |
|
|
|
|
||||
столкновения. |
(Эти понятия бу |
|
|
|
|
|||||
дут |
подробно |
обсуждаться |
в |
|
|
|
|
|||
гл. IV.) |
до столкновения со |
|
|
|
|
|||||
Области |
|
|
|
|
||||||
ответствует т] = |
—оо (рис. 3.5). |
|
|
|
|
|||||
Каковы будут значения рі и р2 |
|
|
|
|
||||||
при парном взаимодействии ча |
|
|
|
|
||||||
стиц, |
если из этой области (ц = |
|
|
|
|
|||||
= —оо) проследить в бесконеч |
|
|
|
|
||||||
ное прошлое, которое представ |
Р и с . |
3.5. Связь между знаком р и |
||||||||
ляет |
собой |
состояние «до |
до |
|
областью столкновения. |
столкновения», т. е. область, находящуюся далеко за пределами зоны взаимодействия двух ча
стиц? В этой области две частицы движутся независимо одна от
136 Гл. I l l . Формальное развитие теории кинетических уравнений
другой как свободные, поэтому их импульсы равны рі и р2. Сле довательно,
Fi feO Fl [рП |t~ = F, (p() F, (pj) - F1(Pl) Fl (Pa). (3.64)
Подставляя это соотношение в (3.63), а затем в (3.59) (и удержи
вая только однородный вклад К 2,в)> получим обычное уравнение Больцмана:
^ + |
^ I |
р Э - |
|
- Л ( Р і) ^ і (Р2)]- |
(3.65) |
Уравнение Больцмана подробно рассматривается в гл. IV, где представлены три других его вывода. Более прямой из них не опирается на последовательность ББКГИ-уравнений. Но весь ма примечательно, что анализ Боголюбова дает независимый (т. е. опирающийся только на уравнение Лиувилля) метод, кото рый позволяет получить больцмановский интеграл столкновений
вявном виде. В данный момент важно, что такой метод существует.
Вглаве IV будет дана всесторонняя трактовка уравнения Больц мана, которая обеспечит более глубокое его понимание.
3.3. Групповые разложения
Совсем недавно были предложены другие методы разложения для исследования ББКГИ-уравнений. Сюда относятся работы Ростокера и Розенблюта (1960), Гернсея (1960), Ленарда (1960), Балеску (1960), Сэндри (1963) и Дюпри (1961). Однако эти работы
впротивоположность предыдущим больше связаны с выводом определенного кинетического уравнения, чем с получением реше ния для Ѵ-частичной функции распределения. (Кинетическим называется любое уравнение в переменных (t, |, х) для неизвест ной функции fi.) Тем не менее все эти методы объединяет то, что
вкаждом из них решение ищется в виде некоторого разложения. В частности, метод Боголюбова связан с разложением по обратным степеням параметра ѵ — удельного объема.
Существуют другие способы разложения, непосредственно опирающиеся на свойства Fz, F3, . . ., F N в предельном случае, когда k-1-> 0 . Что это за свойства? Очевидно, что, когда объем, приходящийся на одну частицу, становится большим, а следова тельно, и расстояние между частицами также возрастает, тогда фаза какой-либо частицы становится менее связанной с фазой любой другой частицы. Это утверждение о независимости состоя ний, записанное через функцию распределения Fs, выглядит так:
^ = П Л ( 0 - |
(3.66) |
!=Т |
|