Метода термодинамика(лекции)
.pdfней излучение. Нетрудно найти температуру, при которой плотность энергии
излучения |
σT 4 |
становится равной плотности |
энергии |
теплового движения |
||||
частиц 3kTn 2 . |
Из равенства σT 4 |
= 3kTn 2 |
следует |
T = 3 3kn σ . Для |
||||
достаточно |
разреженной плазмы, |
когда |
плотность |
частиц |
n ≈ 1021 м-3, |
|||
T ≈ 3 104 |
К. При температурах полной |
ионизации |
(T ≈ 105 |
К) плотность |
энергии излучения в плазме становится преобладающей, что приводит к нежелательному следствию– трудностиадиабатическойизоляциитакойплазмы.
Электростатический потенциал плазмы. Рассмотрим модель плазмы в предположении ее полной ионизации, т.е. состоящей из одинакового числа положительно и отрицательно заряженных частиц (N+ = N− = N / 2 ). При этом свободной энергией излучения пренебрегаем. Внутренняя энергия плазмы E складывается из кинетической энергии хаотического движения ее
частиц Eид = 3NkT 2 =CVT (внутренняя энергия идеального |
газа) и |
|
|
N |
|
энергии электростатического взаимодействия Eq |
= ∑ej ϕ(rj ) 2 |
|
|
j =1 |
|
N |
|
|
E = Eид + ∑ej ϕ(rj ) |
2 , |
(4.64) |
j =1 |
|
|
где ϕ(rj ) – потенциал поля всех остальных зарядов в точке нахождения j -го заряда. Энергия одного положительного (отрицательного) заряда в произвольной точке (например, r0 = 0 ), равна eϕ+(0) (−eϕ−(0) ). Тогда за суммарную энергию кулоновского взаимодействия можно принять
E |
q |
= |
1 Neϕ |
(0) − |
1 Neϕ |
(0) . |
|
(4.65) |
|||
|
|
2 |
|
+ |
|
2 |
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Согласно методу самосогласованного поля, распределение плотности |
|||||||||||
положительных и отрицательных зарядов определяется формулой |
|||||||||||
n (r ) |
= n exp ±e ϕ[n(r )] kT |
, |
(4.66) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где n0 = N /V , ϕ(r ) – |
потенциал в |
точке |
r , |
который |
порождается |
остальными зарядами. После подстановки (4.66) в уравнение Пуассона имеем
212
|
|
|
|
|
|
∆ϕ(r ) = −4πρ(r ) = 4πen0 |
eϕ kT |
−e |
−eϕ kT |
|
(4.67) |
|||||||||
|
|
|
|
|
|
e |
|
. |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Для |
|
достаточно |
разреженной |
изотропной |
|
плазмы |
(eϕ << kT , |
|||||||||||||
|
1 |
|
|
2 ∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
∆ = |
|
|
r |
|
|
) уравнение упрощается |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
∂r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r ∂r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
∆ϕ = χ2ϕ, |
|
χ2 |
= |
|
8πe2n0 |
. |
|
|
|
(4.68) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
kT |
|
|
|
|
|
|
Общее решение |
ϕ(r) =C e−χr r +C eχr |
r |
с учетом его ограниченности во |
|||||||||||||||||
|
|
|
|
|
|
|
1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
всем пространстве (ϕ(r → ∞) < M ) принимает вид |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
ϕ(r) =C |
1 |
e−χr |
|
r . |
|
|
|
|
|
(4.69) |
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Поскольку |
вблизи |
точки |
r = 0 , |
где |
находится положительный или |
|||||||||||||||
отрицательный заряд, потенциал ϕ(r → 0) = ±e r , |
то C1 |
= ±e . Потенциал, |
||||||||||||||||||
создаваемый всеми остальными зарядами в окрестности |
точки r = 0 , |
|||||||||||||||||||
определяется разностью между ϕ(r) и ±e / r |
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
ϕ (0) |
= lim e ±e−χr 1 = |
χe . |
|
(4.70) |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
± |
r →0 r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя (4.70) в (4.64) и (4.65), находим энергию кулоновского взаимодействия и внутреннюю энергию плазмы
E |
q |
= −Ne2χ, |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
E = |
3 |
NkT −Ne |
2 8πe2 |
|
N |
(4.71) |
||
2 |
kT |
V |
|
|||||
|
|
|
|
|
|
Отрицательное значение Eq обусловлено преимуществом сил притяжения,
поскольку каждый заряд находится в окружении противоположно заряженной среды (дебаевская экранировка).
4.10.3. Термодинамические потенциалы и параметры. При независимых переменных температуре T и объема V термодинамическим потенциалом плазмы есть свободная энергия F(T,V ) (3.40)
213
|
|
F = −T ∫ |
|
E |
dT = −T ∫ |
E |
+E |
|
|
|
||||||||
|
|
|
|
ид |
|
q |
dT = Fид +Fq , |
(4.72) |
||||||||||
|
|
|
2 |
|
2 |
|
|
|||||||||||
|
|
|
|
|
|
T |
|
|
|
|
|
T |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
3/2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2πmkT |
|
|
|
|
|
|
|
|
||||
где F |
= −NkT ln |
eV |
– свободная энергия идеального газа, а |
|||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||
ид |
|
|
|
N |
|
|
h |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F = −2 Ne2 |
8πe2N |
– |
|
свободная |
энергия |
кулоновского взаимодействия. |
||||||||||||
q |
3 |
kTV |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При температурах T ≥ 105 К необходимо к (4.72) добавить свободную энергию излучения (4.58)
F ′ = F −αVT3 .
Уравнение состояния ( p = −∂F∂V ), энтропия (S = −(∂F∂T )V ) и
теплоемкость (CV = −T (∂2F∂T 2 )V ) плазмы определяются формулами
|
|
p = NkT − |
1 Ne2 |
8πNe2 |
, |
|
|
|
|
(4.73) |
|
|
|
|
V |
3 |
kTV 3 |
|
|
|
|
|
|
S = (C |
) |
lnT +Nk lnV − |
1 Ne2 |
|
8πNe2 |
, |
(4.74) |
||||
|
V ид |
|
|
3 |
|
|
kT 3V |
|
|
||
C |
= (C |
) |
+ 1 Ne2 |
8πNe2 |
, (C ) |
|
= |
3 Nk . |
(4.75) |
||
V |
V ид |
2 |
kT 3V |
V |
ид |
|
2 |
|
|
Давление и энтропия плазмы меньше, чем у идеального газа, что обусловлено преобладанием в ней сил притяжения. При повышении температуры плазмы приходится затрачивать энергию не только на увеличение кинетической энергии хаотического движения ее частиц, но и на увеличение средней потенциальной энергии кулоновского взаимодействия между ними. Поэтому теплоемкость плазмы больше теплоемкости идеального газа. Кроме явления экранировки к специфическим свойствамплазмыотноситсятермическоеуравнениесостояния
|
|
A |
|
|
Ne2 |
8πNe2 |
|
|
|
|
+ |
|
|
, |
A = − |
|
|
. |
(4.76) |
|
|
|
|||||||
pV = NkT 1 |
|
|
|
|
|||||
|
|
|
|
|
|
3NkT |
kT |
|
|
|
|
V |
|
|
|
|
Здесь уже нельзя выделить, как в уравнении Ван-дер-Ваальса, разную
степень взаимодействия между частицами при изменении их плотности. 214
С помощью формул (4.73) и (4.75) находим следующие условия устойчивости разреженной плазмы
|
|
|
|
|
|
|
|
|
|
NkT |
|
|
1 |
|
|
2 |
8πNe |
2 |
|
|
|
|
|
|||||||
|
|
∂p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
= − |
|
|
|
|
+ |
|
Ne |
|
|
|
|
|
|
|
< 0, |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
||
|
∂V T |
|
|
|
|
V |
|
|
|
2 |
|
|
|
kTV |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(4.77) |
||||||||||||
|
|
T |
= |
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
> 0. |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
CV |
|
|
(C |
) |
|
|
+ |
1 |
Ne |
2 8πNe2 |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
kT 3V |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
V ид |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Отсюда следует, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
NkT |
|
|
|
|
ид |
|
|
|
T |
|
|
|
|
|
T |
|
||||||
|
∂p |
|
|
|
|
|
∂p |
|
, |
|
|
|
|
|
|
(4.78) |
||||||||||||||
|
|
|
|
|
> − |
|
|
|
|
= |
|
|
|
|
|
|
|
|
< |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
∂V T |
|
|
|
V |
|
|
|
∂V T |
|
|
CV |
|
|
|
|
(CV )ид |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
т.е. при одинаковых условиях плазма менее устойчива, чем идеальный газ.
4.10.4. Собственные продольные колебания . Вторая особенность плазмы также связана с дальнодействующим характером кулоновских сил взаимодействия. Речь идет о том, что флуктуация плотности электронов в плазме не релаксирует, как плотность в обычном газе, а колеблется с определенной частотой, зависящей только от концентрации электронов. Причина колебаний в том, что изменение плотности электронов
вкаком-либо месте плазмы связано с появлением объемного заряда, поле которого, действуя на смещенные электроны, приводит к появлению восстанавливающей силы, пропорциональной их смещению. Переменный характер силы вызывает вибрацию электронов с определенной частотой.
Выделим мысленно в плазме с концентрацией электронов n прямоугольный параллелепипед длиной dx и сечением S (объем dV = Sdx ) (рис. 4.12). Вследствие большой массы скорость движения ионов много меньше, чем у электронов, поэтому их можно считать неподвижными. Пусть
внекоторый момент времени электроны выделенного объема сместились по отношению к ионам на величину ξ(x). Тогда возникающий здесь объемный
заряд и его плотность становятся равными
215
|
Рис. 4.12 |
|
||
dq =enS ξ(x) −enS ξ(x +dx) ≈ |
|
|||
≈ −en dξ |
Sdx = −en dξ dV; |
|
||
dx |
|
|
dx |
|
ρ = |
dq |
=en dξ . |
(4.79) |
|
|
||||
|
dV |
dx |
|
Учитывая связь локального заряда с напряженностью электрического поля E
уравнением div E = = 4πρ, или эквивалентным ему для одномерного случая
dΕ dx = −4πen (dξ dx) получаем
|
|
E = −4πenξ . |
|
(4.80) |
|||
Из уравнения движения электрона (mξ = eE ) следует |
|
||||||
|
ξ + 4πne2 |
ξ = 0 , |
или ξ + ω2ξ = 0 . |
(4.81) |
|||
|
m |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
Эти формулы и определяют продольные ( E |
ξ ) электрические колебания в |
||||||
плазме с циклической частотой |
|
|
|
|
|
|
|
|
|
ω |
0 |
= |
4πne2 . |
|
(4.82) |
|
|
|
|
m |
|
|
|
|
|
|
|
|
|
|
|
Так, при |
концентрации плазмы |
n ~ 1016 |
м–3 её собственная |
частота |
|||
ω0 = 5 109 |
с– 1, что соответствует дециметровым волнам. Впервые наличие |
колебаний в плазме было установлено Дж. Релеем (1906) и независимо И.
Ленгмюром (1929 г.), получившим формулу для ω0 . Последняя была названа ленгмюровской частотой колебаний плазмы.
216
5.ТЕОРИЯ МАЛЫХ ФЛУКТУАЦИЙ
5.1.Определение и значение флуктуаций
Флуктуации – случайные отклонения наблюдаемых значений физических величин от их средних значений. Для макроскопических систем наблюдаемые значения физических величин с очень большой точностью совпадают с их статистическими средними, а сколько-нибудь значительные флуктуации встречаются редко. Так, относительная флуктуация аддитивной величины в системе из N частиц обратно пропорционально N . Этот результат является
следствием отдельных макроскопических малых частей системы. Исследование флуктуаций имеет принципиальное значение, так как
позволяет установить границы применимости термодинамических понятий и закономерностей. Кроме того, флуктуации являются существенной особенностью многих физических явлений. Например, с точки зрения термодинамики (рассматриваются только средние значения физических величин!), импульс макроскопического тела, находящегося в среде с температурой T , в состоянии равновесия равен нулю. Согласно же статистической физике, в соответствии с законом о равнораспределении энергии по степеням свободы средняя кинетическая энергия такого тела
должна быть равна 3NkT / 2 , а его третья часть совпадает с Nmvx2 / 2 . Таким образом, средний квадрат каждой декартовой компоненты скорости тела равен vx2 ≡ (∆vx )2 = kT / m . Это объясняет броуновское движение
взвешенных частиц, которое послужило одним из первых объектов изучения флуктуаций. Хаотические отклонения плотности частиц – причина таких явлений как критическая опалесценция, голубая окраска неба, шумы радиоаппаратуры и т.п. В работах А. Эйнштейна и М. Смолуховского (1905-
1906 гг.) построена количественная теория флуктуаций, которая оказалась в
217
|
|
|
|
|
1 |
|
e− |
(ξ−ξ0 )2 |
|
|
w(ξ)dξ = |
|
|
2k a dξ, |
(4.2) |
||||
|
|
|
|||||||
|
|
|
|
2πk a |
|
|
|||
|
|
|
|
|
|
|
|
||
а величина k a = (ξ −ξ0 )2 ≡ (∆ξ)2 |
|
есть |
дисперсия |
(средний квадрат |
отклонения); ее еще называют центральным моментом второго порядка. Если в разложении S(ξ) учитывать последующие члены, то распределение вероятностей флуктуации становится негауссовым и может быть нессиметричным. В этом случае требуется, помимо дисперсии, задание моментов более высокого порядка.
|
Из |
формулы следует, что отклонения |
физического параметра ξ от |
||||
|
|
гораздо чаще встречается в интервале |
|
|
|
|
, чем вне его. |
ξ = ξ |
|
ξ −ξ |
≤ (∆ξ)2 |
||||
0 |
|
|
0 |
|
|
|
Величина (∆ξ)2 дает представление о масштабе флуктуаций.
Знание среднего квадрата отклонения (ξ −ξ0 )2 позволяет найти дисперсию для любой функции случайной величины f(ξ). Учитывая малую величину отклонения, имеем ∆f = f(ξ) − f(ξ0 ) (∂f(ξ0 )∂ξ)(ξ −ξ0 ). Откуда следует
(∆f )2 = ∂f(ξ0 ) ∂ξ 2 (∆ξ)2 .
Заметим, что приведенная формула (5.2) применима к флуктуациям в системах с постоянной энергией.
5.2.2. Флуктуации в квазизамкнутых системах . Произвольную квазизамкнутую систему можно рассматривать как малую часть замкнутой системы или как подсистему, погруженную в термостат с постоянной температурой T0 . Считаем, что флуктуации происходят только в подсистеме, тогда как термостат остается равновесной системой. Состояние
подсистемы определяется некоторым внешним параметром λ. При переходе
220