CC BY
74
ФИЗИКА
J
DOI: 10.18698/1812-3368-2016-1-27-35 УДК 536.75
ОПИСАНИЕ ФЛУКТУАЦИЙ СКОРОСТИ БРОУНОВСКОЙ ЧАСТИЦЫ ПРИ ВОЗДЕЙСТВИИ ПУАССОНОВСКОГО СЛУЧАЙНОГО ПРОЦЕССА
А.Н. Морозов
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]
Проведено описание движения броуновской частицы при воздействии на нее пуассоновского процесса. Найдено стационарное решение для характеристической функции флуктуаций скорости броуновской частицы для случая, если воздействие на нее описывается пуассоновским процессом с нормально распределенными скачками. Вычислены первые четыре момента функции распределения флуктуаций скорости броуновской частицы, а также эксцесс и мера Кульбака для этого распределения. Полученные результаты применены для нахождения характеристической функции флуктуаций напряжения на электролитической ячейке. Установлено, что мера Кульбака флуктуаций напряжения на электролитической ячейке обратно пропорциональна интенсивности пуассоновского процесса и числа ионов в малом объеме электролита. Решено уравнение для характеристической функции распределения флуктуаций скорости броуновской частицы, учитывающее флуктуации коэффициента вязкого трения.
Ключевые слова: броуновское движение, флуктуации скорости, пуассоновский процесс, характеристическая функция, вязкое трение, мера Кульбака, электролитическая ячейка.
FLUCTUATIONS OF THE BROWNIAN PARTICLE VELOCITY INFLUENCED BY A RANDOM POISSON PROCESS
А.N. Мorozov
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]
Brownian motion of a particle influenced by a Poisson process is described. The stationary solution for the characteristic velocity fluctuation function of a Brownian particle is found when influence on it is described by the Poisson process with normal jump distribution. The first four moments of the Brownian particles velocity distribution function as well as kurtosis and Kullback's measure for the distribution have been calculated. The results obtained have been used to find the characteristic function of voltage fluctuations in an electrolysis cell. Kullback's measure for voltage fluctuations in an electrolysis cell is inversely related to the Poisson process intensity and the number of ions in a small volume of electrolyte. The equation for characteristic function of velocity fluctuation distribution of a Brownian particle taking into account fluctuations of the viscous friction coefficient is solved.
Keywords: Brownian motion, velocity fluctuation, Poisson process, characteristic function, viscous friction, Kullback's measure, electrolysis cell.
Описание диффузии и броуновского движения как пуассоновских случайных процессов предложено в работах [1, 2], там же получены уравнения, описывающие поведение соответствующих характеристических функций. Задача нахождения стационарных распределений для стохастической системы взаимодействующих частиц рассмотрена в работе [3], а функции распределения флуктуаций неравномерности вращения Земли получены в работе [4]. Цель настоящей работы — нахождение решений уравнения для характеристической функции флуктуаций скорости броуновской частицы и определение первых четырех моментов распределения этих флуктуаций.
Рассмотрим одномерное броуновское движение частицы, если воздействие частиц среды на броуновскую частицу описывается общим пуассоновским процессом Ш (Ь) со скачками, распределенными по нормальному закону. В этом случае пуассоновский процесс Ш (Ь) задается характеристической функцией вида [2]
1
9« ,t) = exp . Yo kT
exp ( - 2D« ) - 1) Vt t
где Б с = 2---дисперсия пуассоновского процесса Ш (Ь), харак-
ут т
теризующая воздействие частицы среды на броуновскую частицу при единичном соударении; 70 — коэффициент вязкого трения; к — постоянная Больцмана; Т — абсолютная температура среды; ут — интенсивность пуассоновского процесса; т — масса броуновской частицы.
Представим уравнение, описывающее одномерное броуновское движение, в виде дифференциального уравнения Ито [5]
¿у = -7оувь + ¿ш (ь). (1)
Здесь У — скорость броуновской частицы. Решение уравнения (1) име-
г
ет вид У (*) = У С (Ь, т) ¿Шс (т), где
С (Ь,т) = ехр [—7о (Ь - т)] . (2)
Одномерную характеристическую функцию флуктуаций скорости движения У броуновской частицы можно представить как [2, 6]
г
9 (Л; t) = exp < vT
exp ( -1D« G2 (t,T) Л2) - 1
dr\ . (3)
Из формулы (3) можно определить первые четыре момента функции
распределения флуктуаций скорости броуновской частицы
дд (А)
Di =
Do =
D3 =
гдА д2д (А)
= 0;
А=0
(гдА)2 д3д (А)
(гдА)3
А=0
/kT
G2 (t, r) dr = —;
m
= 0;
(4)
А=0
D4 =
д4д (А)
(гдА)4
= 3v?D|
А=0
G2 (t, r) dr I + / G4 (t, r) dr
=з(-Y
m
1 + Y0
(5)
где й (¿, т) — функция, определяемая по (2).
Формулы для моментов (4) и (5) позволяют рассчитать эксцесс функции распределения
D4 3D2
к4 =
D2 D2
= 3 Y0,
(6)
выражение для которого совпадает с выражением, полученным в работе [7].
Для экспериментального определения интенсивности пуассонов-ского процесса ут может быть использована электролитическая ячейка, описанная в работах [8-10]. Как показано в работе [10], экспериментально наблюдаются долговременные изменения меры Кульбака Н флуктуаций напряжения на электролитической ячейке, которая, в свою очередь [7], связана с эксцессом простым соотношением
н = К4. 16
(7)
При расчете эксцесса функции распределения флуктуаций напряжения и на электролитической ячейке необходимо учитывать, что эти флуктуации возникают вследствие броуновского движения достаточно большого числа ионов в малом объеме электролита, находящегося в отверстиях в тонкой пленке, которая разделяет сосуды электролитической ячейки [8, 10]. Кроме того, необходимо учесть низкочастотную фильтрацию флуктуаций напряжения, связанную с наличием собственной электрической емкости электролитической ячейки. Для описания флуктуаций напряжения вместо уравнения (1) необходимо
t
2
t
использовать систему уравнении
вУЕ = -ЪУъ<И + dWщ (*); (8)
ви = -вивь + влу^г. (9)
N
Здесь У^ = У — сумма скоростей ионов в малом объеме электроли-
г=1
N
та; N — число ионов; = ^^ W£i — сумма независимых пуассонов-
г=1
ских процессов воздействующих на каждый ион в электролите; в = 1/(0В) — верхняя частота флуктуаций напряжения, снимаемого с электролитической ячейки; В и О — сопротивление и емкость электролитической ячейки; А = h/(pN) — коэффициент (к — толщина пленки; ц — подвижность всех ионов в электролите, которая для простоты полагается одинаковой).
Характеристическая функция пуассоновского процесса (£) имеет вид
' 1 ~ .2
дщ (т ^) = ехр ^ех^- ^ Ю- ^ ^тt
Решение системы уравнений (8) и (9)
г
и Ои (Ь,т) dWSÍ (т),
вл
где Ои т) = -- {ехр [-в (t - т)] - ехр [-70 ^ - т)]}. Тогда од-
То - в
номерная характеристическая функция флуктуаций напряжения и на электролитической ячейке принимает вид
ди (Xu; t) = exp < vTN
exp( -1 De GU (t,T) XU) - 1
dr\ . (10)
В первом приближении из формулы (10) при условии, что в ^ 70, можно определить второй (Юи2) и четвертый (Юи4) моменты функции распределения флуктуаций напряжения на электролитической ячейке:
д2ди (Xu) Du2 = ч2
(гдХи)
= vTNDe f G2u (t,T) dT = в kTk2
Xu=0 Yo mv2N'
t
Du4 =
d4gu (Au)
(idAu)
Л=0
= 3v2n2d|
GU
(t, т) dr
1
+ N
G4
U
(t, r) dr
= 3
в kTh2 \ 2
1+
в
vT N
ч7о у
В соответствии с формулой (6) получим выражение для эксцесса функции распределения флуктуаций напряжения на электролитической ячейке
KU 4 — 3
в
vT N
(11)
Согласно выражению (11), эксцесс для флуктуаций напряжения на электролитической ячейке, а следовательно, и мера Кульбака (7), пропорциональна верхней частоте в фильтрации сигнала электролитической ячейкой и обратно пропорциональна числу ионов N в малом объеме электролита.
Проведенное описание основывалось на предположении, что в правой части уравнения (1) второе слагаемое, характеризующее случайное воздействие частиц среды на броуновскую частицу, представляет собой производную общего пуассоновского процесса Wg (¿). Однако в работе [7] рассмотрен более общий случай, когда уравнение (1) имеет вид ¿V = —VdW7 (¿) + dWg (¿). Здесь W7 (¿) — процесс, описывающий флуктуации коэффициента вязкого трения и задаваемый характеристической функцией
д7 (^7, ¿) = ехр [(ехр (Ш7) — 1) ут¿], (12)
где Д7 = 70г0 — дисперсия пуассоновского процесса W7 (¿), т0=1/^т — постоянная времени, характеризующая среднее время между очередными соударениями частиц среды с броуновской частицей.
Для указанного случая в работе [7] в предположении, что процесс (¿) является винеровским процессом, а в разложении экспоненты ехр (Ш7) из формулы (12) сохранены первые три слагаемых, было получено уравнение
A^ - iV dg<A - Ag (A)—0.
dA2 270 dA
Yom
В общем случае решение уравнения (13) имеет вид [11]
g (A) — Aa+1/2Za+1/2 (¿VbA),
(13)
(14)
4
2
t
где а = —; Ь = 2--; ^а+1/2 (гл/Ьл) — цилиндрическая функ-
7о 7о т V /
ция [12].
Характеристическая функция (14) должна удовлетворять условию [5]
д (А)|А=о = 1. (15)
Представим решение (14) в виде [11]
д (Л) = Аа+1/2 (01^+1/2 (гЛЛ) + 02Nа+l/2 (гЛЛ)) , (16)
где ^а+1/2 (г^Ьл), N„+1/2 (гл/Ьл) — функции Бесселя первого и второго рода; 01, 02 — произвольные константы. Применение условия (15) для выражения (16) дает следующие значения констант:
01 = 0; 02 = -Г (а I 1/2) (V
где Г (а +1/2) — гамма-функция. Тогда окончательно получаем решение (16):
а+1/2
vT kT
д (X) =
п
гУЬХ
N
а+1/2
iVbx^j.
Г (а +1/2)
Первые четыре момента функции распределения флуктуаций скорости броуновской частицы имеют вид
дд (Л)
D1 =
гдХ
= 0;
X=0
D2 =
д 2д (X)
(гдХ)2
X=0
2vT kT
2vT — y0 m '
(17)
D, =
д 3д (X)
d4 =
д 4д (X)
(гдХ)
3 = 0; (гдХ) X=0 12
k2T2
(18)
А=0 (2рт - З70) (2ит - 70) т2
В этом случае второй (Ю2 (17)) и четвертый (Ю4 (18)) моменты совпадают со вторым (к2) и четвертым (к4) кумулянтами [13]. Эксцесс функции распределения приобретает вид
к4 =
k4 — 3k2 D4 — 3D2 670
k 22
D2 D2
„2 Ю2 2РТ - З70 '
а при условии уТ ^ 70 запишем формулу к4 = 3^0/иТ, совпадающую с выражением (6) и с формулой, полученной в работе [7] для первого приближения.
Следовательно, проведенное описание броуновского движения при воздействии пуассоновского случайного процесса позволило получить в общем виде выражения для характеристических функций флуктуаций скорости движения броуновской частицы и флуктуаций напряжения на электролитической ячейке, а также уточнить эти выражения для случая пуассоновских флуктуаций коэффициента вязкого трения. Показано, что эксцесс функции распределения для всех рассмотренных случаев обратно пропорционален интенсивности пуассоновского процесса, воздействующего на броуновскую частицу.
ЛИТЕРАТУРА
1. Морозов А.Н. Описание диффузии и броуновского движения как пуассоновских случайных процессов // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 1999. № 2. С. 85-90.
2. Бункин Н.Ф., Морозов А.Н.Стохастические системы в физике и технике. М.: Изд-во МГТУ им. Н.Э. Баумана, 2011. 366 с.
3. Калинкин А.В. Стационарное распределение для стохастической системы частиц, взаимодействующих комплексами // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2014. № 4. С. 3-17.
4. Марков Ю.Г., Синицын И.Н. Одно- и многомерные распределения флуктуаций неравномерности вращения Земли // ДАН. 2009. Т. 428. № 5. С. 616-619.
5. Пугачев В.С., Синицын И.Н. Стохастические дифференциальные системы. М.: Наука, 1990. 632 с.
6. Морозов А.Н. Метод описания немарковских процессов, задаваемых линейным интегральным преобразованием // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2004. № 3. С. 47-56.
7. Морозов А.Н.Стационарные распределения флуктуаций скорости броуновской частицы в среде с флуктуирующим коэффициентом вязкого трения // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2014. № 3. С. 26-38.
8. Морозов А.Н.Необратимые процессы и броуновское движение: Физико-технические проблемы. М.: Изд-во МГТУ им. Н.Э. Баумана, 1997. 332 с.
9. Коротаев С.М., Морозов А.Н., Сердюк В.О., Сорокин М.О. Проявление макроскопической нелокальности в некоторых естественных диссипативных процессах // Известия вузов. Физика. 2002. № 5. С. 3-14.
10. Морозов А.Н.Предварительные результаты измерений меры Кульбака флуктуаций напряжения на электролитической ячейке // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2011. № 2. С. 16-24.
11. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Наука, 1976. 576 с.
12. Никифоров А.Ф., Уваров В.Б. Специальные функции математической физики. М.: Наука, 1984. 344 с.
13. Малахов А.Н.Кумулянтный анализ случайных негауссовых процессов и их преобразований. М.: Сов. радио, 1978. 370 с.
REFERENCES
[1] Morozov A.N. Description of diffusion and Brownian motion as accidentai Poisson processes. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 1999, no. 2, pp. 85-90 (in Russ.).
[2] Bunkin N.F., Morozov A.N. Stokhasticheskie sistemy v fizike i tekhnike [Stochastic Systems in Physics and Technology]. Moscow, MGTU im.N.E. Baumana Publ., 2011. 366 p.
[3] Kalinkin A.V. Stationary Distribution for a Stochastic System of Complexes Interacting Particles. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2014, no. 4, pp. 3-17 (in Russ.).
[4] Markov Yu.G., Sinitsyn I.N. One- and Multidimensional Distributions of Fluctuations in the Earth's Irregular Rotation. Dokl. Akad. Nauk [Proc. Russ. Acad. Sci.], 2009, vol. 428, no. 5, pp. 616-619 (in Russ.).
[5] Pugachev V.S., Sinitsyn I.N. Stokhasticheskie differentsial'nye sistemy [Stochastic Differential Systems]. Moscow, Nauka Publ., 1990. 632 p.
[6] Morozov A.N. Method of Describing Non-Markovian Processes Defined by Linear Integral Transformation. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2004, no. 3, pp. 47-56 (in Russ.).
[7] Morozov A.N. Stationary Fluctuations Distributions of Brownian Particle Velocity in a Medium with Fluctuating Viscous Friction Coefficient. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2014, no. 3, pp. 26-38 (in Russ.).
[8] Morozov A.N. Neobratimye protsessy i brounovskoe dvizhenie: Fiziko-tekhnicheskie problemy [Irreversible Processes and Brownian Motion: Physical and Technical Problems]. Moscow, MGTU im. N.E. Baumana Publ., 1997. 332 p.
[9] Korotaev S.M., Morozov A.N., Serdyuk V.O., Sorokin M.O. Manifestation of the Macroscopic Nonlocality in Some Natural Dissipation Processes. Russian Physics Journal, 2002, no. 5, pp. 431-444.
[10] Morozov A.N. Preliminary Results of Recording the Kullback Measure of Voltage Fluctuations on Electrolytic Cell. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2011, no. 2, pp. 16-24 (in Russ.).
[11] Kamke E. Differentialgleichungen: Lösungsmethoden und Losungen, I, Gewohnliche Differentialgleichungen [Handbook of ordinary differential equations]. Leipzig, B.G. Teubner, 1977.
[12] Nikiforov A.F., Uvarov V.B. Spetsial'nye funktsii matematicheskoy fiziki [Special Functions of Mathematical Physics]. Moscow, Nauka Publ., 1984. 344 p.
[13] Malakhov A.N. Kumulyantnyy analiz sluchaynykh negaussovykh protsessov i ikh preobrazovaniy [Cumulant Analysis of Non-Gaussian Random Processes and Their Transforms]. Moscow, Sov. Radio Publ., 1978. 370 p.
Статья поступила в редакцию 23.09.2015
Морозов Андрей Николаевич — д-р физ.-мат. наук, профессор, заведующий кафедрой
"Физика" МГТУ им. Н.Э. Баумана.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул.,
д. 5.
Morozov A.N. — Dr. Sci. (Phys.-Math.), Professor, Head of Physics department, Bauman
Moscow State Technical University.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005
Russian Federation.
Просьба ссылаться на эту статью следующим образом:
Морозов А.Н. Описание флуктуаций скорости броуновской частицы при воздействии пуассоновского случайного процесса // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2016. № 1. C. 27-35. DOI: 10.18698/1812-3368-2016-1-27-35
Please cite this article in English as:
Morozov A.N. Fluctuations of the Brownian particle velocity influenced by a random Poisson process. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2016, no. 1, pp. 27-35. DOI: 10.18698/1812-3368-2016-1-27-35
