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