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

книги из ГПНТБ / Кудрин, Л. П. Статистическая физика плазмы

.pdf
Скачиваний:
23
Добавлен:
22.10.2023
Размер:
17.6 Mб
Скачать

ftA F = ^ РД F 'j 1 —

3 V 2

) - 3/2

In 2

 

■/

\

+ ~r~ ~ T ^ ) ~i— • (31-33)

С учетом квантовой поправки уравнение ионизационного равно­ весия (31.22) модифицируется следующим образом:

2

1

£

 

ехр | — р [/ — 4с3 (2лпгР)

2]

Пе Па

— —

 

 

 

 

2яД2Р

 

 

 

+

У З

^1 —0,1 - ^ J .

(31.34)

Таким

образом,

в

результате учета квантовой

поправки в

ионизационном уравнении возник дополнительный экспонен­ циальный множитель, приводящий к уменьшению степени иони­

 

 

 

 

 

зации. Физически знак по­

 

 

 

 

 

правки

в

 

выражении

 

 

 

 

 

(31.33)

обусловлен

сле­

 

 

 

 

 

дующими

 

эффектами,

 

 

 

 

 

учитываемыми тремя сла­

 

 

 

 

 

гаемыми в множителе при

 

 

 

 

 

X//D. Два слагаемых обус­

 

 

 

 

 

ловлены

квантовыми

по­

 

 

 

 

 

правками к прямому элек-

 

 

 

 

 

трон-ионному и электрон-

 

 

 

 

 

электронному

взаимодей­

 

 

 

 

 

ствию, третий

член

опи­

 

 

 

 

 

сывает

обменное

элек-

 

 

 

 

 

трон-электронное взаимо­

 

 

 

 

 

действие

и также

приво­

 

 

 

 

 

дит

к эффективному

от­

Рис.

33. Зависимость

плотности

атомов

п а

талкиванию.

Хотя сделан­

ОТ ПЛОТНОСТИ ИОНОВ

tli при

постоянной

ная

поправка в принципе

 

температуре.

 

 

неверна

при

k/lD~ 1 и

дует

использовать

для расчета

 

формулу (31.34) не сле­

ионизационного

равновесия в

плотной плазме, тем не менее при т]кл —1 и Т (103—104)°К зна­ чения KJId таковы, что па заметно меняется с введением кванто­ вой поправки. Не исключено, что при этом производная (дпа/дпе) р сильно увеличится и весьма вероятно, что ее обраще­ ние в нуль окажется невозможным.

Если уравнения ионизационного равновесия, вытекающие из выражений (31.31) и (31.32), качественно дают изотерму типа 1 (рис. 33), то с введением квантовой поправки изотерма изме­

нит свой вид и скорее будет описываться

кривой типа 4. При

Г ^ Ю 3

°К в принципе может существовать

(при т)иЛ» 1)

крити­

ческая

температура Г,ф, ниже которой изотерма па(пе)

может

иметь перегиб. При этом два участка изотермы 2, удовлетво­ ряющие неравенству (31.26), описывают различные фазы: менее и более ионизованную. Если это так, то необходимо выполнение

320

условия равновесия фаз, т. е. в точках В и А необходимо равен­ ство температур, давлений и химических потенциалов:

=

Рл = Рв,

1^

=

^ -

(31.35)

Поскольку рассматривается система,

в

которой

нейтральная

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

= п %. По­

этому точки А и В, описывающие

две сосуществующие

фазы,

лежат на одном уровне.

 

 

 

 

 

 

Если расслоение на фазы действительно существует, то при

строгом решении задачи

должна получиться изотерма 3,

везде

' удовлетворяющая неравенству (31.26). Исследование точной об­ ласти устойчивости плотной кулоновской системы, даже далекой от вырождения, должно базироваться на правильном учете квантовых эффектов в термодинамических функциях, а также проводиться с учетом неидеальности нейтральной компоненты плазмы и взаимодействия атом — заряд. Только что рассмотрен­ ная интерполяция в область г)|(Л~ 1 довольно произвольна, и делать отсюда вывод о расслоении на фазы довольно опасно [5].

