Базовые выражения, используемые в электрокинетике (обзор) Текст научной статьи по специальности «Физика»

ISSN 0868-5886

НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2014, том 24, № 4, c. 13-21

МАТЕМАТИЧЕСКИЕ МЕТОДЫ И МОДЕЛИРОВАНИЕ В ПРИБОРОСТРОЕНИИ

УДК 541.13

© Н. Н. Князьков, Б. П. Шарфарец, Е. Б. Шарфарец

БАЗОВЫЕ ВЫРАЖЕНИЯ, ИСПОЛЬЗУЕМЫЕ В ЭЛЕКТРОКИНЕТИКЕ (ОБЗОР)

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

Кл. сл.: электрокинетические явления, уравнение Навье—Стокса, уравнение Нернста—Планка, уравнение теплопереноса, уравнение Пуассона

ВВЕДЕНИЕ

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

ЭЛЕКТРОСТАТИКА

Приведем некоторые базовые выражения, касающиеся электростатики и закона сохранения зарядов. Основной задачей электростатики является отыскание векторов напряженности электрического поля E и электрической индукции (электрического смещения) D в каждой точке электрического поля, созданного данной системой источников поля — электрических зарядов [1, с. 248].

Уравнения электростатики являются частным случаем полной системы уравнений Максвелла для электромагнитного поля [1, с. 501], когда его

характеристики — электрическая индукция (электрическое смещение) D и напряженность электрического поля E — не зависят от времени. Содержательные уравнения в системе уравнений Максвелла имеют в этом случае вид [1]:

V- D ) = ре1 (в СИ),

V- D = V • (^ ) = 4пре1 (в СГСЭ);

V х E = 0 (в СИ и СГСЭ). (2)

Исходя из (2), поле E является потенциальным E = ^р. (3)

Здесь £0 = 8.85-10-12 Кл2/(Н• м2) = 8.85-10-12 Ф/м — электрическая постоянная; £ — относительная диэлектрическая проницаемость среды, показывающая во сколько раз сила взаимодействия между двумя постоянными зарядами в данной среде меньше чем в вакууме; р — скалярный потенциал; ре1 = —--объемная плотность заряда; ^ —

дУ

величина заряда, сосредоточенного в объеме дУ.

Для вектора D справедливо представление

D = ££^ (в СИ); D = £E (в СГСЭ).

Нормальная составляющая вектора D к границе раздела двух диэлектриков (при отсутствии заряда на границе) непрерывна (не претерпевает разрывов) в отличие от нормальной составляющей вектора E . Это следует из (1) при ре1 = 0 на границе.

Выражение (1) является формулировкой теоремы Остроградского—Гаусса в дифференциальной форме.

При заданном распределении свободных зарядов потенциал р электростатического поля в однородном и изотропном диэлектрике в £ раз меньше чем потенциал р0 в вакууме и удовлетворяет уравнению Пуассона, что следует из (1) и (3)

Лр = - р (в СИ),

££

Лр = -^^ (в СГСЭ). (4)

Приведем также дифференциальную форму закона сохранения заряда — уравнение непрерывности для объемной плотности зарядов [1, с. 504]

V-] +

дРе,

а

= 0.

(5)

Здесь ] — плотность токов проводимости.

Отметим, что в иностранной литературе, например в [2], применяются несколько отличные обозначения: через £ обозначается не относительная диэлектрическая проницаемость среды, как в отечественной литературе, а совокупная диэлектрическая проницаемость среды, равная произведению электрической постоянной £0 на относительную диэлектрическую проницаемость среды, обозначаемую в этом случае через £г.

УРАВНЕНИЕ НАВЬЕ—СТОКСА ПРИ ВОЗДЕЙСТВИИ ВНЕШНЕГО ЭЛЕКТРИЧЕСКОГО ПОЛЯ

Закон сохранения импульса при движении вязкой несжимаемой жидкости, к которой приложено внешнее электрическое поле с напряженностью Е, в дифференциальной форме записывается в виде уравнения Навье—Стокса (см., например, работу [2, с. 141])

Р"а V-V)V = ~УР + ^Лу + РШ + Ре,Е . (6)

Здесь V, р — поля вектора скорости и давления в жидкости; р — ее плотность, п — вязкость; g — вектор ускорения силы тяжести; ре1 — по-прежнему объемная плотность электрического заряда.

