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

Физика / 4.___________ ______

.pdf
Скачиваний:
9
Добавлен:
05.02.2016
Размер:
303.68 Кб
Скачать

В.М.Клименко. Статистична фізика

104

 

 

Статистична фізика

§ 80. Імовірність та флуктуації

Наведемо деякі визначення поняття ймовірності.

1. Імовірність wі деякого і-го стану системи визначається границею відношення часу tі, на протязі якого система знаходиться в даному стані, до повного часу Т спостереження за системою при необмеженому збільшенні Т

wi = lim

1

ti .

(1)

 

T→∞ T

i

 

2. Якщо фізична величина А при N вимірюваннях мала Ni число значень Аі, то ймовірність wi того, що величина А має значення Аі дорівнює границі відношення Ni/N при необмеженому зростанні N і вона співпадає з

імовірністю wі і - го стану системи

 

wi = lim

Ni

.

(2)

 

N→∞ N

 

3. Імовірність dw(A) того, що фізична величина має значення в

інтервалі (A; A + dA) пропорційна dA

 

dw(A) = f(A)dA,

(3)

де f(A) називається густиною ймовірності або функцією розподілу ймовірності.

Умова нормування ймовірності для дискретних станів

Wi = 1

(4)

i

 

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

f (A)dA = 1

(5)

сума ймовірностей для всіх елементарних значень параметра стану А є ймовірність достовірної величини і дорівнює 1.

Середнє статистичне значення величини А позначається A і

визначається для дискретних станів так

 

A = Aiwi ,

(6)

i

а для неперервної зміни величини А як

A = Adw ,

де сума й інтеграл беруться по всім можливим станам системи.

Флуктуації це випадкові відхилення термодинамічних параметрів від рівноважних значень. Якщо термодинамічна система складається з невеликого числа частинок, то розподіл частинок, випадково відхиляючись

δρ ,δp ,δT ~

В.М.Клименко. Статистична фізика

105

 

 

від рівноважного, може суттєво змінювати величини термодинамічних параметрів. Для кількісної оцінки флуктуації деякої величини А

використовують квадратичну флуктуацію, яку ще називають дисперсією σ2A .

Дисперсія є середнє значення квадрата відхилення А від її середнього значення:

σ2A = (A)2 = (A − ∆A)2 = A2 2AA + A2 = A2 A2 0.

Величина σA = σ2A називається абсолютною флуктуацією, а δA = σAA

відносною.

Природа флуктуацій термодинамічних величин визначається хаотичним тепловим рухом структурних частинок середовища. Чим більше число N цих частинок, тим менший вплив їх флуктуацій на значення термодинамічних величин. Доведено, що для хімічно однорідного ідеального газу, при сталому об'ємі, флуктуації густини ρ, тиску Р та температури Т обернено пропорційні кореню з числа частинок газу

1N .

Звідси видно, що для середовища з великим числом частинок флуктуаціями термодинамічних величин можна нехтувати.

§ 81. Розподіл Максвелла для частинок за швидкостями

1. Закон розподілу частинок (атомів) ідеального газу за швидкостями визначає їх стаціонарний розподіл за швидкостями в умовах термодинамічної рівноваги й відсутності зовнішнього силового поля.

Розподіл частинок ідеального газу за абсолютними значеннями швидкостей V дослідив Максвелл. Він визначив, що число частинок, швидкості яких лежать в інтервалі (V; V + dV) представляється виразом

dn = n0dP,

де n0 загальне число частинок, а їх частка

dP = dn = f (V)dV n0

є ймовірністю такого розподілу із густиною

mV2

f (V) = Ae 2kT V2 .

В цьому виразі m маса частинки, k стала Больцмана, Т температура, V= Vx2 + Vy2 + Vz2 . Стала нормування А визначається так

 

3

 

 

m

 

 

2

 

A = 4π

 

 

.

 

 

2πkT

 

В.М.Клименко. Статистична фізика

106

 

 

Закон Максвелла можна застосовувати і для скінчених малих величин

n та V

 

 

 

 

n = n f (V)V.

 

 

2. Максимум густини розподілу f(V)

