Влияние реакции диссоциации/рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 3. Оценка возможности возникновения гравитационной конвекции Текст научной статьи по специальности «Химические технологии»

УДК 303.732.4+514.84+515.1+530.1

01.00.00 Физико-математические науки

ВЛИЯНИЕ РЕАКЦИИ ДИССОЦИАЦИИ/РЕКОМБИНАЦИИ МОЛЕКУЛ ВОДЫ НА ПЕРЕНОС 1:1 ЭЛЕКТРОЛИТА В МЕМБРАННЫХ СИСТЕМАХ В ДИФФУЗИОННОМ СЛОЕ. ЧАСТЬ 3. ОЦЕНКА ВОЗМОЖНОСТИ ВОЗНИКНОВЕНИЯ ГРАВИТАЦИОННОЙ КОНВЕКЦИИ

Коваленко Анна Владимировна к.э.н., доцент

Scopus Author ID: 55328224000 SPIN-код автора: 3693-4813 [email protected]

Уртенов Махамет Али Хусеевич д.ф.-м.н., профессор Scopus Author ID: 6603363090 SPIN-код: 7189-0748 [email protected]

Письменский Александр Владимирович к.ф.-м.н., доцент Scopus Author ID:13004856800 SPIN-код: 9932-7747 [email protected]

Кубанский государственный университет, Россия,350040, Краснодар, Ставропольская, 149,

