Математическое моделирование и численное исследование гидродинамики в экспериментальной электрохимической ячейке с вращающимся мембранным диском

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

Математическое моделирование и численное исследование гидродинамики в экспериментальной электрохимической ячейке с вращающимся мембранным диском

Введение

электрохимическая ячейка гидродинамика

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

Рисунок 1. Гидродинамическая и электрическая схемы экспериментальной ячейки с ВМД [4]: 1- верхняя полуячейка с раствором (катодная камера); 2 - мембрана МК-40; 3 - нижняя полуячейка с раствором (анодная камера); 4- капилляр для подвода раствора;5- капилляр для отвода раствора; 6 - поляризующие электроды; 7 - капилляры Луггина-Габера; 8 - галваностат П5848; 9 - мультивольтметр (иономер И-130); 10 - электроды сравнения ; 11 - шкив.

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

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

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

электрохимическая ячейка гидродинамика

В настоящей работе рассматривается задача описания гидродинамики при вращении трубки внутри цилиндрической ячейки вокруг центральной оси (рисунок 2).

Рисунок 2. Схематическое сечение ячейки с ВМД.

Указаны реальные размеры ячейки в мм.

1.1 Область

При численном решении используется осевая симметрия задачи, поэтому описывается половина сечения цилиндра, где и определяются уравнения и граничные условия (см. рис.3).

Рисунок 3. Исследуемая область и ее границы: 1 - ось симметрии, 3 - катионообменная мембрана, 2, 4, 6, 7, 8 - стенки ячейки, 5 - свободная поверхность для открытой ячейки или стенка для закрытой ячейки

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

1.2 Уравнение

Для моделирования течения жидкости используются уравнения Навье-Стокса. В уравнении Навье-Стокса (1) u обозначает скорость (м\с), с ? плотность (кг/м³), з ? динамическую вязкость (Pa·s), и p ? давление (Pa).

(1)

Для постоянного 3D осесимметричного потока размерность уравнений (1) можно уменьшить до 2, где u - радиальная скорость, v - вращательная скорость, и w - осевая скорость (м\с). Модель включает все три скоростных компонента, однако циркулирующий поток является уже 2D.

(2)

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

1.3 Граничные условия

Опишем граничные условия на каждой из границ.

На границе №1:

На границе, соответствующей оси вращения, используется условие симметрии в тангенциальном направлении границы (не в нормальном направлении). Радиальная скорость равна нулю: .

На границах №3-4:

Границы 3, 4 считаются «скользящими стенами» и используется граничное условие (щ - угловая скорость, r - радиус, -линейная скорость оси симметрии, - линейная скорость произвольной точки на границах 3, 4): , где , .

На границах №2, 6, 7, 8:

На данных границах задается условие прилипания, то есть все скоростные компоненты равны нулю: u=(0,0,0).

На границе №5:

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

.

1.4 Начальные значения и свойства жидкости

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

. Плотность и кинематическая вязкость имеют следующие значения: .

1.5 Метод численного решения

Для решения использовался метод конечных элементов с неравномерной сеткой, причем количество элементов 20623.

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

Течение раствора в ячейке зависит от того, является она закрытой или открытой.

Закрытая ячейка

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

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

а) б)

Рисунок 4. Течение раствора при угловой (60 оборотов в минуту) в момент времени : а)-линии тока раствора, б)-направления движения(нормализованные)

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

Открытая ячейка

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

а) б)

в) г)

Рисунок 5. Течение раствора при угловой (120 оборотов в минуту) в момент времени : а)-линии тока раствора, б)-направления движения(нормализованные) и : в)-линии тока раствора, г)-направления движения(нормализованные).

а) б)

в) г)

Рисунок 6. Линии тока раствора раствора при угловой (60 оборотов в минуту) в момент времени: а), б), в) , г)

Выводы

В работе теоретически, на основе математического моделирования и численного решения, исследованы закономерности гидродинамики в электрохимической ячейке с ВМД с учетом её реальных геометрических размеров и конструкционных особенностей и показано, что:

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

2 Несмотря на то, что характер течения во всей ячейке достаточно сложный, под ВМД имеется область, где линии тока раствора близки к логарифмическим спиралям. Граница этой области для открытой ячейки колеблется.

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

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

Литература

электрохимическая ячейка гидродинамика

1. Левич В.Г. Физико-химическая гидродинамика. М.: Физматгиз, 1959, 700с.

2. Исаев Н.И., Золотарева Р.И., Иванов Э.М. // Журн. физ. химии, 1967. Т.41.С.849

3. Завгородных Л.А., Бобрешова О.В., Кулинцов П.И., Аристов И.В. //Электрохимия.2006. Т.42.С.68.

4. Заболоцкий В.И., Шельдешов Н.В., Шарафан М.В.//Электрохимия. 2006, т.42., №11, С.1-7.

5. Патент №78577 РФ, МПК G01N 27/40, 27/333 «Устройство для одновременного измерения вольтамперных характеристик и чисел переноса ионов в электромембранных системах». Шарафан М.В., Заболоцкий В.И. №2008122083/22 от 02.06.2008, опубл. 27.11.2008, Бюл. № 33.

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