Мембраны и мембранные технологии, 2022, T. 12, № 6, стр. 438-451

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

О. А. Рыбалкина a*, К. В. Солонченко a, Д. Ю. Бутыльский a, В. В. Никоненко a, Н. Д. Письменская a

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

* E-mail: olesia93rus@mail.ru

Поступила в редакцию 21.05.2022
После доработки 18.07.2022
Принята к публикации 08.08.2022

Полный текст (PDF)

Аннотация

Впервые исследовано влияние параметров пульсирующего электрического поля (ПЭП) на значения средних за период плотностей тока при электродиализном обессоливании растворов дигидрофосфата натрия. Показано, что в случае растворов дигидрофосфата натрия закономерности воздействия ПЭП в целом совпадают с известными для растворов сильных электролитов. С использованием визуализации электроконвективных течений в обедненном растворе у поверхности анионообменной мембраны установлено, что наблюдаемая разница в поведении мембранных систем обусловлена ослаблением электроконвекции в фосфат содержащих растворах. Высказана гипотеза, что еще одной причиной наблюдаемых отличий является воздействие ПЭП на депротонирование анионов ${{{\text{H}}}_{{\text{2}}}}{\text{PO}}_{4}^{ - }$ при их попадании в объем анионообменной мембраны.

Ключевые слова: пульсирующее электрическое поле, электродиализ, фосфаты, электроконвекция

ВВЕДЕНИЕ

Электродиализ (ЭД), является перспективным методом для извлечения анионов фосфорной кислоты из промышленных [1, 2] и животноводческих [3] стоков, а также из фильтратов мусорных полигонов [4]. Эти же анионы присутствуют в молочной сыворотке, деминерализация которой с применением ЭД стала уже коммерческой технологией [5]. Известно, что использование классического ЭД, когда процесс осуществляют при постоянном токе, зачастую приводит к генерации H+, OH ионов с участием фиксированных групп [6] и к фаулингу (засорению) ионообменных мембран (ИОМ) [7] органическими компонентами перерабатываемых жидких сред, а также к образованию органических и неорганических осадков на их поверхности. Эти явления приводят к деструкции ИОМ, увеличивают энергозатраты и снижают производительность ЭД процессов [8, 9].

Для борьбы с такими нежелательными явлениями чаще всего используют реверсивный электродиализ (ЭДР) [1012]. Этот метод известен с 70-х гг. прошлого века [12]. Он заключается в периодическом изменении направления электрического поля (переполюсовке электродов) и одновременном изменении направления гидравлических потоков в ЭД аппарате. В этом случае камеры концентрирования ЭДР становятся камерами обессоливания и наоборот. Такой прием способствует частичному или полному удалению осадков с поверхности ИОМ. В то же время применение ЭДР приводит к значительным потерям ценных целевых компонентов.

Ю.В. Карлин и В.Н. Кропотов [13] в 1995 г. предложили использовать в ЭД пульсирующее электрическое поле (ПЭП), в котором электрические импульсы постоянного тока чередуются с паузой, в течение которой ток равен нолю. Исследования последних лет позволяют заключить, что воздействие ПЭП в какой-то мере сходно с наблюдаемым в случае ЭДР. Во время паузы (когда ток равен нулю) обедненный и обогащенный диффузионные пограничные слои (ДПС) в растворе у поверхности ИОМ частично релаксируют (“размываются”). Приповерхностные концентрации электролитов со стороны обогащенного ДПС – уменьшаются, т.е. выходят из “опасной зоны”, в которой достигается локальное превышение произведения растворимости некоторых солей, вызывающее осадкообразование [14]. Приповерхностные концентрации электролитов со стороны обедненного ДПС увеличиваются, вызывая ослабление генерации H+, OH ионов на границе ИОМ/обедненный раствор. Локальные сдвиги pH у поверхности ИОМ по сравнению с объемом раствора снижаются, что ведет к ослаблению электростатических взаимодействий некоторых отравляющих веществ, например, белков [15] с поверхностью ИОМ, а также к растворению некоторых (чувствительных к pH) осадков [9]. Кроме того, колебания приповерхностной концентрации электролитов создают электрическую неоднородность, которая способствует более быстрому (и более интенсивному), чем в случае непрерывного протекания электрического тока, развитию электроконвекции [16]. Последняя перемешивает обедненный раствор, обеспечивая дополнительную доставку целевого компонента к поверхности ИОМ и препятствуя осадкообразованию. Формирование и релаксация концентрационных профилей вблизи ИОМ является динамическим процессом. Поэтому соотношение длительности периодов импульса электрического поля и паузы (скважность), а также амплитуда этих импульсов во многом определяют количество переносимых зарядов, энергопотребление и эффективность противодействия фаулингу ИОМ [1720].

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

Заметим, что механизмы воздействия ПЭП на массоперенос и энергозатраты в основном исследовались применительно к растворам сильных электролитов [9, 24, 25]. Вместе с тем в последнее время появляется все больше исследований, согласно которым оптимальные параметры ПЭП во многом зависят от природы веществ, находящихся в перерабатываемом растворе. Например, Lemay и соавторы показали [23], что при деминерализации сладкой молочной сыворотки высокочастотный режим ПЭП, в котором длительность импульса тока и паузы в цикле составляет 0.1 с – 0.1 с (5 Гц), препятствует осадкообразованию солей кальция и магния у поверхности катионообменной мембраны (КОМ) и требует почти такого же времени, включая продолжительность пауз между импульсами тока, что и режим непрерывного постоянного тока для достижения конечной скорости деминерализации 70%. Такой режим ПЭП был признан более предпочтительным по сравнению с режимом 1 с – 1 с, при котором электроконвективные вихри в растворе у поверхности КОМ успевали “рассосаться” в период паузы. В случае безреагентного подкисления казеинат содержащих растворов, осуществляемого ЭД с биполярными мембранами, режим ПЭП 10 с–50 с в комбинации с увеличением скорости прокачки рабочего раствора обеспечивает минимальный фаулинг поверхности катионообменного слоя биполярной мембраны [26].

Несмотря на быстро возрастающее количество публикаций и расширение объектов исследований, открытым остается вопрос и об энергозатратах, требующихся для реализации ЭД в режиме ПЭП [27]. Некоторые исследователи считают, что применение ПЭП для достижения заданной деминерализации растворов приводит к снижению потребления энергии [28, 29] или, по крайней мере, остается на том же уровне [23], что и в случае классического ЭД. Согласно другим работам [30, 31] суммарные энергозатраты (и время, требующееся для деминерализации), значительно увеличиваются.

Заметим, что многие из перерабатываемых с применением ПЭП жидких сред (молочные сыворотки, бычий сывороточный альбумин, прошедшие предварительную обработку жидкие дижестаты (остатки после анаэробного сбраживания биоразлагаемого сырья) продуктов жизнедеятельности человека и животных [3235]) содержат кислотные остатки фосфорной кислоты. Известно, что поведение фосфатов в системах с ионообменными мембранами в значительной мере отличается от поведения сильных электролитов [36, 37]. Вместе с тем воздействие ПЭП на массоперенос в растворах, которые содержат только соли фосфорной кислоты, пока никем не исследовалось.

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

