Электроосмос в тонкой щели Текст научной статьи по специальности «Физика»

Физика

УДК 532.546

Электроосмос в тонкой щели В. В. Кадет, П. С. Корюзлов

Кафедра подземной гидромеханики РГУ нефти и газа им. И. М. Губкина Ленинский пр-т, д. 65, Москва, Россия, 119991

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

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

1. Введение

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

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

Впоследствии, на основании этих данных были предприняты попытки теоретического моделирования электроосмотического течения в пористой среде [5].

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

Статья поступила в редакцию 4 апреля 2008 г.

2. Постановка задачи

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

к\\\\ччч\ч\ччччччто

Течение

\ч\ччччччч\\\\\ччта

3

Н/2

X

1чч\\\учччч\\ччччччч\\\\\чччч\\чччч\та

Рис. 1. Расположение координатных осей и направление вектора напряжённости

электрического поля

На границе твёрдой и жидкой фаз возникает двойной электрический слой (ДЭС) [6], вследствие чего стенки канала оказываются равномерно заряженными, с дзета-потенциалом £. В непосредственной близости к межфазной поверхности находится неподвижный, так называемый адсорбционный слой ДЭС, или слой Гельмгольца, толщина которого не более диаметра гидратированных ионов [6]. За ним располагается диффузная часть ДЭС (слой Гуи), которая после наложения на систему внешнего электрического поля начинает двигаться в результате взаимодействия поля с заключённым в ней объёмным зарядом. В качестве толщины этой части ДЭС обычно принимают I = 1/к, где к — параметр Дебая-Хюккеля [7].

Объёмная плотность заряда ре(у) есть при этом функция поперечной координаты у и определяется уравнением Пуассона [8]

¿2Ф(У) Ре{у)

&У2 £о£г

где £ о — диэлектрическая проницаемость вакуума, ф(у) — электрический потенциал ДЭС.

Рассмотрим уравнение Навье-Стокса для такого течения:

дУ -> ->

р!й + р(у )У = Р + (1)

Вследствие малого раскрытия канала силой гравитации можно пренебречь, и массовая сила у будет представлять собой только действие приложенного электрического поля напряжённостью Е к заряду с объёмной плотностью ре(у). Согласно [9], вследствие малости скорости у инерционный член (У • VУ)У оказывается пренебрежимо малым. Здесь р — плотность электролита, Рх = Ере(у) — плотность массовых сил в проекции на ось Х. При этих предположениях уравнение (1) в проекции на ось X можно представить в виде

ди д 2и (2)

рж = »ду2 -Е£г£0^~. (2)

Начальными и граничными условиями являются: состояние покоя до момента включения внешнего электрического поля, симметрия относительно центра потока и условие прилипания на стенке.

3. Распределение потенциала в сечении канала

Как видно из (2), поле скоростей определяется распределением потенциала ДЭС, поэтому, прежде всего, проанализируем поведение потенциала ф(у). Запишем уравнение Пуассона-Больцмана [6] в безразмерном виде для щелевого микроканала:

ф =8ЬФ(у). (3)

На рис. 2 приведён график распределения потенциала ДЭС в канале: зависимость безразмерного потенциала Ф = (ге^)/(кьТ) (вертикальная ось) от безразмерной координаты у = ку (горизонтальная ось)

¥

2,0 Н

1,5 -1,00,5 -

0 _ 0 2 4 6 8 10 У

Рис. 2. Распределение потенциала ДЭС