Данная статья является продолжением предыдущих работ авторов [Влияние реакции диссоциации/ рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 1. Математическая модель // Научный журнал КубГАУ, 2016. -№07(121) и Влияние реакции диссоциации/ рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 2. Асимптотический анализ // Научный журнал КубГАУ, 2016. -№08(122)] и посвящена оценке возможности возникновения гравитационной конвекции за счет рекомбинации ионов водорода и гидроксила. В статье представлено решение краевой задачи, являющейся математической моделью электродиффузии четырех сортов ионов одновременно (двух ионов соли, а также ионов водорода и гидроксила) в диффузионном слое в электромембранных системах с идеально селективной мембраной, с уравнением теплопроводности и уравнением Навье-Стокса. В статье показана возможность возникновения гравитационной конвекции за счет экзотермической реакции рекомбинации молекул воды в глубине раствора. В статье рассматриваться реакция рекомбинации ионов водорода и

UDC 303.732.4+514.84+515.1+530.1 Physic and mathematics

THE INFLUENCE OF THE REACTION OF DISSOCIATION / RECOMBINATION OF MOLECULES OF WATER ON TRANSPORTATION OF ELECTROLYTE 1: 1 IN MEMBRANE SYSTEMS IN THE DIFFUSION LAYER. PART 3. EVALUATION OF THE POSSIBILITY OF THE EMERGENCE OF GRAVITATIONAL CONVECTION

Kovalenko Anna Vladimirovna Cand.Econ.Sci., associate professor Scopus Author ID: 55328224000 SPIN-code: 3693-4813 [email protected]

Urtenov Makhamet Ali Khuseevich Dr.Sci.Phys.-Math., professor Scopus Author ID: 6603363090 SPIN-code: 7189-0748 [email protected]

Pismenskiy Alexander Vladimirovich Cand.Phys.-Math.Sci., associate professor Scopus Author ID: 13004856800 SPIN-code: 9932-7747 [email protected]

Kuban State University, Krasnodar, Russia

This article is a continuation of the previous works of the authors [The influence of reaction dissociation / recombination of molecules of water on transportation of electrolyte 1:1 in the membrane systems in the diffusion layer. Part 1. Mathematical model // Scientific journal of Kuban State Agrarian University, 2016. No. 07(121) and The influence of the reaction of dissociation / recombination of molecules of water on transportation of electrolyte 1: 1 in membrane systems in the diffusion layer. Part 2. Asymptotic analysis // Scientific journal of Kuban State Agrarian University, 2016. - №08(122)] and devoted to assessing the possibility of gravitational convection due to the recombination of hydrogen and hydroxyl ions. The article presents the solution of a boundary-value problem, which is a mathematical model of electrodiffusion for the four types of ions at the same time (two ions of salts and hydrogen and hydroxyl ions) in the diffusion layer in electro-membrane systems with ideal selective membrane, with the heat transfer equation and the Navier-Stokes equation. The article shows the possibility of the emergence of gravitational convection due to the exothermic reaction of recombination of water molecules in the depth of the solution. The article considered the reaction of recombination of hydrogen ions and hydroxyl, although the main results can be applied, after

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

Ключевые слова: ГРАВИТАЦИОННАЯ КОНВЕКЦИЯ, ЭЛЕКТРОКОНВЕКЦИЯ, ДИГИДРОФОСФАТ, РЕАКЦИЯ ДИССОЦИАЦИИ, РЕАКЦИЯ РЕКОМБИНАЦИИ МОЛЕКУЛ ВОДЫ, ЭЛЕКТРОМЕМБРАННЫЕ СИСТЕМЫ, ДИФФУЗИОННЫЙ СЛОЙ, РЕАКЦИОННЫЙ СЛОЙ, ИДЕАЛЬНО СЕЛЕКТИВНАЯ МЕМБРАНА

Введение

Реакции диссоциации/рекомбинации молекул воды в растворе ИаС1 происходят лишь в сверхпредельных токовых режимах. Однако существуют растворы, в которых подобные реакции происходят и при допредельных токах. Примерами являются многие амфолит-содержащие растворы: виноматериалы, соки, молочные продукты, продукты микробиологической переработки биомассы (аминокислоты, анионы многоосновных карбоновых кислот), муниципальные сточные воды (анионы фосфорной кислоты) и др. Молекулы амфолитов содержат одновременно и кислотные, и основные группы, и поэтому в водных растворах диссоциируют и как кислоты с отщеплением водородных ионов Н+, и как основания с отщеплением гидроксильных ионов ОН-. Эти реакции также сопровождаются тепловыми эффектами.

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

appropriate modifications, and to amfolit-containing solutions, such as wine, juices, dairy products, microbiological processing of biomass (amino acids, anions of polybasic carboxylic acids), municipal effluent (anions of phosphoric acid), etc.

Keywords: DIHYDROGEN PHOSPHATE, ELECTROCONVECTION, DISSOCIATION REACTION, RECOMBINATION REACTION OF MOLECULES OF WATER, ELECTRO-MEMBRANE SYSTEMS,

(ГК). В то же время при электродиализе в тех же условиях растворов, не содержащих амфолитов, например растворе ИаС1, этот эффект не наблюдается, или наблюдается в значительно меньшей степени.

Это позволяет выдвинуть гипотезу, что основной причиной наблюдаемой ГК являются тепловые эффекты, возникающие наряду с другими эффектами концентрационной поляризации. Это поглощение теплоты на границе раствор/анионообменная мембрана (АО) в результате диссоциации (депротонирования) аниона-амфолита:

Н 2 РО- + Q1 ® Н ++ НРО\-, и поглощение теплоты на некотором удалении от межфазной границы в результате реакции рекомбинации (протонирования) с образованием молекулярной формы амфолита - фосфорной кислоты:

Н + + Н2 РО4- + Q2 ® Н3 РО4.

На то, что две эти реакции разнесены в пространстве, указывают вычисления, сделанные в рамках Ш трехслойной модели [7]. Согласно этой модели протоны появляются в обедненном диффузионном слое (ДС) вследствие их доннановского исключения из мембраны и рекомбинируют с однозарядными ионами (гидрокарбонатами) на некотором расстоянии от мембраны внутри ДС.

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

Поскольку рекомбинация ионов Н + и ОН- происходит в глубине раствора и при этом выделяется тепло в узком реакционном слое (РС), то это может привести к ощутимому неравномерному разогреву раствора, и как, следствие к ГК. Для проверки этой гипотезы рассчитаем тепло выделяемое при рекомбинации ионов Н + и ОН-, а затем, используя уравнения теплопроводности и Навье-Стокса рассчитаем течение раствора.

1. Оценка тепла выделяющегося при рекомбинации

Количество тепла выделяющегося при образовании 1 моль воды равно 56,6 кДж/моль [4]. Следовательно, мощность Q источников тепла выделяющегося в каждой точке x пропорциональна количеству молей образующихся молекул воды KrC3C4, за вычетом количества дислоцирующих молекул KdCHlO, с коэффициентом пропорциональности 56,6кДж/моль, т.е.: Q = 56,6(KrC3C4 -KdCHlO).

Из формул (6) и (7) (см. [1]) получаем: KrC3C4 -KdCH2О = - — у3.

2 —х

3 —

Следовательно Q = -56,6 • 103 —у3.

—х

Поскольку вне РС ]3 постоянна, то Q = 0 вне РС. Количество тепла Qw, которое выделяется в РС с единичным сечением за единицу времени:

И Х2 ,х2 — 1

Qw = | Qdx = | Qdx =- 56,6 • 103 | —]з—х = - 56,6 • 103 Г7з ГХ2 ) - ]з (хХ)),

