
книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений
.pdf3.3. Групповые разложения |
137 |
Для случая s = 2 |
|
F2 (1, 2) = F! (1) Fi (2). |
(3.67} |
Если состояния двух частиц не зависят одно от другого, то гово рят, что эти частицы не коррелированы. Мера влияния состояния одной частицы на состояние другой частицы (отличной от первой) определяется корреляцией между частицами. Равенство (3.67) справедливо в предельном случае, когда корреляция отсутствует. Отклонение функции F2 (1, 2) от факторизованной формы
Fi (1) Fi (2) называется корреляционной функцией с62. Таким образом,
F2 (1, 2) = Fi (1) Fi (2) + |
(1, 2). |
Если <ё 2 мала по сравнению с Fi (1) F2 (2), состояния частиц 1 и 2 являются независимыми. Ясно, что это будет иметь место в случае «большого» удельного объема ѵ при условии, что газ однородно заполняет конфигурационное пространство. Этот результат можно сформулировать более точно с помощью г0 — радиуса потенциала межмолекулярного взаимодействия ср и средней тепловой скоро
сти У кТ /т (к — постоянная Больцмана, Т — температура, m — масса частицы). Если ѵ — средний объем, приходящийся на одну частицу, то при ѵ, большом по сравнению с г„, каждая из частиц «не подозревает» о существовании другой. В этом случае соответ ствующим малым параметром будет Ei — (г®/у).
Классический механизм взаимодействия между частицами осу ществляется через потенциал <р. Если температура такова, что кинетическая энергия значительно превышает потенциальную, то последняя не будет играть роли в изменении фаз рассеиваю щихся частиц. В этом случае малым параметром является е2 = = сро/кТо, где фо — «напряженность» потенциала.
Существование в различных предельных случаях малых корре ляций между частицами было использовано Сэндри для приведе ния ББКГИ-уравнений к безразмерному виду. Получающаяся в результате теория одновременно и формальна и иллюстративна. Введем сначала обозначение, в котором последовательность ББКГИ-уравнений (3.19) примет следующий вид:
dt -j- L$FS= n0k sFs+1. |
(3.68) |
Здесь щ — ѵ 1 и оператор Âs определяется равенством
(3.69)
138 Гл. II I . Формальное развитие теории кинетических уравнений
Кроме того, напомним, |
что |
|
|
|
|
|
|
|
L |
|
Ks~&s, |
|
|
k s = |
У, |
-PL. |
è r |
Ѳ - - 2 2 Ѳ » . |
|
|
|
1=1 |
m |
|
|||
|
|
|
|
i<j |
(3.70) |
|
Л |
д |
|
д |
д |
д |
|
|
|
|||||
&іі |
^ дхі фгі |
dpі "г" dxj фг-' |
<Эр^ |
|
||
Запишем систему ББКГИ-уравнений (3.68) через эти новые |
опе |
|||||
раторы: |
|
|
|
|
|
|
|
dF$ |
'-KsFs — ë sFa= n 0k sFs+l. |
(3.71) |
|||
|
dt |
|
|
|
|
|
При такой записи уравнений выявляется только действие параметра п0. Однако рассматривать предельный случай большого размерного удельного объема физически бессмысленно. Корректно только лишь сравнение ѵ с некоторым известным постоянным объемом.
Для того чтобы устранить все размерные параметры в урав нении (3.71), мы должны привести к безразмерному виду все переменные, входящие в это уравнение. С этой целью введем без размерные (отмечаемые значком ~ ) переменные:
х = г0х,
|
Р = |
РоР, |
|
|
|
фи = фофгі |
( либ° GiJ = |
Т5- GU ) > |
|
||
|
t — — t |
|
(3.72) |
||
|
г? |
p |
|
||
|
1 s |
|
|
||
где |
S~ P 30S ’ |
|
|
||
mC, |
mC2 = |
JcT0. |
(3.73) |
||
Po = |
Расстояние r0 — это радиус потенциала взаимодействия. Определе
ние безразмерной силы G и потенциала ср обеспечивает согласован ность введенных переменных.
Подставляя эти новые переменные в уравнение (3.71), получим
( I Ц = ( « л Ш â â « - <3 -74)
Опуская значки «тильда» (все переменные и операторы являются безразмерными), придем к ББКГИ-уравнениям в форме Сэндри:
( + K s — a ë s) Fs = aßÂsF s+I,
а = |
Фо |
(3.75) |
п0г30. |
||
|
кТп |
|
|
3.3. Групповые разложения |
139 |
||||
|
|
|
а ^ fJ k T |
|
|
|
Сильное |
ёзаимодейстдие; |
|
|
|
|
|
низнал |
температура |
€-і |
(<+f) |
(«-1, о |
(<+€+> |
|
|
|
|
||||
|
|
1 0 ,0 |
(і,і) |
а,« -1) |
|
|
Слабое |
ВзаимодейстВие-, |
|
|
|
|
|
Высокая |
температура |
( («, 0 |
(е,1) |
(€,«-4 |
|
|
|
|
|
е |
1 |
ДальнодейстВу- |
|
|
|
|
КороткодейстВу- |
|||
|
|
|
ющие силы; |
|
юицие силы, |
|
|
|
|
разреженная |
|
плотная среда |
|
|
|
|
среда |
|
|
|
Р и с . 3.6. Упорядочение |
параметров (а, |
ß) в групповых |
разложениях. |
Эти уравнения имеют два характерных физических параметра: q>JkTQ и п3г\. Область, представляющая физический интерес, полностью исчерпывается, если позволить этим параметрам при нимать значения, равные по порядку величинам е-1, 1 и е (г — бесконечно малая). Возникающие в результате девять возможных случаев представлены на рис. 3.6.
Задача 3.2. Что можно сказать относительно |
если |
= О, |
|||||
и относительно |
если %3 = 0? |
Описать словами |
«физику» |
||||
каждой ситуации. |
|
|
|
|
|
|
|
Корреляции более высокого порядка определяются по следую |
|||||||
щей схеме: |
|
|
|
|
|
|
|
F2(1, |
2) = F1{1)F2(2) + |
« 2(1, 2), |
|
|
|
|
|
F 3 (1,2,3) = F i |
(1) F i (2) F i |
(3) + F y (1) <ë2(2, 3) + |
F t (2) %2(3,1) + |
||||
|
+ |
Л ( 3 ) а д , |
2) + «a(l, 2, 3), |
|
|
|
|
F,( 1, 2, 3, |
4) = Fi(l) Fi(2) Fi(3) Fi(A) + F i(l) Fi(2)%2(3, 4) + |
||||||
|
+ Fi(2)Fi(3)%z (A, i)-\-Fi (3) Fi (4) %2(1, |
2) + |
|||||
|
+ ^ (4 ) ^ ( 1 ) ^ (2,3)+ ^ ( 1 ) ^ ( 3 ) |
(4,2)+ |
(3.76) |
||||
|
+ |
Fi (2) Fi (4) «2 (3,1) + Fi (1) <ë3(2, 3, 4) + |
|
||||
|
+ + i(2)^з (3,4, l) + Fi (3)«з (4, 1,2) + |
|
|
||||
|
+ ^ ( 4 ) ^ ( 1 , |
2, 3) + « 2(l,2)«a(3, |
4) + |
|
|
||
|
+ |
<ë2(2, 3) <&2(4, 1) + |
(1, 3) «2 (2, 4) + |
|
|
||
|
+ |
^4(1, 2, 3, 4). |
|
|
|
|
140 Гл. II I . Формальное развитие теории кинетических уравнений
Эти формулы обладают следующим свойством:
s
Fs -> II Fl (1), при «г -> о, 2 < і < S.
і=і
Кроме того, если из данной группы частиц выделяется некото рая подгруппа, то функция распределения этой группы частиц соответствующим образом разлагается на множители. Предполо жим, например, что в группе из четырех частиц (1, 2, 3, 4) частицы 1 и 2 находятся каждая в области влияния другой, но далеки от частиц 3 и 4. Каждая из последних двух частиц, подобно первым двум, находится в области влияния другой. Такое состояние характеризуется следующими соотношениями:
F4 (1, 2, 3, |
4) |
= F2 (1, |
2) F2 (3, |
4), |
||
|
^4 = |
0, |
= |
0, |
|
|
(1, 4) = |
(2, |
3) |
= |
<ёг (1, |
3) = |
(2, 4) = 0. |
Подставляя эти выражения в последнее из равенств (3.76), полу чим
Fk{i, 2, 3, 4) = |
Л(1)Л(2)Л(3)Л(4) + Л (1)Л (2)^2(3, 4) + |
|
||
|
+ Fi |
(2>)Fі(4)^2(1, 2) + <ё2(і, 2)«2(3, 4) = |
|
|
= |
ІЛ(1) |
Л(2) + |
%(1, 2)Н/’1(3)Л(4) + %2(3, 4)] |
= |
- |
F2(1, |
2)F2(3, |
4). |
(3.77) |
Рассмотрим также другое часто встречающееся представление корреляционных функций. Оно соответствует так называемому разложению Кирквуда, хорошо известному химикам-теоретикам. Корреляционные функции в этой схеме обозначаются через
%(1, 2), h (1, 2, 3), А (1, 2, 3, 4) и т. д. Они определяются сле дующим образом:
F2 (1,2) = Fi (1)F i (2) 11 + * (1,2)1,
^з(1, 2 ,3 )= F 1(\)Fl(2)Fi(3) [1 + f (i, 2 )+ ?(2,3)+?(3,4)+4(1, 2,3)].
Предположим, что рассматриваемая система бесконечно близка к состоянию с нулевыми корреляциями. В таком случае «все I- частичные корреляционные функции имеют порядок е и выше, так что
F s = f l Fi ( I) + 0 (е) = F™ + eF<15 + . . . , |
(3 .78) |
і=і |
|
где Fs0) записывается в форме произведения (3.66). |
(разрежен |
Поскольку F101 описывает равновесное состояние |
ного газа), то представленное выше разложение справедливо для
3.3. Групповые разложения |
141 |
состояний, близких к равновесию. Рассмотрим теперь некоторые из предельных случаев, изображенных на рис. 3.6.
Случай а: ß = О (г), а = О (1).
Если положить ß = п0г30 = е, то система уравнений (3.75) сведется к следующей (удержим только члены наинизшего порядка по г):
^ + |
= 0, |
(3.79а) |
( 4 + Я2) (В Д ) - |
«Ѳ2е д = 0, |
(3.796) |
( ± . + к , - а Ѳ а) Ц Fi(t) = 0. |
(3.79л) |
Имеем N уравнений для одной неизвестной функции Fr. В этом случае нетривиальное решение возможно, только если уравне ния совместны. То же самое справедливо относительно всех раз ложений ББКГИ-системы около нулевых корреляций. Данное условие существования нетривиального решения следует из пред положения о форме ведущего члена в разложениях (3.78), соглас-
S
но которому F™ = [J Fi (I). Совместность уравнений нулевого і=і
порядка является необходимым условием факторизации Es0’.
Задача 3.3. Установить, являются ли совместными уравне
ния (3.79). |
|
|
|
Случай б: а = О (s), ß ^ |
О (1). |
и удерживая в (3.75) |
|
Полагая а = О (е) в разложении (3.78) |
|||
только члены наинизшего порядка по 8, |
придем к |
следующим |
|
уравнениям: |
|
|
|
( ± |
+ É . ) g F l{l>= 0. |
(3.80) |
Эта согласованная система уравнений имеет решение в виде свободно-молекулярного течения:
F i { x , Р, £) = ф [ х — ( ^ - ) с р] * |
(3-81) |
где ф — любая функция от х — (р/m) і и р. Таким образом, реше ния уравнений (3.80) постоянны вдоль траекторий свободных частиц р = ро и X — (p/те) t = х0.
Полученный результат имеет силу, когда ß ^ О (1), так что решение пригодно для областей (8, е) и (е, 1) (рис. 3.6). Следова
142 Гл. I I I . Формальное развитие теории кинетических уравнений
тельно, г„ может быть порядка гс0_1/3 при условии, что напряжен ность потенциала много меньше кинетической энергии кТ 0. Пер вые два уравнения системы (3.75) до членов О (е) включительно имеют вид:
( ± |
+ К 1) F, (1) = eccßÂ^ (1) F, (2), |
(3.82а) |
|
( 4 t + K z)[F i (1) F, (2) + |
£%>] - еаѲ2Л (1) Fx(2) = |
|
|
|
= eaßÂ2F t (1) Fi (2) Fx(3). |
(3.826) |
|
Для s = 3 будем иметь |
|
|
|
( ~ |
) (F1FiFl + |
е З Д + ... -f е^з) - |
|
— eaQ3FiFiFi = |
eaßA3FlF1FiF1. |
(3.82в) |
Получающуюся в результате систему уравнений в принципе очень легко решить. Во-первых, ІД получается из (3.82а). Когда Ft найдена, оставшиеся уравнения все являются неоднородными уравнениями следующего вида:
+ |
4ßt>). |
(3.83) |
При таком методе решения функции |
распределения |
находятся |
в соответствии со следующей схемой (до членов О (е) включительно):
Я |
|
а д л |
Яа1»/ |
F1F1F1F1 |
Fl1 |
Полезно, между прочим, отметить одно важное свойство корреля ционных функций %і (1,2, . . ., I). Рассмотрим уравнение, опре деляющее
F2 (1, 2) = F, (1) Fi (2) + «2 (1, 2). |
(3.84) |
Но функция Fi в свою очередь определяется через F2-
j ^ 2(1» 2) d2 — VFi(i)\ |
j Fi(2)d2 = V. |
(3.85) |
Сравнивая последние два уравнения, заключаем, что
j « 2(1, 2)d2= j « 2(1, 2)di = 0. |
(3.86) |
|
3.3. Групповые разложения |
|
143 |
||
Для 'ёз мы находим |
|
|
|
|
|
j dl d3 [F, (2) « г (1 , 3) + Fi (3) « 2 (1, 2) + |
s (1, 2, 3)] = 0, |
||||
j di d3 [Ft (1) |
(2, 3) + |
Fi (3) <ёг (1, 2) + |
(1, 2, 3)] = |
0, |
|
dl dl [Fi (1) <ë2(2, 3) + |
Fi (2) <вг (3, |
1) + <ë3 (1, 2, 3)] = |
(3.87) |
||
0 |
|||||
или, используя соотношение (3.86), |
f dl dl<ë3 = 0. |
|
|||
j dld3F>a = |
j dl d3%a= |
(3.88) |
Аналогичные утверждения справедливы для %ѣ, s > 3 . Эти усло вия могут служить дополнительными уравнениями в любом данном методе разложения. Они определяют ограничения, нало женные на корреляционные функции. Так, например, из усло вия (3.86) следует, что не может быть четной функцией по всем из своих аргументов.
Задача 3.4. |
Доказать, |
что ^ cës dk = 0, |
к = 1, |
2, . .., s. |
Случай в: а = |
О (е), ß = |
О (е-1). |
введен |
Розенблютом |
Этот предельный случай впервые был |
и Ростокером (1960) при исследовании соответствующих урав нений для высокотемпературной плазмы. Плазмой называется состояние вещества за пределами атомарной газообразной фазы. Когда атомарный газ нагревается выше некоторого предела, появляется механизм, который непрерывно поглощает кинетиче скую энергию хаотического движения. Это подобно энергии диссоциации в жидкости с полукристаллической структурой, кото рая является как бы каналом для непрерывного поглощения энергии при переходе из жидкого состояния в газообразное. При переходе в плазменную фазу поглощающим каналом является механизм ионизации, и в новой фазе частицами будут положи тельные ионы и отрицательные электроны. Законы, управляющие поведением плазмы, так же сильно отличаются от законов газо образной фазы, как законы газообразной фазы отличаются от зако нов жидкой фазы. Силы взаимодействия в плазме — кулоновские. Они отличаются от всех других сил взаимодействия между части цами очень большим радиусом действия. Эти две характеристики — дальнодействие и высокая температура — дают возможность пред
положить, что ч>о/кТо ~ е и |
г30п 0 ~ |
е-1, так что (роГІп0ІкТ0 ~ 1. |
|
Для кулоновского |
закона |
взаимодействия |
|
|
|
|
(3.89) |
где е — заряд каждой |
частицы, а |
е0 — диэлектрическая прони |
цаемость вакуума. Радиус г0 и напряженность ф0 взаимодействия
144 Гл. II I . Формальное развитие теории кинетических уравнений
■связаны соотношением
Фо |
w |
r |
■ |
|
(3'90> |
Го |
|
||||
В силу того что (p0rln0/kT 0 ~ |
1 , мы можем оценить соответствую- |
||||
щий радиус взаимодействия гр: |
£0кТ0 |
|
|
||
|
|
|
(3.91) |
||
|
|
е2по |
|
||
|
|
|
|
||
Характерная частота со7 определяется |
равенством |
|
|||
|
С \2 |
е2п0 |
(3.92) |
||
|
|
8 |
$П1 |
||
|
|
|
Расстояние гр — это эффективный радиус действия кулоновского потенциала в плазме. Более точные вычисления (впервые прове денные для электролитов Дебаем и Гюкелем (1923)) показали, что заряженная частица, помещенная в нейтральную плазму, вызы вает перераспределение заряда, стремящееся экранировать проб ную частицу, так что вне гр (радиус Дебая) кулоновский потенциал в е_1 (ln е = 1) раз меньше своего значения в вакууме.
Частота сор связана также с основным коллективным колеба нием, которое будет совершать возмущенная плазма, выведенная из состояния локальной зарядовой нейтральности. Установлено, что время образования около пробного заряда дебаевского экра нирующего облака имеет порядок сор1.
Ознакомившись с самыми элементарными понятиями физики плазмы, вернемся теперь к разложению ББКГИ-уравнений, когда
а ~ е и ß ~ 8 _1, |
и к |
концепции |
предельного случая нулевых |
|
корреляций между частицами. |
До настоящего момента мы просто |
|||
считали, что надо положить |
~ |
= О (1)), чтобы формаль |
||
но описать этот предел. |
Однако мы получим более последователь |
|||
ный формализм, если положим |
~ sh~1<ë k. В этом случае |
|||
имеет порядок е3, |
имеет порядок е и т. д. Но обосновано ли физи |
чески такое упорядочение? Иными словами, если 'ëz мала, то будет ли «з еще меньше?
Хотя точное доказательство этого факта получить трудно, некоторые простые оценки, по-видимому, подтверждают правиль ность сделанного нами предположения. Их можно получить,
вычисляя объем |
допустимых |
областей фазового пространства, |
|
где |
отлична |
от нуля *). |
Рассмотрим элементарную задачу |
о трех шариках на прямой проволоке длины L. Областью взаимо действия между любыми двумя из них является г0 (в данном слу
чае это диаметр |
шариков). Чтобы вычислить |
(1, —), зафикси |
руем шарик 1 в |
начале координат, тем самым мы определим пло |
1) Можно рассматривать 4gh как ^-частичное распределение. В этом случае мера равна отношению объема в Г-пространстве, который включает состояния с ^--частичными корреляциями, ко всему объему Г-пространства.
3.3. Групповые разложения |
145 |
скость в трехмерной конфигурационной части пространства Гв. Области, в которой %2 (1, —) отлична от нуля, соответствует заштрихованная площадь на рис. 3.7, а. Когда шарик 2 находится в интервале —г0 < х 2 < + г0, то при ré 2 (1, —) фО шарик 3 может находиться где угодно на прямой. То же самое справед
ливо, |
если шарики 3 |
и 2 поменять местами. Отношение А 2 пло |
|||
щади, |
для |
которой |
%2 (1, —) ^4=0, к |
полной площади равно |
|
(4rJL) |
[1 — |
(r0/L)] для (r0/L) -< |
и равно единице в любом другом |
||
случае. Чтобы определить, где |
(1, —, |
—) фО, поместим опять |
|||
шарик |
1 в точку (0, |
0, 0). Состояниям |
с ненулевой функцией |
|
(1, |
—, —) соответствует заштрихованная область на рис. 3.7, б. |
||
Относительная площадь, на которой %3 (1, —, |
—) конечна, равна |
|||
А |
3 = |
8Гд/L2. Если |
мы обозначим r0/L=== е, то |
в пределе, когда |
£ |
0 , |
А 3 = 2гА2. |
В более общем случае А 2 ~ А ^ 2 (1 — АV2). |
|
|
Наши вычислениясильно упростились, поскольку мы прене |
брегли кинетической энергией частиц. Вообще говоря, для любого заданного прицельного параметра (в классическом случае) дейст вие рассеивающего потенциала сводится к нулю, когда энергия падающей частицы становится достаточно велика. Можно утвер
ждать, что в пределе при ср0ІкТ0 |
0 эффективная область взаимо |
действия уменьшается. С другой |
стороны, когда ц>0/кТ0 оо, |
эффективный радиус взаимодействия становится равным г0. Про стейший способ включения этих соображений в предшествующий анализ состоит в замене г0 на г0 [1 — ехр (—ф0/&Го)1; тогда А 2 ~ ~ Л (1 — Ф), А з ~ .42, где
= ( т ) [ і —exp (-= J^ )] •
10-01243
146 Гл. II I . Формальное развитие теории кинетических уравнений
Если и сравнительно однородны в областях, где они отличныі
от нуля, |
то приведенные выше доводы дают возможность пред |
|||
положить, |
что |
|
||
(где |
ß ~ |
r0/L и а ~ ц>0/кТ0). В области высоких |
температур |
|
(а ~ |
О (е)) |
эта зависимость сводится к следующей: |
|
|
|
|
|
^ ~ T = i ß ^ ’ |
(3-94> |
вто время как в области низких температур
икорреляции более высокого порядка могут быть значительными.
Это простое эвристическое рассуждение показывает, что в неко торых предельных случаях корреляционные функции можно упо рядочить по величине в соответствии со степенью рассматриваемой корреляции. Пределы задаются посредством введенных ранее коэффициентов а и р . Интересно отметить, что даже в таком очень простом представлении мы встречаемся с трудностью, которая появляется в сингулярном пределе а ~ е, ß ~ в-1 (высокая температура, дальнодействие), когда упорядочение согласно сте пени корреляции становится сомнительным.
Мы пытались показать, что упорядочение по относительным величинам корреляционных функций в соответствии со степенью корреляции подсказывается некоторыми простыми физическими соображениями. Соглашаясь с этим, мы формально приходим к следующей|последовательности корреляционных функций:
= [г<ё? + г Ч ? '+ ...], |
|
|
|
%3 = е[еЩ1>+ £г'ё(3а'+ |
...] , |
, |
(3.95) |
« 4==е2[е^1) + 82^ 2,+ |
...]. |
|
|
Полагая а ->- еа, ß -> e-1ß и используя (3.95) и (3.76) в урав нениях (3.75), получим следующую систему уравнений (до чле нов О (в)):
( 4 - |
+ ^ i ) ^ i = ctßÂ1[E1(l)F 1(2) + 8 ^ 1)(l, 2)1, |
(3.96а) |
||||||
( 4 і + К |
г |
~ |
(1) Fi (2) + |
e ^ 1’] = |
|
|
|
|
= |
aßÂ2 [Fi (1) Fi (2) Fi (3) |
+ eF, (1) |
(2, 3) |
+ |
|
|||
|
|
|
+ eFi (2) |
(3, |
1) + BF, (3) |
(1, |
2)1. |
(3.966) |