Кинетические явления в средах с микроструктурой как немарковские процессы Текст научной статьи по специальности «Математика»

НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н.Э. ЬАУМЛНЛ

НАУКА и ОБРАЗОВАНИЕ

Эл JVa ФС 77 - 48211. Государственная регистрация №(I4212Ü0025. ISSN 1994-0408

электронный научно-технический журнал

Кинетические явления в средах с микроструктурой

как немарковские процессы

# 05, май 2012

Б01: 10.7463/0512.0366037

Скрипкин А. В.

УДК 519.62

Россия, МГТУ им. Н.Э. Баумана [email protected]

1. Введение

Описание кинетических процессов проводят, как правило, с использованием дифференциальных уравнений второго порядка (например, уравнения диффузии) [1]. Случайные воздействия, оказываемые на кинетический процесс, в этом случае учитывают путем добавления к соответствующему дифференциальному уравнению случайной силы или источника [2]. При указанном подходе описываемый кинетический процесс будет относиться к классу марковских процессов [1].

Описание кинетических процессов в средах с микроструктурой, а также при наличии частиц или протяженных объектов, характерные размеры которых относятся к области микро- и нанометров, сопряжено с рядом особенностей. Кинетический процесс в таких случаях может проявлять ярко выраженные наследственные свойства [2]. В частности, диффузия микрочастицы вызывает соответствующее движение окружающих ее частиц среды, которое, в свою очередь, начинает оказывать влияние на движение самой микрочастицы. Характер этого влияния в данный момент, очевидно, зависит от всей предыстории движения частицы, которое становится теперь немарковским процессом.

Типичные немарковские процессы широко распространены в природе. Примером немарковского процесса может служить фликкер-шум, наблюдаемый в процессах, имеющих различную физическую природу [3]. В частности, экспериментально наблюдаемые флуктуации кинетических коэффициентов (например, флуктуации коэффициента электропроводности) имеют спектральную плотность, характерную для фликкер-шума [4].

Отметим также, что при воздействии марковского случайного процесса на динамическую систему её отклик представляет собой немарковский случайный процесс. Сумма двух марковских процессов представляет собой немарковский процесс. Немарковскими будут также процессы, образованные при интегрировании марковского или нахождении скользящего среднего от процесса с независимыми значениями [5].

Любая фильтрация (как линейная, так и нелинейная) приводит к преобразованию марковского процесса в немарковский. Результатом большинства динамических измерений выступает немарковский процесс, так как в таких измерительных системах всегда осуществляется фильтрация входного воздействия.

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

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

2. Стохастические дифференциальные уравнения

Кинетические процессы, происходящие в физических системах, традиционно описываются с помощью соответствующих дифференциальных уравнений. Так, при описании движения диффундирующей частицы обычно используется подход, основанный на применении стохастического уравнения Ланжевена [6]. Указанный подход позволяет воспользоваться хорошо разработанной теорией стохастических дифференциальных систем [1], с помощью которой можно определять все необходимые статистические характеристики флуктуаций скорости движения частицы.

Рассмотрим движение шарообразной частицы радиуса Я и массы М в среде (жидкости или газе) с кинематической вязкостью V и плотностью р. Уравнение движения такой частицы имеет вид

= ^+(0(/). (1)

ш

где V ^) - скорость частицы, Fc ) - сила сопротивления, ^) - сумма остальных внешних заданных сил.

Силу сопротивления Fc(t) = Fc0(t) для линейной вязкой среды обычно записывают в

виде

Ко^) = (г), (2)

где у - коэффициент вязкого трения, который для шарообразной частицы равен у = 6пр\Я. Считая спектральную плотность случайного процесса ^ ^) равной

Ом = 2уМкТ, (3)

где к - постоянная Больцмана, Т - температура вязкой жидкости, в которой находится броуновская частица, спектральную плотность флуктуаций скорости V (^) запишем как [6]

а

О"(И) = ^' (4)

где а = 2укТ/М , р = у/М .

Если случайный процесс ^ (г) представляет собой производную от процесса с независимыми приращениями, то скорость V(t) частицы описывается марковским случайным процессом, что позволяет определить для неё любые Ь -мерные характеристические функции, а, следовательно, и любые многомерные функции распределения [1].

Вообще говоря, теория стохастических дифференциальных систем дает возможность определять необходимые статистические характеристики случайного процесса 2(^) в том случае, когда этот процесс удовлетворяет дифференциальному уравнению

Ш2 = а(2, t )Ш + Ъ(2, t )ШЖ (г), (5)

где Ж(^) - процесс с независимыми приращениями. Процесс 2в этом случае является марковским процессом, статистические характеристики которого можно определить путем решения дифференциального уравнения для Ь -мерной характеристической функции

gL (V.., ь Ь; tl,..., ч) [1].

Но существуют процессы, которые невозможно описать с помощью дифференциального уравнения (5). В частности, если случайный процесс 2 (^), описывается с помощью линейного интегрального преобразования

z (t) = J G(t, t) dW (t), (6)

0

где G(t, t) - непрерывная функция переменной t , то этот процесс может не являться марковским процессом. Здесь считается, что интеграл (9) представляет собой интеграл Ито

[7].

Из выражения (6) следует, что начальное условие для процесса Z (t) имеет вид

Z (t )| t=0 = ^ (7) а, следовательно, его одномерная характеристическая функция g1 (Х; t) принимает начальное условие [1]

gi(X; t )| t=o = 1. (8)

Если интегральное преобразование (6) имеет ядро G(t, т), допускающее сведение уравнения (6) к конечномерной системе уравнений типа (5), то задача нахождения статистических характеристик процесса Z(t) может быть решена стандартными методами теории стохастических дифференциальных систем [7]. Но в общем случае такое преобразование оказывается невозможным, что приводит к необходимости разработки адекватной методики нахождения статистических характеристик процесса Z (t), описываемого уравнением (6).

3. Описание немарковского процесса, задаваемого линейным интегральным

преобразованием

Рассмотрим общий метод нахождения L -мерной характеристической функции 1,..., XL; t1,..., tL ) процесса Z (t), описываемого линейным интегральным преобразованием (6) [16, 20].

Начнем с построения одномерной характеристической функции. Будем считать, что преобразование (6) представляет собой интеграл Ито [7]. Разобьем интервал (0, t) на N равных интервалов, продолжительностью At = t/N . Моменты времени tn, соответствующие

окончаниям интервалов, удовлетворяют условию: t0 = 0 < t1 <... < tN_1 < tN = t. Тогда интеграл Ито в (6) можно заменить средним квадратическим пределом

N _1

Z (t) = limS G(t, tn )AW (tn), (9)

n=0

где AW (tn ) = W (tn+1)_ W (tn) - независимые приращения процесса W (tn).

Одномерная характеристическая функция g1 (Ь; t) процесса 2 (t) выражается формулой

й(Ь; t) = ^(/12(0)), (10)

где операция (...) представляет собой нахождение математического ожидания. Подставляя в

формулу (10) выражение (9) имеем

й(Ь;t)= ^ехр^/Ь Шт |1О(t,^)АЖ(^)| = ^Цт Пexp(;ЬО(t,^)АЖ(^. (11)

Так как приращения АЖ(^) являются независимыми, то последнее выражение в (11) приобретает вид

N-1

й (Ь; t) = Ит П (exp(/ЬО(t, ^ )А Ж (^ ))>. (12)

N ша

п=0

Логарифмирование выражения (12) дает

N-1

1пй(Ь;t)= lim | 1пА(Ь;t;^), (13)

п=0

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

¿(Ь; ^ ^) = (exp(/ЬО(t, tn )А Ж (^ ))>. (14)

Учитывая то, что для процесса с независимыми приращениями характеристическая функция приращений записывается через его одномерную характеристическую функцию, имеем

(15)

п" \ п >

А(Ь;,,) = {exp(/ЬО(/,,)АЖ(,))> =

Тогда

1п А(Ь; t; tn ) = ^(ДО^, tn М^ ))>-Ы^/АО^, tn Щя ))>. (16)

Подстановка выражения (16) в формулу (13) дает

1п (Ь. t )= 1{т I-1 1n(exp(гЬО(t, ^ )Ж (tn+l))) - ^(/ЯО^, ^ )Ж (tn )))^ (1?)

Nш« ^ Аt

п=0

После замены в формуле (17) суммирования интегралом Ито и последующего потенцирования получившегося выражения получим

Й1 (ь; t )= exp

где

г

т); т)Шт

(18)

х(ЬО^, т); т) = 1п А (ЬО(t, т); т), (19)

дт

А (ЬО(t, т); т) = ^(ДО^, т)Ж (т))). (20)

Для случая, если процесс Ж (х) является винеровским [1], имеем

Е

1 & х ) =

ехр

а если - пуассоновским [1], то

Е х ) = ехр

V

1 1

1О2 (х, т)0т

|(Е (Х0(х, т))-1)01

(22)

Аналогично для Ь -мерной характеристической функции можно получить

Еь ь ь ; х1,..., хь ) =

ехр

к

Л

где

X

Х10(х1, т)^; т] = -|Т 1п ¿1х, т^; т^

г=1

Г/

¿1 1^г0(хг, т) т

Л

г=1

(

= ехр

\ 1]Г х,0(хг, т)\Ж (т)

г=1

(23)

(24)

(25)

При нахождении интеграла в выражении (23) необходимо учитывать условие

0(х, т)|х> ^ = 0. (26)

Формулы (18) - (20) являются частными случаями выражений (23) - (25) при Ь = 1. В том случае, если процесс Ж (х) является винеровским, то

Еь (^1,..., ь ь хь )=ехР

1 Ь т1п (х1, хк

Ь к 10(хг, т)о(хк, т)0т

г ,к=1

(27)

а если пуассоновским, то

Еь (ь^..^ ь ь хь )=ехР

Л

От

(28)

Е х,, т)|-1

