Численное решение краевой задачи модели переноса бинарного электролита в приближении закона Ома

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

Численное решение краевой задачи модели переноса бинарного электролита в приближении закона Ома

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

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

Чубырь Наталья Олеговна

Хромых Анна Александровна

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

Барсукова Виктория Юрьевна

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

Аннотация

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

Ключевые слова: математическое моделирование, система уравнений нернста - планка и пуассона, численные методы

Summary

Numerical solution transport model binary electrolyte in approximation of ohm's law

Kovalenko Anna Vladimirovna

Urtenov Mahamet Khuseevich

Chubyr Natalia Olegovna

Khromikh Anna Aleksandrovna

Uzdenova Aminat Machametovna

Barsukova Victoria Jurevna

Kovalenko Anna Vladimirovna

The article is devoted to the numerical solution of boundary value problem of the binary electrolyte model of transport in membrane systems in the approximation of Ohm's law. Different numerical methods are offered. The main regularities of transfer are established

Keywords: mathematical modeling, system nernst - planck and poisson equations, numerical methods

Введение

В работе [1] нами было найдено асимптотическое представление решения краевой задачи модели переноса бинарного электролита в мембранных системах в приближении закона Ома. Показано, что канал обессоливания , при каждом фиксированном разбивается на различные области , которые определяются знаком функции обобщенной концентрации : электронейтральности (), пространственного заряда (), промежуточного слоя (окрестность нулей функции ) и асимптотика решения имеет разный вид в различных областях. При «мягких» (допредельных) токовых режимах области и , значительно меньше области электронейтральности . При жестких (запредельных) режимах область пространственного заряда является уже макроскопической и сопоставимой с областью электронейтральности , но промежуточный слой остаются значительно меньшей, чем эти области. Система уравнений с частными производными, описывающая асимптотическое представление решения в каждой из областей меняет свой тип и в них встречаются члены, зависящие от знака искомой функции. Данная статья посвящена разработке метода численного решения таких краевых задач.

1. Краевая задача асимптотического представления решения

1.1 Система уравнений

Асимптотическое представление решения краевой задачи модели переноса бинарного электролита в мембранных системах в приближении закона Ома описывается следующей системой уравнений в безразмерных переменных [1]:

,

,

где - функция Хэвисайда [2]:

Определим дифференциальный оператор:

Тогда система уравнений (1-4) запишется в виде

(5)

Будем искать стационарное решение системы уравнений (5), (6). Поскольку оператор меняет свой тип, то удобно находить стационарное решение системы уравнений (5), (6) методом установления:

(8)

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

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

Постановка краевых условий зависит от моделируемого режима работы электродиализного аппарата.

Как показано в [3] при гальваностатическом режиме естественными являются следующие краевые условия.

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

2) Начальное условие

, (10)

3) Условия согласования граничных условий:

11)

4) Условия согласования граничных и начальных условий:

(12)

Для функции граничные условия и начальное условие будет согласовано, если взять , например, в виде , где мало, например .

2. Явный метод численного решения

2.1 Дискретизация области и уравнений

Область

разбивается с шагом , и по с шагом

.

Для простоты в дальнейшем рассматривается частный случай

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

Вводятся массивы четыре одномерных и шесть двумерных массива соответствующей размерности

, :

При фиксированном вычисляем массив по формуле

.

Аналогично вычисляются массивы . Через будем обозначать значения функции на прошлом слое по , а через на текущем слое по времени.

Вычислим связь , и , из разностных уравнений для

используя явную схему.

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

,

,

,

,

.

Аналогично заменяются конечными разностями и производные от функции : , причем:

Подставим эти выражения в уравнения системы, тогда получим:

,

если ,

если

, ,

где дискретный оператор дивергенции, который вычисляется по формуле:

.

Разрешим эти уравнения относительно , тогда получим:

при :

,

(14)

Выше используется формула

Присваивание значений на границах

(17)

2.2 Алгоритм численного решения

1 шаг. , присвоение начальных значений

, где

, где

2 шаг. , переход на следующий слой

a) Рассчитывается , по формулам (13) - (16), при

.

b) Рассчитывается массивы по формулам:

c) Рассчитываются ,, , , ,, , по формулам (17)

d) Выводим графики ,

3 шаг. Проверка достижения заданного конечного времени

Если , то

, ,

и переход к шагу 2, иначе выход из алгоритма.

3. Различные обобщения

3.1 Полунеявная схема

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

После подстановки в уравнение (7), получится система линейных уравнений для , которую можно решать, например, с помощью продольно-поперечной прогонки [4]. При этом для решения уравнения (8) используется явная схема.

3.2 Неявная схема. Метод итераций

Если все производные заменить конечными разностями на текущем слое, то получим систему нелинейных алгебраических уравнений относительно , . Часть уравнений относительно являются «квазилинейными» алгебраическими уравнения вида:

, (18)

где некоторые известные матрицы, причем функция Хэвисайда применяется к матрице покомпонентно. Для решения уравнения (18) можно использовать следующий метод итераций:

, , (19)

Уравнения (19) можно решать, например, продольно-поперечной прогонкой, причем в качестве начального приближения можно рассматривать значения функции на предыдущем слое, т.е. .

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

4.Некоторые результаты численного решения

Ниже приводятся решения уравнений при следующих данных:

, , ,

Рисунок 1 Графики решений краевой задачи при : а) - , b) - , с) - , d) - , e) - векторное поле плотности тока , , ,

Из рис. 1 следует, что в средней части канала имеется «плато» с практически постоянным значением обобщенной концентрации, равным начальному условию (рис.1а). Это «плато» с увеличением времени постепенно размывается. Вблизи анионообменных и катионообменной мембран распределение обобщенной концентрации практически линейное. Функция при фиксированном времени является практически линейным (рис.1b) по переменным . Из (рис.1c,d,e) следует, что поле плотности тока практически перпендикулярно поверхности мембран.

С увеличением времени значения обобщенной концентрации уменьшаются, при этом сохраняются указанные выше закономерности (см. рис.2a-e).

Рисунок 2 Графики решений краевой задачи при : а) - , b) - , с) - , d) - , e) - векторное поле плотности тока , , ,

Результаты проведенного выше исследования позволяют сделать следующие выводы:

1) Предложен эффективный метод численного решения краевой задачи модели переноса бинарного электролита в мембранных системах в приближении закона Ома.

2) Установлены основные закономерности переноса, среди, которых необходимо особо отметить вывод, что поле плотности тока практически перпендикулярно поверхности мембран. Это позволяет построить различные приближенные аналитические решения.

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

Литература

1 Уртенов К.М. Математическое моделирование тепломассопереноса в электродиализных аппаратах водоподготовки / К.М. Уртенов, А.В. Коваленко, Т.Л. Шапошникова - М.: Финансы и статистика, 2010. -- 214c.

2 Волков И.К. Интегральные преобразования и операционное исчисление: Учеб. для вузов / И.К. Волков, А.Н. Канатников, B.C. Зарубина, А.П. Крищенко -- 2-е изд. -- М.: Изд-во МГТУ им. Н. Э. Баумана, 2002. -- 228 с.

3 Узденова А.М. Математические модели электроконвекции в электромембранных системах/ Узденова А.М., Коваленко А.В., Уртенов М.Х. -- Карачаевск: КЧГУ, 2011. -- 154с.

4 Самарский А.А. Численные методы/ А.А. Самарский, А.В. Гулин, - М.: Наука, 1989.

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