Математическое моделирование электроконвекции в капилляре. Переходный режим

Разработка математической модели переноса ионов бинарной соли при электроосмотическом течении в капилляре. Исследование краевой задачи для связанной системы уравнений Нернста, Планка, Пуассона и Навье-Стокса. Основные закономерности переноса ионов соли.

Рубрика Физика и энергетика
Вид статья
Язык русский
Дата добавления 20.05.2017
Размер файла 3,4 M

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

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

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

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

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

РИНЦ SPIN-код автора: 3693-4813

Scopus Author ID: 55328224000

Сокова Анна Александровна

РИНЦ SPIN-код автора: 3258-0297

Уртенов Махамет Али Хусеевич, д.ф.-м.н., профессор

РИНЦ SPIN-код: 7189-0748

Scopus Author ID: 6603363090

Кубанский государственный университет,

Краснодар, Россия

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

Ключевые слова: КАПИЛЛЯР, ЭЛЕКТРООСМОС, ЭЛЕКТРОКОНВЕКЦИЯ, УРАВНЕНИЕ ПУАССОНА, УРАВНЕНИЯ НАВЬЕ-СТОКСА, УРАВНЕНИЯ НЕРНСТА-ПЛАНКА

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

В работах [3, 4, 5, 6] электроосмотическое течение исследуется с помощью уравнений, описывающих концентрацию ионов, Нернста-Планка, и уравнения Пуассона-Больцмана, которое описывает распределение потенциала двойного электрического слоя (ДЭС) при равновесных условиях [2]. Но, как отмечается в [2], уравнение Пуассона-Больцмана не может быть использовано для описания перекрывающегося ДЭС в очень узких наноканалах, так как предположение электронейтральности не выполняется. В связи с этим многие исследователи для описания явлений, связанных с электрокинетическим переносом, в том числе и электроконвекции, рассматривают систему уравнений Нернста-Планка и уравнение Пуассона для потенциала [1, 2, 7, 8]. Механизмы переноса ионов электроосмотическим течением в наноканале рассматривались в работах [1, 7]. Было установлено, что возле границы мембрана/электрод, к области ДЭС, примыкает область с обедненным раствором и в этой области наблюдаются электроконвективные вихри. В работе [2] исследуется сверхпредельный перенос в узком и длинном канале с параллельными стенками, один конец канала погружен в резервуар, другой конец заделан катионообменной мембраной. Считается, что стенки канала заряжены. При расчетах используются усредненные уравнения Нернста-Планка. Установлено, что электроконвекция вносит основной вклад в сверхпредельный перенос в достаточно широких каналах. Предложенная модель описывает сверхпредельный перенос посредством электроконвекции.

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

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

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

Рисунок 1. Схема капилляра

Уравнения. Массоперенос с учетом электроконвекции в электромембранных системах описывается электродиффузионными уравнениями Нернста-Планка-Пуассона и уравнениями Навье-Стокса, с учетом пространственной силы. Векторная запись этой системы для бинарного электролита, в случае отсутствия химических реакций, имеет вид:

(1)

(2)

(3)

(4)

(5)

(6)

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

Уравнения Нернста-Планка (1) описывают поток растворенных компонентов, обусловленный миграцией в электрическом поле, диффузией и конвекцией; (2) - уравнение материального баланса; (3) - уравнение Пуассона для потенциала электрического поля; (4) - плотность тока в растворе электролита, обусловлена движением заряженных компонентов; уравнения Навье-Стокса (5) и неразрывности для несжимаемой жидкости (6) описывают поле скоростей, формируемое под действием пространственной электрической силы .

Замечание 1. Если подставить (1) в (2), то уравнения (2) запишутся в виде: .

Замечание 2. Плотность силы электрического поля имеет вид:

, (7)

где плотность распределения пространственного заряда, а - напряженность электрического поля. Формулу (7) можно записать и по-другому, с использованием уравнения Пуассона (3):

.

Краевые условия.

1) Граничные условия на входе ():

, (8)

(9)

(10)

Условие (10) означает открытую границу, т.е. раствор может, как входить, так и выходить через нее.

2) Граничные условия на боковых стенках (, ):

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

, (11)

. (12)

3) Граничные условия на мембране ():

