Научный журнал
Современные наукоемкие технологии
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

АЛГОРИТМ РАСЧЕТА ВОЛЬТ-АМПЕРНОЙ ХАРАКТЕРИСТИКИ В ДИФФУЗИОННОМ СЛОЕ ДЛЯ МЕМБРАННЫХ СИСТЕМ В ГАЛЬВАНОДИНАМИЧЕСКОМ РЕЖИМЕ

Чубырь Н.О. 1 Уртенов М.Х. 2 Коваленко А.В. 2 Узденова А.М. 3
1 ФГБОУ ВО «Кубанский государственный технологический университет»
2 ФГБОУ ВО «Кубанский государственный университет»
3 ФГБОУ ВО «Карачаево-Черкесский государственный университет имени У.Д. Алиева»
Электромембранные системы функционируют либо в потенциодинамическом режиме, когда задается скачок потенциала в системе, или в гальванодинамическом режиме, когда задается плотность тока. Для моделирования переноса в потенциодинамическом режиме используется система уравнений Нернста – Планка совместно с уравнением Пуассона для вычисления потенциала. В предыдущих работах были предложены формула и алгоритм вычисления вольт-амперной характеристики (ВАХ) в потенциодинамическом режиме, устойчивый относительно ошибок округления, рассчитаны ВАХ и проведено сравнение с некоторыми экспериментальными ВАХ, снятыми в потенциодинамическом режиме. Однако большинство экспериментальных ВАХ измерено для гальванодинамического режима, для которого в настоящее время нет математической модели переноса и, соответственно, алгоритма расчета ВАХ. Таким образом, появляется проблема вывода уравнений и краевых условий, для моделирования переноса произвольной соли в диффузионном слое в гальванодинамическом режиме и алгоритма расчета ВАХ. Решению этих проблем посвящена данная работа. В работе выведено новое уравнение для плотности тока, которое совместно с уравнениями Нернста – Планка и известной формулой для напряженности электрического поля позволяет построить математическую модель и алгоритм расчета вольт-амперной характеристики в диффузионном слое возле ионообменной мембраны для одномерного нестационарного переноса произвольной соли в гальванодинамическом режиме. Показано, что для бинарной соли полученные уравнения значительно упрощаются.
вольт-амперная характеристика
диффузионный слой
мембранная система
ионообменная мембрана
электродиализ
математическая модель
уравнение Нернста – Планка – Пуассона
1. Ньюмен Дж. Электрохимические системы. М.: Мир, 1977. 463 с.
2. Сарапулова В.В., Небавская К.А., Невакшенова Е.Е., Козмай А.Э., Письменская Н.Д., Ларше К., Систа Ф. Эволюция электрохимических характеристик мембраны АМХ-SB после контакта её поверхности с красным вином // Сорбционные и зроматографические процессы. 2016. Т. 16. № 5. С. 672–685.
3. Коваленко А.В., Бостанов Р.А., Лайпанова З.М., Уртенов М.Х. Моделирование гальванодинамического режима для одномерного нестационарного переноса бинарного электролита // Фундаментальные исследования. 2015. № 12–2. С. 273–277.
4. Харкац Ю.И. О механизме возникновения запредельных токов на границе ионообменная мембрана/электролит // Электрохимия. 1985. Т. 21. № 7. С. 974–977.
5. Gil V.V., Andreeva M.A., Jansezian L., Han J., Pismenskaya N.D., Nikonenko V.V., Larchet C., Dammak L. Impact of heterogeneous cation-exchange membrane surface modification on chronopotentiometric and current-voltage characteristics in NaCl, CaCl2 and MgCl2 solutions. Electrochimica Acta. 2018. V. 281. P. 472–485. DOI: 10.1016/j.electacta.2018.05.195.
6. Чубырь Н.О., Коваленко А.В., Уртенов М.А.Х. Численные и асимптотические методы анализа переноса 1:1 электролита в мембранных системах. Краснодар, 2018. 106 с.

