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

Математическое моделирование переходных процессов в разрядном контуре и разрядной камере электрогидравлической установки

В. М. Косенков, В. М. Бычков

Институт импульсных процессов и технологий НАН Украины, пр. Октябрьский, 43-А, г. Николаев, 54018, Украина, e-mail: v.m.kosenkov@gmail.com

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

Ключевые слова: разряд в воде, разрядный контур, математическая модель, переходные процессы.

УДК 519.87:537.528

ВВЕДЕНИЕ

Электрические разряды в воде находят широкое применение в научных исследованиях и в практических приложениях [1-4]. Для осуществления разрядов в воде используют разнообразные электрогидравлические установки, одним из основных элементов которых является разрядный контур. От эффективности передачи энергии заряда конденсаторной батареи каналу разряда зависит качество выполнения различных технологических процессов. Поэтому изучению влияния параметров разрядного контура на преобразование энергии в нём посвящено много исследований (например, [5-8]). В этих работах высказывалось предположение, что канал разряда представляет собой прямой круговой цилиндр, который в процессе расширения не изменяет свою форму. Исследования [5, 9] об изучении трансформации формы канала разряда свидетельствуют, что это допущение справедливо только для каналов большой длины по сравнению с его максимальным радиусом. В остальных случаях оно не всегда оправдано. Большое влияние на форму и размеры канала разряда оказывают геометрические параметры высоковольтных электродов [10-12]. Однако при выполнении этих исследований не учитывали связь формы канала с процессами в разрядном контуре, считая временную зависимость выделения энергии в канале неизменной. Поэтому изучение влияния электродной системы и разрядной камеры на расширение канала разряда и процессов в разрядном контуре актуально в научном и практи-

ческом отношении. Следовательно, объектом рассмотрения должна быть система, включающая в себя разрядный контур и заполненную водой разрядную камеру.

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

ПОСТАНОВКА ЗАДАЧИ

Разрядный контур импульсных электрогидравлических установок содержит следующие основные элементы: конденсаторную батарею, разрядник, электродную систему и канал разряда между электродами. В зависимости от назначения электрогидравлических установок их разрядные камеры, в которых размещается электродная система, имеют различные формы и размеры. В ряде приложений они имеют цилиндрическую форму и соответственно обладают осевой симметрией. Разрядную камеру такого типа выбрали для моделирования в данной работе, чтобы упростить разработку математической модели и анализ процессов в разрядном контуре. Схема замещения разрядного контура и разрядной камеры представлена на рис. 1. Согласно этой схеме для осуществления электрического разряда между электродами 1 и 6 конденсаторную батарею C заряжают до заданного напряжения Щ. По мере заряда и увеличения напряжения

© Косенков В.М., Бычков В.М., Электронная обработка материалов, 2015, 51(2), 66-72.

на конденсаторной батарее происходит пробой газа в разряднике 5 и на высоковольтный электрод 1 подается напряжение, приводящее к электрическому пробою слоя воды между электродами 1 и 6. По образовавшемуся электропроводному каналу 4 в воде, заполняющей разрядную камеру 3, начинает протекать электрический ток, приводящий к разогреву вещества в канале с образованием низкотемпературной плазмы. То есть канал разряда 4 представляет собой проводник, сопротивление которого существенно изменяется в процессе разряда при разогреве плазмы и её расширении. В процессе разряда элементы разрядного контура, за исключением канала разряда, практически не изменяют своих характеристик. Поэтому общая индуктивность L и активное сопротивление Rsh этих элементов остаются постоянными. Для электрической изоляции стержня высоковольтного электрода 1 от корпуса разрядной камеры используется изолятор 2.

C L Кзи

_11--СГЬ

Рис. 1. Схема замещения разрядного контура и осесиммет-ричной разрядной камеры.

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

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

[5]:

Lft+(Rck+Rsh)i+C\dt_^ (i)

где L - индуктивность контура, Гн; t - сила тока в контуре, А; Rck, Rsh - соответственно сопротивление канала разряда и остальных элементов разрядного контура, Ом; C - емкость конденсаторной батареи, мкФ; t - время, с; U0 - напряжение заряда конденсаторной батареи, В.

Электрическое сопротивление канала разряда можно определить с помощью уравнения [5]:

_ A-IKy-1)

ck P .V -

где Рк - давление в канале разряда, Па; Ук - объем канала разряда, м3; 1к - расстояние между электродами, м; у = 1,26 - эффективный показатель адиабаты плазмы в канале разряда; А - искровая постоянная В2-с/м2 [5].

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