В общем случае, когда интегральное преобразование (6) не сводится к конечномерной системе дифференциальных уравнений, формулы (26), (27) и (28) описывают немарковский случайный процесс.

4. Движение микрочастицы в вязкой среде

При описании движения шарообразной микрочастицы с учетом увлечения частиц среды вместо формулы (2) необходимо использовать выражение [8]

2

М:р:/Ле^потаа^и.щ/doc/366037.htm1

193

Fe(t) = -2прЯ3

1 ШV(t) + Зу V(t) +1 V (т) Шт 3 Л Я2 Я\ п 0 Шт л/Т-т

(29)

Здесь момент времени Т = 0 принят за начало движения частицы. Подставляя (29) в (1), получим уравнение

' М+2 ПРЯ3 ] ШШ

V 3 ) ж

+ 6прvЯV (Т) + 6рЯ2 л/Лу |

ШV(т) Шт Шт д/т-т

= Fо(t) + ^ Ц).

(30)

Введя замены

ШV(т) л_ 6пруЯ ¡1 Fо(t) , ~} = ^у(т)

м + -прЯ3 м + 2 прЯ3 м + 2 прЯ3

2(т) = ^, Л = "У" , 5 = Я,/—, F0(t) = 1 у , ^(Т) = ^ , (31)

Ш ' прЯ 3' м +2 прЯ3' м + 2

3 3 3

придем к выражению

^ 5

2 (Т) + ^ + (т)Шт = Fо(t) + ^), (32)

в котором учтено, что

V (Т) =1-^ Шт, (33)

о Шт

т.к. V (0) = 0. Выражение (32) описывает случайный процесс 2(Т), являющийся в этом случае немарковским случайным процессом [17]. Таким образом, использование для силы сопротивления выражения (29) вместо (2) приводит к необходимости применения интегральных уравнений, а, следовательно, и теории немарковских процессов [2]. Решение интегрального уравнения (32) имеет вид [20]

2 (Т) = } (8(Т -т) - Я(Т, т) )(~0 (т) + )Шт, (34)

0

где резольвента

Я(Т, т) = | Кп (Т, т). (35)

п=1

В последнем выражении члены ряда определяются с помощью рекуррентного соотношения

Кп+1(Т , т) = ¡К^, s)Kn ^, тШ, (36)

т

где

Г В Л

ЗД т) = Л|^1 + (37)

Вычисление по формуле (36) приводит к выражениям для первых членов ряда

К2 (Т, т) = - Л2 ((Т - т) + 4В(Т - т)1/2 + пВ2) (38)

к 3(х, т) = Л2

1 (х - т)2 + 4в(х - т)3/2 + 3 пВ2(х - т) + 2 пВ3 (х - т)

1/2

(39)

Из выражений (37) - (39) видно, что ряд (35) является знакочередующимся, причем резольвента Я(х, т) является функцией только от разности аргументов: Я(х, т) = Я(х -т). Видно также, что абсолютные значения членов ряда (за исключением первого) возрастают при росте х - т таким образом, что |КИ (х - т)| ^ да при (х - т) ^ да. Для не слишком больших х (х < 10с ) и А << 1 с-1 ряд (35) с большой степенью точности можно заменить суммой

Я(х -т) = Л -пЛ2В2 +-Л== - 4ВЛ'л/х-г.

л/х-т

Будем далее считать, что ^0(х) = 0. В этом случае (34) примет вид

(40)

7 (х) = |[5(х - т) - Я(х - т)]£(т)0Г.

(41)

Используя метод описания немарковских случайных процессов, изложенный выше, для одномерной и ь-мерной характеристических функций случайного процесса 2(х), задаваемого интегральным уравнением (41), для не слишком больших х (х < 10с ) и Л << 1 с-1 получим

Е1(^;х) = ехр

-1 аА2 2

(

л2в21п--4лвс4! + (с2 -8Л2в2)х+—Л2ВСл/х7 + 8л4в2х2

5х 3

Еь А ь; х1,..., хь) = ехр

ь ь

-а^^Аг А (х,, хк)

г=1 к=1

(42)

(43)

Здесь а - интенсивность случайного процесса £(х), 5х - малая положительная величина, а функция Qгk (хг, хк ) имеет вид

Ош (хг, хк ) = 2 С X - ЛВС + & -^^Ч )+ + 4 Л2ВС(д + ф! -v(х, - хк )3)- 4л3В2 ^ +

(44)

+ Л2В2(1 + 8Л2(х, -хк))1п

у[хк +л1х~г

- хк + 5х + л/57

+ 4 Л4 В2 л/хкхг3.

В выражениях (42) и (44) С = - Л + пЛ В2.

Многомерные характеристические функции еь (А1,..., Аь; х1,..., хь), задаваемые выражением (43), позволяют найти моменты любого порядка случайного процесса 7(х) [1]. В частности, для корреляционной функции получим

(45)

7 ($¿2 (х2)) = 25О12(х1, х2).

Для нахождения спектральной плотности случайного процесса 2(Т) найдем преобразование Лапласа уравнения (34) при F0(t) = 0, что позволяет записать его в изображениях

Р)

(

1 + Л

1+в^-П

р 4Р

\

=« р)

(46)

где 2(р) и р) - изображения функций 2(Т) и ) соответственно.

Ввиду того, что спектральная плотность О~ = о процесса ) постоянна и равна

О~ =

12пруЯкТ М + 2 прЯ3

(47)

то в соответствии с формулой (46) спектральная плотность процесса 2(Т) при Т ш« принимает вид

/ш

О2 (ш) =

/ш + ЛВл/П • /ш + Л

О~

или

О2 (ш) =

ш

ш2 +.

О~.

(48)

(49)

- ЛВл12пш3 + Л2 В2 пш + Л2 + Л2 ~

Из последнего выражения может быть найдена спектральная плотность флуктуаций скорости броуновской частицы V(). Согласно выражению (31) спектральная плотность

О2 (ш)

° (ш) =

или

ш

Оv (ш) =

О~

ш2 +.

(50)

- ЛВ^ 2пш3 + Л2 В2 пш + Л2 Вл/2лш + Л2 Сравнение формул (50) и (4) показывает, что использование для силы сопротивления выражения (29) вместо (2) приводит к существенному различию спектральных плотностей флуктуаций скорости частицы, в особенности, в полосе низких и средних частот. На рис. 1 показаны в сравнении графики спектральных плотностей, даваемые формулами (4) и (50).

2

Рис. 1. Графики спектральных плотностей, задаваемые формулами (4) (кривая 2) и (50) (кривая 1) при р = 1000 кг/м3, V = 10-6 м2/с, Я = 10-4 м, M = 4-10-9 кг, T = 300 К

Пусть теперь на сферическую частицу вдоль оси Х кроме случайной силы и силы сопротивления действует возвращающая сила /^(х) = -кХ(х). В этом случае уравнение (32) примет вид

7 (х) + Л Л 1 + -В= 2 (т)От = -кХ (х) + ~(х),

4Г-

где

к =

к

2

М + -прЯ3 3

(51)

(52)

Найдем в этом случае спектральные плотности процессов 7(х) и ¥(х) при х ^ да. Учтя, что преобразование Лапласа функций Х(х) и 2(х) связаны соотношением Х(р) = 7( р)/ р2, для уравнения (51) в изображениях имеем

Р)

, к лв4Л л

1 + —т +-+ —

=,( р).

(53)

р 4Р р

Выполнив процедуру, проведенную при выводе формул (49) и (50), для спектральных плотностей процессов 7 (х) и У(х) получим

07 (Ш) =

ш

4

ш4 +

ЛВ>12пш7 + Л2В2пш3 + Л2Вл12пш5 + (Л2 - 2к)ш2 - кЛВ^12пш3 + к

0Е

2 5

(54)

(ш) =

ш

4

ш4 +

ЛВл12пш7 + Л2В 2пш3 + Л2Вл12пш5 + (Л2 - 2к )ш2 - кЛВ^] 2пш3 + к

О,.

2 5

(55)

Из (55) легко найти спектральную плотность для Х(Т) в виде

ОХ (ш) =

1

3 , 7,2 \

ш4 + лв-42пш7 + л вУЩ + л в]2пш5 + (л - 2Л:)ш2 - клв^2пш3 + к: V л/ „ V

Сравним последнее выражение со спектральной плотностью классического

1/2

осциллятора, которая получается из (56) при условии, что В = 0 с , а в формулах (31) (для

2

коэффициента Л), (47) и (52) слагаемое —прЯ3, стоящее в знаменателе, обращается в нуль.

В этом случае (56) переходит в формулу

О™с (ш) =

ш4 + (Л'2 - 2к')ш2 + к'

2 +^2 %

(57)

,, 6пруЯ ,, к 12пруЯкТ

где теперь Л =-, к =—, О~ =-;-.

М М 5 М2

На рис. 2 и 3 изображены спектральные плотности, задаваемые выражениями (56) и

(57) при у = 10-6 м2/с, р = 103 кг/м3, к = 10-2 Н/м, Я = 10 мкм (рис. 2), Я = 100 мкм (рис. 3).

Плотность частицы при этом принимается равной плотности среды.

Рис. 2. Графики спектральных плотностей, задаваемые выражениями (56) (кривая 1) и (57) (кривая 2) при у = 10-6 м2/с, р = 103 кг/м3, к = 10-2 Н/м, Я = 10 мкм

1

Рис. 3. Графики спектральных плотностей, задаваемые выражениями (56) (кривая 1) и (57) (кривая 2) при V = 10-6 м2/с, р = 103 кг/м3, к = 10-2 Н/м, Я = 100 мкм

