Научная статья на тему 'АЛГОРИТМ ЧИСЛЕННОГО ИССЛЕДОВАНИЯ ПРИКАТОДНЫХ ПРОЦЕССОВ ЭЛЕКТРОВЫДЕЛЕНИЯ ЦИНКА'

АЛГОРИТМ ЧИСЛЕННОГО ИССЛЕДОВАНИЯ ПРИКАТОДНЫХ ПРОЦЕССОВ ЭЛЕКТРОВЫДЕЛЕНИЯ ЦИНКА Текст научной статьи по специальности «Математика»

CC BY
3
1
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЭЛЕКТРОВЫДЕЛЕНИЕ ЦИНКА / ПРОГНОЗИРОВАНИЕ / ЧИСЛЕННАЯ ОПТИМИЗАЦИЯ / СТАДИЙНЫЕ ЭЛЕКТРОХИМИЧЕСКИЕ РЕАКЦИИ / МЕТОД РУНГЕ-КАТТА / ZINC ELECTROEXTRACTION / PREDICTING / NUMERICAL OPTIMIZATION / STEPWISE ELECTROCHEMICAL REACTIONS / RUNGE-KUTTA METHOD

Аннотация научной статьи по математике, автор научной работы — Кашапов Наиль Фаикович, Чебакова Виолетта Юрьевна, Чебакова Евгения Вячеславовна, Кормушин Константин Владимирович

Предмет исследования: в работе представлен чис ленный алгоритм нахождения констант скоростей гетерогенных процессов для последующего их исследования. Цель исследования: Цифровое прототипирование электроосаждения цинкового покрытия путем выщелачивания обедненной цинковой руды. Методы исследования: алгоритм представлен на примере задачи прогнозирования выделения цинка при параллельной реакции выхода водорода при электролизе раствора гидроксида натрия, содержащего тетрагидроксоцинкат натрия и верифицирован на натурных экспериментах. Так, в данной статье методом прямого поиска находятся константы прикатодных стадийных электрохимических реакций получения порошка цинка. Сама же система кинетических уравнений, описывающая стадийные прикатодные реакции, решается классическим методом Рунге-Кутта. Основные результаты исследования: моделирование учитывает наличие параллельного процесса выделения водорода и происходит в приближении пористой пленки при образовании осадка. Численный алгоритм, разработанный авторами, также позволяет прогнозировать выход вещества по току.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по математике , автор научной работы — Кашапов Наиль Фаикович, Чебакова Виолетта Юрьевна, Чебакова Евгения Вячеславовна, Кормушин Константин Владимирович

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ALGORITHM FOR THE NUMERICAL INVESTIGATION OF CATHODE PROCESSES OF ZINC ELECTROEXTRACTION

Subject of research: the paper presents a numerical algorithm for detecting the rate constants of heterogeneous processes for their further investigation. Purpose of research: Digital prototyping of electrodeposition of zinc coating by leaching depleted zinc ore. Methods of research: the algorithm is presented on the example of predicting the dezincification during parallel hydrogen release reaction in electrolysis of a sodium hydroxide solution, which contains sodium tetrahydroxozincate, and verified using fullscale experiments. So, in this article, there are detected the constants of cathode stepwise electrochemical reactions for zinc powder obtaining using direct search method. The kinetic equation system itself, which describes cathode stepwise reactions, is solved using classical Runge-Kutta method. Main results of research: the modeling makes allowance for the presence of parallel process of hydrogen release and occurs in the approximation of porous film during the settling-out. The numerical algorithm, developed by the authors, make it also possible to predict the current yield of

Текст научной работы на тему «АЛГОРИТМ ЧИСЛЕННОГО ИССЛЕДОВАНИЯ ПРИКАТОДНЫХ ПРОЦЕССОВ ЭЛЕКТРОВЫДЕЛЕНИЯ ЦИНКА»