Давление в канале разряда можно определить из уравнения баланса энергии, в котором учитывается изменение внутренней энергии плазмы, работы расширения канала разряда и джоулева нагрева плазмы от протекающего через нее электрического тока [5]:

Л d (PV)+P*dV _ t2 R

Y -1 dt k dt

(3)

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

В результате анализа существующих математических моделей гидродинамических процессов в разрядной камере [6, 7, 9-12] выбрана модель, в которой используют предположение о потенциальном безвихревом движении идеальной жидкости. Гидродинамическая модель включает в себя двумерное линейное дифференциальное волновое уравнение гиперболического типа в частных производных. Для принятых допущений волновое уравнение в цилиндрической системе координат имеет вид [12]:

ö2y _ c2 öt2 _ С

1 (r sy)+а2ж

r ör ör öz 2

(4)

где у - потенциал скорости, м2/с; с - скорость звука в воде, м/с; г, z - пространственные координаты, м.

Скорость звука в воде определяли по формуле [12]:

\2 / „ \2

с2 _ ^-(х - !)f -2(х +1)

öy

ör

öy

Hz

(5)

(2)

lk ' k

где с0 - скорость звука в невозмущённой жидкости, м/с; х - эмпирический коэффициент, равный 7, 15 для воды.

Давление в жидкости определяли с помощью линейного интеграла Коши-Лагранжа:

Р=Ро-Р^ 3 (6)

где р - плотность жидкости, кг/м ; Р0 - начальное давление в разрядной камере, Па.

Составляющие вектора скорости жидкости и границы канала разряда определяли по градиенту потенциала:

V =

г дг '

V =

^ дг

(7)

(8)

1_д_( г =0 г дг \ дг Гдх2

(9)

ду дп

= 0.

(10)

На границе канала разряда выполняются условия на контактном разрыве [9]:

ик = иг,

р = р

гк

(11) (12)

где и - скорость жидкости на границе канала, м/с; Р^ - давление жидкости на границе канала, Па.

МЕТОД РЕШЕНИЯ ЗАДАЧИ

Непрерывную пространственно-временную область решения задачи (см. рис. 1) преобразовали в дискретную, задавая шаг изменения времени А/ и пространственных координат Дг и Дг. Номера временных слоёв дискретизации обозначили дискретной переменной п, а неподвижные пространственные слои по координатам г и х -дискретными переменными 1 и 1 соответственно. Схема дискретизации разрядной камеры показана на рис. 2.

Решение уравнений (1) и (3) выполняли по двухшаговому методу Эйлера-Коши [10, 11], имеющему второй порядок точности.

Для расчета гидродинамических процессов в канале разряда использовали уравнение Лапласа [12]:

1+1

1-1

Оно получено как частный случай уравнения (4), с учетом того, что скорость звука в плазме, заполняющей канал разряда, намного больше, чем в воде.

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

На границе раздела двух сред значение потенциала находили из условия, что стенка камеры абсолютно жесткая. Поэтому градиент потенциала по нормали к поверхности стенки п равен нулю:

1-1 1 1+1

Рис. 2. Схема дискретизации разрядной камеры.

На этапе предиктор вычисляли сопротивление канала разряда по формуле (2):

й„ = Л1к(у -1).

ск р "V"

к к

(13)

Затем определяли силу тока, используя дискретный аналог уравнения (1):

1 "+1 = + А[и0 - (Щк + Ял)1" ]-

ь

А/

С ■ ь

(14)

(Л" + А/ ■ 1"),

О п

где Л - интегральная сумма на временном слое п, вычисленная по методу трапеций.

Из дискретного аналога уравнения (3) вычисляли давление в канале разряда:

1

р"1 = (А ■ (у -1)*