Из графиков видно, что с увеличением размера частиц формы кривых становятся похожими, при этом классическому случаю соответствует большее амплитудное значение спектральной плотности. Этот эффект становится наиболее заметным при больших размерах частиц. Видно также, что максимум спектральной плотности в классическом случае соответствует более высокой резонансной частоте. Этот факт является следствием того, что при неклассическом описании к массе частицы добавляется некоторая «эффективная» масса 2

— прЯ. 3

Таким образом, проведенное описание движения броуновской частицы шарообразной формы в неограниченной вязкой среде позволило установить, что флуктуации её скорости представляют собой немарковский случайный процесс. Полученные резонансные кривые для механического осциллятора, размещенного в вязкой безграничной среде, по своей форме отличается от классической. Последний результат может иметь существенное значение для устройств демпфирования колебаний.

5. Движение частицы на малых временных интервалах

Для рассмотрения движения микрочастицы на малых временных интервалах, на которых ее взаимодействие со средой оказывается весьма существенным, проведем преобразование Лапласа исходного уравнения (32).

Перед проведением преобразования Лапласа перепишем выражение (32), представив его в виде

м

ду (г)

сИ

+

ь

¡Г« -т)

су (т) Ст

Ст = 7 (Г) + (Г),

(58)

где функция у(х - т) имеет смысл ядра коэффициента трения и определяется выражением

4 -зг^ , 6рЯVпу

у(г -т) = 6пруЯ + —лрЯ5(г -т) + ,- •

3 -т

Преобразование Лапласа соотношения (58) приводит к равенству

Мр2Х(р) + у(р)р21(р) = Р(р) +у (р).

(59)

(60)

где р - параметр преобразования, Х(р) = ¿(р)/р2, р), (р), /(р) - преобразования Лапласа координаты броуновской частицы X(г), внешней детерминированной силы Я0(г), случайного воздействия %(г) и ядра у (г) соответственно. Полагая далее движение броуновской частицы свободным, при котором внешняя сила 7(г) равна нулю, из (60) найдем

4 (р) (61)

Х( р) =

р (М + у(р))

Последнее выражение позволяет записать формулу для нахождения зависимости координаты броуновской частицы X (г) от времени в виде свертки

I

X (г) ^ Ох (г-т)%у (т)Ст,

0

где для преобразования Лапласа искомого ядра ОХ (г) имеем выражение

(62)

ОХ ( р) -;-,

р (М + у(р))

в котором согласно (59) преобразование функции у (г) определяется как

ч 4 л3 6прЯ24у 6пруЯ у(р) = -прЯ + — +--

3 4р р

(63)

(64)

Рассматриваемому интегральному преобразованию (62), как следует из вышеизложенного, соответствуют следующие одномерная g1(Л; г) и многомерная

gL (Л,..., Л; г1,..., гь ) характеристические функции

gl (Л;г ) = ехр

2^Л2 ¡ОХ (г -т) Ст

(65)

gL (Л... Л; tl,..., tL ) = ехр

1

—V

£ЛЛ ¡ Ох (г1 -т)Ох (Ч -т)Ст

1 ,к=1

(66)

Здесь vw - интенсивность внешнего случайного воздействия. Из (63) и (64) найдем

ър) = ( ! г ) • (67)

ар I р + Оу]р + с I

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

,,4 П3 , 6прЯ 24у 6прvЯ (68)

а = М + -прЯ\ Ь =—--, с = —-----(68)

3 а а