знаходиться в точці V=Vнв, яка

називається найбільш імовірною швидкістю і

 

 

 

V

=

 

 

2RT

=

 

2kT

1.41

kT .

нв

 

 

 

µ

 

 

 

m

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

Дійсно,

 

 

 

 

 

 

mV2

 

 

 

 

 

df (V)

 

 

 

 

 

 

2 2mV

 

 

 

 

 

 

 

dV

 

= Ae

2kT 2V V

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kT

Прирівняємо похідну 0 і знайдемо Vнв

 

 

 

 

 

 

 

 

 

mV2

 

 

 

 

 

 

 

 

 

1

 

нв

= 0

 

 

 

 

 

 

 

2kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

=

 

2kT =

2RT ,

 

 

 

 

нв

 

 

m

 

µ

 

 

 

 

 

 

 

 

 

 

 

 

що й треба було довести.

Вираз dP = dn = f (V)dV значно спрощується, якщо ввести змінну n0

u =

V

, тоді V = uV = u

2kT = u

2RT , dV =

2kTdu ,

 

 

 

нв

m

µ

 

 

 

m

 

Vнв

 

 

 

 

 

dn = n f (u)du ,

f (u) = 4

eu2 u2 = 2.26eu2 u2 .

 

 

 

π

 

 

 

 

 

Графік густини розподілу f(u) по u представлено на мал.39. Максимум

густини такого розподілу приходиться на u=1.

 

 

 

 

3. Середня арифметична швидкість

 

 

 

визначається як

i=N

 

 

 

 

 

 

 

1

 

 

 

 

 

V

=

Vi ,

 

 

 

 

 

 

N i=1

для дискретних значень V та

V = Vf (V)dV,

0

для неперервних значень V.

Розрахунки середньої швидкості частинок ідеального газу за розподілом Максвелла дають

V =

8kT

=

8RT

1. 60

kT .

(1)

 

πm

 

πµ

 

m

 

 

 

В.М.Клименко. Статистична фізика

107

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дійсно,

 

 

 

 

 

 

 

 

mV2

 

 

 

 

 

 

 

 

 

 

mV2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Vf (V)dV = VAe

 

 

V

2dV =

1 AAe

 

V2dV2

 

V

2kT

2kT

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

2

0

 

 

 

 

 

Зробимо заміну змінної

 

mV2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x =

, V

2

=

 

2kT

 

x ,dV

2

=

2kT

dx

 

 

 

 

2kT

 

 

m

 

 

 

m

 

i тепер

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m 3/ 2

 

 

3/ 2

 

 

2kT

2

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xe

 

 

 

 

 

 

V =

2π

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2kT

 

 

 

 

 

 

 

 

m

 

0

 

 

 

 

 

 

 

V =

8kT

1.60

kT .

(2)

 

 

 

πm

 

 

m

 

 

 

Ми проінтегрували

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I = xexdx

 

 

 

по частинам:

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I = −x dex =xex

0

+ exdx =0-ex

0

=1.

0

 

 

 

0

 

 

 

4. Експериментальна перевірка розподілу Максвелла була зроблена лише через 60 років після його встановлення О. Штерном та іншими дослідниками. Зокрема, найбільш простий та надійний дослід провів Б. Ламмерт. Вузький молекулярний пучок М вирізувався системою діафрагм в парах металу, які виходили з молекулярної пічки посудини, в якій випаровувався рідкий метал. Цей пучок проходив через два диски з вузькими радіальними прорізами, паралельно осі, на якій знаходились диски. Напроти пучка розміщувався уловлювач. Прорізи були повернуті один відносно другого на

кут ϕ. При нерухомих дисках пучок не потрапляв в уловлювач, а при обертанні з кутовою швидкістю ω вони проходили вільно. Це означало, що час повороту дисків на кут ϕ дорівнював часу прольоту частинок відстані L

між дисками: t = ωϕ = VL , де V середня швидкість молекул при заданій температурі Т. Звідси маємо

В.М.Клименко. Статистична фізика

108

 

 

V = L ωϕ .

Середню квадратичну швидкість

с =

V2 ,

 

