Научная статья на тему 'Моделирование катодонаправ-ленного стримера в неоднородном электрическом поле'

Моделирование катодонаправ-ленного стримера в неоднородном электрическом поле Текст научной статьи по специальности «Физика»

CC BY
317
90
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ГАЗОВЫЙ РАЗРЯД / ПОЛОЖИТЕЛЬНЫЙ СТРИМЕР / ВОЗДУХ / ДРЕЙФОВО-ДИФФУЗИОННОЕ ПРИБЛИЖЕНИЕ / GAS DISCHARGE / POSITIVE STREAMER / AIR / DRIFT-DIFFUSION APPROXIMATION

Аннотация научной статьи по физике, автор научной работы — Пек Борис Эдуардович, Самусенко Андрей Викторович, Стишков Юрий Константинович

В статье рассмотрено решение задачи о распространении стримера в неоднородном электрическом поле в дрейфово-диффузионном приближении. В качестве газа, в котором развивается процесс, выступает смесь кислорода и азота, моделирующая сухой воздух. Рассматривается три типа частиц: электроны, положительные и отрицательные ионы. Система уравнений включает транспортные уравнения для каждого типа частиц, уравнение Пуассона для электрического потенциала. Особенностью работы является анализ влияния мощности и типа разных источников ионизации на динамику стримера. В качестве источников ионизации рассматривается источник, создающий однородную концентрацию электронов в пространстве, а также фотоизлучение головки стримера. В последнем случае система уравнений дополняется уравнением на концентрацию фотонов. Выявлен начальный участок движения стримера, осуществляемый при помощи электронов, оставшихся от лавинного этапа разряда. Для этого участка характерен рост скорости, а затем спад при переходе стримера через максимум концентрации электронов в лавине. Показано, что при отсутствии фотоионизации либо внешнего источника ионизующего облучения, на этом рост стримера прекращается. Получены данные о влиянии мощности ионизующего излучения, испускаемого головкой стримера, на скорость стримера, радиус его головки, напряжённость поля на переднем фронте стримера. Показано, что увеличение коэффициента поглощения ионизующего излучения головки воздухом затрудняет развитие стримера. Показано, что качественная картина развития стримера не зависит от рассмотренных типов источников излучения и слабо зависит от интенсивности излучения. Библиогр. 10 назв. Ил. 7. Табл. 1.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Simulation of a cathode-directed streamer in non-uniform electric field

Solution of the problem of streamer development in a non-uniform electric field in the drift-diffusion approximation is considered. Models with different types of an ionization source are compared. The data of ionization power and its source type influence on streamer velocity, its head radius, electric field intensity are analyzed.

Текст научной работы на тему «Моделирование катодонаправ-ленного стримера в неоднородном электрическом поле»

Б. Э. Пек, А. В. Самусенко, Ю. К. Стишков

МОДЕЛИРОВАНИЕ КАТОДОНАПРАВЛЕННОГО СТРИМЕРА В НЕОДНОРОДНОМ ЭЛЕКТРИЧЕСКОМ ПОЛЕ

Введение. Основной метод, используемый для компьютерного моделирования развития стримера в газах, - решение системы уравнений в частных производных для концентраций частиц и электрического поля в осесимметричном приближении [1, 2]. В [3] этот метод был описан для моделирования развития электронных лавин в неоднородном электрическом поле; были выявлены эффекты значительного отклонения профиля лавины от функции Гаусса вплоть до появления локального максимума концентрации электронов у поверхности активного электрода. При увеличении напряжения в той же системе становится возможен следующий после лавины этап разряда - лавинно-стри-мерный переход [4] с последующим развитием стримера. Моделирование этих явлений описано в данной работе.

Использована упрощённая модель воздуха, включающая лишь основные реакции, что позволило выявить влияние отдельных механизмов ионизации на структуру и динамику развития стримера. В отличие от [3] кинетические эффекты не учитываются, применяется дрейфово-диффузионное приближение, поскольку оно отражает все основные закономерности развития лавин и стримеров, приведённые в [2]. Система уравнений в частных производных решалась в программном пакете Сот8о1 методом конечных элементов.

