Если, например, тело представляет собой шар с радиусом R, центр которого в момент t находится в точке r0, последнее выражение примет вид
. S = ∫δ(r − r0 (t) − R 2 )dV . |
(7.6) |
Сила, действующая на тело, может быть найдена как интеграл тензора потока импульса по поверхности:
F i = ∫Πik dSk , |
(7.7) |
где интегрирование производится по поверхности тела. С учетом (7.3) это можно переписать в виде
F i = ∫ |
∂Πik |
(7.8) |
|
∂x |
k dV , |
||
|
|
|
|
где интегрирование ведется по объему данного тела и вблизи него (там, где имеет место взаимодействие). Для непроницаемого жесткого шарика интегрирование ведется по тонкому приповерхностному слою.
Рассмотрим равновесный идеальный газ с концентрацией n и температурой T. Плотность потока x-компоненты импульса через площадку, перпендикулярную оси y, определяется формулой mnvxvy, где n=δ(x,y,x). Речь идет о тепловом движении молекул. Вообще, тензор энергии-импульса запишется в виде:
Πik = ∑mδ j vij vkj . |
(7.9) |
j |
|
Усредняя по времени, достаточно большому, чтобы за него произошло достаточно много столкновений (при этом предполагаем справедливой эргодическую гипотезу), получаем
Πrik = ∑m < δjvijvkj >= mn < vivk >, |
(7.10) |
||
j |
|
|
|
где n – средняя плотность. |
|
|
|
С учетом того, что < vi vk >= 0 , если i ¹ k и< vi vk |
>= |
1 |
< v 2 > |
|
|||
3
при i = k , где
111
< v 2 >= |
3kT |
– |
(7.11) |
|
m |
||||
|
|
|
средний квадрат скорости, для тензора потока импульса (7.10) будем иметь
Π r ik = nkTδik |
= pδik , |
(7.12) |
||
где |
|
|
||
p = |
1 |
nm < v 2 |
>= nkT – |
(7.13) |
|
||||
3 |
|
|
|
|
давление идеального газа, δik |
– символ Крон екера. |
|
||
Таким образом, тензор потока импульса равновесного идеального
газа определяется его давлением. |
Подставляя (7.12) в (7.8), получа- |
|||||||
ем |
|
|
|
|
|
|
|
|
F i = ∫ |
∂p |
|
= k ∫ |
∂nT |
|
|||
|
|
dV |
|
i |
dV . |
(7.14) |
||
∂x |
i |
∂ |
||||||
|
|
|
|
|
x |
|
||
Впростейшем случае возмущение потока телом можно рассматривать как эффект непроницаемости поверхности тела для проникновения в него газа. Тогда граничные условия для импульса на поверхности сведутся просто к прерыванию потока на границе.
Вэтом случае легко показать, что сила, действующая на тело, будет пропорциональна градиенту концентрации и размерам тела.
Рассмотрим теперь идеальный газ, в котором мгновенно имеет место небольшой градиент концентрации вдоль оси x:
n = n0 + ax . |
(7.15) |
При этом предположим, что функция распределения молекул по скоростям остается максвелловской. В этом случае на тело, помещенное в газ, действует силапропорциональная градиенту плотности
F = −kT |
∫ ∂x |
∫ |
∂x |
∂n dV = −kTa |
|
dV = −kTaV = −kTV ∂n , (7.16) |
|
x |
|
|
|
где интегрирование проводится по объему взаимодействия V. Знак минус говорит, что речь идет о силе, приложенной средой к телу: эта
112
сила направлена против градиента давления. Разделив на объем взаимодействия, отвечающий частице, получаем плотность силы
f = −kT |
∂n . |
(7.17) |
x |
∂x |
|
|
|
Если имеется ансамбль таких частиц, то формула (7.17) характеризует среднюю силу, действующую на частицы, находящиеся в единице объема.
7.2. Подвижность
Итак, градиент концентрации приводит к наличию силы, действующей на помещенный в газ неподвижный предмет. Под действием этой силы тело, если оно будет отпущено, приобретет некоторую скорость, величина которой определяется подвижностью.
Как показывает эксперимент, при малых (по сравнению с тепловыми) скоростях движения частицы в среде, сила сопротивления среды F пропорциональна скорости и направлена в противополож-
ную сторону: |
|
v = −bF , |
(7.18) |
При рассмотрении подвижности b – |
коэффициента пропорцио- |
нальности между силой и скоростью – |
нужно выделять два случая: |
большая (по сравнению со средней длиной свободного пробега) броуновская частица и малая молекула (например, того же сорта, что и все остальные).
Если броуновская частица движется в направлении действия на нее силы, то возникает дополнительная сила сопротивления. Природа этой силы – в асимметрии числа столкновений частицы с молекулами среды.
Если бы броуновская частица с площадью сечения столкновения S покоилась, то за единицу времени с ним фронтально столкнется Z=vn частиц. Столько же столкновений было бы и с противоположной стороны. Если частица движется со скоростью много меньшей тепловой, то количество ударов спереди и сзади по ходу движения уже не будет в среднем совпадать.
113
Получим давление усреднением ударов о стенку.
Молекула газа, столкнувшись со стенкой и упруго от нее оттолкнувшись, передает ей импульс 2mvx , где vx – проекция скорости
частицы до столкновения на нормаль к поверхности. Какое число молекул сталкивается с поверхностью в единицу времени? Половина из общего числа молекул имеет проекции скорости, направленные в сторону поверхности. Концентрация молекул, скорость кото-
рых находится в интервале [ vx , vx + dvx ] определяется распределе-
нием Максвелла (предполагаем, что газ находится в состоянии термодинамического равновесия):
|
m |
1 |
|
|
2 |
|
|
|
|
|
|
|
|||||
2 |
|
|
||||||
|
|
|
mvx |
|
|
|||
dn = n |
|
|
|
exp |
− |
|
dvx . |
(7.19) |
|
2πkT |
|
|
|
2kT |
|
|
|
Поток этих молекул через единицу площади рассматриваемой частицы (количество ударов молекул с данной проекцией скорости о
единицу площади за единицу времени) равен vx dnx . Тогда импульс,
переданный молекулами, имеющими нормальную к площадке скорость в указанных пределах, единице площади за единицу времени определяется выражением
|
|
m |
1 |
|
|
|
2 |
|
|
2 |
|
|
|||||
2 |
|
2 |
|
|
mvx |
|||
dp = 2mvx |
dn = 2mn |
|
|
|
vx |
exp |
− |
|
|
|
|
||||||
|
|
2πkT |
|
|
|
|
2kT |
|
Интегрируя, находим
m
p = 2mn
2πkT
1 |
∞ |
|
|
2 |
|
|
∫0 |
|
|
vx |
|
2 |
|
|
2 |
|
|
|
− |
mvx |
|
= |
|
|
exp |
|
dvx |
||
|
2kT |
||||
|
|
|
|
|
|
dvx . (7.20)
nkT . (7.21)
Таким образом, получаем давление идеального газа, как и должно, быть в рассматриваемом случае.
Если площадка движется со скоростью ux вдоль положительной нормали к поверхности, то величина переданного импульса возрастает. Из-за роста относительно скорости увеличивается число столк-
114
новений в единицу времени для молекул по сравнению со случаем неподвижной стенки. Учитывая оба фактора, получаем
|
|
m |
1 |
∞ |
|
|
|
|
|
2 |
|
|
|
||||
|
2 |
|
|
2 |
|
|
|
||||||||||
|
|
∫(vx + u x ) |
|
|
mvx |
|
|
|
|||||||||
pu |
= 2mn |
|
|
|
|
|
exp |
− |
2kT |
dvx ≈ |
|||||||
|
|
2πkT |
|
0 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
1 |
|
2ux kT |
|
|
|
|
|
|
|||
|
|
m |
|
|
|
|
|
|
|
m |
|||||||
|
|
|
|
|
2 |
|
|
|
|
|
|||||||
≈ nkT + 2mn |
|
|
|
|
|
|
= nkT 1 + 4ux |
|
|
||||||||
|
|
|
|
m |
|
||||||||||||
|
|
2πkT |
|
|
|
|
|
|
2πkT |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1
(7.22)
2 .
Здесь предполагается, что скорость движения стенки много меньше тепловой скорости, поэтому членом, пропорциональным ux2, можно пренебречь. Если площадка движется против нормали, то ux будет отрицательной, как и соответствующий член в скобках в последнем выражении. Если тело имеет эффективное сечение площадью S, то согласно (7.1) на него будет действовать сила
|
|
|
1 |
|
|
|
|
|
1 |
|
|
||
|
|
|
m |
|
|
|
mkT |
|
|
|
|||
|
= 8u |
2 |
nkTS = 8u |
|
2 |
|
|
||||||
F |
|
|
|
|
nS |
|
|
. |
(7.23) |
||||
|
|
|
|||||||||||
x |
|
x |
2πkT |
|
x |
|
2π |
|
|
|
|
||
Здесь найдена сумма давлений на тело «по» и «против» хода движения. Это позволяет определить подвижность:
|
u x |
|
1 2π |
|
1 |
|
|
||
b = |
= |
2 |
|
|
|||||
|
|
|
|
|
. |
(7.24) |
|||
|
|
|
|||||||
|
Fx |
|
8nS mkT |
|
|
|
|||
Тогда |
|
|
|
|
|
|
|
|
|
|
|
Fx |
= bux . |
|
|
|
(7.23а) |
||
Здесь F – внешняя сила, обратная по знаку силе сопротивления. Это теоретически доказывает линейность силы сопротивления
среды к скорости движения тела в среде.
115