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

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

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

5.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 для всех

Ф и ф — это тождественное обращение в нуль ядра Т . Итак, К —симметричное ядро.

288 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

б) Положительность собственных значений и сумматорная инвариантность собственных векторов х)

Собственные значения ѵ„и собственные функции фп оператора £ определяются обычным образом:

І :і1’п = ѵп'Фп

(5.85)

(при условии, что £ имеет дискретный спектр). Матричные эле­ менты, соответствующие этим собственным векторам, являются

диагональными элементами, что обусловлено симметрией £ (см.

разд. 2.3):

<£пт— ^ dlltyn£tym= Ѵт ^ ^рфпфт = ' ’тбпт*

(5.86)

Диагональные элементы удовлетворяют неравенству (5.79) и ра­ венству (5.80), так что

0^ С £ пп = ѵ„ j dptyn .= v„ II фп [I2.

(5.87)

Следовательно,

 

vn > 0.

(5.88)

Равенство нулю будет только в том случае, если фп — любой из пяти скалярных сумматорных инвариантов. Если начать нумеро­

вать собственные значения и собственные векторы оператора £ с этих нулевых величин, то для произвольной силы взаимодействия (!)

первыми пятью собственными функциями всегда будут

Фі = а

'h = PS*

f s = %

0

Фі = y£z' f s = ч г-

где а, ß, б, у, X —произвольные постоянные. Кроме этого пяти­ кратного вырожденного нуля, все остальные ѵп положительны (опять для произвольной силы взаимодействия), т. е.

ѵ„ > 0, п > 5.

(5.89)

Эти свойства ѵп и фп для п ^ 5 являются следствием трех законов сохранения, а сами законы, конечно, не зависят от частного вида силового поля. Положительность остальных собственных значений вытекает из оЖ’-теоремы Больцмана.

Чтобы убедиться в этом, рассмотрим задачу с начальными данными для пространственно-однородного (д/дх = 0) случая,

1) Метод собственных функций в анализе линеаризованного оператора столкновений первоначально был введен Вэнь-Ченем и Уленбеком (1952).

5.2.Свойства линейного болъцмановского оператора столкновений 289

когда уравнение (5.72) принимает вид

 

 

4 г = - ^ -

(5-9°)

Будем

искать решение этого уравнения,

обращающееся в ф0 (|)

при t

= 0:

 

 

Ф (§, 0) ^ фо (I).

(5.91)

Если ф = ап (t) подставить в предыдущее уравнение, то мы получим ап = ехр ( —vnt). Любая суперпозиция этих решений также является решением уравнения (5.90). Чтобы построить решение, обращающееся в ф0 при t = 0, полагаем

Ф(1. t) = 2 апехр ( — vnt) (1),

(5.92)

так что

 

<Ы!) = 2 апФп(1).

(5.93)

Вследствие ортогональности последовательности {фп} из равен­ ства (5.93) имеем

ап = j Фо (1) фд (І) dp,

(5.94)

что вместе с (5.92) дает искомое решение.

Параметры ѵй1 являются временами релаксации для собствен­ ных состояний ф„. Будем считать, что начальное значение функции распределения ф (g, if) задается (невырожденной) собственной функцией ф„. Тогда соответствующим решением уравнения (5.90) будет

Ф (I, t) = фп (1)

ехр {—vnt).

(5.95)

Система остается в состоянии ф„

и затухает

экспоненциально

во времени. Время релаксации этого собственного состояния — время, за которое ф0 уменьшится в е раз по сравнению со своей первоначальной величиной, — равно ѵй1 • Так как все собственные значения больше либо равны нулю, то ни одно из собственных состояний системы не растет со временем.

Собственные частоты {ѵ„} интересны с формальной точки зре­ ния. Эти параметры связаны с микроскопическим описанием систе­ мы. В частности, они описывают эволюцию во времени функции распределения ер (до членов О (ф)). Если какие-либо собственные частоты отрицательны, то в системе может развиться микронеустойчивостъ. Эта неустойчивость связана с функцией распределения. Простым примером такого явления служит непрерывный рост числа частиц в некотором интервале скоростей. Понятно, что такой процесс может не вызвать никакого заметного изменения в любой из макроскопических переменных.

1 9 - 0 1 2 4 3

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,)

5.2. Свойства линейного больцмановского оператора столкновений 291

Г

5«?

<2

Жесткий потенциал

ѵ(0)

Мягкий потенциал

Р и с . 5.3. Частота столкновений, потенциал взаимодействия и спектр опера­ тора столкновений для жесткого и мягкого потенциалов.