Постановка задачи. Система уравнений состоит из транспортных уравнений для положительных и отрицательных ионов и электронов и уравнения Пуассона:

дпе

т ‘

дп+

~дГ

дп_

~дГ

ё1у - ё1у

—D еУПе —

^;оп (Е) пе (Е) пе Сеепеп+ + ^рЬ + Я

—Б+Уп+ + ц+п+Е = У;оп (Е) Пе — СееПеп+ — Сцп+п- + й'рЬ + 5ип,

—О—'Чп— — ц-п-Е = уай (Е) пе — Сцп+п-,

(1)

Дф =------(п_|_ — п_ — пе); Ё = —Уф.

е0

Здесь пе, п+ и п— - концентрации электронов, положительных и отрицательных ионов соответственно, ф - электрический потенциал. Функция источника для электронов учитывает рождение за счёт ударной ионизации с частотой У;оп (Е), гибель при прилипании с частотой уай (Е), гибель при электрон-ионной рекомбинации с коэффициентом сее. Функция источника для положительных ионов учитывает рождение за счёт ударной ионизации и гибель при электрон-ионной и ион-ионной рекомбинации (для последней задан коэффициент с;;).

Функция источника для отрицательных ионов учитывает их рождение за счёт прилипания электронов к нейтральным частицам и гибель при ион-ионной рекомбинации. Зависимости частоты ионизации и прилипания от напряжённости поля определены путём аппроксимации экспериментально полученных кривых, приведённых в [2, 5-8].

© Б. Э. Пек, А. В. Самусенко, Ю. К. Стишков, 2010

Коэффициент электрон-ионной рекомбинации задан постоянным:

10~13 м3/с, так как он слабо меняется в исследуемом диапазоне напряжённости электрического поля [5].

Для моделирования выбрана одна из простейших систем электродов с неоднородным полем: сфера-сфера. Радиус внутренней сферы - 1 мм, внешней - 2,01 мм. Поле слабонеоднородное: отношение максимального поля к среднему в отсутствие объёмного заряда равно двум.

Расчётная область разделена на две. В области А (рис. 1) решаются все уравнения системы (1), а в области В - только уравнение Пуассона и уравнение на концентрацию фотонов. Разделение на две области позволяет сэкономить расчётные ресурсы. Оно возможно, поскольку заряженные частицы (электроны и ионы) присутствуют практически только в области А, но объёмный заряд возмущает электрическое поле на расстоянии порядка нескольких радиусов стримера.

Граничные условия, обозначенные на рис. 1, представлены в таблице.

Граничные условия

р Электрический потенциал ф Концентрации электронов пе, ИОНОВ П+ И П— Концентрация фотонов прь

1 Ф = 0 — + |х+п+Е = 0 прь = 0

2 Ф = 0 - прь = 0

3 Ш=Ег,=0 - =0

4 Ф = —10 кВ - прь = о

5 - — В+ Уп+ + |х+п+Е = 0 — Ве\7пе — |хепеЕ = 0 — В-Уп— — |х-П—Ё = 0 -

6 Ф = —10 кВ — Ве\7пе — |хепеЕ = 0 — В-Уп— — |х-П—Ё = 0 ПрЬ = 0

7 ЄФ = 0 оп о^=0 £ЇН=0; ^ = 0 ОП ’ оп ’ оп Оп ~ и

В качестве начальных условий (НУ) вблизи катода задано начальное облако электронов: начальная концентрация электронов в облаке задана на уровне 1016 1/м3, в остальном пространстве начальная концентрация электронов задана равной нулю. Концентрации других частиц в начальный момент времени предполагаются равными нулю. Напряжение равно 10 кВ и неизменно во времени.

Для выявления влияния разных источников ионизации на структуру и динамику стримера рассмотрены три качественно разные модели: без фотоионизации и внешнего

