книги из ГПНТБ / Райзер Ю.П. Лазерная искра и распространение разрядов
.pdfном расположении светового луча разряд теряет стабильность и гаснет. Это связано с влиянием архимедовой силы, под действием которой плазма всплывает и выходит из области интенсивного поля. Возможно, здесь сказывается и влияние конвективных га зовых потоков. Вертикальная подача питающего излучения снизу вверх надежно стабилизирует разряд и снимает верхнюю по давлению границу его устойчивого существования. В этом случае плазма, всплывая вверх, попадает в область более интен-
|
|
|
Рис. 7.5. Зависимость от дав |
|
|
|
ления мощности, необходимой |
|
|
|
для поддержания непрерывно |
|
|
|
го оптического разряда гори |
|
|
|
зонтальным лучо!^ |
0 |
10 |
20 |
30 |
|
|
|
р, ат м |
сивного светового поля и снова разрастается, смещаясь вниз по лучу, так что разряд становится устойчивым.
На рис. 7.5 показаны измеренные пределы для существования устойчивого горения при горизонтальной подаче лазерного из лучения в аргоне и ксеноне. Нижние ветви соответствуют обычно му порогу, определяемому потерями, верхние — пределу устой чивости. Область устойчивого горения заключена между верхней и нижней кривыми. При вертикальной подаче излучения верхних кривых нет, остаются только нижние, и разряд горит при любом высоком давлении. Видно, что при повышении давления порог (по потерям) понижается, причем сначала резко, а потом очень мед ленно. Такое поведение порога очень похоже на то, которое на блюдалось в опытах [90] по поджиганию разряда в импульсном режиме (см. подраздел 23.1). И причины эффекта здесь, по-види мому, те же самые. При низких давлениях преобладают теплопроводностные потери и пороговая мощность уменьшается с ростом давления в соответствии с формулой (7.3). При высоких же дав лениях преобладают потери на излучение, которые возрас тают с давлением примерно так же, как и поглощение лазер ного излучения. Поэтому соотношение тепловыделения и потерь становится слабо зависящим от давления.
В работе [4] определялась температура плазмы. По штарковскому уширению — линии водорода измерялась плотность электронов. В аргоне при давлениях 4—16 атм полуширина линии составляла 100—130 А, чему соответствует плотность электронов
5-1017 i /см3. При р = 2 атм 3,5«Ю17 1 /см,3. В предположе нии существования термодинамического равновесия этим значе
232
ниям отвечает температура 23 000°. Легко видеть, что измеренная плотность электронов соответствует почти полной однократной ионизации газа и диапазон возможных вариаций температуры
сучетом различия температур электронов и ионов весьма огра ничен. Непосредственная оценка также показывает, что отрыв температур невелик. Таким образом, температура аргоновой плазмы при 2 атм составляет примерно 23 000°. Измеренная тем пература согласуется с расчетом, согласно которому температура плазмы оптического разряда при условии ее заметной прозрачно сти должна быть близкой к той, при которой ионизация (одно кратная) почти полная и коэффициент поглощения максимален. Плазма светится ослепительным белым светом; непрерывных источников света столь большой яркости раньше не было.
Вработе Фрэнцена [5] были воспроизведены опыты по полу чению непрерывного оптического разряда примерно в такой же постановке, при мощности лазера 0,6 кет. В статье Стюэрдинга [371 были сообщены сведения, почерпнутые в порядке личной информации от Конрада и Сциклауса, о том, что последние получили непрерывный оптический разряд в атмосферном воздухе
спомощью мощного газодинамического лазера на С02. Мощность была 100 кет, луч слегка фокусировался линзой с фокусным рас стоянием 15 ж в пятно площадью 10 см2, поток в фокусе составлял 10 кет!см2. Плазма двигалась от фокуса по расходящемуся лучу
клинзе, а затем останавливалась, как и в опытах [3, 4]. Длитель ность работы лазера и горения разряда не указывается. При таких больших диаметрах луча теплопроводностный механизм выноса энергии из разряда [1, 2] значительно уступает конвек
тивному, связанному с тепловым расширением плазмы12 [37].
28. Высокочастотный индукционный разряд
Основы современной техники индукционных разрядов высоко го давления были заложены в 40-х годах работами Г. И. Бабата [6]. Г. И. Бабат, создав высокочастотные ламповые генераторы с мощ ностями порядка сотни киловатт, получил безэлектродный разряд в воздухе при атмосферном давлении и вводил в разряд огромные по тому времени мощности в десятки киловатт.
Индукционные разряды нашли серьезное практическое при менение. На их основе созданы безэлектродные плазмотроны, которые имеют ряд важных преимуществ по сравнению с дуговы ми. Главное — в безэлектродных плазмотронах генерируется сте рильно чистая плазма, не загрязненная продуктами разрушения электродов. Кроме того, плазменная струя получается более ста
1 Примечание при корректуре.
2 Смит и Фаулер [39], фокусируя линзами с / ^ 14 см луч пятикпловаттного лазера непрерывного действия, получили оптический разряд в атмос ферном воздухе. Экспериментальный порог стабильного поддержания плаз мы — примерно 2 кет, что неплохо согласуется с предсказанной величиной.
233
бильной, и практически не ограничено время работы, тогда как в дуговых плазмотронах большой мощности быстро разрушаются электроды. Температуры плазмы обычно 6000—10 000° в за висимости от мощности и рода газа. Безэлектродные плазмотроны сейчас используют в плазмохимии, металлургии, для выращива
ния кристаллов тугоплавких веществ, |
для физических исследова |
|||
ний, например для изучения воздействия |
высоких температур и |
|||
плазмы |
на материалы. |
из |
главных вопросов, на |
|
28.1. |
Температура плазмы. Один |
|||
который должна ответить теория разряда, |
это какова температура |
|||
|
|
Рис. 7.6. Схема индукционно |
||
|
|
го разряда в трубке, |
вставлен - |
|
|
|
ной в длинный соленоид |
||
|
|
Справа — распределение темпера |
||
|
|
туры |
по радиусу, |
г0 — радиус |
|
|
|
разряда |
|
плазмы, |
от чего и как она зависит. Помимо непосредственного ин |
тереса, который представляет этот вопрос, его принципиальное значение состоит в том, что по известному распределению темпе ратуры в плазме можно вычислить любые другие характеристики разряда: его омическое сопротивление, мощность, вкладываемую в разряд, распределения полей, самоиндукцию.
Довольно полное представление об энергетических и электро динамических закономерностях индукционного разряда можно получить, рассматривая идеализированную задачу об одномерном стационарном статическом разряде, который горит в неподвижном газе внутри диэлектрической трубки, помещенной в бесконечный соленоид. Стационарность состояния обеспечивается тем, что все выделяющееся в разряде тепло отводится от плазмы путем тепло проводности в охлаждаемые стенки трубки (рис. 7.6). Такая зада ча была строго сформулирована в работе В. Н. Сошникова и Е. С. Трехова [7].
В условиях цилиндрической симметрии и стационарности про цесса радиальное распределение температуры Т (г) описывается
следующим уравнением |
баланса энергии: |
|
|
_ _ L -L rJ + а< |
£2> — Ф = 0, |
/ = — A,dT/dr, |
(7.4) |
где Ф — потери на излучение плазмы.
В соленоиде бесконечной длины магнитное поле имеет только осевую, а электрическое — азимутальную составляющие: H = H Z, Е = Е9. В случае мегагерцевых частот с хорошим приближением можно пренебречь токами смещения. Имея также в виду, что Е, Н ~ ехр (—iat), запишем уравнения Максвелла прямо для ком
плексных амплитуд, которые |
для краткости обозначим |
Е, Н : |
|
- dHldr = |
<зЕ, |
J- Л - г Е ^ icoHie. |
(7.5) |
234
Фактически, это не два, |
а четыре |
|
|
||||||
уравнения для четырех |
функций по |
ом-см |
|
||||||
ля, ибо действительные и мнимые |
|
|
|||||||
части Е и |
Н или же амплитуды и |
|
|
||||||
фазы являются |
функциями |
незави |
|
|
|||||
симыми. |
|
|
|
|
|
|
|
|
|
В равновесной плазме с выравнен |
|
|
|||||||
ным давлением проводимость о яв |
|
|
|||||||
ляется функцией одной |
только тем |
|
|
||||||
пературы, о = а (Т), |
причем высо |
|
|
||||||
кочастотная проводимость при ат |
|
|
|||||||
мосферном |
давлении |
практически |
|
|
|||||
совпадает с проводимостью для по |
|
|
|||||||
стоянного тока. На рис. 7.7 показа |
|
|
|||||||
ны кривые 0 (Т) для |
воздуха и ар |
|
|
||||||
гона. |
|
|
условия |
к |
системе |
|
|
||
Граничные |
|
|
|||||||
(7.4), (7.5) таковы. На оси, при г = |
О |
|
|
||||||
вследствие |
симметрии |
/ |
(0) = |
0, |
|
|
|||
Е (0) = |
0. У внутренней поверхности |
|
|
||||||
трубки |
при |
г = R |
поток |
тепла |
и |
|
|
||
температура |
подчиняются |
условию |
|
|
|||||
теплопередачи в стенки. Для просто |
|
|
|||||||
ты можно положить и так: |
Т (Л) |
= |
Рис. 7.7. |
Электропроводности |
|||||
= Тот, |
где |
Гст — какая-то |
неболь |
воздуха и |
аргона при атмос |
||||
шая по сравнению с температурой |
ферном давлении |
||||||||
плазмы температура |
у |
стенки. Маг |
|
|
нитное поле в непроводящей среде у стенки определяется «ампервитками» соленоида:
при г = R |
H{R) = H0 — ~ I Qn, |
(7.6) |
где / 0 — амплитуда тока |
в катушке, а п — число витков на еди |
|
ницу длины [8]. |
|
так как |
Амплитуды / 0 и Н 0 можно считать действительными, |
фаза тока произвольна.
Уравнения (7.4), (7.5) интегрировались численно [7], причем для упрощения решалась обратная задача. На самом деле в соот ветствии с постановкой задачи следует задавать ампер-витки соленоида и радиус трубки R. В результате решения должны опре делиться температура и магнитное поле на оси Т (0), Н (0). Авто ры же, задаваясь произвольно значениями Т (0), Н (0), начинали интегрировать уравнения от точки г = 0. На некотором радиусе температура Т (г) очень резко падала. Радиус R, соответствующий экстраполяции Т к нулю, и считался радиусом трубки, а по ампли туде поля Н в этой точке определялись ампер-витки.
На рис. 7.8 представлены результаты расчета для нескольких вариантов, которые иллюстрируют поведение температуры и по-
235
лей в разряде. Профиль температуры имеет характер плато с не большим провалом в середине и резким спадом у краев. Появле ние провала связано с влиянием потерь на излучение: тепло вы деляется только в периферийном слое плазмы из-за скинирования, а излучает весь нагретый объем.
Многочисленные измерения распределений температур и элек тронных плотностей по радиусу, выполненные различными мето дами, подтверждают такое поведение температуры и также свиде тельствуют о существовании небольшого провала на оси. На рис. рис. 7.9 и 7.10 для примера приведены результаты двух экспери ментальных работ.
Результаты численных расчетов многих вариантов и их разум ное согласие с экспериментом показывают, что изложенная выше постановка задачи правильно отражает существо процесса. Одна ко проникнуться пониманием природы закономерностей явления можно только путем пусть даже приближенных, но аналитиче ских решений. Значительный прогресс в этом направлении был достигнут в работе В. А. Груздева, Р. Е. Ровинского и А. П. Со болева [11], в которой рассматривалась та же самая задача и те же уравнения, но только без учета потерь на излучение (кстати, расчеты [7] показали, что последние существенны только при
236
Рис.7.9. Радиальные распределения температур (1, |
2, 4) |
и |
электронных |
концентраций (3, 5) в аргоне (1—3) и ксеноне (4, 5) [9]] |
|
|
|
Статический разряд в трубке радиуса 35 мм, на частоте 11,5 |
Мгц: |
1 ,2 |
— р = 1 атм |
мощности — 4,7 и 7,2 кет; з — 7,2 кет; 4,5 — р = |
1 атм., 6 |
кет |
Рис. 7.10. Изотермы в аргоне [10]
Разряд |
в потоке |
газа, р = |
1 атм. |
Радиус трубки 11 |
мм, |
трубка |
охлаждается |
водой; |
||
частота |
28 Мгц; |
мощность |
2,5 |
кет, |
расход |
80 смг/сек. |
1 — Т = |
9750°, 2 — |
9650°, |
|
|
|
3— 9400°, |
4 — 9000°, |
5 — 8700°, |
в — 8400° |
|
|
больших |
мощностях, |
когда температуры превышают примерно |
10 000°). |
В работе [И] |
был найден второй интеграл системы (7.4), |
(7.5), который в практически важном случае тонкого скин-слоя непосредственным образом определяет температуру плазмы через ампер-витки соленоида.
Опуская потери Ф в уравнении (7.4) и выражая член джоулева
тепла |
через |
дивергенцию |
потока электромагнитной энергии S |
|||||
(S = |
Sr) |
по |
общей формуле |
(6.26), получим первый |
интеграл |
|||
системы |
(7.4), |
(7.5) |
|
|
|
|
||
|
|
|
J |
+ S = 0; |
/ |
= / г> 0 , |
S C O , |
(1.7) |
где постоянная интегрирования, равная нулю, определяется гра ничными условиями на оси.
237
Уравнение (7.7) выражает тот очевидный факт, что в условиях стационарного процесса без потерь электромагнитный поток в каждой точке компенсируется тепловым, так же как и в целом поток электромагнитной энергии, поступающий в разряд от соле ноида, полностью отводится теплопроводностью в стенки трубки.
Но уравнения Максвелла (7.5), в которых опущен ток смеще ния, позволяют представить поток S в дифференциальной форме:
S |
с2 |
/Л№\ _ |
с2 |
, гг ,2 |
(7.8) |
|
32it2a |
\ dr / |
64я2б |
dr ' 1 |
|||
|
|
Если подставить это выражение, а также дифференциальное выражение (7.4) для / в уравнение (7.7) и умножить все на а, уравнение можно проинтегрировать, и мы получим второй интег рал системы, причем постоянная интегрирования определяется граничными условиями у трубки: при Т = Тст Н = Н0. Будучи отнесенным к точке г = 0, интеграл дает связь максимальной тем пературы плазмы (на оси) Тк = Т (0) с величиной магнитного поля на оси
Тк |
СЧ11 |
[ Я (О)]2 |
(7.9) |
|
J o{T)k{T)dT |
||||
64я2 |
н 2 |
|||
тст«о |
|
|||
|
1 о J |
|
В том практически важном случае, когда скин-слой у поверх ности плазменного столба оказывается тонким по сравнению с его радиусом (фактически с радиусом трубки R), т. е. поле не прони кает в глубь плазмы и Я (0) ж 0, уравнение (7.9) дает непосред ственную связь температуры плазмы Гк с током в соленоиде
Т.к |
(7.10) |
J а(Т)Х (Т) dT = с2Fill64it2 = (1>/2)2. |
|
Тст»0 |
|
Температура в этом случае не зависит ни от частоты поля, ни от радиуса трубки.
В работе [11] развит метод последовательных приближений для нахождения распределений температуры Т (г) и полей и вы числения мощности и в общем случае, когда Н (0) 0. Приведем численный пример. Для разряда в аргоне при атмосферном давле
нии и при |
10п = 13,3 а-в/см по расчету [11] |
получается Тк |
— |
|||
= |
8000°; на частоте |
12 Мгц и при радиусе трубки R = |
3,75 |
см |
||
в |
единицу |
длины |
разряда вкладывается |
мощность |
W |
= |
= 0,21 квт/см; радиус поверхности, где температура снижается
до такой величины, что проводимость уже очень мала, |
Т = |
4500°, |
||||
составляет |
г0 = 0,91 R. По |
известной формуле |
|
|
||
|
б = |
C / V 2ПОЮ= |
5 ,0 3 //Оом-гсм-^Мгц СМ |
|
(7.11) |
|
для толщины скин-слоя [8] |
проводимости схк = |
о (Тк) |
соответст |
|||
вует 6К = |
0,45 |
см. При |
/ 0п = 33 а-в/см |
получается |
Гк = |
238
= 10 000°, W = 1,1 квт/см, r0 = 0,98 R, 5 = 0,3 см. Эти резуль таты неплохо согласуются с экспериментом.
Обращает на себя внимание большое сходство соотношений (7.10) и (6.59). Последнее справедливо для такого же статического разряда, но поддерживаемого электромагнитной волной, которая описывается приближением геометрической оптики (6.28). Причи на сходства состоит в том, что в обоих случаях существует интег рал сохранения потока энергии (7.7) и (6.57) и уравнения Макс велла позволяют приближенно выразить поток электромагнитной энергии в дифференциальной форме. Дифференциальные выраже ния (7.8) и (6.28) соответствуют предельным случаям малых частот, когда несуществен ток смещения, и больших частот, когда спра ведливо геометрико-оптическое уравнение (6.28).
Физический смысл определяющего соотношения (7.10) стано вится наглядным, если связать температуру не с ампер-витками,
ас плотностью потока электромагнитной энергии S0, поступающей
вразряд от соленоида. Это можно сделать, но только приближен но. Уподобим разрядную плазму однородному проводнику радиу са г0 и проводимости п = пк, который отделен от стенок трубки непроводящим зазором,— это называют моделью «металлического цилиндра». (Из рассмотрения графиков Т (г) рис. 7.8 следует,
что такое приближение вполне приемлемо.) Если скин-слой тон кий (б г0), поле спадает в глубь проводника по закону ехр (—ж/б), а поток энергии в проводник выражается простой формулой [8]
|
сН1 |
1 м у/2 |
I |
1 |
(7.12) |
|
16л |
\2ла J ~ 32я2 |
аб ’ |
||
|
|
||||
‘S'oem/cjvt2 = |
3,16-10 |
2 (/г,Кав/см) |
Мгц/^ом-'см-1) 12 |
||
(мощность на единицу длины разряда |
W = 2 |
nr0S0) |
|||
Заменяя в (7.10) |
Н 0 через S0, получим |
|
|
||
тк |
|
|
|
(7.13) |
|
J |
c ( T ) X ( T ) d T ^ ± c KbKS0. |
Допустим, что температура Тк и степень ионизации невелики, так что а — Ne ~ ехр (— 1/2кТ). Вычисляя интеграл (7.13), так как это сделано в подразделе 24.4 при выводе формулы (6.60), найдем
4 ^ „ 2 кТ\ |
(7.14) |
|
что можно интерпретировать так же, как и формулу (6.60). Джоулево тепло выделяется главным образом в слое, где проводимость не более чем в е раз меньше максимальной величины ок на оси. Перепад температуры в этом слое А Т дается формулой (6.56). Толщина слоя, где поле диссипируется, порядка бкСледователь но, радиальный поток тепла, который выносит к трубке энергию
239
S0, поступающую в разряд от соленоида, |
порядка ХКАТ/6К. При |
|
равнивая эту величину S0, и получим (7.14) с точностью до чис |
||
ленного коэффициента. |
|
Питаев- |
Как было показано в работе Б. Э. Мейеровича и Л. П. |
||
ского [12J, в рассматриваемом случае 6К |
г0 и кТк<^_ I, |
распре |
деление температуры и поля в переходном слое от непроводящего газа к плазме, т. е. в области границы разряда, имеют универсаль ный характер, например Т = ТК/ (г/6к). В этой работе выведено и численно проинтегрировано уравнение для безразмерной темпе ратуры и построены профили всех величин в переходном слое. Найдена также точная связь температуры плазмы Тк с потоком электромагнитной энергии S0 в разряд. Эта связь, естественно, дается той же формулой (7.14), но содержит точное значение численного коэффициента— 3,14 вместо 4, что характеризует сте пень приближения модели металлического цилиндра в данном случае.
В заключение этого раздела мы хотели бы еще раз подчеркнуть, что температура разрядной плазмы определяется условием балан са потоков энергии в самой зоне диссипации поля и без рассмотре ния этого процесса, скажем, если в буквальном смысле не выходидить за рамки модели «металлического цилиндра», определить проводимость эквивалентного разряду проводника в принципе невозможно. Что же касается радиуса разряда г„, то его легко свя зать с выделяющейся в проводнике мощностью W условием тепло
отвода этой мощности через «непроводящий зазор». |
Интегрируя |
уравнение теплопроводности в «зазоре», получим |
приближенно |
^ X d T ^ Q K- Q CT^ ^ l n ~ . |
(7.15) |
^ст
28,2. Влияние частоты и порог режима. Вследствие резкого характера зависимости о (Т) даже для небольшого повышения температуры требуется значительное увеличение тока в индукто ре, что сопровождается существенным увеличением мощности, вкладываемой в разряд. Сказанное следует из формул (7.10), (7.13) и усугубляется появлением при высоких температурах потерь на излучение, которые не были учтены. Практически в ин дукционных разрядах атмосферного давления температуры не поднимаются выше примерно 10 000°. Примечательно, что ток в индукторе, необходимый для поддержания данной температуры, по формуле (7.10) не зависит от частоты. Этот вывод, однако, спра ведлив только в случае достаточно высоких частот. На низких частотах, так же как и при небольшой проводимости (температуре), толщина скин-слоя перестает быть малой по сравнению с радиусом разряда, и формулы (7.10), (7.13) теряют силу.
В работе Р. Е. Ровинского и А. П. Соболева [13] был рассмот рен практически важный вопрос: на какой частоте вообще выгод нее всего работать, если, скажем, иметь в виду получение плазмы
240