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

книги из ГПНТБ / Беккер Р. Теория теплоты

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
15.46 Mб
Скачать

ние

зависит

только

от величины

интервала

времени

t'-t:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л (/) Л (О - Ф ( / ' — 0-

 

 

 

(81.2)

 

Введенная, таким

образом,

корреляционная

функция

Ф

(s) при s = 0 равна

Л 2 , при возрастании s она

быстро

 

 

 

стремится

к нулю. Исходя

из до­

 

 

 

пущений

о т, необходимо

потре­

 

 

 

бовать,

 

чтобы

функция

O(s)

 

 

 

практически уже была равна ну­

 

 

 

лю для значений s, еще малых по

 

 

 

сравнению с т.

 

 

 

t, t'

 

 

 

Поэтому

в

плоскости

 

 

 

(рис.

115) подынтегральная фун­

 

 

 

кция

в уравнении (81.1)

сущест­

Рис.

115. К расчету двой­

венно отличается от нуля лишь в

узкой

полоске вдоль

диагонали.

ного интеграла в уравне­

Если

в

качестве

новых

перемен­

нии

(81.1).

 

 

 

 

ных

интегрирования

ввести

 

 

 

 

 

 

 

S =

 

f—t;

 

 

 

 

 

 

U =

t'

+

t,

 

 

 

 

 

то получим

dtdt'

dsdu

и

 

 

 

 

 

 

 

 

 

 

^du

J O ( s )

ds,

 

 

 

 

причем ввиду быстрого стремления функции Ф к нулю с ростом s интегрирование можно производить по полоске (s от со до -j-co; и от 0 до 2т).

Таким образом, просто получаем

 

+ о

 

- i =

0(s)ds=A2a,

(81.3)

X

J

 

где малое по сравнению с т время а может быть обозна­ чено как «время когерентности» функции A (t). Оно яв­ ляется мерой ширины «полосы когерентности» на рис. 115. Интегралы (79.4)

о(0 = о 0 в - р ' + е " в ' J A ( S ) e M d |

о

430

основного уравнения v = — . A (t) можно преобразо­ вать, и не используя скоростей Gx , непосредственно при помощи корреляционной функции (81.2). После возведе­ ния в квадрат и усреднения по многим частицам из урав­ нения (79.4) следует:

 

 

t t

 

 

v*(f) = vl e-W

+er-W [\АЦ)А(V)

em+t)dgd£'.

 

 

о 0

 

 

Используя

обозначения

Л (g) Л (£') =Ф(£ ' — £ ) и

£ = s ,

=

приводим

интеграл в пашем прибли­

жении к виду

 

 

 

 

2/

+ 0 0

 

 

+оо

± - j > d a \

G>(s)ds = ± ( e m -

1) f 0(s)ds,

О

СО'

 

ОО

следовательно,

 

 

 

 

При i->oo правая часть уравнения должна стремить­ ся к выражению kT/tn, в связи с чем для корреляцион­ ного интеграла получим:

Г <t>{s)ds = —$,

 

 

(81.4)

J

Ш

 

 

 

что можно было бы установить еще из уравнений

(81.3)

и (78.5). Корреляция скорости v(t)

вытекает непосред­

ственно из уравнения (79.4) в результате

умножения на

Оо и усреднения по многим частицам:

 

 

