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

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

И. Н. Карпович

Национальный университет водного хозяйства и природопользования, ул. Соборная, 11, г. Ровно, UA - 33028, Украина, e-mail: karpivan@ukr. net

Выполнены расчеты потоков пара и пленочной жидкости под влиянием градиентов напряженности неоднородного электрического поля и давления в зависимости от радиуса и конусности капилляров для воды и этиленгликоля.

Ключевые слова: конический капилляр, жидкая пленка, неоднородное электрическое поле, массоперенос, пористое тело.

УДК 532.685:537.228

ВВЕДЕНИЕ

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

В теоретических исследованиях [1-4] неоднородное электрическое поле (НЭП) рассматривается как эффективный метод управления тепло-массопереносными процессами при взаимодействии различных поверхностей с жидкостями и газами. Особая роль отводится внешним полям высокой напряженности и неоднородности, в частности, полю коронного разряда (ПКР), способного заметно интенсифицировать процессы массообмена в дисперсных системах. Основные выводы развитой теории находятся в согласии с экспериментальными исследованиями влияния ПКР на испарение жидкостей из кварцевых капилляров, которые являются моделью пористых структур. Так, в работе [5] показано, что применение ПКР позволяет повысить скорость испарения воды из капилляров в 3-5 раз. Это объясняется тем, что внутри капилляра, помещенного в ПКР, неоднородность поля способствует добавлению к диффузионному потоку молекул пара электроконвективного потока. Вне капилляра разрядное поле действует как осушитель, понижая давление пара над устьем капилляра. Этот

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

Влияние ПКР на процесс испарения полярных жидкостей проявляется особенно заметно при высоких значениях относительного давления пара (р/р > 0,80). В этих условиях на стенках капилляра за отступающим мениском жидкости образуется полимолекулярная пленка. Если градиент напряженности поля совпадает с потоком испарения, под воздействием силового поля [3, 6] эффективная толщина смачивающей пленки увеличивается за счет подпитки ее из основного столбика капиллярной жидкости. Вытягиваемая толстая пленка в сильно неоднородном поле становится неравновесной. Это может быть вызвано адсорбцией ионов на поверхности капиллярной жидкости в ПКР, которая способна существенно изменять свойства поверхности раздела фаз [4]. Возникающие флуктуации толщины пленки приводят к ее смыканию и образованию перемычек, которые защемляют пузырьки паровоздушной смеси. После образования пузырьков наблюдается их продвижение в глубь капилляра, заполненного жидкостью, то есть жидкость по смачивающей пленке, которая покрывает участок капилляра между менисками пузырька, перемещается в область большей неоднородности поля ближе к устью капилляра. Таким образом, эффект диспергирования жидкого столбика с последующей электромиграцией полярной жидкости по стенкам капилляра к поверхности испарения обеспечивает высокую степень интенсификации процесса испарения.

© Карпович И.Н., Электронная обработка материалов, 2016, 52(2), 82-87.

В сложной и разнообразной геометрии поро-вого пространства реальных капиллярно-пористых материалов преобладающую часть занимают конические и тупиковые поры [7]. Заслуживают внимания явления массообмена в условиях, когда поры заполнены жидкостью частично. Такие процессы имеют место, например, при увлажнении и осушении почвы, изготовлении и эксплуатации керамических изделий, пористых электродов, элементов тепловых труб, строительных материалов и конструкций. Механизм явлений, определяющих закономерности массо-обмена, с учетом как основных особенностей пористой среды, так и внешних воздействий, которые, в частности, могут создаваться силовым полем, изучен недостаточно. Поэтому цель настоящей работы - расчет потоков пара и пленочной жидкости под влиянием градиента напряженности внешнего НЭП и градиентов расклинивающего и капиллярного давлений в коническом капилляре.

ОПИСАНИЕ МОДЕЛИ