где V задается равенством (5.96), а К — равенством (5.70). Опера­

тор ограничен, т. е. для любогоij)из Ь2 (см. стр. 287) существует постоянная М, такая, что

||е%ЧИІ<ЛП|гИ|.

(5.98)

Кроме того, если мы положим

 

Г ^Kt (1, St)

(5.99)

J Уѵ Ш V(gi)

19*

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

5.2.Свойства линейного больцмановского оператора столкновений 293

Двойной индекс у полинома Эрмита Я{п> указывает на то, что это полиномиальный тензор ранга п в трехмерном пространстве. Символом і обозначена последовательность из п индексов (ц, і2, . . ., іп)- Полиномы порождаются весовой функцией

 

 

со (х) = -^|зТ 2" ехР (

г г ) •

 

 

 

(5.105)

ПОЛИНОМ Я (П) =

 

 

, in

задается равенством

 

 

 

(п)

_/

л\n

1 /

д

д

д

)

:

(—1)"

Vnco

Я И, *2, . . . , т

— \

)

ш \

дх^

дх І2

дх-^

)

 

или через порождающую функцию

 

 

 

 

(5.106)

 

 

 

 

 

 

 

 

 

СО

 

 

 

 

 

 

ехр [ - у ( а 2- ж 2) ] =

2

н и!...,іп(*)УнУіг •••

г/in,

(5.107)

 

 

 

 

п=О

 

 

 

 

 

 

 

 

 

 

 

 

 

® І~У і

Несколько первых полиномов имеют вид

Я(0) = 1,

Н^ = хь

(5.108)

//(?> —.X■Т ■_iS• •

Яць = е д а * —

+ жД-й -f xk8u).

Нормировка Я (п) задается равенством

f Я (п)

■Я (и)

rix —г? 1 Аіи •••’1"

(5.109)

J Л Н,

. , . , 1 п л Л, •

. • J n W W X

— Ai!

. . . , 3 n>

где

дИ, • • .. __ ^

 

 

 

 

 

если {гД является перестановкой {/Д, и

 

 

 

ді.і,

• • ■>Щ_

о

 

 

в противном случае. Таким образом,

Я™

образуют ортогональ­

ную с весом dp, систему функций (см. (5.74)).

В следующем разделе мы коснемся в общих чертах метода иссле­ дования нелинейного уравнения Больцмана, который в значи­ тельной мере опирается на только что рассмотренные свойства

полиномов Эрмита.

 

Задача 5.5.

Установить, является ли вполне непрерывным

на интервале

(0,1)

при нулевых граничных условиях оператор

d2/dx2.

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 обозначена

296 Гл. V. Решение и дальнейшие свойства уравнения Больцмана

симметричная сумма градиентов. Например,

V

а<2)

да$

 

 

2

дхк 1 dxj ' дхі

(5.120)

 

 

Символом б обозначен дельта-тензор Кронекера1). Столкновительный член / (П) имеет вид

/<"> =

g (ѵ) g (Ѵі) Г П) (v, vj) dv dv1?

(5.121)

где

/ (ТІ>_ j 7айЙ[Я<П)],

(5.122)

[Я(П)] = Я<п>(y') + Я <П) (ѵО -Я<"> (v) — Я (П) (vt).

Уравнение (5.119) тензорное порядка п. Полезно самим убедить­ ся, что каждый член этого уравнения является тензором /г-го- порядка. Примером может служить член

 

 

 

Д(п)=: J _ C26a(n-l,

 

 

 

дхі

если его записать более

подробно:

 

 

 

_д_

 

 

 

 

дхі:С % ка '^ г Ѵ ...= ^ г о ...

 

Задача 5.6. Доказать следующие тождества:

( 1)

д

Я<П>= 6;Я<П-1>,

 

дѵі

 

(2)

і7іЯ

‘п, = Я{п+1> +

біЯ ‘п' 1>

 

(3)а гбгЯ (П>= а Я (П),

(4)

Vf Sд- Я<п>= пН(П) 2бЯ(П-2).

1 дѵі

Выражение с бгпредставляет сумму по всем членам, в кото­ рых индекс і привязан к 6. Например, (^-компонента 6гЯ а) равна

б4Я (1>= бгьуг + бггуй.

Метод решения состоит в аппроксимации рядов (5.113) конеч­ ным числом членов. Так как старший член разложения пред­ ставляет собой локальное распределение Максвелла (см. (5.115)), то, чем ближе состояние газа к равновесному, тем меньше членов; ряда требуется для описания состояния газа.)*

*) В предыдущих обозначениях 6 = 6; компоненты 6 (или 6) равны 6^..

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