источника (5рь = 0, 5ЦП = 0), с внешним источником ионизации (й'рь = 0, 5^ = 0)

и фотоионизацией за счёт излучения стримера (5рь = 0, 5ип = 0). Далее фотоионизация и ионизация внешним источником называются для краткости дополнительными источниками ионизации (в противовес ударной ионизации, основному источнику).

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

а ___ ^шах

^ип - 7

Т

где птах - максимальная фоновая концентрация свободных электронов, которую может обеспечить внешний источник ионизации, Т - характерное время, за которое источник выводит концентрацию на фоновый уровень. В приведённых результатах задавались значения т = 0,01 нс; птах = 107 м~3. Такой источник обеспечивает постоянную однородную концентрацию электронов в области вокруг стримера.

Для учёта фотоионизации излучением стримера в систему уравнений добавлялось стационарное уравнение диффузии для фотонов высоких энергий (около потенциала ионизации молекулы кислорода и выше), способных привести к фотоионизации молекул воздуха:

АпрЬ = 1 пръ + 9"^'юп ) ne, (2)

где I - коэффициент поглощения ионизующего излучения в воздухе; пр^ - концентрация фотонов. Безразмерный коэффициент д (коэффициент мощности фотоионизации) характеризует отношение частоты производства ионизующих фотонов к частоте ударной ионизации. В [10] указано, что он практически не зависит от напряжённости поля. Изменяя этот коэффициент, мы изменяем мощность ионизующего излучения головки и, таким образом, можем изменять влияние фотоионизации на развитие стримера.

Для учёта фотоионизации в функции источника для электронов и положительных ионов присутствовало слагаемое, описывающее появление ион-электронной пары при поглощении фотона 5рь:

$рь = clnph, (3)

где с - скорость света в вакууме.

Уравнение (2) является приближённым способом описания фотонов. Обоснование применимости такого метода для условий данной задачи (атмосферное давление, размеры расчётной области не более 1 мм) можно найти в статье [9].

Источник й'р^ (фотоионизация излучением головки) даёт неоднородную концентрацию затравочных электронов - при удалении от головки она падает. Источник 5цп (ионизация внешним источником), напротив, даёт однородную концентрацию затравочных электронов во всей расчётной области.

Анализ результатов.

Фаза лавины и лавинно-стримерный переход. Проанализируем основные этапы развития стримера. На рис. 2, 3 показано изменение концентрации электронов в процессе лавинно-стримерного перехода и дальнейшего прорастания встречного положительного стримера к катоду. Здесь и ниже приведены результаты для модели с фотоионизацией излучением стримера с коэффициентами д = 10~6 и I = 103 м_1. Облако электронов, заданное начальными условиями, движется от катода к аноду. Количество электронов в нём постепенно нарастает под действием ударной ионизации, а радуис облака увеличивается из-за диффузии. Объёмный заряд пока слишком мал, чтобы заметно

3,48 не ,,.Ь

3,82 не

\ 4,14 нс

\

і і к

4,50 не

4,80 не

5,06 не

і

0,8

0,6

0,4

0,2

0,2 мм

М

М

Рис. 2. Концентрация электронов в области в разные моменты времени; а - старый б - новый локальные максимумы концентрации в лавинной фазе:

концентрации нормированы на максимум, К - катод, А - анод, М - локальный максимум концентрации

Рис. 3. Концентрация электронов

вдоль оси симметрии в разные моменты времени; а - старый б - новый локальные максимумы концентрации в лавинной фазе:

М - локальный максимум концентрации

I

Анод

Длина, мм

і

Катод

изменить электрическое поле - это так называемая фаза лавины. Она подробно описана в [3], здесь мы не будем на ней останавливаться. В [3] также подробно описан эффект «втягивания» лавины в область слабонеоднородного поля. Этот эффект наблюдался и в названной модели: в определённый момент помимо основного («старого») максимума концентрации возникает второй («новый»), непосредственно на поверхности электрода (рис. 2, 3, £ = 3,48 нс). Появление нового максимума вызвано первыми немногочисленными электронами, достигшими анода. Из-за диффузии отдельные электроны достигают анода раньше, чем центральная часть лавины. Эти электроны из-за высокой напряжённости поля вблизи анода и, как следствие, высокой частоты ионизации приводят к быстрому повышению концентраций непосредственно у поверхности анода. Здесь образуется новый максимум, который быстро обгоняет старый по величине, и в момент £ = 3,82 нс старый максимум становится незаметным на его фоне. Тем не менее, поскольку скорость ионизационных процессов значительно выше скорости

а

0,2 мм

4,14 нс +1

+5

кВ/см 300 200 100 0

Рис. 4- Контурные графики интенсивности ионизации У}оп • пе (а), распределения объёмного заряда (б) и напряжённости электрического поля (в) в межэлектродной области в моменты времени 4,14 и 5,06 нс