Для анализа массопереноса рассмотрим модель элемента порового пространства в виде конического капилляра, частично заполненного жидкостью (рис. 1а). Мениски смачивающей жидкости имеют различную кривизну и разделены пузырьком воздуха. Поверхность капилляра на участке переноса между менисками покрыта смачивающими пленками, толщина которых к(х) зависит от координаты х: к1 = к(0) и к2 = к(!), где 1>>г - длина пузырька. НЭП высокой неоднородности может создаваться системой электродов «игла-пластина» (рис. 1б). Вычисления двухфазных потоков под действием внешнего НЭП проведем с учетом условия локального равновесия, которое обеспечивается равенством химических потенциалов молекул в пленке и паре в каждом сечении капилляра [8]:

J = -1 т V U - I „ V E,

р = RT- ш р

где и - влагосодержание материала; Е - напряженность электрического поля; Х„, Хе - коэффициенты изотермической и электровлагодиффу-зионной массопроводности соответственно.

Если поры заполнены жидкостью частично, перенос вещества происходит в виде пара и жидкости. Движущей силой диффузионного переноса пара является градиент парциального давления пара Vp. Вязкое течение в пленке происходит под действием градиента капиллярного и расклинивающего давления VP. Используя уравнение Фика для диффузионного потока и уравнение гидродинамики тонких пленок [9] для пленочного потока, получим плотность потока массы в коническом капилляре без наложения поля:

777 БУс 2 р, к3

= Jl + J 2 =-^:гVpV р,

КТ 3 п г

где Б - коэффициент диффузии пара; с - концентрация молекул пара; р, - плотность жидкости; ^ - ее вязкость.

Для описания электроконвективных явлений используют приведенную ниже систему уравнений (4), предложенную в монографии [3]:

ЗУ 7 Vp ^ 2 7

— + (иУ)и = g ———^2и +

Зт

+

Р,

Vи = 0;

е0 Е2

2р,

^е, +

Р, Р, 1

2Р,

е oV

Зе,

ЗР,

р,

где Р - давление в пленке жидкости; Т, рц, р = р(х) - температура, давление насыщенного пара и давление пара соответственно; К - газовая постоянная; у- молярный объем жидкости.

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

VI = ; е0е

V7 + ^ = 0; Зт

I = Р еЬЕ + Реи- ПЧРе; Vx А = 0,

где и - средняя скорость движения жидкости; х - время; g - ускорение свободного падения тела; р е - объемная плотность заряда; в0 - диэлектрическая постоянная; в, - диэлектрическая проницаемость жидкости; / - плотность тока; Ь - коэффициент подвижности ионов в жидкости (газе).

Уравнение Навье-Стокса системы (4), кроме обычных факторов (давление, гравитация, вязкость), содержит и плотность электрических сил:

I = РеЕ - 0,5е0Е^е, + 0,5еoV Е2Р,

Зе, ЗР,

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

у

неполном насыщении (капилляр заполнен жидкостью частично) наложение НЭП вызывает электродиффузию пара в свободной от жидкости части капилляра и электромиграцию жидкости в пленке, которая покрывает стенки капилляра между менисками капиллярной жидкости. Плотность потока массы, вызванного действием НЭП, учитывая [2, 3], можно записать:

J = j3 + J = V E + 2 Р 1h fe

e 3 4

.

кТ 3 п г

где ре - дипольный момент молекулы; к - постоянная Больцмана.

Дифференцирование соотношения (1) по координате х позволяет найти зависимость УР от Ур и подставить в (3) вместо УР его значение, выраженное через Ур.

Тогда плотность общего потока массы в коническом капилляре при наличии внешнего НЭП в выбранной системе координат (рис. 1) запишется в виде:

7 Dvcw 2p,RTh3 Vp Dpc J =--Vp —--— + e

RT

2p^fe

3nr

3nvr

= J + J2 + J3 + J4.

—VE + kT

янной остается только сумма слагаемых в уравнении (7), в то время как величина каждого из них меняется по длине капилляра. Продифференцированное по координате х уравнение (7) решалось на ПК методом итераций. Для этого проводилась квазилинеаризация уравнения согласно [13] и находилось численное решение краевой задачи, то есть распределение р = р(х). Точность решения контролировалась путем проверки условия

J - Jn)/JJ<S,