0 МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ I И ИН1<Ф>01>1А|А|1ДИ01Н1 1Н11ЫЕ ТЕХНОЛОГИИ

УДК УДК 519.168 DOI 10.18822/byusu20240111-18

АЛГОРИТМ ЧИСЛЕННОГО ИССЛЕДОВАНИЯ ПРИКАТОДНЫХ ПРОЦЕССОВ ЭЛЕКТРОВЫДЕЛЕНИЯ ЦИНКА

Кашапов Наиль Фаикович,

доктор технических наук, профессор, профессор кафедры Физики перспективных технологий и материаловедения,Казанский (Приволжский) федеральный университет, Заместитель руководителя, «Федеральный исследовательский центр Казанский научный центр Российской академии наук» Казань, Россия E-mail: kashnail@gmail.com

Чебакова Виолетта Юрьевна,

кандидат физико-математических наук,

доцент Института вычислительной математики

и информационных технологий,

Казанский (Приволжский) федеральный

университет

Казань, Россия

E-mail: vchebakova@mail.ru

Чебакова Евгения Вячеславовна,

Химический институт им. А. М. Бутлерова, Казанский (Приволжский) федеральный университет Казань, Россия

E-mail: evgeneachebakova@mail.ru

Кормушин Константин Владимирович,

ассистент кафедры Материаловедения, сварки и производственной безопасности, Казанский национальный исследовательский технический университет им. А. Н. Туполева-КАИ Казань, Россия E-mail: KVKormushin@kai.ru

Исследование выполнено за счет гранта Российского научного фонда (проект № 23-29-00099).

Предмет исследования: в работе представлен численный алгоритм нахождения констант скоростей гетерогенных процессов для последующего их исследования.

Цель исследования: Цифровое прототипирование электроосаждения цинкового покрытия путем выщелачивания обедненной цинковой руды.

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

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

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

ALGORITHM FOR THE NUMERICAL INVESTIGATION OF CATHODE PROCESSES OF ZINC ELECTROEXTRACTION

Nail F. Kashapov

Doctor of Engineering Science, Professor, Professor of the Department of Physics of Advanced Technologies and Materials Science, Kazan Federal University, Deputy Head of «Kazan Scientific Center of the Russian Academy of Sciences» Kazan, Russia

E-mail: kashnail@gmail.com Violetta Yu. Chebakova

Candidate of Science of Physics and Mathematics,

Associate professor at the Department

of Data Analysis and Programming Technologies,

Kazan Federal University

Kazan, Russia

E-mail: vchebakova@mail.ru Evgenia V. Chebakova

Chemical Institute named after A. M. Butlerov, Kazan Federal University Kazan, Russia

E-mail: evgeneachebakova@mail.ru Konstantin V. Kormushin

assistant at the Department of Material Science, Welding and Industrial Safety, Kazan National Research Technical University named after A.N. Tupolev - KAI Kazan, Russia

E-mail: KVKormushin@kai.ru

The study was supported by a grant from the Russian Science Foundation (project No. 23-29-00099).

Subject of research: the paper presents a numerical algorithm for detecting the rate constants of heterogeneous processes for their further investigation.

Purpose of research: Digital prototyping of electrodepo-sition of zinc coating by leaching depleted zinc ore.

Methods of research: the algorithm is presented on the example of predicting the dezincification during parallel hydrogen release reaction in electrolysis of a sodium hydroxide solution, which contains sodium tetrahydroxozincate, and verified using full-scale experiments. So, in this article, there are detected the constants of cathode stepwise electrochemical reactions for zinc powder obtaining using direct search method. The kinetic equation system itself, which describes cathode stepwise reactions, is solved using classical Runge-Kutta method.

Main results of research: the modeling makes allowance for the presence of parallel process of hydrogen release and occurs in the approximation of porous film during the settling-out. The numerical algorithm, developed by the authors, make it also possible to predict the current yield of substance.