в

дрейфа электронов, некоторая относительно малая концентрация электронов остаётся и в отдалении от анода.

Этапы развития стримера. Появлению стримера предшествует описанный выше этап, в конце которого лавина коснулась своей передней частью анода. Электроны быстро уходят на анод, и в прианодном слое постепенно накапливается некомпенсированный положительный заряд (рис. 4б), который через некоторое время начинает существенно влиять на внешнее электрическое поле, вытесняя его из области повышенной концентрации электронов на границу этой области (рис. 4в). Туда же смещается и область повышенной ионизации (рис. 4а). Это приводит к повышению концентрации ионизованных частиц перед границей, и область повышенной концентрации заряженных частиц вновь изменяет свои границы. От анода начинает прорастать положительный стример навстречу лавине - в этом смысле его можно назвать «отражённым» стримером. При этом область повышенной концентрации электронов переходит в состояние плазмы - поле оттуда вытеснено, объёмного заряда практически нет.

Это соответствует сложившимся представлениям о структуре сформировавшегося стримера, для которого характерно наличие плазменного канала, чьё поле вытеснено на внешнюю границу канала [4]. Максимум поля при этом оказывается на передней границе эллипсовидной головки канала (см. рис. 4в). Если лавина развивалась в поле, образованном электродами, практически не влияя на него, то с момента образования стримерного канала проводящая область прорастает за счёт вновь образующейся плазмы и максимум поля смещается от поверхности электрода на внешнюю границу головки стримера, ближе к противоэлектроду. Этот максимум становится новым центром ионизационных процессов, в которых в качестве затравочных электронов участвуют электроны начальной лавины. В результате ионизации возникает новая плазменная область с высокой концентрацией электронов, удлиняющая канал стримера, а максимум поля вновь смещается в направлении к противоэлектроду. Так движется волна ионизации.

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

расстоянии от головки. В качестве затравочных электронов могут выступать оставшиеся от первоначальной лавины и вновь образовавшиеся в результате фотоионизации газа излучением стримера или ионизации газа внешним источником.

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

В сформировавшемся виде стример имеет эллипсоподобную плазменную головку, окаймлённую тонким слоем объёмного заряда (см. рис. 4б), и цилиндрический «плазменный» канал. Электрическое поле в канале и головке стримера существенно понижено, а перед головкой повышенно, зона максимальной ионизации локализуется непосредственно перед головкой стримера (см. рис. 4а, і = 4,14 нс). Полученная в модели структура катодонаправленного стримера имеет свой традиционный вид (рис. 4б и 4в при 4,14 нс), описанный в литературе [2]: квазинейтральная плазменная область с пониженным электрическим полем в канале стримера и область повышенной напряжённости поля перед головкой стримера.

В головке стримера и на границе стримерного канала находится нескомпенсиро-ванный положительный электрический заряд, образованный положительными ионами. Такая структура с ярко выраженным стримерным каналом, замкнутым на электрод, и чётко ограниченной головкой характерна для положительных стримеров. Отрицательные же стримеры имеют диффузный вид, со слабовыраженной структурой головки [1].

