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

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

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

4.2. Концепция столкновений

177

динат, связанной с центром масс (т. е. когда мы ведем наблюде­ ние, будучи «привязанными» к центру масс), показана на рис. 4.5.

Смысл понятий до и после столкновения определяется через потенциал взаимодействия Т (г) либо через радиус взаимодейст­ вия г0. До столкновения — это любой момент времени, когда частицы приближаются одна к другой и г > г0, или, что эквива­ лентно, (г) — 0. После столкновения — это любой момент времени, когда частицы уда­ ляются одна от другой и Т (г) =

=0 ( г > г0).

Относительная скорость г до

ипосле столкновения обладает очень важным свойством. Обо­ значим относительную скорость до и после столкновения через У. Так как относительный им­

пульс р равен рг, то из (4.3) с учетом того, что Н является кон­ стантой движения, следует ра­ венство

 

-іц F2 = - ip F '2,

(4.8а)

 

 

 

так

что

 

Р и с .

4.5. Ориентация

относитель­

 

V = Г ,

(4.86)

ного

радиуса-вектора г до и после

 

столкновения в системе

координат

где

V —величина относитель­

 

центра масс.

 

ной

скорости до столкновения,

 

 

 

а У'

— после столкновения. Таким образом, величина относитель­

ной скорости после завершения столкновения сохраняется. Она не является константой движения, а, вернее, принимает свое первоначальное значение после столкновения.

Задача 4.2. Доказать следующую теорему. Угол рассеяния в системе координат, связанной с центром масс, равен углу рас­ сеяния в системе координат, связанной с какой-либо из двух рас­ сеивающихся частиц. В системе координат, связанной с данной частицей, последняя находится в состоянии покоя.

Ответ. В системе координат, связанной с частицей, угол рассеяния Ѳр есть угол поворота вектора относительной скоро­ сти V при рассеянии, т. е.

Ѵ-Ѵ'

COS Ѳр = у2

Здесь мы использовали полученный

ранее результат: V = V '.

В системе координат, связанной с центром масс, угол рассеяния

12-01243

178

Гл. IV. Уравнения Больцмана,

Крука Бхатнагара Гросса и др.

Ѳс

— это угол поворота скорости

— R при столкновении, т. е.

 

 

cosѲс =

 

 

.

 

 

 

I S x - R I I S i - R ' l

 

(Отметим изменение в обозначениях:

= г^) Но

 

 

— R = £, _ ^НІ1±^2|2_ = .— rü2

