книги из ГПНТБ / Либов Р. Введение в теорию кинетических уравнений
.pdf5.2.Свойства линейного болъцмановского оператора столкновений 287
комбинация |
сумматорных инвариантов (1, |
£2): |
|
|||
|
|
|
Х фф > 0, |
|
(5.79) |
|
Хфф = 0, |
только если ф = а + |
ß-£ + |
(5.80) |
|||
Рассмотрим |
разность |
|
|
|
|
|
|
0 = |
X фя|) |
j |
dpidp (ф.Йл|ф — фКфі). |
(5.81) |
|
Меняя местами | |
и |
во втором интеграле, |
получим |
|
||
|
0 = |
Сdpdpj [фК (І, |
| t) % — ф Д (It, 1) ф], |
|
||
|
|
, |
. „ |
„ |
|
(5-82) |
|
0 = ] dyidii^K d, 1 0 - K & , i)] ф(i) ф (go. |
|
В последнем выражении подразумевается, что ядро К (х, у) дей ствует только на функции от у. Если мы определим ядро
Т (1, Іі) = К (1, Іі)—К (|і, і), |
(5.83) |
то (5.82) можно записать в виде
І!)С(1, 1 0 = о,
G(i, іо = ^ (і)я И іо
для всех допустимых ф и ф . Функции ф и ф называются допусти мыми, если
II ф II2 = |
!фф < |
оо и II ф II2 = І щ < оо, |
|
где 7 == ^ dp,, а II ф II |
называется нормой функции ф. (|| |
ф || также |
|
называют Т/2-нормой |
функции |
(ф). Множество всех |
функций |
с конечной Т,2-нормой образует ^-пространство.) |
|
Предположим, что ядро Т |
не равно тождественно нулю. Это |
||||
означает, что существует такая функция G, что TG отлична от нуля |
|||||
в некоторой области пространства |
1, l t. Тогда существует под |
||||
область D этой области, в которой TG не меняет знак. |
Определим |
||||
G = G в D и |
G — 0 вне D. Для |
такой |
функции |
G интеграл |
|
^ TG =7^=0. Это |
противоречит |
равенству |
(5.84). Следовательно, |
единственный способ удовлетворить условию j Тфф = 0 для всех
Ф и ф — это тождественное обращение в нуль ядра Т . Итак, К —симметричное ядро.
290 Гл. V. Решение и дальнейшие свойства уравнения Больцмана
в) Дальнейшие свойства % и ѵп
До сих пор мы лишь отметили тот факт, что последователь ность {ѵ„} содержит пятикратное вырожденное собственное зна чение (vj = 0), соответствующее пяти сумматорным инвариантам.
Остальные собственные значения оператора X все лежат на поло жительной части действительной оси.
Недавно Трэд (1963) доказал четыре других свойства спектра
(т. е. последовательности собственных значений) оператора X . Для уяснения этих свойств необходимо сначала ввести концепцию жесткого потенциала, который определяется через частоту столк новений (ср. с (5.70)):
(5.96)
Потенциал является жестким, если ѵ (£) изменяется от минималь ного значения до бесконечности, и мягким, если ѵ (£) изменяется от максимума до нуля. Отметим, что потенциал жесткого ядра, рассмотренный в задаче 5.4, включается в класс жестких потенциа лов, так как для него ѵ (|) возрастает как £ для больших Для взаимодействий по закону обратной степенной зависимости (см.
задачу 4.39) G ~ r~N и при N > |
5 функция ѵ (5) является моно |
тонно возрастающей функцией |
изменяющейся с ростом \ от |
своего минимального значения ѵ (0) до бесконечности. Следова тельно, взаимодействия такого рода относятся к классу жестких потенциалов. Для N < 5 частота ѵ (1) монотонно убывает от своей максимальной величины ѵ (0) до нуля, когда | неограниченно растет (рис. 5.3). Для особого случая максвелловских молекул N = 5 и V = ѵ0 не зависит от £.
Для широкого класса пороговых потенциалов Трэд (1963)
установил следующие свойства оператора X .
1. Для жестких потенциалов частота vj = 0 изолирована от
остального спектра частот оператора X . Кроме ѵ4 = 0, спектр X содержит континуум значений, принимаемых частотой столкнове ний V (I), т. е. все значения от минимума ѵ (0) до бесконечности. Для мягких потенциалов континуум собственных значений про стирается от максимума ѵ (0) до нуля.
2. Запишем X в виде
£ = ѵ ( £ ) - 5 Г ,
(5.97)
0Г= j d ^ K d , 1,)
292 Гл. V. Решение и дальнейшие свойства уравнения Больцмана
то легко установить следующее. Если (г|/, Я') —соответственно
собственная функция и собственное значение оператора (/ |
— Ж'), |
|
т. е. если |
• |
|
(/ —Й Г ) ^ '= ^ Ѵ > |
(5.100) |
|
то |
|
|
где |
= |
(5.101) |
|
|
|
|
|
(5.102) |
И > |
vX’. |
(5.103) |
X = |
Это преобразование лучшевсего годится для случая усеченного максвелловского закона взаимодействия (см. разд. 4 ниже), при котором V не зависит от
3.Оператор Ж является вполне непрерывным. Ограниченный
линейный оператор Ж называется вполне непрерывным, если он обладает следующим свойством. Для любой ограниченной после
довательности {ф„} функций |
пространства |
Ь 2 (это означает, что |
|| ф„ К^ R для некоторого |
R и любого |
п) последовательность |
{Ж ^п} содержит по крайней мере одну сходящуюся подпоследова тельность. Из этого факта и спектральной теоремы х) следует,
что Ж имеет дискретный точечный спектр с единственной предель ной точкой в начале.
Физический смысл непрерывной части спектра £ заключается в следующем. Предварительно отметим, что собственные значения VI входят в решение задачи с начальными значениями в виде чле нов ехр (—vit). Для жестких потенциалов, где ѵ отделена от нуля
(0 < V< |
ѵ), все вклады экспоненциальной формы |
убывают по |
крайней мере так же быстро, как ехр ( — ѵ(0) t). С другой стороны, |
||
для мягких потенциалов ѵ может быть как угодно близка к нулю |
||
и вполне вероятно, что некоторые свойства функции |
не будут |
|
убывать со временем по экспоненциальному закону. |
|
|
4. |
Для усеченного максвелловского взаимодействия G ~ г~ь |
для г ^ r0; G — 0 для г > г0; спектр оператора £ является чисто дискретным с единственной точкой накопления ѵ = const. Кроме того, для максвелловских молекул собственные функцйи операто
ра £ будут полиномами Эрмита Я[п) с соответствующими собствен ными значениями:
ѵ(іП) = (Н\п) I £ I Н\п)). |
(5.104) |
0 Краткое изложение этих вопросов можно |
найти у Лорча (1962). |
Более полное представление см. в книге Стоуна |
(1932). |
294Гл. V. Решение и дальнейшие свойства уравнения Больцмана
5.3.Моментный метод Трэда
В1949 году Трэд в своей докторской диссертации (1949) раз работал более общий метод решения уравнения Больцмана. Он известен сейчас как моментный метод Грэда. Причина такого названия вскоре станет ясна. Моментный метод отличается от метода Чепмена — Энскога, поскольку построенные решения не
имеют функциональной формы |
(х, | |
| р, u, Т). |
|
||
Исследование |
начинается с |
определения |
безразмерной ско |
||
рости ѵ: |
|
|
|
|
|
с = |
І —u(x, t) = v V R T |
(u= / |
2 s) |
(5.110) |
и безразмерного распределения g для функции ^ , нормированной
на полную массу системы М = J Цр d \ dx:
SF = |
----= ------------------- |
(RT)3/2 |
^72-со (у). |
(5.111) |
|
(RT)3/2 |
|
' |
Отсюда сразу следуют три свойства функции g:
j |
gdx = 1, |
|
j g v d v = 0, |
(5.112) |
|
j |
gv2dx = 3. |
|
Функция распределения g, |
разложенная по полиномам |
Эрмита |
с весовой функцией со, определенной равенством (5.105), имеет вид
оо
g(*, ѵ,і) = 4 ) Ц |
І «in) (*, t)н г (V) |
(5.113) |
п—0 |
і |
|
или, что то ясе, можно выразить jF через локальное распределе ние Максвелла
Функция # -n) является тензором порядка п и полиномом сте пени п. Нижним индексом і обозначена совокупность из п индек сов. Некоторые примеры таких функций даны в (5.108).
В развернутом виде (5.113) будет выглядеть так:
g = «в [ a ^ H ^ + a ^ H ^ + ^ r alfHir + ~ a b kHbk+ . . . ] . (5.115)
Разложение справедливо (в смысле средней квадратичной схо димости), если функция co-1/2g- квадратично интегрируема, т. е.
5.3. |
Моментный метод Грэда |
295 |
|
если существует интеграл |
|
|
|
И |
|
(5.116) |
|
|
|
|
|
Обращая ряды (5.114), |
получим |
|
|
аГ = \ gH ? ' d y = - |
f |
(5.117) |
|
J |
P |
J |
|
Так как функции Н\п) — полиномы по скорости |, то коэффи
циенты а<и) являются моментами функции распределения fp. Каждый из этих моментов будет некоторой гидродинамической переменной. Коэффициенты a(jn> разложения функции распределе ния по полиномам Эрмита являются гидродинамическими пере менными. Например,
а<0) = 1,
а? = О,
і —■ |
I)' -рЬц |
(5.118) |
Tr а<2>= О, |
||
|
р |
|
Hjj ■ |
2Qi |
|
) Укт |
|
{повторяющиеся индексы означают свертку).
Чтобы построить решение JP, мы должны получить выражения для коэффициентов а-'1’. С этой целью мы сначала подставим раз
ложение (5.114) для jF |
в уравнение Больцмана и над полученным |
||||
выражением |
произведем операцию j |
І7|п> dv. Полагая |
|
||
|
|
|
C = y~RT, |
|
|
получим |
|
|
|
|
|
Да(п) |
I ди |
\ |
|
да\п+1> |
|
Dt |
V дхі |
) |
|
дхі |
|
|
|
+ C a ^ JL- ln (рСп+1) + |
С I а171"1»+ |
|
|
|
+ [ } § + « |
|
+Ь ( Ц ) |
+ |
|
|
+ ( - ^ ^ б + 2 н ) а<"-2>+ ^ ( 2 С 2)6а(П-3’=/<">. (5.119) |
||||
Для |
простоты записи |
нижние индексы у аіП) и / <П) |
опущены; |
||
|
|
|
|
V |
|
D I D t |
—индивидуальная производная, а символом 2 обозначена |