В [2] обосновано введение параметров стримерного канала, которые можно назвать интегральными, поскольку они характеризуют его как цельную проводящую структуру. Они позволяют наглядно представить динамику развития стримера и используются и в других работах [1]. Мы определим интегральные параметры следующим образом:

а) длина стримера, которую мы будем понимать как расстояние от поверхности анода до максимума электрического поля перед головкой;

б) скорость роста стримера, которую мы определим как производную по времени от его длины;

в) эффективный радиус головки, который мы определим как удвоенное расстояние между максимумом напряжённости и точкой перед головкой, в которой поле падает в 2,25 раза по сравнению с максимальным.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Последнее определение связано с тем, что радиус головки влияет на динамику стримера прежде всего тем, что задаёт характерный масштаб падения поля перед головкой [2]. В упрощённой аналитической модели из того же источника предполагается, что поле перед головкой можно приближённо считать полем от сферического проводника: оно спадает по закону обратных квадратов. Легко убедиться, что использованное нами определение даёт точный результат - радиус сферы - для случая сферического заряженного проводника; для стримера же, головка которого имеет сложную форму, это будет лишь некоторый характерный пространственный масштаб головки.

Интегральные параметры позволяют описать влияние различных источников ионизации на динамику стримера. На рис. 5 приведены зависимости длины стримера от

Время, нс

т

Анод

Длина стримера, мм

т

Катод

<—► без фотоинозации

1--- с фотоионизацией (g = 10-6, l = 103 м-1)

2 -. с фотоионизацией (g = 10-6, l = 104 м-1)

3--- с фотоионизацией (g = 10-6, l = 105 м-1)

4 -- с фотоионизацией (g = 10-9, l = 103 м-1)

5 .... с внешним источником

(т = 10-11 с, n = 107 м-3)

max

0,2 0,4 0,6 0,8 1

Длина стримера, мм

Рис. 5. Интегральные параметры стримера - зависимости длины стримера от времени (а), скорости стримера (б) и диаметра головки стримера (в) от его длины в моделях:

без фотоионизации, с внешним источником ионизации и с фотоионизацией за счёт излучения стримера (4 примера)

в

времени, скорости стримера от его длины и диаметра головки стримера от его длины для разных механизмов ионизации.

В отсутствие фотоионизации после прохождения головкой стримера области расположения начального пакета электронов стример останавливается: в определённый момент волна ионизации практически прекращает движение, заряженный слой, окаймляющий стример, утончается; напряжённость поля растёт, пытаясь компенсировать нехватку затравочных электронов перед головкой увеличением скорости ионизации; задача расходится. В моделях, учитывающих фотоионизацию излучением стримера или внешнего источника, скорость стримера, его диаметр, вплоть до длины 0,5 мм, остаются теми же, что и в отсутствие этих источников. Тем самым показано, что на начальном этапе динамика стримера определяется не фотоионизацией, а в качестве затравочных электронов выступают остаточные электроны лавины. С этой позиции можно объяснить и первый максимум скорости на этом этапе развития стримера (рис. 5б): он соответствует максимуму концентрации электронов в основной части лавины. До достижения головкой стримера положения этого максимума скорость стримера растёт, а диаметр головки сужается (рис. 5в). Когда волна ионизации проходит через этот максимум, концентрация остаточных электронов быстро идёт на спад и стример притормаживается.

Этот начальный этап движения стримера можно назвать «несамостоятельным», поскольку затравочные электроны вокруг головки не образовались в результате влияния стримера, а остались от электронной лавины. Такой источник затравочных электронов является «невозобновляемым»: когда стример уйдёт из области, занятой электронной лавиной, затравочных электронов больше не будет.

