Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Жданов Явления переноса в газах и плазме 2008

.pdf
Скачиваний:
188
Добавлен:
16.08.2013
Размер:
6.05 Mб
Скачать

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

′ ′

(1.17)

f (v)f (v1 )dvdv1 = f (v)f (v1 )dv dv1 .

Можно показать, что произведения элементов объема в пространстве скоростей для прямых и обратных столкновений равны (см. гл. 4). Поэтому написанное выше условие переходит в соотношение

f (v)f (v1 )= f (v)f (v1) .

 

Логарифмирование этого соотношения дает

 

ln f (v)+ ln f (v1 )= ln f (v)+ ln f (v1) .

(1.18)

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

ln f = α0 m + α1

mv + α2

1

mv

2

 

 

1

βm(v u)

2

 

 

 

 

 

= α0

 

 

.

(1.19)

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

Скорость u просто характеризует выбор инерциальной системы отсчета и совпадает с макроскопической скоростью газа как цело-

го. В равновесном состоянии удобно полагать u = 0 .

 

В общем случае функция распределения принимает

вид

f = C exp[−β(v u)2 ] .

(1.20)

Константа C определяется из условия нормировки (1.16), которое можно переписать в виде

n = Cexp[−β(v u)2 ]dv = Cexp[−βv2 ]dv =

 

exp(−βvx2 )dvx

exp(−βvy2 )dvy

 

 

π 3 2

= C

 

 

 

 

 

 

exp(− βvz2 )dvz = C

.

 

−∞

−∞

−∞

 

β

 

 

 

 

 

 

 

 

(1.21)

21

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

exp(ax2 )dx = (πa)12 .

−∞

В результате

 

β

3 2

(1.22)

C = n

 

.

 

π

 

 

Для определения β вычислим среднюю кинетическую энергию

молекул, полагая для простоты u = 0 ,

 

 

 

 

 

 

1

mv2 fdv

 

 

 

 

 

 

 

 

 

 

 

 

EK =

2

 

 

 

=

mC

v2 exp(−βv 2 )dv =

 

 

fdv

 

2n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mC

3

 

m

 

 

 

 

= 2π

v4 exp(− βv2 )dv =

 

.

(1.23)

 

 

n

4

β

 

 

 

 

 

0

 

 

 

 

 

 

 

Здесь использовано представление (1.15) и известные значения интегралов из Приложения 1.

Сравнивая (1.23) и определение EK

(1.13), находим

 

β =

m

 

.

(1.24)

 

 

 

2kT

 

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

 

M

 

m 3

2

 

 

m

 

2

 

 

f

 

(v)= n

 

 

 

exp

 

v

 

.

(1.25)

 

 

 

2kT

 

 

 

 

2πkT

 

 

 

 

 

 

 

Легко показать, что если отказаться от условия u = 0 , то соответствующее распределение принимает форму

 

M

 

m

32

 

m

 

2

 

 

 

f

 

(v)= n

 

 

 

exp

 

(v u)

 

 

,

(1.26)

 

 

 

2kT

 

 

 

 

2πkT

 

 

 

 

 

 

 

что соответствует

введению

собственной

скорости

молекул

c = v u или переходу в систему отсчета, движущуюся с макроскопической скоростью газа как целого.

22

Используя (1.25) и (1.15), можно определить относительную долю молекул, абсолютные скорости которых лежат в некотором уз-

ком интервале значений между v

и v

+ dv ,

 

 

 

 

 

dn(v)

 

2

 

m

32

 

 

m

 

2

 

 

 

 

 

= 4πv

 

 

 

 

exp

 

 

v

 

dv

(1.27)

 

n

 

 

2kT

 

 

 

 

 

2πkT

 

 

 

 

 

 

 

Вид распределения dn/ ndv , описываемого выражением (1.27), для двух различных температур T1 и T2 > T1 представлен на рис. 1.3.

1 dn

n dv

T1

T2 > T1

vB (T1 ) vB (T2 )

v

Рис. 1.3

Максимум этого распределения лежит при значении скорости, которая носит название наиболее вероятной скорости молекул. Ее величина определяется условием

dn(v)

= 0 при

v = vВ

,

 

dv

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

 

 

1

mvВ2

= 0

 

 

 

 

2kT

 

 

 

 

 

 

 

 

 

или

23

vВ =

2kT

=

2RT

.

(1.28)

m

 

 

 

M

 

С ростом температуры максимум распределения смещается вправо. Определим еще среднее значение модуля скорости частиц газа или так называемую среднюю тепловую скорость молекул

 

dn(v)

 

m

3

2

 

 

 

 

m

 

 

 

 

 

 

v

3

 

 

2

 

 

 

 

4π

 

 

 

 

 

exp

 

v

 

dv

 

 

 

 

 

 

 

 

 

 

 

