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

книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений

.pdf
Скачиваний:
36
Добавлен:
25.10.2023
Размер:
13.16 Mб
Скачать

3.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)

Соседние файлы в папке книги из ГПНТБ