ОБЪЕКТЫ И МЕТОДЫ ИССЛЕДОВАНИЯ

2.1. Мембраны и растворы

Эксперименты были выполнены с использованием гетерогенной анионообменной мембраны МА-41П (ОАО Щекиноазот). Ее характеристики подробно описаны в [38]. Мембрана МА-41П выбрана, потому что она (а) обладает более высокой химической стабильностью в фосфат содержащих растворах по сравнению с гомогенными мембранами, например, AMX или AMX-Sb (Astom, Япония) [39]; (б) высокая доля ее проводящей поверхности (0.26 в набухшем состоянии), способствуют развитию электроконвекции в растворах сильных электролитов (NaCl) [38]. В качестве вспомогательной использовалась гетерогенная катионообменная мембрана МА-41 (ОАО Щекиноазот). Ее характеристики подробно описаны в [40].

0.02 М растворы NaCl (pH 5.5) и NaH2PO4 (pH 4.4), готовили из кристаллических солей ч. д. а (производитель Вектон, Россия) и дистиллированной воды (электропроводность 1.0 ± 0.1 мкСм/см, pH 5.6).

2.2. Установка для исследования параметров ПЭП

Эксперименты проводили с использованием экспериментальной установки и проточной четырехкамерной ячейки (рис. 1), подробно описанных во многих статьях, например, в [41]. Принципиальная схема этой установки представлена на рис. 1. 0.02 М раствор NaH2PO4 прокачивали через все камеры лабораторной электродиализной ячейки со средней линейной скоростью, V, равной 0.4 cм/с. Межмембранное расстояние, h, составляло 0.66 см. Длина пути обессоливания, L, была 2.0 см.

Рис. 1.

Принципиальная схема установки для изучения массообменных характеристик процесса электродиализного обессоливания 0.02 М раствора NaH2PO4: АОМ* – исследуемая анионообменная мембрана; КОМ, АОМ – вспомогательные катионообменная (МК-40) и анионообменная (МА-41П) мембраны; 1 – емкость с исследуемым раствором; 2 – капилляры Луггина, соединенные с микроемкостями, в которые погружены измерительные Ag/AgCl электроды – 3; 4 – поляризующие платиновые электроды; 5 – электрохимический комплекс Autolab-100; 6 – комбинированный электрод (LE438), соединенный с рН-метром FEP20 FiveEasy™ Plus (Mettler-Toledo) – 7; 8 – кондуктометрическая ячейка (LE703), соединенная с кондуктометром FEP30 FiveEasy™ Plus (производитель Mettler-Toledo) – 9.

Объем растворов, циркулировавших через каждый тракт электродиализной ячейки равнялся 1000.0 см3. Электропроводность, pH и температуру этого раствора контролировали через равные промежутки времени с использованием комбинированного электрода (LE438), соединенного с рН-метром FEP20 FiveEasy™ Plus (Mettler-Toledo), и кондуктометрической ячейки (LE703), соединенной с кондуктометром FEP30 FiveEasy™ Plus (производитель Mettler-Toledo).

В случае получения вольтамперных характеристик (ВАХ) скорость развертки постоянного тока составляла 0.02 мА/с. Скачок потенциала регистрировали с использованием капилляров Луггина (5), соединенных с микроемкостями, в которые были помещены измерительные Ag/AgCl электроды (6). Расстояние от поверхности анионообменной мембраны (АОМ) до кончика капилляра Луггина составляло 0.08 см. Электрическое поле, непрерывное постоянное (НЭП) и пульсирующее (ПЭП), задавали с использованием поляризующих платиновых электродов, соединенных с электрохимическим комплексом Autolab-100 (производитель Metrohm Autolab BV, Недерланды).

В случае НЭП плотность тока, iDC, задавалась, а скачок потенциала, ΔφDC, равнялся измеряемому значению скачка потенциала между капиллярами Луггина в заданный момент времени. Плотность тока определяли как i = I/S, где I – задаваемая сила тока, S – площадь поляризуемой поверхности ИОМ. Пульсирующее электрическое поле представляло собой чередование периодов, каждый из которых состоял из импульса постоянного скачка потенциала длительностью Ton и паузы (Δφ = 0) длительностью Toff. Частоту ПЭП определяли как, f = 1/T, где T = Ton+ Toff (рис. 2).

Рис. 2.

Зависимость скачков потенциала и соответствующих им плотностей тока от времени в режиме ПЭП, когда задается скачок потенциала. Пунктирными линиями указаны средние за период значения $\Delta {{\varphi }_{{av}}}$ и ${{i}_{{av}}}$.

Значения коэффициента заполнения периода (скважность) ПЭП, α, находили как: α = Ton/T. Среднее значение скачка потенциала при применении ПЭП определяли, как Δφav = Δφα, где Δφ – это скачок потенциала, задаваемого в момент импульса электрического поля (в период Ton). Например, чтобы в режиме ПЭП (α = 1/2) обеспечить среднее значение скачка потенциала Δφav = ΔφDC = 0.4 В, задаваемое значение Δφ должна равняться 0.8 В.

Диапазон параметров ПЭП (частоты от 0.01 до 20 Гц и скважности 3/4, 2/3, 1/2, 1/3 и 1/4) был аналогичен исследованному в работах [24, 25], выполненных в растворах NaCl.

Среднюю за период плотность тока для ПЭП с задаваемым значением Δφ рассчитывали по формуле:

(1)
${{i}_{{av}}} = \frac{1}{Т}\int\limits_0^t {idt} .$

2.3. Визуализация электроконвективных течений

Для визуализации электроконвективных течений (ЭК) у поверхности анионообменной мембраны (АОМ) и параллельного получения хронопотенциограмм использовали экспериментальную установку, электродиализную ячейку и методику, описанную в [37]. Необходимость задания плотности тока и регистрации скачка потенциала была предопределена особенностями визуализации ЭК вихрей. Принципиальная схема этой установки и электродиализной ячейки аналогична представленной на рис. 1. Однако электродиализная ячейка имеет другие параметры: h = 0.32 см, L = 0.53 см, расстояние между исследуемой мембраной и кончиком капилляра Луггина – 0.15 см. В 0.02 М растворы NaCl или NaH2PO4, прокачиваемые со средней линейной скоростью 0.07 см/с через камеру обессоливания (КО) и камеру концентрирования (КК) добавлен 10 μM Родамина 6G (RG6), диаметр частиц которого равен 16 Å. Известно [42], что в нейтральных и кислых растворах, R6G диссоциирует, образуя анионы Cl и катионы R6G+. Последние флуоресцируют в диапазоне длин волн 540–630 нм, придавая раствору, обогащенному катионами R6G+ (и Na+) светло-серый цвет на видео, регистрируемых с использованием CMOS камеры оптического микроскопа SOPTOP CX40M. Раствор, в котором эти ионы отсутствуют, имеет черный цвет.

