CC BY
47
ЭЛЕКТРИЧЕСКИЕ ПРОЦЕССЫ В ТЕХНИКЕ И ХИМИИ
КОНДЕНСАЦИЯ ПАРОГАЗОВОЙ СМЕСИ В УСЛОВИЯХ КОРОННОГО РАЗРЯДА
М.К. Болога*, Ф. П. Гросу**, А. А. Поликарпов*, О.В. Моторин*
*Институт прикладной физики АН Молдовы, ул. Академией, 5, г. Кишинев, MD-2028, Республика Молдова, [email protected] **Государственный аграрный университет Молдовы, ул. Мирчешть, 44, г. Кишинев, MD-2049, Республика Молдова, [email protected]
Введение
Вопросы влияния электрических полей на процессы конденсации достаточно подробно представлены в работе [1-3]. Основное внимание уделяется влиянию поля на границу раздела конденсатной пленки, ее деформацию и разрушение. Физическая трактовка явлений осуществлялась
на представлениях идеальных диэлектриков (силы f ~ E2 • Уб) [4] без учета возможностей появления
в среде свободных электрических зарядов (р^0). Такой подход можно считать оправданным, поскольку рассматривались высокоомные жидкие диэлектрики типа фреона и практически однородные поля. Однако отсутствуют данные о влиянии электрических полей на конденсацию паров сравнительно электропроводящих жидкостей, например воды, когда следует считаться с ее электропроводящими свойствами и возможностью электрической зарядки среды в целом. Этот пробел в определенной степени восполнен в работе [5], в которой приведены результаты экспериментальных исследований влияния электрического поля на процесс конденсации водяного пара из паровоздушной смеси.
Опыты проводились на экспериментальной установке, представляющей собой конденсатор пара в виде вертикальной трубы, вдоль оси которой натянута медная проволока с насечками в виде колец, на ее лаковой изоляции среднее расстояние между насечками составляло 5 мм. Насечки наносились специально, чтобы вызвать коронный разряд, поскольку в этом случае следует ожидать максимальный эффект воздействия поля на конденсацию. На проволоку подавалось постоянное высокое напряжение, а труба служила заземленным электродом. С нижнего конца трубы подавались пар и воздух, которые, смешиваясь, образовывали рабочую среду, из которой конденсировался пар. Конденсат стекал в мерный стакан, находящийся под трубой. Снаружи труба-конденсатор охлаждалась потоком воздуха от вентилятора.
Результаты исследований обрабатывались в виде графических зависимостей относительной массы конденсата от массовой концентрации воздуха GE/G0 = /1(с) при различных силах тока I и от силы тока при различных скоростях и внешнего обдува GE/G0 = /2(0. По этим двум зависимостям устанавливалась и третья - GE/G0 = /3(и). Опыты показали, что с ростом концентрации воздуха относительная масса конденсата уменьшается, а увеличение силы тока и скорости обдува трубы приводят к усилению процесса конденсации. Эти результаты, вполне понятные с физической точки зрения, неудобны в целях инженерных расчетов. Поэтому представлялось необходимым их критериальное обобщение, что и предопределило основную цель предпринятой работы.
Частные зависимости
Известно, что обобщенные зависимости ищут в виде степенных функций от обобщенных параметров, называемых критериями подобия [6]. Показатели этих зависимостей можно находить по частным зависимостям в логарифмических координатах:
, Ge lg_E
G0
lg / C); f ~ c2; I, и = const;
igf2(I); f2 ~Im; иC2 =const; (1)
lg /3(u); /3 ~ up; 1, C2 = const© Болога М.К. , Гросу Ф. П., Поликарпов А.А., Моторин О.В., Электронная обработка материалов, 2011, 47(4),
118-122.
Поскольку /1 (с) является убывающей функцией от концентрации воздуха, то целесообразно перейти к возрастающей функции/2(с2) от концентрации пара с2 =1-с, что и отражено в первой строке формул (1). С учетом этого по данным работы [5] на рис. 1 построены зависимости (1) в логарифмической сетке, в которой происходит спрямление кривых. По тангенсам углов наклона этих прямых (рис. 1) находим показатели степеней: п = 1,14; т = 0,15;р = 0,42.
а) ОъЮо = /¡Ы
б) ОЕ/О0 = /2(1)
в) ОъЮо = /зОО
Рис. 1. Зависимости /1(с2),/2(1) и/3(у)
Рис. 2. Схема насечек
Критерии подобия
Для построения обобщенной зависимости необходимы критерии подобия, соответствующие каждой из независимых переменных системы (1). Концентрация пара с2, будучи безразмерной величиной, сама является критерием подобия. Для скорости обдува критерием подобия, естественно, будет число Рейнольдса:
ае=и , (2)
V
где У0 - кинематическая вязкость воздуха комнатной температуры. Охлаждению воздуха
способствует и естественная конвекция вокруг трубы, однако этот эффект мал по сравнению с вынужденным обдувом, и он учитывается массой конденсата в отсутствие поля.
Наиболее сложным является учет влияния электрического поля на конденсацию. Сначала необходимо выяснить физическую сущность механизма взаимодействия поля и среды. В работе [4] данное взаимодействие трактуется с позиций идеальных диэлектриков, в частности путем учета деформации границы раздела пар-жидкость поляризационными силами. Однако, не оспаривая значимость такого эффекта, существеннее является коронный разряд [7] с вытекающими отсюда последствиями униполярной зарядки среды и явления электрического ветра. Оба эти фактора способствуют турбулизации паровоздушной смеси и интенсификации как процесса теплообмена, так и массопереноса. Веским аргументом в пользу вышесказанного предположения служит тот экспериментальный факт, что существенное влияние электрического поля на эти процессы наблюдается одновременно с разрядным током в высоковольтной цепи, причем этот ток однозначно коронной природы, поскольку возникает от областей очень малых размеров («царапин»). Следовательно для характеристики воздействия электрического поля на массо- и теплоперенос необходим критерий, учитывающий интенсивность электрического ветра. Такой критерий установлен нами ранее [8] и имеет вид
' Ж3 У
кц2
(з)
где j - плотность тока коронного разряда; у - плотность среды (смеси газов); п - коэффициент динамической вязкости; I - характерный размер коронирующей области; к - коэффициент подвижности ионов знака коронирующего электрода.
Как отмечено в [8], показатель степени q заключен в пределах 1>д>0,5, причем убывает от единицы до 0,5 по мере развития электрического ветра от ламинарного режима до турбулентного. Выражение (3), названное впоследствии «электрическим» числом Рейнольдса [9], успешно применено для обработки экспериментальных данных по теплоотдаче коронирующей проволоки к различным газам (воздух, углекислый газ, аргон, гелий) при различных давлениях (от 0,1 до 20 ат). Попытаемся использовать это число и для обобщения данных по конденсации в условиях коронного
разряда, поскольку физическая сущность рассматриваемых взаимодействий во многом едина. Конкретизируя параметры, входящие в формулу (3), прежде всего, заметим, что плотность тока будет зависеть, помимо электрического напряжения, и от геометрических факторов. На рис. 2 приведена схема насечек в виде заштрихованных колец.
Поскольку именно от них исходит коронный разряд, то ясно, что характерным размером в формуле (3) должна служить ширина насечки I = 5, которую в расчетах примем порядка ширины острия лезвия (10 ук). При этом плотность электрического тока определится цепочкой равенств:
7' = V1 = N; ^ = пй°5; м=7^ ^7' = %с10ьь' (4)
Паровоздушная смесь по критерию (3) характеризуется термодинамическим параметром, равным отношению у/п2 , который, как можно показать в рамках кинетической теории идеального газа, является величиной, зависящей от давления р, температуры Т и размеров молекул. Действительно, полагая смесь идеальным газом, для плотности будем иметь
у = -^. М, (5)
ЯТ
где М- молярная масса смеси.
С другой стороны, вязкость идеального газа дается формулами:
1 _ „ 1 3ЯТ 1 1 ЯТМ
п = —• У • и •А = —• Шс, П • .-----------р=--- =----;----. -------
1 з ' з 0 У м 42 • л^2 • п пй*2 ыА V 6
где й* - диаметр молекулы; т0 - ее масса; остальные обозначения общепринятые. С учетом этих равенств получим
у (~Л2 АТ Л2
— = 6 р • п2
пй*ЫА
ЯТ
параметр который при данных р и Т зависит лишь от размеров молекул. Для воздуха и водяного пара диаметры молекул приблизительно одинаковы, поэтому дробь в левой части можно считать постоянной и отнести лишь к одной компоненте, удобнее - к воздуху.
Вместе с тем на практике имеем дело с реальными газовыми смесями, для которых вязкость сложным образом зависит как от частных коэффициентов вязкости, так и других физических параметров [10]. Поэтому в расчетах желательно применение более адекватных величин из таблиц экспериментальных данных для смесей. Однако для плотности, как показывают оценки, с достаточной точностью приемлема формула (5), в которой молярная масса определится
общепринятым выражением:
т, + т2
М =----------1---2-------- , (6)
(т1 / М!) + (т2 1М2)
где М\ и М2 - молярные массы; индексы «1» и «2» относятся к воздуху и пару соответственно. Приняв во внимание определения концентраций с1 =т^т; с2 =т2/т; т=т1+т2, получим
Т = Р ,_М^ = _П_. т • = М-!. (7)
' ЯТ 1+у '•с, 1 + у ••с1 ' М 2
Как и следовало ожидать, плотность смеси падает с ростом плотности пара.
Коэффициенты вязкости п [10] и подвижности к, как уже было отмечено, сложно зависят от значений этих коэффициентов для отдельных компонентов, поэтому в первом приближении определим их как среднестатистические:
п = с1П: + С2п = П:(1 -п'с2); п' =1 -(VП:); (8)
к = с1к1 + с2к2 = к1(1 - к'с2); к' = 1 - (к2 / к1).
Эти формулы применимы лишь в силу близости искомых коэффициентов пи к для воздуха и пара. Из них следует, что наличие воздуха, у которого как вязкость, так и подвижность больше, чем у пара,
приводит к увеличению обоих этих коэффициентов, что способствует ухудшению процесса конденсации.
Подставив последнюю формулу (4) в (3), с учетом (7), (8) получим
ЯеЕ =
І'^ Л Ч " 1 "
V пй0 ^к1П1 у _ (1 + У ' С2)(1 - к ' С2)(1 - П ' с2)2 _
(9)
Обобщенная зависимость
Таким образом, имеются три критерия подобия: с2; ЯеЕ и Яе. Ищем обобщенную зависимость в стандартном степенном виде:
ОЕ/О0 = а-сп2 -ЯетЕ-Яер . (10)
Подставив в данное выражение формулу (9), найдем
—
-Е = а-е’Е-Яе"-—
(1 + у' С2Х1 - к' С2)(1 - п^)2
Яе р,
(11)
где выделен отдельно критерий, относящийся только к воздуху:
ЯеЕ1 —
1-І|-52 У1
пй0 ІУіП!2
(12)
Так как выражение в квадратных скобках (11) для данного состава смеси есть функция только от концентрации пара с2, то в первоначальном упрощенном варианте обобщения данных опыта эту функцию объединим с множителем сП и в этом же виде будем искать окончательную зависимость от
с2. Тогда в (10) вместо зависимости от (9) останется зависимость от (12), причем для простоты за показателем степени тд сохраним обозначение т, и тогда искомая критериальная зависимость примет вид
^ = ас2п Яе^-Яер = ас’Пп —0
(ий V
пй0 Ьк{П
\у0 У
(13)
где постоянная а подлежит нахождению на основе опытных данных. Так как ЯеЕ~ I, Яе~и, то е и р считаем найденными (см. выше). Тогда для уравнения (13) находим
Г' IГ' „ „1,14 т> „0,15 т> 0,42
СЕ 1С() — а • с2 • Яея1 • Яе Коэффициент а можно найти как среднее арифметическое выражения
а — (ОБ/О0)(с214^еЕ1^Яе0 42 )-1.
(14)
(15)
Ч
Ч
Обработка экспериментальных данных приводит к значению а — 0,018 и окончательной обобщенной зависимости:
0Е/00 * 0,018 • с^Яе^-Яе0,42. (16)
Следует заметить, что малое значение показателя степени «электрического» числа Рейнольдса свидетельствует об интенсивном перемешивании среды внешним электрическим полем, ибо, как отмечалось, показатель степени этого критерия содержит еще один показатель (д^-0,5), уменьшающийся при развитии электроконвекции.
Данную работу следует рассматривать как одну из первых попыток обработки экспериментальных данных по конденсации с помощью концепций коронного разряда. Предстоят
уточнения роли параметров среды в процессе конденсации при наличии поля. Результаты позволяют сделать общее заключение — о состоятельности данных концепций.
ЛИТЕРАТУРА
1. Болога М.К., Смирнов Г.Ф., Дидковский А.Б., Климов С.М. Теплообмен при кипении и конденсации в электрическом поле. Кишинев: Штиинца, 1987. 320 с.
2. Gidwani A., Molki M., Ohadi M.M. "EHD-enhanced Condensation of Alternative Refrigerants in Smooth and Corrugated tubes," HVAC&R Research, 8(3), p. 219-237 (2002).
3. Enhancement of Condensation Heat Transfer on a Finned Tube Using an Electric Field Ren Cai C. Chu, Shigefumi Nishio, Ichiro Tanasawa, Journal of Enhanced Heat Transfer, v. 8, issue 4, 2001, p. 215-229.
4. Болога М.К., Коровкин В. П., Савин И. К. Двухфазные системы жидкость-пар в электрическом поле. Кишинев: Штиинца, 1992. 243 с.
5. Bologa M.K., Grosu F.P., Polikarpov A.A. and Motorin O.V. Condensation of a gas-vapor mixture in an electric field. Surface Engineering and Applied Electrochemitry. 2009, 45(2), 125-127.
6. Incropera F.P., De Witt D.P., Bergman T.L., Lavine A.S. “Fundamentals of Heat and Mass Transfer”, John Wiley & Sons, 6th Edition, 2007.
7. Капцов Н.А. Электрические явления в газах и вакууме. М.: Гостехиздат, 1950.
8. Болога М.К., Гросу Ф. П., Кожухарь И. А. Электроконвекция и теплообмен. Кишинев: Штиинца, 1977. 320 с.
9. Grosu F.P. and Bologa M.K. Electroisothermal convection and its role in the process of heat exchange. Surface Engineering and Applied Electrochemitry. 2008, 44(3), 187-194.
10. Голубев И. Ф. Вязкость газов и газовых смесей. М.: Физматгиз, 1959. 375 с.
Поступила 21.12.10
Summary
The experimental data for condensation of vapor from a gas-vapor mixture in the presence of an electric field are discussed. It is assumed that a corona discharge or, to be more precise, its accompanying phenomena such as unipolar charging of the medium and the electric wind, which produces intense mixing both of the medium and condensate film, is the main factor influencing the intensification of condensation in an electric field. Criteria of similarity of the process, among which the main attention is drawn to the "electric" Reynolds number, which characterizes the interaction of the electric field with the medium, are considered. We also discuss the parameters of vapor-air mixtures, such as density, viscosity and ion mobility. A simplified method of generalization of experimental data on condensation in the presence of an electric field, based on the concepts of corona discharge, is proposed.