Электромембранные технологии в настоящее время считаются критически важными технологиями. Они используются для решения таких актуальных проблем, как создание безотходных технологий, включая очистку сточных вод, например, гальвонического производства, содержащих ионы тяжелых металлов (цинк, свинец, кадмий и др.) [1]. Важнейшей проблемой являяется очиска воды от радиоактивных изотопов. Электромембранные технологии могут быть использованы для дезактивации воды, если радиоактивные вещества находятся в ней в растворенном состоянии в виде ионов, например для дезактивации слабоминерализованных сбросных вод ядерных реакторов, воды контурного охлождения. При этом требуется ее глубокое обессоливание. Одной из наиболее важных интегральных характеристик переноса ионов соли в мембранных системах является вольт-амперная характеристика (ВАХ), которая применяется при разработке оптимальной конструкции и выборе технологического режима работы, оценке соответствия теоретических и экспериментальных ВАХ [2, 3]. Именно используя ВАХ, вводятся фундаментальные понятия таких критических плотностей токов, как предельная, экзальтационная и Харкаца, и так далее [1]. В работе выведено новое уравнение, заменяющее уравнение Пуассона, позволяющее построить математическую модель переноса и алгоритм расчета вольт-амперной характеристики для произвольной соли в гальванодинамическом режиме. Статья является продолжением и обобщением результатов работы [4] на общий случай произвольной соли.

Рассмотрим отдающий диффузионый слой (ДС). Пусть H – ширина диффузионного слоя, x = 0 соответствует глубине раствора, где выполняется условие локальной электронейтральности, а x = H – границе ионообменная мембрана/раствор.

Перенос ионов произвольной соли в ДС в потенциодинамическом режиме, описывается уравнениями [2, 4]:

chub01.wmf (1)

chub02.wmf (2)

chub03.wmf (3)

chub04.wmf (4)

где chub05.wmf

В связи с тем, что уравнение Пуассона (3) позволяет вычислить значение потенциала, то (1)–(3) используется для моделирования переноса в потенциодинамическом режиме [2]. Нашей целью является вывод дифференциального уравнения для плотности тока, которая должна заменить уравнение (3), формулы и алгоритма расчета ВАХ.

Материалы и методы исследования

Путем тождественных преобразований из (1)–(4) выводятся:

1) формула для вычисления E, которое заменяет формулу (4) для вычисления I;

2) новое дифференциальное уравнение для I.

В итоге получаем математическую модель переноса ионов произвольной соли в гальванодинамическом режиме.

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

Результаты исследования и их обсуждение

1. Преобразование (1)–(4)

1) Вывод формулы для вычисления I.

Из (1) и (4) имеем

chub06.wmf

следовательно (см. [2]):

chub07.wmf (5)

Обозначим

chub08.wmf, RОм(C) = 1/χ(C),

тогда получим

chub09.wmf (6)

Заметим, что RОм(C) проводимость раствора [2].

2) Вывод нового дифференциального уравнения для I.

Подставим (6) в уравнение Пуассона (3), тогда получим уравнение для плотности тока I:

chub10.wmf

или

chub11.wmf

или

chub12.wmf (7)

3) Преобразование уравнений (1), (2).

Дифференцируем (1) по x и получим

chub14.wmf

chub15.wmf

chub16.wmf

chub17.wmf

Заменим chub18.wmf исходя из (3):

chub19.wmf

или

chub20.wmf

Откуда

chub21.wmf (8)

или

chub22.wmf (9)

Уравнение (9) – это нелинейное параболического типа дифференциальное уравнение для парциальных концентраций Ci, зависящее еще и от плотности тока I.

В итоге получили систему уравнений (7), (8), (9) для Ci, ji, i = 1,..., n, I. Она позволяет вместе с соответствующими краевыми условиями однозначно определить Ci, ji, i = 1,..., n, I и являются математической моделью переноса произвольной соли в гальванодинамическом режиме.

