Материал: Лебедев-Степанов ВВедение в самосборку ансамблей наночастиц 2012

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

nv = -DÑn - nbÑU ,

(7.47)

что является обобщением закона Фика на случай диффузии частиц во внешнем потенциале. Применение к (7.47) закона непрерывности (7.29) дает нестационарное уравнение диффузии во внешнем потенциале, обобщающее выражение (7.30):

∂n = Ñ(DÑn) + Ñ(nbÑU ) .

(7.48)

t

 

Аналогично, рассматривая дрейф газа при наличии градиента температуры (данное явление называется термодиффузией), с учетом (4.42) и (7.40), а также соотношения Эйнштейна, можно получить выражение для скорости дрейфа частиц:

v = -

D

ln e−0.5 nl3ÑT .

(7.49)

 

 

T

 

Коэффициент пропорциональности между скоростью дрейфа и величиной − ÑT по определению называется коэффициентом тер-

T

 

модиффузии DT , так что с учетом (7.49) имеем

 

D = D ln e−0.5 nl3 .

(7.50)

T

 

С учетом того, что при обычных температурах тепловая волна де Бройля λ очень мала по сравнению с расстоянием между частицами, коэффициент термодиффузии (7.50) отрицателен. Это означает, что молекулы невырожденного газа стремится концентрироваться в местах с более высокой температурой.

7.7. Аналогии с электричеством

Рассмотрим аналогичные величины (табл. 7.1).

Несмотря на очевидные параллели, соответствие между электричеством и диффузией не является точным. Так, концентрации отвечает и потенциал, и заряд. Аналогия становится более точной с обобщением теории диффузии: диффузионные явления связаны, во-

121

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

 

 

 

Таблица 7.1

 

Аналогии законов переноса

 

Параметр

Электричество

Диффузия

Обобщенная диф-

 

 

 

фузия

 

 

 

 

Заряд

Q

N (число частиц)

N (число частиц)

 

 

 

 

Плотность заряда

ρ

n (объемная кон-

n (объемная кон-

 

 

центрация)

центрация)

 

 

 

 

Плотность тока

J

Nv

Nv

 

 

 

 

Потенциал

ϕ

N

μ

 

 

 

 

Силовое поле

E = – Ñϕ

– Ñn

– Ñμ

 

 

 

 

Плотность мощ-

jE

– nv Ñn

– nv Ñμ

ности

 

 

 

 

 

 

 

Коэффициент

λ (проводимость)

D(коэффициент

b (обобщенная под-

переноса

 

диффузии)

вижность)

 

 

 

 

Закон проводи-

j = – λÑϕ (закон

Nv = – DÑn (закон

nv = – nbÑμ

мости

Ома)

Фика)

 

 

 

 

 

7.8. Уравнения движения в вязкой среде

Движение потока непрерывной среды описывается уравнением вида

ρ

dv

= -Ñp - f (ρ, v, (Ñv), Ñ2 v,...) ,

(7.51)

 

dt

122

где ρ – плотность среды, v – скорость потока, p – давление, член f описывает диссипацию энергии в гидродинамическом потоке. Согласно уравнению Навье-Стокса [25] указанная диссипативная функция для вязкой изотропной жидкости

f (ρ, v) = -hÑ

2

 

h

(7.52)

 

v - x +

Ñ(Ñv),

 

 

 

3

 

где η =νρ и ξ – коэффициенты динамической вязкости и второй

вязкости соответственно; ν – коэффициент кинематической вязкости. Для несжимаемой жидкости функция (7.52) приобретает вид

 

f (ρ, v) = -ηÑ2 v .

(7.53)

Уравнение Навье–Стокса имеет вид

 

ρ

dv

= -Ñp + ηÑ2 v ,

(7.54)

 

 

dt

 

условие несжимаемости

 

 

 

(Ñ × v) = 0 ,

(7.55)

где ρ и η – плотность и вязкость жидкости, v и p – скорость и дав-