V" - V"-1 х(М" (/) - рп к к )+р^п), А/

(15)

где N (/) = (Г )2 ■ Я"к - мощность энергии, выделяемой в канале разряда, Вт.

Начальную скорость перемещения границы канала разряда при заданном начальном давлении в канале Рк0 вычисляли по формуле [13]:

р - р

тт _ 1 к 0 ■* 0 ^ /,п —

Р ■С

(16)

Потенциал на первом и втором временных слоях полагали равным нулю, а на границе канала, на втором слое, определяли из интеграла Коши-Лагранжа (6):

у2,1 = ¥1,1 - р -^ /Р0. (17)

На третьем и последующих слоях времени предварительное значение потенциала жидкости на границе канала определяли по давлению, вычисленному по формуле (6):

^ = ¥",1 - (рГ1 - р,) ■А/ / Р0. (18)

Предварительно поле потенциала во внутренних точках жидкости определяли из уравнения (4) по схеме «крест» [14], имеющей второй поря-

V

док точности по временной и пространственной координатам и простой алгоритм реализации:

~И+1 ».2 2Г 1 г п , п

у«,} = А •с [2 А 2 (Г+1 • %+1, j + Г • %+1, у

2• Аг •г

(19)

-г-+1 • ^ у -2 • г • ^ у - г-1 • ^ у+г-1 • <1, у+

п ~ п , п + п ) , V, У+1 - 2 • V, у + V, у-1] + . п п-1

+г • у)+-¡¡2-]+2 • V, У - V, у.

Шаг по времени определяли из условия устойчивости Куранта-Фридрихса-Леви [11]:

А =

к • min(Аz; Аг)

(20)

72 • с '

где к < 1 - коэффициент запаса устойчивости схемы (число Куранта).