v = v n =

 

2πkT

 

 

 

 

 

 

2kT

 

 

 

.

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

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

v =

8kT

=

 

8RT

.

(1.29)

πm

 

Заметим, что

 

 

πM

 

 

 

 

 

 

 

v =1,128 vВ ,

 

v кв =1,225 vВ .

(1.30)

1.7. Частота столкновений молекул со стенкой

Получим еще одну важную для дальнейших приложений величину, а именно: определим среднее число соударений молекул с единичной поверхностью стенки в единицу времени (или частоту столкновений молекул с единицей поверхности стенки). Направим ось X перпендикулярно стенке. Выделим на стенке малую площадку s . Предположим, что молекулы, имеющие одинаковую

скоростьv ,

движутся в направлении стенки и построим на пло-

щадке

s ,

как на основании, косой цилиндр с образующей v t

(рис. 1.4.). За время t через поверхность стенки

s пройдет чис-

ло

молекул, находящихся в цилиндре, объем

которого равен

( vx

t )

s , поскольку ( vx t ) высота цилиндра.

 

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

ность потока частиц, скорости которых лежат между v

и v + dv ,

определяется выражением

 

vx+ f M (v)dv ,

(1.31)

24

t

 

 

v

S

x

 

 

 

v

Рис. 1.4

где vx+ составляющая скорости молекул в направлении положи-

тельных значений оси X .

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

 

 

 

 

 

 

 

 

 

m 3 2

 

 

 

 

 

 

 

 

 

mv2

 

 

 

ν =

v+ f M (v)dv = n

 

 

 

 

v

 

 

exp

 

x

dv

 

×

 

 

 

 

 

2kT

 

 

 

x

 

 

 

 

 

2πkT

 

 

x

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mvy2

 

 

 

 

mv

2

 

 

 

 

 

 

kT 1 2

 

×

exp

 

dv

y

exp

 

 

z

 

dv

 

 

= n

 

 

.

 

 

 

 

 

 

 

 

 

 

 

2kT

 

 

 

 

 

 

 

 

2kT

 

 

 

 

 

 

z

 

 

2πm

 

 

 

 

 

−∞

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая определение средней тепловой скорости

молекул v

(1.29), это выражение можно представить в виде

 

 

 

 

 

 

 

 

 

 

 

 

1

 

8kT 1 2

 

 

 

n v

 

 

 

.

 

 

 

 

 

(1.32)

 

 

 

 

 

 

ν =

 

 

n

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

πm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя (1.31), можно также строго обосновать вывод выражения (1.11) для давления газа p . Давление на стенку сосуда, в кото-

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

mvx+ (mvx+ )= 2mvx+ .

25

Плотность потока импульса молекул по направлению к стенке со скоростями между v и v + dv определяется умножением mvx+ на

величину (1.31)

mvx+vx+ f M (v)dv = m(vx+ )2 f M (v)dv .

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

p = 2m(vx+ )2 f M (v)dv =

 

m 3 2

 

 

mv2y

 

 

 

= 2mn

 

 

exp

 

dv

y

exp

 

2kT

 

2πkT

 

 

 

 

 

 

 

 

−∞

 

 

 

 

−∞

 

 

mv2

 

 

 

mv2

 

z

 

2

 

 

x

 

2kT

dvz vx

exp

2kT

dvx .

 

 

0

 

 

 

 

Производя интегрирование с использованием значений интегралов Приложения 1, получаем

p =

2

n

mv2

,

(1.33)

3

2

 

 

 

 

что совпадает с выражением (1.11).

1.8. Распределение Больцмана

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

n = n0 exp(mgz kT ) ,

(1.34)

где g - ускорение силы тяжести, z координата, отсчитываемая

вверх по направлению от поверхности Земли. Если температуру Т предполагать постоянной, такое же распределение с высотой имеет давление p . Распределение (1.34) является частным случаем об-

щего закона распределения Больцмана

26

n = n0 exp(U kT ) ,

(1.35)

где U потенциальная энергия молекулы во внешнем силовом поле, n0 значение плотности газа при U = 0 . Барометрической фор-

муле (1.34) соответствует при этом потенциальная энергия молекул в поле силы тяжести U = mgz . Если ввести плотность n (1.35) в

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

f MB = n0

 

m 3

2

 

 

m

v2

 

U

 

 

 

 

 

 

exp

 

 

.

(1.36)

 

 

2kT

 

 

 

2πkT

 

 

 

 

 

kT

 

Полученное выражение справедливо и для неоднородного силового поля U (r). При этом для числа молекул, попадающих в пределы объема dr = dxdydz , расположенного в точке r( x, y, z ) , и имеющих скорости в интервале dv = dvx dvy dvz , имеем

 

 

 

m 3