для скорости используется условие прилипания, для анионов используем условия непроницаемости, концентрацию катионов считаем равной емкости мембраны:

(13)

(14)

4) Начальные условия:

.

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

Основные закономерности течения раствора. Ниже на рис. 2 приведены линии тока раствора.

..

a) b)

c) d)

Рисунок 2. Линии тока жидкости и направления течения при a) t=10c, b) t=60c, c) t=100c, d) t=200 c. Цветом выделена концентрация анионов.

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

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

a) b)

Рисунок 3. Поверхность концентрации катионов при

a) t=10c, b) t=200 c.

a) b)

Рисунок 4. Поверхность концентрации анионов при

a) t=10c, b) t=200 c.

..

a) b)

Рисунок 5. Распределение концентрации a) катионов и b) анионов вблизи границы раствор/ катионообменная мембрана.

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

..

a) b)

Рисунок 6. Распределение концентрации а) катионов и b) анионов в вихревой области. Вихревые области выделены серыми прямоугольникоми.

a) b)

Рисунок 7. Распределение концентрации а) катионов и b) анионов в областях обессоливания. Области обессоливания выделены серыми прямоугольниками.

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

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

Работа выполнена при финансовой поддержке грантов РФФИ 13-08-96525 р_юг_а и РФФИ 13-08-96519 р_юг_а.

Литература

1. Fundamentals and modeling of electrokinetic transport in nanochannels// Hung-Chun Yeh, Moran Wang, Chih-Chang Chang, and Ruey-Jen Yang. Israel Journal of Chemistry. 2014. V. 54. pp. 1533 - 1555.

2. Overlimiting Current in a Microchannel// E. Victoria Dydek, Boris Zaltzman, Isaak Rubinstein, D. S. Deng, Ali Mani and Martin Z. Bazant. Physical Review Letters. 2011. V. 107. pp. 118301 - 118305.

3. Poisson-Boltzmann-Nernst-Planck model// Qiong Zheng, Guo-Wei Wei. The Journal of Chemical Physics. 2011. V. 134. pp. 194101-191017.

4. Numerical simulation of electroosmotic flow in microchannels with sinusoidal roughness// Dayong Yang, Ying Liu. Colloids and Surfaces A: Physicochemical and Engineering Aspects. 2008. pp. 28-33.

5. OnApplicability of Poisson-Boltzmann Equation for Micro- and Nanoscale Electroosmotic Flows// Moran Wang, Shiyi Chen. Communications and computational physics. 2008. pp. 1087-1099.

6. Thermal effects on mixed electro-osmotic and pressure-driven flows in triangle microchannels// Q. Liao, X. Zhu, T.Y. Wen. Applied Thermal Engineering. 2009. V. 29. pp. 807-814.

7. Numerical study on transient induced-charge electro-osmotic flow in a cavity// Y.K. Suh. - Colloids and Surfaces A: Physicochemical and Engineering Aspects. 2011. pp. 111-121.

8. Numerical analysis on electroosmotic flows in a microchannel with rectangle-waved surface roughness using the Poisson-Nernst-Planck model// Sangmo Kang Yong Kweon Suh. Microfluid Nanofluid. 2009. pp. 461-477.

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

...