v(t)v(t') =

~v~2e-m,

 

 

или, в несколько иной записи,

 

 

 

 

о (0 о (О = о 5

в _ р к " г

' .

 

(81.5)

Отсюда особо простым способом определяется сред­

нее квадратичное смещение: имеем

х=

[w(|)rf|,

следо-

 

 

 

о

 

вательпо,

li

х2 = j ' j" v(l)v(l')d£,d%.

6 6

43!

После усреднения и использования уравнения (81.5) при £—^'—s; %-{-%'= и для больших времен ( 6 ^ 1 ) по­ лучаем:

~x*

= —\du

f e-®sXds=~t=2

 

— mBt=2kTBt.

 

2 J

J

P

m

0- c o

Сучетом kTB — D приходим к хорошо известному

выражению х2 = 2Dt.

Если же условие В t^> 1 не выполнено, то в соответ­ ствии с (81.5) необходимо строго рассчитать величину

t t

x*^v2 f f,

о о

Это делается без затруднений и приводит к резуль­ тату

x 2 = 2 D

1 (1-

Это выражение совпадает с полученным выше по ме­ тоду Ланжевена и уже детально обсужденным выраже­ нием (77.6)

82. ДИФФЕРЕНЦИАЛЬНОЕ УРАВНЕНИЕ ФОККЕРА — ПЛАНКА

а) Без внешнего силового поля

Еще в § 79 статистическим методом была рассчитана ве­

роятность

W(v,

t, do) того, что скорость частицы ко вре­

мени i лежит в интервале v,

dv, если к моменту времени

/ = 0 она

имела

скорость v0.

Для решения подобной за­

дачи изберем теперь совсем другой путь, выводя диффе­ ренциальное уравнение для w (v, i) из основного урав­

нения v = — fiv-r-A(t). Пусть

 

q>(0,Ti)*l

(82.1)

представляет вероятность того, что по истечении време­ ни х скорость частицы с начальной скоростью v лежит

между

У + Т ] и v-\-r\-\-dr\.

Тогда для распределения ско­

ростей

к моменту

времени

t-\-x справедливо выражение

 

w{v,t-\-x)=

j w(v

т 1 , ^ ф ( у r|,r])dr).

(82.2)

432

Пусть теперь время т настолько мало, что функция ф(у, г)) значительно отличается от нуля только для та­ ких т], значение которых мало по сравнению с гл. В та­ ком случае можно разложить:

 

w (v, t + т) =

+~~

 

 

 

 

 

-dw

 

1

.д*а>]

 

 

 

j " w (v, f)

ч dv

 

— ч'

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

2

dv*

 

 

 

 

 

X Гф (v, Т))

г|

дф

 

i

1 2 д2 ф

dr\.

 

(82.3)

 

 

 

aw

 

т]

 

Если

ограничиться

 

 

 

 

2

 

dv2

 

по ц,

то в

квадратичными

 

членами

уравнение (82.3) войдут лишь следующие интегралы:

Ф (v, ту) dt) = 1;

т]ф (v, т)) dr) =

т);

 

Л ^ ^ ' ^

 

*1 = 5";

 

^ Л

 

ф (У, Л) dr\ = T i

 

;

|

^ д ф

К

п)

=

 

 

 

 

 

2

 

 

2

 

 

 

dy]

 

 

 

 

 

 

 

 

 

 

 

 

 

да

 

 

 

 

 

 

 

 

 

dv2

 

 

 

dv2

 

 

 

 

 

Эти средние значения можно заимствовать из преж­

них рассмотрений.

Согласно

обозначениям

уравнения

(78.2а) Г| представляет

 

собой

изменение скорости

за

время

т, т.е. т] = U J + I — V j .

 

 

Тем самым

при заданном

v =

= Vj

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т) = — г/Вт и

<Эг|

 

 

•Вт.

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

 

 

 

 

 

 

 

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

2

G

1

= ^ - ВТ, m

следовательно, 6\) 2 /ди=0 .

Таким образом, уравнение (82.3) принимает вид:

w(v,t + т) а ф , / ) =

ауВт + а б т —

Н

Вт—

 

 

 

 

да

dv2

и в пределах т->0 переходит в уравнение

Фоккера-

Планка:

 

 

 

 

 

 

dw

а

, о dw . akT д2 ш

(82.4)

— =

 

Rw + By

hp

dt

 

'

dv

m dv2

 

 

28-480

 

 

 

 

 

433

Если записать уравнение (82.4) в виде

 

dw

о

д l

 

, кТ dw }

,

 

 

, а п

 

^ 6

\wv-\

 

 

 

 

 

(82.5)

 

 

 

 

 

 

 

 

 

 

 

 