Keywords: zinc electroextraction, predicting, numerical optimization, stepwise electrochemical reactions, Run-ge-Kutta method.

ALGORITHM FOR THE NUMERICAL INVESTIGATION OF CATHODE PROCESSES OF ZINC ELECTROEXTRACTION

Nail F. Kashapov, Violetta Yu. Chebakova, Evgenia V. Chebakova, Konstantin V. Kormushin

ВВЕДЕНИЕ

Электрическое цинкование применяется для получения покрытий, обладающих антикоррозионными свойствами, в автомобилестроении и авторемонтном производстве. В работе [1] гальваническое цинкование из сернокислого электролита рассматривается как процесс получения антикоррозионных покрытий. В работе [2] предложено рациональное использование цинковых покрытий для защиты труб от коррозии, предположение выдвинуто на основе анализа различных способов нанесения покрытий. Однако, кроме защитных свойств, данные покрытия обладают достаточной твердостью и износостойкостью, что позволяет использовать электрическое цинкование для восстановления изношенных поверхностей деталей машин. Так, в [3] отмечено следующее преимущество электрического цинкования перед горячим: отсутствие термического воздействия на детали, приводящего к изменению структуры и физико-механических свойств, осаждение покрытия с заданной толщиной, что позволяет минимизировать последующую механическую обработку. При этом покрытие, нанесенное таким способом, по механическим свойствам подобно чистому цинку, т. е. является относительно мягким и более пористым. Также следует отметить, что гальванические не содержат интерметаллических соединений (фаз) и держатся на защищаемом материале только за счет адгезии, т. е. налипания, тогда как получаемые методом горячего цинкования представляют из себя систему железо-цинковых сплавов с постепенно уменьшающимся содержанием железа по мере приближения к внешней стороне покрытия [4]. В работе [5] предлагаются усовершенствования процесса осаждения цинка на восстанавливаемые детали в авторемонтном производстве, позволяющие повысить надежность и прочность сцепления цинка с основой. Также получают гальванические покрытия и из сплавов, содержащих цинк, так, в [6, 7] исследуют свойства цинк-никелевых покрытий, полученных как из щелочных, так и из кислых электролитов. Гальванические способы восстановления имеют несколько разновидностей, из них наиболее распространены: кислые (сульфатные), слабокислые (хлоридные и хлорамонийные) и щелочные (цианидные и цинкатные). От состава и типа электролита зависит качество получаемого осадка, его структура, физико-механические свойства, адгезия с подложкой, виды и скорость протекания электрохимических процессов. Так, в работе [8] обозначено более 20 факторов, влияющих на выбор электролита.

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

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

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

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

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

При выщелачивании руды, в зависимости от концентрации оксида натрия в растворе, цинк в электролите может содержаться как в качестве оксида цинка 7иО, так и в качестве соединения [7и(ОИ)4]2". В работе [10] показано, что 20 % концентрации раствора гидроксида натрия КаОИ хватало, чтобы молярное отношение ионов гидроксильной группы ОН- к содержащемуся цинку было достаточным для существования [7и(ОИ)4]2", при этом, если концентрация ионов гидроксильной группы ОН- мала, то соединение [7и(ОИ)4]2" переходит в оксид цинка 7иО. Таким образом, процесс выщелачивания цинкосодержащей руды в 20 % растворе КаОИ можно описать общей реакцией растворения вида

2 №0Н+ Ъл+ 2 Н20 ^ [2п (ОН)4] + Н2.

Непосредственно по поводу протекания приэлектродных реакций выделения цинка существуют различные мнения. Так, в обзоре [10], опирающемся на работу [11], приведено, что в процессе электроосаждения цинка соединение [7и(ОИ)4]2" подвержено модификации в гидрированный 7и(ОИ)2, при этом

Кашапов Н.Ф., Чебакова В.Ю., Чебакова Е.В., Кормушин К.В.

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