ление как функции координат (радиуса-вектора) и времени.

Однако получить решение системы уравнений (7.54) – (7.55) для тела конечных размеров, движущегося в вязкой жидкости, в аналитическом виде без использования каких-либо приближений пока никому не удалось даже для самого простого случая, когда движущееся тело является твердым шаром.

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

123

При малых числах Рейнольдса Re = ulν (малая частица, движу-

щаяся с небольшой скоростью) из уравнения Навье–Стокса для несжимаемой жидкости можно получить:

В приближении Стокса подвижность частицы радиуса R равна

b =

1

.

(7.56)

 

6πRη

 

 

Коэффициент диффузии с учетом соотношения Эйнштейна

D =

kT

.

(7.57)

6πRη

 

 

 

Бассе (Basset) получил более точное аналитическое решение для вышеупомянутой задачи [25, с. 133]. Силу сопротивления, полученную с использованием этого решения, можно представить в виде суммы трех слагаемых.

 

1 du

 

3ν

 

3 ν

 

t

du dτ

 

 

F = 2πρR 3

 

 

 

+

 

 

u +

 

 

 

∫

 

 

 

 

 

. (7.58)

 

 

 

 

3

 

 

 

 

 

 

 

 

3 dt R

 

 

R π

−∞

 

 

 

 

 

 

 

 

 

 

dτ t − τ

 

В этой формуле три слагаемых. Первое (сила Стокса) пропорционально ускорению тела, второе – скорости движущегося тела. Наконец, третье слагаемое зависит не только от (векторного) значения скорости шара в данный момент, но и ее значений во все предшествующие моменты времени. Это третье слагаемое обычно называют силой Бассе, которая зависит от предыстории движения частицы

.

7.9. Броуновское движение наночастиц

Рассмотрим уединенную малую броуновскую частицу (наночастицу) плотности ρ и объема V в жидкости. Поведение такой части-

цы (для простоты будем рассматривать одномерное движение) обычно описывают уравнением Ланжевена вида

124

m

du

= −βu + F (t) ,

(7.59)

 

 

d t

B

 

 

 

 

где m = ρV – масса частицы, u – ее скорость, β = b−1 –

коэффици-

ент трения (при описании зависимости силы трения от скорости здесь используется приближение Стокса), FB (t) – флуктуирующая

сила Ланжевена (или ланжевеновский источник), с которой молекулы жидкости действуют на рассматриваемую частицу вследствие микроскопического теплового движения, причем стохастические свойства силы Ланжевена удовлетворяют соотношениям

FB (t) = 0 , FB (t ),FB (t') = 2βkTδ (t − t') .

(7.60)

При этом коррелятор тепловых флуктуаций скорости частицы

имеет вид

 

 

 

 

 

 

|t − t'|

 

u (t ),u (t') =

 

kT

 

 

 

−

 

 

 

 

exp

 

,

(7.61)

 

 

 

 

 

m

 

 

 

 

τ

 

 

где время релаксации скорости

 

 

 

 

 

 

 

 

 

 

τ =

m

.

 

 

 

(7.62)

 

 

 

 

 

 

 

 

 

β

 

 

 

 

Заметим, что разделение силы на две части – силу Стокса и силу Ланжевена – имеет смысл, если время корреляции ланжевеновского

источника τL много меньше характерных времен релаксации скоро-

сти. В нулевом приближении по параметру τL τ−1 время корреляции источника можно считать равным нулю. Такой случайный источник называется δ-коррелированным. Наличие дельта-функции δ (t − t') указывает, что время корреляции τL много меньше харак-

терного времени релаксации скорости τ .

Легко убедиться, что, согласно (7.61), средний квадрат флуктуаций скорости, как и требуется, отвечает теореме о равнораспределении поступательной кинетической энергии по степеням свободы:

u 2 =

kT

.

(7.63)

 

m

125

Источник: https://studfile.net/preview/16708784/