ЭЛЕКТРОХИМИЧЕСКИЙ ПОТЕНЦИАЛ

Движущей силой диффузии служит пространственная вариация химического потенциала данного вещества. В отсутствие электрического поля химический потенциал / растворенного вещества для условий, при которых вкладом гидростатического давления можно пренебречь, равен [5, с. 79]

/и = /и0 + RT 1п—, С0

где /и0 — стандартный химический потенциал, зависящий от природы растворителя при постоянной концентрации примеси частиц с0 и отсутствии электрического поля; с — реальное поле концентрации частиц (понятие концентрации ниже уточняется в силу ее неоднозначного определения). Для ионов, движение которых зависит не только от концентрации, но и от электрического потенциала, соответствующим понятием является электрохимический потенциал / [5, с. 79]:

с

/ = /и0 + RT 1п--+ ZFр.

с0

Электрохимический потенциал — мера работы, необходимой для переноса 1 моля ионов из раствора с данной концентрацией и данным электрическим потенциалом в бесконечно удаленную точку в вакууме. Эта работа складывается из затрат на преодоление сил химического взаимодействия /и0 + RT 1п с и работы по переносу зарядов в электрическом поле ZFр [5, с. 80].

Иногда работу по переносу ионов относят не к молю, а к одному иону [5, с. 80]

с

/ = /и0 + kBT 1п--+ Zeр .

С0

Здесь учтено, что kB = R / Ыл и е = F / Ыл . Введены обозначения: Ыл « 6.022141-1023 моль1 — число Авогадро; F — число Фарадея, F = еЫл« « 96485.33Кл - моль-1; Z — заряд иона в единицах заряда протона (безразмерная величина, равная валентности иона с учетом знака его избыточного за-Дж

ряда); R « 8.31446 - — универсальная газомоль - К

вая постоянная; kB = 1.3806488 -1023 Дж/Кл — постоянная Больцмана; T — абсолютная температура в К; е = 1.602176565 -10-19 Кл — элементарный электрический заряд (заряд протона).

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

В работе [2, с. 146] дано определение электрохимического потенциала отдельно для концентраций положительно и отрицательно заряженных ионов (о чем будет говорить соответствующий

Б

спаренный знак ±)

(r) = + kBT ln

c± (r)

± Zep(r).

(7)

Знак " +" соответствует положительным, а знак " -" — отрицательным ионам. Ионы в примеси здесь полагаются обратно заряженными Z+ = Z_ = Z .

В условиях термодинамического равновесия электрохимический потенциал постоянен: Д± (r) = const [2, с. 146]. Исходя из этого, после воздействия на обе части равенства (7) оператором V получается равенство

f

kß^

T ln

c± (r)

Л

= + ZeVp(r),

-о У

которое при малых градиентах температуры VT сводится к более простому равенству

1п С^) = ^ р(г),

С0 квТ

откуда получается равенство (статистическое уравнение Больцмана), связывающее концентрацию ионов с электрическим потенциалом

c± (r) = Со exp

V kT

\