Подобные документы

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

    курсовая работа [2,3 M], добавлен 30.08.2012

  • Основные свойства жидкости. Отсутствие идеальной модели и трудности формулировки общей теории жидкости. Явления переноса: диффузия, теплопроводность и вязкость, их характеристика. Отличия явлений переноса в жидкостях от аналогичных явлений в газах.

    реферат [40,2 K], добавлен 05.06.2009

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

    контрольная работа [1,0 M], добавлен 24.02.2013

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

    дипломная работа [3,4 M], добавлен 26.11.2011

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

    статья [105,9 K], добавлен 30.09.2012

  • Вязкость, движение частиц в вязких средах. Электропроводность и ее виды. Удельная и молярная электропроводность растворов электролитов. Числа переноса и методы их определения. Проверка концентрации кислоты методом потенциометрического титрования.

    курсовая работа [743,5 K], добавлен 17.12.2014

  • Групповая скорость. Парадокс. Вектор Пойнтинга. Проблемы определения скорости переноса энергии. Скорость переноса энергии ТЕ и ТМ волн. Фазовая скорость это скорость движения силового свойства поля.

    реферат [95,4 K], добавлен 02.03.2002

  • Исследование свойств теплопроводности как физического процесса переноса тепловой энергии структурными частицами вещества в процесс их теплового движения. Общая характеристика основных видов переноса тепла. Расчет теплопроводности через плоскую стенку.

    реферат [19,8 K], добавлен 24.01.2012

  • Продольное удержание плазмы в Газодинамической ловушке, поперечные потери, удержание быстрых ионов и микронеустойчивости. Диагностики: двухсеточный зонд, пироэлектрический болометр, 45 анализатор энергий ионов. Результаты измерений и их интерпретация.

    дипломная работа [2,5 M], добавлен 19.02.2013

  • Конвективный теплообмен - распространение тепла в жидкости (газе) от поверхности твердого тела или к ней. Смысл закона Ньютона, дифференциального уравнения Фурье - Кирхгофа и критериального уравнения Навье – Стокса. Теплоотдача при конденсации паров.

    реферат [208,1 K], добавлен 15.10.2011

  • Типы электрохимических цепей и электродов. Сущность метода потенциометрии. Определение растворимости малорастворимой соли на примере хлорида серебра с использованием концентрационной цепи с переносом. Нормальный элемент Вестона, специфика его устройства.

    курсовая работа [3,0 M], добавлен 06.04.2015

  • Уравнение теплового баланса. Теплота, подведенная теплопроводностью и конвекцией, к элементарному объему. Общий вид дифференциального уравнения энергии Фурье-Кирхгофа. Применение ряда Тейлора. Дифференциальное уравнение движения жидкости Навье-Стокса.

    презентация [197,5 K], добавлен 18.10.2013

  • Понятие процесса переноса тепла и вещества, потенциалы переноса. Температурное поле, примеры одномерного и двухмерного полей. Стационарный и нестационарный процесс теплопередачи. Характеристика параметров материала: плотность, пористость, влажность.

    контрольная работа [203,4 K], добавлен 21.01.2012

  • Определение плотности потока формамида через плазматическую мембрану Chara ceratophylla толщиной 8 нм. Расчет коэффициента проницаемости плазматической мембраны Mycoplasma для формамида. Уравнение Фика для диффузии. Расчет потенциала Нернста для ионов.

    контрольная работа [286,8 K], добавлен 09.01.2015

  • Расчет энергии иона. Количественная интерпретация данных о рассеянии быстрых ионов. Метод спектроскопии обратно рассеянных ионов низких энергий. Форма энергетических спектров двухкомпонентных материалов. Спектр кремния с анатомами на поверхности.

    контрольная работа [86,3 K], добавлен 14.11.2011

  • Явления переноса в газах. Число столкновений и средняя длина свободного пробега молекул в газах. Диффузия газов и внутреннее трение. Вязкость и теплопроводность газов. Коэффициенты переноса и их зависимость от давления. Понятие о вакуумном состоянии.

    презентация [2,7 M], добавлен 13.02.2016

  • Потенциал действия и его фазы. Роль ионов Na K в генерации потенциала действия в нервных и мышечных волокнах: роль ионов Ca и Cl. Восстановление от радиационного поражения. Основные методы регистрации радиоактивных излучений и частиц. Их характеристика.

    контрольная работа [17,3 K], добавлен 08.01.2011

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

    курсовая работа [665,8 K], добавлен 29.11.2012

  • Ионная природа мембранных потенциалов. Потенциал покоя, уравнение Нернста. Стационарный потенциал Гольдмана-Ходжкина. Уравнение электродиффузии ионов через мембрану в приближении однородного поля. Механизм генерации и распространения потенциала действия.

    реферат [158,6 K], добавлен 16.12.2015

  • Силы и коэффициент внутреннего трения жидкости, использование формулы Ньютона. Описание динамики с помощью формулы Пуазейля. Уравнение Эйлера - одно из основных уравнений гидродинамики идеальной жидкости. Течение вязкой жидкости. Уравнение Навье-Стокса.

    курсовая работа [531,8 K], добавлен 24.12.2013

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.