Для построения математической модели, описывающей прикатодные процессы, сформулируем следующие предположения:

1)Na2[Zn(OH)4] в электролите полностью диссоциирует на ионы

N [(0Н)4] Д 2Мв+ + [Хп (ОН)4]2- .

Это справедливо при концентрации ги-дроксида натрия в растворе от 20 %.

2) Ионы натрия №+ на электродах не разряжаются, накапливая с течением времени концентрацию в катодном пространстве и уменьшая ее в анодном.

3) На катоде происходит одновременное выделение водорода 2Н20 + 2е д Н2 + 20Н-

и цинка \_Хп{ОН)А] - + 2е д Хп + 40Н- в следующих стадийных реакциях

[Хп О0Н)Л]2- Д гп2+ + 40Н-, Хп2+ + 2е %Еп ,

Н2О->Н+ +ОН-, Н+ + ОН-4Н2О ,

Л к

Н+ + едН , Н + НдН .

4) Несмотря на то, что данные конкурирующие реакции прямо не связаны между собой, скорости стадийных реакций

:вв„

Л = Аехр I

перенапряжения

ЛТ

где А и В некоторые константы, i = 1,5 экспоненциально зависят от перенапряжения на электродах, которое в нашем случае на катоде будет:

ЛТ

- - 1п(С2п2+Сн+ КРНСП) .

перенапряжения

пР

5)Структура металлического осадка определяется рядом факторов, таких как: состав раствора, наличие в нем примесей, предварительной подготовки поверхности, плотности тока и т. д. В работах [12-14] отмечен губчатый и пористый характер покрытия из осажденного цинка. Поэтому для моделирования осаждения цинка используем приближение пористой пленки. В случае пористой пленки считается, что она не оказывает сопротивления реагентом, подходящим к поверхности раздела сред, и не влияет на скорость реакции. В случае непористой пленки сопротивлением пограничного слоя можно пренебречь, считая основным сопротивление слоя выделившегося продукта. При этом, если процесс контролируется диффузией через пограничный

слой, то скорость процесса

О С

со = — ЬС , 8

где 5 - толщина слоя продукта, S - площадь рабочего электрода. Толщину слоя продукта в приближении решения уравнения Яндера можно считать пропорциональной времени:

58 К 82 2К

— = — и 82 = 2К.

дг г

6) Процесс ассоциации воды

Н+ + ОН- д Н20

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

Таким образом, математическая модель состоит из системы, описывающей электрохимические стадийные прикатодные реакции:

дС,

н2о

дЬ

-к3СН20

дСи

дг

дСн

дь

^зСн2о ~ Сгп^н/(Сгп2+Сн+)Сн+ - к5Сн+ С2пСн/(Сгп2+Ся+) - к6Сн+Сн

ас.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

■н2

дЬ

= к6Сн+Сн+

дС

- Сгп^н/(.Сгп2+Сн+)Сгп(он)2~ ~ (1)

дС2п

дЬ

■ = к2 Сгп2+

дС

гп(ону

дг

^гп(он)2"

■Шх

-(С,

гп(он)

2- - С0) к7

Здесь

2-

О

га(0Н)г- - это диффузия ионов

Zn(ОН)2- через пористую пленку осадка,

С0 - концентрация данных ионов Хп(ОН)2 в объеме электролита. Данная запись кинетических уравнений, описываемых задачами Коши, подразумевает положительность констант скоростей реакций к1 > 0 , i = 1,...,7, коэффициента пропорциональности К > 0 и значений концентраций С] > 0,

} = { И20, Zn, И, Zn2+, СИ+, И2, Zn(OИ)2- }.

Расчет массы осажденного цинкового осадка рассчитывается исходя из концентрации цинка и объема пленки:

Мп = С2пЬ2КТ*б5,38 ,

(2)

где - Т конечное время расчетов. Начальное приближение скоростей