0 х\ х\

где х1 - левая граница, а х2 - правая граница реакционного слоя, И -ширина диффузионного слоя (в статьях [1] и [2] он был обозначен 3).

3 1

Используя формулу для ]ъ из [2] получаем Qw = 56,6 10— ^.

^

2. Оценка повышения температуры в реакционном слое

Выделение значительного количества тепла в РС может привести к ощутимому локальному повышению температуры в РС, что в свою очередь может вызвать гравитационную конвекцию.

2.1. Предварительная оценка повышения температуры

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

Количество тепла, которое необходимо сообщить однородному изолированному телу, чтобы повысить его температуру на АТ за единицу

времени равно Q = стЛТ = cpVЛT, где c - удельная теплоемкость, m - масса тела, р - плотность тела, V - его объем. Если бы РС имел единичную глубину и длину, то его температура повысилась бы за счет количества

тепла Qw, выделяемого в таком РС на ЛТ, где ЛT = = , где d

cpV cpd

ширина РС равная [2] d »у^Ь . Поэтому, например, при d = 10-4 см,

с = 4,2^^, имеем ЛТ »1,3К . Температуру можно выразить через ^:

г

ЛТ » 0.36^ .

Таким образом, повышение температуры, как и следовало ожидать, пропорционально плотности тока ^. Так, например, при ^ = 0.4 А, получаем ЛТ »1.4К. Известно [5], что локальное повышение температуры на 1К ведет к интенсивной гравитационной конвекции.

2.2 Определение распределения температуры с помощью нестационарного уравнения теплопроводности с сосредоточенным мгновенным источником тепла

Очевидно, что полученная выше оценка повышения температуры РС является грубой «прикидкой», поскольку из-за теплообмена РС с окружающей частью ДС, а также последнего с окружающим раствором распределение температуры будет неравномерным. Для оценки распределения температуры с помощью нестационарного уравнения теплопроводности с сосредоточенным мгновенным источником тепла, как и выше, предположим, что все тепло выделяется в РС. Так как, процесс рекомбинации происходит во много раз быстрее возникновения и развития ГК, можно предположить, что выделение тепла происходит сразу в начальный момент времени ? = 0.

Поскольку, РС значительно меньше ДС, то предполагаем, что он стягивается в точку х1 = вИ [2], а ДС бесконечно большим, причем вне диффузионного слоя поддерживается постоянная температура Т0.

Для нахождения распределения температуры Т (¿, х) в ДС воспользуемся одномерным нестационарным уравнением теплопроводности с внутренним точечным источником тепла, который выделяется в начальный момент времени (1) с начальным условием (2)

(1)

ЭТ Э ,, ЭТ , „ ,

ср— = — (к—) + г,х - хх), Эt Эх Эх

Т( 0, х) = То.

(2)

(х - х\ ) 2

2

Решение задачи (1), (2) имеет вид [6]: T(t,x) = Т0 +—% е 4а 2t

При х = х1, получаем Т(t,x1) = Т0 +

а

2срл/яа2^

, при t ® 0 +, что

объясняется тем, что все тепло выделяется в этой точке. Если х ф хь то T(t,x) ® Т0, при t ® 0 + . Таким образом, изменение температуры происходит только в некоторой окрестности точки х = х1 (рис. 1).