(p(r)

= c0 exp

+ZF RT

P(r) I. (8)

ЭЛЕКТРИЧЕСКИИ ПОТЕНЦИАЛ

ясно, что потенциал будет зависеть только от нормальной к этой плоскости координаты (например г) и в этом случае он равен [2, с. 147]

Р( г ) = ^ZT Arth

th

4kßT у

exp

г

Ad

Y

(11)

При получении (11) использовались краевые условия c± (да) = c0, (р(да) = 0, <p(surf) = д , где <p(surf) — потенциал на поверхности стенки; д — дзета-потенциал (см., например, [4, 5, 7]). В (11)

ad =

SS0kBT

р ( Ze)

e )2 c

ss0RT 2Z2 F2 C о

(12)

т. н. дебаевская длина, или радиус экранирования

~ С

[5, т. 2, с. 88]. Здесь С0 = —— — концентрация

NА

в моль/м3 (см. ниже). Как видно, длина экранирования существенно зависит от концентрации ионов: в разбавленных водных растворах длина экранирования (т. е. область диффузного двойного слоя) достигает десятков нанометров, а в концентрированных растворах — десятых долей нанометров [там же, с. 89].

В случае выполнения приближения Дебая (или модели Гуи), когда тепловая энергия превосходит электрическую энергию [2, с. 147]

Zeg ■ kBT,

(13)

Электрический потенциал получают из уравнения (4). Сошлемся на работу [2]. Для вычисления величины объемной плотности заряда справа в (4) используется выражение [2, с. 147]

Pel = Ze (С+ (r) - С- (r)) .

(9)

Здесь с+ (г) (с-(г)) — концентрация, подразумевающая число положительно (отрицательно) заряженных ионов в единице объема в точке г . Концентрация с+ (г) находится из выражения (8).

После подстановки (8) и (9) в уравнение Пуассона (4) получено следующее нелинейное уравнение для электрического потенциала (уравнение Пуассона—Больцмана) (см. [2, с. 147]):

Zec,

Ap(r) = 2-0sh -p(r)

SSn

Ze

V kßT

(10)

получают линеаризованное приближение уравнения (10) в виде уравнения Дебая—Хюккеля

1

AP(rЬттP(r) .

An

(14)

Например, в случае зависимости только от координаты г при г > 0 и граничном условии р = д

(точное решение уравнения (10) имеет вид (11)) решение (14) таково [2, с. 148]:

p( z ) = gexp

ad

(15)

В случае цилиндрического капилляра радиусом а уравнение (14) преобразуется к виду

2

Уравнение (10) в общем случае необходимо решать численно, однако в ряде случаев существует аналитическое решение. Один из них — это случай однородно заряженной плоскости. Тогда

1 5

2 +--

5r r 5r

P(r) = 7ГР(г) , r e[0,a], (16) ad

p(a ) = g.

(17)

c

0

z

С 2

Член —у в уравнении (16) равен нулю вслед-дz

ствие предположения о цилиндрической симметрии задачи. Решением этой краевой задачи является [2, с. 149]

р( г ) = Я

Ь(Г/Лр ) ^ (а / Л ) '

(18)

Здесь 10 — модифицированная функция Бесселя нулевого порядка.

Отметим, что по известному распределению потенциала можно легко найти распределение плотности электрического заряда из уравнения (4),

а именно

Ре! = -££0ЛР .

(19)

Для потенциалов (15) и (18) справедливы соответственно выражения [2, с. 148, 149]

Ре! = -££0 ЛР = -££

д 2р "Ж2"

££0Я -^ехр

Лр

(20)

Ре, =■

££

0? 10 (г / Л )

Л 10 (а / Л )

с± (z) = с0ехр

+—— Р( z)

1 + ехр (-z / Л )

kвT

(22)

с± (г) = Са ехр

+-- р(г)

kвT '

1+¿е?10 (г/ Л)

kBT 10 (а / Л )

(23)

др

[8, с. 81], V = Угег , и Vp = —, и, в силу осевой

дг

симметрии, = Уг (г) . Тогда стационарное уравнение (6) без учета силы тяжести записывается в виде

-др+пЛу2 + РЕ = 0. дг

Далее вывод приводится подробно вследствие допущенных в [7] неточностей. После выделения из последнего уравнения плотности заряда

др дг

1

^д2

- + -

1 д

Л

дг г дг

Ре! =■

Е

и подстановки его в уравнение (4) получается уравнение

д2 1 д -2 +--

дг г дг

Л

Р =

1

££0 Ег

др гь'

1

д2 1 д

""Т +--

дг г дг

Уг

У У

решением которого является следующее поле ско-

(21) ростей:

Использование приближения Дебая (13), (14) позволяет в ряде случаев легко найти из (8) выражения для концентраций с± (г). Для плоскости это выражение выглядит так [2, с. 148] (здесь внешняя экспонента заменяется первыми двумя членами ряда):

■--

др а2 - г2

(г) = е (р-Р )-Л ' г 1 ^ *1=а} дг 41

(24)

При интегрировании выражения (24) были приняты следующие краевые условия: у2 = 0 при г = а, Уг и р при г = 0 конечны.

В работе [2, с. 163] выражение (24) представлено в законченном виде путем подстановки в него выражения для потенциала (18), полученного

,10 (г / Лр ) : 10 (а / ) '

из уравнения Дебая—Хюккеля, р( г ) = ?

Для цилиндрического капилляра использование (18) и замена экспоненты первыми двумя членами разложения дает следующее выражение:

Уг ( г ) = Е2 ^ ? 1

I0(г/Лр) ,

—/-7 -1

10 (а / )

др а2 - г2 дг 41

(25)

Из выражения (24) видно, что при отсутствии внешнего электрического поля (Ег = 0), поле скоростей совпадает с течением Пуазейля в круглой трубке

X г ) = -

др а2 - г2

у (г)

дг 41

В работах [2, с. 162] и [7, с. 220] решается уравнение баланса импульса, следующее из уравнения Навье—Стокса для стационарного течения в цилиндрической трубке с осью вдоль оси Ог при предположении Е = Егег. Тогда, как известно

др

а при равенстве градиента давления нулю — = 0

дг

течение определяется как осмотическое, если учесть, что Р\г=а = ? Уг (г)= Ег 1 (Р- ?) .

У

г

г

с

0

с

0

Замечание. По оценкам [2, с. 159] на практике при наличии электроосмоса имеет место следующий порядок величин:

- д -потенциал обычно д = 100 мВ;

- скорость течения в процессе электроосмоса \го = 1 мм/с;

- подвижность для различных ионов представлена в достаточно обширной справочной литературе.

УРАВНЕНИЕ ДИФФУЗИИ ДЛЯ РАСТВОРОВ С УЧЕТОМ КОНВЕКЦИИ

Пусть в растворе присутствует n примесей, включая растворитель. Тогда масса а-й примеси в единице объема будет равна ее парциальной плотности ра. Таким образом, плотность раствора, очевидно, равна

р=Ё р«.

(26)

а г dlff а

J = J dlff + J conv.

а а а

(27)

Плотность потока массы за счет конвекции равна

JConv = Сру . (28)

Диффузионная плотность потока массы при умеренных концентрациях примеси определяется уравнением Фика [2, с. 94]

Jdlff = -D pVC .

а а г а

(29)

выраженное в терминах безразмерного коэффициента концентрации С а (см., например, [2, с. 94]):

dt

+ v-VC = D АС .

а а а

(30)

Уравнение неразрывности для раствора в целом остается таким же, как и для однородной (гомогенной) среды [2, 6 и др.]

^ + V•( ру ) = 0.

ел у '

Безразмерная концентрация Са определяется так [2, с. 94]

С = р р

Плотность потока массы примеси , очевидно, равна сумме конвекционной плотности потока массы JConv и диффузионной плотности потока массы Jd

Для слабых растворов (т. е. при умеренных концентрациях примесей, когда и справедливо уравнение Фика (29)) эквивалентная плотность смеси р приблизительно равна плотности растворителя р*. Уравнение сохранения (30) может быть

Na

переписано в терминах концентраций с а = ,

где N а — число молекул примеси а в объеме V . Из короткой цепочки равенств видно

С = N0.=р^^т. = рр_ С

а -г г *тг а '

V та PV та

т. е. концентрации с а и С а — линейно зависимые функции. Здесь та — масса одной молекулы примеси а . Отсюда уравнение (30) сразу переписывается в терминах концентраций са [2, с. 94]

dEa dt

+ v-Ve = D„Ас .

а а а

(31)

Отметим, что в неподвижной жидкости, когда у = 0, конвективный перенос примесей отсутствует и уравнения (30) и (31) переходят к классическим уравнениям диффузии.

УРАВНЕНИЕ ЭЛЕКТРОДИФФУЗИИ

Источником диффузии является наличие ненулевого градиента химического потенциала [5, с. 79]. Для ионов, движение которых зависит не только от концентрации, но и от электрического потенциала, дополнительно к диффузии возникает миграция ионов под воздействием электрического поля. В этом случае источником движения частиц является наличие ненулевого градиента электрохимического потенциала [5, с. 79-80]. При наличии диффузии, конвекции и миграции плотность потока массы (22) должна учесть и составляющую плотности потока миграции заряженных частиц а-компонента раствора:

J

Здесь Ва, [Ва ] = м2/с — парциальный коэффициент диффузии примеси а .

Используя (27)-(29), легко получить уравнение неразрывности (или непрерывности, или сохранения массы) для а-й примеси в дифференциальной форме,

J dlff + J CDI

а а а

+j :gr.

В работе [7, с. 245] представлена суммарная плотность потока массы J а (уравнение Нернста—Планка)

J =-uaFCа Vp- DaVCa + Сa V..

(32)

Здесь 5« , 5« ] = —2—; V — вектор скорости сме--1 м с