Интерес к

изучению

плотной плазмы неслучаен.

Один из

аспектов практического

использования такой плазмы

состоит

е следующем.

Г. Бете предложил использовать в качестве рабо­

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

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

тогда как ионы можно еще считать практически свобод­ ными [1]. Тогда, как показано в § 22, электронный газ образует кристаллическую решетку Вигнера и, если температура плазмы невелика, электроны осциллируют около равновесного положе­ ния, двигаясь в некоторой потенциальной яме (22.1).

Гамильтониан электронной подсистемы слабо отличается от оператора энергии, рассмотренного в § 22, который записан в предположении о наличии равномерного компенсирующего поля ионов. Тогда

j l

е2

2

За

е2

(31.36)

+ -----

Г2 _

2

 

ог3

 

Л>

 

11 Зак. 635

 

 

 

 

321

где р и т — соответственно

импульс и

масса электрона; г0 —

среднее расстояние между

заряженными

частицами; а — неко­

торый численный коэффициент, эффективно учитывающий кор­ реляцию электронов (а ^ 1 ). Гамильтониан (31.36) является осцилляторным, поэтому сразу же можно выписать значение собственной частоты колебаний

соо = е2/а тг\ — <о~р/За,

(31.37)

где о)р= (4лм,.е2/ т ) — частота ленгмюровских

колебаний для

электронов. Эту частоту запишем через

параметр задачи Т11;л:

Рйю0 = (2а)~'1г i]Jj ф10)~'и ,

(31.38)

где

 

 

/ 0 = те*/Н2 = 2 Ry 27

эв.

 

Если температура плазмы такова, что возбуждается лишь основное колебание со0 и вклад возбужденных состояний атома несуществен, то свободная энергия единицы объема рассматри­ ваемой системы имеет вид [1]

_ 1п ГJ L ( - И - У ‘

- фпа - пеIn Г—

-

я LЯв\2яЛ*р/ J

с L n e \ 2 n h * $ J

J

Г е2В

пе<х( —— ■ЗпеIn (рЙсо0), (31.39) \ Го

где па, ne=tii — соответственно равновесные плотности атомов, электронов и ионов; / — потенциал ионизации атома; е — осно­ вание натурального логарифма; М — масса атома, равная при­ мерно массе иона. Записав закон действующих масс ц„= m + Ц<? и выписав явные выражения для этих величин, после дифферен­ цирования выражения (31.39) по tij получим уравнение иониза­ ционного равновесия:

= (■2сф/0

Г]’/«

ехр

р/ — Y ay]«

(31.40)

*КЛ

 

 

Это уравнение представляет собой аналог известной формулы Саха в рассматриваемой модели. Величина —^-^-)а 11кл +

3 ^ р-1 имеет смысл эффективного снижения потенциала

ионизации, которое оказывается значительным в случае т)кл!Э>1. Оказывается, что функция (дпа/дпР) р имеет максимум.

Исследуем рассматриваемую систему на устойчивость, со­ гласно критерию (31.26). Вычислив производную от па по плот­ ности электронов пе при неизменной температуре, убеждаемся, что система термодинамически устойчива, если

ат)кл < 45/8 ж 6.

(31.41)

322

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

Разумеется, рассмотренная только что модель очень груба и на основании ее вряд ли можно делать категоричные утвер­ ждения об образовании двух фаз *. Совершенно очевидно, что если и невозможно точное термодинамическое рассмотрение плот­ ной плазмы, то изученная модель может быть, конечно, улуч­ шена. Необходимо учесть вандерваальсову неидеалыюсть ато­ мов, а также взаимодействие атомов с заряженными частицами. Отметим, что в работе [1], где рассматривается плотная неиде­ альная плазма, неправильно записан критерий устойчивости, что приводит авторов к ошибочному введению критической тем­

пературы Гкр— 0,1 /;

неправильно утверждается, что при г)—2

плазма неустойчива;

система рассматривается при

|3_13>Ясоо и

не учитываются возбужденные состояния осциллятора и т.