Из рис. 1 следует, что гравитационная конвекция, начинаясь вблизи выхода из КО, достаточно быстро распространяется на весь КО.

а)

б)

в)

г)

д) е)

Рисунок 1. Распределение температуры в поперечном сечении канала и линий тока раствора при средней скорости протока 10'6 м/с при разных значениях времени: а) и б) при г = 10-4 с, в) и г) г = 10-3 с, д) и е) г = 10-2 с.

2.3 Определение распределения температуры с помощью нестационарного уравнения теплопроводности с сосредоточенным постоянно действующим источником тепла

Поскольку реакция рекомбинации молекул воды происходит постоянно, то, в РС тепло выделяется постоянно и ГК, очевидно, будет более интенсивной.

В случае сосредоточенного постоянно действующего источника тепла в точке х1, краевая задача для температуры имеет вид:

дТ д ,, дТ „

= — (к—) + х - х!) дг дх дх

с начальным условием Т( 0,х) = Т0

Решение этой задачи имеет вид [6]:

T(t,x) = Т0 +

кср

С - (х - ху )2 — е л "2 а

Л

4а 2Г , х х1|

С С1 х - х1 ^ ^

2а2

ег/

V V

2аТг

-1

У У

(5)

у

Как видно на рис. 2, распределение температуры, как и в случае с мгновенным источником тепла, приводит к быстрому развитию ГК. Уже для времени t = 0.1с (рис. 2г) хорошо заметно отклонение течения от Пуазейлевского и его сосредоточение в области повышенной температуры. При t = 1с и в последующие моменты времени наблюдается развитое вихревое течение (рис. 2е). Причем скорость этого течения возрастает со временем. Например, профили скорости на рис. 3 показывают, что скорость вихря при t = 10с почти в 2 раза выше, чем при t = 1с.

2.4 Определение распределения температуры с помощью стационарного уравнения теплопроводности с сосредоточенным источником тепла

Если предположить, что в глубине раствора и на межфазной границе поддерживается стационарное состояние, то можно ожидать, что распределение температуры со временем стационарно, тогда вместо уравнения (3) можно использовать стационарное уравнение:

(6)

г) г)Т

-Э- (к Э-) + /о,х) = 0. Эх Эх

Предположим, как и выше, что:

1) Все тепло выделяется в РС.

2) Так как, РС значительно меньше ДС, то предполагаем, что он стягивается в точку х1, а ДС значительно больше, причем вне диффузионного слоя поддерживается постоянная температура.

Тогда для определения температуры имеем уравнение (6), где ¡(г,х) = <Ъ(х - X!), 2 = 56,6 103

кТ" + <Ъ(х - х1) = 0 с краевыми условиями Т( 0) = Т0, Т(И) = Т0

3 Ез т д дТ . , .

3 3 ^, т.е. —(к—) + <8(х - х1) = 0 или дх дх

а)

б)

в)

г)

(7)

(8)

д)

е)

Рисунок 2. Распределение температуры в поперечном сечении КО и линий тока раствора при средней скорости протока 10-6 м/с при разных значениях времени: а) и б) при t = 10-2 с, в) и г) при t = 10-1 с, д) и е) t = 1 с.

Замечание 1. Можно было бы предположить, что ДС бесконечно большой, и ставить краевые условия Т(-<*) = т0, Т(+¥) = т0. Однако уравнение (7), как несложно убедиться, с такими условиями решения не имеет.

Общее решение (7) имеет вид (9), где с(х) - функция Хэвисайда.

в Graph: Velocity field, у component Cm/s)

■■■ Г = 10с .

Г = 0.1с

-----¿_ _________

рррШ 1 т = ю~2 с" t = 10" с

Уф =

Рисунок 3. Динамика распределения продольной (у) составляющей скорости в поперечном сечении канала (посередине)

T(x) = -Q(xc(x - xY) + C\x) + C2: k

(9)

Найдем теперь С1 и С2 используя условия (6), тогда:

б.

T( x) = ( xC( x - xi) - x) + To k