~ моль

си в целом; С«, —— , — соответствующая кон-м

центрация а-компоненты; Ра — коэффициент диф-

, , . г п моль - с м2 - моль

фузии (см выше); иа , [иа]=-=-, —

кг Дж - с

подвижность ионов; га — заряд иона вида а в единицах заряда протона (безразмерная величина, равная валентности иона с учетом знака его избыточного заряда).

Результирующий ток в растворе электролита равен [7, с. 246]

I = F^Zа~Jа . (33)

а

Условие сохранения массы имеет вид [7, с. 246]

дСа

дл

= 5« + Я, =

= ^-(-г^Сс^р- DаVCа + CаV) + Яа . (34)

Здесь Я« — источник, обусловленный гомогенной химической реакцией.

Наконец практически во всем объеме, занятом раствором, кроме тонкого двойного электрического слоя, принимается условие его электрической нейтральности [7, с. 247]

I г«Са = 0.

(35)

Приведем физический смысл уравнений (32)-(35) [7, с. 247]. Эти четыре уравнения дают замкнутое описание процессов переноса в растворах электролитов: согласно уравнению (32) компоненты в растворе могут двигаться за счет миграции (первое слагаемое справа), диффузии (второй член) и конвекции (третий член). Уравнение (33) просто отражает тот факт, что электрический ток обусловлен движением заряженных компонентов. Уравнение (34) дает условие материального баланса, а уравнение (35) — условие электронейтральности.