д.

В качестве задачи

можно предложить читателю

убедиться

в

том. что электронный кристалл Вигнера неустойчив. Для этого достаточно вычислить производную (dP/dV) р.

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

торые,

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

 

СПИСОК ЛИТЕРАТУРЫ

 

 

 

 

 

1. Алексеев

В. А., Велихов Е. П., Лопанцева Г.

Б. Докл. № 5М-74/102 на

 

симпозиуме «Производство электроэнергии с

помощью МГДГ». Зальц­

 

бург,

1966.

и теор.

физ.»,

1967,

т. 52,

2. Ключников Н. И., Тригер С. А. «Ж. экспернм.

3.

с. 276.

 

и теор.

физ.»,

1968,

т. 55,

Ключников Н. И., Тригер С. А. «Ж. эксперим.

4.

с. 1248.

 

физика. М., «Наука», 1964.

Ландау Л. Д., Лифишц Е. М. Статистическая

5.

Норман Г. Э., Старостин А. Н. «Теплофизика

высоких

температур»,

1968,

6.

т. 3, с. 410.

 

р. 75.

 

 

Berlin

Т.,

Montroll Е. J. Chem. Phys., 1952, v. 20,

 

 

7.

Dyson

F.,

Lenard A. J. Math. Phys., 1967, v. 8,

p.

423.

 

 

 

 

* Такое утверждение не следует даже из условия

(31.41).

 

 

11*

Г л а в а д в е н а д ц а т а я

ОБ ЭКСПЕРИМЕНТАЛЬНОМ ИССЛЕДОВАНИИ НЕИДЕАЛЬНОСТИ ПЛАЗМЫ

§ 32. О ДИАГНОСТИКЕ ПЛОТНОЙ ПЛАЗМЫ С ПОМОЩЬЮ РЕЗОНАНСНОГО РАССЕЯНИЯ у-КВАНТОВ

Рассмотрим экспериментальные возможности диагностики плотной системы кулоновских частиц. Плотная плазма — весьма трудный объект для экспериментального исследования. То, что она неудобна и для теоретического изучения, видно из предыду­ щего изложения. Дело в том, что неидеальная плазма обладает многими «неприятными» свойствами, сопутствующими неидеаль-

ности.

Так,

плазма достаточно высокой

плотности

п= 1020-^-

Ч-1021

см~г при относительно низких температурах ~ 1

эв (Р —

—500-^2000

атм) оптически непрозрачна.

Хорошо же

разрабо­

танные оптические методы измерения температуры и плотности в низкотемпературной плазме эффективны лишь в случае малых оптических толщин x = L/l, где L — характерный размер плаз­ менного слоя; I — средний пробег излучения в плазме. С увели­ чением оптической толщины спектроскопические измерения ста­

новятся возможными

лишь для плазмы-излучателя

умеренной

плотности

(~ 5 -1 0 19

см~3), обладающей четким линейчатым

спектром.

Примером

такой плазмы является

ртутная

плазма.

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

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

Однако в реальных экспериментах невозможно получить до­ статочно плотную плазму, однородную по температуре и плот-

32 4

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

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

Частично ионизованный газ может быть неидеальным в двух

отношениях.

Если плотность заряженной компоненты такова,

что энергия

кулоновского взаимодействия на среднем расстоя­

нии между

частицами сравнима с кинетической энергией или

превышает ее, то плазма неидеальна по заряженной компонен­ те. При давлениях 1000—2000 атм становится существенной неидеальность, определяемая вандерваальсовыми силами межмо­ лекулярного и межатомного взаимодействий. Уравнение состоя­ ния для нейтральной компоненты известно лишь в случае простейших газов. Существенный вклад в неидеальность при заметной степени ионизации плазмы дает взаимодействие ато­ мов с заряженными частицами. При больших давлениях прояв­ ляется неидеальность, обусловленная перекрытием волновых функций отдельных электронных состояний атомов.

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