де V2 = V2f (u)dV розрахуємо за

допомогою теореми Больцмана про

0

рівнорозподіл енергії за ступенями свободи для вільної частинки. Згідно цієї теореми середня кінетична енергія теплового руху частинки з масою m дорівнює

m2V2 = 32 kT

і звідси

с= 3mkT = 3RTµ 1.73 kTm .

§82. Барометрична формула

Барометрична формула

mgh

P = P0e kT

визначає залежність тиску атмосфери Р від висоти h над Землею. Щоб знайти цю формулу приймемо, що на поверхні Землі атмосферний тиск P0 , а

на при поверхневих висотах температура є сталою (величина зміни Т з висотою мала порівняно з Т на поверхні Землі). Газ атмосфери можна описати рівнянням КлапейронаМенделєєва, визначивши з нього, наприклад, густину повітря

ρ = RTPµ .

Нехай на висоті h тиск Р. При збільшенні висоти на dh тиск зменшиться на величину тиску стовпа повітря dP=-ρgh. П і с л я підстановки виразу для густини одержимо

dP = − PRTµg dh = − PmgkT dh .

З останнього виразу можна одержати диференціальне рівняння для тиску з розділеними змінними Р та h

dPP = − mgkT dh .

Після інтегрування маємо

lnP = −mgkT h + C.

В.М.Клименко. Статистична фізика

109

 

 

Сталу С знайдемо з умови, що на висоті h = 0 тиск P = P0 і остаточно

mgh

Р = P0 e kT .

Одержаний вираз називається барометричною формулою. За цією формулою можна досить точно вимірювати висоту h по значенню тиску Р, що вимірюється барометром. Проградуйований у значеннях висоти, цей барометр називається альтиметром.

Приклад 1. Розрахувати масу атмосфери Землі, користуючись барометричною формулою.

Дано.

µ = 0.029 кг, R =8.31

Дж

 

 

 

, R

 

= 6.4 106 м, P =1.01 105

Па, g =9.8

м

.

моль К

 

 

Обчислити: М ?

 

 

З

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

с2

 

 

 

 

 

 

 

 

Розв’язок.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dm = m0ndV = m0 n0e

m0g

x 4πRЗ2 dx

 

 

 

 

 

kT

 

 

 

 

 

Зробимо заміну змінної

 

 

 

 

m0g

 

 

 

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y =

x, dx =

 

dy.

 

 

 

 

 

 

 

m0g

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

kT

 

 

 

 

P0

 

 

 

 

 

 

dm = m0n0e

y 4πRЗ2

 

dy = 4πRЗ2

eydy .

 

 

m0g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

P0

 

 

 

 

 

 

 

 

P0

 

 

0= 4πRЗ2

P0

 

 

 

M = −4πRЗ2

 

eydy = −4πRЗ2

ey

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

0

 

 

 

 

 

 

 

 

 

g

 

g

 

 

Остаточно

M =5.3 1018 кг.

Для порівняння маса Землі складає MЗ =5.98 1024 кг, тобто в 1ю2 106 разів більша за знайдену масу атмосфери.

§ 83. Розподіл Максвелла - Больцмана

1. Розподіл Максвелла для частинок ідеального газу за значеннями кінетичної енергії має вид

dn = 2n (kT)

3

e

Ek

Ek

dEk ,

kT

2

 

 

 

 

 

 

π

 

де dn число частинок із загальної їх кількості n, кінетичні енергії яких лежать в інтервалі (Ek ;Ek + dEk ).

Цей розподіл можна одержати, якщо в розподілі по швидкостям провести заміну змінної

nx,y,z

В.М.Клименко. Статистична фізика

110

 

 

 

 

 

 

 

 

 

 

 

 

Ek =

mV2

,

V2 = 2Ek , dV = dEk

=

dEk .

 

2

 

 

 

m

 

 

 

 

 

 

mV

 

2mEk

Підставивши ці вирази в розподіл Максвелла, одержимо

 

 

 

 

m

3

 

 

Ek

 

2E

 

dE

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dn = 4π

 

 

 

 

e

 

kT

 

k

 

k

,

 

 

 

 

 

 

 

 