ALGORITHM FOR THE NUMERICAL INVESTIGATION OF CATHODE PROCESSES OF ZINC ELECTROEXTRACTION

Nail F. Kashapov, Violetta Yu. Chebakova, Evgenia V. Chebakova, Konstantin V. Kormushin

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

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

Решение задачи по поиску констант скоростей прикатодных электрохимических процессов и прогнозированию выхода строится на минимизации функционала:

п т 2

\ \ / / - эксперимент расчетные \

ЬЬ\С1,У ) .

1=1У=1

Здесь I отвечает за суммирование по времени, V же отвечает за суммирование по веществам, концентрации которых известны из экспериментов на определенные моменты времени Т1, I = 1,...,п.

В идеале желательно иметь данные как по водороду, так и по массе осаждаемого осадка, но, как правило, это - масса осаждаемого вещества в нескольких временных точках, и функционал приобретает вид:

п 2

М(К,кх,^,к7) = Ь{Мэптримет:ы - мрпГные) (3)

1=1

Численный алгоритм исследования катодных процессов

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

Данный алгоритм опирается на расчет значений функционала (3), от которого ищется минимум. Значение данного функционала вычисляется путем численного решения системы кинетических уравнений (1) методом Рунге-Кутта четвертого порядка и

U (t) - (ch2o, Сн+, Сн, CHi, СZn(OH)1-CZni*, Czn) ,

соотношения (2). Из семейства явных методов Рунге-Кутта четвертого порядка был выбран метод, называемый классическим методом Рунге-Кутта четвертого порядка. Точности четвертого порядка в целом достаточно для соответствия точности снятия экспериментальных данных. Для описания его реализации к данной задаче введем следующие обозначения:

T

2 C с

Zn(OH)

К = (kj, k2,k3, kq.K) ,

F (t K U\- (/l tK'U) 'f2 t'K'U) (tK'U) tK'U)'T ' ' ' ' f5 (t,K,U), f6 (t,K,U),f7 (t,K,U))

Здесь (°)T - операция транспонирования, а функции f¡ (t,K,U),f2 (t,K,U),

f, (t,K_,U) (t,K,U), _ _ f5 (t, k, U) , f6 (t, K, U) , f7 (t, K,U)

соответствуют правым частям системы (1):

frit,К, 0) = -к3СН2О,

f2(t,K,Ü) = к3СН2о - к5 CZnCH/(CZn2+CH+)CH+, f3(t,R,0) = к5Сн+ CZnCH/(CZn2+CH+) - к6Сн+Сн+; fA(t,K,V) = к6Сн+Сн+;

f¡(t,K, I/) = *j2Kt ow>2~ ~ ~ CrnCH/(CZn2+CH+)CZn(OH)2-

f6(t,K,D) = fci CZnCH/(CZn2+CH+)CZn(OH)2- - k2CZn2+ f7(t,K,0) = k2C2„ 2+

Так, система (1) перепишется в векторном виде:

Up-- F (t,K,U), где t0 < t < T.

Здесь T - максимальное значение времени, на котором было зафиксировано экспериментальное значение. Покроем отрезок времени (t0,T) сеткой

an = (t0,tl,t2,-,t = t0 + ¡h,-,tn = T} ,

где h - t¡+1 -t¡ - const. Значения в половинных узлах сетки будем обозначать как:

t 1 - t¡ + h / 2.

¡+—

2

В узлах сетки an введем сеточные функции U¡ -U (t¡) и F¡ (t¡,K,U¡) соответствующие численному решению разностной задачи:

U+i -и1 + hd¡.

Здесь

di -(d?}+2d2)+2df3)+d4)), dW- Fi (ti ,K,Ui ),

d,(2)- F

t j,K,Ui + h*df') /2

i+—

V 2 ,

Кашапов Н.Ф., Чебакова В.Ю., Чебакова Е.В., Кормушин К.В.

dP = FAt i,K,Ui + /2

I A+,