Здесь предлагается метод измерения температуры в плотной оптической непрозрачной плазме с помощью ядерной резонанс­ ной у-флуоресценции [5]. Если ввести в плазму радиоактивные ядра, обладающие достаточно жестким у-излучением, можно получить информацию из любого объема оптически непрозрач­ ного газа. Например, для плазмы с плотностью п ~ \ О21 смгъ и атомной массой ядер А —150 при характерном размере плазмен­ ного слоя L » 1 см ослабление интенсивности у-квантов с энер­

325

гией Ev = 3004-500 кэв составит около 5%. Для такого излуче­ ния плазма является практически прозрачной.

Трудность задачи заключается в необходимости измерения малого эффекта, обусловленного относительно небольшим доп­ плеровским уширением у-линии. Например, для энергии излуче­ ния £^ = 300 кэв (при Л = 150, р-1~1 эв) отношение A /£v~10-6,

где А — допплеровская ширина линии.

Можно показать, однако, что использование явления ядерной резонансной у-флуоресценции, методы которой получили широ­ кое развитие в ядерной физике при измерениях естественной ширины и спинов короткоживущих состояний ядер [4], позволяет преодолеть эту трудность. Ядерная резонансная у-флуоресцен- ция состоит в возбуждении ядерного уровня под действием из­ лучения, обусловленного «разрядкой» ядра того же изотопа с той же энергией возбуждения. Резонансный характер явления обеспечивает высокую избирательную способность процесса, что и позволяет измерить малые эффекты, вызываемые допплеров­ ским уширением у-линии излучателя или рассеивателя.

Наблюдение же резонансного рассеяния у-квантов сущест­ венно осложняется сдвигом линии испускания относительно ли­

нии

поглощения на величину

 

 

EeK = 2EH = tfylM(*,

(32.1)

где

Еп — энергия отдачи ядра; Еу — энергия

у-кванта; М

масса ядра; с — скорость света. Энергию смещения для восста­ новления условия полного резонанса можно компенсировать не­ сколькими способами: 1) нагреванием источника либо рассеи­ вателя; 2) механическим смещением источника относительно рассеивателя; 3) использованием отдачи ядра в предшествую­ щем процессе распада ядра источника; 4) использованием ско­ рости, приобретаемой ядром источника при его образовании в результате ядерной реакции. Для измерения температуры плаз­ мы с помощью резонансного рассеяния у-квантов в принципе можно применить все указанные выше методы компенсации £ см. При этом плазма может быть использована как в качестве ре­ зонансного рассеивателя, так и в виде источника резонансного у-излучения.

Рассмотрим возможную схему опыта с плазмой в качестве резонансного рассеивателя на конкретном реальном примере ксе­ ноновой плазмы. Условие компенсации энергии отдачи (по спо­

собу /)

обеспечивается тем, что сама плазма уже

является до­

статочно

нагретым рассеивателем. В этом опыте

может быть

использована резонансная пара ядер 13Ч—131Хе.

 

 

Схема возможного опыта представлена на рис. 34. Коллими­

рованный пучок у-квантов

источника

1311—131Хе* с энергией

Еу =364

кэв

рассеивается в

некоторой

области

плотной ксено­

новой (|3|Хе)

плазмы. Такой

плазмой

может

быть, например,

326

плазма дуги высокого давления либо СВЧ-нагревателя. Форми­ рование направленного пучка и выделение резонансно рассеян­ ных под углом 0 у-квантов осуществляется системой свинцовых коллиматоров. Выделенные из локальной области плазмы кван­ ты регистрируются стандартной аппаратурой с амплитудным дискриминатором.

Резонансная пара 1311-—131Хе хорошо удовлетворяет требо­ ваниям, предъявляемым к источнику и рассеивателю в такого

рода сложных опытах. На рис. 35 приведена схема распада 13Ч. Изотоп 1311 обладает оптимальным периодом полураспада по отношению к p-излучению (7’1/2=8 дней), обеспечивающим до­ статочно высокую активность источника, а основным у-перехо-

дом при разрядке

возбужденных состояний 131Хе* является

переход 5/2+ -^-3/2+ с

энергией £ v = 364

