К теории скрытых предельных диффузионных токов.
Сокирко А.В., Харкац Ю.И.
В работе проведено теоретическое исследование параллельно про-
текающих электродных процессов в условиях, когда продукт одной из
электродных реакций вступает в гомогенную реакцию с реагентом другой
электродной реакции. В результате предельный ток второй реакции зависит
от тока первой реакции и наблюдается явление, названное Кемулей и
Грабовским "скрытым предельным током". Проведенные аналитические и
численные расчеты охватывают широкую область параметров задачи и
включают , как частный случай, ре- зультаты Кемули и Грабовского для
предела бесконечно быстрой константы скорости гомогенной реакции.
Явление скрытых предельных токов, широко распространенное в
вольтамперометрии, связано с особенностями протекания параллельных
электрохимических реакций в том случае, когда продукт одной из реакций
необратимо взаимодействует в растворе с реагентом дру гой реакции /1-3/.
В результате такого взаимодействия диффузионный поток реагента во второй
реакции зависит от скорости протека ния первой реакции, т. е. процессы
уже не являются независимыми.
Пусть, например, в растворе присутствуют два вещества А и В,
которые восстанавливаются на электроде, причем для определенности будем
считать, что В восстанавливается при более отрицатель- ных потенциалах.
A + ne --> pA , (1) B + me --> rB . (2)
Eсли в растворе протекает реакция pA +qВ --> A B , (3)
то по терминологии работы /1/ наблюдается скрытый предельный ток первого
рода вещества В, определяемый как разность предельных токов в отсутствии
вещества А i и в его присутствии i .
Если же в растворе протекает реакция
то наблюдается скрытый предельный ток второго рода i =i -i . Проведенный
в /1,2/ теоретический анализ скрытых предельных
токов относится к случаю, когда реакции (3) или (5) протекают бесконечно
быстро по сравнению со скоростью диффузионного транспорта реагентов и
продуктов реакции.
В настоящей работе развивается теория скрытых диффузионных
токов, не использующая указанного выше предположения. Мы ограничимся
расмотрением скрытых предельных токов пер-
вого рода ( задача о скрытых предельных токах второго рода решается
аналогичным образом) и кроме того, будем для простоты полагать p=q=1.
Считая, что рассматриваемые процессы протекают в условиях избытка
фонового электролита, когда эффекты электромиграции несущественны *),
уравнения, описывающие распределение концентрации вещества В и продукта
реакции (1) А в диффузион- ном слое, можно записать в виде: @+4L
-------------------------------- *) Взаимодействие двух кинетически не
связанных электродных процессов за счет миграционного транспорта ионов
приводит к эффекту экзальтации миграционного тока, теория которого
изложена в /4-7/
Здесь D и D - коэффиценты диффузии соответствующих веществ, K -
константа скорости реакции (3), - - координата, перпендику-
лярная поверхности электрода. В качестве граничных условий будем
считать, что заданы кон-
центрации C в об'еме раствора
@+1L
и нулевая концентрация вещества А @+1L
Здесь д - толщина диффузионного слоя, которая считается прибли- женно
одинаковой для веществ А и В. При достижении предельного тока по
веществу В выполняется условие C (0)=0. Кроме того, на электроде задано
значение j =-j , где j - поток вещества А, т.е. @+3L
Перейдем к безразмерным переменым : @+6L
Тогда система уравнений (6)-(7) и граничные условия к ним могут быть
записаны в виде : @+7L
Нашей целью является нахождение связи между Y(0)=dY/dt... - безразмерным
предельным током вещества В и л- безразмерным пре- дельным током
вещества А при различных значениях параметра в .
Система (12)-(15) в общем случае не имеет аналитического ре-
шения. Ниже будут изложены результаты исследования ряда предельных
случаев, для которых можно получить приближенные решения задачи
(см.Таблицу), а также результаты численного решения на ЭВМ.
1. Случай малых скоростей реакции ( в << 1 ).
\---------------------------------------------- Будем считатать, что
поток л достаточно велик (л >>1), так
что лвЬ 1. При этом функция Х имеет характерный масштаб ХЬл , а функция
Y¦1. Это позволяет искать решение задачи в виде разложения по параметру
1/л, первые члены которого имеют вид X=лX +X , Y=Y , где X , X , Y -
функции с масштабом порядка единицы.
Из (12)-(13) с использованием граничных условий (14)-(15)
можно получить линейную связь между и : @+1L
Подставляя и , выраженные через , , и , и рассматривая члены порядка л,
получим . Используя (20)-(22), уравнение (7) можно привести к виду @+2L
Решением этого уравнения служит линейная комбинация функций Эйри @+1L
где ; и - константы. Определив и из граничных условий (15)-(16) и
дифференцируя (19), получаем исkомое выражение для предельного потока
вещества В: @+2L
В случае достаточно малых в, когда q<<1, можно воспользоваться ведущими
членами асимптотических разложений функций Эйри /8/. В результате (20)
приводится к виду @+1L
Формулу (21) можно получить и непосредственно из исходной системы
(12)-(15), представляя решение в виде разложения по малому параметру в :
X = X + вX , Y= Y + вY, где X , X , Y , Y функции порядка 1 и считается,
что вл<<1. Поведение X(t), Y(t), oпреде- ляемых (19) и (17) показана на
рис.1 Таким образом, выражение (20) применима при в<<1 и любых л.
2. Случай больших потоков вещества А ( л >>1 ). \
----------------------------------------------
В отличиe от случая малых скоростей реакций в<<1, влЬ1 при котором
функции X(t), Y(t) изменялись плавно во всем интер- вале 0єtє1, в случае
вЄ1 и вл >>1 при значениях t-1 функция Y(t) изменяется весьма резко, а
на всем остальном интервале из- меняется плавно, принимая значения,
близкие к нулю. Другими словами , у границы диффузионного слоя t=1
Формируется пограничный слой.
Поскольку в рассматриваемом случае Y(0)=е<<1 связь X с Y (17) может быть
записана в виде: @+1L
7. Сокирко А.В., Харкац Ю.И. Диффузионно-миграционные токи в параллельно
протекающих процессах электроосаждения металлов и восстановления
анионов. - Электрохимиия ( в печати ).
@:
Подпись к рисунку. \ -----------------
Области различных видов решений в плоскости параметров (о,л)
для z1=2, z2=z3=1: I - область линейных профилей концентраций (формулы
(21)),
II - область нелинейных профилей концентраций (см. (30)-(33)), III -
область полной депрессии тока j1 ((49)-(52)); кривые 1, 2 отвечают
формулам (55), (56), соответственно.
3. Случай малых потоков вещества А ( л << 1 ) .
\ ---------------------------------------------
При л << 1 и в достаточно малых, так что вл << 1 поток ве-
щества дается формулой (21). Ниже мы рассмотрим случай вЬ1,в>>1.
Поскольку характерный масштаб изменения функции X(t)Ь л <<1,
будем искать X,Y в виде разложений по параметру л, первые члены которых
имеют вид: @+1L
где X , Y , Y - функции с масштабом изменения порядка 1. Подставляя (25)
в (13) получим с учетом граничных условий Y =t. При этом уравнение (12)
принимает вид d X /dt =вtX. Его решени- ем, удовлетворяющим (15), служит
линейная комбинация функций Эйри
@+2L
Используя теперь линейную связь X и Y (17) находим поток В: @+4L
Используя ведущие члены асимптотических разложений функций Эйри, можно
убедиться, что при в<<1 (27) переходит в (21).
При больших скоростях реакции в>>1 решению задачи соответству-
ет существование пограничного слоя вблизи t=0. При этом Х(1)<<1 и
решение задачи находится по схеме, аналогичной описанной выше для случая
погранслоя при t=1. Получающийся результат совпaдает с выражением для
потока вещества В, следующим из (27) при в >>1. @+2L
Таким образом, формула (27) справедлива при л<<1 и любых зна-
чениях в .
4. Большие скорости реакции ( в >>1 ) . \
------------------------------------
При расмотрении этого случая будем различать две ситуации л < 1 и л > 1.
Пусть л<1. При больших скоростях реакции в>>1 практически
все вещество А успевает прореагировать с веществом В вблизи электрода. У
поверхности электрода формируется пограничный слой, в котором Y(t) имеет
резкий максимум, так как Y(0)=0, а за максимумом Y(t) экспоненциально
убывает, в силу того, что все вещество А прореагировало около электрода.
Именно такой предельный случай и соответствует рассмотрению скрытых
предельных токов, проведенному в /1/ , где считалось, что в пределе в=..
поток вещества В на электрод равен 1-л.
Будем считать, что при в >>1 поток вещества В на электрод
отличается от 1- л на малую величину е : Y(0)= 1 + е - л .
Тогда с учетом (15) связь Х и Y дается выражением
@+1L
где t =е/(1+е)-е . В пределе в--> ,е=0, Y(t)=t, X(t)=0 Будем искать
решение в виде
@+3L
При t>>t в (30) можно пренебречь Х и получить уравнение для Х
@+2L
решение которого, удовлетворяющее (15), дается выражением: @+3L
Отметим, что при t -t второе слагаемое в выражении для Х эк-
споненциально мало по сравнению с первым и им можно пребречь.
Соответствующая Х функция Y находится из (30) и равна: @+2L
В области tєt, Y Ь t , а X Ь е-t . Поэтому, пренебрегая в (30) Х по
сравнению со вторым слагаемым, приходим к следующему уравне- нию для Y @+1L
Решением этого уравнения служит линейная комбинация функций Эйри @+2L
Связь S и S находится с помощью граничных условий (15) и усло- вий
сшивки решений. Условия сшивки для Y и Y имеют вид: @+2L
Определив коэффиценты S ,S ,S из (29) и (31), получаем уравнение для е :
@+2L Определяемое уравнением (32) значение е имеет при значениях л не
близких к единице порядок в . В частном случае л<<1 из (32),(29) можно
получить формулу (28). Легко убедиться,что при 1-л<<1 Аве Y(0)=1-л << 1.
Таким образом, (29) дает правильное решение при малых потоках вещества А
л<<1 и при всех значениях л<1 и в>>1.
Обратимся теперь к ситуации в>>1,л>1. В отличиe от рассмотрен-
ного выше случая л>>1, пограничный слой возникает здесь не на правой
границе области t=1, а в некоторой промежуточной точке внутри области
0<t<1.
Будем искать решение задачи в виде первых членов разложения
по малому параметру в : @+1L где X , X , Y , Y Ь1. В нулевом приближении
имеем X Y =0, везде, кроме окрестности точки t =1-1/л, этому уравнению и
граничным условиям удовлетворяют функции @+3L
В окрестности t=t сосредоточена реакционная зона, где веще-
ство В вступает в реакцию с веществом А . Лишь небольшая часть вещества
В доходит до электрода: Y(0)=е<<1. Для нахождения е будем искать решение
в виде, следующем из (33) и (34): @+3L где Y , X Ь в << 1 . Схема
нахождения решения полностью анало- гична описанной выше. Достаточно
далеко от точки t=t можно пребречь в (30) Y или X , получить в каждой из
двух областей уравнение Эйри и найти коэффиценты в линейной комбинации
функций Эйри из граничных условий и условий сшивки функций Y и Y и их
произ- водных при t=t . В результате получим @+2L При получении (36)
считалось, что л-1 >> 1, л<<в т.е. точка t расположена не слишком близко
к концам интервала t >>е, 1-t >> е Если же считать, что 1-t << е, то
можно убедиться, что Y (0) дается формулой (24). В случае t <<1, Y(0) - 0.
Профили концентраций X(t), Y(t), для случая в>>1, лЬ1 пока-
заны на рис. 3.
5. Численное решение задачи. \ ------------------------
Для получения полной картины поведения предельного потока вещества В при
произвольных значениях параметров в и л система (12)-(15) решалась
численно.
Схема решения строилась следующим образом. Из соотношения
(17) при t=0 следует @+1L Подставляя (37) в (17), имеем @+1L Используя
(38), преобразуем (12) к виду @+1L Уравнение (39) с граничными условиями
X(1)=0, X(0)=з решалось численно методом линеаризации и прогонки
/9/-/10/. Из решения находилось значение X(0)=-л и Y(0)- предельный
поток вещества В, которые являлись функцией з , и, таким образом, в
параметрическом виде искомая зависимость Y(0)(л,в). Расчитанные
зависимости полного потока УJ=л+Y(0) на электрод А и В для ряда значений
параметра в представлены на рис. 4. Семейство кривых УJ показывает
степень влияния величины скорости реакции в на результирующую кинетику
процесса. Рассмотренному в /1-3/ пределу в = . соотве- тствует УJ =1 при
л<1 и л при л>1. Скрытый предельный ток вещества В , равный разности
токов при в=0 и в=0, дается соотношением ДJ=1-Y(0) и показан как функция
л на рис. 5.
6.Заключение. \ ------------
Проведенные расчеты зависимости величин скрытых предельных токов от
параметра л для широкого набора значений в дают возможность оценивать по
соответствующим экспериментальным зависимос- тям значение константы
скорости гомогеннной реакции К. Для ее нахождения нужно сопоставить
экспериментальную кривую зависимо-
сти суммарного тока от состава раствора с расчитанными кривыми и выбрать
то значение в, которое соответствует соответствует тео- ретической
кривой, наиболее точно описывающей экспериментальные данные.
Авторы благодарят Л.Г.Феоктистова за полезное обсуждение работы
@:
Литература.
1.Kemula W.,Grabowski Z.R. Roszn. Chem., 1951, v.25, p. 350-366.
2.Феоктистов Л.Г., Жданов С.И., Изв. АН СССР, отд. хим. наук
1963, с. 45-52.
3.Галюс З. Теоретические основы электрохимического анализа. М.,Мир,
1974, 552 с.
4.Харкац Ю.И., Электрохимия, 1978, т.14, N 11, с. 1716-1720. 5.Харкац
Ю.И., Электрохимия, 1978, т.14, N 12, с. 1840-1844. 6.Kharkats Yu.I., J.
Elektroanal. Cnem., 1979, v. 105, N 4
p. 97-114.
7.Тополев В.В., Харкац Ю.И., Электрохимия, 1983, т. 19, N 4,
с. 515-520.
8.Справочник по специальным функциям. Ред. М.Абрамовиц, И.Сти- ган, М.,
Наука, 832 с.
9.Самарский А.А. Теория разностных схем. М., Наука, 1977, 656 с.
10.Дулан Э., Миллер Дж., Шилдерс У., Равномерные численные методы
решения задач с пограничным слоем. М., Мир, 1983, 198 с.
@:
Подписи к рисункам:
.
Рис. 1. Зависимость X(t), Y(t) и в X(t) Y(t) при X(0)=1 , в=1 Рис. 2.
Зависимость X(t), Y(t) и в X(t) Y(t) при X(0)=8 , в=10 Рис. 3.
Зависимость X(t), Y(t) и в X(t) Y(t) при X(0)=1.5,в=1000 Рис. 4.
Зависимость УJ от л для значений параметра в :
1 - 0 ; 2 - 1 ; 3 - 10 ; 4 - 100 ; 5 - 1000 ; 6 -
Рис. 5. Зависимость скрытого предельного тока от л для значений
параметра в: 1 - 0 ; 2 - 1 ; 3 - 10 ; 4 - 100 ; 5 - 1000 ; 6 - @:
Сводная таблица приближенных аналитических формул для зависи-
мости потока вещества В Y(0) от л для различных областей значе- ний л и
параметра в.
@+3L
------------------------------------------------ | в | в<<1 | в Ь 1 | в
>> 1 |
| л | | | |
|----------------------------------------------- | л<<1 | (21) | (27) |
(28) |
|----------------------------------------------- | | | Численные ре- |
лє1 (29),(32) |
| лЬ1 | (21) | шения на ЭВМ. | лЄ1 (36) |
|----------------------------------------------- | | | | (36) л <<в |
| л>>1 | (20) | (24) л >>в>1 | (24) л >>в |
------------------------------------------------ Стрелка показывает, что
соответщюх ШЭ@ А
Д` А "пА ;лицы для различных областей