Замечание 1. Отметим, что уравнение потока (32) не выполняется для концентрированных растворов [7, с. 248].

Замечание 2. В области правомерности уравнения потока (32) оказывается, что подвижности и коэффициенты диффузии ионов связаны. Эта связь осуществляется уравнением Нернста— Эйнштейна [7, с. 258]

которое строго применимо при бесконечном разбавлении.

Покажем связь соотношений (32) и (36) с распределением Больцмана (8) [7, с. 215].

Поток каждого из растворенных компонентов описывается уравнением Нернста—Планка (32). В случае равновесного состояния (V = 0, 5« = 0) имеем

-ZаUаFCаVр- О^С а = 0.

После подстановки в последнее уравнение значения Ра из (36) получается уравнение

УСа

Са

RT

Vр,

откуда и следует распределение Больцмана (8).

Замечание 3. Покажем, что концентрация, оп-

N

ределенная выше са = , где Nа — число молекул примеси а в объеме V, и концентрация С а, фигурирующая в выражениях, начиная с уравнения (32), равны между собой с точностью до постоянного множителя. Сначала отметим, что наличие Nа частиц в объеме V равносильно наличию в этом объеме Nа / Nл молей. По определению,

Са =

число молей N с

аа

V

НлУ Мл

Поэтому уравнение сохранения массы справедливо для обеих концентраций при условии, что

источник Яа

моль

м3 - с

в (34) заменяется на

Л т ~ число частиц га = NлЯа , -3-, а именно

м3 - с

О = Я1и

(36)

дс

= -V - (-zauaFсaVр- Рсуса + с«V) + г« . (34а)

ЭЛЕКТРОЛИТИЧЕСКАЯ ДИССОЦИАЦИЯ

Электролитической диссоциацией называют распад молекул растворенного вещества на положительно и отрицательно заряженные ионы в результате взаимодействия с растворителем. Это явление обусловлено совместным действием двух причин [1, с. 386, 387] (далее в этом пункте все выражения взяты из этого источника): теплового движения и взаимодействия полярных молекул растворенного вещества, состоящих из взаимосвязанных ионов, с полярными молекулами раствори-

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

а =

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

В условиях динамического равновесия между процессами диссоциации и молизации удовлетворяет уравнению

1 - а

—— = const • n0.

а

const

F

j = —Z+ no+(u+ + u_) E =

N A

(37)

где Z+ — валентность положительных ионов в растворе; п0+ — число положительных ионов

в единице объема электролита; и+ и и — подвижности соответственно положительного и отрицательного ионов, т. е. средние скорости движения этих ионов под воздействием электрического поля, напряженность которого равна единице.

Основные выражения таковы:

- закон Ома

] = иЕ,

где и — удельная проводимость среды;

- закон Гаусса

V • (££0Е) = ре1;

- уравнение непрерывности

■ + V• \ = 0,

SPel

St

(38)

(39)

(40)

Если п0 ^ 0, то а ^ 1, т. е. в слабых растворах а «1 и почти все молекулы диссоциированы. С увеличением концентрации раствора убывает. В сильно концентрированных растворах

где ре1 — удельная плотность электрического заряда.

Подставляя (38) и (39) в (40), легко получаем уравнение

SPi

а

+ — Pi = 0,

а

Отметим, что вода также диссоциирует. Константа диссоциации воды при температуре +25 °С равна 1.86 10-16. Диссоциация воды — процесс эндотермический (т. е. поглощающий тепло), поэтому с повышением температуры, в соответствии с принципом Ле-Шателье, диссоциация усиливается.

ЗАКОН ОМА ДЛЯ ТОКА В ЖИДКОСТЯХ

Плотность тока в жидкостях равна сумме плотностей токов положительных и отрицательных ионов

] = + ]+ .

Зависимость плотности тока ] в жидкостях от напряженности Е электрического поля, приложенного к электродам, имеет вид [1, с. 388]

St ss0 решение которого имеет вид

Pel = Po exP

(41)

££

где х = —0 — время релаксации заряда. и

Для хороших проводников х ~ 0(10~19) с, для хороших изоляторов х ~ 0(103) с.

В зависимости от соотношения величин х и ( различают область, в которой необходимо изучать процесс.

При х >> Л (изоляторы) имеем приближение электростатики, при х << Л (проводники) и х ~ Л имеем приближение электрического тока: стационарного (х << Л), либо переменного (х ~ Л).

ЭЛЕКТРОФОРЕЗ

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

n

0

n

0

t

т

в тепловом движении. Слой ионов на частице и неподвижные и диффузные противоионы и образуют имеющий нейтральный заряд двойной электрический слой. После приложения статического электрического поля напряженностью E частицы начинают двигаться в нем в направлении, определяемом знаком заряда поверхности частицы. Пусть заряд частицы Ze, где Z целое число (отрицательное или положительное), e — элементарный заряд (электрона). Тогда электрическая сила Fel, действующая на частицу, равна

Fei = ZeE .

При движении на частицу действует сила сопротивления Стокса

Fdrag = -bn1Uep .

Здесь uep — вектор скорости электрофореза, которая определяется из баланса действующих на частицу сил

Fdag + Fei =-6n71üUeV + ZeE = 0.

Отсюда

Ze bnia

■E.

(42)

r =

kBT bnDi

стоянии покоя в неподвижной жидкости. В [2] приведена версия уравнения этого процесса

4 з du

— na р — = -bniau + F 3 Ир dt ' ex

(43)

Здесь рассматривается одномерное движение частицы плотностью рр и радиусом а со скоростью и в поле вязких сил Стокса (первый член справа) и постоянной силы Fext, приложенной в начальный момент времени t = 0 . Решение (43) имеет вид

К„

u = -

bnia

f f 1 - exp

91

2ра

w

//

(44)

В работе [2] показано, что, согласно (44), время, характеризующее стабилизацию процесса т =

2Ppa2 9i

-, при р , примерно совпадающей с плот-

ностью воды, равно т = 5 мкс. Тогда из (44) следует, что через время t ~ 2т = 10 мкс частица практически достигает стационарной скорости

F„

Замечание. Фигурирующий в (42) радиус а иона не принимается равным его реальной величине а « 0.05 нм, а вследствие процессов гидратации принимается равным а « 0.2 нм. Это является следствием того факта, что ионы присоединяют к себе слой молекул воды (растворителя) [2, с. 145]. В этой же работе приведена таблица значений ионной подвижности различных ионов в водных растворах.

В работе [5, т. 2, с. 80] приведено выражение для радиуса иона с учетом гидратации

Здесь kB — постоянная Больцмана (kB «

«1.38 •Ю-23 Дж• К-1); Т — абсолютная температура (К); D, м2 •с-1, — коэффициент диффузии; г), Па• с, — динамическая вязкость растворителя.

СКОРОСТЬ ПЕРЕХОДНОГО ПРОЦЕССА ПРИ РАБОТЕ С МИКРОЧАСТИЦАМИ

