Математическое моделирование процесса обессоливания в ЭДА с учетом электроконвекции и геометрической неоднородности поверхности мембраны

Размещено на http://www.allbest.ru/

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

Узденова Аминат Магометовна

Коваленко Анна Владимировна

Уртенов Махамет Хусеевич

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

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

В данной статье, являющейся продолжением работы [8] предлагается моделирование, и исследование электроконвекции в слое раствора прилегающей к катионообменной мембране камеры обессоливания электродиализного аппарата (ЭДА).

Экспериментально доказано [4], что поверхность мембраны является геометрически неоднородной. Кроме того, в настоящее время ведутся исследования по целенаправленной модификации поверхности мембраны и ее влияния на массоперенос. Установлено, что в зависимости от модифицирующего агента и микроструктуры мембраны наблюдается увеличение или уменьшение скорости массопереноса [1]. Теоретически и экспериментально установлено важная роль электроконвекции для лучшего понимания механизмов сверхпредельного тока и интенсификации массопереноса в электродиализных аппаратах [2,3]. В связи с этим возникает проблема исследования взаимодействия электроконвекции и геометрической неоднородности поверхности мембраны.

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

Будем моделировать геометрическую неоднородность поверхности катионообменной мембраны треугольными выступами и кавернами.

а)

б)

Рисунок 1 - Исследуемая область: а) мембрана с треугольными выступами и прилегающий слой раствора, б) мембрана с треугольными кавернами и прилегающий слой раствора. Раствор прокачивается слева направо с некоторой линейной скоростью

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

Пусть и - ширина и длина рассматриваемой области, соответственно, - начальная (линейная) скорость прокачивания раствора, - соответствует условной межфазной границе мембрана/раствор, - внешняя («свободная») граница раствора, - соответствует входу, а - выходу (рис. 1).

Для моделирования электроконвекции в данных условиях, будем использовать связанную систему электродиффузионных уравнений [6] и уравнений Навье-Стокса [7] в приближении Буссинеска, с учетом пространственной силы. Векторная запись этой системы для бинарного электролита, в случае отсутствия химических реакций, имеет вид:

(1)

(2)

(3)

(4)

(5)

(6)

где - градиент, - оператор Лапласа, - плотность силы электрического поля, - скорость течения раствора электролита, - характерная плотность раствора, - давление, - концентрации катионов и анионов в растворе, соответственно, - зарядовые числа катионов и анионов, - коэффициенты диффузии катионов и анионов, - потенциал электрического поля, - диэлектрическая проницаемость электролита, - постоянная Фарадея, - газовая постоянная, - абсолютная температура, - время, - коэффициенты кинематической вязкости. При этом - неизвестные функции, зависящие от времени и координат , . В системе уравнений (1)-(6) уравнения (1)-(4) описывают электрохимические поля, а уравнения Навье-Стокса (5), (6) описывают движение раствора в приближении Буссинеска под действием пространственной электрической силы. Электромембранные системы используются, как правило, в двух разных режимах работы: потенциостатическом режиме, когда постоянным поддерживается падение потенциала в цепи и гальваностатическом режиме, когда постоянным поддерживается ток в цепи. В этой работе мы будем рассматривать потенциостатический режим, которому соответствует условие:

, (7)

означающее, что величина падения потенциала в рассматриваемой области постоянна.

2. Математическая постановка задачи

После ряда преобразований уравнения (1)-(6) запишутся в виде

, (8)

, (9)

, (10)

, (11)

, (12)

Наряду с условием (7) будем использовать следующие граничные условия:

1) На поверхности катионобменной мембраны будем считать граничную концентрацию катионов равной фиксированному заряду внутри мембраны:

математический электроконвекция обессоливание электродиализный

(13)

Кроме того, предположим катионообменную мембрану идеально селективной, т.е. непроницаемой для анионов:

(14)

Для скорости используем условие прилипания:

(15)

2) На внешней границе слоя раствора, , будем считать концентрацию ионов постоянной:

(16)

Для скорости течения раствора будем использовать условие исчезновения нормальной скорости и вязкого напряжения на внешнем «свободном» краю слоя раствора:

(17)

3) На входе в рассматриваемую область будем считать, что концентрации, потенциал, скорость течения раствора распределенными линейно, причем их распределение соответствует предельной плотности тока и при выполняется условие электронейтральности ():

, (18)