2πkT

 

 

 

 

 

 

 

 

 

m

2mEk

 

і остаточно

 

 

 

 

 

 

 

 

 

Ek

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

Ek

 

 

 

 

dn = 2n (kT)

 

 

e

 

kT

 

dEk ,

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

2. Зваживши на те, що P=nkT, з барометричної формули одержимо розподіл Больцмана для частинок у потенціальному полі

Eп

n = n0e kT ,

де Еп = mgh потенціальна енергія частинки в полі тяжіння Землі. Обидва розподіли можна об'єднати в один розподіл Максвелла-Больцмана, який визначає число частинок в одиниці об'єму, що знаходяться в

потенціальному полі з енергією Еп і мають швидкість в інтервалі

[V,V+dV]. Цей розподіл Максвелла - Больцмана має вигляд

 

 

 

 

 

 

 

3

 

 

E

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

dN = N0

 

 

 

 

 

e

 

kT dVx dVx dVz dxdydz ,

 

 

 

 

 

 

 

mV2

2πkT

 

 

 

 

 

 

 

де E = Eп +

повна енергія частинки. Потенціальна енергія частинки

 

2

 

 

 

= E

 

 

(x,y,z) , а V2 = V2

+ V2

+ V2 .

залежить від координати E

п

п

 

 

 

 

 

 

 

 

 

 

x

y

z

§ 84. Статистичний розподіл частинок за напрямком руху

Хаотичність теплового руху n частинок ідеального газу означає їх рівнорозподіл руху за напрямками швидкостей, а саме: усі напрямки теплового руху рівно ймовірні. Якщо в простір частинок увести декартову систему координат із 3-ма осями x,y,z і 6-ма напрямками по ним (±), то в кожному з них рухається

= 16 n

частинок. Наприклад, якщо частинки мають середню швидкість теплового руху V, то в додатному напрямку осі ОХ через перпендикулярну до неї площину S за час dt пройдуть всі частинки, що знаходяться на відстані dL = V·dt і їх загальне число становить

dN+ = nx dL S = 16 n V S dt .

В.М.Клименко. Статистична фізика

111

 

 

§ 85. Ефективний діаметр та переріз розсіювання

Ефективний діаметр dеф частинок ідеального газу є найменша

відстань між центрами двох частинок, які здійснюють центральний пружний удар. За порядком величини

dеф 1010 м.

Таблиця 5. Значення ефективного діаметра деяких молекул

 

Газ

 

Азот

 

Аргон

Кисень

Водень

 

Гелій

 

 

d·1010 м

 

3,75

 

3,64

3,61

2,74

 

2,18

 

Ефективний газокінетичний поперечний переріз σеф дорівнює площі

круга з радіусом

dеф :

σеф = πdеф2 .

Це означає, що

дві частинки, які

одночасово проходять через σеф , обов'язково співударяються.

§ 86. Число зіткнень частинки за одиницю часу

Число співударянь частинки із середньою швидкістю V при концентрації n за 1с дорівнює

Z = 2nσефV .

Щоб одержати цей вираз, розглянемо траєкторію, що її опише частинка ідеального газу за час dt. Вона являє собою деяку ламану лінію, у вузлах якої відбулися зіткнення з іншими частинками. ЇЇ довжина становить dL=Vdt, де V є середня швидкість частинок. Побудуємо на цій ламаній, як на осі, циліндр з основою, рівною ефективному перерізу σеф . Приймаючи до уваги

властивість ефективного перерізу, можна стверджувати, що розглядувана частинка обов’язково співудариться з усіма нерухомими частинками

dN = ndLσеф = nVσефdt ,

які знаходяться в циліндрі з концентрацією n. Але

частинки рухаються і

тому замість швидкості V потрібно взяти відносну середню квадратичну

швидкість Vв, яку розрахуємо так

r

r

 

 

 

 

 

 

r

 

r

 

 

 

 

V =

= V2 + V2

 

 

 

 

 

 

 

(V

V )2

2V V .

 

 

в

1

2

 

1

2

 

r

1

 

2

 

 

 

 

 

 

 

r

r

 

 