В дальнейших расчетах и рассуждениях будем использовать модули векторных величин. Давление в пленке складывается из капиллярного Рк = - c/r, обусловленного ее кривизной, и расклинивающего П, обусловленного действием поверхностных сил [8]. Толщину жидкой пленки определяем на основании экспериментальных изотерм расклинивающего давления n(h), полученных в работе [10] для пленок воды и эти-ленгликоля на поверхности цилиндрических кварцевых капилляров. Применяя аппроксимацию [11] n(h) = A/h3 (A = const), удовлетворительно описывающую известные экспериментальные данные, для давления в пленке получим:

р = - О-П=- О - А

r r h

Задавая радиус капилляра

r

где a = (r2 - r1)// - конусность капилляра, можно выразить h(x) через P(x). Замена P(x) на p(x) позволит найти решение уравнения (7) в виде

р=р(х).

Краевые условия уравнения (7) с учетом [12] имеют вид:

Р1 = Р

При использовании условия локального равновесия, как было показано в работе [9], посто-

где n - номер итерации и 5 = 0,01 - заданная относительная погрешность.

Из-за различной кривизны менисков (r1 < r2), ограничивающих пузырек паровоздушной смеси в капилляре, градиент давления пара будет положительным (p1 < p2). Поэтому первый и второй члены уравнения (5) отрицательны: первый определяет диффузионный поток молекул пара, а второй - обратный поток жидкости в пленке под действием градиента давления, обусловленного конусностью капилляра. Третье слагаемое, которое соответствует электродиффузионному потоку дипольных молекул пара, и четвертое, характеризующее электромиграцию жидкости под действием НЭП, при VE > 0 будут положительными.

В проведенных вычислениях использованы табличные значения физических констант, которые отвечали заданной температуре. На участке l зависимость E(x) принималась линейной, то есть VE = const. При линейной зависимости диэлектрической проницаемости s1 от плотности диэлектрика р1 в изотермических условиях выражение (5) преобразуется к виду f = peE + 0,5е0(е1 -1)VE2. В связи с отсутствием корректных значений объемной плотности заряда ре в капиллярной жидкости расчеты выполнены при ре = 0. Учитывая эти условия, определялись как отдельные составляющие плотности потока, так и общий поток массы J.

Расчеты выполнены для двух полярных жидкостей с различной летучестью и вязкостью: воды и этиленгликоля. Зависимости J и отдельных составляющих потока от радиуса и конусности капилляров получены для значений VE от 3-107 до 3-108 В/м2 при среднем значении E = 106 В/м.

Для воды принято (при T = 293 К): D = 0,252-10-4 м2/с; pe = 6,2-10"30 Кл-м; с = 17,29-10-3 кг/м3; pj = 998 кг/м3; ц = 10-3 Н-с/м2; v = 18,06-10-3 м3/кмоль; S! = 78,3; с = 72,75-10-3 Н/м; ps = 2,38-103 Н/м2.

Для этиленгликоля: D = 0,52-10-5; pe = = 7,4-10-30; с = 3,5-10-3; pj = 1113; ц = 19-10-3; v = = 55,77-10-3; S! = 38,7; с = 46,1-10-3;ps = 0,068-103.

Г

%

НР

Рис. 1. Расчетная схема элемента порового пространства (а) и схема расположения капилляра (б) в неоднородном электрическом поле.

Отношение 1/г менялось от 10 (для больших г) до 1000 (для малых г). В этих условиях, как было показано авторами [14], в капиллярах радиусов г > 1 мкм могут оставаться устойчивыми толстые (100-300 нм) метастабильные пленки жидкости. При увеличении длины пузырьков толстые пленки становятся неустойчивыми, что приводит к резкому снижению скорости массопереноса в НЭП.

При неполном насыщении капиллярно-пористых тел основным механизмом массопере-носа является перемещение жидкости в пленках, которые покрывают стенки капилляров. В цилиндрических капиллярах радиусами 1-100 мкм пленочное течение может осуществляться под влиянием градиента силового поля [10], причем вклад потока пара в общий поток массы здесь не превышает 1%. В конических капиллярах, как известно, соотношение потоков выглядит по-другому: даже при отсутствии внешних воздействий имеет место поток 32 жидкости в пленках, обусловленный разностью капиллярных давлений менисков с различными радиусами. Вклад потока 32 в общий поток массы с увеличением конусности капилляров а возрастает, и при достаточно больших значениях а становится доминирующим [12].

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Принимая поток 32 основным (уравнение 7), влияние НЭП на массообмен в капиллярах с различной конусностью можно оценить отношением абсолютных величин 3 /1 32 |. На рис. 2 показаны зависимости этого отношения от конусности капилляра, рассчитанные для х = 1/2. Величина каждого из слагаемых (отдельных потоков) в уравнении (7) меняется по длине капилляра. Постоянным остается только суммарный поток 3.