В этих экспериментах значения средней за период плотности тока определяли, как ${{i}_{{av}}}$ = iα, где i – это плотность тока, задаваемого в момент импульса электрического поля (в период Ton). Например, чтобы в режиме ПЭП (α = 1/2) обеспечить среднее за период значение плотности тока ${{i}_{{av}}}$ = iDC = 5A/см2 задаваемое значение i должна равняться 10 A/см2.

Среднее за период значение скачка потенциала находили по формуле:

(2)
$\Delta {{\varphi }_{{av}}} = \frac{1}{Т}\int\limits_0^t {\Delta \varphi dt} .$

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

3.1. Вольтамперные характеристики и механизмы транспорта фосфатов

На рис. 3 представлена вольтамперная характеристика мембраны МА-41П в 0.02 М растворе NaH2PO4. Штрихпунктирная линия – это плотность тока, равная значению предельного тока $\left( {i_{{{\text{lim}}}}^{{{\text{Lev}}}}} \right),$ рассчитанному по уравнению Левека:

(3)
$i_{{{\text{lim}}}}^{{{\text{Lev}}}} = \frac{{FDC}}{{h\left( {{{T}_{1}} - {{t}_{1}}} \right)}}\left[ {1.47{{{\left( {\frac{{{{h}^{2}}V}}{{LD}}} \right)}}^{{1{\text{/}}3}}}} \right].$
Здесь F – постоянная Фарадея, D и t1 – коэффициент диффузии электролита и электромиграционное число переноса противоиона в бесконечно разбавленном растворе, С – молярная концентрация электролита в растворе. Число переноса противоионов в мембране, T1, считалось равным единице.

Рис. 3.

Вольтамперная характеристика мембраны МА-41П в 0.02 М NaH2PO4 (а) и фрагмент той же кривой, из которой вычтена омическая составляющая скачка потенциала (б). Штрихпунктирная линия соответствует предельному току, $i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ = 1.62 мА/см2, рассчитанному по уравнению (5). Точечными линиями обозначены значения скачков потенциала ΔφDC и соответствующих им плотностей токов, iDC, с которыми сравниваются данные, полученные в режиме ПЭП при Δφ = const = 1.61 В (обозначен пунктирной линией).

Форма ВАХ является типичной для анионообменных мембран в растворах NaH2PO4 [39,41]. Кривая имеет слабо выраженное наклонное плато I (рис. 3б) в области $i_{{\lim }}^{{{\text{Lev}}}}$. Оно соответствует достижению минимальных концентраций электролита в обедненном растворе, граничащем с МА-41П, и максимальному парциальному току анионов ${{{\text{H}}}_{2}}{\text{PO}}_{4}^{ - }$ в мембране [36]. Рост тока сверх $i_{{\lim }}^{{{\text{Lev}}}}$ обусловлен частичным депротонированием анионов ${{{\text{H}}}_{2}}{\text{PO}}_{4}^{ - }$, попадающих в МА-41П, результатом которого является превращение однозарядных анионов в двухзарядные анионы ${\text{HPO}}_{4}^{{2 - }}$ и доннановское исключение протонов в обедненный раствор. Механизм этого явления, которое для краткости далее обозначается термином “диссоциация кислоты” представлен на рис. 4б.Хорошо выраженное плато II (рис. 3а) регистрируется примерно при двукратном превышении $\left( {i_{{{\text{lim}}}}^{{{\text{Lev}}}}} \right).$ Оно соответствует ситуации, при которой подавляющее большинство анионов ${{{\text{H}}}_{{\text{2}}}}{\text{PO}}_{4}^{ - }$ уже подверглись депротонированию и переносятся через МА-41П в виде двухзарядных анионов. Поэтому дальнейшее увеличение скачка потенциала не вызывает заметного приращения плотности тока. Следующий за плато II заметный рост плотности тока в основном вызван генерацией H+, OH ионов на границе мембрана/раствор с участием фиксированных групп МА-41П [36], а также электроконвективным перемешиванием обедненного раствора. Как показано в работах [37, 43] электроконвекция в растворах NaH2PO4 развивается гораздо слабее по сравнению с растворами NaCl (рис. 4а).

Рис. 4.

Схема переноса анионов и генерации протонов в системах АОМ/раствор NaCl (а) или раствор NaH2PO4 (б). Цифрами обозначены: генерация H+, OH ионов с участием фиксированных групп АОМ (1); электроконвекция (2); генерация H+ ионов благодаря депротонированию анионов фосфорной кислоты (“диссоциация кислоты”) (3). Пояснения в тексте. Адаптировано из [43].

Причиной ослабления ЭК в фосфат содержащих растворах является более существенное, чем в растворах NaCl, снижение области пространственного заряда (ОПЗ). Действительно, в растворах NaCl протоны, поступление которых от поверхности АОМ в обедненный раствор вызывает снижение ОПЗ, генерируются только благодаря диссоциации воды с участием фиксированных групп мембраны. В фосфатсодержащих растворах к этому механизму добавляется “диссоциация кислоты”. Знание этих различий в поведении мембранных систем в растворах сильных электролитов (NaCl) и амфолитов (анионы фосфорной кислоты) позволяет объяснить экспериментальные результаты, представленные на рис. 5–8.

Рис. 5.

Частотные зависимости значений средних за период ПЭП плотностей тока, полученные для мембраны МА-41П в 0.02 М растворе NaH2PO4 при скачке потенциала в момент импульса электрического поля Δφ = 1.61 В. Маркеры соответствуют значениям ${{i}_{{av}}}$, полученным при скважностях 1/4 (розовый цвет); 1/3 (голубой цвет); 1/2 (красный цвет); 2/3 (фиолетовый цвет); 3/4 (зеленый цвет). Точечные линии тех же цветов – это плотности тока для каждой из скважностей ПЭП (табл. 1), которые соответствуют значениям ΔφDC на ВАХ (режим НЭП), равных $\Delta {{\varphi }_{{av}}}$. Штрихпунктирная линия – предельный ток, рассчитанный по уравнению (3). Пояснения в тексте.

Рис. 6.

