Автокорреляционные функции в обобщенной модели Больцмана-Энскога Текст научной статьи по специальности «Математика»

Автокорреляционные функции в обобщенной модели Больцмана-Энскога

Н. Г. Иноземцева1,а, И. И. Масленников2,6, Б. И. Садовников2

1 Университет «Дубна». Россия, 141980, Московская область, Дубна, Университетская ул., д. 19.

2 Московский государственный университет имени М. В. Ломоносова, физический факультет, кафедра квантовой статистики и теории поля. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2.

E-mail: а nginozv@mail.ru, ь ilyamaslennikov@mail.ru

Статья поступила 18.09.2012, подписана в печать 28.09.2012.

Цель работы состоит в исследовании асимптотических свойств временных автокорреляционных функций для нелинейной обобщенной модели Больцмана-Энскога, содержащей дальнодействующую компоненту взаимодействия между частицами. На основе анализа нелинейных особенностей кинетического уравнения Больцмана-Энскога непосредственно выявляется роль нелинейных эффектов при приближении к состоянию равновесия. Показано, что автокорреляционные функции обладают степенной асимптотикой а эффекты, связанные с «включением» дальнодействующей компоненты, приводят к изменению коэффициента при . Полученные результаты устанавливают замкнутое выражение для определения коэффициентов в асимптотическом разложении автокорреляционных функций скорости и термодиффузии.

Ключевые слова: модель Больцмана-Энскога, кинетические уравнения, автокорреляционные функции.

УДК: 533.7. PACS: 51.10.+у; 05.20.-у; 05.20.Dd.

Введение

Проблема изучения асимптотического поведения временных автокорреляционных функций рассматривалась в рамках различных подходов к описанию явлений переноса. Один из них, основанный на последовательной схеме кластерных разложений в цепочке ББГКИ [1-3], позволил связать воедино «аномальные» явления в классической статистической механике твердых сфер, а также учесть эффекты высших порядков по плотности в системе. Достоинством другого подхода, основанного на анализе нелинейных особенностей кинетических уравнений типа уравнений Больцмана [4, 5] либо Больцмана-Энскога, является возможность непосредственного выделения роли нелинейных эффектов при приближении статистических систем к состоянию равновесия. Следует отметить, что интерес к одно-частичным кинетическим уравнениям при описании эффектов высшего порядка по плотности был в значительной степени стимулирован работой H.H. Боголюбова [6], впервые установившего и строго доказавшего существование точных микроскопических решений для этого класса уравнений.

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

Как было показано в работе [6], эта модель обладает точными микроскопическими решениями и может быть использована для изучения высших приближений по плотности. Отметим, что модель рассматривалась ранее в работе [7], где были определены поправки к скорости звука с учетом дальнодействия в системе, а также собственные функции линеаризованного оператора Больцмана-Энскога с дальнодействующей компонентой. Здесь мы покажем, как, исходя из результатов работы [7], найти замкнутое выражение для опреде-

ления коэффициентов в асимптотическом разложении автокорреляционных функций скорости и термодиффузии.

1. Автокорреляционные функции в обобщенной модели Больцмана-Энскога с дальнодействующей компонентой

Кинетическое уравнение модели может быть записано в виде