Поток 32 всегда направлен в узкую часть капилляра, то есть для избранной нами конфигурации внешнего поля 32 < 0. С увеличением конусности а абсолютная величина потока 32 возрастает. При 3 /| 32 |= -1 основным механизмом переноса является поток в пленках под дей-

ствием капиллярных сил. Например, для воды при г1 = 1 мкм и а > 10-5 можно считать 3 « 32 (рис. 2).

Рис. 2. Зависимости 3 / | 32 | от конусности капилляров а

для воды (1,3) и этиленгликоля (2,4): 1 — г1 = 1 мкм, I = 0,5 мм; 2 — г1 = 10 мкм, I = 1 мм; 3 — г1 = 40 мкм, I = 1 мм; 4 — г1 = 50 мкм, I = 1 мм; Е = 106 В/м; УЕ = 108 В/м2.

С уменьшением конусности абсолютное значение потока 32 стремится к нулю, следовательно, 3 /1 32 х. В этих условиях вклад электромиграционного потока в пленках 34, вызванного действием НЭП, становится подавляющим - он в несколько раз превышает пленочный поток 32 под действием градиента капиллярных сил. Отсюда, в частности, следует, что существенного ускорения массопереноса можно достичь как раз для таких жидкостей, когда градиент капиллярного давления мал (плохое смачивание, низкое межфазное натяжение, большие радиусы пор). Полем можно удерживать и перемещать в капилляре несмачивающую его жидкость при условии, что отрицательные значения второго слагаемого компенсируются положительными значениями четвертого слагаемого в уравнении (7). В то же время, меняя знак поля, можно препятствовать вхождению смачивающей жидкости в капилляры или ускорять этот процесс. Таким образом, применение НЭП создает возможность управления скоростью капиллярного массообмена в нужном для практики направлении.

Из рис. 2 видно, что с уменьшеним радиуса капилляров более весомое значение приобретает конусность. Если при Т\ = 40 мкм для воды влияние конусности проявляется при а > 10-3, то при Т\ = 1 мкм - уже при а > 10-6 Это объясняется увеличением вклада потока J2 в пленках в общий поток J при уменьшении среднего радиуса капилляров [15]. При других одинаковых условиях влияние поля более заметно на перемещение пленок воды, полярность которой (в1 = 78,3) по сравнению с этиленгликолем (в1 = 38,7) выше.

Расчеты показали, что для некоторого значения конусности а = а* суммарный поток массы J = 0, то есть поток пара J1, электродиффузионный поток J3 и электромиграционный поток J4 уравновешиваются обратным пленочным потоком J2. Такие циркуляционные потоки реализуются не в полностью насыщенных жидкостью пористых средах [16]. Поток массы под действием градиента капиллярного давления уравновешивается потоком, вызванным действием силового (в данном случае гравитационного) поля. Локальная конусность капилляра а*, которая соответствует условию J = 0, зависит от градиента напряженности поля и радиуса капилляра. Составленная программа позволяет найти ее значение.

В реальных пористых телах капилляры имеют переменный радиус. Поэтому пузырек воздуха, который перемещается в противоположном к градиенту напряженности направлении, может задержаться в месте расширения поры, если значения r, а и EVE таковы, что выполняется условие J = 0. Это означает, что количество влажного воздуха или пара, а следовательно, и влагосо-держание пористого тела могут не меняться, несмотря на наличие градиента силового поля.

ЗАКЛЮЧЕНИЕ