2

 

1

mv2

 

 

dn = n

 

 

 

 

 

exp

 

 

 

+U (r) dvdr .

(1.37)

 

 

 

 

2

 

0

 

2πkT

 

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Эта формула носит название распределения Максвелла-Больцмана.

1.9. Равновесное распределение по скоростям и внутренним состояниям молекул

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

27

Вращательное движение молекул иногда описывается с помощью некоторых классических моделей молекул твердых тел («шероховатые» и «нагруженные» сферы, сферо-цилиндры и т.д. [2]). Более общим и эффективным оказывается использование полуклассического метода описания, когда функция распределения fi = f (v,Ei ) предполагается зависящей от скорости поступатель-

ного движения молекулы v и от дискретных значений ее внутренней энергии Ei . Конкретный вид спектра этих значений энергии

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

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

равновесия приводит в этом случае к соотношению

 

fi (v)f j (v1)= fk (v)fA (v1) ,

(1.38)

где индексы i, j и k,A отвечают дискретным значениям энергии сталкивающихся молекул до и после столкновения соответственно. Из этого соотношения следует, что величина ln fi должна быть

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

 

 

 

1

 

2

 

 

ln fi = α0 m + α1

mv + α2

 

 

mv

 

+ Ei .

(1.39)

2

 

 

 

 

 

 

 

 

Из соотношения (1.38) следует, что коэффициенты α0 , α1 и α2

оказываются одинаковыми для любых i .

В отличие от случая одноатомных молекул здесь существенно то, что в столкновениях сохраняется полная энергия молекул, равная

сумме ее кинетической энергии mv2 2 и внутренней энергии Ei . Коэффициенты α0 , α1 и α2 находятся из условий, определяющих плотность n и среднюю энергию единицы объема газа nE . Если u = 0 , эти условия записываются как

28

n = fi dv , n E

= ∑∫

1

 

mv 2

 

 

 

 

 

 

 

+ Ei dv

(1.40)

 

 

i

 

 

 

 

 

i

 

 

2

 

 

 

 

 

 

 

и соответствующая функция распределения принимает вид

 

 

 

m 32

1

 

 

 

 

 

m

 

2

 

Ei

 

f (v, Ei

)= n

 

 

 

exp

 

 

 

 

 

 

v

 

 

 

(1.41)

 

Q

 

2kT

 

 

 

 

2πkT

 

 

 

 

 

 

 

kT

 

или

 

 

 

 

 

 

 

 

 

E

 

 

 

 

 

 

 

 

 

f (v, Ei )= f M (v)

1

 

 

 

 

i

 

 

 

 

 

 

 

 

exp

 

 

 

 

.

 

 

 

 

(1.41')

Q

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Величина

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q = exp

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.42)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

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

[3].

Функция распределения f (v,Ei ) записывается, таким образом, в

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

EJ ,v = EJ + Ev ,

где J соответствует вращательному, а v колебательному квантовому числу. При этом статистическая сумма Q записывается в виде произведения

Q = Q rot Qvib .

Статистические суммы для вращательных (rot) и колебательных (vib) степеней свободы определяются выражениями

Qrot = (2J +1)exp

 

EJ

 

 

 

 

E

v

 

 

,

Qvib = exp

 

 

. (1.43)

 

 

J

 

 

kT

 

v

 

kT

 

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

29

gJ = 2J +1 определяет число состояний с одной и той же энергией

иносит название кратности вырождения или статистического веса.

Определение средней энергии единицы объема газа nE (1.40)

подразумевает, что она представляет собой сумму вкладов от поступательных или транспортных (tr) и внутренних (int) степеней свободы молекул

 

 

 

 

 

n E =n Etr

+n Eint

,

 

 

(1.44)

где

 

 

 

 

 

 

Ei

exp(Ei / kT )

 

 

n E tr =

3

nkT =

3

p

 

 

 

 

, n Eint

= n

i

 

 

.

(1.45)

 

 

 

 

2

2

exp(Ei / kT )

 

 

 

 

 

 

 

i

Введем безразмерные переменные

 

Ei

 

 

E int

 

εi exp(− εi )

 

εi =

,

ε =

 

=

i

 

.

kT

kT

exp(− εi )

 

 

 

 

 

 

 

 

 

 

 

i

 

Тогда для средней энергии единицы объема газа имеем

3

 

n E = nkT

 

+ ε(T ) .

2

 

 

(1.46)

(1.47)

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

При одновременном учете вращательных и колебательных степеней свободы выражение (1.47) принимает вид

 

 

3

+ εrot (T )

+ εvib (T )

 

 

n E

= nkT

 

,

(1.48)

2

 

 

rot

 

vib

 

где усреднение в выражениях для ε производится независимо по вращательному и колебательному спектрам молекул.

30