вектор

 

 

В записаному виразі Vв = V1

- V2

 

відносної швидкості двох частинок. З достатньою

точністю можна покласти, що

 

V2

=

V2

= V2 , а

 

r r

 

 

 

 

 

 

 

1

 

 

 

 

2

 

r r

 

 

 

 

 

 

 

 

 

 

 

 

 

V1V2 = 0 тому, що кут α між векторами V1 i V2

змінюється в межах від 00

до 1800 і

 

= 0. В

cosα

результаті

маємо,

що

Vв = V

2

 

 

 

і

 

число

В.М.Клименко. Статистична фізика

112

 

 

співударянь за час dt дорівнює dN = 2nVσефdt , а число співударянь за одиницю часу дорівнює

Z = dN

= 2nVσеф.

dt

 

§ 87. Середня довжина вільного пробігу

Середня довжина вільного пробігу λ є середня відстань між двома послідовними співударяннями частинок, що рухаються рівномірно й прямолінійно. Величину λ можна обчислити, виходячи з того, що на довжині траєкторії dL=Vdt відбувається Zdt співударянь і тоді середня довжина вільного пробігу між співударяннями становить

λ =

dL

=

Vdt

=

1 .

 

Zdt

 

2nVσефdt

 

2nσеф

Приймаючи до уваги, що n = kTP , можна записати

 

1

 

kT

λ =

2 nσеф

=

2 πdеф2 P .

§ 88. Явища переносу

Явища переносу (дифузія, теплопровідність, внутрішнє тертя) полягають у виникненні напрямленого переносу в газах маси (дифузія), внутрішньої енергії (теплопровідність) та імпульсу направленого руху відповідно.

Явища переносу виникають при просторовій неоднорідності концентрації n, температури Т і швидкості напрямленого руху частинок u. Вирівнювання значень n, T, u в об'ємі відбувається за рахунок теплового руху частинок, що характеризується середньою швидкістю V.

1. Експериментально визначено, що при одновимірній дифузії однорідного газу, число частинок N, що переносяться за час dt через елементарний поперечний переріз S дорівнює

N = −D dndx Sdt ,

де D коефіцієнт дифузії.

2. Експериментально визначено, що для ідеального газу в одновимірному стаціонарному випадку кількість тепла dq, що переноситься в напрямку ОХ нормально до площадки dS за час dt описується рівнянням Фур'є:

В.М.Клименко. Статистична фізика

113

 

 

dq = −χ dTdx Sdt ,

де χ коефіцієнт теплопровідності.

3. Сила внутрішнього тертя F, що діє на площадку S поверхні шару Ох при просторовій неоднорідності швидкості напрямленого руху u=u(x) OX задається рівнянням Ньютона

F = −ηdudx S,

де η коефіцієнт в'язкості.

§ 89. Дифузія

Для визначення коефіцієнта дифузії D, розглянемо одновимірний рух частинок, як показано на Мал.42. Нехай концентрація частинок є функцією координати х: n=f(x), а число частинок, що рухаються в напрямкові осі ОХ дорівнює

nx = 16 n = 16 f (x) .

В додатному напрямку ОХ через S ОХ, що має координату х, за 1с пройде

dN+ = 16 V S f (x − ∆x)

частинок, де величина х задає координату точки, починаючи з якої частинки без співударяння пройдуть через S. Такою величиною є середня довжина

вільного пробігу х = λ. Аналогічно в протилежному напрямкові пройде

 

 

dN=

1

V S f (x + ∆x)

 

 

 

 

 

6

 

 

 

частинок, де також х = λ. З визначення похідної

 

df

= f'(x) = lim

f(x +

xxf(x)

 

dx

 

x

 

 

 

x0

для достатньо малих х, що задовольняють певній точності обчислення похідної, можна записати

f '(x) f (x + ∆x) f (x)

x

і звідси

f (x + ∆x) = f (x) + f '(x)x ,

або в загальному виді значення функції в точці х±∆х можна записати через значення функції та її похідної у точці х таким чином

f (x ± ∆x) = f (x) ± f '(x)x ,

причому

Соседние файлы в папке Физика