mo"

 

то можно заметить, что dwldt=Q,

если w = Ce

l k T . Рас­

пределение Максвелла,

естественно,

стационарно.

 

Полное интегрирование уравнения (82.4) удается пу­

тем изящной замены переменных. Сначала

вместо

v в

качестве переменной вводим

величину p — ve$'.

В

этом

случае

имеем:

 

 

 

 

 

 

 

 

 

 

 

 

 

w{v,t)

= w{pe-V,t)

 

=

Y{p,t),

 

 

 

следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

dw

_ dY

в/.

(Pw _

dW_^t

 

 

 

 

 

do

dp

'

dv2

 

 

dp2

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt ~~ dt

dp\

dt Jv

 

~~

dt

 

P

dp

'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

образом, из уравнения

 

(82.4)

получаем:

 

 

 

 

dt

К

V

 

 

.Эр2

 

 

 

 

 

 

где q = р —

 

 

 

 

 

 

 

 

 

 

 

 

\

m

 

 

 

 

 

 

 

 

 

 

 

 

Теперь подставим У—%е^{, откуда

Тем самым

имеем: dt

dt

1

 

 

 

dt

4

dp2 '

Наконец, путем

замены

е1^1

dt — dft вводим новую шка­

лу времени,

т.е.

полагаем:

 

Таким образом, уравнение (82.4) приведено к при­ вычной форме диффузионного уравнения

д%(Р, О)

Э 2 х ( Р . О )

3d

ар2

434

Стандартное решение этого диффузионного уравне­ ния имеет вид:

1

( Р - Р . ) 3

Х(Р,#;Ро) = -Цг е

«*> .

4я<70'

 

При f>->-0 это решение переходит в %(р, 0; ро)=6(р — —ро).

Если с помощью выражений

w = Xe*>; , = i L p . о

1

( ^ _ , ) :

т

р = ve^; р 0 = о0

снова перейти к первоначальным переменным, то полу­ чим:

m(v—v„ e ~ P Q 2

' O ;

1 2zikT ( l e ^ ) J

Таким образом, заново получен закон (79.7) прибли­ жения к распределению Максвелла, на этот раз исходя из дифференциального уравнения (82.4).

б) При наличии силового поля

Пусть теперь на частицу воздействует дополнительно внешняя сила К(х) с потенциалом Ф(х), следовательно, К(х) = —дФ/дх. Аналогично предыдущим параграфам отыщем вероятность f(x, v, t)dxdv нахождения частицы в интервале х, dx и v, dv. Если в данном случае выра­ жение

 

 

 

<p(x,v,t,i\)dldr\

 

(82.6)

представляет

собой

вероятность

изменения

х или v

в пределах

времени т на |, d\ и соответственно

т), dx\, то

применима

аналогия с уравнением

(82.2):

 

f (х, v,t

- f

т) =

j ' \ f (х — I, V — Т|, t) ф — g, V —

 

 

-'+

 

 

 

 

 

-x\,l,i\)dtdi\.

 

(82.7)

Так же как в (82.3) мы должны разложить подынтег­ ральную функцию по % и г) до квадратичных членов. По-

28*

4.35

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

 

О =

- 6

У +

 

 

 

+A(t),

 

 

 

 

 

 

т

 

 

 

 

т. е. после

интегрирования

 

по t

от /т до

( / + 1 ) т

[см. (78.2а)]:

 

 

 

 

 

 

 

 

 

- v. = -

ро, т +

 

т + G .

(82.8)

К тому же при т, меньшем

времени торможения 1/|5,

имеем:

 

 

 

 

 

 

 

 

 

 

 

x . ^ - x ^ v . x .

 

(82.9)

Используя замены £ =

X j

+ i

 

r

i = U j + i — v j t

из уравне­

ний (82.8) и

(82.9)

находим

теперь средние

значения

 

 

 

 

(усреднение производится по многим частицам с равны­ ми Vj = vl), ограничиваясь линейными членами по т:

Ч = ( - Р * - I - ^ L ) г,

\ т j dv

 

 

i f

= G2 =

 

рт, I = vx,

 

= т.

 

 

 

 

 

 

от

 

 

 

di>

 

 

 

Остальные

средние

значения, так же как £2 , |т) и т. п.,

при линейном приближении не вносят

никакого

вкла­

да

в т.

 

 

 

 

 

 

 

 

 

 

 

Следовательно, при разложении в ряд подынтеграль­

ного выражения (82.7) и выполнении

интегрирования

остается

 

 

 

 

 

 

 

 

 

f(*,v;

t + x) =

f { x , v - t ) - - ^ i + ± ^ L ^ - f ^ - .

 

 

 

 

 

 

dv

 

2 dv2

 

dv

 

 

Вводя указанные выше средние значения, произведя

деление на т и перейдя

к Игл т-*-0,

 

получим

уравнение

Фоккера — Планка в виде

 

 

 

 

 

 

^

^ -

P

^ +

^ ^

- ^

-

A

W

-

f .

(82.10)

at

 

Ov \

in

dv j

dx

 

m

dv

 

 

 

 

Для проверки этого уравнения прежде всего устано­

вим, что известное

распределение

 

 

 

 

 

 

 

 

 

 

0= + Ф (X)

 

 

 

 

 

 

 

f = Ce

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

436

действительно стационарно, т. е. обращает в нуль пра­ вую часть уравнения (82.10).

Если

из уравнения

(82.10) нужно получить

данные

об изменении во времени плотности р, то вначале

напра­

шивается

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

 

 

p(x,t) =

§f(x,v,f)dv.

(82.11)

Если выполнить это интегрирование в уравнении (82.10), то производные по v не дадут никакого вклада, так как \(х, v, t) при |и|->оо очень быстро стремится к нулю. Получим

 

=

— f vf(v, x)dv =

— [pv),

 

 

dt

дх J v

'

dx K

'

 

т. е. хотя

и правильный,

но, к сожалению, крайне три­

виальный

результат. Он выражает лишь

закон

сохране­

ния массы.

 

 

 

 

 

На основании

прежних результатов

можно

ожидать

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

 

 

 

 

К (х)

kT

dp

 

 

 

 

• р

дх

 

 

 

 

mfi

m(3

 

 

где B = l/mB — подвижность и D = £77mB—коэффициент диффузии; следовательно, надеемся получить дифферен­ циальное уравнение

=

? _ Ш £ !р

_ _ £ _ . * - ) .

(82.12)

dt