dP= Fa (t1+1,K,U1 + h*dP) .

Сам алгоритм расчета минимума происходит итерационно по каждой из констант (которые далее в описании считаются независимыми координатами).

Первый шаг метода заключается в исследующем поиске вокруг базисной точки g(K,k1,...,ki,...,k1) по вектору приращений H(h,h1,.,h7), где hA - это шаг по координате kA, i = 1,...,7. Для осуществления данного этапа по каждой координате отдельно рассчитываются значения функционала (3) в трех точках: непосредственно в точке g(K,k1,...,ki,...,k7), и точках, отстоящих от нее на шаг по исследуемой координате. Затем точку, в которой функционал принимает минимальное значение, принимаем за новую базисную точку. Данные итерации продолжаются до тех пор, пока базисная точка не перестанет меняться. При этом в алгоритм расчета закладывается, что точка, в которой значение координаты или любой из концентраций при расчете становится меньше нуля, исключается из рассмотрения.

Например, для координаты, отвечающей за коэффициент пропорциональности K в точках

g(K,kv...,ki,...,k7),

gi (K + h,ki,...,ki,...,k7),

g2 (K - h,ki,...,ki,...,k7),

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

После расчета по последней координате получается новая точка g, которая может как совпадать с первоначальной точкой g' = g, так и отличаться от нее g' ф g .