Выражение (67) не позволяет в общем случае найти обратное преобразование Лапласа. Однако, если Ь > 2у[с , то (67) может быть представлено в виде

'У+.Ж + У + У 4 ^

11 ар[р+Ь4р+с)=а [ р+Тр +7р-Х"+Тр

Х2 У

(69)

где

х1 =1 (-Ь + 4ьГ-4~с), х2 =1 (-Ь Ь2 - 4с ), (70)

1 Ь Ь2 + Ьл/Ь2 - 4с - 2с - Ь2 + Ьл/Ь2 - 4с + 2с

У1 = ", У2 = -—, У3 =-2 I 2-, У4 =-2 / 2-. (71)

с с 2сЧЬ2 - 4с 2с2ЫЬ2 - 4с

Заметим, что необходимое условие Ь > удовлетворяется для широкого диапазона

частиц и сред, соответствующих реальным условиям. В частности, при движении частицы

радиусом Я = 1 мкм в воде (при нормальных условиях плотность р = 1000 кг/м ,

кинематическая вязкость V = 10-6 м2/с) неравенство Ь > 2у[с удовлетворяется при М < 5,2 -10-16 кг (легкие пористые частицы). В дальнейшем в качестве тестовой частицы (имея в виду последующее сравнение с результатами экспериментальных исследований) будем брать частицу плотностью р0 = 100 кг/м и радиусом Я = 0,1 мкм, движущуюся в воде.

Обратное преобразование функции (67) при учете (69) приводит к следующей формуле для ядра Ох (г - т) исходного преобразования (62):

Ох (г -т) =1

а

У2 + у3 + у4 у1 + 2, +

д/п(г -т) (72)

+ х1 у3 ехр(х12 (г - т))егГс(- г -т)+ х2 у4 ехр(х22 (г - т))егГс(- х24Т-х)

Здесь функция егБс(х) представляет собой дополнительный интеграл ошибок. Из определения переменных у2, у3 и у4 следует, что у2 + у3 + у4 = 0. Тогда окончательно получим

Ох (г - т) = 1 [у + х Уз ехр(х1 (г - т))ег&(- х л/Г-т)+

+ х2у4 ехр(х22(г -т))ег&(- х2л/г-т)]. Найденное выражение (73) принципиально позволяет определить средний квадрат смещения частицы ^х2(г)^ согласно формуле

X 2(г )\ = д 2 ^) , (74)

^ ' д(^)2 ,=0

где характеристическая функция g1(Л; г) задается выражением (65).

При больших г - т, в связи с тем, что функция егБс(х) убывает быстрее возрастания функции ехр(х2), ядро Ох (г - т), как следует из (72), сводится к постоянной величине

Ох (г - т)|г-т^ш = ^ (75)

В этом случае из соотношений (65) и (74) получим

„2, л у,2 1 2 у0кТ 2кТ

Х (г) г = V,Лг = Vг =г =-г = 2^г, (76)

1 г а ас у 0 у 0

что совпадает с известной формулой Эйнштейна.

Для приближенного нахождения величины ^х2(гна малых временных интервалах

(порядка 10-3 -10-2 с) в выражении (73) можно разложить функции ехр(х) и ег&(х) в ряд. Ввиду того, что величины х1 и х2, вообще говоря, принимают большие значения (для тестовой частицы ~104 с"12), то следует удержать значительное число слагаемых разложения функций ехр(х) и егГе(х). Оставив 30 слагаемых ряда, можно найти, что для значений

г - т < 10-3 с функция Ох (г - т) хорошо аппроксимируется зависимостью

\0,5

О^-)!,-„„ = . (77)

3,5 • 10(г-т)0 а

Использование найденной аппроксимирующей зависимости для ядра Ох (г -т) позволяет найти, после подстановки выражения (77) в формулу (65) и последующего применения (74), что для малых г среднее значение квадрата смещения броуновской частицы задается квадратичной по времени функцией

(х2 (г }) = Шаг2, (78)

где параметр а « 210 с"1.

Проведенные в работах [9, 10] экспериментальные исследования броуновского движения частицы (параметры которой близки к рассматриваемой нами выше тестовой частицы) в вязкой среде на малых временных и пространственных интервалах показали, что для времен порядка 10-3 -10-2 с средний квадрат смещения хорошо описывается полученной выше зависимостью (78) (при этом в экспериментах, в которых частица двигалась в воде, получалось значение а = 180 с"1), а на больших временных (более 10-1 с) и пространственных интервалах удовлетворительной является классическая формула (76), что также следует из рассматриваемой выше теории.

Хорошая согласованность полученных результатов, базирующихся на формуле (29) для силы сопротивления, с соответствующими экспериментами на малых и больших временных интервалах, говорит о правомерности использования выражения (29) при учете взаимодействия движущейся микрочастицы и частиц среды, приводящего к увлечению среды микрочастицей.

Таким образом, проведенное описание движения броуновской частицы шарообразной формы в неограниченной вязкой среде, учитывающее ее взаимодействие с окружающими частицами среды, позволило установить, что флуктуации её скорости представляют собой немарковский случайный процесс. Найденные статистические характеристики броуновского движения имеют существенное отличие от аналогичных характеристик, получаемых при классическом описании броуновского движения, использующего методы теории стохастических дифференциальных систем. Наиболее ярко отличие в движении броуновской частицы проявляется на малых временных и пространственных интервалах, что обусловлено принципиальной необходимостью учета на таких интервалах взаимодействия броуновской частицы и частиц среды.

6. Диффузия в пространстве, окружающем микрочастицы жидкости.

Испарение микрокапель

Рассмотрим неподвижную каплю жидкости сферической формы с первоначальным радиусом Я0, вокруг которой происходит диффузия частиц пара (рис. 4). Начало сферической системы координат поместим в центр капли. Концентрацию насыщенного пара жидкости капли обозначим через п0, а концентрацию пара у поверхности будем считать

некоторой произвольной функцией времени и обозначим как пЯ (г). Также будем предполагать коэффициент диффузии Б постоянным. Массу одной частицы жидкости

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

Очевидно, ввиду симметрии рассматриваемой задачи концентрация частиц пара будет зависеть только от расстояния до центра капли г, а также времени г. Уравнение диффузии в сферических координатах для такого случая принимает вид

дп(г, г) Б д2(гп(г, г))

дг

дг2

(79)

Граничные и начальные условия записываются согласно сказанному выше с помощью выражений

п(г,г )| г=Я = пя (г),

(80)

п(г, г )| г=0 = По. (81)

Заметим, что в условии (80) радиус частицы Я является величиной, в общем случае зависящей от времени.

Рис. 4. Диффузия пара над поверхностью частицы сферической формы

Принимая за положительное направление потока частиц с поверхности q(г) направление вдоль координаты г, и считая, что этот поток пропорционален разности концентраций пЯ (г) и п0, получим

q(г) = -а(пя (г) - по ) + ^ ч (г).

(82)

В последнем выражении учтены флуктуации потока частиц, выражаемые последним слагаемым ^ (г) . Величина а - коэффициент аккомодации.

Очевидно, поток испаряющихся частиц q(г) удовлетворяет также общей формуле

/ Ч тлдп

q(г) = - Б — дг

(83)

г=Я

Для потока массы qm (г), равного количеству массы частиц, испарившихся (конденсировавшихся) с единицы площади в единицу времени, получим

qm (г) = -та(пк (г) - «0 ) + (г), (84)

/ Ч тлдп

qm(г) = -тБ—

дг г=Я

В выражении (84) случайный массовый поток ^ (г) = т^ (г) . Приравнивая соотношения (84) и (85), получим

(85)

- тБ — дг

= -та(пк (г) - «0 ^^ (г). (86)

г=Я

Рассмотрим прежде всего задачу в предположении постоянства радиуса частицы: Я(г) = Я . Такое приближение физически обосновано для случая частиц большого радиуса, а также при относительно коротком времени наблюдения процессов диффузии, в течение которого изменением размеров капель можно пренебречь.

Решим систему уравнений (79) - (81). Для этого введем вспомогательную функцию Г (г, г) , определяемую выражением

п(г, г) = ^^ . (87)

г

Подстановка последнего выражения в (79) - (81) приводит к системе уравнений для функции Г(г,г)

дГМ = Б^ГМ (88)

дг дг2

Г (г, г )| г=я = Я0«я (г), (89)

Г (г, г )| г=0 = «0 г, г > я,. (90)

Из записанных выражений (88) - (90) видно, что функция Г (г, г) формально соответствует одномерному уравнению диффузии. Поэтому решение поставленной задачи имеет вид [11]

(г - Я,)2

Г (") = I Т^ехр

4Б(г - т)

пЯ {т)ёт. (91)

24ю 0 (г -т)3

Заметим, что в записанном соотношении (91) уже учтено, что между потоком, связанным с начальной концентрацией частиц пара, и потоком испаряющейся жидкости почти сразу устанавливается динамическое равновесие, поэтому слагаемое, связанное с наличием начальной концентрации, в (91) опущено.

Найдем производную

дГ (г, г)

дг

принимая во внимание соотношение (85). Получим

г=Я

дГ (г, г)

дг

А_

тБ

qm (г) + ^ (г)

(92)

Производная последнего выражения по координате г определяется из уравнения (91) и

равна

дГ(г, г)_ Я

дг

1

(г -т)3/2 2Б(г -т)

(г - Я0)2

5/2

ехр

(г - Я0)2 4Б(г -т)

пЯ (т)йт.

Вычисление полученного интеграла (93) по частям дает

дГ (г, г) _ Я

дг

1 .171

4пй 0 Vг-т

^ехр

(г - Я0)2 4Б(г -т)

йпЯ(т)

ёт

йт.

Сравнивая (92) и (94), получим выражение для массового потока qm (г) в виде

(93)

(94)

. . ти г

qm (г) = I

тБ \ ( Я 1

Я * ^л/пБ л/г-т Введем замены с помощью формул

5п(г) = пЯ (г) - п0,

+1| йпя(т> йт.

] йт

г (г) =

йЬп(г)

йг

(95)

(96)

(97)

С использованием обозначений (96) и (97), из выражений (84) и (85) окончательно получим

ГI Б 1 Б

т11 л--,-+а+--

01 V п Vг-т Я

л

г (т)йт = Е (г).

(98)

Таким образом, случайный процесс г (г) описывается с помощью линейного интегрального оператора Вольтерра первого рода, ядро которого имеет вид суммы слагаемого абелевого типа, постоянной величины и величины, зависящей от радиуса сферической частицы. Описание случайного процесса г(г) с помощью интегрального оператора вида (98) указывает на немарковский характер рассматриваемого процесса. Укажем также на то, что величина 5п(г) также в общем случае является немарковским случайным процессом, так как представляет собой интеграл по времени от немарковского процесса г (г).

Преобразование Лапласа уравнения (98) приводит к уравнению

г=Я

Г О а О )

Л — + — +-

р р Я0 р )

^(P) = £qm (P),

где p - параметр преобразования. Из последней формулы находим

P

(99)

Др) =

m(Л/Dp +р)

* ^ (P),

где введена замена

О

В = а +— .

Я,

Воспользовавшись определением (97), получим

т(р) = ■ 1

(100)

(101)

^ (р).

(102)

т^Ор +р) т

Так как интеграл в левой части (98) представляет собой свертку, то решение для функции 5п(г) будем также искать в виде интегрального оператора Вольтерра первого рода, имеющего вид свертки

г

(103)

5«(г) = | К (г -т)^ (т)А.

Здесь К (г - т) - некоторая функция, которую следует найти.

В связи с тем, что массовый поток qm (г) и величина 5п(г) связаны соотношением (84), из (103) получим

г

(104)

qm (г ) = | 7 (г -т)^ (т)А,

где

7 (г -т) = 5(г - т) - аmK (г -т). (105)

Выполняя преобразование Лапласа выражения (103) и приравнивая полученный образ 5«(р) выражению (102), для образа ядра интеграла (103) К (р) получим

к (р ) =

1

(106)

mQDp+в).

Выполнив теперь обратное преобразование Лапласа выражения (106) [9, с. 210], будем иметь

К (г-т) =

1

1 -_Р_

\т]п(г -т) 4О

ехр

в2(г -т)

О

егБс

Рл/Г-

т

4О

(107)

где егТе(х) - дополнительный интеграл ошибок. Ввиду того, что функция егТе(х) убывает быстрее с ростом х по сравнению с ростом величины ехр(х2), то при (г -т) ^ го ядро К (г - т) стремится к нулю.

Рис. 5. График функции К (г-т), задаваемый выражением (107)

График функции К (г -т) схематически представлен на рис. 5. При этом принималось, что значения соответствующих параметров (в данном случае - массы частицы т, коэффициента диффузии Д коэффициента а и радиуса сферической капли Я0) равны единице. Это также было обусловлено тем, что быстрый рост экспоненты в (107) и одновременное очень быстрое убывание функции ег1Ъ(х) затрудняло расчеты с использованием реальных параметров. Однако, очевидно, что для таких параметров функция К (г - т) будет вести себя аналогичным образом.

Аппроксимация построенного графика степенной функцией с хорошим приближением дает для значений (г -т)> 0,05 с обратно пропорциональную зависимость функции К (г - т) от разности г - т:

К (г-т)=006

г-т

(108)

Применяя метод, изложенный ранее, для одномерной g1 (А; г) и ¿-мерной gL(А1,...,АL;г1,...,гь) характеристических функций случайного процесса 8п(г) получим соответственно выражения

gl

-(А; г ) =

ехр

1 '

- аА21К2 (г -т)й?т

(109)

gL А L ; ^^ гL ) =

= ехр

(

1

--а 2'^к

к=1 0

V

2

2 4 | К2 (гк - т^т + 2 2 А к А, | К г - т)К (г, - т)Л

| -"Лч к,1=1, 0 к <1

(110)

Рис. 6. Графики функции g1(А; г) при г = 5-10-6 с (кривая 1) и г = 1 с (кривая 2)

Схематический график функции g1 (А; г) для различных г изображен на рис. 6. Из него видно, что одномерная характеристическая функция g1 (А; г) имеет вид гауссовой кривой, ширина которой уменьшается с увеличением г.

Полученные характеристические функции позволяют найти любые статистические характеристики случайного процесса 5п(г) . В частности, для его математического ожидания

5п(г)) и момента второго порядка (8п(г1)8п(г2)) имеем выражения

5п(г)):

(А; г)

/ЗА

= 0,

А=0

8п(г1 )5п(г2 )) = -д 2 ^ (А-А 2; г- г2)

/ЗА1ЗА 2

А1=0, 0 А2=0 0

а| К(г1 - т)К(г2 - т)ёт.

(111)

(112)

В соотношении (112) считается, что гу < г2.

Последняя формула позволяет найти дисперсию процесса 5п(г), которую обозначим через Д(г) . В этом случае г = г1 = г2, поэтому

д(г) = (§п2 (г)} = а | К2 (г - т)^т.

(113)

Статистические характеристики флуктуаций потока массы частиц qm (г) задаются, очевидно, формулами, аналогичными (109) - (113), в которых ядро К (г - т) заменено на ядро

7 (г -т).

Заметим, что интегрирование в полученных выше интегралах проводится фактически в пределах от нуля до г - 5г, где 5г - малый положительный параметр, равный по порядку величины времени свободного пробега частиц пара. Введение этого параметра физически обосновано тем, что не имеет смысла говорить о взаимодействии частиц для меньших времен, следовательно интервал (г - 5г; г) не оказывает влияния на статистические характеристики соответствующих случайных процессов. Отсутствие интегрирования на интервале (г - 5г;г) исключает рассмотрение области сингулярного поведения интегралов при (г - т) ^ 0.

Найдем теперь одномерную функцию плотности вероятности р(5п) флуктуаций 5п(г), используя найденное выражение для одномерной характеристической функции ^(А; г) (4.74). Согласно определению [1] , запишем

1 ад

р(5п) = — | & (А; гУ5«ЫХ, (114)

-ад

то есть

1 ад 1 2

р(5п) = — | ^ )Х в''5«Х ёА. (115)

2П -ад

Интегрируя последнее выражение, получим [12, с. 344]

р(5п) = , 1 ехр -—«— . (116)

^ ' Т^паС) Я 2А(г)) ' '

