Введение
Электродиализный аппарат имеет периодическую структуру, состоящую из чередующих камер обессоливания и концентрирования, а также двух электродных камер. Изменение концентрации в камере концентрирования можно учесть в граничных условиях. Таким образом, основной задачей является моделирование переноса в камере обессоливания. Пусть соответствует условной межфазной границе анионообменная мембрана/раствор, соответствует условной межфазной границе катионообменная мембрана/раствор, – входу, а – выходу из камеры обессоливания, – заданная скорость прокачивания раствора. Будем считать, что скорость течения раствора имеет форму параболы Пуазейля или моделируется по-другому.
Для математического моделирования явлений переноса для бинарного электролита в камере обессоливания используется краевая задача для системы уравнений Нернста – Планка и Пуассона [4], которая достаточно сложна для аналитического и численного решения.
2. Постановка задачи. Нами в работе [5], для 1:1 электролита с одинаковыми коэффициентами диффузии катионов и анионов, была выведена упрощенная модельная задача с функцией Хэвисайда [1] и показано, что она достаточна адекватна. В данной работе предлагается асимптотическое решение соответствующей краевой задачи, которая имеет безразмерный вид [1, 2]:
, (1)
, (2)
, (3)
,
, (4)
, , ,
, , (5)
где – искомая напряженность, , – концентрации катионов и анионов, – обобщенная общая «концентрация», – плотность тока, – функция тока для плотности тока , т.е. , , – безразмерный малый параметр, равный удвоенному квадрату отношения Дебаевской длины к ширине канала, , – число Пекле, – заданная скорость протока электролита в камере обессоливания, – функция Хэвисайда, – ток в цепи, – длина канала, – кососимметрическое скалярное произведение.
При запредельных режимах функция меняет знак в области по крайней мере, для небольших , причем дважды. Из этого следует, что функция имеет в разных частях области различные асимптотические представления.
Нами показано, что область , где значение функции положительно является областью электронейтральности, а отрицательно – пространственного заряда.
Рисунок 1. Область знакопостоянства функции
Особенностью задачи является то, что функция из уравнения (1) находится независимо от функций и .
3. Асимптотическое решение в области электронейтральности
3. 1. Асимптотическое разложение
Используем в области электронейтральности () для асимптотического решения следующие разложения:
, .
Подставляя данные разложения в систему уравнений и приравнивая коэффициенты при равных степенях , получим уравнения, решая которые можно найти все коэффициенты разложения. Для нахождения решения -го приближения получаем следующую систему уравнений:
(6)
(7)
где
3. 2. Алгоритм решения начального и первого приближения
Для начального приближения имеем систему уравнений:
,
,
которая после ряда преобразований принимает вид:
(8)
(9)
Из системы уравнений (9)–(10) видно, что вначале из уравнения (9) находится решение , далее, используя условие , вычисляем , затем находим .
Для приближения первого порядка получаем систему уравнений:
,
которая после некоторых преобразований преобразуется к виду:
(10)
(11)
Из системы уравнений (10)–(11) видно, что первое приближение находится по такому же алгоритму, что и начальное приближение.
4. Асимптотическое решение в области пространственного заряда
4. 1. Асимптотическое разложение
Сделаем замену: , и используем в области пространственного заряда () для асимптотического решения следующие разложения тогда
,
,
тогда для -го приближения имеем систему уравнений
(12)
(13)
Система уравнений (12) может быть записана в виде:
,
где .
Заметим, что для получаем систему линейных уравнений с нулевым определителем . Для однозначного определения - го приближения используется условие разрешимости - го приближения, которое имеет вид
(14)
4. 2. Алгоритм решения начального и первого приближения
Система уравнений для нулевого приближения имеет вид
(15)
(16)
Система уравнений для первого приближения имеет вид
(17)
(18)
Система уравнений для второго приближения для имеет вид:
(19)
Уравнение (15) можно записать в виде
(20)
Из этого уравнения не удается однозначно найти и , поэтому необходимо использовать условие разрешимости уравнения для , которое, как следует из (14), имеет следующий вид:
(21)
Решая систему уравнений (20) и (21), получаем однозначное выражение через :
, (22)
с учетом которого уравнение (18) преобразуется к виду:
(23)
Уравнение (23) является квазилинейным уравнением параболического типа.
Из уравнения (19) имеем следующее условие разрешимости для первого приближения
(24)
Решая систему уравнений (21) и (24), выражаем напряженность электрического поля через плотность тока для первого приближения
, (25)
где , .
Из уравнений (18) и (19) можно получить выражение для нахождения :
,
в которое из (21) подставляем и получаем уравнение для . Соответствующее уравнение не приводится из-за ее громоздкости.
5. Алгоритм численного решения начального приближения
Для численного решения вводится в рассмотрение дифференциальный оператор , который определяется следующим образом [3]:
,
Тогда для функций и получается система уравнений
(26)
(27)
с соответствующими краевыми условиями (4)–(5).
Для численного решения использовался метод установления, для дискретизации применялись явная и неявная разностные схемы. На рисунке 2 приведены графики решений краевой задачи для функции обобщенной концентрации и плотности тока .
а)б)
Рисунок 2. Графики решений краевой задачи: а) , б)
Заключение. Краевая задача (11)–(15) относится к классу «жестких» задач [1] из-за наличия малого параметра и поэтому при численном решении при значениях малого параметра порядка требуется шаг дискретизации порядка , что еще приемлемо по времени решения.
Однако реально меняется в пределах от до . При таких значениях малого параметра численное решение становится затруднительным. Предлагаемое асимптотическое решение позволяет провести численный анализ краевой задачи при произвольно малых значениях параметра .
Рецензенты:
Осипян Валерий Осипович, д.ф-м.н., доцент, профессор кафедры информационных технологий Кубанского государственного университета, г. Краснодар.
Семенчин Евгений Андреевич, д.ф.-м.н., профессор, заведующий кафедрой высшей алгебры и геометрии Кубанского государственного университета, г. Краснодар.
Криштоп Виктор Владимирович, д.ф.-м.н., профессор, заведующий кафедрой «Физика», Дальневосточный государственный университет путей сообщения, г. Хабаровск, профессор Университета Kwangwoon University, Korea.
Библиографическая ссылка
Хромых А.А., Чубырь Н.О., Коваленко А.В., Уртенов М.Х. НАХОЖДЕНИЕ ВЫСШИХ АСИМПТОТИЧЕСКИХ РАЗЛОЖЕНИЙ КРАЕВОЙ ЗАДАЧИ МОДЕЛИ С ФУНКЦИЕЙ ХЭВИСАЙДА // Современные проблемы науки и образования. – 2013. – № 4. ;URL: https://science-education.ru/ru/article/view?id=9842 (дата обращения: 07.10.2024).