4) Алгоритм расчёта ВАХ.

а) Задается некоторый темп роста α функции I(t, H), начиная с некоторого значения I0: chub23.wmf

б) Решается (7), (8) с соответствующими граничными условиями и I(t, x) и Ci, i = 1,..., n.

в) Рассчитывается E по формуле (6).

г) Находится скачок потенциала по формуле chub24.wmf

2. Упрощение модели в стационарном случае

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

chub25.wmf

chub26.wmf

chub27.wmf

поэтому алгоритм рсчета ВАХ становится проще.

Алгоритм расчета ВАХ

а) Задается набор значений Il, l = 1,..., m.

б) Решается уравнение (6) с соответствующими краевыми условиями и I = Il, находятся Ci, i = 1,..., n.

в) Определяем E:

chub28.wmf

г) Определяем скачок потенциала по формуле chub29.wmf

Набор chub30.wmf представляет собой вольт-амперную характеристику.

3. Упрощение модели при выполнении условия электронейтральности chub31.wmf

Прт выполнения условия chub32.wmf уравнения (1), (2), (6), (7) упрощается, так как I не зависит от x, хотя потоки могут зависит от x. Действительно, из (2) следует

chub33.wmf

или

chub34.wmf, т.е. I = I(t)

где I(t) задается.

Кроме того, соотношение chub35.wmf позволяет найти одну из концентраций если известны chub36.wmf концентрации, например если известны Ci, i = 1,..., n – 1, то chub37.wmf. Следовательно, число уравнений для определения концентраций должно быть n – 1, поэтому из (6) получаем

chub38.wmf (10)

где E определяется по формуле (6).

Уравнение (7) выполняется тождественно.

Алгоритм расчета ВАХ

а) Задается некоторый темп роста α плотности тока I(t), начиная с некоторого значения I0: chub39.wmf.

б) Решается система уравнений (10) и определяется Ci, i = 1,..., n – 1, затем вычисляется Cn по формуле chub40.wmf.

в) Рассчитывается напряженность электрического поля по формуле (6).

г) Находится скачок потенциала по формуле chub41.wmf

4. Упрощение модели п. 3 для бинарного электролита

Для бинарного электролита, можно ввести эквивалентную концентрацию C: chub42.wmf, тогда из (10) получаем

chub43.wmf

или с использованием эквивалентной концентрации получим

chub44.wmf,

chub45.wmf

Откуда следует

chub46.wmf

или

chub47.wmf (11)

где chub48.wmf – коэффициент диффузии электролита [1, 2].

Вычислим

chub58a.wmf

chub58b.wmf,

chub50.wmf

или

chub51.wmf (12)

Вычислим j1 и j2 потоки. Несложно показать

chub52.wmf (13)

chub53.wmf (14)

Алгоритм расчета ВАХ

а) Задается некоторый темп роста α плотности тока I(t), начиная с некоторого значения I0: chub54.wmf.

б) Из (11) находим С, и затем chub55.wmf, chub56.wmf.

в) Рассчитываются потоки по формулам (13), (14).

г Рассчитывается E по формуле (12).

д) Находится скачок потенциала по формуле chub57.wmf

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

Заключение

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

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта 19-08-00252 А «Теоретическое и экспериментальное исследование вольт-амперных характеристик электромембранных систем».


Библиографическая ссылка

Чубырь Н.О., Уртенов М.Х., Коваленко А.В., Узденова А.М. АЛГОРИТМ РАСЧЕТА ВОЛЬТ-АМПЕРНОЙ ХАРАКТЕРИСТИКИ В ДИФФУЗИОННОМ СЛОЕ ДЛЯ МЕМБРАННЫХ СИСТЕМ В ГАЛЬВАНОДИНАМИЧЕСКОМ РЕЖИМЕ // Современные наукоемкие технологии. – 2019. – № 10-1. – С. 92-96;
URL: https://top-technologies.ru/ru/article/view?id=37703 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674