кэв, 7’1/2=2,8- 10~п сек

и коэффициентом внутренней конверсии

а = 0,02. Интенсивность

у-излучения для всех остальных переходов составляет лишь 20% обшей интенсивности.

Если источник 1311 находится в твердом или жидком состоя­ нии, то интенсивность межатомных столкновений достаточно ве­ лика, т. е. Тет'Ст, где т — время жизни уровня. При этом ядра 131Хе * находятся в термическом равновесии и «забывают» об энергии отдачи при p-распаде, а линия испускания остается до­

статочно узкой.

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

327

монохроматическими у-квантами описывается формулой Брейта — Вигнера:

2/ i + l

л+2

ГУ

(32.2)

С +з (^)

2

( £ - £ а)2 + (Г/2)2 ’

2/о + 1

 

где 11 и /о — полные угловые моменты возбужденного и основ­ ного состояний ядра; Е и й — энергия и длина волны падающего монохроматического излучения; Еа— резонансная энергия; Гу

и Г — естественная и полная ширина линии соответственно. В нашем случае Г7— Г, поскольку коэффициент внутренней кон­

версии мал (0,02), и выражение

(32.2) можно переписать в виде

«Урез ( Е) =

 

(32.3)

 

4 ° маКС Еа)‘*+ (Г/2)2

ГД6 О^макс сечение процесса

при Е = Еа (полный резонанс).

Для рассматриваемого примера 0°макс= 2,8-10~20 см2/атом. Энер­ гия кванта, падающего на неподвижное ядро, Е и энергия кван­ та Е' относительно ядра, движущегося в направлении кванта со скоростью и, связаны соотношением

Е' = Е[\ + (и/с)],

(32.4)

328

которое с учетом максвелловского распределения для компо­ ненты скорости ядер в направлении источника приводит к сле­

дующему распределению у-квантов относительно рассеивателя:

^ £,= ^ Крхе ^ У К '

где допплеровская ширина линии поглощения

 

Аа =(£'/с)(2/Мра),/%-

 

 

 

(32.6)

ра==11Ша\ Та— температура

рассеивателя;

М — масса

ядра.

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

 

 

 

 

 

 

°Рез (Е, Г/А) =

f арез (,Е') Р (F ) dE’ =

o°aKc V (х, t),

(32.7)

где

 

 

 

 

 

 

х = (Е — Еа)/Г/2;

t = (Д/Г)2, ¥ (х, t) =

(я/2/)‘/г ехр (—х2/40-

После интегрирования с учетом упрощений, связанных с тем,

что Д /Г»1, выражение для

сечения рассеяния

принимает вид

Od (Е) = Q

Смаке

ехр [ - ( - Ц

^

) 2] •

(32.8)

Если учесть, что линия испускания также имеет допплеровскую ширину и распределение у-квантов источника описывается вы­ ражением

Nd (E) = — ^ — ехр

(

Е - Е еу -

(32.9)

У п Ае

Ае J .

 

 

 

где Ае= (Е/с) (2/Мре) 1/2— допплеровская ширина линии испу­ скания; Ее— энергия центра линии испускания; pes=ljkTe\ Те— температура источника, то получим формулу для среднего се­ чения рассеяния с учетом допплеровской ширины линий испу­ скания и поглощения:

-------- трг,

f oD (Е) Nd (Е) dE

Срез. D ( Е )

Т~Г

 

J Nd (Е) dE

сти

г V

л

ехр

(fi. -

Еар 1

(32.10)

2

макс

(Де +

Д2а)'/2

\ 2+

д« J '

 

 

 

 

Здесь Еа—Ее= Е 2/Мс2.

Если Да сравнимо с £ см, то резонансное рассеяние наблюдае­ мо и сильно зависит от температуры. Формулу (32.10) удобно переписать в виде

Y я

Мс2

 

 

X

Срез.D(£)

Мс2 2k (Та+ Те)

 

 

X ехр

0

Мс2

gexp(— g2), (32.11)

" м а к с А

 

Mc2-2k(Ta +

Te)

 

 

 

329

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