(10)

2.4. Определение распределения температуры с помощью уравнения теплопроводности с распределенным источником тепла

Полученные выше в п.2.2 оценки можно уточнить с учетом того, что в нашем случае внутренние источники тепла непостоянны:

/(I, х) = /(I, х) = 0 = -56, 6 • 103 —/3 ,

ах

тогда вместо (7) можно использовать уравнение:

^-(А = 56,6 • 103 — /3. (11)

дх дх ах

Интегрируя один раз (11) получим (12), где а постоянная интегрирования.

дТ п

к— = 56,6 103 • /3 + а , (12)

дх

Поскольку ДС занимает небольшую часть КО, то можно предположить, что в ядре потока температура постоянна, т.е. при х = 0 поддерживается постоянная температура Т (0) = Т0. Интегрируя уравнение (12) еще раз от 0 до произвольного х, получим:

3 х

Т(х) = Т0 + ах + • | Jз(s)ds . (13)

к к 0

Из [2] следует, что /3 = — %, при х е [0, х:], где х: = (0-л1Ш\ 1п Ш)И и

р

/3 = 0, при [х2,И], где х2 = (0 + д/Ш| 1пШ)И . Поэтому при хе [0,х1 ], получаем

3

ГГ^У 1 гг, 1 У 56,6 • 10 т . Г 1 1

Т(х) = Т0 +-(а + —--¡ж)х, а при [х2,И]:

к р

а 56 6 103 56 6 103 х2 56 6 103 х

Т(х) = Т0 + Vх + —Ч--I /3(Э— + —--I /3(Э— + —6--I /3(Э—

к к Г1 к к

0 х1 х2

Так как /3 = 0, при [х2 ,И], последний интеграл равен нулю, интеграл

х2 х1 1 а 566 •Ю3 1

|/3(8— мал, а | /3(8— = — %хь то Т(х) » Т0 +-х + —--— 1жх1.

Ч 0 р к к р

3) Поскольку отрезок [х1 ,х2 ] мал, то предыдущие формулы можно

продолжить на него по непрерывности.

Остается определить постоянную интегрирования а исходя из

краевого условия при х = И. Постановка краевого условия при х = И

затруднено тем, что х = И является межфазной границей, в которой

происходят сложные процессы. Наиболее существенными из них являются

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

1) Предположим процессы диссоциации молекул воды и джоулев нагрев раствора уравновешивают друг - друга и температура при х = И остается неизменной, поэтому Т(И) = Т0. Тогда определяя а из условия

3 3

Т(И) = т0 + аи +56'6 10 •—% ей = Т0, получаем а = -56'6 10 • % е,

к к F F

56 6 10З

следовательно, функция Т(х) имеет вид Т(х) = Т0 + —'-(1 -е)1мх, при

кБ

3

хе [0,еИ] . Т(х) = Т0 + —6-е^(к -х), при хе [еи,и].

кБ

Таким образом, функция Т(х) на отрезке [0,еИ] линейно растет от

56 6 10 3

Т0 до Ттах = Т0 + —'-(1 -е)е1м, а на отрезке [еи,и] линейно падает до Т0.

кБ

2) Джоулев нагрев раствора преобладает, по крайней мере, вначале, над охлаждающим эффектом диссоциации молекул воды. Этот случай рассматривается ниже в п.З.

3. Учет джоулевого нагрева раствора

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

^ та Л мА -,-г, Л-3 Вт л 2 Дж раствора Qдж = I Аф = 10—-• 1В = 10 10 3—- = 10 —.

см см см • с

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

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

смоделировано также как выше в виде точечного источника в некоторой

точке х2 > х1.

а)

б)

в) г)

Рисунок 4. Распределение температуры и линии тока раствора при разных временах: а) и б) при t = 10-4 с; в) и г) при ? = 10-3 с.

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

получим

ТО,х) = Т0 +

Q