Таким образом, флуктуации 5п(г) подчиняются нормальному распределению Гаусса, зависящему от времени. Графики функции р(5п) при различных г изображены на рис. 7. Видно, что с течением времени график плотности вероятности р(5п) «размывается» по оси 5п, а при малых значениях времени флуктуации величины 5п(г) сконцентрированы преимущественно в области малых значений 5п , стремясь к 5 -функции при г ^ ад .

0,0006

0.0005

0,0004

0.0003

0,0002

0,0001

0.0000

/7(5«), м2

/ / * / s \ \\ \\

\\

4 \\ \\ \ \

/ - ■ --— / i \ \ s

Ьп,

м'

-4000

-2000

2000

4000

Рис. 7. Графики функции p(dn) при t = 5-10-6 с (кривая 1) и t = 1 с (кривая 2)

Воспользовавшись соотношением (102), найдем спектральную плотность флуктуаций величины 5n(t)

(117)

°Ъп(ш) ш2+ ^2Бш+р2) ' Из последней формулы видно, что при ш ^ 0 значение спектральной плотности 0Ъп (ш) стремится к постоянной величине

GSn (о)

а

а

™2г>2

m 2р2

Í

m

а + -

D

R

2 '

(118)

Для спектральной плотности потока массы G (о) получим выражение

Gm (®) =

R2 Do+ R р/ 2 Do + D2 R2 (Do + pV 2 Do +p2)

а.

(119)

С, (со), кг2слм"

610 510 410 310 2-10 1-10

-8

-й

-8

2

1

ш, С

О 0,5-104 1,0-104 1,5-104 2,0-104 2,5-104 3,0-104

Рис. 8. Графики спектральной плотности массового потока (ш), соответствующие выражению (119), для воды (кривая 1) и этилового спирта (кривая 2)

Формула (119) показывает, что значение спектральной плотности Оч (ш) при ш — 0 стремится к постоянной величине

( ч О2а О2ст

Оа (ш) = —тт = 7-^г, (120)

^ ^> Я2в2 (Я>а + О)2

а при ш —> ад спектральная плотность массового потока Оч (ш) принимает значение мощности флуктуаций величины ^ (г)

Мш) =о. (121)

Графики, соответствующие выражению (119), для воды и этилового спирта изображены на рис. 8.

Заметим, что статистические характеристики для величин 5п(г) и qm (г) в случае испарения с плоской поверхности жидкости в полупространство находятся из полученных выше выражений при Я0 = ад .

7. Испарение микрокапли. Случай изменяющегося радиуса

Рассматриваемый до сих пор квазистационарный случай, при котором пренебрегалось изменением радиуса капли, очевидно, может быть использован лишь как первое приближение. Поток частиц через поверхность капли приводит к изменению массы капли и

ее радиуса. Особенно важно учитывать изменение радиуса в случае маленьких (субмикронного размера) капель, так как даже в течение небольшого промежутка времени их радиус может измениться значительно относительно первоначального. Ясно, что квазистационарное приближение в таком случае окажется неприменимым. Изменение радиуса капли необходимо учитывать при длительном наблюдении над ее поведением, даже в случае крупных (десятки и сотни микрон) капель аэрозоля. Наконец, при расчете времени «жизни» капли (до ее исчезновения в результате «схлопывания») учет изменения радиуса частицы, очевидно, необходим принципиально.

Однако решение задачи об испаряющейся капле жидкости при учете изменения ее радиуса осложняется ее математической сложностью (поставленная задача, как известно, относится к так называемому классу задач Стефана, аналитическое решение которых известно в очень ограниченном числе случаев [13]). В связи с этим получение статистических характеристик физических величин, описывающих поставленную задачу, может быть осуществлено только с использованием численных методов.

При проведении численного расчета использовались формулы, полученные для квазистационарного случая, при учете изменения радиуса капли. Изменение радиуса капли за время ёг можно найти, зная массовый поток qm (г). В самом деле, за время ёг масса капли М изменится на некоторую величину ёМ, определяемую равенством

ём = -аm (г) • 4пЯ2 (г) • ёг. (122)

С другой стороны общую массу капли можно определить согласно выражению

М = р- 4 пЯ3 (г), (123)

где р - плотность жидкости. Находя производную ёМ из последней формулы и

ёг

приравнивая ее той же производной из (4.87), получим

ёя(г)_ qm (г)

(124)

ёг р

Полученная зависимость радиуса капли от текущего массового потока через ее поверхность использовалась на каждом шаге итерации при численном расчете.

При численном расчете количество шагов итерации было равным 30000, шаг по времени Дг = 0,1 с, масса молекул, составляющих каплю, m = 10-22 кг, коэффициент

диффузии О = 2 • 10 5 м2/с, начальный радиус частицы Я, = 50 мкм, плотность жидкости

3 23 3

р = 1000 кг/м , концентрация насыщенного пара п0 = 10 м" .

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

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

При численном расчете был выбран также относительно большой шаг итерации по времени, что требовалось необходимостью устойчивости расчета при росте шага итерации. Однако, и в этом случае статистические характеристики рассматриваемой модели не претерпевают принципиальных изменений. В частности, они сохраняют свой вид при кратном изменении шага по времени.

Рис. 9. Характерный вид измененения радиуса капли Я с течением времени

Рис. 10. Характерный вид массового потока дт () в зависимости от времени

Рис. 11. График спектральной плотности флуктуаций массового потока 0„ (о)

На рис. 9 показан характерный график зависимости радиуса частицы от шага итерации к, а на рис. 10 изображен соответствующий график массового потока для этого случая. График на рис. 11 показывает спектральную плотность массового потока Од (о),

изображенного на рис. 10. Найденное уравнение кривой аппроксимации для построенного графика, подчиняющееся степенному закону, имеет вид

Оа (о) = 10-7о-103.

Чт ^ '

(125)

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

По результатам 3000 независимых реализаций были найдены те моменты итерационного процесса, при которых происходит исчезновение капли. На рис. 12 показан график распределения количества частиц Му (к), исчезнувших в к-ый шаг итерации. Видно,

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

Рис. 12. График распределения количества частиц N (к), исчезнувших на к-ом шаге

итерации

Полученные результаты могут иметь существенное значение при исследовании капельных аэрозолей, находящихся в атмосфере. Рассмотренная модель, учитывающая флуктуации массового потока через границу раздела жидкой и газообразной сред, может быть использована при метеорологических исследованиях, например, при прогнозировании времени исчезновения тумана; при экологическом мониторинге загрязнения городов и промышленных территорий объектами, находящимися в состоянии аэрозоля; при исследовании химических взаимодействий веществ, одно из которых находится в аэрозольном состоянии и др.

8. Теплопроводность вокруг сферической микрочастицы

Рассмотрим неподвижную сферическую частицу радиуса Я, температуропроводности Хм и объемной теплоемкости Су, в центр которой поместим начало сферической системы координат (рис. 4, в котором теперь вместо концентрации следует брать температуру).

Температуру поверхности частицы будем считать некоторой функцией времени ТЯ (^). Среду вне сферической частицы (при г > Я) считаем однородной с постоянными параметрами: плотностью р, теплопроводностью к и температуропроводностью х, причем X << хм. Начальная температура во всем пространстве равна некоторой постоянной

величине Т0. Изменением радиуса сферической частицы вследствие изменения температуры пренебрегаем.

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

дТ(г, X)_ хд2(гТ(г, X))

(г > К) (126)

дх г дг с соответствующими граничным и начальным условиями

Т (г, Х)| г=К = Тк (X), Т (г, Х)| х=0 = Т0. (127)

Для потока тепла чт (х) через поверхность сферической оболочки радиуса К справедливо общее соотношение

дТ (г, X)

Чт (х ) = -к

(128)

дг

С другой стороны, тот же тепловой поток чт (X) при учете его флуктуаций через поверхность оболочки радиуса К может быть определен с помощью выражения

ЧТ (X) = - ¿Ш + £ (X), (129)

где ^ (X) - случайный тепловой поток (источник Ланжевена), свойства которого зависят от характера источника флуктуаций, причем ^ (X= 0. Отметим, что соотношение (129)

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

Во многих случаях можно считать, что поток ^ (X) представляет собой белый шум с

интенсивностью V. Граничную частоту о такого шума можно оценить с помощью

2р

характеризующих задачу параметров: радиуса частицы К и коэффициента температуропроводности материала частицы %м , согласно выражению

огр =Ку. (130)

В частности, для медной частицы с коэффициентом температуропроводности Хм = 10-4 м2/с и радиусом 10 мкм для величины Огр получим значение порядка 106 с-1.

Величину интенсивности V случайного теплового потока можно оценить по формуле

к 1

V = -^квТ2, (131)

где кв - постоянная Больцмана. Для рассматриваемой выше частицы, помещенной в воду,

3 2 4

получим оценку V = 10- Дж /м с.

Решение поставленной задачи (126) - (127) будем искать аналогично задаче (79) - (81) с помощью введения вспомогательной функции / (г, 1), определяемой выражением

/ (г, 1) = гТ (г, 1). (132)

Решение такого рода задач представляет собой, как известно [4], сумму двух слагаемых. Первое слагаемое связано с наличием первоначального распределения температуры (и не зависит от соответствующих изменений температуры на границах), второе - обязано присутствием граничных функций температур или потоков. Из физических соображений ясно, что в рассматриваемой задаче, из-за наличия термодинамического равновесия в начальный момент времени, неслучайная составляющая потока тепла через границу среды, вызванная начальным распределением температуры, будет иметь нулевое значение. Тогда решение задачи (126) - (127) запишем с помощью выражения

/ (г,1) =

Я г г - Я

—^ ехР

(г - Я)2

х -т)

ТЯ (т)йт.

2л/ПХ ¡С -Т)3

Для получения уравнения, определяющего величину теплового потока воспользуемся соотношением, следующим из (128) и (129),

д/ (г,1)

дг

Я

=--Чт С) + Тя (X).

