
папа Жужа / Конспект лекций
.pdf
Одна из форм записи функции распределения Максвелла имеет вид
|
æ |
m |
ö |
3 2 |
− |
m υ2 |
|
|
|
|
|
|
2 k T |
|
2 |
|
|
||||
F(υ) |
ç |
|
÷ |
e |
|
|
4 p υ |
|
, |
(19) |
|
|
|
|
|||||||
F(υ) = ç |
|
÷ |
|
|
|
|||||
|
è |
2 pk T ø |
|
|
|
|
|
|
|
где υ – скорость на длине свободного пробега; m – масса одной молекулы; k – постоянная
υ Больцмана; Т – температура.
График функции F(υ) показан на рис. 10. Как и следовало ожидать, F(υ) = 0 при υ = 0 и υ = ¥, т.е. в газе нет неподвижных молекул и движущихся с бес-
конечно большими скоростями.
Найдем наиболее вероятную скорость υВЕР, определяющую максимум кривой распределения. Для этого следует взять произ- водную и приравнять её к нулю (постоянные множители при этом вынесутся за знак производной):
|
|
æ |
|
m |
|
|
|
ö3 2 |
|
|
d |
æ |
− |
m υ2 |
|
ö |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
ç |
|
|
|
|
÷ |
4p |
|
çe |
2 k T υ2 |
÷ = 0 |
, |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
ç |
2pk T |
÷ |
|
|
|
d υ ç |
|
|
|
|
÷ |
|
|
||||||||||||
|
|
è |
ø |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
æ m |
|
|
|
|
|
é − 2 k T æ |
|
|
m |
è |
|
ö 2 |
ø |
|
|
|
|
||||||||||
ö |
|
|
|
|
|
|
|
|
|
− 2 k T |
|
ù |
|
||||||||||||||
|
|
|
|
3 2 |
|
|
|
|
|
m υ2 |
|
|
|
|
|
|
|
|
|
|
|
|
m υ2 |
|
|
|
|
ç |
|
|
÷ |
|
|
4p êe |
ç |
- |
|
|
|
|
|
2 υ÷ |
υ + e |
|
|
2 υú = 0 |
, |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
ç |
|
÷ |
|
|
|
ê |
|
|
ç |
|
2k T |
|
|
÷ |
|
|
|
|
|
ú |
|||||||
è |
2pk T ø |
|
|
|
|
|
è |
|
|
|
ø |
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
ë |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
û |
|
|
|
|
æ |
|
|
m |
ö3 2 |
|
|
|
− |
m υ2 |
æ |
|
|
m υ |
2 ö |
|
|
|
|
|||||||
|
|
|
|
4 p υ e |
2 k T |
|
|
|
|
|
|
||||||||||||||||
|
|
ç |
|
|
|
|
÷ |
|
|
|
|
|
ç |
2 - |
|
|
÷ |
= 0 |
. |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
ç |
|
|
|
÷ |
|
|
|
|
|
|
|
|
ç |
|
|
k T |
÷ |
|
|
|
|||||
|
|
è |
2 pk T ø |
|
|
|
|
|
|
|
|
è |
|
|
ø |
|
|
|
|
Корнями последнего уравнения будут
υ = 0 |
, |
υ |
2 |
= ∞ |
, |
υ3 = |
2k T |
. |
1 |
|
|
|
m |
||||
|
|
|
|
|
|
|
|
Первые два корня – это минимумы функции F(υ), а третий корень – максимум (наиболее вероятная скорость):
|
|
|
|
|
|
|
υBEP = |
2kT |
= |
2RT |
. |
(20) |
|
|
m |
|
M |
20

