Перенос ионов в трехслойных ионообменных мембранных системах при интенсивных токовых режимах — страница 7

  • Просмотров 375
  • Скачиваний 13
  • Размер файла 136
    Кб

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

одинаковы. Величина шага определяется автоматически быстротой изменения интегрируемых функций. Точка wi становится точкой разбиения исходного отрезка на подотрезки, если не выполняется хотя бы одно из условий , (1) где – интегрируемые функции, M – наперед заданная константа. Кроме того, вводится замена переменных: ; ; , (2) где C1 – концентрация противоионов; CА – концентрация коионов; Е – напряженность электрического поля.

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

типа NaCl через отдающий противоионы диффузионный слой толщины . В новых переменных краевая задача записывается в виде: (3) Для повышения надежности вычислительных итераций также использовался метод продолжения по параметру и предложенная А.Н. Тихоновым регуляризация метода Ньютона. В качестве параметра продолжения был выбран малый безразмерный параметр . В диссертационной работе получено решение задачи (3) для значений

малого параметра вплоть до 10-7, в то время как использование метода параллельной стрельбы с шагом постоянной длины позволяет получить решение только для (К.А. Лебедев). Таким образом, за счет модификации метода удалось понизить величину малого параметра, для которого метод дает устойчивое решение, на два порядка. Показано совпадение найденных решений при <10-5 с асимптотическими решениями, полученными М.Х. Уртеновым, а при

>10-5 с решениями, полученными К.А. Лебедевым методом параллельной стрельбы с постоянным шагом. Четвертая глава посвящена исследованию строения двойного электрического слоя (ДЭС) на межфазной границе. Рассматривается перенос ионов сильного электролита типа 1:1 с учетом пространственного заряда как в диффузионном слое, так и в фазе мембраны. Математическая модель представляет собой совокупность следующих уравнений: 