Проведенные численные расчеты не позволили выявить влияние взаимодействия потоков пара и жидкости в капилляре. Учет влияния неаддитивности потоков, а также капиллярного антидавления [17], вызванного гистерезисом краевого угла (отличием наступающего и отступающего краевых углов менисков), позволит уточнить найденные решения. Тем не менее полученные результаты, в частности, зависимость обратного потока жидкости в пленках под действием градиента капиллярного давления от конусности капилляра, находятся в соответствии с данными работы [18], где исследовался механизм массообмена в коническом капилляре, погруженном в жидкость. Значения суммарного потока J в капиллярах с достаточно малой конусностью, рассчитанные по уравнению (7),

согласуются с экспериментальными данными [9], полученными при перемещении под действием НЭП смачивающих пленок воды в цилиндрических капиллярах.

В дальнейших исследованиях межфазного массообмена необходимо учитывать фактор электризации, объемную и поверхностную проводимость капиллярной жидкости, а также изменение свойств твердой поверхности в электрическом поле [4, 19]. Интересно оценить роль капиллярных волн на поверхности мениска и флуктуаций толщины смачивающей пленки [20, 21], способствующих ее неустойчивости. В связи с этим заслуживает внимания исследование явлений генерации и интенсивного развития автоколебаний столбиков капиллярной жидкости в высокоградиентном НЭП.

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

ЛИТЕРАТУРА

1. Гросу Ф.П., Болога М.К., Болога Ал.М. ЭОМ. 2013, 49(1), 20-29.

2. Панасюк А.Л., Панченко М.С., Старов В.М., Чураев

Н.В. Инженерно-физический журнал. 1978, 35(1), 93-100.

3. Болога М.К., Гросу Ф.П., Кожухарь И.А. Электроконвекция и теплообмен. Кишинев: Штиинца, 1977. 320 с.

4. Жакин А.И. Успехи физических наук. 2013, 183(2), 153-177.

5. Карпович И.Н., Чураев Н.В., Панченко М.С.

Инженерно-физический журнал. 1981, 41(6), 1049-1056.

6. Тамм И.Е. Основы теории электричества. М.: Наука, 1989. 616 с.

7. Аксельруд Г.А., Альтшулер М.А. Введение в капиллярно-химическую технологию. М.: Химия, 1983. 264 с.

8. Дерягин Б.В., Чураев Н.В. Докл. АН СССР. 1972, 207(3), 572-575.

9. Дерягин Б.В., Нерпин С.В., Чураев Н.В. Коллоидный журнал. 1964, 26(3), 301-306.

10. Карпович И.Н., Чураев Н.В., Панченко М.С. Коллоидный журнал. 1984, 46(1), 114-118.

11. Чураев Н.В. Свойства смачивающих пленок жидкостей. В кн. Поверхностные силы в тонких пленках и устойчивость коллоидов. М.: Наука, 1974. 81-89.

12. Киселева О.А., Старов В.М., Чураев Н.В. Коллоидный журнал. 1977, 39(6), 1164-1167.

13. Беллман Р., Калаба Р. Квазилинеаризация и нелинейные краевые задачи. М.: Мир, 1968. 186 с.

14. Викторина М.М., Дерягин Б.В., Ершова И.Г., Чураев Н.В. Докл. АН СССР. 1971, 200(6), 1306-1309.

15. Дерягин Б.В., Чураев Н.В. Смачивающие пленки. М.: Наука, 1984. 160 с.

16. Глобус А.М. ДАН СССР. 1965, 160(1), 81-84.

17. Зорин З.М., Соболев В.Д., Чураев Н.В. Коллоидный журнал. 2001, 63(2), 200-207.

18. Мигун Н.П., Шнип А.И. Инженерно-физический журнал. 2002, 75(6), 145-150.

19. Болога М.К., Гросу Ф.П. ЭОМ. 2012, 48(5), 93-103.

20. Карпович И.Н., Панченко М.С. Инженерно-физический журнал. 2006, 79(5), 20-26.

21.Алиев И.Н., Мильвидский А.Р., Наумов И.А.

Инженерно-физический журнал. 2002, 75(5), 86-87.

Поступила 16.12.14 После доработки 14.09.15 Summary

Calculations have been made of the flows of vapor and a liquid film under the influence of the tension gradients of the inhomogeneous electric field and pressure depending on the radius and taper of the capillary for water and ethylene glycol.

Keywords: tapered capillary, liquid film, non-uniform electric field, mass transfer, porous body.