Найдем значение функции распределения в максимуме, под-
ставив формулу (20) в уравнение (19):
F(υ |
|
|
m |
|
m |
|
æ |
|
m 2k T |
ö |
|
2k T |
|
|||||||||
BEP |
) = |
|
expç |
- |
÷ |
4p |
, |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
2pk T 2pk T |
ç |
|
|
|
|
÷ |
|
m |
||||||||||||
|
|
|
è |
|
2k T m ø |
|
|
|||||||||||||||
|
|
|
|
4 |
|
|
|
|
|
|
|
|
. |
|
|
|
|
|||||
|
|
|
F(υBEP ) = |
|
|
|
m |
|
~ |
m |
|
|
|
(21) |
||||||||
|
|
|
e 2 πk T |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
Из формул (20) и (21) следует, что при увеличении темпера- туры или уменьшении массы молекулы максимум кривой смеща- ется вправо и становится ниже (рис. 11). Однако площадь под кривой из условия нормировки (S = 1) сохраняется.
Зная функцию распределения F(υ), можно найти среднюю
(арифметическую) скорость:
F(υ) T1 < T2 m1 > m2
0 υ
Рис. 11
Тогда
|
|
|
|
|
|
|
|
|
|
|
υ = |
8k T |
= |
|
8 R T |
. |
(22) |
||||
π m |
|
π M |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
Так, например, для T = 300 К средние ско- рости молекул кислорода и водорода равны со- ответственно 445 м/с и 1782 м/с.
Функции F(υ) можно придать другой вид, удобный при расчетах, если ввести относи-
тельную скорость: u = υ/υBEP.
υ = u υBEP = u |
|
2 k T |
|
, |
d υ = |
2k T |
d u . |
|
m |
m |
|||||||
|
|
|
|
|
|
Если подставить получившиеся выражения для υ и dυ в формулу dN/N = F(υ) d υ, то можно получить функцию распреде- ления F(u) для относительной скорости:
d N |
|
m |
|
|
m |
|
æ |
|
mu2 2k T ö |
4pu2 |
2k T |
|
|
2k T |
|
d u |
|
|||||||
= |
|
|
|
expç |
- |
|
|
÷ |
|
|
|
, |
||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
N 2pk T 2pk T |
ç |
|
|
÷ |
|
|
|
|
|
m |
m |
|||||||||||||
è |
|
2k T m ø |
|
|
|
|
|
|
||||||||||||||||
|
|
d N |
= F(u) d u , |
где |
F(u) = |
4 |
|
|
e−u2 u2 . |
|
|
|
(23) |
|||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
N |
|
π |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
У 70% всех молекул скорость отличается от наиболее веро- ятной не более чем на 50% (рис. 12). А скорости, превышающие наиболее вероятную более чем в 5 раз, наблюдаются у одной из 12 млрд молекул.
21

Распределение Максвелла по- |
F(u) |
70% |
|
||||
зволяет объяснить существование и |
|
||||||
|
|
|
|
|
|
||
рассеяние атмосферы планет. Что- |
|
|
|
|
|
|
|
бы покинуть Землю, молекула долж- |
|
|
|
|
|
|
|
на иметь скорость, превышающую |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
0 0,5 |
1 1,5 5 u |
||||||
вторую космическую (11,2 км/с). Эта |
|||||||
скорость в 25 раз превышает наибо- |
|
|
Рис. 12 |
||||
лее вероятную скорость для молекул |
|
|
|
|
|
|
кислорода. Поэтому число покинувших Землю молекул кислоро- да очень мало. Однако легкие газы (водород, гелий) в основном рассеялись и остались «тяжелые» газы с небольшой скоростью молекул (азот, кислород, аргон, углекислый газ). Атмосферы со- хранились у тех планет, у которых сильное тяготение (высокая вторая космическая скорость) и низкая температура (низкая скорость самих молекул). Атмосферы состоят в основном из «тяжелых» газов – азот, кислород, аммиак, метан и т.п.
6.СРЕДНЯЯ ДЛИНА СВОБОДНОГО ПРОБЕГА
ИЧИСЛО СТОЛКНОВЕНИЙ
Минимальное расстояние, на которое сближаются центры двух молекул при столкновении, называется эффективным диа-
метром молекулы d, а величина s = pd 2 - эффективным сечением молекулы (рис. 13). Средняя длина свободного пробега мо-
лекулы – это прямолинейный участок пути, проходимый молеку- лой между двумя последовательными соударениями.
|
Пусть имеется газ с концентрацией n и движется только од- |
||
|
|
|
на молекула, а остальные неподвижны. За 1 с |
|
|
|
эта молекула пройдет путь, равный её сред- |
|
|
|
ней скорости áυñ, и столкнётся со всеми мо- |
d |
|
σ |
лекулами, центры которых расположены в |
|
|
|
ломаном цилиндре длиной áυñ и площадью |
|
|
|
|
|
|
|
основания, равной эффективному сечению |
|
d |
(рис. 14). Умножив объем цилиндра на кон- |
|
|
|
|
центрацию n, получим число столкновений: |
|
Рис. 13 |
n' = pd 2 áυñ n. |
22