Коснемся теперь скорости переходного процесса при работе с микрочастицами. В работе [2] рассматривается ускоренное движение сферической частицы, находящейся в начальный момент в со-

и ~ и = .

6жт]а

Отметим, что при получении этой оценки были приняты некоторые допущения. Например, выражение для силы сопротивления движению частицы ^ = 6лщ (сила Стокса) справедливо только в стационарном режиме движения частицы с постоянной скоростью. В случае, если частица обтекается жидкостью, уравнение для скорости частицы определяется системой уравнений движения, описывающей возмущенное наличием включения несущее течение жидкости, включающей также и уравнение для скорости самого включения, которая, как правило, не совпадает с линиями тока и амплитудой скорости несущего течения [3, с. 28-30].

СВЯЗАННЫЕ СИСТЕМЫ УРАВНЕНИЙ ДЛЯ ПОЛНОГО РЕШЕНИЯ ЗАДАЧИ

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

- уравнение Навье—Стокса (6);

-уравнение материального баланса — (30)

(случай без электродиффузии) и (34) (при ее наличии);

- уравнение Пуассона (4);

- в случае наличия переменного теплового поля

t

u

сюда еще необходимо добавить уравнение тепло-переноса, которое при малых возмущениях может быть записано в линеаризованном виде (см., например, [9])

1 еТ' =ДТ '+а 0Т0 еР'

Ь St

к0 St

Здесь абсолютная температура процесса равна Т = Т'+ Т0, где Т0 — равновесная температура, а Т' — ее возмущение, при этом полагается Т' ■ Т0; х0 — коэффициент температуропроводности жидкости; к0 — коэффициент теплопроводности жидкости; р' — избыточное давление;

1 (еУч

VI ет

а0 ^ —| | — термодинамический коэффици-

ент расширяемости при постоянном давлении.

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

ВЫВОДЫ

СПИСОК ЛИТЕРАТУРЫ

1. Яворский Б.М., Детлаф А.А. Справочник по физике. М.: Наука, 1971. 939 с.

2. Bruus H. Theoretical microfluidics. Oxford University Press, 2008. 346 p.

3. Дружинин О.А. Исследование динамики вихревых потоков и волн в дисперсных и стратифицированных средах. Дис. ... д-ра физ.-мат. наук. Нижний Новгород: ИПФ РАН, 2004. 300 с.

4. Духин С.С., Дерягин Б.В. Электрофорез. М.: Наука,

1976. 332 с.

5. Рубин А.Б. Биофизика. Т. 2. М.: Наука, 2004. 469 с.

6. Нигматулин Р.И. Основы механики гетерогенных сред. М.: Наука, 1978. 336 с.

7. Ньюмен Дж. Электрохимические системы. М.: Мир,

1977. 464 с.

8. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидродинамика. М.: Наука, 1988. 736 с.

9. Шарфарец Б.П., Князьков Н.Н., Пашовкин Т.Н. О математической постановке задачи движения вязких сжимаемых теплопроводящих жидкостей в термоупругой трубке // Научное приборостроение. 2013. Т. 23, №. 4. С. 85-90.

Институт аналитического приборостроения РАН, г. Санкт-Петербург

В работе дан обзор основных выражений, позволяющих математически моделировать сложные задачи электрохимии вообще и электрокинетические задачи в частности в условиях связанности различных взаимодействующих в них полей.

Контакты: Шарфарец Борис Пинкусович, sharb@mail.ru

Материал поступил в редакцию: 24.09.2014

THE BASIC EXPRESSIONS USED IN THE ELECTROKINETIC PHENOMENA

(REVIEW)

N. N. Knyaz'kov, B. P. Sharfarets, E. B. Sharfarets

Institute for Analytical Instrumentation of RAS, Saint-Petersburg

The paper presents an overview of the main expressions, allowing to carry out mathematical modeling of complex problems of electrochemistry in General and electrokinetic problems in particular in conditions of interdependence presented in different physical fields.

Keywords: electrokinetic phenomena, Navier—Stokes equation, equation of Nernst—Planck, equation of heat transfer, Poisson equation

Contacts: Sharfarets Boris Pinkusovich, Article arrived in edition: 24.09.2014

sharb@mail.ru