[(г»' - v)

(Т X

= па2

дt дг

(г>'-г>)<Х>0

х \f{t,r,v)f(t,r+ao■,v'*)

п Г9Ф(|г'-г|) , д .

+ ^ V (1)

т J ог от

гдер(^, г') = г, о- — единичный вектор;

V* = V + — V1* = V1 — <г((®' — ®)о-); а —

диаметр области «жесткого» соударения; Ф(|г' — г|) — дальнодействующая компонента потенциала взаимодействия двух частиц в системе.

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

Cs(t) = lim V~l ( 5>("*(0)) } ,

N

n = — = const,

(2)

где V — объем системы взаимодействующих частиц; (...) — усреднение по равновесному большому каноническому ансамблю; — скорость к-й частицы в момент времени 8 =(■/], А),

jv(v) = mvxvy, jx(v) = -(mv2 - 5ß l)vx, ß =

где Т — температура; m — масса частицы.

1

КьГ

Для произвольного взаимодействия в системе выражение (2) может быть преобразовано к виду

Q(t) = п2

düo ф{у o)js(vo)

ArAv 4>{v)js{v)4!{t, г, v), (4)

где 4>{t,r,v) — отклонение функции распределения f{t,r,v) от равновесного значения (¡>{v).

Будем использовать систему нормировки, в которой (¡>{v) = (/Зт/2-тг)3/2 exp(^l/2/3mo2), при этом f(t,r,v) = cj>(v)(l +4'(t,r,v)). Заметим, что следствием преобразования формулы (2) к виду (4) является ограничение на функцию 4>{t,r, v)

Ф(0, г, о) = [пфШ ~lN(r - r0)S(v ^v0) + f (г), (5)

где /(г) — произвольная функция координат; первое слагаемое в (5) описывает 5-образное отклонение от равновесной по скоростям в начальный момент времени. Выражение (5) содержит значительный произвол; величина оо может быть выбрана в интервале (—оо, оо). Единственным ограничением на N(r — ro) является условие нормировки

drAf(r-ro) = 1.

(6)

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

Таким образом, для вычисления автокорреляционных функций согласно (4) необходимо, как и в случае уравнения Больцмана-Энскога, найти решение уравнения (1) с начальным условием (5).

С учетом соотношения = +

имеем для функции уравнение

d^it г v) <7Ф ,

— ' ' +г>——nA(vW(t, г, v) = па1 at or

х T[V(t, r, v)V(t, г, dü' der +

1 дф д

п

т ф(и) dv дг п 1 д(фЧ') д т ф(и) dv dr

dr' Av' г', г/)Ф(|г - г'|) +

dr' di/ ф{у')ЪЦ, г', г/)Ф(|г - г'|),

(7)

где операторы А(г>) и Т определены следующими соотношениями:

Х(г>)Ф0, г, v) = а2 [(ü' - v)cr] г, v*) +

+ r+aer)v'* - 4>(t, г, v) - 4>(t, r-atr)v'] ф(у') dü der,

(8)

T[V(t,r,v)V(t,r,v')] =

= 4>(t, r, v*)4>(t, r+acr, v'*) - 4>(t, r, v)4>(t, r arr, v').

(9)

Отметим, что эффекты конечных размеров области «жесткого» соударения рассматривались ранее в работах [5, 8, 9], где было показано, что смещение пространственных аргументов в (8) и (9) на ±ает приводит к поправкам высшего порядка по плотности системы. Здесь при рассмотрении влияния дальнодействующей

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

(I gradФ|N -1

т. е. размеры области «жесткого» соударения значительно меньше расстояний, на которых существенно изменяется потенциал Ф. В уравнении (7), таким образом, будем использовать операторы Ао(г>), То:

Ло(®)\1>{t,r, v) = а2

[(»' — v)(T] х

[{v'—v)tr\ x

(г>'-г>)<х>о

x г, г>*)+Ф(^, r, v'*)-4f(t, r, v)-4?(t, r, v')\<j>(v') dv' der

(10)

— линеаризованный оператор Больцмана;

To[V(t,r,v)V(t,r,v')] =

= r, v*)V(t, r, v'*) - r, v)V(t, r, v'). (11)

Нелинейное уравнение (7) естественно исследовать методом теории возмущений. Действительно, для асимптотически больших значений t возмущение г, мало (f(t,r,v)^Kj>(v), t-¥ оо).

В области малых t нелинейная часть уравнения (7), соответствующая «жесткому» бинарному соударению, также мала по сравнению с линейными слагаемыми, поскольку

%№(t, г, v)V(t, г, г»')] [(v'-v)tr]<l>(v') dv' der = 0.

Отдельного рассмотрения требует последнее нелинейное слагаемое в (7). При малых t оно содержит производную от a(v - ®о) ■ Однако, как легко видеть, выражение для автокорреляционных функций (4) содержит усреднение по ®о- При усреднении (7) по ®о нелинейный член исчезает; таким образом, для

этих целей теорию возмущений можно использовать и в окрестности точки t=0. Представляя искомое решение уравнения (7) в виде (12):

Ф = ф® + + ...)

(12)

получим с учетом соотношений (10), (11) и предположения Ф® линеаризованное уравнение для Ф«»

dt ' ' дг п 1 дф д m ф dv дг

<1 г' &о' ф(ю')¥0)У, г', г»')Ф(|г ^ г'|) = 0.

(13)

Для следующего члена разложения (12), очевидно, имеет место неоднородное уравнение

• V-

dt дг

п 1 д д

m ф dv дг

л

иА0Ф(1) -

dr' di/ ф{у')¥1)Ц, г', г)')Ф(|г ^ г'|) =

= па

(v'^v)cr] Т [Ф(0)(t, г, v)V(t, г, г»')]ФЫ) df der+

12 ВМУ. Физика. Астрономия. № 1

п 1 д{ф¥^>) д т ф

dv дг

йг' <!©' г', г/)Ф(|г - г'|).

(14)

Отметим, что начальные условия для можно

выбрать в виде

Ф(0)(0,г, = Ф(0, г, V),

Ф(1)(0,г, ©) = 0.

(15)

Однородное уравнение (13) не содержит явной зависимости от аргумента поэтому его решение с начальным условием (15) может быть легко найдено:

¥0)(t, г, v) = e^w"iv)t¥0)(0, г, v),

где

№ф{и)т](г, v) =

где

Q(t, г, v) = па2

— v)(T] х

(г>'-г>)<Х> О

х ?o[e^iv)t4>(0,r,v)e^iv,)t4>(0,r,v')] dvd<т + m ф dvy '

дг

Cf\t) = n2

Cf\t) = n2

Cs{t) = C«\t) + Cf\t),

d®o Ф{роШ®о) x d r (k> |j(D)^D)e"f»wf(0, r, v),

dr d®

3. Асимптотика автокорреляционных функций

Отметим, что вследствие соотношения | dг>o Ф(щ) х х М®о) = 0 выражение (19а) можно записать в форме

Cf\t) = n

Поскольку функция Ы(г - го) нормирована условием (6), а (20) содержит интегрирование по г, можно опустить в Wф{v) все слагаемые, содержащие оператор

д/дг, после чего найдем = пЦ^, е"А°{ >'¡¡(V)) —

стандартное выражение для автокорреляционной функции в линеаризованной модели Больцмана, экспоненциально убывающее с ростом I [5].

Переходя к исследованию выражения (196) описывающего нелинейные эффекты в нашей модели, заметим, что структура С}г, (18) позволяет выделить в СР(0 два слагаемых:

(16)

cf\t) = cf(o + cf(o,

(21)

где

T](r, v) -

dr' dv' ф{^)г]{г', г/)Ф(|г - г'|). (17)

п 1 дф д т ф ду дг,

Поскольку зависимость правой части неоднородного уравнения (14) от I может быть определена с помощью (16), (17), его решение с нулевым начальным условием можно представить в виде

I

Ф(1)(^, г, V) = г, V) dt',

о

cif(0 =

паi

dk

dvф{v)is{v)

e(t-t')nAo(v) fa! x

(2-тг)3 dv' ф(ч)')%

X - »') + nfkN;uW) + ntkNkU{v)

(21a)

где fa, Nk — фурье-образы функций /(r), N(r — ro) соответственно; оператор Hk{v) определен соотношением

nk{v)r){v) =

п 1 дф

dv' r)(v')<&{k),

(22)

где Ф(й)

ненты;

т ф dv

фурье-образ дальнодействующей компо

dr' dv' ф^е-ЪМ'Щ0, г', г/)Ф(|г - г'|). (18)

Cf = *

Й2 т

d®o <М®оШ®о)

dr dv

e(t-t')nAo(v) fa! х

С учетом (16), (18) получим следующее разложение для автокорреляционных функций:

(

(19)

(19а)

' ф(ю) dv д_ дг

dr' dv' -J_^(r-ro)5(c-co)+n/(r)\

\ФШ /

(21b)

Применяя при интегрировании по г в (216) соотношение J7(r)g(r)dr = $f(k)g(—k)dk, представим выраже-

ние для Cf\t) в виде

е\%W'-fiQ(t>,r,v)dt'. (19b)

п 1

1

dk

dvф{v)jii{v)

e(t-t')n\o(v)

dv' фЫ)-гг^г X >(o)

д

drd® - r0)js(v). (20)

x

x {'-щШ28(v^v') + nfkNZ~js(v') + nfZ~Nkjs(v)) x

x (-/АФ(А)). (23)

Асимптотическое поведение С®(О, С®(0 при t ^ оо определяется структурой особенностей функции

(^¡2Чр) = I ер1С^(1)(И в конечной области комплексной о

р-плоскости. Непосредственное вычисление С®(р) дает

СГ(Р) = П[

х Т0

М®)

р - /гЛ0(®)

(24а)

где

(25)

V) = - V) + пЬЩф') +

где {ф{р} — ортонормированный базис собственных функций оператора Ао, соответствующих элементарным гидродинамическим решениям линеаризованного уравнения Больцмана с волновым вектором к [5];

8 = /ЗпФ(0)х1 —.

',(1,2)

не являются взаимно

Отметим, что функции ортогональными.

Разложим величину р^(г>, V1) по собственным функциям 5^(1»), оставляя лишь член с функциями (28):

ли)

(29)

При этом для определения ведущей сингулярности в точке р = 0 достаточно рассмотреть вместо (24а), (24Ь) выражения

С®

(р) =

М®)

т

(к/ ф(ч>')

чр - /гЛ0(®) Ак

(2-тг)3

1

ф(ю)

9®

х

р+/й - тгк (V)(г»')—) —.«Ао (®) —.«Ао ).

(24Ь)

Поскольку среди собственных значений оператора Ло(®) есть ненулевое, очевидно, что для функций (24а), (24Ь) точка р =0 является крайней правой особой точкой, и таким образом поведение величины (24а), (24Ь) в окрестности р =0 определяет искомые ведущие члены в асимптотике автокорреляционных функций при t-¥oo.

Все собственные функции оператора Ло(г»), соответствующие нулевому собственному значению, ортогональны /¿(г>)> скалярное произведение определено посредством (25); вследствие чего структуры

в (24а), (24Ь) не дают сингулярной особенности в точке р = 0, и достаточно исследовать выражение типа

is(v) -/гЛ0(®)'

Аг

(2-тг)3

р + 5А(®) -5_а(®')/

где введено обозначение

5а(®) = - 7га(®)) - /гЛ0(®).

Собственные функции оператора 5а(®), соответствующие собственным значениям, обращающимся в нуль при —0. были найдены в работе [7]. В нулевом приближении по \к\ имеем

-1/2

фм'2)(®) = (2а (а т V»2 - ¿2

ф(/) _ шО) к '

/ = 3,4,5,

¡¡(у) -/гЛ0(®)'

-0-А

(30а)

С<2)

ь2

/¿(р) -/гЛ0(®)'

с:1£

Щ3

■0-А

Ф(А)

(ЗОЬ)

где 4°

собственные значения 5а(®), соответствующие собственным функциям Фоа(®)- Следует подчеркнуть, что р^ не содержат величины поскольку одночастичные токи /¿(®) ортогональны всем используемым нами функциям Фоа(®)-

Величины г^ для малых значений \к\ были вычислены в работе [7]:

,(1.2)

а

2 - 82\к\ + \к

2 Фт + 2Д?)

(31)

,(3)

~ \к\2От;

,(4,5)

Д„

где Ит, Д? — коэффициенты сдвиговой вязкости и термодиффузии, определяемые из уравнения Больцмана.

Заметим, что структура особенности при р=0 определяется интегрированием в (30а), (ЗОЬ) в области малых Щ, причем ведущий член в (30а), (ЗОЬ), определяющий асимптотику корреляционных функций, соответствует тем слагаемым в суммах по (/,/), для которых ^к+^к ~ Щ2 ■ Поскольку числитель подынтегрального выражения в (ЗОЬ) содержит й, в окрестности точки

(28)

р = 0 С^(р)>>С^(р) и поэтому С^(0 » при

t—¥oo. Следовательно, для изучения асимптотического поведения С®(0 достаточно рассмотреть выражение

13 ВМУ. Физика. Астрономия. М 1

cf\t)~J h{v)

rtAo(f)

dv' <p(v')

ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2013. № 1 Г dk

(2-тг)3'

х ]. (32)

где Y1 * означает суммирование по значениям индексов (/,/), определяемых условием + ~

const ~k2. Полагая в (32) гЦ} + = k2{\f + А^), найдем

Ф) = -r3/2[js(v),

г dQk 4-тг

d®'

(Af + A^

47Г3/2

. (33)

Таким образом, в нелинейной модели Больцма-на-Энскога с дальнодействующей компонентой автокорреляционные функции также обладают степенной

асимптотикой . Постоянные , X_k могут быть вычислены в рамках теории Энскога; эффекты, связанные с «включением» дальнодействующей компоненты Ф(|г — г'|), приводят к изменению величины коэффициента при ^3/2 согласно (29), (33).

Список литературы

1. Боголюбов H.H. Проблемы динамической теории в статистической физике. М.; J1., 1946.

2. Ernst М.Н., Dorfman J.R. 11 Physica. 1972. 61. P. 157.

3. Dorfman J.R., Cohen E.G.D. 11 Phys. Rev. 1972. A6. P. 776.

4. Ubbink J.Т., Hange E.H. 11 Physica. 1973. 70. P. 297.

5. Иноземцева HS., Садовников Б.И. // ТМФ. 1977. 31. С. 260.

6. Боголюбов H.H. // ТМФ. 1975. 24. С. 242.

7. Бочков С.Н., Иноземцева Н.Г, Садовников Б.И. 11 ОИЯИ. р 17-81-10. Дубна, 1981.

8. Dorfman J.R., Cohen E.G.D. 11 Phys. Rev. 1975. A12. P. 292.

9. Sadovnikov B.I., Inozemtseva N.G. 11 Physica. 1978. 94A. P. 615.

The autocorrelation functions in the generalized Boltzmann-Enskog model N.G. Inozemtseva1,0 , I.I. Maslennikov2'ft, B.I. Sadovnikov2

1 International University «Dubna», Universitetskaya str. 19, 141980 Dubna, Moscow Region, Russia.

2 Department of Quantum Statistics and Field Theory, Faculty of Physics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia.

E-mail: a nginozv@mail.ru, b ilyamaslennikov@mail.ru.

The asymptotic behavior of autocorrelation functions for the nonlinear generalized Boltzmann-Enskog model with a long-range interaction between particles is examined. On the basic of this examination of nonlinear peculiarities of Boltzmann-Enskog kinetic equation the role of nonlinear effects near the equilibrium state is established. It is shown that the asymptotic behavior of autocorrelation functions is characterized by t^2 and effects of the long-range interaction lead to the change of the coefficient in the t^2. The obtained results give a closed expression for determing the asymptotic expansion of the autocorrelation functions of the velocity and thermodiffusion.

Keywords: Boltzmann-Enskog model, kinetic equations, autocorrelation functions. PACS: 51.10,+y; 05.20.-y; 05.20.Dd. Received 18 September 2012.

English version: Moscow University Physics Bulletin 1(2013).

Сведения об авторах

1. Иноземцева Наталья Германова — докт. физ.-мат. наук, профессор; e-mail: nginozv@rnail.ru.

2. Масленников Илья Игоревич — аспирант; e-mail: ilyamaslennikov@mail.ru.

3. Садовников Борис Иосифович — докт. физ.-мат. наук, профессор; тел.: (495) 932-80-10; e-mail: sadovnikov@phys.msu.ru.