Математическое моделирование и численное исследование гидродинамики в экспериментальной электрохимической ячейке с вращающимся мембранным диском
Исследование гидродинамики экспериментальной электрохимической ячейки с вращающимся мембранным диском. Анализ течения раствора в открытых и герметично закрытых ячейках. Закономерности гидродинамики с учетом реальных размеров экспериментальной ячейки.
Рубрика | Производство и технологии |
Вид | статья |
Язык | русский |
Дата добавления | 13.05.2017 |
Размер файла | 833,5 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на 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 обозначает скорость (м\с), с ? плотность (кг/м3), з ? динамическую вязкость (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
Подобные документы
Описание режимов работы ситчатой и колпачковой тарелок ректификационной колонны. Экспериментальное определение гидравлического сопротивления сухой и орошаемой тарелки. Расчет гидродинамики тарельчатых колонн и сравнение с экспериментальным результатом.
лабораторная работа [265,5 K], добавлен 15.12.2014Исследования процессов взаимодействия образцов конструкционных материалов ЯЭУ с жидкометаллическими теплоносителями. Моделирование взаимодействия реакторных сталей на установке ЭУ "ВД". Использование метода вращающегося диска для натриевого теплоносителя.
дипломная работа [1,2 M], добавлен 03.01.2014Проект установки для изучения течения и процессов теплоотдачи в сложных пространственных каналах. Определение расчётных параметров течения в экспериментальной установке на четырёх участках. Разработка методики определения расхода воздуха по его нагреву.
дипломная работа [1,8 M], добавлен 06.06.2013Индицирование рентгенограммы поликристаллов, анализ параметров и типа решетки для кубической сингонии, средней категории. Ячейки Бравэ и систематические погасания. Оформление экспериментальной рентгенограммы. Программы для автоматического индицирования.
курсовая работа [90,0 K], добавлен 15.12.2015Исследование процесса проектирования подсистемы моделирования работы гибкой производственной ячейки и графического представления результатов на экране. Анализ формирования параметров оборудования, путем сопоставления с необходимым коэффициентом загрузки.
курсовая работа [1,3 M], добавлен 09.04.2012Сравнение фонтанирования и псевдоожижения. Разработка метода расчета коэффициента гидравлического сопротивления топочного устройства и технологической схемы экспериментальной установки. Изучение движения газовзвеси в экспериментальной топочной камере.
курсовая работа [900,1 K], добавлен 31.07.2015Технико-экономическое обоснование и выбор типа установки электрохимической защиты газопровода. Расчет установки электрохимической защиты, эксплуатация протекторной станции. Техника безопасности и мероприятия по охране окружающей среды при эксплуатации.
курсовая работа [750,2 K], добавлен 07.03.2012Расчет параметров электрохимической обработки детали. Изучение процессов на поверхности твердого тела при вакуумном ионно-плазменном напылении порошка борида циркония. Анализ показателей температурных полей при наплавке покрытия плазменно-дуговым методом.
курсовая работа [2,5 M], добавлен 06.12.2013Численное исследование силового взаимодействия газовой струи и несжимаемой жидкости через контактную поверхность. Физико-математическое моделирование кислородно-конвертерного процесса. Влияние управляющих параметров (давления и температуры в газопроводе).
дипломная работа [2,5 M], добавлен 18.02.2011Принципы построения комбинированной гидродинамической модели аппарата методом декомпозиции функции отклика системы на возмущение идентификацией простейших типовых гидродинамических моделей. Разработка химического реактора с учетом его гидродинамики.
контрольная работа [304,4 K], добавлен 02.12.2015