Поле потенциала, вычисленное по формуле (19), обеспечивает малую погрешность величины потенциала, но относительно большую погрешность его производных по времени и пространственным координатам, что приводит к большим дисперсионным погрешностям полей давления и скорости, определяемых по формулам (6)-(8). Чтобы уменьшить эту погрешность до необходимого минимума, выполняли регуляризацию поля потенциала х^1 по неявной разностной схеме [14]: п+1 ,,2 2(, г 1 г ~п+1

V,у =А с {к1 •[7"М"2"~(г+1 • ¥'+1,у +

2•Аг •г- ^

~п+1 ~п+1 ~ ~ п+1

+г • ¥/+1, у - г +1 • V,- у - 2 • г • V,- у -

-п+1 -п+1 -п+1 ч

г 1 • V; ,■ + г , • V- , ■ + г • V- 1 ) + - ,-1 гг, у ,-1 4-1, у , 4-1, у

~п+1 ~ ~п+1 ~п+1 V-, у+1 - 2 • V-, у + V-, у-1

(21)

Ае

+к2(г+1 •^+1, у + г • <+1, у- г+1 • < у

2•Аг •г

п. п п п п ч

-2 • г • V; ; - г 1 • V,- ; + г 1 • V ■ 1 ■ + г • V ■ 1 ■) + г 4, у I-1 4 , у 1-1 4-1, / г 4-1, у

> у - -1 ^ , у

п ~ п п

[I [I [I

Ч у+1 -2•У, у +У, у-1 + 2 п п-1

+-¡¡2-]} + 2 • V, у - V, у ,

где к1 и к2 - весовые коэффициенты (к1= 0,06; к2 = 0,94).

Определяли скорость границы канала разряда во всех точках пересечения координатных линий разностной сетки с его поверхностью из условия на контактном разрыве [13]:

иИ+1 = 0 5 • (ип + ип ) +

г (к,у ^г», у^^гй+1,

г ¡к, у' г Ук+1, у>

Р«+1 - 05 • (рп + рп )

1 к,у V-1 к, у т -1 к+1, у )

(22)

&, у й+1, у >

Р •с

Интегрируя по времени скорость канала, определили его радиус на п+1 слое. Интегрируя

площадь поперечного сечения канала по его длине, определили объём канала на п+1 слое.

На этапе «корректор» вычислили сопротивление канала разряда

( 1 л \

Rn:l = 2 А • 42(у -1)

1

1

рпуп V кУ к

р п+1тр п 1 ь ' ь

(23)

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

А/ '

Ь

и 0 - (+ Rл )-

А/ ( ^ + рИ+1

1 5"1 +М-—

Ь • С

(24)

Затем определили давление в канале разряда: Ркп+1 = ^(А/ • (у -1) • (Ы"(/) + #"+1 {/) -

к

уп - V«-1 ^

_ Пп к г к _ рп+1 * к к 1 к

А/

ур,п+1 - у,п

)+• у п),

А/

(25)

скорость границы канала:

^п+к1, у = 0,5 • (иПпи, у + ^к+1, у) +

Рп+ -0,5• (Рпу + Рп+, у)

(26)

¡k ,у ¡k+1, у'

Р •с

и потенциал на границе канала разряда на п+1 временном слое:

Уu = VIу -(Ркп+1 -Р,)•А//Ре. (27)

По полученным потенциалам вычислили поле скоростей Уг, Уz и давление в разрядной камере. Полученные характеристики использовали для определения поля скорости звука на п+1 временном слое по формуле (5).

По скорости движения границы канала разряда вычислили его новое положение на п+1 -временном слое ^у^1). Затем определили расположение границы канала по отношению к узлам дискретной сетки. При новом цикле действий исключали из процесса вычисления по формулам (19) и (21) точки, которые оказались внутри канала разряда. Благодаря этому приему удалось выполнить расчет гидродинамических процессов в разрядной камере с учетом подвижности границы канала разряда на неподвижной сетке. В точках сетки, расположенных внутри канала, потенциал вычисляли, решая краевую задачу для уравнения (9) на каждом временном слое по методу последовательной верхней релаксации [14]. Таким образом, задача сформулирована в физически и геометрически нелинейной постановке.

ТЕСТИРОВАНИЕ МАТЕМАТИЧЕСКОЙ МОДЕЛИ

Тестирование математической модели выполнили на основе результатов моделирования, приведенных в работе [6].

+

+

Rck, Ом

1,0 0,8 0.6 0,4 0,2 0,0

—---- 1 2 -

0 50 100 150 t, мкс

Рис. 3. Сопротивление канала разряда.

80 40 О -40

]

/Л ----2

' \

/

О 50 100 t, мкс

Рис. 5. Сила тока в канале разряда.

Rk, мм 40 30 20 10 О

-]

--2

О 100 200 300 1, мкс

Рис. 7. Радиус канала разряда.

В [6] был выбран радиус камеры 200 мм и длина 100 мм. Расстояние между электродами -100 мм, радиус электродов принят равным радиусу камеры. Напряжение заряда конденсаторной батареи - 50 кВ, индуктивность - 4 мкГн, емкость конденсаторной батареи - 16 мкФ, сопротивление контура - 0,04 Ом, искровая постоянная А [5] принята равной 105 В-с/м2. Полученные результаты представлены на рис. 3-8, на которых цифрой 1 отмечены данные из [6], цифрой 2 - расчет по модели, описанной в данной работе. Как следует из этих рисунков, предлагаемая математическая модель удовлетворительно согласуется с данными работы [6], несмотря на сделанные упрощающие предположения.

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

РЕЗУЛЬТАТЫ ИССЛЕДОВАНИЯ

Рассмотрен пример влияния дифракции ударных волн, генерируемых электрическим разрядом в воде, вблизи круговой границы торцов электродов на переходные процессы в разрядном контуре со следующими параметрами: С = 2 мкФ, Ь = 1 мкГн, и0 = 50 кВ, 1к = 25 мм, Яаи = 0,01 Ом. Радиус разрядной камеры составлял 200 мм, а её длина - 75 мм. Использовали

N, МВт 800 п 600 400 200 0

Л i

!j \ j ----2

г\ -

k ____^

О 50 100 t, мкс

Рис. 4. Мощность выделения энергии.

U, кВ

105

0 50 100 t, мкс

Рис. 6. Напряжение на канале разряда.

Рк, МПа 400300200-

0 100 200 300 t, мкс

Рис. 8. Давление в канале разряда.

электроды цилиндрической формы с плоскими торцами. Выполнено сравнение характеристик двух разрядных контуров, отличающихся радиусами электродов, равными 200 и 10 мм (рис. 9-12).

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

Радиус канала возле торцов электродов увеличивается быстрее, чем посередине. В целом скорость расширения канала становится больше. Вследствие этого быстрее увеличивается объём канала и уменьшается давление плазмы (рис. 11, 1 - электрод 200 мм, 2 - 10 мм).

После этого момента времени (« 10 мкс) форма канала разряда начинает отличаться от прямого цилиндра (рис. 10).

В результате произошедших изменений сопротивление канала начинает увеличиваться быстрее при электродах радиусом 10 мм, чем при электродах большого радиуса (рис. 12, 1 - электрод 200 мм, 2 - 10 мм). Из формулы (2) следует, что давление плазмы уменьшается в большей степени, чем увеличивается объём канала разряда. Это соответствует адиабатическому процессу расширения плазмы с показателем адиабаты, большим единицы (у = 1,26), и данным работ [15, 16].

Рис. 9. Поле давления через 10 мкс от начала разряда. Rl-, мм 1 1

Рис. 10. Профили канала через 15 мкс после начала разряда (1 и 2), 24 мкс (3 и 4), 35 мкс (5 и 6). Кривые с чётными номерами получены, когда радиус электрода равен 200 мм, а нечётные - 10 мм.

Рь. МПа

Rtk, Ом 0,20] 0,150.100,050,00 -I

_ _, — —

10 20 30 40 t, мкс

Рис. 11. Давление в канале.

ВЫВОДЫ

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

ЛИТЕРАТУРА

1. Гулый Г. А. Научные основы разрядно-импульсных технологий. Киев: Наукова думка, 1990. 208 с.

2. Golovashchenko S.F., Bessonov N.M., Ilinich A. M. Two-step Method of Forming Complex Shapes from Sheet Metal. J Mater Process Tech. 2011, 211, 875-885.

3.

4.

5.

6.

7.

0 10 20 30 40 t, мкс

Рис. 12. Сопротивление канала.

Rohatgi A., Stephens E.V., Soulami A., Davies R.W., Smith M.T. Experimental Characterization of Sheet Metal Deformation during Electro-hydraulic Forming. J Mater Process Tech. 2011, 211, 1824-1833.

Rohatgi A., Stephens E.V., Soulami A., Davies R.W., Smith M.T., Soulami A., Ahzi S. Electrohydraulic Forming of Sheet Metals: Free-forming vs. Coni-caldie Forming. J Mater Process Tech. 2012, 212, 1070-1079.

Кривицкий Е.В., Шамко В.В. Переходные процессы при высоковольтном разряде в воде. Киев: На-укова думка, 1979. 208 с.

Косенков В.М. Влияние длины канала высоковольтного разряда в воде на эффективность пластического деформирования цилиндрической оболочки. ЖТФ. 2011, 81(10), 133-139.

Dubovenko K.V. Allowance for the Interaction between the Underwater Electric Discharge Channel Plasma and the Shock Wave Reflected from the

Chamber's Wall. Surf Eng Appl Electrochem. 2013, 49(1), 28-35.

8. Блащенко А.Д., Половинка В.Д., Тертилов Р.В. Двухконтурный генератор импульсных токов с регулируемым включением контуров. ЭОМ, 2013, 49(1), 97-101.

9. Барбашова Г.А., Иванов А.В. Расчёт двумерного нестационарного течения жидкости при подводном электрическом взрыве проволочки конечной длины. Гидромеханика. 1986, 53, 16-19.

10. Barbashova G.A. Electric Circuit Inductance Dependence of a Hydrodynamic Load on a Thrown Plate during High-speed Impact Welding of Two Plates. Surf Eng Appl Electrochem. 2012, 48(3), 260-263.

11. Барбашова Г.А., Косенков В.М., Цуркин В.Н. О влиянии гидростатического давления на гидродинамические процессы при электроразряде. Сб.

науч. тр. ИИПТНАН Украины. 2000, 3, 36-42.

12. Косенков В.М., Каменская Л.А. Расчет расширения канала электрического разряда в жидкости, описываемой в потенциальном приближении. Акустический вестник. 2001, (2), 47-53.

13. Годунов С.К., Забродин А.В., Иванов М.Я., Край-ко А.Н., Прокопов Г.П. Численное решение многомерных задач газовой динамики. М.: Наука, 1976. 400 с.

14. Калиткин Н.Н. Численные методы. М.: Наука. 1978. 512 с.

15. Kosenkov V.M., Bychkov V.M., Zhekul V.G., Poklonov S.G. Electric Discharges of Big Length in Water and their Effect on the Energy of the Plastic Deformation of a Cylindrical Shell. Surf Eng Appl Electrochem. 2013, 49(3), 241-246.

16. Косенков В.М., Бычков В.М., Жекул В.Г., Поклонов С.Г. Экспериментальное исследование влияния электрических разрядов большой длины в воде на энергию деформации цилиндрической оболочки. Письма в ЖТФ. 2012, 38(9), 103-110.

Поступила 06.12.13

Summary

A mathematical model of the electrodynamic processes in the discharge circuit and the discharge chamber of an electro hydraulic assembly is developed. In addition, the algorithm of the solution of the obtained system of equations on the base of numerical methods is proposed. The simulation of the discharge channel expansion in water is realized in a geometrically nonlinear model, but on the finite difference scheme for a stationary grid. Both the mathematical model and the algorithm developed have been tested after the example of similar problems solution a well as by the comparison with the experimental data. The influence of the finite diameter of the electrodes on transient processes in the discharge circuit and the dynamics of the discharge channel expansion have been analyzed.

Keywords: discharge in water, discharge circuit, mathematical model, transient processes.