к

(133) Чт(1 ^

(134)

Определив производную 1 (г1) из формулы (133), получим

дг

д/ (г, 1)

дг

Я

24пх\

1 _ 2(г - Я)2 (1 -т)3/2 " 4х(1 -Т)5/2

" (г - Я)2 "

ехр _ 4х(1 -т) _

ТЯ (т)йт

(135)

Интегрирование по частям позволяет записать последнее выражение в виде

д/ (г,1)

дг

Я

урХ о V

'г 1

1^= ехр -т

(г - Я)2

4х(1 -т)

йТЯ (т)

йт

йт.

(136)

Нахождение из формулы (136) величины

д/ (г, 1)

дг

и ее подстановка в (134)

приводит к следующему соотношению для теплового потока Чт (1)

Чт(1) =

к

■Пс 0У^^-Т Я

йТЯ (Т) Т.

йт

(137)

Введя замены

1 (1) =

йТЯ (1)

л

£ =—£

(138)

г=Я

и воспользовавшись выражением (129), получим искомое соотношение, связывающее случайный тепловой поток ^ (х) через границу сферической оболочки и производную

температуры тк (X) на поверхности по времени, имеющее вид интегрального уравнения Вольтерра второго рода

7 С 1 I-^

1_^

1 (X) + 3к

7 (г) dт = iqт (X ). 039)

о ^ х-т к

Полученное уравнение (139), ядро которого представляет собой сумму слагаемого абелевого типа и постоянной величины, не может быть сведено к конечной системе стохастических дифференциальных уравнений. Таким образом, случайные процессы 1 (X) и qт (X), а также флуктуации температуры Тк (X), представляющей собой интеграл по времени от функции 1 (X), являются немарковскими [2].

Заметим также, что соотношение (139) справедливо при описании одномерного процесса теплопроводности в полупространстве над поверхностью плоского слоя конечной толщины. Действительно, рассмотрим сферическую оболочку радиуса К и толщины И, такой, что И << К. Теплоемкость С такой оболочки определяется равенством С = су ■ 4жR2 И . Вводя С

величину с =-= с И, представляющую собой теплоемкость оболочки, отнесенную к

единице площади ее поверхности, для теплового потока получим равенство

а (х) = - '4П2И ^^^ + г (X) = -с5 ^^^ + г (X) . Тогда для интегрального т 4жЯ2 dX ат 5 dX ат

стохастического уравнения, описывающего процесс распространения тепла в среде,

ограниченной бесконечным плоским слоем, вместо соотношения (139) получим выражение

1 (X) + -К-1 -Л-1 (т)dт = # (X) . (140)

сн п 0 V1 -т

Общий вид решения (139) можно представить с помощью интегрального оператора

1 (X) = 4 (X) + |к(X -т)4 (т^т, (141)

о

где К(X -т) - резольвента для уравнения (139), которая может быть записана с помощью бесконечного ряда [14]

ад

К(X -т) = 2 ип (X -т), (142)

п=1

в котором

(V - т) = и1 (V - я)ип (я - т) ds

- т) = -

3к

1 + 4ПХ

Л

^л/Пх ^л/Т-т я

(143)

(144)

У

Рис. 13. Схематический график функции К (т - т)

Расчет по формуле (143) показывает, что для реальных сред в случае не слишком больших значений V -т ряд (142) с хорошей степенью точности может быть определен с помощью суммы первых его членов. На рис. 13 изображен схематический график функции К(V -т), в случае, если сохранены первые 10 слагаемых (для простоты все постоянные приняты за единицу). Хорошо видно, что построенный график представляет собой убывающую кривую. Аппроксимируя его степенной функцией, получим в хорошем приближении зависимость

К (V -т) = 0,61(Т -т)-057. (145)

В частном случае теплопроводности над плоскостью (при Я = го ) ряд (142) принимает

вид

1 го

К (V-т) = — 2 (-1)4 (V-т) V-т п=1

п/2

где

а =-

1

' 3к ^

Г(п /2)

V

я4х

(146)

(147)

В последнем выражении Г( х) - гамма-функция.

Для одномерной V) и многомерной gL(А1,...,А1;^) характеристических

функций случайного процесса Z (V) получим соотношения

и

gl(Л;t) = ехр

gL ; tl,•••, tL) = ехР

9ул

2 í

2с1 R2

^у^ о

| К 2(> -т)ёт

9у

"2 п 2 ^

сул г ,¿=1

£ ЛЛ {К(^ -т)К-т)^т

(148)

(149)

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

Найденные формулы (148) и (149) позволяют определить моменты любого порядка для процесса 2^). Для математического ожидания {2)), момента второго порядка

7(^)2^2)) и дисперсии а2^) = ^2, получим

(2 «)> =

д^(Л; t)

гдЛ

= 0,

Л=0

2 (tl) 2&)) =

д 2 g 2(Л,Л2; t1, t2)

гдЛI дЛ2

9у }

Л=о, су П2

Л2=0

а2^) = ( 2 2(t)) = ( 2 ^) 2 (t)) = | К- т)От

| К(^ -т)К(^ -т>Г

с2 п

(150)

(151)

(152)

Уравнение (148) позволяет найти также одномерную плотность вероятности р(2) флуктуаций величины 2 ^) с помощью определения

1 ад

р( 2)=П t )е-Л2^Л •

Получим

Р( 2) = -

1

^ехр

,2 Л

2а ^)

(153)

(154)

что соответствует нормальному распределению Гаусса. Из выражения (152) видно, что дисперсия величины 2(^ растет с течением времени, из чего следует «размывание» плотности вероятности р(2) при увеличении t•

Найдем теперь спектральные плотности флуктуаций величины 2 ^), температуры поверхности Тп ^) и теплового потока qT ^). Для этого проведем преобразование Лапласа исходного интегрального уравнения (139). В результате имеем

Р) +

3к

1 4%

г( Р) = 4 (Р)

(155)

А/Х1ТР ЯР

где символами со «шляпкой» обозначены образы соответствующих функций, р - фурье-образ переменной 7.

С помощью определения спектральной плотности установившегося случайного процесса (при 7 = го ), а также при учете того факта, что спектральная плотность белого шума равна его интенсивности, для спектральной плотности флуктуаций Z (7) получим

02 (а) =

а

с2Я2а2 42кс¥Я 3/2 к2 л/2к 1/2 к2

—-+-а +—а +—¡=а +—г

9 з4х х Я

V

Из выражения (156) находим спектральные плотности флуктуаций ТЯ (7) и дТ (7):

(156)

(а) =

с!Я2а2 л/2кс„К 3/2 к 2 >/2

V'

' а3/2 +--а + -

к

I- 1С I 1С I I—

3л1х х

а12 +-

к

к

(157)

а* (а) =

с2 К 2а2 л/2кс„К

+ ■ а3/2

з4Х

к2 л/2к 1/2 к2

—а + —/=а + —г

V х Я2 у

V

(158)

9

6 5 4 3

2 [

О

ч/

к.

п

со, с

100

200

300

400

500

Рис. 14. Графики спектральной плотности флуктуаций температуры поверхности медной частицы От (а) при Я = 10мкм (кривая 1), Я = 12,5мкм (2), Я = 15мкм (3)

Графики на рис. 14 и 15 иллюстрируют найденные зависимости (157) и (158) для различных значений радиуса Я. При построении графиков параметр V рассчитывался согласно выражению (131). Из рис. 14 хорошо видно, что с увеличением радиуса Я графики спектральной плотности ОТ& (ш) флуктуаций температуры ТЯ ^) располагаются ниже

аналогичных графиков, соответствующих меньшим по размеру частицам. Последнее означает уменьшение дисперсии флуктуаций температуры сферической частицы с ростом ее радиуса. Спектральная плотность ОТ. (ш) в области низких частот оказывается обратно пропорциональной радиусу.

Рис. 15. Графики спектральной плотности флуктуаций теплового потока через поверхность медной частицы О (ш) при Я = 10мкм (1), Я = 12,5мкм (2), Я = 15мкм (3)

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

9. Теплопроводность вокруг микронитей

Рассмотрим бесконечно протяженное цилиндрическое тело малого радиуса Я, изготовленное из хорошо проводящего материала с объемной теплоемкостью Су и температуропроводностью хс, находящееся в неограниченной теплопроводящей среде с теплопроводностью К и температуропроводностью % (рис. 16). Будем считать, что

температура поверхности цилиндра в данный момент времени всюду одинакова и является некоторой функцией времени: Т (Я, 7) = ТЯ (7). Заметим, что в случае реальных сред такое

условие соответствует цилиндрическим телам не очень большого радиуса Я при рассмотрении явления теплопроводности на расстояниях, меньших или порядка Я от поверхности цилиндра. Будем также считать, что в начальный момент времени температура во всем пространстве, включая температуру материала цилиндра, была одинакова и равна нулю [19].

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

теплопроводности вида

