ЭЛЕКТРИЧЕСКИЕ ПРОЦЕССЫ В ТЕХНИКЕ И ХИМИИ
М.К. Болога , Ф.П. Гросу , И В. Кожевников
ОСОБЕННОСТИ ЭЛЕКТРОГИДРОДИНАМИЧЕСКИХ ТЕЧЕНИЙ В МНОГОЭЛЕКТРОДНЫХ СИСТЕМАХ
* Институт прикладной физики АНМ, ул. Академией, 5, г. Кишинев, MD-2028, Республика Молдова **Государственный аграрный университет Молдовы, ул. Мирчешть, 44, г. Кишинев, MD-2049, Республика Молдова
Введение
Одно из наиболее ярких в научном и прикладном аспекте проявлений взаимодействия электрических полей с диэлектрическими средами - это электрогидродинамические (ЭГД) течения, управление и оптимизация которых сопряжены с немалыми трудностями, и их преодоление, безусловно, заслуживает внимания.
С этих позиций в работе приводятся результаты экспериментальных и теоретических исследований ЭГД течений в трех- и четырехэлектродных системах применительно к реализации ЭГД насосов диэлектрических жидкостей.
1. Результаты экспериментальных исследований. Наибольший интерес с практической точки зрения, а именно организации сквозной прокачки диэлектрической жидкости - представляет исследование электрогидродинамических течений в узких каналах при малых расстояниях между электродами, обеспечивающих приемлемые напряжения для развития ЭГД явлений. В этих условиях применение пузырьков воздуха [1] для визуализации течений невозможно ввиду их влияния на распределение электрического поля и объемного заряда в межэлектродном пространстве (МЭП).
Структуры ЭГД течения исследовали с помощью теневого метода; экспериментальная установка и методика исследований описаны в [2]. Электродами служили лезвие (э - эмиттер) и две проволоки диаметром 1 мм (к - коллектор), которые расположены параллельно друг другу (рис. 1), на одинаковом расстоянии, симметрично относительно кромки лезвия, аналогично системе игла-кольцо. Электроды установлены в канале шириной 5j = 6,6 мм, являющемся участком замкнутого трубопровода. Одновременно с измерением максимальной скорости сквозного течения жидкости в трубопроводе вне поля фотокинокамерой проводилась регистрация структуры ЭГД течения, визуализированного теневым методом. В качестве рабочей среды использовалось очищенное трансформаторное масло.
а б в г
Рис. 1. Структуры ЭГД течений при различных напряжениях, кВ: а - 10, б - 15, в (d1 = 9,5 мм). Рабочая среда- трансформаторное масло
20, г - 20,
© Болога М.К., Гросу Ф.П., Кожевников И.В., Электронная обработка материалов, 2007, № 6, С.29-34.
29
Теневые картины структуры ЭГД течений при симметричном расположении центральной струи относительно проволок коллектора и при геометрических параметрах S = 5 мм (зазор между образующими проводов), d = 1,5 мм (расстояние между эмиттером и коллектором) приведены на рис. 1, 2. По мере повышения напряжения между электродами до - 20 кВ максимальная скорость течения в трубопроводе сначала возрастает с 0,5 см/с (при U = - 15 кВ) до 0,7 см/с (при U = - 20 кВ), а затем снижается практически до 0 (при U = - 25 кВ). Наблюдается усиление флуктуирования центральной струи между проводами коллектора. Видимо, это одна из причин понижения интенсивности сквозной прокачки. Другая причина состоит в том, что заряженная эмиттером диэлектрическая жидкость не успевает полностью разрядиться на коллекторе (ввиду ограниченности поверхности). В результате за МЭП возникает сила кулоновского взаимодействия между заряженной средой и электродами коллектора, направленная против потока жидкости и возрастающая с увеличением напряжения, как следствие возникают возвратные течения в виде вихрей, направленных к коллектору (рис. 1,а,б,в).
Рис. 2. Структуры ЭГД течений при различных расстояниях d1 между четвертым электродом и
коллекторами, мм: а -щ б - 9,5; в - 4,5; г - 1,5; U = -25 кВ
Об этом же свидетельствуют эксперименты с дополнительным заземленным электродом 01 мм (ДЭ, рис. 1,г), расположенным симметрично и на расстоянии d1 от проводов коллектора, который влияет как на структуру ЭГД течения (рис. 1,в,г), так и на интенсивность сквозной прокачки. Дополнительный электрод способствует “вытягиванию” заряженной среды из приколлекторной области, что хорошо заметно при сравнении структуры ЭГД течений (рис. 1,в и г). По мере приближения центрального дополнительного электрода к боковым (уменьшения расстояния d\) наблюдаются стабилизация центральной струи (рис. 2) и вырождение вихрей за коллектором
Рис. 3. Зависимость максимальной скорости течения в трубопроводе от напряжения при различных расстояниях между электродами, мм: 1 - d = 1,7, S = 5,0; d1: 1 - без дополнительного электрода, 2 -1,5, 3 - 4,5, 4 - 9,5; 5 - d = 3, Si = 5,0, d1 = 1,5; d = 1,7, S = 1,8, dj: 6 - без дополнительного
электрода, 7 - 6,9, 8 - 2,5, 0,8; и= 3,9-10-12 См-м-1
30
Сужение канала 5i и одновременно зазора S между проволоками коллектора в отсутствие дополнительного электрода приводит к уменьшению амплитуды колебаний центральной струи, интенсивности возвратных течений как в МЭП , так и за ним, что положительно сказывается на скорости течения среды (рис. 3, кривые 1 и 6). В то же время на коллекторах не происходит полной нейтрализации заряженной среды. Это подтверждается данными по воздействию четвертого электрода на интенсивность сквозной прокачки (рис. 3, кривая 7), которая практически остается постоянной при изменении расстояния между коллекторами и дополнительным электродом в пределах 0,8-2,5 мм (рис. 3, кривая 8). Видимо, по мере приближения четвертого электрода к коллекторам (d1 < 2,5 мм) увеличивается сопротивление, оказываемое им потоку жидкости, что компенсирует повышение скорости прокачки.
2. Теоретическая интерпретация экспериментальных результатов. Адекватными представляются две математические модели для описания экспериментально наблюдаемого ЭГД течения. Одна из них - модель диффузорного течения типа Галлея [3], другая - модель плоской затопленной струи, аналогичная рассмотренной Ландау [3]. Хотя математические трудности на пути решения обеих задач преодолимы, результаты оказываются весьма громоздкими и мало обозримыми, поэтому на начальном этапе ограничимся более прозрачной с физической точки зрения моделью.
В соответствии с представленной экспериментальной картиной ЭГД течения в системе электродов O, Ob O2, O3 (рис. 4) плоская струя от лезвия О пробивается вверх, сквозь толщу жидкости. При этом движущей является кулоновская сила объемной плотностью
f = PE, (1)
в области заземленных коллекторов О1 и О2 возникают два примерно симметричных вихря (рис. 1, 2).
На начальных стадиях течения, то есть при малых напряжениях, вихри незначительны по размерам и течение диффузорное радиальное через всю дугу окружности радиуса R, измеряемую центральным углом 2а (см. рис. 4). По мере развития течения, то есть увеличения напряжения, вместе со скоростью течения увеличиваются и размеры вихрей “О1” , “О2”, заслоняя тем самым свободный просвет для сквозной прокачки жидкости.
Итак, усиление электрической прокачки жидкости приводит к двум противоборствующим тенденциям: увеличению расхода за счет роста скорости и одновременно уменьшению расхода из-за сокращения свободного просвета канала, измеряемого центральным углом 20 (рис. 4). Так как при начальных напряжениях расход Q ~ 0, а при больших также стремится к нулю, ввиду полного перекрытия сквозного прохода жидкости возросшими вихрями, то ясно, что должен существовать оптимальный (максимальный) расход при некотором промежуточном режиме (напряжении) работы
31
ЭГД насоса. Из простых физических соображений попытаемся дать количественное объяснение этим представлениям, определяя первоначально расход по формуле
Q = у [ dz [ и(R) • Rd9 = 2ybRu • 9, (2)
J 0 J—9
где у — плотность среды, b - длина кромки лезвия; и — средняя на дуге 29 скорость. Для ее определения вспомним, что “электрический ветер” от острий в газах (при хорошем разряде) и в жидкостях имеет по существу одну и ту же природу и движущую силу [4—6], определенную формулами
j = кРЕ ^ f = j / К (3)
где к — подвижность носителей зарядов знака электрода-лезвия, j определяемая равенством
I
j =
плотность электрического тока,
(4)
3 • п-d •b
Здесь полный ток I отнесен к площадям трех проводников-коллекторов диаметром d0 и длиной b. Ранее нами показано [5], что среднюю скорость электрического ветра можно определить формулой
и=с
(г, Л ( ji3 л
г
V1 2 J
V
' г
(5)
l = dr
0 ’
куу J
где l - характерный размер, в качестве которого, очевидно, следует принять v — кинематическая вязкость, 1 > m > 1/2и убывает по мере перехода ламинарного режима течения в турбулентный [5]; C (Г1 / Г2) — коэффициент пропорциональности, зависящий от симплексов
подобия, порождаемых граничными условиями.
Угол 9 можно определить формулой
9 = а — 5 = а — 2r0/ R, (6)
где r0 — радиус боковых вихрей согласно рис. 4, для оценки которого полагаем, что вихрь круговой и
частицы, участвующие в вихревом движении, подвержены действию лишь центростремительных сил, являющихся одновременно и кулоновскими, поэтому априори предположим, что эти силы равны:
уи
= РЕ = j / к
(7)
откуда
2г = Л,0
2к уи
(8)
r
0
где
Подставляем (4) в (5), затем (5) и (6) в (2):
(г Л ( тЛ‘
v
Q = 2YbR I d C
Г
И 2
I
VI* J
a —
2 (г Л ( г Л2т
0'w v с2 £i_ 1
6nbd0ку v
RI d2
г
И 2 J
VI* J
или
Q = A (agm — B ■ g3m—1),
V
f г Л
A _ 2ybR—C d0
г
И 2 J
_____I ; T _ ЗпЬкуу2
; g = I ’ I* _ d2
1 * 14a
6пЬку•v
2 С T ^
B _ « ’ C2
RI*d0
г
V1 2 J
Исследуем (9) на экстремум:
dQ
dg
= Aamg"
1 B ( 3m — 1)
am
g
(9)
(10)
32
Приравнивая нулю это выражение, находим экстремальное значение безразмерной силы тока:
£о =
am
B (3m -1)
i
2 m-1
(11)
Находим вторую производную
dg
Подставив сюда (11), получим
Щ = Aam (m - 1)g~ (1 - B (3m ~(1)(3”’~ 2)
dg am (m -1)
2m-1
g
d 2Q dg2
= -Aam (2m -1) gn
Поскольку 0,5 < m < 1, то d2Q / dg < 0 и функция (9), то есть расход, имеет максимум при
g = g0, даваемом выражением (11), причем
Qmax
aA (2m -1) 3m -1
am
B (3m -1)
m
2 m-1
(12)
§ 5 6
Приведем численный расчет Qmax, полагая a = arctg----= arctg—-— = 1,21 рад,
2d1 2 * 1,5
b = 10 *10~3 м; R = 3,18 *10-3 м; к = 10-8 м2/В * с; d = 10-3 м; v = 22 *10-6 м2/с; у = 0,9 *103 кг/м3;
A = C(Г1/Г2)* 1,26*10-3 кг/с; I* = 0,4*10-9 A; B = C(Г1 /Г2)*1,7*10_1; m = 0,7 е(0,5;1,0). Получим
Qr = C -0'75 • 2,86.10-3 (кг/с).
В экспериментах и ~ 1 см/c=10"2 м/с; S ~ 1 см2 = 10-4 м2 и Q^^x ~ yuS ~ 10-3 кг/с, то есть
по порядку величины, то же самое и эмпирический коэффициент пропорциональности ~1 -необходимый и верный признак состоятельности теоретических концепций.
Из (9) следует, что если сила тока достигает критической величины
gкр |п
a 12 m-1
то работа насоса полностью прекращается.
Подставив численные данные в (13), найдем
g^ ~ 5/C"5,
то есть критические токи такого же порядка, как и I*, которые назовем “приведенными”. Связь между g0 и g^ дается формулой
(13)
g кр g0
3m -1 ]2 m-1
m
(14)
которое при m ~ 0,7 дает ~ 3-кратное значение g0.
Следует заметить, что формулы (9)-(14) носят степенной характер с переменным показателем степени, меняющимся в узких пределах (0,5-1,0) и в сочетании с одним-единственным эмпирическим коэффициентом пропорциональности делают их удобными и гибкими для расчетных целей.
Таким образом, показано, что неустойчивость струи жидкости, направленной от центрального электрода (лезвия), и вихреобразование в трехэлектродной системе можно частично подавить, используя дополнительный электрод, установленный за МЭП, или уменьшив зазор между боковыми электродами, тем самым снизив их расфокусирующее действие на центральную струю течения, что в целом приводит к увеличению интенсивности сквозной прокачки диэлектрической жидкости. Экспериментально обнаружено и теоретически подтверждено наличие максимума скорости прокачки в зависимости от силы тока (напряжения) в ЭГД насосе. На основе простой физической модели ЭГД
33
течения получены расчетные формулы для определения расхода жидкости через ЭГД насос и его оптимальных режимов. Теоретические результаты согласуются с экспериментальными.
Статья выполнена при поддержке гранта РФФИ 06-08-91020-Мол/А и Академии наук Молдовы 06.13CRF, а также Институционного проекта 06.408.003F.
ЛИТЕРАТУРА
1. Стишков Ю.К., Остапенко А.А. Электрогидродинамические течения в жидких диэлектриках. Л., 1989, 176 с.
2. Болога М.К., Кожевников И.В. Влияние изоляционного покрытия на процессы зарядообмена между электродом и жидкостью // Электронная обработка материалов. 2006. № 1. С. 32-36.
3. Ландау Л.Д., Лифшиц Е.М. Гидродинамика.
4. КапцовН.А. Электрические явления в газах и вакууме. М., Гостехиздат, 1950. 836 с.
5. Болога М.К., Гросу Ф.П., Кожухарь И.А. Электроконвекция и теплообмен. Кишинев: Штиинца, 1977. 320 с.
6. Остроумов. Взаимодействие электрических и гидродинамических полей. М.: Наука, 1979.
Поступила 14.08.07
Summary
Electrohydrodynamic flow in the system “needle-two wires” have been investigated using the Shadow method with simultaneous recording of transversal pumping intensity. Not-stable electrohydrodynamic flow appears as a fluctuation of central jet (flow from the blade) between electrodes that results in decreasing of efficiency of three electrode system. The possibility of electrohydrodynamic flow stabilization by means of decreasing interelectrode gap and using of additional, fourth electrode, which is situated behind the first two electrodes, is shown. Theoretical model of the pumping process, considering the flow as a diffuser flow from the blade, limited by two lateral eddies near the upper electrode-collectors, is proposed. As electrohydrodynamic flow develops, the average speed of diffuser flow as well as eddy dimensions are increasing. These two factors influence on the maximum flow rate of electrohydrodynamic pump. The maximum flow rate was found theoretically and the optimum parameters of electrohydrodynamic pump were established.
34