Как следует из (17) и (18) граничные условия в точке не согласованы. Однако при численном решении это достаточно быстро сглаживается.

4) На выходе из рассматриваемой области будем использовать «мягкие» условия на концентрации и потенциал:

, (19)

5) Начальные условия при примем согласованными с остальными граничными условиями:

, (20)

Нами были проведены численные эксперименты для раствора NaCl в широком спектре таких параметров, как начальная концентрация, скорость прокачивания раствора, межмембранное расстояние, длина, скачок электрического потенциала и определены основные закономерности распределения электрохимических (концентрации, напряженности электрического поля, и т.д.) и гидродинамических полей. Ниже представлены некоторые результаты численных экспериментов при следующих входных параметрах: ширина слоя раствора мм, длина мм, начальная скорость прокачивания раствора м/с или м/с, начальная концентрация раствора моль/м³, температура раствора K, начальная плотность раствора кг/м³, коэффициент кинематической вязкости м²/с, коэффициент диффузии катиона и аниона, соответственно, м²/с, м²/с, разность потенциалов электрического поля в канале В.

3. Основные закономерности процесса электроконвекции

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

3.1. Изменение полей концентраций, потенциала и скорости.

Рисунок 2 - Поля концентраций катионов, рассчитанные спустя 100 секунд с момента включения электрического тока при скорости прокачивания раствора м/с и скачке потенциала В. а) гладкий канал, б) канал с кавернами и в) со вставками без учета электроконвекции и г), д), е) с учетом электроконвекции, соответственно.

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

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

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

Рисунок 3 - Линии тока раствора электролита, рассчитанные спустя 100 секунд с момента включения электрического тока при скорости прокачивания раствора м/с и скачке потенциала В а) гладкий канал, б) канал с кавернами и в) со вставками с учетом электроконвекции, соответственно.

Рисунок 4 - Линии тока раствора электролита, рассчитанные спустя 100 секунд с момента включения электрического тока при скорости прокачивания раствора м/с и скачке потенциала В. а) гладкий канал, б)канал с кавернами и в) со вставками с учетом электроконвекции, соответственно.

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

3.2. Вольтамперная характеристика

Сравнение вольтамперных характеристик (рис. 5) показывает, что электроконвекция во всех трех случаях приводит к росту тока над предельным значением.

Рисунок 5 - Вольтамперные характеристики для области с гладкой мембраной, с мембраной с кавернами и выступами с учетом и без электроконвекции при скорости прокачивания раствора м/с

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

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

Из проведенного выше исследования можно сделать вывод, что при проведении процесса обессоливания в ЭДА, при больших падениях потенциала, модификация поверхности мембраны, особенно с образованием каверн, перспективно для интенсификации процесса, вследствие возникновения и развития электроконвекции.

ПРИМЕЧАНИЯ

1. Nikonenko V., Pismenskaya N., Belova E., Sistat Ph., Larchet Ch., Pourcelly G. // Adv. Colloid Interface Sci. 2010. V. 160. P. 101.

2. Rubinshtein I. Theory of concentration polarization effects in electrodialysis on counter-ion selectivity of ion-exchange membranes with differing counter-ion distribution coefficients // J. Chem. Soc. Faraday Trans. 1990. Vol. 86. №10.P.1857.

3. Rubinshtein I., Maletzki F. Elektroconvection at an electrically inhomoheneous permselective membrane surface //J.Chem. Soc., Faradey Trans.II. 1991. V.87. №13. P.2079-2087.

4. Sata, T. Studies on anion exchange membranes having permselectivity for specific anions in electrodialysis - effect of anion exchange membranes on permselectivity of anions // J. Membr. Sci. - 2000.-Vol.167.- №1.- P 1-31.

5. Лопаткова Г.Ю. Влияние свойств поверхности ионообменных мембран на их электрохимическое поведение в сверхпредельных токовых режимах: Дис. … канд. хим. наук. Краснодар. - 2006. -180 с.

6. Ньюман Дж. Электрохимические системы. М.: Мир, 1977. 463 с.

7. Роуч П. Вычислительная гидродинамика. М.: Мир, 1980. 616 с.

8. Узденова А.М., Коваленко А.В., Уртенов М.Х. Моделирование электроконвекции в электромембранной системе при наличии вынужденной конвекции // Экологический вестник, 2011, №1. С. 77 - 83.

Размещено на Allbest.ru