В действительности движется не одна, а все молекулы, по-
этому в последнюю формулу должна входить не средняя ско- |
||||||||||
|
|
рость относительно стенок сосуда, а ско- |
||||||||
|
d |
рость относительно других молекул. Мож- |
||||||||
|
|
|
|
|
|
|
|
|
||
но доказать, что |
uOTH = |
2 u . |
Поэтому |
|||||||
|
|
|||||||||
|
|
|||||||||
|
|
среднее число столкновений за 1 с равно: |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ν = |
|
π d 2 υ |
|
n |
. |
|
|
Рис. 14 |
2 |
|
(24) |
Средняя длина свободного пробега l равна отношению длины пути, пройденного молекулой, к числу испытанных ею на этом пути столкновений:
λ = |
υ t |
= |
υ |
= |
|
1 |
= |
|
1 |
|
. |
(25) |
||
νt |
ν |
|
|
πd 2 n |
|
|
σn |
|
||||||
2 |
2 |
Таким образом, длина свободного пробега молекул тем меньше, чем больше их концентрация и эффективное сечение.
Оценим порядок величин длины свободного пробега и числа соударений (т.е. для формулы (25) необходимо найти d и n). Для нахождения d вспомним, что в жидкостях молекулы располага- ются достаточно плотно друг к другу. Один моль воды (18 г) за- нимает объём 18 см3 = 18 × 10 −6 м3. Разделив объём одного моля на число молекул в одном моле NA, получим приблизительно объём одной молекулы и её линейный размер:
|
|
V |
|
18×10−6 |
|
|
|
|
|
|
|
||
|
= |
= |
= |
|
−30 3 |
d = 3 30 ×10−30 » 3×10−10 |
|
||||||
VМОЛ |
M |
|
|
30 × 10 |
м. |
||||||||
6 ×1023 |
|||||||||||||
|
NA |
|
|
м , |
|
|
|
При испарении воды размер молекулы d не изменяется, но теперь один моль любого газа (и водяного пара) при нормальных условиях занимает объём 22,4 ×10−3 м3. Разделив число молекул в одном моле NA на объём одного моля газа, получим концентра-
цию n:
n = 6 ×1023 / 22,4 ×10 –3 » 3 × 1025 м −3.
Подставим найденные значения d и n в формулу (25): l = 1 / ( 2 × p × 9 ×10 –20 × 3 ×1025) » 10 −7 м.
23