dT (r, t ) _ / 3 T (r, t ) +1 dT (r, t

dt

= X

dr r dr

(159)

J

Рис. 16. Теплопроводность в пространстве, окружающем цилиндрическую поверхность Начальные и граничные условия задаются соответственно выражениями

Т (г, 0| 7=0 = 0, (160)

Т(г, 7)| г=я = Тя (7). (161)

Тепловой поток через поверхность цилиндрического тела дТ (7) определяется, как и ранее, общим соотношением

а^ (162)

(t ) = -3T (r,t )

It

r=R

dr

где функция (t) представляет собой случайный тепловой поток, обусловленный указанными ранее причинами, статистические свойства которого определяются характером источника флуктуаций, причем среднее значение (t= 0. С другой стороны, такой поток

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

Чт (^ = - С^ё^ (). (163)

2 dt

Следует подчеркнуть, что выражение (163) справедливо для случая предполагаемого нами наличия высокой тепловой проводимости материала цилиндра по сравнению с проводимостью среды, то есть при условии х << Хс.

Функция ^ (/), как и ранее, во многих случаях может быть представлена в виде

белого шума с интенсивностью V, ограниченного сверху некоторой предельной частотой ш , которую с точностью до постоянной можно оценить с помощью формулы

ш =%. (164)

тах я 2 4 '

Для оценки величины штах будем считать цилиндр изготовленным из меди (для которой коэффициент температуропроводности хс = 10-4 м2/с) и имеющим радиус Я = 10 мкм. Получим, что штах = 106 с-1. Таким образом, белый шум, характерный для флуктуаций теплового потока через цилиндрическую поверхность, оказывается широкополосным по спектру.

Оценку величины интенсивности V можно произвести по формуле

К 1

V = Я3 квТ2. (165)

При температуре Т0 = 300 К для медного цилиндрического тела радиусом Я = 50 нм, помещенного в воду, оценка, даваемая выражением (165), приводит к величине интенсивности v = 104 Дж2/(м4-с).

Можно показать, что если в некоторый момент времени т на каждой единице длины рассматриваемой цилиндрической поверхности мгновенно выделилось тепло Q(т), то температура среды при г > Я в момент времени t определится выражением [11]

Т(г, t) = ^(тО(г, Я, t -т), (166)

К

где О (г, Я, t -т) - функция влияния мгновенного цилиндрического источника тепла (функция Грина), определяемая согласно

С г2 + Я2 ^ Г Яг ^ 4х(-т)) 01Ш-т),

где /0( х) - модифицированная функция Бесселя нулевого порядка.

О(г, Я, t -т) = -—1-- exp

4пx(t - т)

Яг (167)

Отметим, что Q(т)S(t -т) = 2пЯдт (7) . Тогда в случае произвольного теплового потока через поверхность цилиндра дт (7), температура на расстоянии г от оси симметрии цилиндра будет определяться с помощью интеграла

т (г, 7) = Ж-

Г г2 + Я2 ^ /

л

4х(7 -г)

Яг

2х(7 -*).

йт •

(168)

Найдем производную

дТ (г, 7) дг

воспользовавшись выражением (168) и тем фактом,

что = !,( г):

д2

дТ (г, 7)

дг

Я г дт(т)

4хК0(7 -г)2

ехр

4х(7 -т)

Я!,

Яг

2х(7 -г).

- г!,

Яг

2х(7 -т).

(169)

йт.

Из полученного соотношения (169), при учете формул (162) и (163), находим

^ (7 -т)2 (т)йт = ^Тт (7),

С,Я

2

(170)

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

2 (7) = ^^ (171)

Ж

и ядра интегрального уравнения (170)

¿(7 -т) +

К(7 -т) = Я2

4Х(7 -т)

ехр

Я

2 Л

2Х(7 -т)

!

Я2

2Х(7 -т)

} ( -!

У

Я2

2Х(7 -т)

(172)

Таким образом, процесс распространения тепла в пространстве, окружающем цилиндрическую поверхность малого радиуса Я, при наличии случайного теплового потока через нее, описывается стохастическим интегральным уравнением Вольтерра второго рода (170), что означает немарковский характер флуктуаций величины 2 (7), а, следовательно, и

флуктуаций температуры поверхности цилиндра тЯ (7) и потока дт (7).

Заметим, что при больших значениях радиуса Я выражение (172) может быть записано с помощью приближенной формулы [15]

К (7-т) = 3(7-т) + -

4х

Яу1 - т) '

при подстановке которой в уравнение (5.59) вместо получим интегральное уравнение вида

(173)

суя

2 1(')^ЖТ1 (т)Л=^(174)

описывающее явление теплопроводности в полупространстве над плоской поверхностью.

Отметим также на уже отмечающийся ранее факт, заключающийся в том, что для случая, когда основным источником теплопроводности является теплоперенос, то верхний предел интегрирования в уравнении (170) следует, вообще говоря, брать равным 1 — 81, где 81 - малый параметр, по порядку величины равный времени свободного движения частиц теплопроводящей среды.

Рассмотрим теплопроводность в среде, окружающей цилиндрическую оболочку, в Я2

случае, когда параметр —<< 1, что, как видно из выражения (164), соответствует случаю

X

очень большой граничной частоты спектра флуктуаций теплового потока через цилиндрическую поверхность (белый шум). Положив также 1 — т >81, после разложения модифицированных функций Бесселя и экспоненты в ряд и удержания первых слагаемых разложения, получим

п 2 1 1 2

1 (1) — — Г—Ц- 1 {т)с1т=— £ (1). (175)

4x0(1 -т)2К) СуЯ^'

Решение полученного интегрального уравнения Вольтерра второго рода (175) в

общем случае можно записать в виде [14]

1 (1) = 7^ 1 [8(1 — т) + Р (1 — т)]£т (г)Л-, (176)

СУЯ 0

где резольвента Р(1 — т) определяется реккурентным соотношением

Р (1 — т) = £ Рп (1 — т). (177)

п=1

Здесь

Р(1 — Г) = —(178) 1 4Х(1 — т)2

г—8г

Р (1 — т) = | Р1(1 — 8)Рп—х(8 — г)^, п > 1. (179)

т+8г

п2

Указанное выше условие _<< 1 делает ряд (177) быстро сходящимся. При условии

X

п2

_<< 81 (цилиндрические поверхности малого радиуса, или тонкие нити), вторым и

X

последующими слагаемыми ряда (177) можно пренебречь по сравнению с первым, и окончательно получим

Я2 1

Р (г-г) =

(180)

4х (г -г)2

Исходное интегральное уравнение (176) для случая резольвенты Р (г -г) вида (180) позволяет найти любые статистические характеристики процесса 2(г). Так, для одномерной g1(Л; г) и многомерных gL (Л[,...,ЛЬ; г1,..., гь ) характеристических функций процесса 2 (г) при условии, что случайный тепловой поток ^ (г) представляет собой белый шум интенсивности у, получим

gl(Л; г) = ехр

1 Я \Л

2,л2 Г

24 хС

V V

J__1

7

(181)

gL (Л,...А; г1,..., tL) =

= ехр

1 Я2у

Е

Л1Лк

24 х2су & .(г, - гк )2

к <1

_1+ч - г++

8г гчк (г, - ч) г, - ч гк (г, - ч)

(182)

Графики функций g1(Л; г), задаваемых соотношением (181), изображены на рис. 17 для различных значений г. Здесь и в дальнейшем в качестве примеров использованы медные цилиндрические микротела малого радиуса, помещенные в воду (Я = 10-8 м, у = 104 Дж2/(м4-с), х = 2 -10-7 м2/с, Су = 1,3 -103 Дж/(Юм3)).

Рис. 17. Графики функций gl(Л;г) для значений г = 10 7 с (кривая 1), г = 10-6 с (кривая 2), г = го (кривая 3)

Выражения для характеристических функций, определяемых уравнениями (181) и (182), позволяют найти любые статистические характеристики процесса 1 (1). Для математического ожидания {2(1)) и дисперсии Б1 (1) процесса 1 (1) из (181) получим

(1 (1)) =

д&(Л; 1)

гдЛ

= 0,

Л=0

В, (1) = 12(1) = —

д&(Л; 1)

дЛ2

я V

1__1

8 ~ 7

(183)

(184)

Л=0 24х Су V*

Из найденной формулы (184) видно, что установившийся процесс теплопроводности (при 1 ^ да ) сопровождается постоянным значением дисперсии В1 (1) скорости изменения

температуры цилиндрической поверхности 1 (1):

В (1 )| = К у 1 (1 )1 ^ = 24х2С2813

(185)

Для корреляционной функции К (11, 12) = ( 1 (11) 1 (12)) при помощи двумерной характеристической функции g2(Л1,Л2; 11, 12), определяемой из (182), получим

К (11, О = —

дЕ 2(Л1,Л2; ^ 12)

дЛ1дЛ2

Л=0,

Л=0

(186)

1 Я2у

1

48 х2С2у(12 —11 )2

1 1 — 1т +2 , 181

— + ^-2-^ +-1п- 2

81 1112(12 — 11) 12 — 11 11(12 — 11)_ Из зависимости (186) следует, что при 11 ^да, 12 ^да (установившийся процесс

теплопроводности), корреляционная функция К (11,12) = К (12 —11) = К (т), принимает вид

К (т) =

1 Я2у 1

48 хС2 г2

1 12, 81 — + — + — 1п— 81 т т т

(187)

Графики функций К (1 — т, 1), задаваемых соотношением (186), изображены на рис. 18 для

различных значений параметра 1.

Найденная корреляционная функция (186) позволяет найти одностороннюю спектральную плотность процесса 1 (1)

I

01 (а, 1) = 21К (1 — т, 1 )со$>атс1т.

(188)

На рис. 19 проиллюстрирован результат численного вычисления односторонней спектральной плотности 01 (а, 1) с помощью формулы (188) для различных значений

времени t. В частном случае малых значений с и не очень больших I из (186) и (188) при учете условия 8t << t получим

^ (с, t )| ^ = 2{ К (t-г, t)

1 -

„2_2Л со г

2

йг =

1 Я2У t

48 /с2 а

_4_

-с

(189)

Рис. 18. Графики функций К- г, t) для значений t = 10 7 с (кривая 1) t = 10-6 с (кривая 2), t = 1с (кривая 3)

<?г(й), 0

Ю"10, К2 /с

4,001 4

3,999 3,998 3,997 3,996 3,995 3,994 3,993

1

си, с

20

40

60

80

100

Рис. 19. Графики односторонней спектральной плотности 02 (с, t) для значений времени

t = 0,01с (кривая 1) и t = 1с (кривая 2)

Для одномерной функции плотности вероятности р(2; t) согласно определению найдем выражение

1 ад

р(2, t) = — [ g1(Л; t)ехр(-/Л2)йЛ = 2п

(190)

V

1

2яБ (t)

ехр

,2 Л

2 В(1 )

График на рис. 20 изображает зависимость (190) для различных 1. Хорошо видно, что с течением времени график плотности вероятности флуктуаций величины 1 (1) «размазывается» вдоль оси 1, стремясь к стационарной гауссовой кривой при 1 ^ да .

р&Л, Ю10, с/К

г,ю'9, к/с

-2,0 -1,5 -1,0 -0,5 0 0,5 1,0 1,5 2,0

Рис. 20. Зависимость плотности вероятности р(1, 1) для 1 = 10—7 с (кривая 1) и 1 = 10—6 с

(кривая 2)

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

10. Заключение

Проведенное описание движения микрочастицы в вязкой среде, а также явлений теплопроводности в средах с микронаноструктурой, как случайных процессов, носящих немарковский характер, показывает их значительное отличие от аналогичных процессов, исследуемых классическими методами. Результаты работы представляют большое значение

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

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

1. Пугачев В.С., Синицын И.Н. Стохастические дифференциальные системы. - М.: Наука, 1990. - 632 с.

2. Морозов А.Н. Необратимые процессы и броуновское движение. - М.:Изд-во МГТУ им. Н.Э.Баумана, 1997. - 332 с.

3. Бочков Г.Н., Кузовлев Ю.Е. Новое в исследованиях 1/Г-шума // Успехи физических наук. - 1983. - Т. 141., вып. 1. - С. 151 - 176.

4. Букингем М. Шумы в электронных приборах и системах: Пер. с англ. - М.: Мир, 1986. - 399 с.

5. Феллер В. Введение в теорию вероятностей и ее приложения: В 2 т.: Пер. с англ. -М.: Мир, 1984. - Т. 1. - 528 с.

6. Климонтович Ю.Л. Статистическая физика. - М.:Наука, 1982. - 608 с.

7. Гардинер К.В. Стохастические методы в естественных науках. - М.: Мир, 1986. -

528 с.

8. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. - М.: Наука, 1986. - 736 с.

9. Bagayev S. N., Orlov V. A., Panov S. V. The Law of Motion of Brownian Particles at Short Observation Times // Laser Physics. 2001. Vol. 11. No. 11. P. 1232-1234.

10. Lucic B., Jeney S., Tischer C. et al. Direct observation of nondiffusive motion of a Brownian particle // Phys. Rev. Lett. 2005. Vol. 95. P. 160601-1-160601-4.

11. Тихонов А.Н., Самарский А.А. Уравнения математической физики. - М.: Наука, 1977. - 736 с.

12. Прудников А.П., Брычков Ю.А., Маричев О.И. Интегралы и ряды. - М.:Наука, 1981. - 800 с.

13. Рубинштейн Л.И. Проблема Стефана. - Рига: Звайгзне, 1967. - 458 с.

14. Волтерра В. Теория функционалов, интегральных и интегро-дифференциальных уравнений. - М.: Наука, 1982. - 304 с.

15. Никифоров А.Ф., Уваров В.Б. Специальные функции математической физики. -М.:Наука, 1984. - 344 с.

16. Морозов А.Н. Необратимые процессы и броуновское движение. - М.: Изд-во МГТУ им. Н.Э.Баумана, 1997. - 332 с.

17. Morozov A.N., Skripkin A.V. Spherical particle Brownian motion in viscous medium as non-Markovian random process // Physics Letters A. 2011. V. 375. P. 4113 - 4115.

18. Морозов А.Н., Скрипкин А.В. Описание испарения сферической частицы жидкости как немарковского случайного процесса с использованием интегральных стохастических уравнений // Известия вузов. Физика. 2010. №11. С. 55 - 64.

19. Морозов А.Н., Скрипкин А.В. Распространение тепла в пространстве вокруг цилиндрической поверхности как немарковский случайный процесс // Инженерно-физический журнал. 2011. Т. 84. №6. С. 1121 - 1127.

20. Морозов А.Н. Метод описания немарковских процессов, задаваемых линейным интегральным преобразованием // Вестник МГТУ, Естественные науки. 2004. №3. - С. 47 -56.

scientific periodical of the bauman mstü

SCIENCE and EDUCATION

EL JV® FS 77 - 4821 1. №0421200025. ISSN 1994-0408 electronic scientific and technical journal

Kinetic phenomena in media with microstructure

as non-Markov processes

# 05, May 2012

DOI: 10.7463/0512.0366037

Skripkin A.V.

Russia, Bauman Moscow State Technical University

[email protected]

The article presents a review of investigation of kinetic processes such as diffusion and thermal conduction in media with microstructure, near to the particles of micron and nanometer size, as well as microfilaments. The flow features of kinetic processes on small spatial and temporal scales are described. It was shown that these processes belong to the class of non-Markov processes and require the use of integral transformations for their description.

Publications with keywords: diffusion, thermal conductivity, medium with microstructure, non-Markov process

Publications with words: diffusion, thermal conductivity, medium with microstructure, non-Markov process

References

1. Pugachev V.S., Sinitsyn I.N. Stokhasticheskie differentsial'nye sistemy [Stochastic differential systems]. Moscow, Nauka, 1990. 632 p.

2. Morozov A.N. Neobratimye protsessy i brounovskoe dvizhenie [Irreversible processes and Brownian motion]. Moscow, Bauman MSTU Publ., 1997. 332 p.

3. Bochkov G.N., Kuzovlev Iu.E. Novoe v issledovaniiakh 1/f-shuma [New in research of 1/f-noise]. Uspekhifizicheskikh nauk, 1983, vol. 141, no. 1, pp. 151 - 176.

4. Buckingham M.J. Noise in Electronic Devices and Systems. Ellis Horwood (John

Wiley), Chichester, 1983. 372 p. (Russ. ed.: Bukingem M. Shumy v elektronnykhpriborakh i sistemakh. Moscow, Mir, 1986. 399 p.).

5. Feller W. An Introduction to Probability Theory and Its Applications . Vol. 1. 3rd ed. John Wiley & Sons, Inc., New York-London-Sydney, 1968. 509 p. (Russ. ed.: Feller V. Vvedenie v teoriiu veroiatnostei i eeprilozheniia. V. 1. Moscow, Mir, 1984. 528 p.).

6. Klimontovich Iu.L. Statisticheskaiafizika [Statistical physics]. Moscow, Nauka, 1982. 608 p.

7. Gardiner C. W. Handbook of Stochastic Methods for Physics, Chemistry and the Natural Sciences. 2d ed. Springer-Verlag, 1985. 442 p. (Russ. ed.: Gardiner K.V. Stokhasticheskie metody v estestvennykh naukakh. Moscow, Mir, 1986. 528 p.).

8. Landau L.D., Lifshits E.M. Gidrodinamika [Hydrodynamics]. Moscow, Nauka, 1986. 736 p.

9. Bagayev S.N., Orlov V.A., Panov S.V. The Law of Motion of Brownian Particles at Short Observation Times. Laser Physics, 2001, vol. 11, no. 11, pp. 1232-1234.

10. Lucic B., Jeney S., Tischer C., Kulik A.J., Forro L., Florin E.-L. Direct observation of nondiffusive motion of a Brownian particle. Phys. Rev. Lett, 2005, vol. 95, pp. 160601-1-160601-4.

11. Tikhonov A.N., Samarskii A.A. Uravneniia matematicheskoi fiziki [The equations of mathematical physics]. Moscow, Nauka, 1977. 736 p.

12. Prudnikov A.P., Brychkov Iu.A., Marichev O.I. Integraly i riady [Integrals and series]. Moscow, Nauka, 1981. 800 p.

13. Rubinshtein L.I. Problema Stefana [Stefan problem]. Riga, Zvaigzne, 1967. 458 p.

14. Volterra V. Theory offunctionals and of integral and integro-differential equations. Blackie & Son, 1930. (Russ. ed.: Volterra V. Teoriia funktsionalov, integral'nykh i integro-differentsial'nykh uravnenii. Moscow, Nauka, 1982. 304 p.).

15. Nikiforov A.F., Uvarov V.B. Spetsial'nye funktsii matematicheskoi fiziki [Special functions of mathematical physics]. Moscow, Nauka, 1984. 344 p.

16. Morozov A.N. Neobratimye protsessy i brounovskoe dvizhenie [Irreversible processes and Brownian motion]. Moscow, Bauman MSTU Publ., 1997. 332 p.

17. Morozov A.N., Skripkin A.V. Spherical particle Brownian motion in viscous medium as non-Markovian random process. Physics Letters A, 2011, vol. 375, pp. 4113 - 4115.

18. Morozov A.N., Skripkin A.V. Opisanie ispareniia sfericheskoi chastitsy zhidko-sti kak nemarkovskogo sluchainogo protsessa s ispol'zovaniem integral'nykh stokhasticheskikh uravnenii [Description of evaporation of a spherical particle of the fluid as a non-Markovian stochastic process with the use of integral stochastic equations]. Izvestiia VUZov. Ser. Fizika, 2010, no. 11, pp. 55 - 64.

19. Morozov A.N., Skripkin A.V. Rasprostranenie tepla v prostranstve vokrug tsi-lindricheskoi poverkhnosti kak nemarkovskii sluchainyi protsess [Distribution of heat in the space around the cylindrical surface of a non-Markovian random process]. Inzhenerno-fizicheskii zhurnal, 2011, vol. 84, no. 6, pp. 1121 - 1127.

20. Morozov A.N. Metod opisaniia nemarkovskikh protsessov, zadavaemykh lineinym integral'nym preobrazovaniem [Method of description of non-Markovian processes defined by a linear integral transformation]. VestnikMGTUim. N.E. Baumana. Ser. Estestvennye nauki [Bulletin of the Bauman MSTU. Ser. Natural science], 2004, no. 3, pp. 47 - 56.