Математическое моделирование электроконвекции в мембранных системах
На правах рукописи
Узденова Аминат Магометовна
Математическое моделирование электроконвекции
в мембранных системах
05.13.18 – математическое моделирование, численные методы
и комплексы программ
АВТОРЕФЕРАТ
диссертации на соискание ученой степени
кандидата физико-математических наук
Краснодар – 2012
Работа выполнена в ФГБОУ ВПО «Карачаево-Черкесский государственный
университет им. У.Д. Алиева»
Научный руководитель: доктор физико-математических наук, профессор
Уртенов Махамет Али Хусеевич
Официальные оппоненты: доктор физико-математических наук,
Демехин Евгений Афанасьевич,
кандидат физико-математических наук,
Письменский Александр Владимирович
Ведущая организация: ФГБОУ ВПО «Российский государственный университет нефти и газа имени И.М. Губкина» г. Москва
Защита состоится «19» марта 2012 г. в 16 ч. 30 мин. на заседании диссертационного совета Д 212.101.17 в ФГБОУ ВПО «Кубанский государственный университет» по адресу: 350040, г. Краснодар, ул. Ставропольская, д. 149, в ауд. 231.
С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Кубанский государственный университет», с авторефератом – на сайте http://www.kubsu.ru/
Автореферат разослан «16» февраля 2012 г.
Ученый секретарь
диссертационного совета Д 212.101.17
кандидат физ.-мат. наук, доцент В.Ю. Барсукова
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ
Актуальность темы. Математическое моделирование электроконвекции привлекало внимание многих ученых. В работах по электроконвекции в мембранных системах выделяют два её типа: объемную электроконвекцию, обусловленную действием электрического поля на пространственный заряд в квазиэлектронейтральном растворе, и электроосмос, связанный с существованием на границе фаз двойного электрического слоя. В настоящее время доказано отсутствие объемной электроконвекции в реальных условиях. В зависимости от структуры двойного электрического слоя для электроосмотического скольжения выделяют два фундаментально различных режима: электроосмос первого и второго родов.
Теория электроосмоса первого и второго родов была разработана С.С. Духиным, Б.В. Дерягиным, Н.А. Мищук, Э. Жолковским, Э. Стодом, М. Воротынцевым др. В работах И. Рубинштейна, Б. Зальцмана, И. Лерман и др. были проведены теоретические и экспериментальные исследования электроконвекции. В указанных выше работах заложены основы математического моделирования и теории электроконвекции в мембранных системах. Однако в них ограничиваются рассмотрением сравнительно простых математических моделей при различных упрощающих предположениях:
1) предположение об электронейтральности раствора,
2) стационарности электроконвекции,
3) принятие концентрации постоянной на границе с мембраной,
4) рассматривалась лишь непроточная ячейка и др.
Таким образом, тему диссертации, посвященную построению и исследованию математических моделей электроконвекции в мембранных системах, свободных от указанных выше ограничений, разработке эффективных алгоритмов численного и асимптотического анализа, установлению основных закономерностей электроконвекции, следует признать актуальной.
Актуальность темы исследования подтверждается также поддержкой, оказанной работе Федеральным Агентством по образованию и науке РФ в рамках темы 1.4.08 («Методы регулярного представления сингулярно возмущенных уравнений и их приложения. Метод модулирующих функций в обратной задаче теории фильтрации» (направление фундаментальных научных исследований «Рациональное природопользование») и гранта РФФИ-Юг (№ 09-08-96529 «Модифицирование поверхности ионообменных мембран с использованием углеродных нанотрубок с целью совершенствования процессов электродиализного обессоливания и концентрирования»).
Объектом исследования являются закономерности электроконвекции в мембранных системах.
Предметом исследования являются математические модели электроконвекции в мембранных системах.
Цель исследования. Разработка математических моделей электроконвекции в мембранных системах, асимптотических и численных методов их исследования.
Цель исследования предопределила следующие задачи исследования:
- Построение математических моделей электроконвекции в мембранных системах.
- Построение эффективных численных и асимптотических методов решения.
- Разработка проблемно-ориентированных программ для проведения вычислительных экспериментов.
- Установление основных закономерностей электроконвекции в мембранных системах.
Научная новизна:
- Разработана иерархическая система новых моделей, адекватно описывающих электроконвекцию в мембранной системе в потенциостатическом и гальваностатическом режимах: общая модель, модель в приближении соленоидальности поля плотности тока, модель в приближении закона Ома. Эти модели разработаны, как для гладких мембран, так и мембран с выступами и кавернами.
- Предложены новые алгоритмы и методы численного и асимптотического решения краевой задачи, отвечающей модели в приближении закона Ома, особенностью которых является оригинальное сочетание численных и асимптотических методов с последовательным уточнением области применимости асимптотических решений. Даны асимптотические оценки электрического и концентрационного полей.
- Впервые установлены основные закономерности электроконвекции в проточных мембранных системах: структура и динамика электрического и концентрационного полей, динамика возникновения и развития электроконвективных вихрей, влияние геометрической неоднородности поверхности мембраны на электроконвекцию.
Научная и практическая значимость
- Предложенные в диссертации уравнение для плотности электрического тока, формула для напряженности электрического поля, модели в приближении соленоидальности плотности тока и закона Ома могут быть использованы для разработки математических моделей в других задачах, описывающих перенос в мембранных, нано- и микрофлюидных системах.
- Методы асимптотического и численного решения краевых задач, предложенные в диссертации, могут быть применены при решении краевых задач для системы квазилинейных уравнений математической физики.
- Установленные в диссертации основные закономерности электроконвекции могут быть использованы научно-исследовательскими группами, проектными организациями для повышения эффективности электродиализных аппаратов водоподготовки, при разработке новых конструкций этих аппаратов и нано- и микрофлюидных устройств.
Основные положения, выносимые на защиту:
- Иерархическая система математических моделей электроконвекции в мембранных системах для потенциостатического и гальваностатического режимов, выведенная на основе предложенного в диссертации метода моделирования процессов переноса в мембранных системах.
- Алгоритмы численного и асимптотического решения краевых задач, соответствующих математическим моделям. Особенностью предлагаемых методов является оригинальное сочетание численных и асимптотических методов, причем в ходе решения последовательно уточняются области применимости асимптотических решений. Приближенные аналитические формулы для расчета электрохимических полей.
- Новое уравнение для плотности электрического тока, устанавливающее соответствие между плотностью тока, напряженностью электрического поля и концентрацией электролита. Формула для напряженности электрического поля, устанавливающая зависимость напряженности электрического поля от плотности тока, градиента концентрацией электролита и скорости конвективного потока электролита.
- Основные закономерности электроконвекции в мембранных системах, а именно структура распределения электрохимических и гидродинамических полей по ширине и длине канала обессоливания, влияние на перенос ионов соли электроконвекции и геометрической неоднородности поверхности мембран при интенсивных токовых режимах, критическая (пороговая) кривая возникновения электроконвекции в координатах «скорость вынужденного течения раствора/падение потенциала», установленные с помощью комплексного исследования с использованием математических моделей и вычислительного эксперимента.
- Комплекс проблемно-ориентированных программ для проведения вычислительных экспериментов по исследованию электроконвекции в канале обессоливания электродиализного аппарата с учетом геометрической неоднородности поверхности мембран и проведение вычислительных экспериментов.
Внедрение. Результаты диссертационного исследования использованы в работе инновационного технологического Центра «Кубань-Юг» при проектировании новых систем водоподготовки, в учебном процессе ФГБОУ ВПО «Карачаево-Черкесский государственный университет им. У.Д. Алиева».
Обоснованность и достоверность результатов, полученных в диссертации, обеспечивается выбором адекватных моделей, применением строгих математических методов, сравнением с результатами других авторов и с литературными экспериментальными данными, в тех случаях, когда это возможно.
Личный вклад автора. Основные результаты работы получены лично автором: построена иерархическая система математических моделей процесса электроконвекции; построены алгоритмы численного и асимптотического решения краевых задач соответствующих математических моделей; разработан комплекс проблемно-ориентированных программ для проведения вычислительных экспериментов по моделированию массопереноса в канале обессоливания электродиализного аппарата с учетом электроконвекции и геометрической неоднородности поверхности мембран; выявлены основные закономерности электроконвекции в канале обессоливания электродиализного аппарата для бинарного электролита.
Апробация работы. Результаты диссертации были доложены:
- на 4-х Международных конференциях: «Ion transport in organic and inorganic membranes» (г. Туапсе, 2008 – 2011 гг.) и 2-х Всероссийских конференциях: «VI-VII Всероссийские конференции «Современное состояние и приоритеты развития фундаментальных наук в регионах»» (г. Анапа, 2009, 2010 гг.);
- на научных семинарах кафедры математического анализа КЧГУ (2008, 2009, 2010 гг.), прикладной математики КубГУ (2009, 2010 г.), кафедры физической химии КубГУ (2011 г.);
- на научных конференциях студентов и аспирантов КЧГУ (г. Карачаевск, 2009, 2010, 2011 гг.) и КубГУ (г. Краснодар, 2009, 2010 гг.).
Публикации. По материалам диссертации опубликовано 18 печатных работ, из них 1 монография, 8 статей, 9 тезисов докладов, в том числе 5 статей в журналах, рекомендованных ВАК для опубликования основных результатов докторских и кандидатских диссертаций.
Структура и объём диссертации. Диссертация состоит из введения, четырех глав, заключения, списка литературы (164 наим.), двух приложений. Работа изложена на 144 стр., в том числе содержит 33 рисунка, 1 таблицу.
СОДЕРЖАНИЕ РАБОТЫ
Во введении обоснована актуальность диссертационной работы, сформулированы цель и задачи, перечислены результаты, выносимые на защиту, сформулированы научная новизна и практическая ценность исследования, определен личный вклад автора, указано содержание работы по разделам.
В первой главе приводятся краткие сведения по электромембранным системам и процессам, необходимые для понимания и содержательной интерпретации математических моделей массопереноса в этих системах. В этой же главе представлены основные законы, описывающие массоперенос в мембранных системах с учетом электроконвекции. Анализируются известные математические модели электроконвекции С.С. Духина, Б.В. Дерягина, И. Рубинштейна, Б. Зальцмана и др., составляющие основу теории электроконвекции в мембранных системах.
Пусть и – ширина и длина камеры обессоливания, соответственно, - средняя скорость прокачивания раствора, - соответствует условной межфазной границе катионообменная мембрана/раствор, - соответствует условной межфазной границе анионообменная мембрана/раствор, - входу, а - выходу из камеры обессоливания, – скорость течения раствора, – характерная плотность раствора, – давление, – концентрации катионов и анионов в растворе, соответственно, – зарядовые числа катионов и анионов, – коэффициенты диффузии катионов и анионов, соответственно, – потенциал электрического поля, - напряженность электрического поля, – диэлектрическая проницаемость электролита, – постоянная Фарадея, – газовая постоянная, – абсолютная температура, – время, – коэффициент кинематической вязкости, – плотность силы электрического поля. В данной задаче – неизвестные функции, зависящие от времени и координат , .
При моделировании электроконвекции в главах 2-4 учитывается тот факт, что мембранные системы используют, как правило, в двух разных режимах работы: потенциостатическом режиме, которому соответствует условие:
, (1)
означающее, что величина падения потенциала в камере обессоливания постоянна, и гальваностатическом режиме, когда ток, протекающий через любое сечение камеры обессоливания, является постоянным, т.е.:
, (2)
где - средний ток.
К системе уравнений, которые будут ниже использоваться для разработки различных математических моделей электроконвекции, должны быть добавлены соответствующие краевые условия, причем они должны быть согласованы со свойствами мембран и величиной или («интенсивным» или «мягким» токовым режимом).
Наряду с одним из указанных выше условий на потенциал или плотность тока, используются и другие условия для потенциала (плотности тока), для скорости и следующие граничные условия для концентраций:
1). На поверхности катионообменной мембраны будем считать граничную концентрацию катионов равной фиксированному заряду внутри мембраны:
(3)
Кроме того, предположим катионообменную мембрану идеально селективной, т.е. непроницаемой для анионов:
(4)
2). На поверхности анионообменной мембраны будем считать граничную концентрацию анионов равной фиксированному заряду внутри мембраны:
(5)
Предположим анионообменную мембрану идеально селективной, то есть непроницаемой для катионов:
(6)
3). На входе в рассматриваемую область будем полагать:
(7)
4). На выходе из рассматриваемой области будем использовать «мягкие» условия на концентрации:
(8)
5). Начальные условия при примем согласованными, по возможности, с остальными граничными условиями.
Во второй главе описана разработанная в диссертации иерархическая система математических моделей электроконвекции в гальваностатическом режиме. Для математического моделирования явлений переноса, как правило, используются уравнения Нернста-Планка и Пуассона, совместно с системой уравнений Навье-Стокса. Система уравнений Нернста-Планка и Пуассона удобна для моделирования потенциостатического режима, так как уравнение Пуассона используется для нахождения потенциала. Однако она неудобна для исследования гальваностатического режима, поскольку нет явного уравнения для плотности тока. Отсюда следует, что систему уравнений Нернста-Планка и Пуассона необходимо преобразовать к виду удобному для решения указанной проблемы. Для этого оказалось необходимым решить следующие две задачи:
- Вывести дифференциальное уравнение для плотности тока , которое должно использоваться вместо уравнения Пуассона.
- Вывести формулу, выражающую напряженность электрического поля через плотность тока и концентрацию, которая должна использоваться вместо уравнения плотности тока.
В §2.1 предлагается математическая модель переноса ионов в канале обессоливания электродиализного аппарата в гальваностатическом режиме в виде системы квазилинейных уравнений с частными производными.
В §2.2 предложена новая математическая модель переноса ионов в канале обессоливания электродиализного аппарата в гальваностатическом режиме в приближении соленоидальности плотности тока. Введено новое понятие «функция тока» для плотности тока:
(9)
В трехмерном случае функции соответствует векторный потенциал для плотности тока.
Указанная выше модель описывается следующей системой уравнений:
, , (10)
(11)
. (12)
Здесь - кососимметричное скалярное произведение. Эта система уравнений дополняется соответствующими краевыми условиями.
В §2.3 для симметричного бинарного электролита приведен алгоритм полной декомпозиции системы электродиффузионных уравнений и получена система декомпозиционных уравнений, включающая новое уравнение для функции тока для общей плотности тока.
(13)
(14)
(15)
Здесь - обобщенная суммарная концентрация, - общая плотность тока, , .
В §2.4 предложена новая математическая модель переноса ионов в канале обессоливания электродиализного аппарата в гальваностатическом режиме в декомпозиционных переменных, приведен переход к безразмерному виду и дана оценка параметров характерных для камеры обессоливания
В третьей главе предлагаются математические модели электроконвекции в потенциостатическом режиме. Математическая модель рассматривается сначала для канала обессоливания, исследование, которого позволяет сформулировать более простую математическую модель электроконвекции в диффузионном слое и определить границы ее применимости.
В §3.1 сформулирована математическая модель электроконвекции в камере обессоливания электродиализного аппарата. Модель строится на основе системы уравнений Нернста-Планка и Пуассона (18)-(20). Здесь (16), (17) – и уравнения Навье-Стокса, описывающие движение раствора с учетом пространственной электрической силы.
, (16)
, (17)
, (18)
, (19)
. (20)
Для решения задачи (16)-(20) применяется метод конечных элементов.
Показано, что при заданном падении потенциала электроконвекция проявляется, когда скорость прокачки не превышает некоторое пороговое значение. При этом течение в камере обессоливания электродиализного аппарата несимметрично: для раствора NaCl электроконвекция возникает вначале только у катионообменной мембраны. Показано, что это связано с различием значений коэффициентов диффузии для катионов и анионов. Численно определена пороговая кривая возникновения электроконвекции: чем больше скорость вынужденной конвекции, тем выше пороговая разность потенциалов.
Рис. 1. Зависимость пороговой разности потенциалов от скорости прокачивания раствора при
Нами было проведено численное исследование устойчивости, которое показало систематическое уменьшение колебаний плотности тока с уменьшением величины возмущающей добавки и, поэтому, можно утверждать, что процесс сверхпредельного массопереноса ионов соли устойчив как при отсутствии, так и при наличии электроконвекции.
В §3.2 приводится упрощенная модель электроконвекции в диффузионном слое, прилегающем к катионообменной мембране.
В §3.3 предлагается математическая модель влияния электроконвекции на перенос ионов соли в камере обессоливания электродиализного аппарата при наличии вынужденной конвекции и геометрической неоднородности поверхности мембран. Геометрическая неоднородность поверхности мембраны моделируется кавернами и выступами. Выявлено, что наличие геометрической неоднородности поверхности может приводить как к увеличению, так и к уменьшению массопереноса. Вклад от неоднородности в плотность тока зависит от скорости прокачивания раствора, приложенного скачка потенциалов, размеров неоднородности.
В четвертой главе проводится численное исследование электроконвекции в гальваностатическом режиме. В §4.1 главы дан вывод модельной задачи в «приближении закона Ома» из системы декомпозиционных уравнений.
В §4.2 приведено найденное нами асимптотическое представление решения краевой задачи из §4.1. Выявлена структура электрохимических полей и показано, что канал обессоливания, разбивается на области электронейтральности и пространственного заряда, причем асимптотическое представление решения имеет разный вид в областях электронейтральности и пространственного заряда. В §4.3 выведены формулы для расчета силы электрического поля, даны относительные оценки этой силы в областях электронейтральности и пространственного заряда.
В §4.4 предлагается следующий алгоритм численного решения краевой задачи, соответствующей математической модели электроконвекции в «приближения закона Ома»:
1. Задаются две начальные функции, определяющие области электронейтральности и пространственного заряда, например, для нахождения и для нахождения , причем , . Пусть уже определены приближения, , .
2. Определяются области и .
3. В области пространственного заряда решается уравнение для с соответствующими краевыми условиями. Полученное решение обозначим и рассчитаем .
4. В области пространственного заряда рассчитывается сила по формуле
,
а в области электронейтральности полагается равным нулю, т.е.: .
5. Решая краевую задачу для системы уравнений Навье-Стокса, находим распределение поля скоростей .
6. Решается уравнение для функции
,
с соответствующими краевыми условиями.
Полученное решение обозначаем .
7. Находится следующие приближения для функций, определяющих области электронейтральности и пространственного заряда:
;
.
8. Определяются области и .
- Расчеты продолжаются до тех пор, пока не выполнится условие остановки:
,
где - некоторая заданная точность, а - площадь фигуры.
В §4.5 проведен численный анализ математической модели электроконвекции в «приближении закона Ома». Установлены основные закономерности электроконвекции в гальваностатическом режиме.
Верификация результатов моделирования электроконвекции в мембранных системах проведено сравнением ширины диффузионного слоя с экспериментальными данными (рис. 2), сопоставлением вольтамперной характеристики с результатами Гнусина Н.П. и др. (ELDIAL) (рис. 3), течения с экспериментальными данными С. Рубинштейна и др. Результаты показывают хорошее согласие, таким образом, предложенные нами математические модели электроконвекции, адекватны.
Рис. 2. Сопоставление вычислений ширины диффузионного слоя (1) с экспериментальными данными (2) Шапошника В.А. и др.[1]
Рис. 3. Вольтамперная характеристика при скорости прокачивания раствора м/с, начальной концентрации раствора моль/м3;
ElConv – расчет по модели (16)-(20); расчет по ELDIAL[2]
Теоретически рассчитанная нами вольтамперная кривая (рис. 3) имеет классический вид. А именно, на графике вольтамперной характеристики можно выделить три области: начальный участок (допредельный режим, ток возрастает с ростом разности потенциалов), плато (предельный режим, рост разности потенциалов приводит к незначительному увеличению тока), сверхпредельный режим (ток возрастает с ростом напряжения). Таким образом, можно сделать вывод, что сверхпредельный режим объясняется появлением электроконвекции, поскольку электроконвективные вихри обеспечивают доставку раствора из ядра потока к границе мембрана/раствор, что приводит к дополнительному переносу ионов, т.е. к интенсификации массопереноса.
По итогам комплексного исследования электроконвекции определены следующие основные закономерности электроконвекции в канале обессоливания электродиализного аппарата:
- Проявление электроконвекции зависит от значений коэффициентов диффузии катионов и анионов. Для раствора NaCl электроконвекция возникает в диффузионном слое, прилегающем к катионообменной мембране при истощении концентрации катионов соли при запредельной плотности тока. Причиной электроконвекции в проточных мембранных системах является действие силы электрического поля на пространственный заряд, которые образуется на межфазной границе мембрана/раствор электролита.
- Процесс электроконвекции быстро выходит на стационарный режим. Поскольку протяженность пространственного заряда увеличивается с разбавлением раствора (с уменьшением концентрации соли) или (и) увеличением плотности тока, то вклад электроконвекции в сверхпредельный массоперенос растет.
- С увеличением линейной скорости протока раствора вклад электроконвекции в сверхпредельный массоперенос уменьшается, с увеличением абсолютной величины падения потенциала вклад электроконвекции в сверхпредельный массоперенос увеличивается. При заданной скорости прокачки электроконвекция проявляется, когда разность потенциалов превышает некоторое пороговое значение. При заданном падении потенциала электроконвекция проявляется, когда скорости прокачки не превышает некоторое пороговое значение. Численно определена пороговая кривая возникновения электроконвекции: чем больше скорость вынужденной конвекции, тем выше пороговая разность потенциалов (рис. 1).
- Сверхпредельный режим объясняется появлением электроконвекции, поскольку электроконвективные вихри обеспечивают доставку раствора из ядра потока к границе мембрана/раствор, что приводит к дополнительному переносу ионов, т.е. к интенсификации массопереноса. Во всех вариантах введения возмущения наблюдается уменьшение колебаний плотности тока с уменьшением величины возмущающей добавки и можно утверждать, что процесс сверхпредельного массопереноса ионов соли устойчив несколько выше и под пороговой кривой, как при отсутствии, так и наличии при электроконвекции.
- Наличие геометрической неоднородности поверхности может приводить как к увеличению, так и к уменьшению массопереноса. Вклад от неоднородности в плотность тока зависит от скорости прокачивания раствора, приложенного скачка потенциалов, размеров неоднородности. В случае каверн наблюдается усиление добавочного массопереноса при совмещении электроконвекции и геометрической неоднородности. В то же время, наличие вставок ослабляет добавочный массоперенос, вызванный электроконвекцией.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ
В итоге проведенных в диссертации исследований можно сформулировать следующие основные результаты и предложения.
- Предложена иерархическая система математических моделей процесса электроконвекции в мембранных системах для потенциостатического и гальваностатического режимов: общие модели на основе связанной системы уравнений Нернста-Планка, Пуассона и Навье-Стокса, модели в приближении соленоидальности и закона Ома.
- Разработаны эффективные алгоритмы численного и асимптотического решения краевых задач, соответствующих математических моделей. Выведены приближенные аналитические формулы для расчета электрохимических полей. Выведены новое уравнение для плотности электрического тока и формула для напряженности электрического поля.
- Установлены причина и основные закономерности электроконвекции в мембранных системах в двумерном случае, а именно структура распределения электрохимических и гидродинамических полей по ширине и длине канала обессоливания, влияние на перенос ионов соли электроконвекции и геометрической неоднородности поверхности мембран при интенсивных токовых режимах. Определена критическая (пороговая) кривая электроконвекции в координатах «скорость/падение потенциала». Проведена верификация математических моделей и алгоритмов численного решения путем сравнения с результатами других авторов.
- Создан комплекс проблемно-ориентированных программ для проведения вычислительных экспериментов по моделированию процесса переноса ионов соли в канале обессоливания электродиализного аппарата с учетом электроконвекции и геометрической неоднородности поверхности мембран.
ОСНОВНЫЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ
В СЛЕДУЮЩИХ РАБОТАХ:
Монография:
- Узденова А.М., Коваленко А.В., Уртенов М.Х. Математические модели электроконвекции в электромембранных системах. – Карачаевск.: Изд-во КЧГУ, 2011. - 156 с.
Статьи в журналах, рекомендованных ВАК для опубликования результатов кандидатских и докторских диссертаций:
- Узденова А.М., Коваленко А.В., Уртенов М.Х. Моделирование электроконвекции в электромембранной системе при наличии вынужденной конвекции // Экологический вестник научных центров Черноморского экономического сотрудничества. – Краснодар. – 2011. - №1. – С. 77 – 83.
- Узденова А.М., Коваленко А.В., Никоненко В.В., Уртенов М.Х. Математическое моделирование электроконвекции в электромембранных системах с вынужденной конвекцией // Конденсированные среды и межфазные границы. – Воронеж. – 2011.-Т.13 - №4. - С. 492-498.
- Узденова А.М., Коваленко А.В., Никоненко В.В., Уртенов М.Х. Математическое моделирование электроконвекции в канале обессоливания электродиализатора с учетом вынужденной конвекции // Экологический вестник научных центров Черноморского экономического сотрудничества. – Краснодар. – 2011. - №4. – С. 86-93.
- Узденова А.М., Коваленко А.В., Уртенов М.Х. Математическое моделирование процесса обессоливания в ЭДА с учетом электроконвекции и геометрической неоднородности поверхности мембраны [Электронный ресурс] // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. - №72(08). - 2011. - С. 1-11. Режим доступа: http://ej.kubagro.ru/2011/08/pdf/18.pdf
- Узденова А.М., Коваленко А.В., Никоненко В.В., Уртенов М.Х. Причины возникновения электроконвекции в электромембранных системах [Электронный ресурс] // Политематический сетевой электронный научный журнал Кубанского государственного аграрного университета. - №73(09). - 2011. - С. 1-14. Режим доступа: http://ej.kubagro.ru/2011/09/pdf/33.pdf.
Публикации по теме диссертации в других изданиях:
- Узденова А.М., Уртенов М.Х. Условие скольжения и его использование в моделях мембранной электрохимии // Ионный перенос в органических и неорганических мембранах. Электромембранные технологии на базе фундаментальных исследований явлений переноса. Материалы российской конференции с международным участием. – Краснодар, 2008. - С. 231.
- Узденова А.М., Уртенов М.Х. Анализ влияния условия скольжения на течение электролита // Ионный перенос в органических и неорганических мембранах. Электромембранные технологии на базе фундаментальных исследований явлений переноса. Материалы российской конференции с международным участием. – Краснодар, 2008. - С. 232-233.
- Uzdenova A.M., Urtenov M. Kh. About a mathematical model of electroosmosis of the second kind // Ion transport in organic and inorganic membranes. International conference. Book of Abstracts. – Krasnodar. 2009. - P. 209.
- Uzdenova A.M., Urtenov M. Kh. Instability influence on processes in electro-membrane systems // Ion transport in organic and inorganic membranes. International conference. Book of Abstracts. – Krasnodar. 2009. - P. 210.
- Узденова А.М., Уртенов М.Х. Условие скольжения в задачах моделирования течения электролита в электромембранной системе // Современное состояние и приоритеты развития фундаментальных наук в регионах: Труды VI Всероссийской научной конференции молодых ученых и студентов. Краснодар: Просвещение-Юг, 2009. - С. 110-112.
- Узденова А.М., Уртенов М.Х. Моделирование электроконвекции в мембранных процессах, связанной с неоднородной электропроводностью поверхности мембраны // Современное состояние и приоритеты развития фундаментальных наук в регионах: Труды VI Всероссийской научной конференции молодых ученых и студентов. Краснодар: Просвещение-Юг, 2009. - С. 106-110.
- Узденова А.М. Алгоритм численного решения одной математической модели электроосмоса второго рода // Прикладная математика XXI века. Материалы IX объединенной научной конференции студентов и аспирантов факультета компьютерных технологий и прикладной математики. – Краснодар, 2009. - С. 115-117.
- Uzdenova A.M., Urtenov M.Kh. Using the COMSOL Multiphysics in mass transfer modeling in desalination channel of the electrodialysis device // Ion transport in organic and inorganic membranes. International conference. Book of Abstracts. – Krasnodar. 2010. - P. 193-195.
- Узденова А.М. Математическая модель стационарного течения электролита в примембранном слое // Прикладная математика XXI века: материалы X объединенной научной конференции студентов и аспирантов факультета компьютерных технологий и прикладной математики. – Краснодар, 2011. - С. 91-93.
- Узденова А.М., Уртенов М.Х. Численное решение задачи массопереноса в канале обессоливания электродиализного аппарата методом конечных элементов при выполнении условия электронейтральности // Современное состояние и приоритеты развития фундаментальных наук в регионах: Труды VII Всероссийской научной конференции молодых ученых и студентов. Краснодар: Просвещение-Юг, 2010. - С. 168-171.
- Узденова А.М., Уртенов М.Х. Анализ различных моделей стационарного переноса бинарного электролита в КО ЭДА // Современное состояние и приоритеты развития фундаментальных наук в регионах: Труды VII Всероссийской научной конференции молодых ученых и студентов. Краснодар: Просвещение-Юг, 2010. - С. 171-174.
- Uzdenova A.M., Kovalenko A.V., Urtenov M.Kh., Nikonenko V.V. Modeling of electroconvection in membrane system // Proceedings of the international conference “Ion transport in organic and inorganic membranes” June 6-11. – Krasnodar. 2011. - P. 205-207.