Пусть средняя скорость молекулы водяного пара равна 600 м/с. Тогда из формулы (25) найдём число столкновений:
ν = λυ
= 10600−7 = 6 ×10 9 c−1.
Естественно возникает вопрос, можно ли считать идеальным газ, в котором молекулы каждую секунду сталкиваются миллиар- ды раз, т.е. «взаимодействуют». Однако длина свободного пробе- га в 100–1000 раз больше размеров молекулы. Другими словами, время столкновения молекул примерно в 100–1000 раз меньше времени между столкновениями. Следовательно, подавляющую часть времени молекулы движутся свободно.
Состояние газа в сосуде называется вакуумом, если длина
свободного пробега сравнима с линейными размерами сосуда или превышает их. Газ в этом случае называют ультраразреженным. Различные степени вакуума (низкий, средний, высокий, сверхвы- сокий) создаются в разных технических устройствах (электро- лампы, радиолампы, электронно-лучевые трубки телевизоров, системы вакуумного напыления металлов и т.п.).
ТРИЗ-задание 11. Вакуум в науке, технике и быту
Использование пустоты (вакуума) вместо вещества – это один из законов развития технических систем и один из изобретательских приё- мов. Наберите в любой поисковой системе в Интернете слова «вакуум- ный», «вакуумная», «вакуумное» и составьте список технических уст- ройств и технологий с использованием вакуума.
7. ЯВЛЕНИЯ ПЕРЕНОСА
Явления переноса объединяют группу процессов, связан- ных с неоднородностями плотности, температуры или скорости упорядоченного перемещения отдельных слоёв вещества. При нарушении равновесия система стремится к нему вернуться, при этом в веществе возникает упорядоченный, направленный пере-
нос массы (диффузия), импульса (внутреннее трение) и внут-
ренней энергии (теплопроводность). Такие явления представ- ляют собой необратимые процессы. Для рассмотрения явлений переноса необходимы два понятия из векторного анализа.
24

Градиент скалярной функции f (вектор) – это производная по направлению скорейшего (максимального) изменения этой
функции. Например, градиент вдоль оси z: df/dz. |
|
|
||||
Поток вектора (скаляр). Пусть име- |
|
|
||||
ется поверхность S (рис. 15), в каждой точке |
dS |
|
||||
которой можно установить нормали (пер- |
|
a |
||||
пендикуляры) dS . Пусть имеется векторное |
S |
dS |
||||
поле a , пронизывающее поверхность S. |
|
|
||||
Элементарный поток вектора dN равен ска- |
|
|
||||
лярному произведению: |
dN |
= (r |
). В |
ча- |
Рис. 15 |
|
|
a,dS |
|
|
стном случае, когда однородное поле пер- пендикулярно поверхности S, имеем N = aS.
7.1. ДИФФУЗИЯ В ГАЗАХ
Диффузия – обусловленное тепловым движением самопро- извольное взаимное проникновение и перемешивание частиц со- прикасающихся газов, жидкостей или твердых тел. Диффузия
происходит в направлении уменьшения концентрации вещества (т.е. оттуда, где вещества больше, туда, где его меньше).
Процесс диффузии может происходить между различными агрегатными состояниями вещества. Например, в полупроводни- ковой технологии p-n-переходы создаются диффузией примеси, находящейся в газообразном или жидком состоянии, в твёрдый полупроводник (нагретый до 1000–1300 °С).
Установлено, что поток молекул, прошедших за время dt че- рез площадку dS, перпендикулярную направлению переноса ве- щества, определяется законом Фика:
dN = −D |
dn |
dS dt |
, |
(26) |
|
||||
|
dz |
|
|
где dN – количество частиц (молекул); D – коэффициент диффу- зии; dn/dz – градиент концентрации в направлении оси z. Знак «минус» в формуле (26) обусловлен тем, что поток молекул на- правлен в сторону убывания концентрации. Знак «минус» ком-
пенсирует отрицательный знак градиента (рис. 16):
25

