М.С. Апфельбаум
ОБ ОДНОЙ МОДЕЛИ РАСЧЕТА ИОНИЗАЦИОННОГО ВЫДЕЛЕНИЯ ЭНЕРГИИ В СЛАБОПРОВОДЯЩЕЙ ЖИДКОСТИ ПРИ ЕЕ ВЗАИМОДЕЙСТВИИ С ПРЕДПРОБОЙНЫМ ЭЛЕКТРИЧЕСКИМ ПОЛЕМ
Институт теплофизики экстремальных состояний объединенного института высоких
температур Российской Академии наук, Ижорская, д. 13,строение 2, г. Москва, 125412, Россия msa@ihed.ras.ru
Введение
Достаточно сильные отклонения от закона Ома для твёрдых слабопроводящих сред в пред-пробойных постоянных по времени приложенных полях плоского конденсатора экспериментально обнаружил Пуль [1]. Позже аналогичный нелинейный рост электрического тока с ростом высокого напряжения как в однородном, так и в сильно неоднородном приложенном электрическом поле, экспериментально фиксировался как в жидких (типа трансформаторного масла), так и газообразных сла-бопроводящих средах.
Теоретически рост проводимости рассматриваемых сред с ростом модуля напряженности электрического поля обосновал Френкель в [2], используя уравнение Аррениусовского типа для объемной скорости ионизации с учетом снижения потенциала такой ионизации сильным электрическим полем. Аналогичное снижение работы выхода в предпробойном поле определено Шоттки. Онзагер [3] получил математически близкое к френкелевскому выражение для зависимости предпробойной скорости диссоциации слабых электролитов от модуля напряжённости электрического поля методами физической кинетики. Постулируя выполнение закона химического равновесия (действующих масс) для рассматриваемых сред типа трансформаторного масла, Онзагер [3] и Френкель [2] получили и зависимости их проводимости от напряженности электрического поля вплоть до её пробойной величины. При этом ими не учитывались влияние дрейфа зарядов и движения жидкости в таком поле. Кроме того, в [2, 3] не учитывается их диффузия, обусловленная градиентом объёмных концентраций зарядов. В работах известных учёных теоретически не исследовались временные изменения предпробойных электродинамических характеристик, наблюдаемые при проведении опытов даже с высоковольтным источником постоянного приложенного напряжения и малой мощности. Отметим, что случаи использования источника переменного или импульсного по времени высокого напряжения с немалой разрядной мощностью в работе не рассматриваются.
С учётом таких макроскопических процессов были выведены электрогидродинамические уравнения предпробойного образования зарядов в слабопроводящих жидких средах типа трансформаторного масла. Показано [4, 5], что для применения формул Френкеля и Онзагера при расчёте установившейся проводимости исследуемых сред достаточно сохранения условия их электрической (плазменной) квазинейтральности в рассматриваемых полях. При этом плотность образующегося в предпробойных полях объёмного заряда должна быть гораздо меньше плазменной плотности самонейтрализующегося фона зарядов разных знаков. Тогда по [4] система электродинамических уравнений зарядообразования может решаться независимо от системы гидродинамических при малости и соизмеримости скорости течений среды по сравнению со скоростью дрейфа ионов в электрическом поле. Отметим также, что в гидродинамических уравнениях движения изучаемых достаточно плотных и теплопроводных сред в сильном электрическом поле, которые рассматриваются в работе, следует учитывать влияние вязкости, например по [6]. В таких полях молекулы жидкостей типа трансформаторного масла и имеющихся в них примесей частично распадаются на ионы. Считается, что свободные электроны в указанных средах до образования в них электрических дуг и электрического пробоя в основном прилипают к нейтральным молекулам. Поэтому дополнительная макроскопическая электронная компонента при описываемом теоретическом моделировании предпробойных явлений в исследуемых средах нами не рассматривается.
© Апфельбаум М.С., Электронная обработка материалов, 2009, № 4, С. 38-46.
В случае плоского высоковольтного конденсатора квазистационарная экспоненциальная вольт-амперная характеристика достаточно легко получается из френкелевской формулы для проводимости и уравнения сохранения заряда в дифференциальной форме. В случае сферического конденсатора распределение установившегося потенциала электрического поля в его межэлектродном промежутке удовлетворяет получаемому из этого закона сохранения заряда и френкелевской проводимости достаточно сложному обыкновенному дифференциальному уравнению первого порядка с интегральным условием равенства разности этих потенциалов приложенному напряжению на обкладках конденсатора.
Из аналитического решения этого уравнения, полученного в [4] и удовлетворяющего такому условию, следуют омическая линейность вольт-амперной характеристики в слабых полях и её квад-ратичность в предпробойных, наблюдаемых при их сильной неоднородности. Полученное решение описывает математически и физически более строго, чем униполярные решения для таких полей в случае разреженных газов типа таунсендовского, переходной от линейного к квадратичному участок вольт-амперной характеристики. В этом решении вместо таунсендовской эмпирической константы напряжения зажигания короны для расчёта таких характеристик используется справочная или аккуратно определяемая по измерениям омических вольт-амперных характеристик, а в некоторых случаях рассчитываемая по составу молекул среды константа её низковольтной проводимости. В полученной нами теоретической высоковольтной вольт-амперной характеристике учтено, в отличие от таунсен-довских, влияние температуры и зарядового состава молекул среды на эти предпробойные характеристики. Отметим также, что константу типа напряжения зажигания короны для рассматриваемых сред использовать при получении расчётных вольт-амперных характеристик менее целесообразно, чем для разреженных газов, ввиду имеющихся трудностей по экспериментальной регистрации такого напряжения в достаточно плотных, вязких и теплопроводных средах. При этом актуальными являются задачи расчёта ионизационного выделения энергии, приводящего к образованию разрядных дуг в таких средах и к их пробоям в окрестностях высоковольтных электродов. Для высоковольтного электрода, близкого к точечному, такая энергия нами ранее уже рассчитывалась [5]. Актуальными являются задачи расчёта такой энергии для диагностики электрических пробоев вплоть до аварий в маслонаполненных трансформаторах. Ниже будут получены аналитическое решение электрогидродинамических дифференциальных уравнений для такой нестационарной задачи в случае плоскопараллельных электродов с межэлектродным расстоянием d<<yfS (здесь S - площадь поверхности пластины конденсатора) и численное решение уравнений электродинамического приближения для аналогичной задачи в случае коаксиальных цилиндрических электродов с межэлектродным расстоянием гораздо меньшим, чем длина их образующих. По результатам решений будут оценены время, за которое происходит выделение энергии в нагретом цилиндрическом слое масла в трансформаторе у его поверхности, приводящее как к пробоям масла, так и к авариям в трансформаторах ввиду пробоев, а также порядок величины такого выделения.
Теоретическая модель и исходные уравнения
Для объемных концентраций зарядов и, возможно, имеющихся примесей выполняются следующие неравенства в случае слабопроводящих сред:
1п±<< Па; (1)
\пр <<па.
Объёмная скорость ионизации с последующим прилипанием свободных электронов к молекулам рассматриваемых жидкостей или их частичной диссоциации - это термодинамическая функция объёмных концентраций молекул среды и, возможно, имеющихся в ней молекул примесей, а также температуры и напряжённости электрического поля. Рекомбинацию зарядов ввиду прилипания свободных электронов в рассматриваемых жидкостях предполагаем только ионной (двухчастичной):
Ж=Ж!(Па,Пр, Т, \Е\)=Щ(Па,Пр, Т, 0)К\Е\) (2)
Жг = Кгп+п- (3)
= 2(К+ ь- )\е\. 880
Выражение, используемое для Кг, в 1903 году получено Ланжевеном. Онзагер, применяя методы физической кинетики, обосновал [3] слабое влияние электрических полей вплоть до полей пробоя на величину ионного коэффициента рекомбинации. В [3] описано достаточно сильное влияние электрического поля на величину константы (коэффициента) диссоциации рассматриваемых слабых
жидких электролитов. При этом для коэффициента рекомбинации в [3] было также получено второе из уравнений (3), а для монотонно возрастающей зависимости константы диссоциации от модуля напряжённости макроскопического электрического поля выражение в виде аналитической и в нуле функции Бесселя первого порядка мнимого аргумента, близкой в предпробойных полях к неаналитической в нуле френкелевской экспоненте от корня этого модуля.
Для коэффициента диффузии свободных зарядов используем соотношение Нернста-Таунсенда-Эйнштейна:
D+ =
квТЬ±
И
(4)
Вид непрерывной скалярной функции f векторного аргумента в (2) используем, с учётом вышеизложенного, по Френкелю [2]:
3
\^е\2 . (5)
1
f (E) = exp(P|E|2); р =
4
nsBgквТ
В этих условиях уравнения зарядообразования выглядят по [4] так:
к Tb
+(V,Vq)-kBTbAq+(E,Va)=, Аф = --
д t
Zlel
0
88
0
(6А)
f+(V,Va)+b2(E,Vq)-ZeA A°-
k^Tb a2^
(6B)
88,
exp(P|E|1// 2)+8|-=0.
0 0
Систему электродинамических уравнений (6), где учтены диффузионный и конвективный переносы ионных компонент в движущейся среде и уравнение Пуассона для потенциала макроскопического электрического поля, следует замкнуть математически гидродинамическими уравнениями для её нейтральной компоненты, близкими к уравнениям для всей смеси при выполнении соотношений (1). При этом дифференциальное уравнение, описывающее неразрывность массы несжимаемых сред, предстанет в виде
divV = 0. (7)
К (7) согласно [6] достаточно добавить только уравнение движения, не рассматривая при этом уравнение сохранения энергии и уравнения состояния. Эти уравнения, например, при теоретическом анализе охлаждения трансформаторов следует решать, используя решения выписанных выше электрогидродинамических уравнений (при этом вместо уравнения энергии можно использовать стандартное уравнение притока тепла [6], не учитывая слабого джоулева нагрева и влияния электрической поляризации, а также диссипации энергии, обусловленной вязкостью жидкости). Поле сил электрической поляризации связанных зарядов только перераспределяет давление рассматриваемой среды, что позволяет использовать вместо импульсного уравнения уравнение переноса вихрей [6]:
X = rot V ; (8)
Р dt "P(X 'V)V-PVAX = [Vq,E].
В начальный момент времени
q(t=0)=V(t=0)=0; a(t=0)=a0.
(9)
(10)
Аналитическое решение для установившегося распределения потенциала высоковольтной сферы выписанных выше уравнений получено в [4] при выведенных из них там же условиях квазинейтральности рассматриваемых слабопроводящих сред. Оно выглядит так:
ф(г ) =
f
I
1
^ 2
v4na0у
88,
Y.
E
exp
1 та0
E2 =■
Л"
sign (ф(г0 ));
4na0 r
= и.
(11)
Аналогично выглядит решение в случае сферического высоковольтного конденсатора с заданной разностью потенциалов и между его обкладками. Вышеупомянутая линейность вольт-
амперной характеристики слабых полей и квадратичность сильно неоднородных предпробойных получаются из (11) соответствующими предельными переходами. Распределение же предпробойного установившегося потенциала сферической симметрии по (11) порядком отличается от лапласовского ввиду влияния поля объёмного заряда. Вместе с тем в случае плоского высоковольтного конденсатора электрическое поле вне приэлектродных погранслоёв нарушения квазинейтральности электрических зарядов остаётся однородным, что подтверждается из закона сохранения заряда. Зависящая от граничных условий на электродах неоднородность пространственного распределения полей образуется при наложении предпробойного напряжения наряду с объёмным зарядом только в таких слоях. Для этих случаев результаты разных экспериментов в поле плоского конденсатора подтверждают формулу расчета квазистационарной предпробойной вольт-амперной характеристики в условиях френкелевской квазинейтральности рассматриваемых сред, которая теоретически получается из соответствующего интеграла одномерного уравнения сохранения заряда в дифференциальной форме. Тогда в случае определения константы такого интегрирования по электрическому току эта зависимость выглядит следующим образом:
U d
exp
I £
(12)
Нестационарные решения: аналитическое для случая плоского высоковольтного конденсатора и численное для цилиндрического.
Независимо от вихревых характеристик предпробойных течений рассматриваемых сред в межэлектродном промежутке плоского конденсатора электродинамическим уравнениям (6) и начальным условиям (10) удовлетворяет достаточно простое частное решение:
I = aSU / d,q = 0,
E
= U/d;
2Л/х
880 [°0т + 880 -(880 -°0т) ^ ,
а=—г-?-; (13)
т| а0Т + 880 +(880 -о0т)e
т= 880 ^exp (0,5р
E
1 / 2
Решения этих уравнений, удовлетворяющие электрохимическим (в случае диссоциации) или ионизационным граничным условиям на пластинах конденсатора, а также интегральному условию независимости от времени разности потенциалов между обкладками пластин, близки за пределами неквази-нейтральных погранслоёв к решению (13). Теоретически это может быть обосновано описанными, например, в [7] методами асимптотических погранслойных разложений по малым параметрам, соответствующим режимам квазинейтральности, при членах уравнений (6) со старшими пространственными производными. Скорости предпробойных течений рассматриваемой слабопроводящей среды при этом могут быть и малыми, и соизмеримыми по сравнению со скоростями дрейфа ионов в электрическом поле. Выражения для таких критериев подобия описываемых явлений выписаны, например, в [4, 5, 8]. Там же проведены оценки величин этих параметров для исследуемых предпробойных режимов. При этом показано, что для достаточно широкого класса рассматриваемых сред условия электрической (плазменной) квазинейтральности могут выполняться и при наложении сильных электрических полей. Тогда независимо от граничных условий на электродах ампер-временную характеристику плоского конденсатора в случае источника постоянного по времени высокого напряжения можно в одномерном приближении рассчитывать по (13), откуда следует, что с ростом приложенного напряжения (до пробоя) усиливаются и временные отклонения предпробойных токов до установления. Из (13), как и следовало ожидать, формула Френкеля [2] для зависимости проводимости от напряжённости поля получается в виде временного асимптотического приближения. Далее следует, что для идеальных диэлектриков известная формула становится практически неприменимой, ибо время релаксации рассматриваемых сред при этом становится достаточно большим. Оно является и временем накопления ионизационной энергии, большая величина которой приводит к электрическим пробоям. Как следует из результатов [8], в условиях квазинейтральности формулу (13) можно получить также при разных величинах подвижностей положительных и отрицательных зарядов. Аналитического аналога формулы (13) в случае цилиндрического и сферического конденсаторов, а тем более
для двумерных конфигураций высоковольтных электродов типа "провод над плоскостью " и трёхмерных типа "игла-плоскость" из уравнений (6) получить не удаётся даже без учёта влияния гидродинамики на электродинамические характеристики. В этих случаях расчёт ампер-временных характеристик до установления можно проводить численно, используя упрощённую для решения внешних по [7] задач систему уравнений (6) с начальными условиями (10). Аналог же формулы (13) для сильных неоднородных электрических полей любой сложной конфигурации может быть получен только при сохранении их лапласовости (при этом соответствующая (13) общая формула, в которой f определяется не обязательно по (5), была выписана в [9]). Для получения простейшего численного решения при цилиндрическом конденсаторе и последующих расчётах нелинейных вольт-амперных и предпробойных ампер-временных характеристик использовалась система нестационарных одномерных электродинамических уравнений, получаемая из уравнений (6А) с применением в одномерных случаях уравнения Гаусса для напряжённости электрического поля цилиндрической симметрии вместо уравнения Пуассона для её потенциала:
дq + Б -да = q -а дt дг 8 -80'
да = ар -ехр(в-Б172)-а2 . (14)
дt 8 -80
дБ + Б = q дг г 8 - 80
При получении системы уравнений (14) в левых частях (6А) сохранены члены, описывающие нестационарные изменения квазинейтральной проводимости. При этом можно анализировать и образование объемного заряда с ростом от нулевого в начальный момент времени до установления (ранее на возможность образования такого заряда указал Остроумов [10]). Исключенные, как и в [11], члены с пространственными старшими производными влияют при выполнении критерия квазинейтральности только на изменение искомых распределений электродинамических характеристи в неквазинейтральных тонких слоях у границ объёмов, заполненных рассматриваемыми средами. При получении численного решения системы уравнений (14) использовались начальные условия (10) для проводимости и плотности объемного заряда, а также условие лапласовости электрического поля цилиндрического конденсатора в начальный момент времени:
Б (г,0) = и . (15)
г - /п(
Г2
В (15), в отличие от аналогичных задач, решённых численно в предыдущих работах [12, 13], г1 и г2 - внутренний и внешний радиусы цилиндрического слоя жидкости типа слоя масла в трансформаторе (у его поверхности), охлаждаемом таким маслом. При выборе шага сетки по радиальной координате считалось, что толщины образующихся погранслоёв нарушения квазинейтральности малы по сравнению с этим шагом. Расчеты производились от внутренней цилиндрической поверхности рассматриваемого слоя к внешней. В расчётах отмечалось, что эта поверхность заряжена положительно. При смене знака её заряда задача расчёта выделения энергии в таком слое сводится к предыдущей со сменой знака по напряженности поля и плотности объемного заряда. Последнее утвеждение обосновано численно в [11], где с использованием обобщений уравнений типа (14) на неизотермический случай рассматривалась аналогичная нестационарная задача об электрическом поле в окрестности остывающей дисперсной сферической частицы. Но, в отличие от неизотермической, при решении описываемой в настоящей работе изотермической задачи нестационарный интеграл сохранения заряда
дБ ^ I
880 --+ аБ =--(16)
0 дt 2пг/
не удаётся использовать для понижения порядка системы уравнений (14).
При расчетах, результаты которых будут представлены ниже, интеграл (16) применялся для контроля численного счета разностных уравнений, получаемых из дифференциальных уравнений (14) в безразмерной форме. Кроме того, по (16) определялась искомая зависимость электрического тока от времени до установления. В (16) / - длина образующей рассматриваемого цилиндрического слоя (/=1,626 м для поверхности трансформатора).
чин:
Переход к безразмерным уравнениям осуществлялся изменением искомых и исходных вели-
~ * ~ Ей ~ цй2 * = - ,Е = — ,ц = ■
и
Л / 2
88пи
в=в(й) '=7
(17)
а =
а
а 0ехр
~
в
2
V У
Тогда получаемая из (14) система безразмерных уравнений примет вид
д Е Е ~
—+- = ц,
д г г
дц Еда ~ ~ — + —— = -Ц а,
д * д г
(18)
да
-1]
—- = ехр \ в
д ~
Начальные условия для безразмерных переменных выглядят следующим образом:
С ~ ^ в
а(г, 0) = ехр
2
V У
ц(г,0) = 0,
(19)
Е (г, 0) = ■
г 1п
гу
V г2 У
При этом интеграл (16) в безразмерном виде выглядит так:
д Е ~~ I
-+аЕ=—,
(20)
д *
где
I = ■
2 л/ а0иехр
в
2
V У
Из аналогичной задачи для плоского конденсатора и её аналитического решения (13) следует, что установление происходит при времени * ~ 1. Поэтому для получения искомого численного решения выбирался шаг по времени А * = 0,02 . Установление происходило примерно за 50 шагов. Для проведения конкретных численных расчетов по радиальной координате выбирался слой у поверхности трансформатора с внутренним радиусом 0,25 м и внешним, равным радиусу поверхности трансформатора (0,42 м). То есть считалось, что охлаждаемый слой у поверхности трансформатора, заполненного маслом, можно аппроксимировать цилиндрическим слоем толщины й (й = г2 - г!).При расчётах из справочников выбирались диэлектрическая проницаемость среды (2,24) и ее низковольтная
1
проводимость (10 11-). Среди сред с такой низковольтной проводимостью и диэлектрической
т
проницаемостью можно указать, например, очищенное трансформаторное масло.
1,2Е+008 6Е+007 0
0 0,004 0,008 0,012
Рис. 1. Предпробойная ампер-временная характеристика цилиндрического охлаждаемого слоя у поверхности маслонстолненного трансформатора (пояснения в тексте)
I, in А
"" и, МУ
Рис. 2. Квазистационарная вольт-амперная характеристика трансформаторного масла для генерируемых трансформатором наиболее высоких нпапряжений
Температура среды, в отличие от предыдущих работ [12, 13], в которых численно решались близкие задачи при цилиндрической симметрии, считалась не комнатной, а близкой к температуре нагретого масла у обмоток трансформатора (400К). Приложенное напряжение выбиралось равным 10 МУ(по порядку величины самого высокого напряжения, генерируемого трансформатором). При счете по пространственной координате шаг равнялся 0,1. Дифференциальные уравнения (18) заменялись соответствующими разностными [4, 9, 14]. Из них на каждом следующем временном шаге расчитывались сначала пространственное распределение проводимости, затем плотности объемного заряда. Для получения распределения напряженности электрического поля в межэлектродном промежутке использовалось уравнение Гаусса для поля цилиндрической симметрии. После выполнения этих расчетов на каждом временном шаге по разностному аналогу интеграла (20) рассчитывались временные изменения полного тока до установления. На каждом шаге пересчитывались соответствующие матричные массивы пространственно-временных распределений искомых переменных. После чего они подставлялись в разностный аналог интеграла (20). Отметим, что результаты расчетов на первых шагах по времени уже описывались в [4, 9]. При этом там же описано отличие результатов по расчетам q и Е для начальных лапласовских высоковольтных полей плоского и цилиндрического конденсаторов. Кроме того, в [3, 9] отличие искомых электродинамических характеристик в случае однородности предпробойного приложенного поля от сильной неоднородности проанализировано и в установившихся q и Е. Отличие же структуры течений в сильном неоднородном приложенном поле от ячеистых течений в плоском конденсаторе описаны ранее. Поэтому в данной работе, исходя из полученных аналитических и численных решений, опишем, как и в [13], сходство искомых характеристик для случая однородного и сильно неоднородного приложенных полей. Оно состоит в достаточно быстрой в обоих случаях временной эволюции и быстром накоплении энергии, приводящих к электрическим пробоям, например в трансформаторах. Из формулы (13) следует, что время выделения такой энергии, например для трансформаторного масла, составляет доли секунды. Отметим, что
счетная устойчивость описанной в настоящей статье разностной схемы не исследовалась. Но это может быть сделано по аналогии с анализом устойчивости [14] для близких одномерных задач.
Изменение <5 обусловлено влиянием предпробойного поля на скорость реакции диссоциации (ионизации) рассматриваемых сред. При этом для исходных данных, при которых проводились численные расчеты, схема которых описана выше, достаточно быстрое установление проводимости происходило за достаточно малое время (тоже доли секунды). Такие же времена установления проводимости получаются и из аналитического решения (13). Кроме того, это сходство касается и ампер-временных характеристик.
На графике рис. 1 представлена временная эволюция предпробойного тока до установления с учетом результатов численных расчетов уравнений (18) и их интеграла (20). Этот график геометрически близок к кривым предпробойных нестационарных токов плоского конденсатора, которые можно построить по формуле (13). На графике рис. 2 представлена квазистационарная вольт-амперная характеристика, аппроксимируемая квадратичной зависимостью, как и в [12], но для более высоких порядков U, генерируемых трансформатором. Для предпробойных квазистационарных режимов пространственное (радиальное) распределение напряженности поля цилиндрической симметрии из закона сохранения заряда получаем:
E
ехр
^Р *1^ I
2
E
J
2п«Лг
(21)
где I - квазистационарный ток.
Интегральный заряд Q вне приэлектродных слоев нарушения квазинейтральности
Q=|qdV.
(22)
V
где q определяется из последнего уравнения (14) и уравнения (21), вычисляется аналитически. При этом следует:
0= XI[ехрС-^^) - ехр()]; Ех=Е(гх\ Е2=Е(г2). (23)
Тогда по (23) для данных, при которых проводились описанные численные расчеты при установлении, по порядку величины для трансформаторного масла можно достичь 0= 10-2 Кулон даже без учёта объёмного заряда, образующегося в неквазинейтральных погранслоях (для дуговой электрон-ионной плазмы такие слои известны как легмюровские [15]). Этот теоретический результат свидетельствует о высокой степени зарядки в высоковольтных установках. В частности, в трансформаторах, что подтверждается экспериментально. Поэтому при разработке конструкций таких установок и их эксплуатации следует предусматривать дополнительные устройства нейтрализации такого заряда или его отводки в связи с возможными электрическими пробоями. Для больших объемов масла, например в трансформаторах, генерирующих напряжение и= 10МУ, порядок накапливаемой (в нагретом слое масла у его поверхности) энергии Ж=0и достигает 0,1 МДж даже без учета энергии, выделяемой в ленгмюровских пограничных слоях нарушения квазинейтральности [15]. Удельная плотность выделяемой энергии в случае трансформаторного масла по порядку величины для приводимых выше исходных данных поверхности трансформатора получается близкой к порядку взрывных (1 МДж/кг) [16]. То есть расчет выделения в жидкости энергии, обусловленного предпробойными электродинамическими и ионизационными процессами по предложенной модели, можно использовать для диагностики пробоев, наблюдаемых в высоковольтных устройствах, например в трансформаторах.
Основные обозначения:
и - приложенное электрическое напряжение; й - межэлектродное расстояние; е - заряд протона; кв - постоянная Больцмана; В - коэффициент диффузии; Ь - коэффициент подвижности; q - объемная плотность заряда; V - скорость движения среды; п - концентрация; Е - напряженность электрического поля; - площадь поверхности электрода; г - радиальная координата; t - время; т - время релаксации заряда; х - завихренность; е - диэлектрическая проницаемость; I - сила электрического тока; Т — абсолютная температура; / — функция скалярного или векторного аргумента; Ж — объемная скорость образования зарядов или нейтралов; Ъ — валентность иона; о — проводимость среды
ЛИТЕРАТУРА
1. Poole H.H. On the dielectric constant and electrical conductivity of mica in intense field // Philosophical magazine. 1916. S. 6. V.2. № 187. Р. 112-129.
2. Френкель Я.И. К теории электрического пробоя в диэлектриках и электронных полупроводниках // Журнал экспериментальной и теоретической физики. 1938. Т. 8. Вып. 12. С. 1291-2001.
3. Onsager L. Deviation from Ohm s law in a weak electrolytes // Journal of chemical physics. 1934. № 2. Р. 599-611.
4. Апфельбаум М.С. Об одной расчётной схеме электрогидродинамических течений // Электрохимия. 1986. Вып. 11. С. 1463-1471.
5. Апфельбаум М.С. Задача о точечном электровзрыве слабопроводящих сред // Электронная обработка материалов. 2000. № 6. С. 31-42.
6. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидромеханика. М., 1989.
7. Васильева А.Б., Бутузов В.Ф. Асимптотические разложения решений сингулярно возмущённых уравнений, Москва, 1973.
8. Апфельбаум М.С.,Козыренко В.Е. Физическое моделирование и расчёт электрогидродинамических течений // Электрохимия. 1991. Вып. 7. С. 855-863.
9. Апфельбаум М.С. Ионизация и течение слабопроводящей жидкости в неоднородном электрическом поле // Электронная обработка материалов. 1988. № 1. С. 60-65.
10. Остроумов Г.А. Взаимодействие электрических и гидродинамических полей. М., 1979.
11. ApfelbaumM.S., Apfelbaum E.M. On electric field distribution in the vicinity of particles in a weakly ionized disperse medium // Chem.Phys, Report. 2000. Vol.18. Р. 2313-2328.
12. Апфельбаум М.С. Одномерные стационарные задачи расчёта предпробойных вольт-амперных характеристик слабоионизованных сред // Электронная обработка материалов. 2005. № 2. С. 50-53.
13. Апфельбаум М.С. Одномерные нестационарные задачи расчёта предпробойных электродинамических характеристик слабоионизованных сред // Электронная обработка материалов. 2008. № 3. С. 36-45.
14. Марчук Г.И. Методы вычислительной математики. М., 1970.
15. Tonks L., Langmuir I. A general theory of the plasma of an arc // Physical review. 1929.V.14. P.876-922.
16. Зельдович Я.Б., Райзер Ю.П. Физика ударных волн и высокотемпературных гидродинамических явлений. М., 1966.
Поступила 23.03.09
Summary
A theoretical model of electrohydrodynamic pre-breakdown phenomena in a weakly ionized media is proposed. Differential equations and initial conditions for describing are written. The 1D analytical and numerical solutions of these equations for calculations current time characteristics are given. The ionization breakdown energy is calculated.