ТЕОРИЯ ЭФФЕКТА ЭКЗАЛЬТАЦИИ МИГРАЦИОННОГО ТОКА С УЧЕТОМ РЕАКЦИИ
ДИССОЦИАЦИИ - РЕКОМБИНАЦИИ ВОДЫ.
@+6
Харкац Ю.И., Сокирко А.В. Институт электрохимии им. А.Н.Фрумкина АН СССР,
117071, Москва, V 71, Ленинский проспект, 31, ( U.S.S.R. )
@+1l
Исследовано влияние процесса рекомбинации ионов ОНФ и НХ на эффект
экзальтации миграционного тока, наблюдающийся при парал- лельном
восстановлении катионов и нейтрального вещества [1,2]. Показано, что в
зависимости от соотношения концентраций разряжающихся катионов и ионов
НХ в растворе возможны качественно раз- личные типы поведения
расматриваемой системы. Обнаружено, что при определенных условиях в
кислых растворах может наблюдаться необычный предельный
диффузионно-миграционный ток разряда катионов, связанный с уменьшением
их концентрации до нуля внутри диффузионного слоя. Для случая
нейтральных растворов получены приближенные аналитические решения и
численные решения на ЭВМ соот- ветствующей системы электродиффузионных
уравнений. Показано, что для реальных значений константы равновесия
диссоциации - реком- бинации воды учет процесса рекомбинации приводит к
относительно небольшой поправке к теории эффекта экзальтации
миграционного тока, развитой в [3-6]. @+1L
1. Введение. \ -----------
Явление увеличения предельного тока реакции восстановления катионов при
одновременном восстановлении нейтрального вещества было экспериментально
обнаружено в 1936 году Кемулей и Михальским [1] и получило название
эффекта экзальтации миграционного тока. Количественное изучение эффекта
экзальтации миграционного тока было проведено Гейровским и Бурешем [2].
Исследовалось необратимое восстановление ионов NaХ из разбавленных
растворов NaCl Сначало записывалась полярографическая волна,
соответствующая разряду ионов NаХ, затем раствор насыщался кислородом из
воздуха и записывалась новая полярографическая волна. Предельный ток во
втором случае оказывался выше суммы тока восстановления кислорода и
предельного тока восстановления NаХ в отсутствие кислорода.
В работах [3-6] была предложена теория эффекта экзальтации,
учитывающая тот факт, что в результате протекающей параллельно к ^ z1+
реакции разряда катионов А1 ^ z1+
А + z1eФ ---^ AЦ (1)
реакции восстановления нейтрального вещества ( например, О2) в
диффузионном слое появляются отрицательно заряженные продукты (например,
ОНФ):
О2 + 2 Н2О + 4 еФ --^ 4 OHФ. (2) Эти вещества могут отсутствовать в
объеме раствора, тем не ме- нее, система электродиффузионных уравнений,
описывающая парал- лельно протекающие процессы, должна включать
уравнение для ука- занных анионов. Зависимость плотности предельного
тока восста- новления катионов i1 от плотности тока восстановления
нейтраль- ного вещества i2 имеет вид /4/: @+6
^ z1+z2
^
i1 = ----- FD1C¦/L + i2 D1z1 / D3z3, (3) \
\ z2
@+6
где D1, z1 и D3, z3 - коэффиценты диффузии и зарядности катионов и
продуктов восстановления нейтрального вещества, соответствен- но, z2 -
зарядность анионов, С¦ - концентация катионов в объеме электролита, L -
толщина диффузионого слоя Нернста, F - число Фарадея.
Как следует из (3), предельный ток разряда катионов линейно возрастает с
увеличением тока восстановления нейтрального вещества.
При построении теории эффекта экзальтации миграционного тока считалось,
что продукты восстановления нейтрального вещества ( в данном случае -
ОНФ ) в объеме раствора отсутствуют. Кроме того, не учитывалось
существования катионов НХ, которые могут появиться в результате
гомогенной реакции диссоциации - рекомбинации воды:
^ <---
OHФ + HХ ---^ H2O. (4) В настоящей работе проводится теоретический
анализ парал-
лельно протекающих процессов восстановления катионов и раство- ренного в
воде кислорода при произвольных рН раствора. Образую- щиеся при
восстановлении О2 анионы ОНФ поступают в диффузионный слой, где вступают
в реакцию рекомбинации с катионами НХ. Если концентрация катионов НХ в
растворе невелика, то можно ожидать, что результирующий процесс будет
весьма сходен с процессом, наблюдающимся при эффекте экзальтации
миграционного тока. В обрат- ном предельном случае,когда концентрация
ионов НХ в растворе велика, можно ожидать, что практически все ионы ОНФ
прореагируют с НХ непосредственно у поверхности электрода, так что
результирующая реакция восстановления кислорода выглядит следующим образом:
О2 + 4 НХ + 4еФ ---> 2Н2О. (5)
При этом мы имеем параллельное восстановление катионов раствора и
катионов НХ, описываемое теорией эффекта корреляционной экзальтации
миграционных токов [7].
При промежуточной концентрации катионов НХ в растворе можно
ожидать перехода от ситуации, описываемой эффектом экзальтации
миграционного тока, к ситуации, соответствующей эффекту корреляционной
экзальтации миграционных токов. @+6 @#5l
2. Общая постановка задачи. \ ---------------------------
@+6
Для простоты описания будем считать все ионы в системе однозарядными.
Полученные результаты легко распространяются на любые системы ионов,
имеющие равные зарядности.
Система электродиффузионных уравнений, описывающих парал-
лельное восстановление катионов АХ по схеме (1) и восстановле- ние
кислорода в реакции (2) с последующей рекомбинацией (4) внутри
нернстовского диффузионного слоя, имеет вид: @+6
^ dc1 dЩ i1 L
^
--- + c1 --- = ----- = j1, (5) \
\ dx dx FD1C¦
@+9
^ dc2 dЩ
^
--- + c2 -- = 0, (7) \
\ dx dx
@+9
^ d}c3 d dЩ L}CТ
^
--- - -- ( c3 -- ) = ---- ( k2c3c4 - k1/(CТ)} ), (8) \
\ dx} dx dx D3
@+9
^ d}c4 d dЩ L}c¦
^
--- + -- ( c4 -- ) = ---- ( k2c3c4 - k1/(CТ)} ), (9) \
\ dx} dx dx D4
@+7
c2 + c3 = c1 + c4. (10) Здесь с1, с2, с3, с4 - концентрации катионов АХ,
анионов,
ионов ОНФ и ионов НХ, обезразмереннные на C¦_- концентрацию ка- тионов
АХ в объеме раствора, Щ=Fж/RT - безразмерный электричес- кий потенциал;
x - безразмерная координата (0єxє) ; i1 - плот- ность тока разряда
катионов, k1- константа скорости рекомбинации НХ и ОНФ ионов; k2 -
константа скорости реакции диссоциации во- ды, остальные обозначения
стандартны. Уравнение (10) выражает условие локальной электронейтральности.
Граничные условия при x=1 ( граница диффузионного слоя с
объемом раствора ) имеют вид: с1(1)=1, с2(1)=1+k-m, с3(1)=m, c4(1)=k,
Щ(1)=0. (11)
Комбинируя уравнения (8) и (9), можно получить: @+6
^ dc3 dЩ dc4 dЩ i2 L
^
D3( --- - c3 -- ) - D4( --- + c4 -- ) = - ---- = -D3j2. (12) \
\ dx dx dx dx F C¦
@+6
Здесь i2 - плотность тока восстановления кислорода, j2 - соот-
ветствующий безразмерный поток.
Ток восстановления нейтрального вещества ( в данном случае
- О2 ) не зависит от скорости восстановления других веществ и
определяется только концентрацией О2 в перемешиваемой части раствора.
Таким образом, значение плотности тока i2 в данной задаче можно считать
заданным внешним параметром.
Значения констант k1 и k2 при нормальных условиях равны
^ . . .
[ 8 ] k1=2,5 10ФЎ сФ{ и k2 = 1,4 10{{ л} / моль} с, поэтому при условии,
что реакция (4) протекает в основном в направлении рекомбинации ионов НХ
и ОНФ, ее можно считать быстрой по сравнению со скоростью диффузионно -
миграционного транспорта ионов. Таким образом, можно считать, что внутри
диффузионного слоя поддерживается равновесие:
c3 c4 = k1 / k2(CТ)} п е}, (13)
где е} обозначает константу равновесия реакции (4), обезразме- ренную на
квадрат концентрации разряжающихся катионов в объеме раствора. При
характерном значении СТ порядка 0,1М и вышеуказанных значениях k1 и k2
значение е} - 1,3 10{а1.
Таким образом, исходная система электродиффузионных уравне-
ний свелась к системе из трех дифференциальных уравнений (6), (7), (12),
двух алгебраических соотношений (10), (13) и гранич- ных условий (11).
Эта система уравнений в общем виде не имеет аналитического решения,
поэтому для ее решения требуются допол- нительные упрощения. Так, в
излагаемом ниже разделе 3 делается предположение о том, что е=0. В
разделах 4-6 считается, что е - малая величина еа1, однако отличная от
нуля. При этом результаты раздела 4 относятся к случаю малых токов j1а1,
а результаты раздела 6 - к случаю больших токов j1>2. В разделе 5
представлено общее решение исходной системы (6), (7), (10) - (13) для
случая нейтрального раствора, а в разделе 7 описаны результаты числен-
ного интегрирования системы уравнений. @+6 @#5l
3. Эффект экзальтации миграционного тока в кислых растворах.
\ ----------------------------------------------------------- @+6
Как отмечалось в разделе 2, значение константы равновесия е}а1, и
поэтому с высокой точностью можно положить е}-0, что позволяет вместо
(13) использовать приближенное равенство с3с4-0. Это означает,что в
любой точке внутри диффузионного слоя в нулевом приближении существуют
либо ионы ОНФ либо ионы НХ, но не оба сорта ионов вместе. Поскольку мы
рассматриваем кислую среду, то нужно положить значение величины m,
отвечающей концентрации ионов ОНФ в объеме раствора, в граничном условии
(11)
равным нулю. Таким образом, весь диффузионный слой разбивается точкой
рекомбинации ионов НХ и ОНФ x=г на две области: 0єxєг и гєxє1. В
области, расположенной правее x=г, можно считать с4>0 и с3-0, а в
области, расположенной левее, наоборот, с3>0 и с4-0. Указанные
приближения позволяют перейти от уравнений (12), (13) к двум уравнениям:
@+6
^ dc3 dЩ i2 L
^
--- - c3 -- = - ----- = -j2, 0<x<г, (14) \
\ dx dx FD3c¦
@+8
^ dc4 dЩ i2 L
^
--- + c4 -- = ----- = j2о, г<x<1, (15) \
\ dx dx FD4c¦
@+6
где о=D3/D4. Отметим, что введенные в (6) и (14) безразмер- ные потоки
j1 и j2 положительны.
Таким образом, исходная электродиффузионная система уравне-
ний (6)-(9) расщепилась на две более простые системы уравнений первого
порядка. В области 0<x<г система уравнений дается формулами (6), (7),
(14) и (10) и совпадает с системой уравнений, используемой при описании
эффекта экзальтации миграционного тока. В области гєxє1 система
уравнений дается формулами (6), (7), (15) и (10) и совпадает с системой
уравнений, используемой при описании эффекта корреляционной экзальтации
миграционного тока.
В точке x=г должны быть непрерывны все концентрации с . Са-
\ й
ма точка г подлежит определению одновременно с нахождением решения
системы электродиффузионных уравнений. В соответствии со сказанным выше,
ее положение должно зависеть как от величин токов i1 и i2, так и от
параметров о и k, т.е. от величины концентрации ионов НХ в растворе.
Следует ожидать, что по мере роста последней величина x=г будет
приближаться к электроду, и при достаточно высоких концентрациях НХ,
когда реакция (1) сменя- ется реакцией (2), г станет равной нулю.
Начнем исследование со случая, когда тoчкa x=г расположена
внутри интервала (0,1). Найдем распределение концентраций с1(x), c2(x),
c3(x) и
c4(x) в каждой из областей 0єxєг и гєxє1. Распределение концентрации
с2(x) в области гєxє1 легко по-
лучается сложением и последующим интегрированием уравнений (6), (7) и
(10) с учетом граничных условий (11):
с2 = (j1 + оj2)(x-1)/2 + 1 + k. (16)
Используя (16), из уравнения (7) находим @+6
^ j1 + оj2
^
Ш = - -------------------------. (17) \
\ 2(1+k) + (j1 + оj2) (x-1)
@+6
Подставляя (17) в (6) и интегрируя получающееся линейное уравнение с
учетом (11), получаем распределение катионов с1(x) в области гєxє1: @+6
@
^ 1-j1(1+k)/(j1+оj2) j1(1+k) j1+оj2
^
с1 = ------------------------- + -------[(x-1)------ + 1]. \
\ (x-1)(j1+оj2)/ (1+k)2 + 1 j1+оj2 2(1+k) (18)
@+6
Распределение концентрации с4(x) в этой области следует из условия
электронейтральности (10) и условия с3(x>г)=0.
При x=г значения концентраций с1 и с2, даваемые формулами (16) и
(18),принимают в силу вышеупомянутого условия с3(г)= =с4(г)=0 равные
значения
c1(г)=G, c2(г)=G. (19)
При нахождении распределения концентраций с (x) в области \ й
0єxєг в качестве граничных условий при x=г будем использовать
соотношения непрерывности
c1(г-0) = c1(г+0) = G, c2(г-0) = c2(г+0) = G. (20)
Два условия (19) представляют собой систему уравнений, из ко- торых
можно определить неизвестные величины г и G:
@+6
^ ФФФФФФФФФФФФФФФФФФФ
G = А (1+k)(1 - j1k/оj2), (21)
@+6
^ ФФФФФФФФФФФФФФФФФФФФ
г = 1 - 2(1+k)[1 - А(1 - j1k/оj2)/(1+k)] / (j1+оj2). (22) Решение
системы уравнений (6), (7), (10), (14) в области (0єxєг) с граничными
условиями (10) может быть найдено способом, анало- гичным
использованному в [3], и дается соотношениями:
c1 = (j1-j2)(x-г)/2 + G, (23)
@+6
^ dЩ j1+j2
^
-- = -----------------, (24)
\
\ dx (j1-j2)(x-г) + 2G
@+5
^ (j1+j2)/(j1-j2)
^
c2 = G[(j1-j2)(x-г)/ 2G + 1] . (25) Распределения концентраций всех
компонентов и потенциала
для полной области координаты 0єxє1 даются формулами (16) -(18), (20) -
(25) и (10) при выполнении условий
оj2 Є kj1, 0 є г є 1, (26)
c (x) Є 0 при 0єxє1. (27)
\ й
Отметим, что выполнение неравенства (26) необходимо как для
существования G ( см. (21)), так и для существования г (см. (22)).
Нетрудно убедиться, что условие гє1 не накладывает никаких
дополнительных ограничений на параметры о, k, j1, j2. Определенная
формулой (22) величина г в общем случае убывает с ростом j1 при
фиксированном значении j2. Условие г=0 приводит к следующей связи j1 и j2:
j1 = оj2 ( 1 - (1+k)( 2k / оj2 - 1 )} ) / k. (28)
Учет неравенства (26) приводит к тому, что j1 в (28) может изменяться в
пределах 0єj1є2. При этом величина оj2 в (28) из- ^ ФФФФ ФФФФ
меняется от значения оj2 =2А 1+k (А 1+k - 1), когда j1=0, до значения
оj2=2k при j1=2. На плоскости (j2,j1) условию 0<гє1 соответствует
область I ( рис. 1 ), расположенная правее кривой j1(j2), определяемой
формулой (28), и ниже прямой j1 = оj2/k, ( см. (26) ), соответствующей
условию G=0.
Рассмотрим теперь подробнее ситуацию, когда реакция восстановления
кислорода идет по схеме (5), и анионы ОНФ в диффузи- онном слое
отсутствуют. При этом во всей области 0єxє1 распределения с1, с2, с4 и Щ
описываются уравнениями (6), (7),
(15) и (10) и граничными условиями (11), решения которых даются
формулами (16), (18), (10). Значения поверхностных концентраций с1(0) и
с4(0) зависят как от параметров о и k, так и от значений токов j1 и j2.
В соответствии с общепринятым определением предельного тока будем
полагать, что предельному току
^ ·
восстановления катионов j1 соответствует условие с1(0)=0. Этот ток
зависит от тока j2 параллельно протекающей реакции восста^ ·
новления кислорода. Зависимость j1 от j2 вытекает из (18) и условия
с1(0)=0, дается соотношением [4]: @#10
------- \ ·
@
^ 2 (1+k) k-оj2/j1
^ · |
j1 = --------(1 - А ------ ) (29) \ ·
\ 1+оj2/j1 1 + k
@+6
^ ·
и показана на рис. 1. Кривая j1(j2) выходит при j2=0 из точки @+5
^ · ФФФФФФФФФФФФ
j1 = 2(1+k)(1-А 1 - 1/(1+k) ), соответствующей формуле Эйкена [ ], и
при значении оj2/k=2 по- ^ ·
падает в точку Q, где j1=2, в которой смыкается с кривой (28) ^ ·
и прямой j1=оj2/k. При j1<j1 и с4(0)>0 поведение системы соответствует
области II на плоскости j1,j2 ( рис. 1 ). При постоянном значении потока
j1 и увеличении потока j2 концентрация с4(0) уменьшается ( см. рис. 2,
кривые 1,2 ). При выполнении условия @
с4(0)=0, которое, как нетрудно убедиться, совпадает с условием г=0 (
кривая 3 на рис. 2 ), мы приходим к правой границе об- ласти II,
описываемой формулой (28). Подчеркнем, что обращение с4(0) в нуль не
означает возникновения предельного тока реакции восстановления
кислорода, а соответствует переключению механизма электродной реакции со
схемы (4) на схему (2), (3). При дальнейшем росте j2 в области I точка
x=г, определяемая формулой (22), смещается вправо ( кривые 4, 5 на рис.
2 ), а в облас- ти 0<x<г распределены ионы ОНФ, генерируемые на электроде.
Возвращаясь к исследованию области I на рис.1, проанализи-
руем характер изменения концентрации катионов с1(x) при измене- нии тока
реакции их восстановления j1 и фиксированном потоке j2.
При j1=0 концентрация с1(x) линейно убывает при возрастании
x в области 0єxєг и продолжает убывать в области гєxє1. При уве-
личении потока j1 концентрация с1(x) в области 0єxєг по-прежнему
изменяется линейно, причем с1(0) и с1(г)=G убывают с ростом j1. При
достаточно больших значениях j1>оj2/(1+2k) на графике зави- симости
с1(x) в области гєxє1 появляется минимум при x=гЦ ( см. рис. 3а),
определяемом соотношением: @#6
^ ФФФФФФФФФФФФФФФФФФФФ
гЦ= 1 - 2(1+k)[1 - А(оj2-j1k) /(j1(1+k)) ] /(j1+оj2).
При дальнейшем росте j1 возможна ситуация (при j1>j2), когда с1(0)
станет меньше, чем G (см. рис. 3б). Условие с1(0)=0, где с1 дается
формулой (23), соответствует общепринятому определе- ^ ·¦
нию предельного тока разряда катионов j1 и приводит к следуюшему
соотношению: ^ ·¦ ФФФФФФФФФФФФФФФФФ
j1 = 2[1+k + J(о+1)А(1 - k/(оJ))(1+k) /(1-J)] / (1+оJ), (30)
^ ·¦
где величина J=j2/j1 ограничена неравенством J<1, которое совместно с
условием (25) реализуется в случае к<о. Зависимость ^·¦
j1 (j2), определяемая формулой (30) при к<о, показана на рис.1. ^ ·¦
Подчеркнем, что кривая j1 (j2) при к<о расположена под пря- мой j2=
kj1/о, соответствующей условию G=0 и, таким образом, ^ ·¦
кривая j1 (j2) служит верхней границей области I, в которой существует
решение с 0єгє1. Она начинается в точке Q на рис. 3.1 и ^ ·¦ ФФФФФФФФФФФ
при j2>>1 асимптотически выходит на прямую j1 =j2+ 2А(1kо)(1+k).
Распределения концентраций всех компонентов, соответствующие условию
с1(0)=0 и формуле (30),показаны на рис. 3в.
В случае к>о концентрация с1(0) не обращается в нуль ни при
каких j1 и j2. При этом увеличение тока j1 при фиксированном j2 приводит
к опережающему снижению значения с1(г)=G по сравнению с с1(0).
Одновременно значение x=гЦ,где с1(x) имеет минимум, ^ ·г
сближается с x=г,и при достаточно больших значениях j1=j1 точка гЦ
сливается с г, и с(гЦ)=G oбращается в нуль. Это происходит ^ ·г Ь
при j1пj1 и гпг:
^ ·г
j1 = оj2 / k, (31) @+6
^ Ь ·г
г = гЦ = 1 - 2 / j1 = 1 - 2k / оj2. (32) В этом случае область I
ограничена условиями G=0 и г=0.
Ситуации, когда с1(г)=G=0, соответствуют следующие из фор-
мул (16) - (18), (23) - (25) распределения концентраций ионов и
электрического поля: ^ Ь Ь
с1=(x-г)/(1-г),c2=(1+k)c1,c3=0,c4=kc1, ^ Ь Ь
dЩ/dx =-1/(x-г) при x>г; (33) ^ Ь
c1=c3=j2(1-о/к)(г-x)/2, c2=c4=0, ^ Ь Ь
dЩ/dx =-(k+о)/[(k-о)(г-x)] при x<г. (34) ^ Ь ·г
Отметим, что в области 0<x<г при j1=j1 анионы соли отсутствуют:
с2(x)=0, т.е. раствор представляет собой бинарный электролит, состоящий
из разряжающихся катионов и генерируемых на электроде анионов ОНФ.
Распределения концентраций (33) и (34) показаны на рис. 3г.
Зависимости предельного тока разряда катионов от тока реак-
ции восстановления кислорода при различных значениях k и о пока@
^ · ·
заны на рис. 4. При j2 тaких, что j1<2,зависимости j1(j2) сле- дуют
формулам теории корреляционной экзальтации миграционного ^ ·
тока. Дальнейшее поведение j1 зависит от соотношения между пара^ ·
метрами k и о. При kЄо (кривые 3,4) j1 линейно зависит от j2 в ^ ·
соответствии с (31). При 0<k<о (кривая 2) j1 определяется фор- мулой (30).
При уменьшении k/о до нуля, точка Q на рис. 4 стремится по ^ ·¦
горизонтали j1=2 к вертикальной оси.При этом кривая j1 переходит ^ ·¦
в прямую j1 =j2+2, соответствующую теории эффекта экзальтации
миграционного тока [3].
В случае, когда в объеме раствора имеется не кислая, а ще-
лочная среда, величину k в граничном условии нужно положить равной нулю:
k=0; m>0. Распределение концентраций во всем объеме диффузионного слоя
дается решением системы (6), (7), (10), (14) с граничными условиями
(11). Можно легко показать, что ^ ·
предельный ток разряда катионов j1 дается тем же соотношением, ^ ·¦
что и для нейтральной среды (k=0): j1 =j2+2,и не зависит от концентрации
ионов ОНФ в объеме раствора.
Таким образом, в рассматриваемой системе мы встретились с
необычным поведением концентрации разряжающихся катионов внутри
диффузионного слоя. Профиль концентрации с1(x) при определенных условиях
имеет минимум, что приводит к возможности обращения с1 ^ Ь
в нуль при некотором значении x=г внутри диффузионного слоя. Такое
поведение с1 обусловлено существенным влиянием на процесс подвода
катионов к электроду двух фaкторов - во-первых, параллельно протекающего
процесса восстановления кислорода и, во-вторых, гомогенной реакции
рекомбинации ОНФ и НХ ионов внутри диффузионного слоя.
Обобщенное условие возникновения в рассматриваемой системе предельного
тока имеет вид: @
с1(x) = 0, 0 є x < 1, (35)
и, как показало проведенное исследование, может реализоваться как у
поверхности электрода x=0, так и "оторваться" от электрода ^ Ь
и перейти внутрь диффузионного слоя в точку x=г<1. Отметим, что при о_1,
т.е. D3_D4, в области x-г, где проте-
кает реакция рекомбинации ОНФ и НХ ионов, происходит резкое из- менение
величины электрического поля Ш и dШ/dx.При этом, как следует из
уравнения Пуассона,в указанной области должны наблю- даться аномалии в
распределении объемного заряда.В частности, при стремлении к нулю
величины G=с1(г), внутри диффузионного слоя формируется резкий пик
объемного заряда. Возможность формирования пика объемного заряда при
прохождении электрического тока ранее была предсказана для систем,
имеющих наряду с подвижными также и фиксированные заряды [10-11]. @#5l
@+6
4. Поправки к диффузионно-миграционному току восстановления
\ ------------------------------------------------------------ @
катионов j1 в случае малых токов j2. \ ------------------------------------
@+6
Исследование, проведенное в предыдущем разделе, существенно использовало
условие, что скорость рекомбинации НХ и ОНФ ионов весьма велики, и,
следовательно, основная реакционная зона расположена вблизи точки x=г.
Кроме того, использовалось приближение малой константы диссоциации
е}_-0. Это привело к тому, что в точке г концентрации с3, с4 стремятся к
нулю, а остальные концентрации испытывают излом. В более строгой
постановке, использующей соотношение (13), снижение концентраций
происходит не до нуля, а до весьма низких значений порядка е, а профили
всех концентраций сохраняют плавность. Хотя при этом все обнаруженные
особенности в распределении ионных компонентов сохраняются, для
получения более полной картины представляет интерес получить
распределение концентраций при условии малого, но конечного значения еа1.
В этом и последующих разделах мы ограничимся только случаем нейтрального
раствора с3(1)=с4(1), откуда с помощью (13) можно @
записать граничные условия (11) в виде: с1(1) = с2(1) = 1,
с3(1)=с4(1)=е, Щ(1)=0. (37)
Кроме того, мы будем интересоваться случаем предельного тока по катионам:
с1(0) = 0. (38) Сначалa найдем решение, соответствующее j2=0. Подставляя
(13) в (12), после несложных преобразований получаем, что в этом случае
отсутствует как перенос катионов НХ, так и перенос ОНФ:
@+6
^ d c3 d Щ
^
---- + c3 --- = 0, (39)
\
\ dx dx
@+9
^ d c4 d Щ
^
---- - c4 --- = 0. (40)
\
\ dx dx
@+7
Складывая уравнения (7) и (39) и вводя суммарную концентрацию анионов,
мы приходим к задаче Эйкена [9]. Распределение концен- траций показано
на рис 5. Потенциал Щ ведет себя логарифмическим
образом в зависимости от расстояния x:
@+6
^ Щ
е = 1 + j1(x-1)/(1+е), (41)
а концентрации с2 и с3 линейно убывают при приближении к электроду:
^ Щ Щ
c2 = e , c3 = е e , (42)
^ Ф
однако не принимают нулевых значений на электроде: с2(0)-Ае, c3(0) -
е~'}. Отметим, что для рассматриваемой системы это является общим
правилом: концентрации всех компонентов, кроме с1, не могут одновременно
обратиться в нуль на электроде в режиме предельного тока в силу условия
(11). Это приводит к тому, что градиент потенциала и сам потенциал, как
следует из (7), также остаются конечными. Таким образом, в
рассматриваемой системе, учитывающей диссоциацию, отсутствует
логарифмическая расходимость потенциала, обычно встречающаяся в
электродиффузионных задачах при расчете предельного тока.
^ -Щ
Поскольку из (13) и (42) следует, что с4 = е е , то практически во всем
диффузионном слое выполняются неравенства
с3 << с4, (43)
с3, с4 << с1, с2. (44)
Значение предельного тока дается формулой Эйкена, записываемой в виде:
^ ФФФФФФФ
j1 = 2(1+е)(1- Ае/(1+е) ). (45)
Найдем теперь решение в случае малых, но отличных от нуля значений j2,
таких, что неравенства (43), (44) по-прежнему выполняются. Пренебрегая в
(12) с3 по сравнению с с4, можно привести (12) к уравнению:
@+6
@
^ d c4 d Щ
^
---- - c4 --- = - о j2. (46)
\
\ dx dx
@+6
Вводя суммарную концентрацию анионов, мы сводим задачу к задаче о
разряде двух сортов катионов с1 и с4 при неподвижном фоне анионов, т.е.
к задаче о корреляционной экзальтации миграционного тока, решенную в
[4,7]. Соответствующая зависимость для предельного тока j1 oт j2 имеет вид:
@+6
^ ФФФФФФФФФФФФФФФФФФФФФФФФФ
j1 = 2(1+е) - оj2/2 - А 4(1+е)(е-оj2/2) + (оj2)}. (47)
Соответствующие формуле (47) зависимости j1(j2) при j2а1 для ряда
значений параметра е приведены на рис.6. При j2=0 формула (47) переходит
в (45). Подчеркнем, что в соответствии с @
(45), (47) предельный ток j1 снижается при учете процесса рекомбинации
ОНФ и НХ ионов на величину, пропорциональную е{'}.
Перейдем к описанию общего решения полной системы уравне-
ний, не использующему предположения о малости е.
@+6 @#5l
5. Общее решение задачи в случае j2_0.
\ -------------------------------------
@+6
Нашей целью является расчет значения предельного тока j1 из системы трех
дифференциальных уравнений (6), (7), (12), двух алгебраических
соотношений (10) и (13), а также граничных условий (38), (39).
Подставляя в (10) и (12) условие (13), можно исключить из них
концентрацию с4:
@+7
^ d c3 d Щ е}
^
---- - c3 --- = - j2 / ( 1 + ----- ), (48)
@
\
\ dx dx о c3}
@+7
c1 + е}/c3 = c2 + c3. (49)
При е--^ 0 эта система уравнений переходит в уравнения из [5],
соответствующие отсутствию учета диссоциации воды, решение которых
приводит к выражению (3) для предельных токов.
Правая часть уравнения (48) нелинейным образом зависит от с3, и поэтому
систему уравнений (6), (7), (48), (49) не удается свести к линейной
путем перехода к новой независимой переменной Щ, как это было сделано в
[4]. Таким образом, поставленная за- дача требует разработки новых
методов интегрирования рассматри- ваемой системы нелинейных
дифференциальных уравнений.
Исключая dЩ/dx из (5) и (48) с помощью (7), получаем: @+7
^ d w
^
--- = j1 c2, (50)
\
\ d x
@+9
^ d v j2 c2
^
--- = ----------, (51)
\
\ d x c3} + е}/о
@+7
где введены
w = c1 c2, v = c2 / c3 (52)
- безразмерные комбинации концентраций. Выражая из (52) с1 и с2 через
с3, w, v и подставляя в (49), можно выразить c3 через w, v:
@+7
^ w/v + е}
^
c3} = --------. (53)
\
\ 1 + v
@+6
Таким образом, с помощью (49), (52) и (53) концентра- ции всех
компонентов раствора выражаются через w,v, и в систему
(50), (53) входят только функции w и v, однако довольно сложным образом.
Граничные условия для w, v приобретают вид: w| =1, v| =1/е, (54)
\ x=1 x=1
@+6
w| =0. (55) \ x=0
Поделив (50) на (51), перейдем к новой независимой перемен- ной v:
@+7
^ d w J 1 1
^
--- - ------ w = J е} ( --- + --- ), (56)
\
\ d v (1+v)v 1+v о
@+7
где J=j1/j2 - отношение токов. При выводе (56) было использо- вано
соотношение (53). Уравнение (56) - линейное неоднород- ное уравнение
первого порядка. Его решение, удовлетворяющее при х=1 граничным условиям
(54), записывается в виде: @+7
^ v 1 1 dv'
^
w(v) = wк(v) [ 1 + Ы J е} ( ---- + --- ) ------ ], (57) \
\ 1+е 1+v' о wк(v')
@+7
где wк(v) - решение однородного уравнения (54), удовлетворяю- щее
граничным условиям (54): @#12
\ J
^ (1+е) v
^
wк(v) = [ ------- ] . (58) \ 1 + v
@+7
Таким образом находится зависимость w(v), с помощью которой, используя
(52) и (53), можно получить распределение концентраций всех компонентов
как функций независимой переменной v при условии, что известно значение
J. Для нахождения J необ-
ходимо использовать граничное условие (55). Интегрируя (51), получаем:
@+7
^ 1 1/е dv 1/е dv е}
^
j2 = j2 Ы dx = Ы ---- c3 + Ы ---- ----. (59)
\
\ 0 vЦ qЯЯЯЯЯЯЯЯЯЯЯ
AЯЯЯ_ЯЯ_ИЮ_ЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯрс?АБЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯЯррp____ЯЯ
AЯЯАЯЮ_ ___ __ЯЯ_ГЯЮ_ ________ ___Я_АПГЯЮ___ГЯЮ_ _ _@< _рс(55), w(vЦ)=0.
Поскольку, как легко видеть из (58), w (v) при положи-
тельных v в нуль не обращается, условием для нахождения vЦ явля-
ется обращение в нуль квадратной скобки в (57):
@+7
@
^ 1/е 1 1 dv' 1
^
Ы ( ---- + --- ) ------ = ---. (60) \
\ vЦ 1+v' о wк(v') 1+е
@+7
Таким образом, в (59) можно найти пределы интегрирования (
из (60), (58)) и подинтегральную функцию ( из (53), (57), (58)).
Произведя интегрирование в правой части (59), мы получаем в правой части
некоторую функцию отношения токов J. Наконец, можно выразить из
получившегося уравнения j1 как функцию j2, получив обобщение формулы (3)
на случай экзальтации миграционного тока, учитывающей рекомбинацию ионов
НХ и ОНФ.
Описанная процедура позволяет получить решение поставленной
задачи в квадратурах. Для определения зависимости j1 oт j2 необходимо
аналитическое вычисление входящих в (57), (59), (60) ин- тегралов. Эта
задача разрешима для интеграла (57) только при не@
которых специальных значениях J. ( Интеграл (57), входящий в (59),
берется только при целых и полуцелых J, а интеграл в пра- вой части (59)
в аналитическом виде не вычисляется.) Ниже будет приведено приближенное
решение задачи, использующее малость па- раметра е. @+6 @#5l
3.5 Приближенное решение в случае 1<J<2.
\ ----------------------------------------
@+6
Для того, чтобы найти приближенное аналитическое решение задачи, можно
воспользоваться малостью параметра е: е<<1. Нашей целью будет нахождение
поправочного члена к правой части выражения (1.46), являющегося функцией
как е, так и j2. Рассматривая (1.46) в качестве нулевого приближения,
можно отметить, что условие J<2 соответствует токам j2>2. Область J>2
или 0<j2<2 соответствует переходной области от j2=0 к j2>2, где
поправочные ^ Ф
@
члены имеют промежуточное значение между членами порядка Ае при j2=0 и
членами порядка е при j2>2, как это будет показано ниже. Отметим, что
условие J>1 или j1>j2 выполняется вплоть до очень больших значений j2Ь
1/еЄ10~, т.е. во всем интересующем нас диапазоне параметров j2. При
v>>vЦ можно записать с точностью до членов порядка е: ^ J
w(v) = [ 1+еJ(1-1/о)] [ v / (1+v) ] , (61) а при v>vЦ ^ J-1 J-1
w(v) = [ 1+еJ(1-1/о)] (1-Jv) v [ v - (vЦ) ], (62) где
^ J-1
(vЦ) = Jе} ( 1 + 1/о ) / ( J - 1 ) - (63) - малая величина. Переходя в
(59) к интегрированию по ln v, заменяем пределы интегрирования на
ln(1/е) и ln(vЦ), которые, в отличие от самих величин 1/е и vЦ, можно
считать величинами по@
рядка единицы. Таким образом, область интегрирования имеет размер
порядка единицы, что подтверждает правомерность сохранения в
подинтегральной функции величин порядка не ниже, чем е, поскольку нас
интересует только первая поправка к (3) по е. Подставляя (61) и (63) в
правую часть (59) и производя интегрирование, по- лучаем уравнение,
задающее в неявном виде зависимость j1(j2,е): @+6
^ 2 е
^
j1 = ----- - ----- f(J), (64) \
\ J - 1 J - 1
@#12
где
\ ФФФФФФФФФ ФФФФФФФФФФФ ФФФФФФФФФФФ
^ | J(1+1/о) |J(1+1/о) |J(1+1/о)
^
f(J) = - --------+1/о-1-(--2arcctg -------- -1)/ -------- -1 \ А А А /о
\ J - 1 J - 1 J - 1
(65) Первое слагаемое в правой части (64) отвечает нулевому приближению
(3), а второе представляет собой малую поправку. Поэтому уравнение (64)
можно разрешить относительно j1 методом последо- вательных приближений.
Сначала находим J в нулевом приближении:
JТ = ( j2 + 2 ) / j2. (66)
Умножая (64) на J-1 и подставляя в правую часть получившегося выражения
JТ вместо J, получаем:
j1 = 2 + j2 - еf(JТ) - (67)
- итоговое выражение, являющееся обобщением ранее полученого результата
(3). Определяемая формулой (67) зависимость j1(j2) по- казана на рис. 7,
где показана также зависимость (37).
@+6 @#5l
7. Численное решение задачи на ЭВМ.
\ -----------------------------------
@+6
Система уравнений (50), (51) с граничными условиями (54), (55) решалась
также численным интегрированием на ЭВМ. Алгоритм решения строился
следующим образом. При фиксированных j1, j2 система (50), (51) с
граничными условиями (54) интегрировалась методом Рунге-Кутта.При
интегрировании могли реализоваться два варианта: либо на некотором шаге
по координате x значения w или v становились отрицательными (детальный
анализ показывает, что это происходит одновременно),либо выполнялось
соотношение w(x=0) >0.При фиксированном j2 реализующийся вариант
определялсь только значением j1. Таким образом, можно говорить о
существовании не- которой функции от j1, корню которой соответствует
w(x=0)=0. Для нахождения корня этой функции использовался метод
половинного деления. Сопоставление приближенных аналитических
зависимостей j1(jк) (67) с численными расчетами для различных значений е
по- казывают их неплохое согласование. На рис.8 приведены профили
концентраций компонентов с , полученные при численном решении \ й
задачи.
@+6
8. Заключение.
\ --------------
@+6
Таким образом, проведенный анализ показал, что эффект диссоциации ионов
ОНФ и НХ проявляется в нейтральных растворах наиболее существенно при
j2<<1, когда предельный ток восстановления ^ Ф
катионов снижается на величину порядка Ае. При больших токах j2> 2
снижение предельного тока восстановления катионов имеет поря- док е.
Численные расчеты в промежуточной области е<j2<2 показа- ли, что
снижение предельного тока имеет порядок промежуточный ^ Ф
между Ае и е и плавно зависит от величины тока восстановления
нейтрального вещества.
В случае протекания рассматриваемых процессов в кислых растворах процесс
рекомбинации ионов ОНФ и НХ проявляется более существенно, и при
достаточно высоких рН результирующий экзальта- ционный ток приближается
к значению, описываемому теорией эффекта корреляционной экзальтации
миграционного тока [4]. При этом обнаружено, что при достаточно больших
концентрациях в растворе НХ может наблюдаться необычный предельный
диффузионно - миграционный ток разряда катионов, связанный с уменьшением
их концент- рации до нуля внутри диффузионного слоя.
@:
@
Литература.
\ ----------
@+6
1. Кemula W., Michalski M.// Roczniki Chem., 1936, v. 16, N 3, p. 535 -541.
2. Heyrovsky J.,Bures M.// Coll. Czech. Chem. C