n1 > n2 |
|
|
|
|
dn |
Dn |
n |
- n (< 0) |
|
||||||
|
|
|
|
dz |
» Dz |
= z2 |
- z1 (> 0) < 0. |
|
|||||||
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||
z1 < z2 |
z |
|
|
|
|
|
2 |
1 |
|
|
|
|
|||
|
|
Умножив обе части формулы (26) на мас- |
|||||||||||||
|
Рис. 16 |
|
|
|
|
||||||||||
|
|
|
|
су одной молекулы, получим выражение для |
|||||||||||
|
|
|
|
|
|
потока массы и градиента плотности: |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
dm = −D |
d ρ |
dS dt |
. |
(27) |
||
|
|
|
|
|
|
|
|
|
|||||||
n1 |
dS |
|
|
n2 |
|
|
|
|
d z |
|
|
||||
|
|
|
|
|
N ' |
|
Рассмотрим процесс диффузии в |
||||||||
|
N '' |
|
|
|
|
газах и найдем коэффициент диффу- |
|||||||||
|
|
λ |
|
λ |
|
z |
зии. Пусть диффузия – стационарная, |
||||||||
|
|
|
|
т.е. градиент концентрации (или плот- |
|||||||||||
|
|
|
|
|
|||||||||||
|
n' |
|
|
n'' |
|
ности) является постоянным и не зави- |
|||||||||
|
|
|
|
||||||||||||
|
|
|
|
сящим от времени. Рассмотрим пло- |
|||||||||||
|
|
Рис. 17 |
|
|
|||||||||||
|
|
|
|
щадку dS, |
перпендикулярную |
оси z |
(рис. 17). Пусть n1 > n2.
Из-за теплового движения молекулы будут переходить через площадку dS и слева направо, и справа налево. Однако из-за раз- ности концентраций возникает диффузионный поток в направле- нии оси z. За время dt через площадку dS пройдут N молекул:
N = N ' − N ''. |
(28) |
Поскольку тепловое движение хаотическое, тепловые ско- рости молекул равномерно распределены по трем осям x, y, z. Причем из 1/3 всех молекул, которые движутся вдоль оси z, одна половина движется в положительном направлении оси z, а другая половина – в отрицательном. Таким образом, в положительном направлении оси z движется 1/6 часть молекул.
За время dt до площадки dS долетят все движущиеся по на- правлению к ней молекулы, заключенные в элементе объёма с основанием dS и высотой áυñdt. Умножив эти элементы объёма на концентрации молекул с одной и другой стороны площадки (n' и n''), получим количество молекул N ' и N '', прошедших через пло- щадку dS в противоположных направлениях:
N ' = |
1 n' υ dt dS , |
N '' = |
1 |
n'' υ dt dS . |
(29) |
|
6 |
||||||
|
6 |
|
|
|
26

Подставив уравнения (29) в формулу (28), получим:
N = |
1 |
(n'−n'') υ dS dt. |
(30) |
|
6 |
|
|
Очевидно, что через площадку dS будут пролетать лишь те молекулы, которые испытали последнее соударение на различ- ных расстояниях от dS. Значит, n' и n'' – это концентрации моле- кул на расстоянии длины свободного пробега λ по обе стороны от площадки. Так как dn/dz – разность концентраций, приходящихся на единицу длины, то на расстоянии 2λ разность концентраций
равна
n'−n''= −2λ dn . |
(31) |
|
|
dz |
|
Подставив уравнение (31) в формулу (30), получим: |
|
|
N = − 1 |
υ λ dn dS dt. |
(32) |
3 |
dz |
|
Сравнивая формулы (26) и (32), находим коэффициент диф- фузии:
D = |
1 |
υ λ |
. |
(33) |
|
3 |
|
|
|
Полученная формула (33) определяет коэффициент само- диффузии, т.е. диффузию молекул газа в среде того же газа.
Когда молекулы обеих компонент смеси неодинаковы по массе и эффективному сечению, коэффициент диффузии опреде- ляется более сложным выражением. Например, в смеси водорода и углекислого газа водород должен диффундировать быстрее, чем углекислый газ, так как средняя скорость молекул водорода почти в 5 раз выше, чем у углекислого газа, и кроме того, λ(H2) > λ(СО2). Это значит, что объём водорода, переносимого в одном направлении, больше объёма углекислого газа, переноси- мого в противоположном направлении. При этом неизбежно воз- никает разность давлений, которая влияет на процесс диффузии.
ТРИЗ-задание 12. Диффузия в науке, в технике и быту
Используя ключевое слово «диффузионный(-ая, -ое)», найдите в Интернете информацию для изучения технических устройств и техноло- гий с использованием диффузии. А что такое «диффузия инноваций»?
27

7.2. ВНУТРЕННЕЕ ТРЕНИЕ
Явлением внутреннего трения (вязкостью) называется возникновение сил трения между слоями жидкости или газа,
движущимися друг относительно друга параллельно и с разными
по величине скоростями. Например, при лами- |
|
нарном течении скорость жидкости или газа в |
|
трубе изменяется по параболическому закону: |
|
скорость максимальна в центре трубы и равна |
|
нулю у стенок (рис. 18). |
Рис. 18 |
Рассмотрим границу раздела двух сосед- |
|
них слоёв газа (жидкости), движущихся с разными скоростями (рис. 19). Более быстрый слой стремится увлечь за собой более медленный слой, действуя на него с силой F1, направленной по течению. Более медленный слой одновременно стремится замед- лить движение более быстрого слоя, действуя на него с силой F2.
Причиной вязкости является наложение двух движений: упорядоченного движения слоёв газа (жидкости) с различными скоростями и теплового движения молекул. При этом молекулы в своём тепловом движении, переходя из слоя в слой, имеют раз- ные скорости упорядоченного движения, и происходит перенос импульса упорядоченного движения молекул (рис. 20). Ось z на рис. 20 указывает направление потока импульса.
Медленный слой |
z |
|
|
|
F1 |
|
z |
||
|
|
|
|
υ1 |
|||||
|
|
|
|
|
|||||
υ1 |
|
|
|
|
|
|
|
dS |
υ2 |
υ2 |
|
|
|
|
|
|
|
υ3 |
|
Быстрый слой |
|
F2 |
|
|
dS |
|
υ 4 |
||
Рис. 19 |
|
|
|
|
|
|
|
Рис. 20 |
Модуль силы внутреннего трения dF, действующей на пло- щадку dS, лежащую на границе между слоями, определяется фор-
мулой Ньютона´ :
|
dυ |
|
|
|
|
||
dF = η |
|
|
dS |
, |
(34) |
||
dz |
|||||||
|
|
|
|
|
|
||
|
|
|
|
|
|
|
28

где h – коэффициент внутреннего трения (динамическая вязкость, или просто вязкость); dυ/dz = grad υ – градиент скорости (изменение скорости движения слоев на единицу длины в на- правлении нормали к поверхности слоя). Жидкости, подчиняю- щиеся уравнению (34), называются ньютоновскими´ (вода, гли- церин и др.), а неподчиняющиеся – неньютоновскими´ (в них вязкость зависит от градиента скорости; примером такой жидко- сти является раствор крахмала в воде). Единицей динамической вязкости является паскаль-секунда (Па×с). Динамическая вязкость газов при нормальных условиях имеет порядок 10−5 Па×с.
Наряду с динамической вязкостью используется также кинематическая вязкость n, определяемая как отношение дина- мической вязкости к плотности среды: n = h/r.
Учитывая второй закон Ньютона (dP/dt = F), формулу (34) можно представить в другом виде, как поток импульса dР через площадку dS:
dP = - h |
dυ |
dS dt |
. |
(35) |
|
dz |
|||||
|
|
|
|
Знак «минус» в формуле (35) обусловлен тем, что импульс передается в направлении убывания скорости (градиент скорости отрицателен).
Если известна плотность газа r, то из молекулярно-кине- тической теории можно найти, что коэффициент вязкости газа
определяется выражением
h = |
1 |
υ l r |
. |
(36) |
|
3 |
|
|
|
Расчет вязкости для жидкости очень сложен, так как дви- жение молекул в жидкости происходит за счет «перескоков» мо- лекул из одного «оседлого» положения в другое. Динамическая
вязкость достаточно хорошо описывается формулой вида
η ≈ AebT , где А и b – эмпирические постоянные, определяемые
свойствами жидкости, Т – температура. Динамическая вязкость жидкостей резко уменьшается с повышением температуры (а у газов увеличивается). Динамическая вязкость обычных, не очень вязких жидкостей (вода, бензин) имеет порядок 10−3 Па×с.
29