В моделях, где учитывается собственная фотоионизация или внешний источник, в определённый момент времени скорость стримера после прохождения начального этапа вновь начинает расти. Динамика стримера зависит от механизма и параметров дополнительных источников ионизации. Этот этап можно назвать «самостоятельным», так как стример далее не зависит от электронов начальной лавины. В моменты времени 4,8 и 5,06 нс (см. рис. 2 и 3) виден локальный максимум (М) концентрации электронов, соответствующий области, в которой произошёл этот переход.

Момент перехода и область, где он происходит, зависит от интенсивности фотоионизации дополнительного источника и той концентрации электронов перед головкой, которую он может обеспечить. Чем больше заданный коэффициент мощности фотоионизации, тем раньше проявляется её влияние: тем раньше начинается вторичный рост скорости стримера.

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

На рис. 6 слева показаны линейные графики изменения напряжённости поля в процессе прорастания стримера к катоду.

Непосредственно перед головкой стримерного канала формируется передний фронт - область максимальной напряжённости поля (около 350 кВ/см), толщиной примерно 50 нм, из-за чего здесь образуется узкая область повышенной ионизации (см. рис. 4а). Поле перед головкой спадает плавно, подобно полю от заряженной сферы.

Напряжённость поля в плазменном канале на два порядка меньше, чем на фронте. Поле в нём слишком низкое, и ударной ионизации в канале не происходит. Здесь ключевую роль играет гибель электронов в результате рекомбинации и прилипания. На рис. 6 видно, что поле в хвостовой части канала вблизи анода примерно на порядок выше поля в остальной части канала. Это вызвано тем, что вблизи анода толщина канала меньше (рис. 4в), и чтобы обеспечить подпитку движущейся головки стримера, т. е. пропустить такой же ток через меньшее сечение канала, требуется повышенное поле в этой области проводящего канала.

Влияние дополнительных источников ионизации на развитие стримера. Были рассмотрены модели с разными значениями коэффициента мощности ионизующего излучения: д = 10~6 и д = 10~9. Как видно на рис. 5, 6, увеличение д приводит к более раннему переходу на самостоятельный этап развития стримера. Скорость стримера для приведённых значений д отличается не более чем в 2 раза, напряжённость поля на поверхности головки - не более чем в 1,4 раза, а эффективный диаметр головки практически не зависит от д. При этом разница в значениях д составила 3 порядка. Таким образом, подтверждается приведённое в [2, 8] утверждение о том, что скорость

02

к

bf :

0,2 0,4 0,6 0,8 1 Длина, мм в

t = 4,8 нс

Т

Анод

Длина стримера, мм

Т

Катод

c d

■*—► без фотоинозации

1--- с фотоионизацией (g = 10-6, l = 103 м-1)

2 -. с фотоионизацией (g = 10-6, l = 104 м-1)

3--- с фотоионизацией (g = 10-6, l = 105 м-1)

4 .. с фотоионизацией (g = 10-9, l = 103 м-1)

5 ...... с внешним источником

(т = 10-11 с, n = 107 м-3)

max

Рис. 6. Серия графиков напряжённости электрического поля вдоль оси симметрии (а) в разные моменты времени для модели с фотоионизацией и огибающие (б) подобных серий графиков для всех моделей; поле в один момент времени (в); схематично можно выделить участки:

а — хвостовая часть канала, Ь — канал, с — головка, ё — передний фронт стримера

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

В модели также варьировался и коэффициент поглощения фотоизлучения I. С одной стороны, увеличивая I, мы увеличиваем мощность источника электронов 5^ (3), что может ускорить стример; с другой стороны, затравочные электроны образуются ближе к фронту, что должно замедлить стример. Так, если характерная длина поглощения I-1 была порядка и больше радиуса головки (1 мм и 100 нм), влияние I на параметры стримера было слабое. Однако уменьшение I-1 до 10 нм привело к почти двукратному падению скорости стримера и существенному росту напряжённости поля перед головкой. Очевидно, это связано с тем, что, уменьшив характерную длину поглощения (т. е. прозрачность газа для данного излучения), мы уменьшили расстояние от фронта, на котором поглощаются фотоны и образуются затравочные электроны, а значит, уменьшили длину развития лавинного процесса перед фронтом стримера. Область локализации фотонов в этом случае сконцентрирована в тонком слое вокруг стримерной головки (рис. 7в).

