Электрический взрыв проводников (именуемый в дальнейшем ЭВП) традиционно является одним из способов получения высокотемпературной плазмы металлов. Это объясняется принципиальной простотой установок для ЭВП, представляющих собой обычно контур, содержащий накопитель энергии, коммутатор, токоподводящие цепи и собственно взрываемый проводник. Поскольку речь идет о плазме металлов, в дальнейшем подразумевается, что взрыв происходит в среде, не содержащей значительных количеств газа (вакууме).
Однако реализация процесса вклада энергии в проводник пропусканием через него электрического тока наталкивается на ряд трудностей. Одна из них связана с геометрическими особенностями системы типа цилиндрического проводника, взаимодействующего с собственным магнитным полем. В области микросекундных времен, как показано в работах [1−2], электрический взрыв, как правило, связан с магнитно-гидродинамическими (МГД) неустойчивостями, приводящими к распаду проводника на жидкометаллические фрагменты, разделенные плазменными областями. В этих областях перенос заряда осуществляется электронами, энергия которых существенно выше тепловой. Энергия электронов передается жидкому металлу и испаряет его. На финальной стадии процесса формируется сильно неоднородная область, состоящая из смеси плазмы, паров металла и неиспарившегося материала проводника. Вследствие шунтирования канала дуговым разрядом, развивающимся в наружных слоях, дальнейший вклад энергии в плазму становится неэффективным.
Развитие техники эксперимента позволило создать импульсные установки с высокой скоростью ввода энергии в проводник. Известны критерии, позволяющие оценить условия, при которых магнитно-гидродинамические неустойчивости при электрическом взрыве проводников не успевают развиться до испарения проводника. При этом для типичных значений диаметров проводников порядка десятых долей миллиметра, время, за которое в проводник вкладывается энергия, близкая к энергии сублимации должно быть менее 1 мкс. При быстром вводе энергии в проводник МГД — неустойчивости перестают определять характер процесса [3] в стадии, предшествующей испарению.
Однако, быстрый вклад энергии в одиночный проводник, осуществляемый, как правило, при помощи высоковольтной формирующей линии, приводит к образованию вокруг плотного ядра сравнительно малоплотной, хорошо проводящей «шубы», которая резко снижает сопротивление плазменного канала, уменьшая эффективность передачи энергии в приосевую область, занятую плотной плазмой. В результате этого, даже применяя источники энергии тераваттного диапазона мощности, при взрыве одиночного проводника удается разогреть шубу до температур не более 500эВ [4]. Образующиеся в приосевой области плотные высокотемпературные образования («плазменные точки») играют не очень значительную роль в энергобалансе плазмы. Вкладываемая в них энергия не превышает нескольких процентов от общего энергосодержания плазмы [11].
Дальнейшим развитием этого направления является использование в качестве нагрузки источника энергии цилиндрических оболочек (лайнеров) или многопроволочных массивов[5]. Результаты, достигнутые этим методом являются весьма высокими (температура 1−2 кэВ, электронная концентрация 1019-Ю20 см"), коэффициент передачи энергии в нагрузку достаточно высок благодаря тому, что на финальной стадии процесса сжатия происходит трансформация кинетической энергии материала лайнера в тепловую. Однако, установки, использующие принцип сжатия массива тонкопроволочных лайнеров являются чрезвычайно дорогими и сложными в силу необходимости использования сверхвысоких напряжений.
В связи с вышеизложенным представляется актуальным поиск альтернативных методов получения горячей плазмы материалов с высокими атомными номерами, применение которой является весьма перспективным с точки зрения технологических и научных применений (рентгеновская литография, спектроскопия, плазменная обработка материалов, создание активных сред для лазеров коротковолнового диапазона).
Таким методом может быть применение сильного магнитного поля, направленного вдоль оси проводника. В ряде работ продольное поле применялось как средство подавления магнитно-гидродинамических неустойчивостей. В ходе этих исследований, выполнявшихся с СПбГТУ [6], а также в численных экспериментах [7], была отмечена возможность нагрева плазмы до высоких температур за счет токов, индуцируемых в среде при ее расширении. Результаты предварительных оценок и расчетов дали основания предположить, что наряду с известными методами получения плотной горячей металлической плазмы может быть использован альтернативный метод, использующий взрыв проволочки в сильном магнитном поле. Успехи технологии получения сильных магнитных полей, в частности, создание ряда установок для получения импульса поля с индукцией 100 Т, позволяют считать, что эксперименты в полях такого масштаба являются вполне достижимыми в настоящее время. 7
Целью данной работы является экспериментальное подтверждение эффекта образования в периферийной области проводника, взрываемого в сильном продольном магнитном поле, плотной горячей металлической плазмы, нагреваемой азимутальным током, и возможности использования электрического взрыва в продольном поле проводников для создания источников плазмы с высоким энергосодержанием.
1. Литературный обзор.
Выводы. Экспериментально подтверждено образование вокруг канала плазмы электрического взрыва в продольном магнитном поле облака хорошо проводящей плазмы, в котором, судя по данным диамагнитных измерений протекает азимутальный ток с высокой плотностью. Поскольку концентрация частиц периферической плазмы значительно меньше, чем в «керне», можно ожидать, что периферия канала взрыва имеет весьма высокую температуру.
5. Численное моделирование периферийного плазменного слоя при электрическом взрыве проводника в сильном продольном магнитном поле с целью интерпретации экспериментальных результатов.
5.1. Постановка задачи.
Как было показано в гл. 2, из решения модельных задач о движении проводящей среды в магнитном поле следует возможность существования режимов, при которых вместо электромагнитного торможения проводника имеет место его ускорение. Результаты расчетов, приведенные в работах [7,59] предсказывают возможность существования режимов нестационарного течения, когда имеет место не только ускорение, но и нагрев плазмы ЭВП в магнитном поле порядка 100−120Т. В то же время, из экспериментальных данных следует, что при электрическом взрыве проводника в продольном магнитном поле значительно меньшей амплитуды (50−70Т) возникает быстро расширяющееся светящееся облако плазмы малой плотности, механизм разогрева которого вероятно аналогичен исследованному авторами расчетных работ, а условия развития определяются спецификой пограничного слоя взрываемого проводника.
Учитывая сложность более или менее полного расчета процесса электрического взрыва, для интерпретации экспериментальных данных имеет смысл выделить в отдельную задачу расчет расширения мало плотного плазменного облака в продольном магнитном поле, используя в качестве начальных условий варьируемые параметры периферической плазмы. Практически, при варьировании начальных параметров плазмы делается попытка решить обратную задачу определения их реальных значений путем сопоставления результатов расчетов и экспериментов.
Другой целью численных расчетов является качественный анализ влияния конкретного вида зависимости электропроводности плазмы от ее термодинамических параметров на характер процесса с точки зрения возможности дополнительного разгона и разогрева наружных слоев плазмы по сравнению с адиабатическим расширением. Несмотря на качественный характер расчетов, их результаты могут быть использованы как для интерпретации экспериментальных данных, так и для выбора условий дальнейших экспериментов с целью достижения высокой скорости и температуры плазмы при электрическом взрыве в сильном магнитном поле.
5.2 Описание расчетной модели.
Для численного моделирования периферийного плазменного слоя с целью интерпретации экспериментальных результатов была применена модифицированная компьютерная программа, использовавшаяся для моделирования течения плазмы с постоянной проводимостью. Основным принципом при разработке расчетной модели было ее максимальное упрощение и возможность использования параметров, надежно измеренных экспериментально. К ним могут быть отнесены результаты, полученные из сопоставления результатов электроннооптического фото-хронографирования, измерений диамагнитного сигнала, а также осциллограмм тока и напряжения на проводнике.
Внутренняя область взрываемого проводника представляет собой плазму с плотностью порядка 1026 м~3. По литературным данным известно, что температура плазмы на поверхности столба при взрыве проводника в режимах, близких к использовавшемуся в экспериментах при отсутствии магнитного поля) имеет порядок 3−5 эВ, а
94 «У Я 3 концентрация частиц в потоке плазмы, вытекающем из ядра — 10 -10 м~
VI
Поскольку существующие модели электрического взрыва не описывают процесс ионизации пара на поверхности проводника, предполагалось, что кратность ионизации плазмы Z, поступающей из ядра проводника соответствует равновесной для заданной начальной температуры [55].
В отличие от программы, применявшейся для решения модельных задач с постоянной электропроводностью, в данной серии расчетов модель была несколько уточнена за счет использования уравнения состояния плазмы более адекватно описывающего процесс расширения по сравнению с моделью идеального газа, учета протекания по периферийной плазме аксиального тока и наличия азимутального магнитного поля.
Упрощенный характер расчетной модели, в которой рассчитывались только параметры внешней мало плотной области плазмы, определил выбор вида уравнения состояния плазмы. В настоящее время существует большое количество широкодиапазонных уравнений состояния для металлов, позволяющих достаточно хорошо описывать их поведение при изменении плотности от твердотельной до величин, характерных для идеального газа [75,76]. Однако в данной работе, в связи с тем, что центральная плотная область плазмы не рассматривалась, плотность плазмы ограничивалась величиной порядка одной десятой твердотельной. Как показывают оценки, в этой области параметров достаточно хорошим приближением является уравнение Ван дер Ваальса в виде N0 а
Р ~ у ~у 2 ,[77] где N — число частиц в объеме V, 0 температура. Параметры, а и Ь определялись по параметрам критической точки для меди по данным работы [64]. Уравнение Ван дер Ваальса достаточно хорошо описывает давление плазмы при начальных параметрах, лежащих выше критической точки.
Величина аксиального тока в плазме задавалась в соответствии с данными эксперимента в предположении, что в шубе протекает половина суммарного тока в цепи.
5.3 Выбор модели электропроводности плазмы.
Известные соотношения для проводимости, полученные при допущении о термодинамическом равновесии среды, перестают быть справедливыми в области малой плотности, где, согласно расчетам, концентрация частиц имеет порядок 1021 м" 3. Действительно, поскольку плотность азимутального тока имеет порядок 108 А/м2, то при указанной концентрации электронов их дрейфовая скорость имеет порядок 10б м/с, что существенно выше тепловой скорости ионов с температурой Т^1кэВ. При таких условиях предположение об отсутствии термодинамического равновесия и снижении электропроводности вследствие раскачки хаотических колебаний в плазме кажется вполне обоснованным [48]. Поэтому расчет параметров плазмы мог проводиться как с использованием нескольких известных моделей проводимости, справедливых в предположении термодинамического равновесия, так и при применении в расчете эмпирических моделей, параметры которых подбираются из условия минимального отличия результатов расчетов и экспериментов.
В качестве первого приближения расчет могла закладываться проводимость, вычисленная по модели Спитцера: х = 0.61 • 109 Г3/2 / X ! 1п Л, где Т — температура плазмы, измеренная в килоэлектронвольтах, Ъ — кратность ионизации, 1пЛ — кулоновский логарифм [44].
Результаты расчетов разлета плазмы с начальной плотностью 60 л кг/м, температурой 5 эВ, скоростью истечения 3 км/с в магнитном поле В2=50Т приведены на рис. 46.
Расчет с использованием указанной модели дает быстрое торможение границы плазмы, не наблюдающееся в экспериментах. Величина расчетного диамагнитного сигнала значительно меньше наблюдаемой в экспериментах.
На следующем этапе расчетов для описания проводимости плазмы была использована аппроксимация модели Беспалова — Полищука [78], которая хорошо описывает поведение проводимости термодинамически-равновесной плазмы в широком диапазоне плотностей и температур. Не останавливаясь на физическом обосновании модели, приведем ход зависимости проводимости от плотности и температуры (рис. 47). Как видно из приведенной диаграммы, при температурах масштаба 1 эВ и плотностях, близких к твердотельной, проводимость плазмы резко падает с уменьшением ее плотности и уменьшается с ростом температуры. При высоких температурах и малых плотностях, характерных для наружных слоев расширяющейся плазмы, поведение проводимости плазмы приближается к спитцеровскому закону, т. е. проводимость мало зависит от плотности и резко растет с ростом температуры. Из приведенных на рисунке линий равной проводиости видно, что в области температур 10−100 эВ электропроводность плазмы еще существенно зависит от плотности.
V м/с
Т эВ
16 000 14 000 12 000 10 000 8000 6000
1 25
2 / 5 0
0 1Е-07 2Е-07 ЗЕ-07 t с
И. м
ДФ Вб
0.001 б)
1 с
Рис. 46. Результаты расчета процесса разлета плазмы с начальной плотностью 60 кг/м3 в магнитном поле В=50 Т с использованием модели проводимости Спитцера. Кривые 1 и 2 (а) — скорость границы и температура, 3 и 4 (б) — радиус плазмы и диамагнитный сигнал. 5 — радиус плазмы по данным хронографирования.
Ом-м)
ТК Р кг/м3
Рис. 47. Зависимость проводимости медной плазмы от плотности и температуры, рассчитанная по модели [78]
Результаты расчетов с применением указанной модели весьма резко расходятся с экспериментальными данными, предсказывая очень быстрый разгон с последующим торможением плазмы за время порядка десятков наносекунд (рис. 48). Была проделана серия расчетов с введением поправочного множителя к исходной модели. Как видно из рис. 48, удовлетворительно описывает эксперимент вариант расчета с проводимостью, сниженной на порядок по сравнению с моделью Беспалова — Полищука. Величина расчетного диамагнитного сигнала при этом по порядку величины близка к измеренной.
Качественным отличием результатов, полученных с использованием модели Беспалова — Полищука, которая учитывает зависимость проводимости плазмы от ее плотности, является значительный дополнительный разгон плазмы в магнитном поле по сравнению с адиабатическим расширением. Физический смысл этого явления заключается в том, что снижение плотности наружных слоев плазмы в процессе расширения ведет к уменьшению их проводимости и увеличению омического нагрева более глубоких областей плазмы. Рост газокинетического давления в этих слоях приводит к дополнительному разгону периферии, что, вероятно, и наблюдается в эксперименте. В пользу такого механизма дополнительного ускорения границы плазмы в модели с электропроводностью, падающей с уменьшением плотности, говорит сравнение профилей давления, температуры и плотности тока, полученные в расчетах. На рис. 49,50 показаны профили, полученные в вариантах расчетов с постоянной проводимостью и с использованием модели Беспалова-Полищука (с поправочным множителем 0.1). Как видно из графика, во втором случае имеется характерная область повышенного давления и плотности тока, которая по мере расширения плазмы перемещается в относительно более глубокие слои плазмы, способствуя разгону наружной области.
V м/с Т эВ к м АФ Вб с t с а) б)
Рис. 48. Результаты расчета процесса разлета плазмы в магнитном поле В=50 Т с использованием апроксимации модели проводимости Беспалова-Полищука. На рис. (а) — кривые 1,2 -скорость границы, 3,4 -температура на границе относятся к соответственно к величине поправочного множителя 0.5 и 0.1. На рис. (б) для тех же значений поправочного множителя — кривые 1,2 -внешний радиус плазмы, 3,4 -диамагнитный сигнал.
Р Н/м
Ум/с
1.0Е + 09 8.0Е+08 6.0Е+08 4.0Е+08 2.0Е+08 О. ОЕ + ОО
4Ь 4 с
4а /
1 / /
Ж' Vx3b Зс
За —
0.0012
0.0016 И м
15 000 12 000 9000 6000
3000 0.002
РН/м2 1.0Е+09
8.0Е+08
6.0Е+08
4.0Е+08
2.0Е+08
0.0Е+00
4Ь
4с
4а
V ,
Л^Д д
0.0012
0.0016
И м
Ум/с 18 000
0.002
Рис. 50. Профили давления и скорости для расчета с использованием модели постоянной проводимости а=5000 (Ом-м)" 1 (вверху) и апроксимации модели Беспалова-Полищука с поправочным множителем 0.1 (внизу). 3 -давление, 4- скорость
Индексы а, Ь и с относятся к моментам времени 20, 40 и 60 не.
Поскольку расчеты с использование моделей, построенных для термодинамически-равновесной плазмы, неточно описывают процесс разлета плазмы в магнитном поле, имеет смысл произвести подбор некоторой постоянной величины проводимости, дающей наилучшее совпадение расчета и эксперимента. юооо
К м
0.001
2Е-07 4Е-07 6Е-07
2Е-07 4Е-07
ДФ Вб 0.0Е+00
— 2.0Е-04
6Е-07 t с t с
Рис. 51. Результаты расчета процесса разлета плазмы в магнитном поле В=50 Т с использованием постоянной проводимости о=1000 (Ом-м)-1 и а=5000 (Ом-м)-1. На рис. (а) — кривые 1,2 — скорость границы, 3,4 -температура на границе относятся к соответственно к величине проводимости а=1000 (Ом-м)-1 и а=5000 (Ом-м)-1. На рис. (б) для тех же значений проводимости — кривые 1,2 -внешний радиус плазмы, 3,4 -диамагнитный сигнал.
Одним из главных критериев совпадения является величина диамагнитного сигнала к моменту достижения плазмой стенок камеры. Результаты расчета для величины проводимости а=1 ООО (Ом-м)" 1 и а=5000 (Ом м)" 1 приведены на рис. 51. Из графика видно, что основные характеристики процесса удовлетворительно описываются расчетом с постоянной проводимостью а=5000 (Ом-м)" 1. Следует заметить, что эта величина проводимости по порядку величины совпадает с проводимостью, полученной из осциллограмм напряжения. При этом, однако, значение скорости границы плазмы ниже измеренного экспериментально, что может говорить о том, что в реальном процессе происходит снижение проводимости плазмы на периферии, приводящее к дополнительному ускорению границы.
Возможность создания мощного источника импульсного источника рентгеновскогоизлучения на основе ЭВП в сильном магнитном поле.
Проведенные эксперименты не дают однозначного ответа на вопрос о перспективности использования ЭВП в продольном магнитном поле, как источника рентгеновского излучения. Используя расчетную модель с постоянной проводимостью плазмы можно проиллюстрировать изменение суммарного выхода излучения из расширяющегося плазменного столба с заданной начальной температурой в режимах, характеризуемых наличием и отсутствием продольного магнитного поля. С этой целью расчетная модель была модифицирована добавлением в уравнение энергобаланся члена, учитывающего потери энергии плазмы на излучение. Поскольку полный учет влияния излучения на энергобаланс является весьма сложной проблемой, для сравнительных оценок использовались литературные данные о мощности излучения из объема плазмы, расчитанные для режима коронального равновесия [55]. Предполагалось, что плазма является оптически прозрачной. Результаты расчета зависимости полной энергии излучения на единицу длины плазменного столба от его начальной температуры при наличии продольного поля В7=50Т и в случае его отсутствия для начальной плотности плазмы 60 кг/м3 и проводимости а=5000(0м-м)~1 приведены на рис52. Как видно из результатов расчета, наличие продольного магнитного поля увеличивает суммарную энергию излучения по сравнению с режимом В2=0, что может быть объяснено дополнительным разогревом периферийной плазмы при расширении поперек магнитного поля, а также уменьшением скорости адиабатического охлаждения плазмы за счет тормозящего действия поля.
О Дж/см
70 -60 -50 -40
30 -20 -10
0 -,
О 50 100 150 200
ТэВ
Рис. 52. Зависимость полной энергии излучения от начальной температуры плазменного столба при наличии продольного поля В2=50Т (1) и в случае его отсутствия (2).
Учет потерь на излучение в суммарном энергобалансе плазмы приводит к снижению расчетной температуры поверхности плазмы в несколько раз по сранению с вариантом без учета излучения. Тем не менее, и в этом случае температура плазмы достаточно высока для того, чтобы являться излучателем мягкого рентгеновского излучения (температура поверхностных слоев плазмы к момету перехода в режим течения с постоянной скоростью на порядок превышает начальную).
Заключение
.
Проведенные экспериментальные и расчетные исследования вакуумного ЭВП в сильном продольном магнитном поле подтвердили возможность получения горячей плазмы на периферии канала разряда за счет генерации в ней азимутального тока при расширении поперек силовых линий магнитного поля.
Экспериментально определенная скорость расширения канала плазмы в исследованных режимах достигает 20 — 30 км/с. Измеренная величина диамагнитного сигнала соответствует линейной плотности азимутального тока 106А/м.
Проведенные численные расчеты на основе данных эксперимента позволяют оценить плотность тока и температуру плазмы в одножидкостном приближении на периферии плазмы. По данным расчетов плотность азимутального тока достигает 109А/м2, что примерно на порядок превышает плотность аксиального тока. При этом расчетная температура периферийных слоев плазмы составляет несколько десятков электрон-вольт (без учета излучения).
На основе сравнения результатов эксперимента с расчетами, проведенными с использованием различных моделей проводимости плазмы, показано, что наблюдаемое в экспериментах ускорение внешней границы плазмы может быть объяснено смещением локального максимума плотности тока внутрь плазмы за счет уменьшения проводимости внешних слоев.
В заключение перечислим основные результаты, полученные автором в данной работе: 1. Впервые проведено экспериментальное исследование электрического взрыва проводника в вакууме в продольном магнитном поле с амплитудой 50 Т и выше. Отработана технология проведения экспериментов по электрическому взрыву в сильном магнитном поле, обеспечивающая получение воспроизводимых результатов.
2. Впервые экспериментально показано, что в условия проведенных экспериментов наблюдается переход наружных слоев плазмы канала электрического взрыва в режим разогрева и ускорения за счет взаимодействия с магнитным полем. При этом граница плазмы движется со скоростью порядка 20 км/с.
3. Впервые применена техника диамагнитных измерений для диагностики плазмы электрического взрыва, что в сочетании с высокоскоростным фотоэлектронным хронографированием позволило измерить азимутальный ток в плазменном облаке.
4. На основании полученных экспериментальных данных о параметрах плазмы была проведена серия численных расчетов, продемонстрировавших возможность ускорения плазмы магнитным полем при возникновении аномального электрического сопротивления в областях с высокой плотностью тока.
5. Показано, что азимутальный ток, генерируемый в расширяющейся плазме, перераспределяет энергию плазмы таким образом, что наружные слои получают основную ее часть, нагреваясь до существенно более высоких температур, чем при электрическом взрыве без продольного поля. Этот факт может быть использован при разработке источника высокотемпературной плазмы для источников рентгеновского излучения и плазменных размыкателей.
Основное содержание диссертации опубликовано в следующих работах:
1. Адамьян Ю. Э. Березкин А.Н. Шнеерсон Г. А. VI Электрический взрыв проводника в магнитном поле. Всесоюзная конференция по физике низкотемпературной плазмы. Тезисы докладов, т.2, стр.387−389, Ленинград, 1983 г.
2. Адамьян Ю. Э., Березкин А. Н., Бочаров Ю. Н. и др. Стенд для создания сильных импульсных магнитных полей. ПТЭ, 1994, вып.5, с. 134−139.
3. Адамьян Ю. Э., Василевский В. М. Кривошеев С.И., Колгатин С. Н.,
Шнеерсон Г. А. Шутов В.Л. Электрический взрыв проводников в сильном магнитном поле. Письма в ЖТФ, 1995, N. 23, с. 43 -46.
4. Адамьян Ю. Э., Березкин А. Н., Василевский В. М. и др. Процесс расширения и нагрева плазмы при электрическом взрыве проволочки в сверхсильном магнитном поле. Тезисы докладов 6 Международной конференции по генерации мегагауссных полей и родственным экспериментам, г. Саров, 5−10 августа 1996.
5. Yu.E. Adamyan, Y.M. Vasilevsky, S.N. Kolgatin, G.A.Shneerson, Forming of high speed radial plasma flow due to its heating by induced azimuthal current in electrical explosion of wires in extra high magnetic field, Proceedings of 12th international conference on high-power particle beams, Haifa, Israel, 1998.
6. Адамьян Ю. Э., Василевский B.M., Колгатин C. H, Шнеерсон Г. А.
Эффект нагрева поверхностного слоя плазмы при взрыве проводника в вакууме в сильном продольном магнитном поле. ЖТФ, 1999, вып 5, с. 121−127