
книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений
.pdf4.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) |