дх { пф

mf> дх j

 

Как установил Крамере1 , уравнение (82.12) действи­ тельно вытекает как приближенное уравнение из выра­ жения (82.10):

df

^ил

dv

 

L\(fV

+ lLJL

+ д Т

V

dt

\

dx I \

dx dv

m$

dx

 

K

(X)

f)~

\ m\°>

m(3 дар

(82.13)

 

 

пф

j

dx

 

Смысл данного преобразования выясняется при рас­ смотрении рис. 116: для того чтобы получить плотность р в точке х0, интегрировали функцию f(x, v) вдоль пунк­ тирной линии (параллельно оси v) согласно уравнению (82.11). Рассмотрим теперь наклонную прямую g, про-

1 Kramers Н. А. — «Physica», 1940, v. 7, p. 284.

437

ходящую через

точку

х0

и имеющую

уравнение

= Р(х0 х). Вдоль этой

прямой dv = —$dx.

При

движе­

нии вдоль этой прямой любая функция

ср(х,

v)

изменя­

ется на

 

 

 

 

 

 

 

 

 

 

dv

дер

дер

 

 

dv

дх

dv

dx

p

 

 

 

И Л И

W

 

V dv

ф .

P

dx

Следовательно,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

a

dx

<P =

P

dv

J

вдоль

g.

 

 

 

 

 

 

 

dv

 

 

 

 

 

 

 

 

 

 

Таким образом, если проинтегрировать уравнение

(82.13)

не вдоль пунктирной

 

вертикали,

а вдоль

наклон­

 

 

 

 

 

 

ной

прямой

g,

то

первое

 

сла­

 

 

 

 

 

 

гаемое в правой части урав­

 

 

 

 

 

 

нения

станет

равным

 

нулю.

 

 

 

\ 1

 

Выясним,

при

каких

обстоя­

 

 

 

\ l

 

тельствах

подобное

интегриро­

 

 

 

у

 

вание

по

 

наклонной

прямой

 

 

 

хЛ

 

может

рассматриваться

в

ка­

 

 

 

 

l\

 

честве

замены

интегрирования

 

 

 

 

 

 

при

постоянном

х.

 

 

 

 

 

 

 

 

 

 

 

 

Функция

f(x,

v,

t)

практи­

 

 

 

 

 

 

чески становится равной

нулю,

Рис. 116.

Замена

интег­

если v превышает максималь­

рирования

по

v

при

по­

ное значение v0.

На

прямой g

стоянном

х0

интегриро­

в точке w0

координата

х отсто­

ванием

вдоль

наклонной

ит

от х0 на

у0 /Р- Но

а0

пред­

прямой

g.

о =

р(х„ —

х).

ставляет

собой

длину

пути

 

 

 

 

 

 

 

 

 

 

 

 

торможения

частицы,

движу­

щейся

со скоростью

v0

И подвергающейся только

воздействию

трения р.

Если

теперь К(х)

и

[(х,

v)

настолько слабо зависят от х, что на расстояниях поряд­ ка пути торможения самой быстрой частицы они еще мо­ гут рассматриваться почти постоянными, то в самом деле можно заменить интегрирование остальных членов в уравнении (82.13) вдоль g интегрированием при по­

стоянном Хо. И в этом случае было бы точно

получено

требуемое уравнение (82.12).

 

При

термическом равновесии v0 имеет

порядок

У kT/tn.

Следовательно, справедливости

уравнения

438

(82.12) можно ожидать лишь в том случае, когда выпол­ нены следующие условия:

Ж г / ^ ±

< < К »

J

L

l /

^ i « / .

(82.14)

dx • V т

&

дх

у

т

р

 

Обсуждение

процессов, для которых данные

условия

не выполняются, можно найти в работе Крамерса и Чан-

драсекара

(см. сноску на стр. 429).

 

8». СПЕКТРАЛЬНОЕ

РАЗЛОЖЕНИЕ

СТАТИСТИЧЕСКОЙ

ФУНКЦИИ'

В предыдущих

параграфах мы охарактеризовали нере­

гулярное

ускорение

A(t)

через среднюю

квадратичную

 

т

 

 

 

 

скорость

^ A { t ) d t ^

и

корреляционную

функцию

A(t)A(t+s)

6

 

 

 

 

=<D(s).

 

 

 

Во многих случаях, в особенности если данное уско­ рение влияет на колебательную систему, целесообразно получить сведения о спектральном распределении. Для этой цели требуются некоторые математические вы­ кладки.

Под статистической функцией A(t) понимаем функ­ ции, подобные скорости v(t) при броуновском движении частицы или силе воздействия окружающих молекул на данную частицу. В частности, предполагаем, что среднее

значение A(t) за большой интервал времени

равно ну­

лю, или в более строгой

формулировке

 

 

Ita (

i f Л ( • ) » ) - о

 

 

 

о

 

 

и

далее, что хотя A(t)

и изменяется

крайне

хаотично,

но

за очень большие

промежутки

времени

характер

функции остается в сущности без изменения, так что при

не слишком

малых

значениях т среднее

квадратичное

значение —

j * A2(t)dt

= А2

не зависит от t\.

 

U

 

 

 

Для дальнейшего

расчета

предположим также, что

A(t) отличается от нуля лишь в пределах

конечного, да-

1 Becker R, — «Z. angew. Phys.», 1954, Bd 6, S. 23].

439

Соседние файлы в папке книги из ГПНТБ