Наличие только внешнего источника ионизации привело к развитию стримера, качественно и количественно весьма схожего с моделью с фотоионизацией при одной из

1//= 10 мкм

0,2 мм

1

0,8 0,6

0,4

0,2

0

Рис. 7. Концентрации фотонов в расчётной области, нормированные на максимум, для моделей с разным коэффициентом поглощения I в момент времени £ = 4,8 нс

б

а

в

комбинаций I и д (см. рис. 5, 6). Это говорит о том, что тип источника затравочных электронов вокруг головки качественно не влияет на картину развития стримера - необходимо лишь его наличие, чтобы вокруг продвигающейся головки постоянно возникали затравочные электроны. При исследованных мощностях излучения качественных отличий в развитии стримеров выявлено не было.

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

В работе в рамках компьютерной модели стримера рассчитана структура короткого катодонаправленного стримера в слабонеоднородном поле.

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

На стадии сформировавшегося стримера ионизация газа идёт преимущественно в тонком полусферическом слое непосредственно перед головкой, а в канале преобладает прилипание электронов.

Ионизация обеспечивает постоянное удлинение проводящего стримерного канала и перемещение окаймляющего головку слоя объёмного заряда на внешнюю границу текущего состояния проводящей головки, т. е. в направлении распространения стримера.

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

Рассчитанная скорость прорастания стримера в воздухе на самостоятельном этапе развития составила около 1000-1600 км/с, что согласуется с данными некоторых экспериментов 1000 км/с) [9], несмотря на различия в длине стримера: в эксперименте измеряется средняя скорость стримера на отрезке длиной как минимум несколько сантиметров [8], а в представленной модели длина стримера достигает лишь миллиметра.

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

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

Из результатов моделирования видно, что радиус головки сформировавшегося стримера на начальных этапах развития монотонно растёт. На масштабах порядка 1 мм и времени порядка 5 нс реакция прилипания электронов и наличие отрицательных ионов практически не влияют на процесс распространения стримера.

Литература

1. Luque A., Ratushnaya V., Ebert U. Positive and negative streamers in ambient air: modeling evolution and velocities // J. Phys. (D). 2008. Vol. 41. 234005.

2. Базелян Э. М., Райзер Ю. П. Физика молнии и молниезащиты. М., 2001. 320 с.

3. Стишков Ю. К., Самусенко А. В. Особенности распространения электронных лавин

в неоднородных электрических полях // Вестн. С.-Петерб. ун-та. Сер. 4: Физика, химия. 2009. Вып. 3. С. 36-44.

4. Райзер Ю. П. Физика газового разряда. М., 1992. 536 с.

5. Dutton J. A survey of electron swarm data // J. Phys. Chem. Ref. Data. 1975. Vol. 4. N 3. P. 577-856.

6. Gallagher J. W., Beaty E. C., Dutton J., Pitchford L. C. An annotated compilation and appraisal of electron swarm data in electronegative gases // J. Phys. Chem. Ref. Data. 1983. Vol. 12. N 1. P. 133-152.

7. Королёв Ю. Д., Месяц Г. А. Физика импульсного пробоя газов. М., 1991. 224 с.

8. Базелян Э. М., Райзер Ю. П. Искровой разряд. М., 1997. 320 с.

9. Winands G. J. J. Efficient streamer plasma generation: Ph. D. thesis. Netherlands, 2007.

10. Luque A., Ebert U., Montijn C., Hundsdorfer W. Photoionization in negative streamers: fast computations and two propagation modes // Appl. Phys. Lett. 2007. Vol. 90. 081501.

Статья поступила в редакцию 20 мая 2010 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.