(£ s \ =

----- tL у.

 

^51

т1+ т2

"4 + m2

VS1

 

Аналогично

 

 

 

 

 

 

 

i ; - R '= - J - v ' ,

 

так

что

 

 

ѴѴ'

 

 

 

 

COS Ѳс

 

 

 

 

 

 

р2 •

 

 

Сравнивая последнюю формулу с предшествующим выражением для Ѳр, получим желаемый результат:

Ѳс = ѳ р .

Соответствующие углы и векторы, использованные в данном доказательстве, показаны на рис. 4.5 и 4.6.

т і

Р и с. 4.6.

 

а — рассеяние в системе координат частицы m2; б

рассеяние в систе-

ме координат ЦМ.

 

В лабораторной системе координат уравнения сохранения импульса и энергии имеют вид

P i

+

Р 2

Р І +

Р а ,

(4.9a)

Рі

,

РІ

_

Р?

Pi

(4.96)

2mi

1 2m2

 

2mj

2m2

 

 

4.2. Концепция столкновений

179

Эти уравнения связывают р(, р' с рь р2. Однако для заданной совокупности начальных импульсов р4, р2 существует континуум значений pj, p', удовлетворяющих уравнениям сохранения. Это не является неожиданностью, поскольку уравнения сохранения не содержат информации относительно потенциала взаимодейст­ вия Т (г). Какую начальную информацию достаточно иметь, чтобы определить конечные импульсы pj, р2? Мы нашли, что звдач

Р и с . 4.7. Рассеяние в системе координат, связанной с относительным век­ тором г.

двух тел эквивалентна задаче одного тела с фиктивной массой ц и радиусом г. Случай рассеяния в системе координат, где начало, этого вектора фиксировано, изображен на рис. 4.7. Единичный апсидальный вектор а делит пополам угол между V и У'. При­ цельный параметр s таков, что угловой момент рассматриваемой совокупности частиц в данной системе координат (эквивалентный моменту в системе координат, связанной с центром масс) равен: [гEs. Прицельный параметр s, а также векторы относительной скорости У и V' характеризуют свойства асимптотических состоя­ ний при столкновениях. Модуль скорости Е и прицельный пара­ метр s сохраняются при данном столкновении, так что эти две. скалярные величины могут характеризовать данное столкновение^

Задача 4.3. Пусть в данный момент времени положение и импульс центра масс системы двух тел определяются векторами (R, Р). Показать, что угловой момент системы равен

L = R X Р + epsE,

где е — единичный вектор нормали к плоскости, в которой лежит вектор г ( t ) .

Как будет показано, р( и р'2 определяются через V' и постоян-. ный вектор Р. Вектор V' имеет величину Е, а его ориентация

12*.

180 Гл. IV. Уравнения Больцмана, Крука Бхатнагара Гросса и др.

I

задается углом рассеяния Ѳ (или углом ф = — (л — Ѳ)). Угол

рассеяния определяется как функционал от потенциала взаимо­ действия Т (г), относительной скорости V и прицельного пара­ метра s, т. е.

Ѳ = Ѳ (s, Т , V).

(4.10)

(Это также будет установлено позднее.) Теперь мы можем заклю­ чить, что для заданного взаимодействия Тіг) и импульса центра

масс Р конечные импульсы р(, р' определяются через s и р4, р2. Поэтому изменение импульсов при парном столкновении можно представить следующим образом:

І(Ри Ра) —►(РІ.

(4.11)

Индекс s указывает, что прицельный параметр есть свойство столкновения.

Существуют два типа столкновений, которым в теории кинети­

ческих

уравнений уделяется особое внимание. Это инверсив­

ные и

реверсивные столкновения.

Инверсивным по отношению

к [(рі,

р2)

(рі, Рг)]5 является столкновение, которое имеет конеч­

ным состоянием (р4, р2). Мы видим,

что если [(pt, р2) ->- (р(, р')]3

удовлетворяет уравнениям сохранения, то [(p', p') ->■ (р1? p2)]s также им удовлетворяет (следует только прочитать уравнения сохранения справа налево). Прицельный параметр s гарантирует, что (p', pj) правильно ориентированы, чтобы дать (рь р2). Отсюда следует, что начальными импульсами инверсивного столкновения являются (pj, Рз)3, т. е. конечные импульсы первоначального ■столкновения (рис. 4.8). Можно применить следующую символи-

4.2. Концепция

столкновений

181

ческую запись:

 

 

Инверсия [(Рі, р2) -*■ (Р;, р;)],

= [(р;, Р;)

(Pl, P2)]s. (4.12>

Второй важный тип столкновения — реверсивное столкновение. Реверсивным по отношению к [(Pl, р2) -ѵ (р;, P;)]s является

столкновение, которое имеет конечным состоянием (—Рі, —р2). Обращаясь опять к уравнениям сохранения, приходим к выводу,

Р и с. 4.9. Столкновение и его реверсия.

что начальными импульсами будут ( —рі, —p')s (рис. 4.9). То есть, если частицы с импульсами ( —Рі\ —рі) взаимодействуют с при­ цельным параметром s, то в результате взаимодействия получают­ ся ( —Рі, —p2)s. Снова можно использовать символическую запись:

Реверсия [(Рі, р2) -> (pi, Pi)]s= [(—РІ, —РІ) —►( —Pi, —p2)ls- (4.13)

Хотя эти утверждения выведены для общего случая, когда гПі ф т2, при изучении уравнения Больцмана для газа из одина­ ковых частиц достаточно рассмотреть только частный случай т1 = т2. Тогда уравнения сохранения импульса и энергии (4.9) будут зависеть только от скоростей І! и | 2 (Р = т0>- Они примут вид:

іі + І2 — 1і + І2,

(4.14)

 

Эти уравнения имеют решения:

 

Іі = 1і +сс(а-Ѵ),

(4.15)

li = 12 — a(a-V),

V ^ S a -S i-

182 Гл. IV. Уравнения Больцмана, Крука Бхатнагара Гросса и др.

(См. рис. 4.7.) Из последних формул явно видно, что существует

целый континуум скоростей

рассеяния ( |4', §'), совместимых

в силу уравнений сохранения

с начальными скоростями ( |1; | 2).

Этот континуум порождается при варьировании единичного апсидального вектора ос. Однако, если заданы s, V и взаимодействие f (г), то динамическое решение для орбит определяет значение ос, (т. е. либо Ѳ, либо ф), которое в свою очередь фиксирует (§', Ij).

Здесь остаются нерешенными две задачи. Во-первых, угол рассеяния Ѳдолжен быть записан как функционал от s, V и потен-

Р и с. 4.10. Координаты в эксперименте по рассеянию.

циала взаимодействия Т . Во-вторых, мы должны определить, какая часть пучка (невзаимодействующих) частиц рассеивается на данном рассеивателе. Мы найдем, что, сформулировав решение второй задачи, мы естественным образом придем к решению первой.

Рассмотрим однородный пучок падающих частиц с энергией Е и интенсивностью I (число частиц/сек см2). (См. рис. 4.10.) Число частиц, рассеянных в элементе телесного угла dQ около Q, пропор­ ционально интенсивности падающего пучка I и элементу телес­ ного угла dQ. Множитель пропорциональности а — это диффе­ ренциальное сечение рассеяния.

Іа (Q) dQ равно числу частиц, отклонившихся в пределах dQ

около О за секунду.

(4.16)

Но такое же число частиц прошло через элемент кольца s ds dcp, так что

Іа dQ = I d(fs ds.

(4.17)

Азимутальный угол cp, который определяет сечение падающего пучка частиц,— это тот же угол, что и в сферической системе

4.2. Концепция столкновений

183

координат с началом, фиксированным в рассеивателе. Таким

•образом, dQ = d cos Ѳdcp. Подставляя это равенство в (4.17), получим

о ( £ , Ѳ ) = - » ( £ , Ѳ ) | ^ . .

(4.18)

Как находится s = s (Е , Ѳ), будет описано ниже. В то время как о связано с числом частиц, рассеянных в определенном направле­ нии, величина

2 = | а с Й 2 = лг02

(4.19)

4я

полное сечение рассеяния, связана с полным числом частиц, рассеянных из падающего пучка. Радиус взаимодействия равен г0,

так что частицы с прицельным параметром s >> г0 не рассеиваются из пучка. Полное поперечное сечение 2 — это заградительная площадь, которую представляет рассеиватель для падающего пучка частиц. Для однородного пучка с площадью поперечного сечения А ^ 2 величина 2 ІА является долей частиц, рассеян­ ных из пучка (по всем направлениям).

В дальнейшем мы вернемся к более детальному описанию того, как вычисляется дифференциальное сечение рассеяния а (Ѳ, Е), что, в свою очередь, связано с более точной оценкой угла рас­ сеяния Ѳ. Этот угол просто вычисляется при исследовании приве­ денного движения фиктивной частицы массы р, управляемого относительным гамильтонианом. Хотя сама частица фиктивная, ее радиус-вектор г — нет. Действительно, г является относитель­ ным вектором между и т 2. Угол между векторами гдо и гПОСле (прежнее обозначение соответственно г и г ' ) равен углу рассея­ ния Ѳ.

Поскольку относительное движение, описываемое вектором г, лежит в плоскости, состояние системы определяется двумя обоб­ щенными координатами, углом ф и модулем г. Относительный

гамильтониан частицы с массой р и импульсом

р,

движущейся

в поле центральной

силы

с потенциалом

Tif),

 

 

 

Н = ^

+ Т(Г),

 

 

 

(4.20)

запишется через эти

переменные

в виде

 

 

 

 

 

Ч =

4

г +

Ѣ

+

Г

( г ) -(4 '21)

Так как Н цикличен по ф, то рфявляется константой движения, обозначаемой через L. Относительный гамильтониан сам по себе также постоянен, поскольку энергия совокупности частиц, изме­ ренная в системе координат, связанной с центром масс, будет

184 Гл. IV . Уравнения Больцмана, Крука БхатнагараГросса и др.

равна

 

 

 

 

 

 

Е = 2Ргц ^ 2цг2 ■Г ( Г ) .

 

(4.22)

Уравнения Гамильтона дают:

 

 

 

Рг

dr

- У І ( Е - Г ) -

L*

(4.23)

dt

[г2?-2

Р

 

 

d<$>

[ХГ^

d$=* L dt

 

(4.24)

 

dt

 

 

 

Выражая в первом уравнении dt через dj>, получим

 

* = 1

 

L dr/r^

■const.

(4.25)

у 2ц [Е ‘V0 (r)]-(L2/r2)

 

 

Типичный случай рассеивания вследствие отталкивания в системе координат, где начало вектора г фиксировано, показан на рис. 4.7.

Р и с . 4.11. (Сечение рассеяния одинаково для всех трех столкновений.)

Это не есть система координат, связанная с центром масс. В послед­ ней г вращается относительно фиксированного центра масс (рис. 4.5). Однако, как мы видели, угол рассеяния Ѳв этой системе координат тот же самый, что и в системе координат центра масс. В каждом случае Ѳ — это угол между гд0 и гПОслег который не зависит от выбора начала координат.

4.2. Концепция столкновений

185

Диаграмма рассеивания (рис. 4.7) симметрична относительно линии апсиды гт іП, так что г = 0 для rmin. Эта симметрия являет­

ся следствием

сохранения

углового момента

L =

p,Fs = цКѴ

и относительной скорости

V — V', так что s = s '.

Поскольку г

обращается в

нуль на апсиде, из уравнения

(4.23)

следует, что

в этой точке знаменатель в интеграле (4.25) равен нулю. Однако сам интеграл остается ограниченным.

Угол я|), изображенный на рис. 4.11, получается при интегри­ ровании (4.25) от Гщіп до оо. В результате будем иметь

(£/r12) dr

(4.26)

у г р ^ - 'Г Д - ^ / г 2)

Угол рассеяния Ѳ связан с -ф, как показано на рис. 4.7:

 

Ѳ +

2ф = л.

(4.27)

Лучше выражать угол ф через

величину, обратную

радиусу,

и =

г-1. Кроме того, угловой момент L = p,Fs связан с энергией

Е = -g -^F2 соотношением

 

 

 

 

 

Е =

 

L2

(4.28)

 

 

 

2jxs2

 

Если

этим воспользоваться,

то

(4.26) примет вид

 

 

j

 

 

s du

(4.29)

 

 

 

 

 

 

 

 

f

 

 

1 - S 2uz - { T

 

{U)IE) = 0.

 

Так как ф определяет вектор а, то (4.29) представляет собой явное решение для а, выражаемое через Т , E n s . Для потенциала в форме

Т (г) = Kr~N = KuN

(4.30)

(4.29) примет вид

 

s du

(4.31)

У 1 — s2«2—Kun/E

 

Вводя безразмерный обратный радиус

(4.32)

ß = su

186 Гл. IV. Уравнения Больцмана, Крука БхатнагараГросса и др.

и безразмерный прицельный

параметр

 

7

/ Е \ U N

(4.33)

 

 

в (4.31), получим

 

 

W = !

p l -р - ф /!,)» ’

(4.34)

 

Задача 4.4. Для притягивающих потенциалов из (4.34) полу­ чаем:

? - * + ( ! ) " ■

Показать, что это уравнение не имеет положительных решений для N ^ 2 и b < 1.

Ответ. Перепишите уравнение в форме

Р2 = 1 + ( т Г ^ -

Как мы далее увидим, Vo d cos Ѳ является важным множите­ лем, входящим в уравнение Больцмана. Вспоминая (4.17), полу­ чим, что

 

Vo d cos Ѳ =

Kscfe.

(4.35)

Через

безразмерный параметр Ъ это сечение выражается следую­

щим

образом:

 

 

 

Vo d cos Ѳ= ( ) 2/N

db.

(4.36)

Полученный результат явным образом демонстрирует зависимость печения рассеяния от скорости V, которая целиком заключена в мультипликативном коэффициенте FC,- 4>/iV.

Очевидно, что случай N = 4 заслуживает особого внимания. Молекулы, взаимодействие которых описывается таким потен­ циалом, называются максвелловскими молекулами. Для таких молекул вес рассеяния Vo d cos Ѳ не зависит от V. Этот факт значительно упрощает вычисления, относящиеся к линеаризован­ ному уравнению Больцмана. С этими вопросами мы столкнемся в гл. V, где линеаризованное уравнение Больцмана будет подробно обсуждаться.

Сейчас обратимся к вычислению а для кулоновского взаимо­

действия, когда Т = К /r. Интегрируя (4.34) при N

= 1, получим

= ~ 1§2Ф-

(4.37)

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