Хронопотенциограммы и визуализация ЭК вихревых структур в обедненном растворе у поверхности мембраны AMX, полученные в режиме НЭП (а), а также в режиме ПЭП при α = 1/2, f = 0.5 Гц (б) и α = 1/2, f = 5.0 Гц (в). Исследования выполнены в 0.02 М растворе NaCl при плотности тока в момент импульса электрического поля 10.11 мА/см2 (${i \mathord{\left/ {\vphantom {i {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}} \right. \kern-0em} {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}$ = 3.0); $i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ = 3.39 мА/см2 рассчитан по уравнению (3).

Рис. 7.

Хронопотенциограммы и визуализация ЭК вихревых структур в обедненном растворе у поверхности мембраны AMX, полученные в режиме НЭП (а), а также в режиме ПЭП при α = 1/2, f = 0.5 Гц (б) и α = 1/2, f = 5.0 Гц (в). Исследования выполнены в 0.02 М растворе NaH2PO4 при плотности тока в момент импульса электрического поля 10.11 мА/см2 (${i \mathord{\left/ {\vphantom {i {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}} \right. \kern-0em} {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}$ = 5.5); $i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ = 1.85 мА/см2, рассчитан по уравнению (3).

Рис. 8.

Частотные зависимости средних за период ПЭП значений плотностей тока, полученные для мембраны МА-41П в 0.02 М растворе NaH2PO4 при средних за период скачках потенциала 1.00 В (а), 1.60 В (б) и 2.19 В (в). Эти значения соответствуют токам iDC = $1.36i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ (а), $1.85i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ (б) и $2.16i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ (в), которые обозначены пунктирной линией черного цвета. Процедура определения этих токов на ВАХ показана на вкладках. Цифры у маркеров обозначают скважность периода. Пояснения в тексте.

3.2. Пульсирующее электрическое поле: Δφ = const (i = const)

Рисунок 5 демонстрирует зависимости средних за период значений плотностей тока от частоты накладываемого ПЭП, полученные для мембраны МА-41П в 0.02 М растворе NaH2PO4 при фиксированных значениях скачка потенциала в момент импульса электрического поля, Δφ = 1.61 В, но разных величинах скважностей периода. Соответствующие им значения $\Delta {{\varphi }_{{av}}}$ при исследованных величинах скважностей суммированы в табл. 1. В случае скважностей 1/4 и 1/3 плотности тока, iDC, отвечающие $\Delta {{\varphi }_{{av}}}$ = ΔφDC на ВАХ (рис. 2), находятся в области допредельных токов (i < $i_{{{\text{lim}}}}^{{{\text{Lev}}}}$). При α = 1/2 эта плотность тока немногим превышает $i_{{{\text{lim}}}}^{{{\text{Lev}}}}.$ Для значений скважностей 2/3 и 3/4 – это область плотностей токов $i_{{{\text{lim}}}}^{{{\text{Lev}}}} \ll i < i_{{{\text{lim}}2}}^{{{\text{exp}}}}.$

Таблица 1.  

Средние значения скачка потенциала за период ПЭП при исследованных скважностях и соответствующие им значения постоянного тока на ВАХ (рис. 3), достигаемые при скачке потенциала 1.61 В в момент импульса электрического поля

α 1/4 1/3 1/2 2/3 3/4
$\Delta {{\varphi }_{{av}}}$, В 0.40 0.53 0.80 1.06 1.20
iDC 0.94 1.21 1.79 2.33 2.56
$i_{{av}}^{{{\text{max}}}}$ 1.18 1.37 1.93 2.33 2.49
$i_{{av}}^{{{\text{max}}}}$/iDC 1.25 1.13 1.08 1.00 0.97
$i_{{av}}^{{{\text{min}}}}$ 0.99 1.26 1.56 1.88 2.27
$i_{{av}}^{{{\text{min}}}}$/iDC 1.05 1.04 0.87 0.80 0.89

На первый взгляд представленные на рис. 5 зависимости в целом подчиняются закономерностям, известным для растворов сильных электролитов, например, NaCl [44]. Значения ${{i}_{{av}}}$ растут с увеличением Ton, т.е. в ряду 1/4 < 1/3 < 1/2 < 2/3 < 3/4. При скважностях 1/4 и 1/3 значения iav растут в области низких частот, но при более высоких частотах (f > 0.5–0.8) достигают максимума и перестают существенно зависеть от частоты. При скважностях 2/3 и 3/4 средняя за период плотность тока принимает максимальное значение в области малых частот и становится минимальной в области больших частот. Скважность α = 1/2, по-видимому, является переходным вариантом: максимальные значения ${{i}_{{av}}}$ наблюдаются в области средних частот 0.45 < f < 2.3.

Такой вид зависимостей ${{i}_{{av}}}$ (f, α), согласно исследованию Ибла [45], определяется природой воздействия ПЭП на обедненный диффузионный пограничный слой (ДПС). Его можно условно разделить на две зоны. Зона ДПС, граничащая с объемом раствора, является стационарной, в то время как зона у поверхности мембраны является осциллирующей. В момент импульса электрического поля концентрация электролита в осциллирующей зоне существенно снижается. Во время паузы имеет место частичная релаксация концентрационного профиля, которая приводит к существенному снижению сопротивления примембранного раствора. В начале следующего периода ПЭП Δφ прикладывается к системе с более низким сопротивлением, чем в случае НЭП. Поэтому достигаемое значение ${{i}_{{av}}}$ превышает iDC при $\Delta {{\varphi }_{{av}}}$ = ΔφDC. Средний за период ПЭП выигрыш в скорости массопереноса по сравнению с НЭП обусловлен высоким вкладом электромиграционного переноса при повторном задании Δφ. Ширина осциллирующей зоны обедненного ДПС и степень его релаксации растут с увеличением Toff, которая определяет значения f и α (рис. 2). Верхний предел ${{i}_{{av}}}$ (скорости массопереноса) достигается, когда толщина осциллирующей зоны ДПС становится пренебрежимо малой по сравнению с его полной толщиной. Нижний предел ${{i}_{{av}}}$ достигается, когда период колебаний осциллирующей зоны оказывается много больше времени релаксации концентрационного профиля [24].

Следует заметить, что превышение ${{i}_{{av}}}$ над iDC (более значительное в области ${{i}_{{av}}} \ll i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ и менее существенное при ${{i}_{{av}}}$ близких к $i_{{{\text{lim}}}}^{{{\text{Lev}}}}$) является характерным и в случае раствора NaCl [44], и в случае раствора NaH2PO4 (рис. 5). Различия в поведении этих систем возникают, если ${{i}_{{av}}} \gg i_{{{\text{lim}}}}^{{{\text{Lev}}}}.$ Действительно, в растворе NaCl (мембрана Nafion-117) [43] удается достичь максимального (до 33%) прироста ${{i}_{{av}}}$ над iDC по сравнению с допредельными токовыми режимами, если α = 1/2; при α > 1/2 максимальные и минимальные значения ${{i}_{{av}}}$ всегда превышают iDC или близки к нему. В растворе NaH2PO4 (рис. 5) ${{i_{{av}}^{{{\text{max}}}}} \mathord{\left/ {\vphantom {{i_{{av}}^{{{\text{max}}}}} {{{i}_{{DC \approx 1}}}}}} \right. \kern-0em} {{{i}_{{DC \approx 1}}}}},$ а $i_{{av}}^{{{\text{min}}}} < {{i}_{{DC}}}$ (табл. 1).

Чтобы понять причины наблюдаемых различий, сравним параметры ЭК вихрей, получаемых в режимах НЭП и ПЭП в растворах NaCl (рис. 6) и NaH2PO4 (рис. 7) при одинаковой задаваемой плотности тока в момент импульса электрического поля, равной 0.65 мА/см2. Как было показано в предшествующих исследованиях [1, 43], при таком токе в режиме НЭП ЭК вихри отчетливо визуализируются в обеих исследуемых системах. В случае ПЭП частоты 0.5 и 5 Гц при α = 1/2 обеспечивают оптимальное увеличение выходов по току сильных электролитов при сохранении тех же энергозатрат, что и в режиме НЭП [17, 23].

Закономерности развития ЭК в системах ионообменная мембрана/раствор сильного электролита хорошо изучены. Результаты этих исследований обобщены в нескольких обзорах, например, в [46]. Наибольшее воздействие на распределение концентраций электролита в обедненном растворе у поверхности мембраны оказывает неравновесная ЭК, которую называют электроконвекцией Рубинштейна–Зальцмана. Она характеризуется образованием кластеров крупных вихрей и развивается при i > $i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ в пороговом режиме после снижения концентрации электролита в примембранном растворе до минимальных значений. Причем, зона примембранного раствора, перемешиваемая ЭК вихрями расширяется с увеличением ${i \mathord{\left/ {\vphantom {i {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}} \right. \kern-0em} {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}.$

Благоприятные для старта ЭК значения скачка потенциала в режиме НЭП достигаются в точке 2 хронопотенциограммы, ХП, (рис. 6а) после включения тока [42]. Однако максимальные (для заданного тока) размеры вихрей формируются только в точке 3 (то есть примерно через 24 с с момента включения тока), обеспечивая стационарные значения регистрируемого на ХП скачка потенциала. О наличии этих вихрей свидетельствуют осцилляции Δφ на ХП, а также каплевидные темные области на кадрах видео. Согласно расчетам, выполненным с использованием 2D модели [1, 42], эти наиболее затемненные каплеобразные области характеризуются наименьшей скоростью вращения ЭК вихрей и поэтому являются наиболее обессоленными.

Из результатов математического моделирования [16] следует, что в случае ПЭП после выключения тока кластеры крупных ЭК вихрей быстро (в течение 0.02 с) затухают, потому что исчезает движущая сила – электрическое поле. Однако остаточные слабые вихри продолжают существовать достаточно длительное время благодаря сохранению электроосмотических потоков жидкости. Эти потоки обусловлены неравномерным распределением концентрации противоионов внутри двойного электрического слоя на границе АОМ/раствор. В результате неравномерное распределение объемной электрической силы подпитывает ЭК вихри [16].

Время релаксации концентрационного поля может составлять несколько секунд. Поэтому подпитка ЭК вихрей химической энергией этого неоднородного электрического поля может длиться достаточно долго. После наложения следующего импульса электрического поля после паузы ЭК может восстановиться значительно быстрее, чем после включения тока в режиме НЭП. Причиной является неоднородность остаточного поля концентраций, которая усиливает тангенциальную составляющую электрического поля, то есть действует аналогично электрической неоднородности поверхности мембраны [47]. Действительно, в случае f = 0.5 Гц и α = 1/2 (Ton = 1 с, Toff = 1 c) через тот же интервал времени с момента начала эксперимента, как и в случае НЭП (рис. 6а) осцилляции Δφ на ХП появляются уже через 0.2 с после включения тока и не исчезают в течение паузы (рис. 6б). Вместе с тем размеры кластеров вихрей (вернее их наиболее обессоленной части) в момент включения тока оказываются меньше, а значение Δφ – выше по сравнению с режимом НЭП. По-видимому, за столь короткое время отдельные ЭК вихри не успевают объединиться в более крупные кластеры. К завершению паузы обедненный ДПС не релаксирует полностью: толщина зоны с неоднородным распределением концентраций электролита у поверхности АОМ оказывается сопоставимой с толщиной обедненного ДПС в момент импульса тока.

Увеличение частоты до значения 5 Гц (Ton = 0.1 с, Toff = 0.1 c), то есть на порядок (рис. 6в), позволяет сохранить примерно те же значения Δφ, что и в случае, представленном на рис. 6б. Однако, распределение концентраций NaCl в обедненном ДПС в момент импульса и в момент паузы является более равномерным по сравнению с рис. 6а, 6б. По-видимому, с одной стороны, быстрая смена импульсов тока и пауз не позволяет сформироваться обессоленному ядру кластеров вихрей. С другой стороны, релаксация обедненного ДПС находится на начальной стадии, поэтому Δφ в момент паузы не достигают нулевых значений (рис. 6в).

Следует подчеркнуть, что задаваемая плотность тока в момент импульса составляет $3.0i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ (NaCl) и $5.5i_{{{\text{lim}}}}^{{{\text{Lev}}}}$ (NaH2PO4). Несмотря на значительное превышение ${i \mathord{\left/ {\vphantom {i {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}} \right. \kern-0em} {i_{{{\text{lim}}}}^{{{\text{Lev}}}}}}$ размеры вихрей в системе с раствором NaH2PO4 оказываются меньше (рис. 7а) а интенсивность их окраски глубже по сравнению с раствором NaCl. Эта разница обусловлена более интенсивным поступлением протонов от мембраны в обедненный раствор NaH2PO4. Действительно, в случае NaCl (рис. 4а) источником этих протонов является генерация H+, OH ионов с участием фиксированных групп, тогда как в растворах NaH2PO4 к этому механизму присоединяется механизм “диссоциации кислоты” (рис. 4б). Катионы H+, поступающие в примембранный раствор, уменьшают (по абсолютной величине) плотность объемного электрического заряда, знак которого изначально формируют противоионы ${{{\text{H}}}_{{\text{2}}}}{\text{PO}}_{4}^{ - }.$

Условия для развития ЭК ухудшаются; вихри вращаются медленнее, а раствор в их ядре обессоливается сильнее по сравнению с NaCl. Меньшая подвижность анионов ${{{\text{H}}}_{{\text{2}}}}{\text{PO}}_{4}^{ - }$ по сравнению с анионами Cl является причиной более высокого сопротивления раствора NaH2PO4 по сравнению с NaCl. Это обстоятельство обеспечивает более значительный омический скачок потенциала (φΩ) в момент включения тока, а также более медленную релаксацию обедненного ДПС во время паузы (рис. 6, 7). В пользу последнего обстоятельства свидетельствуют более быстрое появление осцилляций Δφ и ЭК вихрей в момент включения тока (которое регистрируется на видео); более значительные осцилляции скачка потенциала и более интенсивная окраска примембранного раствора NaH2PO4 во время паузы. В остальном закономерности развития электроконвекции в обеих исследованных системах являются сходными.

Более слабое развитие ЭК в момент импульса тока (скачка потенциала) может быть причиной менее существенного приращения ${{i}_{{av}}}$ по сравнению с iDC в растворах NaH2PO4 по сравнению с NaCl, наблюдаемого при скважностях 1/4, 1/3 и 1/2 (рис. 7). Однако, природа снижения ${{i}_{{av}}}$ по сравнению с iDC при α = 2/3 и 3/4 остается не вполне понятной. Чтобы объяснить наблюдаемый ход кривых, необходимо принять во внимание, что в режиме ПЭП увеличение примембранной концентрации NaH2PO4 благодаря облегчению развития ЭК и/или частичной релаксации обедненного слоя должно способствовать росту ${{i}_{{av}}}$ (как и в случае сильных электролитов). Вместе с тем подавление превращения анионов ${{{\text{H}}}_{{\text{2}}}}{\text{PO}}_{4}^{ - }$ в двухзарядные анионы ${\text{HPO}}_{4}^{{2 - }}$ (рис. 4б), напротив, может способствовать снижению ${{i}_{{av}}}$ по сравнению с iDC. Причем такое снижение можно расценивать как положительный эффект от применения ПЭП, потому что его результатом может стать рост выхода по току по целевому компоненту – фосфору [36]. Действительно, один атом фосфора содержится в любом анионе фосфорной кислоты. Поэтому превращение всех однозарядных анионов в двухзарядные увеличит регистрируемый ток в 2 раза, при этом полезный массоперенос останется неизменным. По-видимому, наблюдаемые нами закономерности воздействия ПЭП на частотные зависимости средних за период плотностей тока (рис. 5) являются результирующей этих двух разнонаправленных факторов. Причем, именно повышение концентрации NaH2PO4 у поверхности АОМ благодаря развитию ЭК, вероятно, способствует сдерживанию превращения однозарядных анионов фосфорной кислоты в двухзарядные. Этот фактор должен в большей степени проявляться при высоких значениях скачков потенциала (плотностей тока) в момент импульса электрического поля. В следующем разделе представлены результаты исследования именно таких интенсивных токовых режимов.

3.3. Пульсирующее электрическое поле в интенсивных токовых режимах: Δφav = const

На рис. 8 показаны зависимости средних за период значений плотностей тока от частоты накладываемого ПЭП, полученные для мембраны МА-41П в 0.02 М растворе NaH2PO4 при средних за период скачках потенциала, $\Delta {{\varphi }_{{av}}}$, равных 1.00, 1.60 и 2.19 В. Соответствующие им значения Δφ, при исследованных величинах скважностей периодов ПЭП суммированы в табл. 2.

Таблица 2.  

Значения скачков потенциала в момент импульса электрического тока, Δφ, обеспечивающих достижение исследованных величин $\Delta {{\varphi }_{{av}}}$ и α

Δφav, В α
1/4 1/3 1/2 2/3 3/4
1.00 4.00 3.00 2.00 1.50 1.33
1.60 6.44 4.83 3.22 2.41 2.15
2.19 8.76 6.57 4.38 3.29 2.92

Как видно из рис. 8а, описывающего случай $\Delta {{\varphi }_{{av}}}$ = 1.0 В $\left( {i_{{{\text{lim}}}}^{{{\text{Lev}}}} < {{i}_{{DC}}} < i_{{{\text{lim}}2}}^{{{\text{exp}}}}} \right),$ ход частотных зависимостей средней за период плотности тока имеет следующий вид. При α = 1/4 значение ${{i}_{{av}}}$ в режиме ПЭП растет с увеличением частоты и достигает максимума, равного 0.4iDC, при f ≈ 0.05 Гц (T ≈ 16 c). Увеличение α до 1/3 приводит к росту численных значений средних за период плотностей тока и к изменению хода кривой: ${{i}_{{av}}}$ медленно нарастает с увеличением частоты ПЭП и достигает постоянных значений, равных 0.68iDC, при f ≈ 0.36 Гц (T ≈ 0.7 c). Увеличение α до 1/2 и 2/3 сопровождается дальнейшим ростом ${{i}_{{av}}}$, причем максимальные значения, равные (0.9–0.95)iDC наблюдаются в диапазоне частот 0.2–2.0 Гц (5 > T > 0.5 c). Использование α = 3/4 дает примерно те же ${{i}_{{av}}}$, что и в случае скважностей 1/2 и 2/3, однако максимальные значения средней за период плотности тока достигаются при самой низкой из исследованных частот (f = 0.01 Гц, T = 100 c), затем немного снижаются и перестают изменяться при f > –0.25 Гц (T < 4 c).

При Δφav = 1.6 В $\left( {{{i}_{{DC}}} \approx i_{{{\text{lim}}2}}^{{{\text{exp}}}}} \right),$ рост ${{i}_{{av}}}$ при увеличении скважности периода ПЭП сохраняется (рис. 8б). Однако, значения всех средних за период плотностей тока находятся в интервале от (0.43 до 0.56) iDC и демонстрируют слабую зависимость от частоты ПЭП. Исключение составляет только кривая для α = 3/4, которая характеризуется небольшим снижением ${{i}_{{av}}}$ при f > 0.05 Гц (T < 20 c). Если же $\Delta {{\varphi }_{{av}}}$ = 2.19 $\left( {{{i}_{{DC}}} > i_{{{\text{lim}}2}}^{{{\text{exp}}}}} \right),$ эти кривые делятся на две группы (рис. 8в). Одна из них (α = 1/4 и α = 1/3) соответствует средним за период плотностям токов 0.38 < ${{i}_{{av}}}$/iDC < 0.42, другая (1/2 < α < 3/4) локализована в окрестностях ${{i}_{{av}}}$ 0.58iDC. Рост ${{i}_{{av}}}$ с увеличением частоты наблюдается только в случае α = 1/3 при f > 0.15 Гц (T < 6.7 c), а также α = 2/3 и α = 3/4 при f > 2.3 Гц (T < 0.4 c).

Таким образом, во всех представленных на рис. 8 случаях выполняется тенденция, характерная и для сильных электролитов [43]: в целом, ${{i}_{{av}}}$ растет с увеличением α. Вместе с тем, явления, связанные с протонированием-депртонированием частиц фосфорной кислоты в растворе и объеме мембраны, вносят свои коррективы в ход получаемых зависимостей.

При $\Delta {{\varphi }_{{av}}}$ = 1.00 В (рис. 8а) высокие значения скачка потенциала в момент импульса электрического поля, характерные для α = 1/4 и α = 1/3 (табл. 2) вызывают развитие электроконвекции [37], увеличивающей плотность тока в момент импульса электрического тока. По-видимому, максимум на кривой α = 1/4 (рис. 8а) соответствует частотам, при которых электроконвекция в момент импульса развивается в полной мере и ЭК вихри не успевают рассосаться в момент паузы. При более высоких частотах ЭК не успевает развиться до такой степени, чтобы обеспечить низкое сопротивление примембранного раствора в момент паузы. Следует заметить, что для высоких значений скачков потенциала в момент импульса тока (Δφ = 4.0 при α = 1/4) характерна интенсивная генерация H+, OH ионов с участием фиксированных групп, результатом которой является превращение анионов ${{{\text{H}}}_{{\text{2}}}}{\text{PO}}_{4}^{ - }$ в ${\text{PO}}_{4}^{{3 - }}$ в режиме НЭП [36, 41]. В том случае, если в режиме ПЭП удается ослабить это явление, можно ожидать, что отношение ${{i}_{{av}}}$/iDC станет меньше 0.5, что и наблюдается в эксперименте.

Снижение Δφ при α = 1/3 приводит к сокращению Toff, в течение которого идет релаксация обедненного ДПС. Поэтому в этом случае побеждает другая известная для сильных электролитов тенденция, а именно, рост ${{i}_{{av}}}$ с увеличением f [43]. Разница заключается в том, что отношение ${{i}_{{av}}}$/iDC в случае раствора NaH2PO4 оказывается ниже по сравнению с раствором NaCl. Мы предполагаем, что в паузах между импульсами электрического поля уменьшается движущая сила, способствующая переносу протонов из объема мембраны в обедненный ДПС. Результатом этого снижения является ослабление превращения однозарядных анионов фосфорной кислоты в двухзарядные и, соответственно, уменьшение регистрируемой плотности тока. При скважностях 1/2, 2/3 и 3/4 значения Δφ оказываются недостаточными для формирования крупных ЭК вихрей [37, 41]. Поэтому максимумы на кривых α = 1/2 и α = 2/3, скорее всего, обусловлены конкурированием тенденций 1 (в области малых частот ${{i}_{{av}}}$ растет с увеличением α [44]) и 2 (в области больших частот ${{i}_{{av}}}$ растет с уменьшением α [44]). В случае α = 3/4 доминирует тенденция 1. Сокращение Toff уменьшает время, в течение которого ослабляется исключение протонов из объема АОМ. Поэтому при скважностях 1/2, 2/3 и 3/4 значения ${{i}_{{av}}}$ мало отличаются от значений iDC.

Слабое влияние частот и скважностей периодов ПЭП на значения ${{i}_{{av}}}$ при $\Delta {{\varphi }_{{av}}}$ = 1.68 В (рис. 8б) можно объяснить тем, что при достигаемых в момент импульса значениях Δφ (табл. 2) конкуренция описанных выше явлений в основном приводит к ослаблению трансформации однозарядных анионов фосфорной кислоты в трехзарядные.

При увеличении $\Delta {{\varphi }_{{av}}}$ до значения 2.19 В при α = 1/4 и α = 1/3 доминирует ослабление превращения анионов фосфорной кислоты в трехзарядную форму в объеме мембраны. В случае скважностей 1/2, 2/3 и 3/4 росту ${{i}_{{av}}}$ способствует формирование кластеров крупных ЭК вихрей при наложении более длительного, чем при α = 1/4 и α = 1/3, импульса электрического поля.

Заметим, что представленная нами интерпретация хода частотных зависимостей ${{i}_{{av}}}$ в режимах ПЭП основываются на проведенных ранее исследованиях, но все же являются гипотезой. Для проверки этой гипотезы требуется осуществить математическое моделирование изучаемых процессов, а также провести эксперименты по электродиализному обессоливанию растворов NaH2PO4 и сравнить массообменные характеристики электродиализа в режимах НЭП и ПЭП. Эти исследования планируется выполнить в дальнейшем.

ЗАКЛЮЧЕНИЕ

Экспериментально установлено, что закономерности воздействия пульсирующего электрического поля (ПЭП) на систему АОМ/раствор NaH2PO4 в целом аналогичны полученным ранее для случая растворов сильных электролитов, в частности, NaCl. При постоянной величине задаваемого в момент импульса ПЭП скачка потенциала (Δφ) значения средней за период плотности тока (${{i}_{{av}}}$) растут с увеличением Ton, т.е. в ряду 1/4 < 1/3 < 1/2 < 2/3 < 3/4. При скважностях 1/4 и 1/3 iav растет в области низких частот и остается практически неизменным при более высоких частотах (f > 0.5–0.8). При скважностях 2/3 и 3/4 средняя за период плотность тока принимает максимальное значение в области малых частот и становится минимальной в области больших частот. Скважность α = 1/2, по-видимому, является переходным вариантом: максимальные значения ${{i}_{{av}}}$ наблюдаются в области средних частот: 0.45 < f < 2.3. Такой вид зависимостей ${{i}_{{av}}}$ (f, α) определяется существенным снижением сопротивления примембранного раствора благодаря частичной релаксации концентрационного профиля в момент паузы по сравнению с достигаемому в период импульса. В начале следующего периода ПЭП электрическое поле прикладывается к системе с более низким сопротивлением, чем в случае непрерывного постоянного электрического поля (НЭП). Поэтому достигаемое значение ${{i}_{{av}}}$ превышает iDC, если ${{i}_{{av}}} \leqslant i_{{{\text{lim}}}}^{{{\text{Lev}}}}.$ Вместе с тем, в случае раствора NaCl удается достичь значительно большего прироста ${{i}_{{av}}}$ над iDC по сравнению с раствором NaH2PO4.

При ${{i}_{{av}}} \gg i_{{{\text{lim}}}}^{{{\text{Lev}}}},$ в системе АОМ/раствор NaH2PO4 значения iDC всегда превышают ${{i}_{{av}}}$. Результаты визуализации электроконвективных течений в обедненном растворе у поверхности АОМ в режимах НЭП и ПЭП позволяют заключить, что наблюдаемая разница в поведении мембранных систем обусловлена ослаблением электроконвекции в фосфат содержащих растворах по сравнению с растворами NaCl. Высказана гипотеза, согласно которой еще одной причиной наблюдаемых отличий является воздействие ПЭП на депротонирование анионов ${{{\text{H}}}_{2}}{\text{PO}}_{4}^{ - }$ при их попадании в объем анионообменной мембраны. Это полезное для увеличения выходов по току явление может приводить к уменьшению регистрируемых значений ${{i}_{{av}}}$. Проверку этой гипотезы планируется выполнить в дальнейших исследованиях, сравнив степени электродиализного обессоливания растворов NaCl и NaH2PO4 в режимах НЭП и ПЭП.

Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 20-38-90054.

Список литературы

  1. Rotta E.H., Bitencourt C.S., Marder L., Bernardes A.M. // J. Memb. Sci. 2019. V. 573. P. 293–300.

  2. Melnikov S., Loza S., Sharafan M., Zabolotskiy V. // Sep. Purif. Technol. 2016. V. 157. P. 179–191.

  3. Wang X., Zhang X., Wang Y., Du Y., Feng H., Xu T. // J. Memb. Sci. 2015. V. 490. P. 65–71.

  4. Oliveira V., Dias-Ferreira C., González-García I., Labrincha J., Horta C., García-González M.C. // Waste Manag. 2021. V. 125. P. 293–302.

  5. Dufton G., Mikhaylin S., Gaaloul S., Bazinet L. // Membranes. 2020. V. 10. № 1. P. 14.

  6. Tanaka Y., Senō M. // J. Chem. Socie., Faraday Transactions 1: Physical Chemistry in Condensed Phases. 1986. V. 82. № 7. P. 2065–2077.

  7. Cifuentes-Araya N., Pourcelly G., Bazinet L. // J. Colloid Interface Sci. 2011. V. 361. № 1. P. 79–89.

  8. Malek P., Ortiz J.M., Richards B.S., Schäfer A.I. // J. Membr. Sci. 2013. V.435. P. 99–109.

  9. Mikhaylin S., Nikonenko V., Pourcelly G., Bazinet L. // J. Membr. Sci. 2014. V. 468. P. 389–399.

  10. Lee H.J., Moon S.H. // Korean J. Chem. Eng. 2004. V. 21. № 3. P. 629–634.

  11. Shi L., Xie S., Hu Z., Wu G., Morrison L., Croot P., Hu H. // J. Memb. Sci. 2019. V. 573. P. 560–569.

  12. Al-Amshawee S., Yunus M.Y.B.M., Azoddein A.A.M., Hassell D.G., Dakhil I.H., Hasan H.A. // Chem. Eng. 2020. V. 380. № 122231.

  13. Karlin Yu.V., Kropotov V.N. // Russ. J. Electrochem. 1995. V. 31. № 5. P. 472.

  14. Zhao D., Lee L.Y., Ong S.L., Chowdhury P., Siah K.B., Ng H.Y. // Sep. Purif. Technol. 2019. V. 213 P. 339–347.

  15. Bazinet L., Geoffroy T.R. // Membranes. 2020. V. 10. P. 1–72.

  16. Uzdenova A.M., Kovalenko A.V., Urtenov M.K., Nikonenko V.V. // Electrochem. Commun. 2015. V. 51. P. 1–5.

  17. Lemay N., Mikhaylin S., Bazinet L. // Innov. Food Sci. Emerg. Technol. 2019. V. 52. P. 221–231.

  18. Gonzalez-Vogel A., Rojas O.J. // Desalination 2020. P. 474. P. 114190.

  19. Dufton G., Mikhaylin S., Gaaloul S., Bazinet L. // Int. J. Mol. Sci. 2019. V. 20. P. 797.

  20. De Jaegher B., De Schepper W., Verliefde A., Nopens I.A // J. Memb. Sci. 2022. V. 642. P. 119975.

  21. Merkel A., Ashrafi A.M. // Int. J. Mol. Sci. 2020. V. 20. № 8. P. 1918.

  22. Sosa-Fernandez P.A., Post J.W., Ramdlan M.S., Leermakers F.A.M., Bruning H., Rijnaarts H.H.M. // Desalination. 2020. V. 484. P. 114424.

  23. Lemay N., Mikhaylin S., Mareev S., Pismenskaya N., Nikonenko V., Bazinet L. // J. Memb. Sci. 2020. V. 603. № 117878.

  24. Sistat P., Huguet P., Ruiz B., Pourcelly G., Mareev S.A., Nikonenko V.V. // Electrochimica Acta. 2015. V. 164. P. 267–280.

  25. Butylskii D., Moroz I., Tsygurina K., Mareev S. // Membranes. 2020. V. 10. № 3. P. 40.

  26. Nichka V., Mareev S., Pismenskaya N., Nikonenko V., Bazinet L. // Membranes. 2022. V. 12. 193.

  27. Martí-Calatayud M. C., Sancho-Cirer Poczatek M., Pérez-Herranz V. // Membranes. 2021. V. 11. № 1. P. 43.

  28. Sosa-Fernandez P., Post J., Ramdlan M., Leermakers F., Bruning H., Rijnaarts H. // Desalination. 2020. V. 484. P. 114424.

  29. Gao Q., Li Q., Lei C., Fu R., Wang W., Li Q., Liu Z. // Materials. 2020. V. 13. P. 2501.

  30. Liu R., Wang Y., Wu G., Luo J., Wang S. // Chem. Eng. 2017. P. 322. P. 224–233.

  31. Monetti J., Ledezma P., Virdis B., Freguia S. // ACS Omega. 2019. V. 4. № 1. P. 2152–2159.

  32. Zhang Y., Van der Bruggen B., Pinoy L., Meesschaert B. // J. Memb. Sci. 2009. V. 332. P. 104–112.

  33. Ghyselbrecht K., Jongbloet A., Pinoy L., Meesschaert B. // J. Environ. Chem. Eng. 2020. V. 8. № 6. P. 104492.

  34. Shi L., Xie S., Hu Z., Wu G., Morrison L., Croot P., Hu H., Zhan X. // J. Memb. Sci. 2019. V. 573. P. 560–569.

  35. Shi L., Xiao L., Hu Z., Zhan X. // Water Cycle. 2020. V. 1. P. 54–62.

  36. Rybalkina O.A., Sharafan M.V., Nikonenko V.V., Pismenskaya N.D. // J. Memb. Sci. 2022. V. 651. P. 120449.

  37. Рыбалкина О.А., Мороз И.А., Горобченко А.Д., Письменская Н.Д., Никоненко В.В. // мембраны и мембранные технологии. 2022. Т. 12. № 1. С. 38–46. [Rybalkina O.A., Moroz I.A., Gorobchenko A.D., Pismenskaya N.D., Nikonenko V.V. // Membranes and Membrane Technologies. 2022. V. 4. № 1. P. 31–38.]

  38. Pismenskaya N.D., Pokhidnia E.V., Pourcelly G., Nikonenko V.V. // J. Memb. Sci. 2018. V. 566. P. 54–68.

  39. Рыбалкина О.А., Цыгурина К.А., Сарапулова В.В., Мареев С.А., Никоненко В.В., Письменская Н.Д. // Мембраны и мембранные технологии. 2019. Т. 9. № 2. С. 131–145. [Rybalkina O.A., Tsygurina K.A., Sarapulova V.V., Mareev S.A., Nikonenko V.V., Pismenskaya N.D. // Membranes and Membrane Technologies. 2019. V. 1. № 2. P. 107–119.]

  40. Демина О.А., Кононенко Н.А., Фалина И.В., Демин А.В. // Коллоидный журн. 2017. Т. 79. № 3. С. 259–269. [Demina O.A., Kononenko N.A., Falina I.V., Demin A.V. // Colloid J. 2017. V. 79. № 3. P. 317–327.]

  41. Rybalkina O., Tsygurina K., Melnikova E., Mareev S., Moroz I., Nikonenko V., Pismenskaya N. // Int. J. Mol. Sci. 2019. V. 20. № 14. P. 3593.

  42. Magut P.K.S., Das S., Fernandv V.E., Losso J., McDonough K., Naylor B.M., Aggarwal S., Warner I.M. // J. Am. Chem. Soc. 2013. V. 135. P. 15873–15879.

  43. Pismenskaya N., Rybalkina O., Moroz I., Mareev S., Nikonenko V. // Int. J. Mol. Sci. 2021. V. 22. № 24. P. 13518.

  44. Зырянова С.В., Бутыльский Д.Ю., Мареев С.А., Письменская Н.Д., Никоненко В.В. Пурселли Ж. // Электрохимия. 2018. Т. 54. № 10. С. 885–892. [Zyryanova S.V., Butyl’skii D.Y., Mareev S.A., Pis’menskaya N.D., Nikonenko V.V., Pourcelly G. // Russ. J. Electrochem. 2018. V. 54. № 10. P. 775–781.]

  45. Ibl N. // Surface Technology. 1980. V. 10. № 2. P. 81–104.

  46. Никоненко В.В., Мареев С.А., Письменска Н.Д., Узденова А.М., Коваленко А.В., Уртенов М.Х., Пурсели Ж. // Электрохимия. 2017. Т. 53. № 10. С. 1266–1289. [Nikonenko V.V., Mareev S.A., Pis’menskaya N.D., Uzdenova A.M., Kovalenko A.V., Urtenov M.K., Pourcelly G. // Russ. J. Electrochem. 2017. V. 53. № 10. P. 1122–1144.]

  47. Rubinstein I., Zaltzman B., Kedem O. // J. Memb. Sci. 1997. V. 125. P. 17.

Дополнительные материалы отсутствуют.