ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2020, том 30, № 3, с. 9-18 ФИЗИКА ПРИБОРОСТРОЕНИЯ -
УДК 541.182+537
© Б. П. Шарфарец, В. Е. Курочкин, 2020
ЭЛЕКТРОФОРЕЗ НА СУММАРНОМ ПОСТОЯННОМ И ПЕРЕМЕННОМ ЭЛЕКТРИЧЕСКИХ ПОЛЯХ.
I. ПЕРЕМЕННОЕ ЭЛЕКТРИЧЕСКОЕ ПОЛЕ, КРИТИЧЕСКИЙ ОБЗОР СОСТОЯНИЯ ВОПРОСА
Проанализированы предложенные ранее механизмы создания переменным электрическим полем нелинейных стационарных течений вблизи дисперсных частиц: инерционный и поляризационный. Такое рассмотрение предпринято авторами в целях изучения в следующих работах эффекта совместного использования электрофореза на постоянном и переменном электрических полях. В процессе рассмотрения существующих моделей возникновения стационарных течений при воздействии переменного поля вследствие наличия поляризационных эффектов обнаружено, что, кроме стационарного течения и переменного течения на частоте задающего переменного поля, возникает еще и переменное течение на кратной частоте. Эти результаты могут оказаться полезными при возможном осуществлении электрофореза на суммарном постоянном и переменном электрическом поле.
Кл. сл.: электрофорез, электрическое поле, дисперсная частица, переменные течения, стационарные течения, уравнение Навье—Стокса, функции тока, течения с удвоенной частотой
ВВЕДЕНИЕ
Изучению возможности увеличения подвижности ионов и дисперсных частиц при капиллярном электрофорезе в настоящее время уделяется повышенное внимание. Этому способствуют выявленные в последнее время эффекты, связанные с нелинейной зависимостью между подвижностью ионов (дисперсных частиц) и приложенным электрическим полем (см., например, обзоры [1-3], а также [4, гл. 10]). Причем в разных случаях возникает различная прямая степенная зависимость между скоростью частицы и напряженностью приложенного электрического поля. Обычно эти зависимости рассматриваются для случая постоянного электрического поля.
Кроме прямого влияния на подвижность частиц вследствие нелинейной зависимости между скоростью частиц и напряженностью приложенного постоянного электрического поля, следует отметить также сравнительно недавно открытый и исследованный эффект возникновения в окрестностях частиц при приложении переменного электрического поля целого спектра разнообразных течений различной интенсивности [5, 6] (см. также библиографию в них). Все эти течения имеют общую особенность: их скорость квадратично зависит от напряженности приложенного электрического поля [5, 6]. В работе [6] рассматриваются нелинейные течения, вызванные поляризацией двойного
электрического слоя (ДЭС), в то время как в работе [5], кроме этого, рассмотрено еще два механизма, обуславливающих возникновение стационарных нелинейных течений: инерционного и электрокапиллярного.
В работе [5, с. 65] приведены данные о том, что скорости электрофореза в открытом объеме при постоянном электрическом поле более чем в 25 раз превосходят скорости стационарных течений при приложении переменного электрического поля при одинаковых величинах амплитуд напря-женностей электрического поля. Однако в условиях капиллярного электрофореза заряженных частиц вполне может оказаться целесообразным совместное использование и постоянного, и переменного электрических полей. Для этого необходимо критически рассмотреть выбор математического формализма и физических моделей рассматриваемой проблемы.
ПОСТАНОВКА ПРОБЛЕМЫ
Настоящая работа посвящена случаю, когда электрофорез осуществляется на переменном электрическом поле. Целью настоящей работы является критический обзор известных в настоящее время посвященных этому вопросу работ, оценка выбора физических и математических моделей, посвященных электрофорезу на перемен-
ном электрическом поле. Существует несколько механизмов образования стационарных течений жидкости в окрестностях частиц. Здесь будет рассмотрено два из них: инерционный и поляризационный механизмы образования стационарных течений в условиях приложения переменного поля. Рассмотрение сопровождается массированными ссылками на работы [5, 6], в которых приведена обширная библиография по анализируемому вопросу.
ИНЕРЦИОННЫЙ МЕХАНИЗМ ОБРАЗОВАНИЯ СТАЦИОНАРНЫХ ТЕЧЕНИЙ ПРИ ВОЗДЕЙСТВИИ ПЕРЕМЕННОГО ЭЛЕКТРИЧЕСКОГО ПОЛЯ
Воздействие на дисперсные частицы переменного электрического поля вызывает их колебательные движения. В этом случае причиной возникновения стационарных течений в окрестности частиц служат инерционные эффекты. Для этих течений, как упоминалось выше, характерна квадратичная зависимость скорости от амплитуды напряженности электрического поля. Данный эффект проявляется как для твердых частиц, так и для капель, совершающих колебательное движение относительно их сферической формы [5, с. 112].
В работе [5, гл. 4] для исследования заявленной тематики за основу берется неусеченная система гидродинамических уравнений Навье—Стокса для несжимаемой вязкой жидкости (см., напримр, [7, § 16]; [8, с. 39, 44]):
( дУ
р I — + у • Уу | = -vp + )Ду , V-у = 0 .
(1)
можно представить через безразмерную функцию тока для сферических координат [5, с. 113] (см. также Приложение 1)
у =- U 0 ЗФ у = Uo дФ
г 2 дв
г sinв дг
(3)
Здесь г = R / а — безразмерное расстояние от центра частицы, где а — радиус сферической частицы; ио — амплитуда скорости колебательного процесса; Ф — безразмерная функция тока,
равная Ф = у , где у — размерная функция и оа
тока, описанная в Приложении 1; (R ,в,ф) — исходные сферические координаты, отсчитываемые из центра частицы.
После применения к (1) оператора ротор Ух получается следующее общее выражение для вычисления функции тока Ф в сферических координатах [5, с. 113]:
Е4Ф- Re • 8ш в
дФ д дФ д дв дг дг дв
Е 2Ф
г ^ш2 в
-21 а| Ае2ф=о.
. ^ ) дг
(4)
Здесь Re = °и0 — число Рейнольдса процесса;
V
у = г) / р — кинематическая вязкость; 8 =
2v
а
1/2
= I — | — глубина проникновения вязкой вол-
ны; Е2 — оператор Стокса, равный в сферических
д2 sinв д ( 1 д ^
(2) координатах [8, с. 125] Е =—- +
Здесь у — вектор скорости течения жидкости; р и р — соответственно ее плотность и давление в ней; ) — динамическая вязкость жидкости.
Отметим, что при рассмотрении этого механизма не принимается допущение о малости числа Рейнольдса, равного отношению инерционного члена у • Уу к вязкому члену )Ду, и тем самым не отбрасывается нелинейный инерционный член в уравнении (1), что привело бы к линеаризации уравнения движения (1) и соответственно к невозможности моделирования нелинейного течения, которое не может возникнуть в задаче, описываемой линейной моделью. Задача решается в терминах безразмерной функции тока для осесиммет-ричных течений в свободном пространстве. Рассматривается сферическая частица. Тогда (размерные) компоненты вектора скорости у (уR,Ув,t)
дЯ2 Я2 дв\ sinв дв
а применительно к (4)
2 д2 8ш в д ( 1 д
Е2 =^ +
дг2 г2 дв I sin в дв
[5, с. 113];
г = а t — безразмерное время; а — угловая частота переменного электрического поля.
Нелинейное по Ф уравнение (4) решается итерациями с помощью разложения
Ф = Ф0 + Яе Ф1 + (Яе)2 Ф2 +... . (5)
Вначале получено из (4) решение Ф0 (г,в,t) = F0 (г,в)епри Яе = 0. РешенияФп, п = 2, 3, 4,... могут быть получены итеративно после подстановки разложения (5) в (4) и последовательного нахождения решений Ф1, Ф2,... .
Следующее решение ищется в виде (г,в,t) = Фа (г,в) + Фt (г,в,t), а стационарное течение подчиняется уравнению
Е4Ф = sin в х
8Яе (Ф0) а 8 Яе (Ф0) а
8в 8г
8г 8в
Е 2Яе (Ф0
г2 sin2 в
(6)
ствия постоянного E0 и гармонического E1 векторов электрической напряженности), а также уравнения непрерывности потока катионов и анионов (см. Приложение 2) и уравнения Пуассона для потенциала электростатического поля (см. [10, §§ 3.1 и 3.2], а также [5, с. 157]; [6])
Р^ = "¥Р + ^ + Я А + Р^г 1 (7)
где ( ) означает усреднение по времени, а Яе (Ф0) — действительная часть Ф0.
В [5, § 4.2] найдено решение (6) для стационарных течений в окрестности осциллирующих жидких капель, в [5, § 4.3] такое решение найдено для осциллирующих жестких сферических частиц.
В заключение данного пункта отметим, что при нахождении стационарных решений, вызванных инерционными свойствами как жидкости, так и включений, нелинейность учитывалась только в уравнении движения Навье—Стокса (1). Краевые условия, принимавшиеся в задаче, оставались линейными.
ПОЛЯРИЗАЦИОННЫЙ МЕХАНИЗМ ОБРАЗОВАНИЯ СТАЦИОНАРНЫХ ТЕЧЕНИЙ ПРИ ВОЗДЕЙСТВИИ ПЕРЕМЕННОГО ЭЛЕКТРИЧЕСКОГО ПОЛЯ
План изложения
Поляризационный механизм достаточно хорошо изучен (см., например, [4, гл. 6-8]; [5, 6]; [9, гл. 9]; [10, гл. 3, 4, 11]; [12]).
Далее рассмотрим последовательно поляризацию незаряженных или слабо заряженных частиц, когда отсутствует необходимость учета механизмов концентрационной поляризации и поверхностной проводимости частиц [5, 6]. После чего коснемся более сложного механизма поляризации частиц, когда этими механизмами пренебречь нельзя [11, 12]. Остановимся и на случае сильно заряженных частиц, когда проявляется качественно новый эффект, связанный с наличием у них объемного заряда, квадратичного по напряженности электрического поля [13].
Поляризация слабо заряженных частиц
Общая постановка задачи [5, гл. 5]; [6]
Система уравнений, описывающая рассматриваемый процесс, состоит из гидродинамического блока, включающего модернизированное для малых чисел Рейнольдса уравнение движения (1) (называемое в этом случае уравнением Стокса) и уравнения непрерывности для жидкости (2) (будем далее рассматривать случай совместного дей-
8Т
V-у = 0,
V-г+^ = 0, 81
Дф = -
4Р
(в СГСЭ).
(8) (9) (10)
Здесь рассматривается случай симметричного электролита с модулем зарядного числа ионов Z и концентрацией зарядов положительного с+ и отрицательного с- ; ре — плотность электрического заряда; г, г' — соответственно диэлектрическая проницаемость жидкости и частицы. В качестве вектора ^ рассматривается плотность потока этих ионов (см. Приложение 2, (П6))
•± ± г~\±гт ± , ±
\ = с Ушр - О ^ с + с у = = ±с± л± Е - О± V- с + с± у,
(11)
где Л — подвижность ионов; О — коэффициент диффузии; Е — вектор напряженности приложенного электрического поля.
Далее будем рассматривать стационарную и нестационарную задачи отдельно.
Нестационарная задача. Диэлектрическая частица в проводящей жидкости [5, § 5 .1.1]; [6]
Принимается, что в декартовых координатах поле Е1 ориентировано вдоль оси Oz: Е! (0,0, Е1),
где Е1 = Е 1в-им. Электрическое поле Е1 изменяет концентрации положительных с+е—м и отрицательных с~е~иМ ионов. Результирующая плотность заряда ре1 в (10) определяется из выражения
Ре1 (R,в, t) = ге (с+ (R,в) - с^ (R,в))е~ш =
= 2е5с (R,в) е-1Ш. (11)
1) Отметим, что в оригинальных работах [5, 6] автор рассматривает однородное уравнение движения (7). В качестве источника фигурирует только переменное поле Е1. Действие этого источника проявляется в [5, 6]
в краевом условии и в распределении поля в отсутствие поляризации.
х
г
Здесь е — заряд протона; 8с = с + - с — разностная концентрация ионов.
Пусть 8с± (Я,в)е= с±(Я,в)е-а - с0 — невязки концентраций катионов и анионов относительно их равновесной концентрации с0 .
Допущение 1
Принимается допущение о малости числа Пекле Ре = Уа / D С 1, где У — характерная величина
амплитуды скорости; а — радиус частицы; D — коэффициент диффузии ионов (принимается одинаковым D = D+ = D ). Как показано в [10], это условие справедливо даже для заряженных частиц в случае слабой поляризации, т.е. при условии
Sc± (R,6) < со.
(12)
-iaSc = DA Sc + 2с0цАф. Окончательно из (10) и (13) получено
A Sc = y2Sc,
(13)
(14)
Sc (R,6 t ) = -у-p, (R,e, t),
где
p0 (R,6,t) = -E1Rcos6e~ia
f Ъ ^ \ ^ cos6 ( (R,6, t) = a3E1-- x R
(16)
(17)
1 -
3k
1 +
2sy
2 Л
v sK /
(1 - ya)2 +1 -(ya)
2y2
1 +
2s
(1 + ya )2 +1 -(ка)
(18)
(p (R,6, t ) =
3aE1e
-y(R-a) /
yR (Ka )2
1 + — yR
\
1 +
2s
-cos6e
(19)
(1 + ya )2 +1 -(Ka)
При Ре С 1 в (11) можно пренебречь конвективным членом с± у . После этого система (7)—(10) распадается на две независимые подсистемы: гидродинамическую (7), (8) и электрическую (9), (10). Задача поляризации таким образом решается с помощью только электрической системы (9), (10).
Для упрощения выкладок принимается также допущение о равенстве коэффициентов диффузии D = D+ = D и подвижности ¡л = /и+ = ¡л- ионов. Тогда (9) сводится к одному выражению
Здесь (р0 — потенциал, вызванный исходным полем Е1, определяемый из условия, что напряженность, соответствующая потенциалу, равна
Е(Я,в,г) = Е(Я,в)е,
E( R,6)e"
= E
где
E,
где у2 = к2 - ia / D ; к = Я4 ; Я — длина Дебая. Подставка решения (14) в выражение для ре1 = 2е8с(Я,в), а затем в уравнение Пуассона (10) и учет граничных условий на границах дает окончательно для диэлектрических частиц
р( Я,в, г ) =
= Р0 (Я,в, г) + ра (Я,в, г) + рр (Я,в, г), (15)
определен выше: E1 (x, y, z, t ) = = ( 0,0, E) e-M, E = const.
Составляющая pd (19) суммарного потенциала, вызванная поляризацией частицы, представляет собой осциллирующий дипольный момент. Наконец еще одно следствие поляризации, вызванной внешним электрическим полем, — появление скачка потенциала pp.
Как видно из (15)-(19), пока все решения исходной системы (9), (10) линейны, являются гармоническими функциями частоты со исходного электрического поля (отсутствуют постоянная составляющая и кратные гармоники).
Допущение 2. Тонкий двойной слой
Выполняются условия тонкого двойного слоя Ka ^ 1, частота поля много меньше частоты максвелл-вагнеровской поляризации ю ^ Dk2 . В этом случае выражения (16), (18) и (19) упрощаются [5, 6]:
Pd (R,6, t ) = -
a3 E1 cos6
--
2 R2
/ D а Л 3aE1 1 -K(R-a)
Pp (R,6, t)=--e
2 s
1 +— Ka
(20)
cos6e-iat, (21)
Sc (R,6, t ) = -c0
2eZ
~kT
Pp (R,6,t) .
(22)
2
x
e
2
s
X
X
2
Б
s
Допущение 3. Малость числа Духина
(
ехР
Du
2кТ
ка
^ 1. Здесь — штерновский
потенциал. Приведенные выше выкладки справедливы при отсутствии у частицы поверхностного заряда, наличие которого, согласно [10], вызывает необходимость учета эффектов, связанных с поверхностной проводимостью и концентрационной поляризацией. Для частиц с тонким двойным слоем относительный вклад этих эффектов характеризуется величиной числа Би. Отсюда следует, что при условии тонкого двойного слоя поляриза-
е^, ,
цию слабо заряженных частиц, т.е. когда -— < 1,
2кТ
можно не только качественно, но и количественно описывать как поляризацию незаряженных частиц [10]. В связи с этим предложенная в [5, 6] теория в случае частиц с тонким двойным слоем пригодна для описания не только незаряженных, но и также слабо заряженных частиц, отвечающих условию Би < 1.
При этом найденное распределение рр (R,в,t) и разностной концентрации ионов 8с(Я,в,/)
в пределах двойного слоя необходимо изменить на величину сферически симметричных доба-вокЯ) и 8с0 (Я), характеризующих распределение соответственно потенциала и концентрации объемного заряда в равновесном двойном слое [5, с. 166]; [6].
Кроме поляризации диэлектрических частиц, в [5, § 5.1.2] рассмотрен также случай поляризации идеально поляризуемой частицы. Техника решения этой задачи идентична технике, приведенной выше. Отличия состоят лишь в граничных условиях.
Отметим в заключение этого раздела, что до сих пор рассматривались только линейные модели, при которых не должно создаваться эффекта возникновения стационарных течений.
Стационарные поляризационные течения, вызванные переменным электрическим полем.
Диэлектрическая частица в проводящей жидкости [5, § 5.2]; [6]
Электроосмотическое скольжение жидкости по поверхности частицы вызвано действием тангенциальной составляющей электрического поля на распределенные в подвижной части двойного слоя электрические заряды. При тонком двойном слое ка ^ 1 полагается, что скорость жидкости на
поверхности частицы совпадает со скоростью электроосмотического обтекания этой частицы.
Допущение 4 [5, с. 170]; [6] Изменение скорости электроосмоса по поверхности частицы за счет выделения в диффузном слое двойного слоя при ка ^ 1 поляризационных
зарядов формально можно учесть с помощью добавки в формуле Смолуховского к равновесному ^-потенциалу величины ф°р = рр (а,в, t), действительной составляющей выражения (21) для скачка потенциала (21) в точке Я = а: рр (Я,в, t)| званного поляризацией двойного слоя
вы-
0 3аЕ1 Р =
1
2 £ 1 +--ка
-cosвcosat.
(23)
Кроме того, меняется и значение электрического поля, вызывающего движение: кроме тангенциальной составляющей поля Е1, добавляется тангенциальная составляющая электрического поля, вызванного скачком потенциала поля ра. В результате тангенциальная составляющая скорости жидкости при Я = а равна [5, с. 171]; [6]
V
в 1я=
4 яг]
(£ + Р)
д(Рр + Ра )
адв
(24)
Здесь рр определяется из (23):
р0 (Я = а,в,t) = -aElCOsвcosat. (25)
ра определяется в (20) при Я = а и с заменой экспоненты ее действительной частью
Ра (a,в,t) = -
а3 Е1 cosв 2
аЕ1
-cosюt =
2
cosвcosюt.
(26)
Радиальная компонента скорости на поверхности частицы равна нулю
= 0.
Я\Я=а
(27)
На бесконечности для радиальной скорости ставится краевое условие
Ук\ = и0 = Ее~1а'.
(28)
Для описания осесимметричных в данном случае течений в окрестности частицы решается система уравнений (7), (8), (24), (27), (28). Для этого в [5, 6] используется функция тока Ф (3). Согласно [5, 6], решение (как отмечает автор, линейной
Б
Е
а
Я=а
задачи) ищется в виде
Ф( г,в, г) = Ф0 (г,в, г) + Ф р (г,в, г), (29)
где Ф0 (г,в,г) соответствует решению задачи об
электрофорезе в переменном электрическом поле без учета поляризации двойного слоя и соответствует члену Ф0 в разложении (5), который раскрыт в [5, с. 119] и представляет собой периодическую функцию времени с частотой а; Ф (г,в,г) — искомая добавка к решению, обусловленная только поляризацией частицы. Функция Ф (г,в,г) ищется в [5, 6] в виде
Ф р (г,в, г ) = Фа (г,в) + Фг (г,в, г), (30)
означает ос-
где Фа (г,в) = (фр (г,в,г)} г, где ( )[
реднение по времени.
В [5, с. 172, 173]; [6] приведены окончательные уравнения и граничные условия для определения функции тока Фа, а также связанной с ней скорости течения:
E 4Ф„ = 0,
Ф.
дФ„
Ф
= 0,
дг
= -a sin2 OcosO.
(31)
(32)
(33)
Решение системы (31)-(33) имеет вид [5, с. 173]; [6]:
Фа (r,0) = ---1jsin2 0cos0 , и соответственно для компонент скорости:
V (r,0) = -^f0(1 -1 j(3sin2 0 -1) ,
1 .
Ve (r,0) = aU0—sinOcosO r
где
a = -
9 sE¡a 1
32U0 Щ 1 + s
H— ка s'
Аналогичные результаты получены в [5, 6] и для идеально поляризуемых сферических частиц. В [5, 6] приведены различия в решениях для этих двух случаев, на которых здесь не останавливаемся.
Замечания по поводу приведенных в [5, 6]
результатов, полученных для нелинейных поляризационных течений
Остановимся на некоторых замечаниях к изложенным выше результатам, касающимся поляризационных течений.
Во-первых, изначально автором рассматривается линейная (как он сам утверждает в [6, с. 363]) задача (7), (8), (24), (27) и (28) (нумерация соответствует настоящей работе). После чего делается вывод о наличии возможности искать ее решение в виде суммы периодической функции времени с частотой а , совпадающей с частотой внешнего электрического поля, и не зависящей от времени функции (см. выражения (29), (30)). Однако известно, что в линейной задаче при наличии периодического по времени источника невозможно получение апериодического решения.
Во-вторых, автор рассматривает линейную однородную систему уравнений движения Навье— Стокса (7), (8). Однако краевое условие (24) по умолчанию предполагает нелинейность задачи в силу того, что, несмотря на наличие исходного поля частотой а , явочным порядком появляются постоянная составляющая а = 0 и кратная частота 2а. Поэтому для справедливости утверждения о линейности рассматриваемой задачи необходимо сторонний источник характеризовать краевыми условиями (24)-(28).
В-третьих, если рассмотреть краевое условие (24), то выясняется неотмеченное автором [5, 6] еще одно качество получаемого решения. А именно: разложим краевое условие (24) на слагаемые
V
4лг]а
fAd^oL +
дО дО
+
дУо + ^ дО дО
Таким образом, на примере диэлектрической сферической частицы показан механизм образования поляризационных нелинейных течений, индуцированных переменным электрическим полем.
Учитывая выражения для pp (23), p0 (25) и pd (26), а также что й = const, делаем вывод о том,
„ ( 8p0 dp Л
что первое слагаемое в нем й I —0 + —- I содер-
v 86 86 )
жит гармоники только исходной частоты со, в то
0 (8p0 dp, Л
время как слагаемое pp I + I содержит как
постоянную составляющую, соответствующую со = 0, отвечающую за наличие стационарного те-
r
r =1
е
R=a
X
R = а
чения, так и гармонику двойной частоты 2а. Таким образом, в процессе присутствует три, а не две составляющие поляризационного течения: стационарное течение, соответствующее а = 0, периодическое по времени течение с частотой а и периодическое по времени течение с двойной частотой 2а.
По мнению авторов, способ решения задачи (7)—(10) в работах [5, 6] является несколько искусственным, поскольку изначально ставится линейная задача, а затем привносятся нелинейные артефакты типа краевого условия (24). Видимо, поэтому в [5, 6] не была обнаружена составляющая решения на кратной частоте 2а , что, впрочем, не умаляет изобретательности автора работ [5, 6] в решении столь сложной задачи.
ВЫВОДЫ
В работе проанализированы два предложенных ранее механизма создания переменным электрическим полем нелинейных стационарных течений вблизи дисперсных частиц. Из существующих рассмотрены два механизма, инициирующих создание стационарных течений: инерционный и поляризационный. Такое рассмотрение предпринято авторами в целях изучения в следующих работах эффекта совместного использования электрофореза на постоянном и переменном электрических полях. В процессе рассмотрения существующих моделей возникновения стационарных течений при воздействии переменного поля вследствие наличия поляризационных эффектов обнаружено, что, кроме стационарного течения и переменного течения на частоте задающего переменного поля, возникает еще и переменное течение на кратной частоте. Эти результаты могут оказаться полезными при возможном осуществлении электрофореза на суммарном постоянном и переменном электрическом поле.
ПРИЛОЖЕНИЕ 1. ФУНКЦИИ ТОКА
Для течений несжимаемой жидкости часто используется метод функции тока, позволяющий свести решение к нахождению единственной скалярной функции — функции тока х, у) . Для двумерных задач используется функция тока Ла-гранжа [8, § 3.1]. Существуют классы трехмерных течений, которые также можно единственным образом описать с помощью одной скалярной функции. Практически очень важным случаем является обтекание тела вращения потоком жидкости, параллельным его оси симметрии (осесимметричные течения). В этом случае используется метод функ-
ций тока, предложенный Дж. Стоксом. Функция тока в общем случае несжимаемой жидкости, когда справедливо полное уравнение движения На-вье—Стокса [8, с. 124]
(ду ^
р I — + V • Vv I = -Ур + ]Ду,
подчиняется динамическому уравнению
г4
уЬ у - у
(ду д ду д^Ь2у д
дz ду ду дz
7
дt
(Ь2у) = о
В случае уравнения Стокса
ду _ р — = -Ур + Т]Ду , дt
когда силы вязкости превосходят инерционные силы (малые числа Рейнольдса), функция тока подчиняется уравнению [8, с. 124], но в уравнении для функции тока нелинейный член выпадает:
д
уЬ> Ь2У) =
дt
В стационарном случае это уравнение максимально упрощается
Ь4у = 0.
Дифференциальный оператор Ь2 (оператор Стокса) в цилиндрических координатах выражается так [8, с. 123]:
Р2 д Ь = у —
д
7-
ду ^ ду
> д2
+-;
дz'
А 2 д
= Д---.
у ду
Фигурирующие в динамических уравнениях независимые переменные у = у (у, 7) таковы: г — ось симметрии; у — расстояние, отсчитываемое по нормали к оси симметрии, что совпадает с азимутально-симметричной цилиндрической системой координат (р, 7) .
В сферических координатах выражение для Ь2 имеет вид [8, с. 125]:
Ь2 =
д2 sin 0 д
дг2 +
1 д
г2 д0 I sin0 д0
Выражение для оператора Ь в криволинейных координатах вращения (q1,q2,ф) представлено
в работе [8, выражение (4.7.17)].
Граничные условия для функции тока приводятся в работе [8, § 4.13]. Кроме того, в [8, § 4.23] приведена схема получения решения уравнения Ь4 у = 0 с помощью решения уравнения Ь2 у = 0.
Отметим, что после получения решения у = = у (у, 7) компоненты скорости течения могут быть определены из соотношения [8, с. 119]
v = —
y
I [«p.
^ ip y
Что дает для цилиндрической системы координат
1 8^
vp = p, p 8z
1 8w
vz =---
P 8P
и для сферической системы координат
1 8цт 1 8у
'а — "
V = —
' D
R2 sin 6 86
R sin 6 8R
8p w • 8t
(П1)
Здесь р — объемная плотность некоторой величины q; ] — плотность потока величины q; и — скорость изменения величины q в единице объема. В случае, если величина q сохраняется (например, энергия), то (П1) записывается в виде
8P + V- j = 0.
8t
(П2)
8Pe + V- j = о
8t J e
(П3)
где ие — добавление величины заряда q в единице объема в единицу времени за счет внешнего источника.
В случае <ге = 0 уравнение (П3) приводится к виду (П2)
p+ V-j = 0, 8t J
(П4)
В более общем случае, если (д2,р) образуют правую систему ортогональных криволинейных координат вращения, составляющие скорости (у1,у2 ) определяются выражением (4.4.3) работы [8].
ПРИЛОЖЕНИЕ 2. УРАВНЕНИЕ НЕПРЕРЫВНОСТИ
Дифференциальная форма общего уравнения непрерывности такова [14]:
Электромагнетизм. В случае, если в качестве величины q рассматривается совокупный электрический заряд некоторого объема жидкости, тогда аналогом объемной плотности р следует рассматривать объемную плотность заряда ре, в качестве плотности потока заряда следует рассматривать вектор плотности электрического тока и (П1) перепишется в виде
а в стационарном случае, когда плотность заряда не зависит от времени pe (x, t) = pe (x), выражение (П4) сводится к уравнению
V-j = 0. (П5)
Далее для других физических явлений рассматриваем случаи сохранения величины q, т.е. о=0.
Теория волн. В теории волн уравнение непрерывности выражает собой закон сохранения энергии в элементарном объеме, в котором распространяются волны любой природы:
- + V- j = 0.
8t
Здесь w — плотность энергии W; j = = j (x, y, z, t) — вектор плотности потока энергии в точке четырех векторов (x, y, z, t).
Гидродинамика. В гидродинамике уравнение непрерывности записывают так (его также называют уравнением неразрывности)
8P + V-( pv ) = 8P + pV- v + v -Vp = 0,
8t 8t
где p = p( x, y, z, t) — плотность жидкости; v = v ( x, y, z, t) — вектор скорости жидкости; j = pv — вектор плотности потока жидкости с направлением, совпадающим с направлением вектора скорости. Для однородной несжимаемой жидкости p = const, и уравнение непрерывности упрощается
V-v = 0.
Движение ионов в жидкости. В работе [10, гл. 3, §§ 1, 2] уравнение непрерывности рассматривается применительно к следующим величинам: в качестве величины q принимается количество в некотором объеме некоторых одинаковых ионов. Тогда в качестве плотности величины q рассматривается объемная концентрация этих ионов c (количество ионов в единице объема). В качестве вектора j рассматривается плотность потока этих ионов (количество ионов, проходящих через единичную площадку, нормальную направлению их результирующей скорости в единицу времени). Вектор j описывается уравнением Нернста—
Планка (см., например, [2]; [15, с. 157]; [16, с. 245])
•+ + гл+гт + , +
Г = с-утщг - С>У-с + с~ у =
= ±с± ¡и± Е - D±V • с± + с± у. (П6)
Здесь с — концентрация рассматриваемого вида ионов; D — коэффициент диффузии этого иона; у — вектор скорости окружающей жидкости; у 1швг — миграционная скорость иона, равная [15, с. 144]
у тбг = иЕ = Е ; 6жт]а
и — миграционная подвижность иона; Z — зарядное число иона (валентность с учетом знака заряда); е — заряд протона; Т] — динамическая вязкость жидкости; а — радиус иона; верхний индекс ± относится соответственно к величинам, характеризующим катионы и анионы.
Окончательно уравнение непрерывности потока ионов записывается стандартно в виде (см., например, [10, с. 76])
—+ V- j = 0. dt
(П7)
где вектор плотности потока ионов ] определяется уравнением Нернста—Планка (П6).
В случае неизменной концентрации ионов последнее уравнение модифицируется к виду (см., например, [10, с. 61])
V-Г = 0.
Работа выполнена в ИАП РАН в рамках Государственного задания 075-00780-20-00 по теме № 0074-20190013 Министерства науки и высшего образования.
3. Мищук Н.А., Баринова Н.О. Теоретическое и экспериментальное исследование нелинейного электрофореза // Коллоидный журнал. 2011. Т. 73, №: 1. С. 74-82. URL: https://eHbrary.ru/item.asp?id=15598772
4. Духин С.С., Дерягин Б.В. Электрофорез. М.: Наука,
1976. 328 с.
5. Мурцовкин В.А. Электрогидродинамические эффекты в поверхностных явлениях. Дис. ... д. ф.-м. н. М.: Инст. физ. химии, 1992. 300 с.
6. Мурцовкин В.А. Нелинейные течения вблизи поляризованных дисперсных частиц // Коллоидный журнал. 1996. Т. 58, № 3. С. 358-367.
7. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидродинамика. М.: Наука, 1986. 736 с.
8. Хаппель Дж., Бреннер Г. Гидродинамика при малых числах Рейнольдса. М.: Мир, 1976. 630 с.
9. Левич В.Г. Физико-химическая гидродинамика. М.Ижевск: Институт компьютерных исследований, 2016. 708 с.
10. Духин С.С., Шилов В.Н. Диэлектрические явления и двойной слой в дисперсных системах и полиэлектролитах. Киев: Наукова думка, 1972. 207 с.
11. Шилов В.Н., Духин С.С. Теория поляризации диффузной части двойного слоя сферической частицы в переменном электрическом поле // Коллоидный журнал. 1970. Т. 32, № 1. С. 117-123.
12. Симонов И.Н., Шилов В.Н. Теория поляризации диффузной части тонкого двойного слоя сферических проводящих частиц в переменном электрическом поле // Коллоидный журнал. 1973. Т. 35, № 2. С. 381-385.
13. Симонова Т.С., Духин С.С. Нелинейная поляризация диффузной части тонкого двойного слоя сферической частицы // Коллоидный журнал 1976. Т.38, №1. С. 7985.
14. Continuity equation // Wikipedia.
URL: http s ://en.wikipedia. org/ wiki/Continuity_equation
15. Bruus H. Theoretical microfluidics. Oxford University Press, 2008. 346 p.
16. Ньюмен Дж. Электрохимические системы. М.: Мир,
1977. 464 с.
СПИСОК ЛИТЕРАТУРЫ
Жуков А.Н., Курочкин В.Е., Шарфарец Б.П. О нелинейности электрокинетических явлений. Обзор // Научное приборостроение. 2020. Т. 30, № 1. С. 17-21. URL: http://iairas.ru/mag/2020/abst1.php#abst2 Курочкин В.Е., Шарфарец Б.П. О подвижности ионов при электрофорезе в однородной жидкости и в пористых средах в случае линейной зависимости скорости от амплитуды вектора электрической напряженности // Научное приборостроение. 2020. Т. 30, № 2. С. 33-40. URL: http://iairas.ru/mag/2020/abst2.php#abst5
Институт аналитического приборостроения, Санкт-Петербург
Контакты: Шарфарец Борис Пинксович, sharb@mail.ru
Материал поступил в редакцию 8.05.2020
ISSN 0868-5886
NAUCHNOE PRIBOROSTROENIE, 2020, Vol. 30, No. 3, pp. 9-18
ELECTROPHORESIS ON THE TOTAL CONSTANT AND HARMONIC ELECTRIC FIELD.
I. HARMONIC ELECTRIC FIELD, A REVIEW OF THE STATE OF THE ISSUE
B. P. Sharfarets, V. E. Kurochkin
Institute for Analytical Instrumentation of RAS, Saint Petersburg, Russia
The previously proposed mechanisms for creating an alternating electric field of nonlinear stationary flows near dispersed particles: inertial and polarizing, are analyzed. Such a consideration was undertaken by the authors in order to study in the following works the effect of the joint use of electrophoresis in a constant and an alternating electric field. In the process of considering the existing models of the appearance of stationary flows under the influence of an alternating field due to the presence of polarization effects, it was found that in addition to the stationary flow and the sinusoidal flow at the frequency of the driving sinusoidal field, an sinusoidal flow occurs at a dual frequency. These results may be useful in the possible implementation of electrophoresis in the total constant and alternating electric field.
variable flows, stationary flows,
Keywords: electrophoresis, electric field, dispersed particle, the Navier-Stokes equation, double frequency flows
REFERENCES
1. Zhukov A.N., Kurochkin V.E., Sharfarets B.P. [On the non-linearity of electrokinetic phenomena. Overview]. Nauch-noe Priborostroenie [Scientific Instrumentation], 2020, vol. 30, no. 1, pp. 17-21. DOI: 10.18358/np-30-1-i1721 (In Russ.).
2. Kurochkin V.E., Sharfarets B.P. [On the mobility of ions during electrophoresis in a homogeneous liquid and in porous media in the case of a linear dependence of the velocity on the amplitude of the electric voltage vector]. Nauchnoe Priborostroenie [Scientific Instrumentation], 2020, vol. 30, no. 2, pp. 33-40. DOI: 10.18358/np-30-2-i3349 (In Russ.).
3. Mishchuk N.A., Barinova N.O. [Theoretical and experimental study of nonlinear electrophoresis]. Kolloidnyj zhurnal [Colloid Journal], 2011, vol. 73, no. 1, pp. 74-82. DOI: 10.1134/S1061933X11010133
4. Duhin S.S., Deryagin B.V. Elektroforez [Electrophoresis]. Moscow, Nauka Publ., 1976. 328 p. (In Russ.).
5. Murtcovkin V.A. Elektrogidrodinamicheskie effekty v po-verhnostnyh yavleniyah. Diss. dok. fiz. mat. nauk [Elec-trohydrodynamic effects in the superficial phenomena. dock. phys. math. sci. diss.]. Moscow, Frumkin Institute of Physical chemistry and Electrochemistry Russian academy of sciences, 1992. 300 p. (In Russ.).
6. Murtcovkin V.A. [Nonlinear flows near polarized particulate matter]. Kolloidny zhurnal [Colloidal journal], 1996, vol. 58, no. 3, pp. 358-367. (In Russ.).
7. Landau L.D., Lifshiz E.M. Teoreticheskaya fizika. T. 6. Gidrodinamika [Theoretical physics. Vol. 6. Hydrodynamics]. Moscow, Nauka Publ., 1986. 736 p. (In Russ.).
Contacts: Sharfarets Boris Pinkusovich, sharb@mail.ru
8. Happel J., Brenner H. Low Reynolds number hydrodynamics. Springer, Netherlands, 1983. 553 p. (Russ. Ed.: Happel Dzh., Brenner G. Gidrodinamika pri malyh chis-lahRejnoldsa. Moscow, Mir Publ., 1976. 630 p.
DOI: 10.1007/978-94-009-8352-6.
9. Levich V.G. Fiziko-himicheskaya gidrodinamika [Physical and chemical hydrodynamics]. Moscow—Izhevsk, Institute of Computer Science, 2016. 708 p. (In Russ.).
10. Dukhin S.S., Shilov V.N. Dielektricheskie yavleniya i dvojnoj sloj v dispersnyh sistemah i polielektrolitah [Dielectric phenomena and double layer in disperse systems and polyelectrolytes]. Kiev, Naukova dumka Publ., 1972. 207 p. (In Russ.).
11. Shilov V.N., Dukhin S.S. [Theory of polarization of the diffuse part of the double layer of a spherical particle in a variable electric field]. Kolloidny zhurnal [Colloidal journal], 1970, vol. 32, no. 1, pp. 117-123. (In Russ.).
12. Simonov I.N., Shilov V.N. [Theory of polarization of the diffuse part of the thin double layer of spherical conducting particles in an alternating electric field]. Kolloidny zhurnal [Colloidal journal], 1973, vol. 35, no. 2, pp. 381385. (In Russ.).
13. Simonova T.S., Dukhin S.S. [Nonlinear polarization of diffuse part of thin double layer of spherical particle]. Kolloidny zhurnal [Colloidal journal], 1976, vol. 38, no. 1. pp. 7985. (In Russ.).
14. Continuity equation // Wikipedia.
URL: https://en.wikipedia.org/wiki/Continuity_equation
15. Bruus H. Theoretical microfluidics. Oxford University Press, 2008. 346 p.
16. Newman J. Electrochemical Systems. Wiley-Interscience, 672 p. (Russ. Ed.: Nyumen Dzh. Elektrohimicheskie sis-temy. Moscow, Mir Publ., 1977. 464 p.).
Article received by the editorial office on 8.05.2020