Расчёт сделан при следующих значениях параметров: объёмная концентрация ионов п = 6,23 ■ 1022 1/м3 (10-4 моль/л), соответствующая дзета-потенциалу ( = 50 мВ раствора NaCl, валентность ионов z = 1, материал стенок канала — стекло, плотность жидкости р = 998 кг/м3, вязкость ц = 0,9 ■ 10-3 Па-с, диэлектрическая проницаемость электролита ег = 80, раскрытие канала Н = 3 ■ 10-7 м.

Зависимость (-потенциала от концентрации ионов предполагалась линейной на интервале от 50 до 150 мВ. Этому диапазону изменений (-потенциала соответствуют концентрации 10-4-10-6 моль/л.

Из приведённого на рис. 2 графика видно, что в пристеночной области, где объёмная концентрация ионов резко возрастает, увеличивается и потенциал ДЭС — его максимум достигается на стенке канала, в неподвижной адсорбционной части ДЭС.

На рис. 3 приведены результаты расчётов для значений дзета-потенциала 50 мВ, 100 мВ и 150 мВ, которые показывают, что с увеличением дзета-потенциала увеличивается толщина диффузного слоя ДЭС и, соответственно, возрастает электрический потенциал, наведённый его ионами. Для £ =150 мВ потенциал остаётся отличным от нуля даже в окрестностях центра канала при его раскрытии 3 ■ 10-7 м.

4. Поле электроосмотической скорости

Полученный выше результат позволяет перейти непосредственно к решению уравнения (2).

Для того, чтобы привести (2) к виду линейного уравнения теплопроводности, введём функцию U(y,t) = u(y,t) - usip(y)/(, где u(y,t) = Vx, us = Ee0erQ/p — скорость Смолуховского [6].

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

ч? 8 -

6 -

4 -

2 -

О 2 4 6 8 10 7

Рис. 3. Профиль распределения потенциала ДЭС при разных дзета-потенциалах:

1 — С = 50 мВ, 2 — С = 100 мВ, 3 — С = 150 мВ

*.0=-о - ^)+£» (^-кьь (- (^-Ц*

(кН л

[' - cos (^ ) *} . (4)

где и = (u/us), Ф = (zеф)/(кьТ), ( = (zeQ/(къТ), у = ку, t = (v/H2)t (здесь n, и z — объёмная концентрация ионов и валентность ионов соответственно, е — заряд протона, къ — постоянная Больцмана, Т — абсолютная температура, к =

1 /2

((2z2e2n,)/(къТ)) — параметр Дебая-Хюккеля, а 1/к — толщина диффузной части ДЭС).

Рассмотрим влияние основных параметров на распределение скорости. График зависимости (4), приведённый на рис. 4, показывает, что картина формирования электроосмотического течения отличается от таковой для пуазейлев-ского течения — движение жидкости начинается с тонкого слоя вблизи стенки канала, то есть в диффузной части ДЭС, в которой концентрация ионов, способных к передвижению под действием внешнего поля, максимальна.

Рис. 4. Профили скорости электроосмотического течения в зависимости от времени

За время порядка нескольких миллисекунд скорость течения в области ДЭС достигает максимума. При этом вдали от стенки она падает практически до нуля, так как в этой области концентрация ионов мала. Течение в слое ДЭС шириной порядка I (см. рис. 1) можно назвать активным вследствие того, что именно оно приводит в движение остальные слои жидкости благодаря вязкостным силам. Профиль скоростей развившегося электроосмотического течения близок к таковому для турбулентного течения — жидкость движется с одинаковой скоростью почти по всему сечению канала, за исключением небольшой области вблизи стенки канала.

Поперечный размер канала влияет, главным образом, на время установления течения и величину средней скорости. Для оценки времени установления учтём, что выход на стационарный режим происходит, когда экспоненциальный множитель в (4) стремится к нулю. Соответственно, с течением времени это слагаемое становится пренебрежительно мало и его можно отбросить. Если принять, что существенное влияние на зависимость от времени этот множитель перестаёт оказывать при значении порядка 0,001, то есть меньше десятой процента от максимальной скорости течения, то оценкой времени установления можно считать

га — з(я>). (5)

С уменьшением толщины щели максимум профиля скорости в канале смещается в сторону стенки канала. Средняя скорость течения, при прочих равных условиях ((-потенциале равном 50 мВ, и напряжённости 100 В/м), больше для канала с большим Н.

Для случая неперекрывающихся ДЭС в плоском микроканале максимальная скорость при развившемся течении наблюдается в области, где электрокинетический потенциал практически равен нулю Ф(у) — 0.

Можно также отметить, что максимальная скорость, достигаемая при развившемся течении, не зависит от распределения ионной концентрации и размеров канала. Этот результат, полученный на основании (5), согласуется с экспериментальным фактом [9].

5. Ограничения, накладываемые на модель в случае переменного циклического внешнего поля

Строгое математическое моделирование транспорта ионов в ДЭС, генерируемого приложенным извне переменным электрическим полем, должно учитывать нестационарные эффекты электропереноса. Однако, согласно работе [10], порядок времени отклика ДЭС на влияние поля составляет 10-8-10-7 с. Это, как минимум, на порядок меньше, чем характерное время эволюции электроосмотического течения (5), которое имеет значение 10-6-10-4 с для раскрытий 10-6-10-5 м. Поэтому при временах больше указанных процессом «релаксации» ДЭС можно пренебречь. Например, даже для частоты поля 500 КГц при раскрытии 10-6 м характерное время составляет 2,3 ■ 10-6 с. Кроме того, в [10] отмечается, что при приложении переменного поля равновесное распределение Больцмана для ионной концентрации справедливо вплоть до частоты 1 МГц. Следовательно, уравнение Пуассона-Больцмана остаётся правомерным и в этих условиях.

6. Двухфазное электроосмотическое течение

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

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

ёж/

+ ^ Г^)

-1

X

н н

Ь - х/

ь -Г, 41Р1Е*1 ^ ^ -

Рг2Н

2

н н

0 0

2(Ь - х/)

-4

J J Рг Рхг ¿у ¿у

00

(6)

где Ь — длина щели, индекс % = I, г относится к свойствам слева и справа от фронта, Рс — капиллярное давление на границе раздела двух жидкостей, произведение Рг^хг(у) определяется из (2) и (3).

Для проверки результатов теоретических расчётов был спланирован и проведён цикл экспериментов. Использовавшаяся экспериментальная установка, детально описанная в работе [11], позволила проследить влияние (-потенциала, раскрытия щели и напряжённости внешнего поля на скорость движения фронта раздела двух жидкостей при электроосмотическом течении.

В качестве вытесняющей жидкости использовались растворы МаС1 различной концентрации (0,01 Н, 0,005 Н, 0,001 Н). Вытесняемой жидкостью было минеральное масло с известными реологическими характеристиками.

После определения дзета-потенциалов жидкостей на данной подложке, для случаев различной напряжённости поля, различных дзета-потенциалов вытесняющих растворов и раскрытия канала были проведены измерения зависимости положения X/ фронта раздела двух фаз от времени. В качестве координаты фронта бралась координата, полученная усреднением минимального и максимального значения X/ на данный момент времени. Фактические данные сравнивались с расчётными, полученными по зависимости (6).

На рис. 5 приведены графики для различных значений напряжённости поля (53 В/см, 75 В/см, 83 В/см). Раскрытие канала составляет 50 мкм, дзета-потенциал вытесняющего раствора равен 97 мВ.

0,7 0,9 1,1 1,3 1,5 1,7 1,9 2,1 2,3 2,5 Время, ч

Рис. 5. Модельные зависимости положения фронта от времени при различных напряжённостях поля: 1 — 83 В/см, 2 — 75 В/см, 3 — 53 В/см. Ромбами, треугольниками и кругами изображены соответствующие измеренные значения

X

Расчётные средние скорости движения: 7,46 ■ 10-6 м/с для напряжённости 53 В/см, 1,05 ■ 10-5 м/с для напряжённости 75 В/см и 1,17 ■ 10-5 м/с для напряжённости 85 В/см.

Сравнение данных эксперимента и модели для различных значений дзета-потенциала вытесняющего раствора приведено на рис. 6. Напряжённость поля постоянная: 75 В/см.

Как и в случае с увеличением напряжённости поля, наибольшая негладкость фронта наблюдалось при максимальных скоростях течения, в данном случае при 0 = 110 мВ.

S ff

о о.

■е<

0,14 0,12 0,1 0,08

g, 0,06

0,04 0,02 0

I/

--И

-Г

0,7

0,9

1,1 1,3 Время, ч

1,5

1,7

Рис. 6. Модельные зависимости координаты фронта от времени при различных

дзета-потенциалах: 1 — 110 мВ, 2 — 93 мВ, 3 — 87 мВ. Квадратами, треугольниками и ромбами изображены соответствующие измеренные значения

Модельный расчёт для канала с раскрытием 10 мкм, напряжённостью поля 75 В/см и дзета-потенциалом 97 мВ показал значение скорости движения фронта, почти в полтора раза меньшее по сравнению с наблюдаемой. Это можно связать с тем, что при таких раскрытиях заметно возрастает влияние поверхностной проводимости канала. Так как в модели это явление не учитывается, то возможным вариантом корректировки могло бы быть определение эффективного дзета-потенциала для канала с таким раскрытием. Измеренная средняя скорость фронта для раскрытия 100 мкм (1,8-10—5 м/с) отличается от модельной всего на 9%.

7. Выводы

Решена задача о течении раствора электролита во внешнем электрическом поле в тонкой щели с диапазоном раскрытия 10-6 м < Н < 10-8 м. Показано, что время установления течения прямо пропорционально квадрату раскрытия щели и обратно пропорционально кинематической вязкости ts ~ Н2 ¡v. Экспериментально подтверждено, что диапазон раскрытия щели, при котором внешнее поле оказывает существенное влияние на движение раствора электролита, лежит в пределах нескольких толщин ДЭС. Получены аналитические выражения для профиля потенциала ДЭС и скорости течения в случае постоянной напряжённости внешнего электрического поля.

На процесс эволюции течения сильно влияют ширина раскрытия канала, электрокинетические свойства раствора и материал стенок канала. Максимальная скорость при развившемся течении наблюдается в области, где электрокинетический потенциал практически равен нулю и не зависит от распределения ионной концентрации и размеров канала. Для каналов с раскрытием большим 10-7 м установившийся профиль течения является фактически однородным (и = const), с уменьшением раскрытия структура течения усложняется. Для использования модели в случае переменного поля накладывается ограничение на его частоту « 500 кГц).

Наибольшей эффективности вытеснения можно добиться путём снижения концентрации вытесняющего электролита при постоянной напряжённости внешнего поля.

Было показано, что при раскрытии меньшим 10 мкм необходимо введение поправки на поверхностную проводимость канала. Эта поправка может быть реализована дополнительным измерением дзета-потенциала для данного раскрытия щели.

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

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

Литература

1. Electroosmotic Capillary Flow with Nonuniform Zeta Potential / M. G. Garguilo, J. I. Molho, J. G. Santiago et al // Anal. Chem. — 2000. — Vol. 72. — Pp. 10531057.

2. Cummings E. B., Griffiths S. K., Nilson R. H. Irrotationality of Uniform Elec-troosmosis // Proc. SPIE Microuidic Devices and Systems. — Vol. 3877, No 2. — 1999. — Pp. 180-189.

3. Fluid Transport Mechanisms in Micro Uidic Devices / J. M. Molho, A. E. Herr, M. Desphande et al // Proc. ASME Micro-Electro-Mechanical-Systems (MEMS). — No 66, 66. — 1998. — Pp. 69-76.

4. Bianchi F., Ferrigno R. Finite Element Simulation of an Electroosmotic-Driven Flow Division at a T-Junction of Microscale Dimensions // Anal. Chem. — 2000. — No 72.

5. Kang Y., Yang C., Huang X. // J. Micromech. Microeng. — 2004. — No 14. — Pp. 1249-1257.

6. Фролов Ю. Г. Курс коллоидной химии (Поверхностные явления и дисперсные системы): учебник для вузов. — М.: Химия, 1982. — 280 с.

7. Ландау Л. Д., Лифшиц Е. М. Статистическая физика. — М.: Наука, 1976. — 470 с.

8. Будтов В. П. Физическая химия растворов полимеров. — СПб.: Химия, 1992.

9. Schlichting H. Boundary-Layer Theory. — New York: McGraw Hill, 1979. — 310 p.

10. Hsu J. P., Kuo Y. C.and Tseng S. J. Dynamic Interactions of Two Electrical Double Layers // J. Colloid Interf. Sci. — 1997. — Vol. 198.

11. Батырбаев Н. Д., Кадет В. В. Современный этап разработки нефтяных месторождений западного Казахстана. Проблемы и решения. — М.: ОАО ВНИ-ИОНГ, 2006.

UDC 532.546

Electroosmosis Flow in Slit MicroChannel

V.V. Kadet, P. S. Koryuzlov

Oil and Gas Underground Hydromechanics Department Gubkin Russian State University of Oil and Gas 65, Leninskiy prospekt, Moscow, 119991, Russia

The process of electroosmotic flow in slit microchannel which thickness is of an order of several EDL thicknesses is investigated in this paper. The model and derived analytical solutions allow to predict velocity profiles and time of developing the steady state flow in cases of both alternating and constant electric field. The results of the simulation are in good agreement with the obtained experimental data.