Затем, если первоначальная базисная точка совпадает с вновь получившейся точкой g'= g, то уменьшаем значения компонент вектора покоординатных приращений и продолжаем поиск вокруг базисной точки. В противном случае, когда g ф g, осуществляется поиск по образцу. В нем за новую базисную точку берутся значения g', g + 2*(g- g'), при которых, опять-таки, значение функционала (3) минимально. Условием выхода является

ограничение на размер компонент вектора приращений И.

При расчете целевой функции использовался численный метод из семейства методов Рунге-Кутта четвертого порядка точности, т. е. его погрешность составляет О(дг4), где 1 - шаг по времени. Для поиска целевой функции использовался метод Хука-Дживса, его точность определяется путем задания ограничения возможного уклонения расчетного значения параметров от значений параметров, определяющих минимум е1. Таким образом, погрешность метода представляет шах1 (Д?4.

Результаты численных расчетов

Численный алгоритм, реализующий поставленную математическую постановку осаждения цинкового порошка, верифицирован на натурном эксперименте, описанном в [16]. В ней исследовалось влияние приложенного напряжения на приэлектродные процессы электроосаждения цинка в растворах гидрок-сида натрия на катоде из нержавеющей стали. Концентрация ионов цинка Хп2+ в приготовленном растворе составляет 10 г/л. Рабочий электрод сделан из нержавеющий стали и имеет площадь в 1 см2. Воспользуемся условиями эксперимента для вывода начальных условий:

1000 грамм/л раствора дает 80 %

С

NaOH\i=0

дает 20 %

Таким образом,

, _ 1000 * 0,2 _

С __

NaOH\t=0

250 грамм/л

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0,8 6,25 моль/л.

Молярная масса цинка 65,38 г/моль, и концентрация в 10 г/л дает приблизительно 0,153 моль/л, и, следуя реакции растворения

2ИаОН + гп + №2 [гп(ОН)А] д 2Ма+ + \гп{ОН)^]2- + Н2

будем иметь C

Zn(OH)2

t=0

; 0,153 моль/л.

При этом 0,153 моль/л цинка связывает всего ~ 0,306 моль/л ионов гидроксильной группы, остальные 5,944 моль/л участвуют в реакции диссоциации МаОН д Ыа+ + ОН-, таким образом, СОН_\==д ~ 5,944 моль/л.

Напрямую начальная концентрация ги-дроксильного иона не участвует в выделенных нами прикатодных стадийных реакциях, но ионное произведение воды показывает, что при обычных условиях и СН+\==д к 0,168-14 моль/л. Цинк связан гидроксидной группой, и значение ионной концентрации возьмем близким к «нулю», с учетом того, что осаждение металлов происходит вокруг зародышей, электрохимические процессы стадии образования

ALGORITHM FOR THE NUMERICAL INVESTIGATION OF CATHODE PROCESSES OF ZINC ELECTROEXTRACTION

Nail F. Kashapov, Violetta Yu. Chebakova, Evgenia V. Chebakova, Konstantin V. Kormushin

Таблица 1. Диапазоны получаемой экспериментально массы Zn и расчетных значений при силе тока в 500 А по экспериментальным данным в t = 1 час.

T, min Диапазон экспериментальных данных [16] Расчетные значения

10 0.4371-0.46 0.315

20 1.22-1.38 0.893

40 2.765-2.928 2.58

60 4.514-4.638 4.635

70 5.508-5.693 5.84

100 8.13-8.54 9.1

10 20 30 40 50 60 70 80 90 100

t, min

Рисунок 1. Расхождение результатов расчетов в зависимости от количества выбранных точек.

зародышей (адатомов) отличаются от стадии их роста, которая нас фактически и интересует Сгт\г=0 ~ 0, концентрация же атомов цинка взята равной «нулю» С2п|10 = 0.

Концентрацию воды возьмем

г

С

и,о

t=0

750 грамм :18_

41,7 моль/л,

л моль

при обычных условиях в воде содержится СН2 |=0~ 4,29 * 10-7 Ммоль/л, далее возьмем СД=0 « 0, с учетом необходимого наличия зародышей.

В работе [16] показаны экспериментальные данные по выходу цинка в %, рассчитанных по приведенной в ней следующей формуле: Ш; - т

гп ~

n7„ =

1.22 х I х t

100 %,

где т0 - первоначальная масса катода в г, т1 - масса катода с осажденным на нем цинком в 1-й момент времени, показатели

снимались через каждые 10 минут, I - сила тока в А, а t = 10 * I/60. Так как в статье [16] не приведены данные по выходу водорода, то будем строить функционал только по известной массе осажденного цинка.

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

Кашапов Н.Ф., Чебакова В.Ю., Чебакова Е.В., Кормушин К.В.

Также было проведено исследование влияния количества взятых точек и показано, что оно не оказывает решающего влияния на погрешность. В работе [15] этот вопрос исследовался для процессов выделения вещества в газовой фазе, когда применялось приближении неизменности и однородности площади электрода. Процесс отложения цинкового осадка на электрод показан на рисунке (1). На нем приведено сравнение рассчитываемых данных при силе тока 500А в двух случаях: расчеты приводились по одной точке 1 = 1 час, и по двум точкам (1 = 1 час и 1 = 10 мин).

ЗАКЛЮЧЕНИЕ И ВЫВОДЫ

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

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

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

1. Манаков, Д. А. Актуальные вопросы защиты от коррозии / Д. А. Манаков, О. В. Немыкина. - Текст: непосредственный // Молодежный вестник ИрГТУ. - 2019. - Т. 9, № 2. - С. 88-91.

2. Проскуркин, Е. В. Цинковые покрытия - основные современные системы защиты труб от коррозии /

E. В. Проскуркин, В. А. Геловани, А. Н. Сонк, И. В. Петров, И. П. Ярема, Д. А Сухомлин. - Текст: непосредственный // Сталь. - 2018. - № 6. - С. 32-37.

3. Юдин, В. М. Исследование электролитического цинкования при высоких плотностях тока / В. М. Юдин, М. Н. Вихарев. - Текст: непосредственный // Труды ГОСНИТИ. - 2008. - Т. 102. - С. 162-164.

4. Полькин, В. И. Цинк для защиты от коррозии / В. И. Полькин. - Текст: непосредственный // Гидротехника. - 2019. - № 1 (54). - С. 86-89.

5. Захаров, Ю. А. Совершенствование технологического процесса гальванического цинкования деталей транспортно-технологических машин и комплексов / Ю. А. Захаров, И. А. Спицын, Е. В. Ремзин, Г. А. Мусатов. - Текст: непосредственный // Модели, системы, сети в экономике, технике, природе и обществе. - 2014.

- № 4 (12). - C. 105-111.

6. Бурляев, Д. В. Цинк-никелевые сплавные покрытия: кинетика электроосаждения, коррозия и селективное растворение. Обзор / Д. В. Бурляев, О. А. Козадеров, П. Волович. - Текст: непосредственный // Конденсированные среды и межфазные границы. - 2021. - Т. 23(1). - С. 3-15.

7. Dutra, C. A. M. Corrosion Resistance of Zn and Zn-Ni Electrodeposits: Morphological Characterization and Phases Identification/ Conceiçao A. M. Dutra , José W. J. Silva , Roberto Z. Nakazato // Materials Sciences and Applications, - 2013. - vol. 4. - P. 644-648

8. Захаров, Ю. А. Выбор, контроль и корректировка электролита цинкования восстанавливаемых поверхностей деталей автомобилей / Ю. А. Захаров, Г. А. Мусатов.

- Текст: электронный // Инженерный вестник Дона.

- 2015. - № 2. - URL: ivdon.ru/ru/magazine/archive/ n2y2015/2861 (дата обращения: 08.08.2021).

9. Петров, В. В. Технология переработки отходов твердосплавного инструмента / В. В. Петров, А. Е. Проценко, К. Р. Сапожник. - Текст: непосредственный // ^рои-тельные материалы и изделия. - 2022 - Т. 5, № 3. - С. 27-34.

10. Zhang, Y. Leaching of zinc from calcined smithsonite using sodium hydroxide / Yucheng Zhang, Jinxia Deng, Jun Chen, Pando Yu, Xianran Xing // Hydromettallurgy. -2013. - vol. 131. - P. 89-92.

11. Brenner, A. Electrodeposition of alloys principles and practice / Brenner A. - New York: Academic Press, 1963, -vol. 1. - 402 p.

12. Abdel-Aal, E. A. Kinetics of sulfuric acid leaching of low-grade zinc silicate ore / E. A. Abdel-Aal // Hydrometallurgy. -2000. - Vol. 55. - P. 247-254.

13. Leung, P. K. Zinc deposition and dissolution in methanesulfonic acid onto a carbon composite electrode as the negative electrode reactions in a hybrid redox flow battery / P. K. Leung, C. Ponce-De-Leon, C. T. J. Low,

F. C. Walsh, // Electrochimica Acta . - 2011. - vol. 56(18). -P. 6536-6546.

14. Gurmen, S. laboratory-scale investigation of alkaline zinc electrowinning/ S. Gurmen, M. A. Emre//Minerals

ALGORITHM FOR THE NUMERICAL INVESTIGATION OF CATHODE PROCESSES OF ZINC ELECTROEXTRACTION

Nail F. Kashapov, Violetta Yu. Chebakova, Evgenia V. Chebakova, Konstantin V. Kormushin

Engineering. - 2003 - Vol. 16(6). - P. 559-562.

15. Кашапов, Р. Н. Кинетика двухфазных газожидкостных сред в процессах электролиза / Р. Н. Кашапов, Л. Н. Кашапов, Н. Ф. Кашапов, В. Ю. Чебакова. - Текст: непосредственный // Теплофизика высоких температур.

- 2021. - Т. 59, № 6. - С. 869-876.

16. Zhang, Y. The Electrowinning of Zinc from Sodium Hydroxide Solutions/ Yucheng Zhang, Jinxia Deng, Jun Chen, Pando Yu, Xianran Xing. //Hydromettallurgy. - 2014.

- vol. 146. - P. 59-63.

i Надоели баннеры? Вы всегда можете отключить рекламу.