2cp^[шa2t

(х - Х1 )2 4а ^ +. QД

( х - х2 )

2

2срл[ка2

ж е 4а2t

график

которой и соответствующие линии тока раствора приведены на рис.4.

е

Из рис. 4 следует, что в начальные моменты времени, ГК в основном определяется джоулевым нагревом раствора.

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

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

Исследование выполнено при финансовой поддержке РФФИ в рамках научных проектов: № 16-08-00128 А "Теоретическое и экспериментальное исследование гравитационной конвекции в мембранных системах с учетом реакции диссоциации/рекомбинации молекул воды", № 16-48-230856р_а "Теоретическое и экспериментальное исследование влияния тепловых эффектов на гравитационную конвекцию в мембранных системах с амфолитами".

Библиографический список

1.Коваленко А.В. Влияние реакции диссоциации/рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 1. Математическая модель / Коваленко А.В., Уртенов М. Али Х., Сеидова Н.М., Письменский А.В. // Научный журнал КубГАУ, 2016. - №07(121). - Режим доступа: http://ei.kubagro.ru/2016/07/pdf/122.pdf

2. Коваленко А.В. Влияние реакции диссоциации/рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 2. Асимптотический анализ / Коваленко А.В., Уртенов М.-Али Х., Сеидова Н.М., Письменский А.В. // Научный журнал КубГАУ, 2016. - №08(122). - Режим доступа: http://ej .kubagro.ru/2016/08/pdf/17.pdf

3. Belashova E.D., Pismenskaya N.D., Nikonenko V.V., Sistat P., Pourcelly G. Current voltage curve of an anion exchange membrane in a monosodium phosphate solution. Modelling and experiment.

4. Волков А.И., Жарский И.М. Большой химический справочник. -Мн.: Современная школа, 2005. - 608 с.

5. Певницкая, М.В. Интенсификация массопереноса при электродиализе разбавленных растворов // Электрохимия. 1992. - Т.28, № 11. - С. 1708-1715.

6. Тихонов А.Н., Самарский А.А. Уравнения математической физики. Учебное пособие для вузов. - 5-е изд. - М.: Наука, 1977. - 735 с.: ил.

7. Nikonenko, V., Lebedev K., Manzanares J.A., Pourcelly G. // Electrochim. Acta. 2003. V.48, N 24. P.3639-3650.

References

1.Kovalenko A.V. Vlijanie reakcii dissociacii/rekombinacii molekul vody na perenos 1:1 jelektrolita v membrannyh sistemah v diffuzionnom sloe. Chast' 1. Matematicheskaja model' / Kovalenko A.V., Urtenov M. Ali H., Seidova N.M., Pis'menskij A.V. // Nauchnyj zhurnal KubGAU, 2016. - №07(121). - Rezhim dostupa: http://ej.kubagro.ru/2016/07/pdf/122.pdf

2. Kovalenko A.V. Vlijanie reakcii dissociacii/rekombinacii molekul vody na perenos 1:1 jelektrolita v membrannyh sistemah v diffuzionnom sloe. Chast' 2. Asimptoticheskij analiz / Kovalenko A.V., Urtenov M.-Ali H., Seidova N.M., Pis'menskij A.V. // Nauchnyj zhurnal KubGAU, 2016. - №08(122). - Rezhim dostupa: http://ej .kubagro.ru/2016/08/pdf/17.pdf

3. Belashova E.D., Pismenskaya N.D., Nikonenko V.V., Sistat P., Pourcelly G. Current voltage curve of an anion exchange membrane in a monosodium phosphate solution. Modelling and experiment.

4. Volkov A.I., Zharskij I.M. Bol'shoj himicheskij spravochnik. -Mn.: Sovremennaja shkola, 2005. - 608 s.

5. Pevnickaja, M.V. Intensifikacija massoperenosa pri jelektrodialize razbavlennyh rastvorov // Jelektrohimija. 1992. - T.28, № 11. - S. 1708-1715.

6. Tihonov A.N., Samarskij A.A. Uravnenija matematicheskoj fiziki. Uchebnoe posobie dlja vuzov. - 5-e izd. - M.: Nauka, 1977. - 735 s.: il.

7. Nikonenko, V., Lebedev K., Manzanares J.A., Pourcelly G. // Electrochim. Acta. 2003. V.48, N 24. P.3639-3650.