CC BY
2
Math-Net.Ru
Е. В. Казаковцева, А. В. Коваленко, А. В. Пись-менский, М. Х. Уртенов, Гибридный численно-аналитический метод решения задач переноса ионов соли в мембранных системах с осевой симметрией, Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2024, номер 1, 130-151
001: 10.14498^^2043
Использование Общероссийского математического портала Math-Net.Ru подразумевает, что вы прочитали и согласны с пользовательским соглашением
http://www.mathnet.ru/rus/agreement
Параметры загрузки:
IP: 109.252.33.182
29 сентября 2024 г., 12:20:43
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2024. Т. 28, № 1. С. 130—151_
ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu2043
EDN: BOXTTX
УДК 517.958:531.72
Гибридный численно-аналитический метод решения задач переноса ионов соли в мембранных системах с осевой симметрией
Е. В. Казаковцева, А. В. Коваленко, А. В. Письменский, М. Х. Уртенов
Кубанский государственный университет, Россия, 350040, Краснодар, ул. Ставропольская, 149.
Аннотация
Целью работы является разработка нового гибридного численно-аналитического метода решения краевых задач с осевой симметрией, например, с вращающимся мембранным диском, основанного на сращивании асимптотического решения вблизи катионнообменной мембраны (КОМ) с численным решением в оставшейся части области. Для это используется следующий метод:
1) предлагается базовая математическая модель переноса ионов соли в электрохимической ячейке с вращающимся катионнообменным мембранным диском на основе общих законов сохранения, представленных уравнениями Нернста-Планка-Пуассона и Навье-Стокса с естественными граничными и начальными условиями. Эта модель не содержит никаких подгоночных параметров или упрощающих предположений. Однако численное решение соответствующей краевой задачи представляет собой значительные вычислительные сложности при реальных концентрациях раствора и больших скачках потенциала и угловой скорости вращения мембранного диска, связанные с большими градиентами концентрации
Математическое моделирование, численные методы и комплексы программ Научная статья
© Коллектив авторов, 2024 © СамГТУ, 2024 (составление, дизайн, макет)
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru) Образец для цитирования
Казаковцева Е. В., Коваленко А. В., Письменский А. В., Уртенов М. Х. Гибридный численно-аналитический метод решения задач переноса ионов соли в мембранных системах с осевой симметрией // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2024. Т. 28, № 1. С. 130-151. EDN: BOXTTX. DOI: 10.14498/vsgtu2043. Сведения об авторах
Екатерина Васильевна Казаковцева А https://orcid.org/0009-0003-0040-0880 старший преподаватель; каф. анализа данных и искусственного интеллекта; e-mail: vivkaterina@mail. ru
Анна Владимировна Коваленко О) https://orcid.org/0000-0002-3991-3953 доктор технических наук, доцент; заведующий кафедрой; каф. анализа данных и искусственного интеллекта; e-mail: savanna- [email protected]
Александр Владимирович Письменский © https://orcid.org/0000-0003-4046-2229 кандидат физико-математических наук; заведующий кафедрой; каф. прикладной математики; e-mail: [email protected]
Махамет Хусеевич Уртенов © https://orcid.org/0000-0002-0252-6247 доктор физико-математических наук, профессор; профессор; каф. прикладной математики; e-mail: [email protected]
и потенциала вблизи КОМ в квазиравновесной области пространственного заряда (ОПЗ);
2) область решения разбивается на две части, одна из которых представляет собой небольшую область возрастания катионов (ОВК), расположенную вблизи КОМ, и оставшуюся основную часть области (ОЧО);
3) в ОВК находится аналитическое решение методом сращивания асимптотических решений;
4) в ОЧО строится упрощенная математическая модель, которая отличается от базовой математической модели таким граничным условием на границе с ОВК, которое позволяет затем произвести сращивание решения соответствующей краевой задачи с решением в ОВК.
Основным результатом является гибридный численно-аналитический метод, который дает возможность проводить численный анализ переноса ионов соли при реальных концентрациях раствора электролита бинарной соли при широком диапазоне изменения скачка потенциала и угловой скорости вращения мембранного диска.
По результатам работы можно сделать следующий вывод: сочетание аналитического (асимптотического) метода решения в области по-гранслоя и численного решения в остальной области, за исключением погранслоя, с их последующим сращиванием позволяет построить эффективный гибридный численно-аналитический метод решения задач переноса ионов соли в мембранных системах с осевой симметрией.
Ключевые слова: асимптотическое решение, гибридный численно-аналитический метод, вращающийся мембранный диск, осевая симметрия, цилиндрическая система координат.
Получение: 22 июля 2023 г. / Исправление: 12 апреля 2024 г. / Принятие: 15 апреля 2024 г. / Публикация онлайн: 17 мая 2024 г.
Введение. Гидродинамические системы с осевой симметрией в настоящее время широко используются и интенсивно исследуются [1—8]. Так, например, в работе [9] найдено точное решение в рамках уравнений Эйлера закрученных осесимметричных стационарных течений идеальной несжимаемой жидкости, а в работе [10] в рамках уравнений Навье—Стокса рассмотрены нестационарные осесимметричные течения однородной вязкой несжимаемой жидкости, в которых осевая и окружная скорости зависят только от радиуса и от времени, а радиальная скорость равна нулю. Произведено расщепление краевых задач для рассматриваемого типа течений на две задачи, каждая из которых содержит две неизвестные функции (давление и одна из компонент скорости). Актуальным является исследование аналогичных задач, но с дополнительными процессами переноса вещества за счет конвекции, диффузии, электромиграции и т.д. Изучению свойств переноса ионов соли в системах с вращающейся жидкостью, в том числе методом вращающегося мембранного диска, посвящены ряд работ [1-8]. Эта работа является продолжением работы [3], где приведена базовая модель, основанная только на основных законах сохранения и естественных граничных условиях и не использующая никаких гипотез для упрощения. Базовая модель служит для проверки исследования основных свойств и вывода упрощенных моделей на основе этих свойств. Достоинствами базовой модели являются общность и отсутствие подгоночных параметров, к недостаткам можно отнести сложность расчета при высоких
скоростях вращения мембранного диска и скачке потенциала, а также реальных концентрациях раствора. Это связано с большими градиентами концентрации и потенциала вблизи катионнообменной мембраны (КОМ). Для преодоления указанных ограничений предлагается разбиение области решения на две подобласти: первая — основная часть области (ОЧО) вдали от катионообменной мембраны; вторая — область возрастания катионов (ОВК) вблизи КОМ, границей которой служит точка локального минимума концентрации катионов около КОМ. Для каждой из этих областей из базовой модели выводятся упрощенные модельные задачи, удобные как для численного решения (в ОЧО), так и для аналитического (в ОВК). Аналитическое решение зависит от того, больше или меньше ток, протекающий через систему, по сравнению с предельным током [11]. В допредельном режиме ОВК совпадает с квазиравновесной областью пространственного заряда (ОПЗ), однако в сверхпредельном режиме ОВК должна содержать наряду с квазиравновесной ОПЗ еще и промежуточную область (слой), необходимую для сращивания решений. Решения в этих подобластях сращиваются, что позволяет искать решение при большой угловой скорости, начальной концентрации и скачке потенциала. Модель без области возрастания катионов назовем упрощенной моделью без ОВК. Аналогичная модель переноса ионов соли для плоского канала впервые была предложена, по-видимому, в работе [12], однако в этой работе не было исследовано решение в ОВК и поэтому не было произведено сращивание.
1. Формулировка базовой модели и модели без ОВК в задачах с осевой симметрией. В допредельном режиме область решения состоит из области электронейтральности (ОЭН) и квазиравновесной ОПЗ. В сверхпредельном режиме между ОЭН и квазиравновесной ОПЗ возникает расширенная ОПЗ. В расширенной ОПЗ электромиграция преобладает над диффузией, а в области электронейтральности миграционный и диффузионный потоки равны. В квазиравновесной ОПЗ миграционный поток равен диффузионному, но они противоположны по направлению, что приводит к тому, что ток становится равным нулю, а также в этой области концентрации и напряженность электрического поля резко возрастают, поэтому их градиенты имеют значительные величины, а это затрудняет численное решение в данной области. Поэтому возникает необходимость в формулировке упрощенных моделей, которые дают решение, близкое к решению базовой математической модели, и позволяют выполнять численные расчеты при большой угловой скорости, начальной концентрации и скачке потенциала.
Приведем краевую задачу базовой математической модели, которая состоит из области исследования, системы уравнений, а также краевых условий.
1.1. Область исследования. Для моделирования течения жидкости и переноса ионов соли в электрохимической ячейке с вращающимся кати-оннобменным диском [2] благодаря осевой симметрии можно описать только половину сечения цилиндрической области (рис. 1), которую обозначим через О.
1.2. Система уравнений. В области О для описания переноса ионов соли используются связанная система уравнений Навье—Стокса с объемной электрической силой / и Нернста—Планка—Пуассона, которые в цилиндрической системе координат с учетом осевой симметрии принимают следующий
Рис. 1. Исследуемое сечение области О и ее границы: 1 — глубина раствора, где выполняется условие электронейтральности; 2 — ось симметрии; 3 — катионнообменная мембрана;
4 — открытая граница [Figure 1. Investigated cross-section of the region О and its boundaries: 1—depth of the solution where the electroneutrality condition is satisfied; 2—axis of symmetry; 3—cation-exchange membrane; 4—open boundary]
вид:
du / du v2 du \ dp P dt P\dr r W dz) dr ^
1 d / du \ u d2u r dr\ dr J r2 dz2
+ fr
dv ( dv uv dv \
Pd~t + РГд~г + )="
dw / dw dw \ dp P dt P\dr W dz) dz
1 d / dv\ v d 2v d d 2 d 2
1 d / dw \ d2w .rdr\ dr J dz2
+u
+ fz;
где
, p dC. Ni = ZiDiCiEr - Di -d^ + C\u) er + CiV+
F _ „ „ dC
+ ( КГ ^DiCiEz -Di d^ + Clw^ez, i = 1, 2; (1)
dCi _ 1 d ( N \ dNi,z . _ (2]
dt rdr г,г dz ' ' '
-2Ф + -2Ф = -1 -Ф - F ,C +zC , (3) dr2 dz2 r dr e
I = F (ziN! + Z2N2), (4)
19 / дФ \ d2 Ф е(д2ф + °2ф )дФ +£1 (дФ V
г \г дг\ дг у дг2 ) дг V дг2 дг2 ) дг г\дг / '
и = 0,
_ /1 д /дФ\ д2Ф\ дФ _ £(д2Ф + д2Ф\дФ +£1дФдФ г \гдг\ дг / дг2 у дг \дг2 дг2 / дг г дг дг1 и, V, и> — радиальная, азимутальная и аксиальная компоненты скорости течения раствора; N1, N2, С1, С2 — потоки и концентрации катионов и анионов
в растворе соответственно; , 22 — зарядовые числа катионов и анионов; I — плотность тока; И2 — коэффициенты диффузии катионов и анионов соответственно; Ф — потенциал электрического поля; Е = —УФ — напряженность электрического поля; е — диэлектрическая проницаемость электролита; Р — постоянная Фарадея; К — газовая постоянная; Т — абсолютная температура; Ь — время; р — плотность; ц — динамическая вязкость; р — давление.
1.3. Краевые условия для базовой модели. Граница 1 является входом для раствора; для концентраций выполняется условие электронейтральности; задается постоянная начальная концентрация СО; для скорости ставится условие отсутствия нормального напряжения (УЙ + (УЙ)Т)п = 0; задается нулевое давление; граница 1 считается также эквипотенциальной (анодом) — Ф(*,г, 0) = Оф.
Граница 2 является осью симметрии.
Граница 3 соответствует вращающейся идеально селективной КОМ и является выходом для катионов, концентрация которых постоянна и определяется емкостью мембраны — С\(Ь, г, Н) = С^т. КОМ предполагается идеально селективной и поэтому для анионов задается условие отсутствия потока — —п ■ М2 = 0; для потенциала Ф(Ь,г,Н) = 0 задается азимутальная скорость в виде г, Н) = шг, а остальные компоненты скорости принимаются равными нулю.
Граница 4 является выходом для раствора; ставится условие выноса ионов
только конвективным потоком — = —и ■ Сг, г = 1, 2; для потенциала за-> / дФ дФ \т п
дается условие — п ■ [г-¡^, = 0; для скорости предполагается условие
отсутствия нормального напряжения— (Уй + (УЙ)Т)п = 0.
В начальный момент времени концентрация ионов соли предполагается постоянной — С\(0, г, г) = С2(0, г, г) = С0, а раствор — неподвижным.
1.4. Краевые условия для модели без ОВК. Все начальные и граничные условия остаются такими же, как в базовой модели, за исключением условия на границе ОВК, для которой задается концентрации катионов — —пУС1\_ =0, где гт — координата минимума концентрации, зависящая от угловой скорости вращения диска и начальной концентрации. Поскольку толщина области возрастания катионов (ОВК) мала, т.е. хт ~ Н, это условие может быть заменено на условие — пУС\ \х—н = 0. Таким образом, модель без ОВК будет определена во всей ячейке. Ниже показано, что ее решение близко к решению базовой модели в ОЧО, но, естественно, их решения отличаются в ОВК.
1.5. Сравнительный анализ численных результатов исследования базовой модели и модели без ОВК. Были проведены исследования для разнообразных исходных данных (концентрация, угловая скорость и т.д.), соответствующих сверхпредельному режиму, когда возникают электроконвективные вихри. В качестве примера ниже приведены результаты вычислений при следующих условиях:1 Со = 0.01 моль/м3, ш = ж/2, йф = 1.5 В, радиус вращающегося мембранного диска г = 1 мм, высота ячейки — 2 мм.
На рис. 2 приведены линии тока жидкости в базовой модели и для модели без ОВК, которые показывают образование электроконвективного вихря,
хЭти данные соответствуют конструкции электрохимической ячейки с вращающимся мембранным диском, приведенной в [2].
причем линии тока раствора практически полностью совпадают.
Surface: Concentration (mol/m3) Streamline: Velocity field
Рис. 2. Линии тока жидкости при Ь = 1000 е в базовой модели (а), в базовой модели около катионнообменной мембраны (Ь), в модели без области возрастания катионов (с), в модели без области возрастания катионов около катионнообменной мембраны (^
[Figure 2. Fluid streamlines at t = 1000 s in the base model (a), in the base model near the cation-exchange membrane (b), in the model without the cation increase region (c), in the model without the cation increase region near the cation-exchange membrane (d)]
Концентрации, потенциал и т.д. зависят от t, г и z. Для того чтобы сделать обозримой их зависимость от z, рассмотрим графики при фиксированных г начиная со значения г = 0.1 мм и до значения г = 0.9 мм с шагом 0.1 мм. На всех приведенных ниже рисунках по оси абсцисс показаны координаты высоты ячейки z, причем z = 0 соответствует глубине раствора, где выполняется условие электронейтральности, а z = 2 мм — условной границе раствор/КОМ. По оси ординат показаны требуемые величины (концентрации, потенциал, компоненты скорости и т.д.).
На рис. 3, а показаны профили концентрации катионов С\ в базовой модели, по которым можно сделать вывод, что концентрация сначала почти
постоянна, а затем линейно убывает, что соответствует поведению в области электронейтральности, далее образуется расширенная ОПЗ и концентрация достигает своего минимума, а потом нелинейно возрастает в квазиравновесной ОПЗ (рис. 3, Ь).
На рис. 3 показаны графики концентрации катионов в базовой модели и модели без ОВК, из которых можно сделать вывод, что они совпадают с большой точностью в области электронейтральности и расширенной ОПЗ, однако около КОМ отличаются. В базовой модели концентрации достигают минимума в некоторой точке хт, расположенной вблизи КОМ, и далее увеличиваются до граничного значения С1,т, в то время как в модели без ОВК концентрации продолжают монотонно убывать, причем в области (хт, 1] практически постоянны, что объясняется разницей в граничном условии для катионов около КОМ.
Графики концентрации анионов С2, показанные на рис. 4 для базовой модели и в модели без ОВК, ведут себя примерно одинаково, включая область около КОМ.
Для анализа распределения пространственного заряда рассмотрим разность концентраций С1 — С2, которая представляет собой плотность распределения пространственного заряда, нормированную на число Фарадея Р. Из анализа графиков разности концентраций С1 — С2 (рис. 5, Ь, с) в базовой модели и модели без ОВК следует их совпадение с высокой точностью в области электронейтральности и расширенной ОПЗ, однако, как показано на рис. 5, а, с базовой модели дополнительно формируется квазиравновесная ОПЗ.
На рис. 6 показано визуальное совпадение с большой точностью графиков потенциала базовой модели и модели без ОВК при небольших С^,т — порядка Со. При этом график разницы потенциалов этих двух моделей (рис. 7) показывает небольшие отличия (максимум которых достигается при г = 0.1 мм и в сечениях от г = 0.3 мм до 0.5 мм). Это связано с тем, что при данных условиях скачок потенциала в ОВК мал.
Если Ск,т значительно больше Со, например, в тысячу и более раз, что вполне реально, то тогда разница будет более существенной и она компенсируется скачком потенциала в ОВК, рассчитанным аналитически.
По графикам радиальной скорости (см. рис. 8) базовой модели и модели без ОВК можно сделать вывод об их высокой степени совпадения.
На рис. 9 и рис. 10 приведены графики азимутальной и аксиальной компонент скорости для базовой модели и модели без ОВК соответственно, которые показывают их практически полное совпадение.
Таким образом, модельная задача приближает решение базовой задачи с высокой точностью всюду за исключением ОВК, где решение должно быть дополнено аналитическими формулами, полученными асимптотическим методом.
2. Аналитическое решение краевой задачи в ОВК. Аналитическое решение в ОВК зависит от величины скачка потенциала, а именно допредельного или сверхпредельного. В допредельном режиме ОВК совпадает с квазиравновесной ОПЗ, а в сверхпредельном режиме эти области очень близки. Из численного исследования базовой модели следует, что толщина ОВК не зависит от радиальной координаты г, радиальная (и) и аксиальная (-ч) компоненты скорости близки к нулю, а азимутальная ( ) компонента скорости
Line Graph: Concentration (mol/m)
Line Graph: Concentration (mol/m5)
Line Graph: Concentration (mol/m3)
0.009
\V1
0.G075 0.007
1
1
S
<s0, 4
0.0035 0.003 0.0025 0.002 0.0015 0.001 §§§§§§§§§ r-0.6mm r=0.8mm
u
Vm
i i
z-coordlnate (m) b
le Graph: Concentration (mol/m3)
со -J
Рис. 3. Концентрация C\ при t = 1000 с и разных сечениях по г в базовой модели (а), в базовой модели вблизи катионнообменной мембраны (Ь), в модели без области возрастания катионов (с), в модели без области возрастания катионов вблизи катионнообменной мембраны (d)
[Figure 3. Concentration С\ at t = 1000 s and different cross sections for r in the base model (a), in the base model near the cation-exchange membrane (b), in the model without the cation increase region (c), in the model without the cation increase region near the cation-exchange
membrane (d)[
Une Graph: Concentration (mol/m5)
0.4mm 0.6mm O.Bmm
Une Graph: Concentration (mol/m )
Line Graph: Concentration (mol/m*)
If
0.007 M
0.004 0.003 0.002 0.001
§§§§§§§§§ 1.3mm D.6mm D.8mm
11
\\ \ \ \
\\ \ \ \ -e- Time-lOOOs. f-0.
\\ \ \ \ —Tima=1000s. r=0. 5mm
\\ \ \ \ Tlme= 1000 s. r-0 6mm
\\ \ \ \ —A— Tlme=i000s, f=0. Time-1000 s, r=0. ^mm
\\ \ \ \
\\ \ \ \
V\ \ \\
« \ W
v\\ v\
\\\
...........X.
0
Une Graph; Concentration {mol/m3!
0,0016 0.0015 0.0014 ООО s, ООО s, ООО s, ООО s, r=0.1mm r=0.2mm
0.0012 ООО S, ООО S, г=0.6шш
TIme= 000 s, r-0.9mm
0.0009 0.0008 0.0007 0.0006 0,0005
\ \ \
N \4 \ \
0
Рис. 4. Концентрация С'з при t = 1000 с и разных сечениях по г в базовой модели (а), в базовой модели вблизи катионнообменной мембраны (Ь), в модели без области возрастания катионов (с), в модели без области возрастания катионов вблизи катионнообменной мембраны (d)
[Figure 4. Concentration С'з at t = 1000 s and different cross sections for r in the base model (a), in the base model near the cation-exchange membrane (b), in the model without the cation increase region (c), in the model without the cation increase region near the cation-exchange
membrane (d)[
шшш
ШШШ
ШШШ
нш ж
НИМ уш М ! IГ
1
ШШШ
ШШШ
нш ж
I I
I* I
1
ШШШ
и~ «г «* «* и* «* «*
¡¡¡¡¡¡¡¡I
ШШШ жн ж
I I ! I I I I I I I I I ! I
(£шЛош) г>хэ
1'
Line Graph: Potential (V)
Graph: Potential (V)
Рис. 6. Графики потенциала при t -[Figure 6. Plots of the potential at t -
- 1000 с и разных сечениях по г в базовой модели (а), в модели без области возрастания катионов (Ь)
1000 s and different cross sections for r in the base model (a), in the model without the cation increase
region (b)[
Рис. 7. Графики разности потенциалов между базовой моделью и моделью без области возрастания катионов при
t = 1000 с и разных сечениях по г [Figure 7. Plots of the potential difference between the base model and the model without the cation increase region at t = 1000 s and different cross sections for r]
Graph: Velocity field, r component (m/s)
Рис. 8. Графики радиальной скорости при t = 1000 с и разных сечениях по г в базовой модели (а), в модели без области возрастания катионов (b) [Figure 8. Plots of the radial velocity at t = 1000 s and different cross sections for r in the base model (a), in the model without the
cation increase region (b)[
Line Graph: Velocity field, phi component (m/s)
Line Graph: Velocity field, phi component (m/s)
Рис. 9. Графики азимутальной скорости при t = 1000 с и разных сечениях по г в базовой модели (а), в модели без области возрастания катионов (b) [Figure 9. Plots of the azimuthal velocity at t = 1000 s and different cross sections for r in the base model (a), in the model without
the cation increase region (b)]
Рис. 10. катионов
0.006 0.0055 -0,005 -0.0045 -0-004 -0.0035 -0.003 -0,0025 -0.002 -0,0015 -0.001 -0,0005 -
-1000 s
1000 s Concentration <C3I (mol/m3) - 1000 S C1-C2 (C1-C2) (moi/m'l
Рис. 11. Графики аксиальной скорости при t = 1000 с и разных сечениях по г в базовой модели (а), в модели без области возрастания катионов (b) [Figure 11. Plots of the axial velocity at t = 1000 s and different cross sections for r in the base model (a), in the model without the
cation increase region (b)]
0 0.0002 0.0004 O.OOOG 0.0000 0.001 0.0012 0.0014 0.0016 0.0010 0 0.0002 0.0004 0.0006 0.0000 0.001 0.0012 0.0014 0.0016 0.0018
^coordinate (m) ^coordinate (m)
a b
Графики аксиальной скорости при t = 1000 с и разных сечениях по г в базовой модели (а), в модели без области возрастания (b) [Figure 10. Plots of the axial velocity at t = 1000 s and different cross sections for r in the base model (a), in the model without the
cation increase region (b)[
Une Graph: (mol/m5) Line Graph: Concentration (mol/m') Line Graph: C1-C2 (mol/m5)
a
OJ
a W
О
to Д a to a
f4
to ^
о to a
a
a w
о ^
to
¡4
a
о
№
§
a
о W a a<
to
почти постоянна в ОВК. Поэтому в уравнении (1) останется только третья компонента:
F дС■
N,z = ™ZtDtCtEz - А , i = 1, 2.
RT dz
Кроме того, толщина ОВК также практически не зависит и от времени t, т.е. наблюдается квазистационарный режим. В результате левая часть уравнения (2) равна нулю, а в правой части в силу независимости от г равна нулю производная потока по z, кроме того, в уравнениях (3) и (4) также остаются составляющие, зависящие только от z. Таким образом, для нахождения неизвестных функций в ОВК получается следующая система уравнений:
F дС
Ni}Z = -RTZiDiCiEz - Diд-Л, i = i, 2; (5)
^ = 0, i = 1,2; (6) д
d2 Ф F
M = - F (^ + ^2); (7)
Iz = F (ZlNUz + Z2N2,z ). (8)
Для асимптотического решения системы (5)—(8) выполним переход в области ОВК (zm ^ z ^ Н) к безразмерному виду с использованием характерных величин: z(u) = z/Н; i(u) = Сг/С0; N(u) = Nl/N0; d(u) = Dl/D0;
Üu) = Iz//о; = Ф/Фо; Eif) = Ez/Eo; e^ = е/Ъ0. Здесь Со
— начальная
концентрация раствора; Н — высота ячейки; Ио — коэффициент диффузии электролита; Ф0 = КТо/Е — тепловой потенциал; 10 = ЕИоСо/Н — характерный ток; N0 = ИоСо/Н — характерный поток ионов; Ео = Фо/Н — отношение теплового потенциала к характерной ширине; Ъо = Н2РСо/Фо — характерная величина, имеющая размерность электрической постоянной.
Параметр е(и является малым параметром, имеющим смысл отношения квадрата толщины области пространственного заряда к квадрату высоты ячейки. Другими словами, это отношение (квадратов) толщин заряженной и незаряженной областей раствора.
Параметр е(и может быть записан в виде £ ^ 0.58 ■ 10-12/Со.
При начальной концентрации Со = 0.01 моль/м3, использованной для численного расчета, получаем е(и ~ 5.8 ■ 10-11. Поэтому в дальнейшем е(и считается малым параметром.
Система уравнений (5)-(8) в безразмерной форме, например, для раствора КаС1, с учетом идеальной селективности КОМ принимает следующий вид (индекс «и» для простоты записи опущен):
От = СЕ - 1Я, *Ст = -С2Е*, е~Е = С1 -С2, (9)
аг аг аг
где Е(г, е) = -ОаФ/Оаг — напряженность электрического поля, а Ф(г, е) — его потенциал; С1(г, е), С2(г, е), 1г — искомые концентрации катионов и анионов соответственно; 1г —ток, равный потоку катионов; е > 0 — малый параметр.
В качестве краевых условий в безразмерном виде при г = гт/Н ставятся условия сращивания с упрощенной моделью, а при = 1 — условия
С1 \z=1 = С1,М,
_СдФ +
dz dz
= 0, Ф\г_1 = 0. (10)
z=1 Z=1
2.1. Решение в допредельном случае. Система уравнений (9) представляет собой сингулярно-возмущенную задачу. Проведя замену
£ = (г - 1)//, Ег = %, е)//е, С1(г, е) = С^, е), С2(г, е) = С2({, е)
можно убедиться, что в допредельном режиме ОВК совпадает с квазиравновесной ОПЗ.
Действительно, эти замены при малых е в начальном приближении сводят систему уравнений к классической системе уравнений Больцмана—Дебая, которая и описывает квазиравновесную ОПЗ [11]
§ = С1Е, § = | = С1 -С2 (11)
с соответствующими краевыми условиями, которая имеет аналитическое решение:
Е{г е) = , (12) / 1 -р ехр(/-40.^ - 1)//е)
, . ейЕ еЕ2 а . ,
С (^, £) = 2^е + ~Г - а, (13)
Е Е2 а
С ^, £) = - 2*; + -Г - а, (14)
где
а = -(С1(гт, е)+С2(гт, е)) ^-{С1(1, е) + С2(1, е)) < 0, Р — некоторое положительное число, которое определяется из условия
Е(0, е) = у/2(С1>т + а),
выведенного из граничных условий (10) с использованием первого интеграла системы уравнений (11) и эквивалентного условию С1(1, е) = С\)т.
При допредельном режиме нет расширенной ОПЗ, поэтому Е(гт - 0, е) — значение численного решения (решение в ОЧО в точке хт) должно быть конечным и, соответственно, Е(хт + 0, е) (значение аналитического решения в ОВК) тоже должно быть конечным. Это условие выполняется, если хт = = 1 - ку/е 11п е\, тогда получаем Е (гт + 0, е) = 4у/Ру/-а, при е ^ 0, если взять к = 1/(2у/-а). Равенства Е(гт - 0, е) = Е(гт + 0, е) можно добиться с использованием уже следующего первого приближения.
Сравнение точек хт, полученных при допредельном режиме в результате численного решения базовой модели, где хт вычисляется как точка, в которой достигается минимум концентрации, и аналитического решения в ОВК показывает их совпадение с высокой точностью: в базовой модели хт ~ 1.976 мм, а в модели ОВК хт ~ 1.999 мм.
Концентрации сращиваются в силу выбора а и формул (13), (14).
2.2. Решение в сверхпредельном случае. В сверхпредельном режиме наряду с квазиравновесным погранслоем возникает и расширенная ОПЗ [11], где концентрации ионов малы (рис. 11) (имеют порядок О (у/е)) и поэтому
процедура сращивания в допредельном режиме не может быть использована и должна быть модифицирована. Возьмем сечение = 0.5 мм и построим графики концентрации катионов и анионов, а также функцию С1 - С2 (см. рис. 11). Из этого рисунка видно, что в сверхпредельном состоянии ОВК состоит из двух частей: квазиравновесной ОПЗ и промежуточной области, которая служит для сращивания решения в квазиравновесной ОПЗ и в расширенной ОПЗ, являющейся частью ОЧО. Причем на правой границе промежуточной области сращивается напряженность электрического поля, а на левой границе — концентрация катионов. Точка гт, где производная концентрации равна нулю, близка к 1. При допредельном режиме гт и гк совпадают, при сверхпредельном режиме гт < гк (см. рис. 11).
При асимптотическом решении краевой задачи для системы (9) в области (гт, 1] в сверхпредельном режиме примем естественное для этого режима предположение об отсутствии в этой области анионов, то есть С2(г, е) = 0, г е (гт, 1].
Тогда получаем систему уравнений
=СЕ - 1г, ^ = С1,
а а
которая имеет аналитическое решение в общем виде:
= 1 46-Рехр(6 (г - 1)/-) Е(г, £)= — 1 -Рехр(26(г - 1)/^), (15)
= 4Ъ2-Рехр(б(г - 1)/—)(1 + Рехр(26(г - 1)/—)) (1б) (1 -Рехр(26(г -1)/^))2 ,
где > 0.
2.3. Сращивание решений и определение констант. Возьмем гк = 1 - к1 — +----^ 1, е ^ 0 + .
Для сращивания напряженности Е( , ) в точке к применим соотношение
Е( к + 0, ) = Е( к - 0, ),
где в Е(гк+0,е) используется решение в интервале (гк, 1], а в Е(гк-0,е) —решение в расширенной ОПЗ Е(г, е) = д/2( 1гг - I)/^/е, которое продолжается в интервал ( т, к). Таким образом получаем, что
1 4Ь—Рехр(-к1 Ь) _ \/2( 1г - I)
Отсюда
Положим
у/ё1 -Р ехр(-2к1 Ь) ^е
к = _11п(-46УР + у/8Р(2Ь2 + 1г - I)ч 1 Ь Ч 2Ру/2( 1г -1) У
гт = 1 -к2у[е11п +-----> 1, е^- 0 + .
Условие сращивания концентрации С1(г, е) в точке хт имеет вид
С1(г 0, е) = С1 (г + 0, £),
где в С1(гт - 0, е) используется решение в расширенной ОПЗ С1(г, е) = = Ьу/ё/\/2( 1гг - I), а в С1 (гт + 0, е) —решение из интервала (^, 1], продолженное в интервал (гт, гк). Откуда получим, что
к2 = \, Р = (-)2.
2 2 Ъ ' Мъ2^2( 1Х - I)) Для нахождения Ь используем условие С1(1, е) = С1т, откуда
4-
Ь=1 С1, т(1 -Р)2
2У уР(1+Р) .
3. Алгоритм гибридного численно-аналитического решения.
1) Численно решаем модель без ОВК и находим С]^(1, е), С2(1, е).
2) Находим численно скачок потенциала для упрощенной модели в ОЧО. Находим скачок потенциала для базовой, используя соотношение
Ф0 = [ Е(г, е)йг = Jo
[■¿т Г^к Г1
/ Е(г, е)йг + Е(г, е)йг + Е(г, е)йг =
'0 -'хт хк
Г Г ¿к Г1 ЙГ\
= Е(г, е)йх + Е(г, е)йг - =
Jo Jzm ]хк С2
= - 1п С2(1, £\ + П Е(г, е)йх+ Г Е(г, е)йг. С2(гк, £) Jo Лт
С учетом того, что Хк ~ ~ 1, получаем
1
С2(гк, е) ' ,}о
С (1 , ) 1
Фо 1п с2; , \ + Е(г, е)йг или Фо и Фовк + Фочо. С2(гк, е) Уо
Здесь первое слагаемое Фовк — скачок потенциала в области возрастания катионов (около КОМ), а второе слагаемое Фочо —скачок потенциала, численно рассчитанный в ОЧО.
3) Находим аналитическое решение в ОВК по формулам (12)—(14) при допредельном токовом режиме и по формулам (15), (16) при сверхпредельном.
4) Используя шаги 1) и 3), получаем решение базовой задачи.
Заключение. В статье предложен новый гибридный численно-аналитический метод, позволяющий находить решение базовой задачи при произвольных скоростях вращения мембранного диска и начальной концентрации. Это достигается путем сращивания численного решения упрощенной модели без
ОВК с асимптотическим решением в области возрастания катионов. Предложенный метод позволяет проводить численный анализ переноса ионов соли в реальных концентрациях раствора электролита бинарной соли при широком диапазоне изменения скачка потенциала и угловой скорости вращения мембранного диска.
Конкурирующие интересы. Заявляем, что в отношении авторства и публикации этой статьи конфликта интересов не имеем.
Авторский вклад и ответственность. Е.В. Казаковцева — разработка численно-аналитического метода и программы численных экспериментов; компьютерное и математическое моделирование; интерпретация полученных результатов; работа с черновиком и переработанным вариантом рукописи. А.В. Коваленко — компьютерное моделирование; обработка и анализ экспериментальных данных; подготовка первичного варианта рукописи. А.В. Письменский — обработка и анализ экспериментальных данных, интерпретация полученных результатов. М.Х. Уртенов — идея исследования; математическое моделирование; разработка численно-аналитического метода; интерпретация полученных результатов; подготовка первичного варианта рукописи; работа с черновиком и переработанным вариантом рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Финансирование. Исследование выполнено за счет гранта Российского научного фонда № 24-19-00648, https://rscf.ru/project/24-19-00648/.
Благодарность. Авторы благодарны рецензентам за тщательное прочтение статьи, ценные предложения и комментарии.
Библиографический список
1. Бурмашева Н. В., Просвиряков Е. Ю. Точное решение уравнений Навье-Стокса, описывающее пространственно неоднородные течения вращающейся жидкости // Тр. ИММ УрО РАН, 2020. Т. 26, №2. С. 79-87. EDN: IAWMLK. DOI: https://doi.org/10.21538/ 0134-4889-2020-26-2-79-87.
2. Заболоцкий В. И., Шельдешов Н. В., Шарафан М. В. Исследование электромассопере-носа хлорида натрия через катионообменную мембрану МК-40 методом вращающегося мембранного диска// Электрохимия, 2006. Т. 42, №12. С. 1494-1500. EDN: HYLEUX.
3. Казаковцева Е. В. Теоретическое исследование квазиравновесной области пространственного заряда в мембранных системах с осевой симметрией // Перспективы науки, 2023. №6(165). С. 58-68. EDN: DRJSOK.
4. Просвиряков Е. Ю. Восстановление радиально-осевой скорости в закрученных осесим-метричных течениях вязкой несжимаемой жидкости при лагранжевом рассмотрении эволюции завихренности // Вестн. Удмуртск. ун-та. Матем. Мех. Компьют. науки, 2021. Т. 31, №3. С. 505-516. EDN: ORVWHT. DOI: https://doi.org/10.35634/vm210311.
5. Чубырь Н. О., Коваленко А. В., Уртенов М. Х., Гудза И. В. Математическая модель стационарного переноса ионов соли в сечении канала при равновесии // Моделирование, оптимизация и информационные технологии, 2022. Т. 10, №3(38). EDN: DIECUK. DOI: https://doi.org/10.26102/2310-6018/2022.38.3.009.
6. Achoh A., Melnikov S., Bondarev D. Electrochemical characteristics of the MF-4SK membrane doped with the hyperbranched phosphorylated dendrimer BOLTORN H20 / Ion Transport in Organic and Inorganic Membranes: Conference Proceedings (Sochi, 22-27 May 2023). Krasnodar, 2023. pp. 15-17. EDN: RPXPHA.
7. Bondarev D., Eterevskova S., Zabolotsky V., et al. Homogeneous anion-exchange membrane with heterocyclic functional groups / Ion Transport in Organic and Inorganic Membranes: Conference Proceedings (Sochi, 22-27 May 2023). Krasnodar, 2023. pp. 31-32. EDN: BPLDVU.
8. Melnikov S. Experimental and theoretical study of ion transport through bilayer ionexchange membranes / Ion Transport in Organic and Inorganic Membranes: Conference Proceedings (Sochi, 22-27 May 2023). Krasnodar. pp. 185-187. EDN: WKHJME.
9. Prosviryakov E. Yu. Non-helical exact solutions to the Euler equations for swirling axisym-metric fluid flows// Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019. vol.23, no. 4. pp. 764-770. EDN: WITBIY. DOI: https://doi.org/10.14498/vsgtu1715.
10. Сизых Г. Б. Расщепление уравнений Навье-Стокса для одного класса осесимметрич-ных течений // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2020. Т. 24, №1. С. 163-173. EDN: RUDNVJ. DOI: https://doi.org/10.14498/vsgtu1740.
11. Заболоцкий В. И., Никоненко В. В. Перенос ионов в мембранах. М.: Наука, 1996. 392 с.
12. Узденова А. М., Коваленко А. В., Уртенов М. Х., Никоненко В. В. Теоретический анализ влияния концентрации ионов в объеме раствора и у поверхности мембраны на массоперенос при сверхпредельных токах// Электрохимия, 2017. Т. 53, №11. С. 14211433. EDN: JVMTDS. DOI: https://doi.org/10.7868/S0424857017110032.
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2024, vol. 28, no. 1, pp. 130-151
d https://doi.org/10.14498/vsgtu2043
ISSN: 2310-7081 (online), 1991-8615 (print)
MSC: 80A30, 35Q60, 78A57
Hybrid numerical-analytical method for solving the problems of salt ion transport in membrane systems with axial symmetry
E. V. Kazakovtseva, A. V. Kovalenko,
A. V. Pismenskiy, M. Kh. Urtenov
Kuban State University,
149, Stavropolskaya st., Krasnodar, 350040, Russian Federation.
Abstract
The aim of this study is to develop a new hybrid numerical-analytical method for solving boundary value problems with axial symmetry, for example, with a rotating membrane disk, based on matching the asymptotic solution near the cation-exchange membrane (CEM) with the numerical solution in the rest of the region. For this, the following method is used:
1) a basic mathematical model for the transfer of salt ions in an electrochemical cell with a rotating cation-exchange membrane disk is proposed based on the general conservation laws represented by the Nernst-Planck-Poisson and Navier-Stokes equations with natural boundary and initial conditions. This model contains no fitting parameters or simplifying assumptions. However, the numerical solution of the corresponding boundary value problem presents significant computational difficulties for real solution concentrations and large jumps in the potential and angular velocity of the membrane disk rotation, associated with large concentration and potential gradients near the CEM in the quasi-equilibrium space charge region (SCR);
Mathematical Modeling, Numerical Methods and Software Complexes Research Article
© Authors, 2024
© Samara State Technical University, 2024 (Compilation, Design, and Layout) 3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Kazakovtseva E. V., Kovalenko A. V., Pismenskiy A. V., Urtenov M. Kh. Hybrid numerical-analytical method for solving the problems of salt ion transport in membrane systems with axial symmetry, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2024, vol. 28, no. 1, pp. 130-151. EDN: BOXTTX. DOI: 10.14498/vsgtu2043 (In Russian). Authors' Details:
Ekaterina V. Kazakovtseva https://orcid.org/0009-0003-0040-0880
Senior Lecturer; Dept. of Data Analysis and Artificial Intelligence; e-mail: [email protected]
Anna V. Kovalenko © https://orcid.org/0000-0002-3991-3953
Dr. Techn. Sci., Associate Professor; Head of Department; Dept. of Data Analysis and Artificial Intelligence; e-mail: savanna- [email protected]
Alexander V. Pismenskiy © https://orcid.org/0000-0003-4046-2229 Cand. Phys. & Math. Sci.; Head of Department; Dept. of Applied Mathematics; e-mail: [email protected]
Makhamet Kh. Urtenov © https://orcid.org/0000-0002-0252-6247
Dr. Phys. & Math. Sci., Professor; Dept. of Applied Mathematics; e-mail: [email protected]
2) the solution region is divided into two parts, one of which is a small cation increase region (CIR) located near the CEM, and the remaining main part of the region (MPOR);
3) in the CIR, an analytical solution is found by the method of matching asymptotic solutions;
4) a simplified mathematical model is constructed in the MPOR, which differs from the basic mathematical model in such a boundary condition at the boundary with the CIR, which then allows us the solution of the corresponding boundary value problem to be matched with the solution in the CIR.
The main result is a hybrid numerical-analytical method that allows one to carry out a numerical analysis of the transfer of salt ions at real concentrations of a binary salt electrolyte solution in a wide range of changes in the potential jump and the angular velocity of the membrane disk.
Based on the results of the work, the following conclusion can be drawn, that the combination of the analytical (asymptotic) method of solving in the region of the boundary layer and the numerical solution in the rest of the region, with the exception of the boundary layer, with their subsequent splicing, makes it possible to construct an effective hybrid numerical-analytical method for solving the problems of salt ion transport in membrane systems with axial symmetry.
Keywords: asymptotic solution, hybrid numerical-analytical method, rotating membrane disk, axial symmetry, cylindrical coordinate system.
Received: 22nd July, 2023 / Revised: 12th April, 2024 / Accepted: 15th April, 2024 / First online: 17th May, 2024
Competing interests. We declare that we have no conflict of interest regarding the authorship and publication of this article.
Authors' contributions and responsibilities. E.V. Kazakovtseva: Development of a numerical-analytical method and a program of numerical experiments; Computer and mathematical modeling; Interpretation of the obtained results; Writing — original draft and review & editing. A.V. Kovalenko: Computer modeling; Processing and analysis of experimental data; Writing — original draft. A.V. Pisemskiy: Processing and analysis of experimental data; Interpretation of the obtained results. M.Kh. Urtenov: Idea of study; Mathematical modeling; Development of a numerical-analytical method; Interpretation of the obtained results; preparation of the initial version of the manuscript; Writing — original draft and review & editing. The authors bear full responsibility for providing the final manuscript for publication. The final version of the manuscript was approved by all authors.
Funding. The research was funded by the Russian Science Foundation (project no. 2419-00648), https://rscf.ru/en/project/24-19-00648/.
Acknowledgments. The authors are grateful to the reviewers for the careful reading of the article, valuable suggestions, and comments.
References
1. Burmasheva N. V., Prosviryakov E. Yu. Exact solution of Navier-Stokes equations describing spatially inhomogeneous flows of a rotating fluid, Trudy Instituta Matematiki i Mekhaniki URO RAN, 2020, vol.26, no. 2, pp. 79-87 (In Russian). EDN: IAWMLK. DOI: https://doi. org/10.21538/0134-4889-2020-26-2-79-87.
2. Zabolotskii V. I., Shel'deshov N. V., Sharafan M. V. Electric mass transfer of sodium chloride through cation-exchange membrane MK-40: A rotating membrane disk study, Russ. J. Electrochem., 2006, vol. 42, no. 12, pp. 1345-1351. EDN: LJVHER. DOI: https://doi.org/10. 1134/S1023193506120123.
3. Kazakovtseva E. V. A theoretical study of the quasi-equilibrium region of the space charge in membrane systems with axial symmetry, Perspektivy Nauki, 2023, no. 6(165), pp. 58-68 (In Russian). EDN: DRJSOK.
4. Prosviryakov E. Yu. Recovery of radial-axial velocity in axisymmetric swirling flows of a viscous incompressible fluid in the Lagrangian consideration of vorticity evolution, Vestn. Udmurtsk. Univ. Mat. Mekh. Komp. Nauki, 2021, vol.31, no. 3, pp. 505-516 (In Russian). EDN: ORVWHT. DOI: https://doi.org/10.35634/vm210311.
5. Chubyr N. O., Kovalenko A. V., Urtenov M. Kh., Gudza I. V. Mathematical model of salt ion stationary transport in the cross section of the channel at equilibrium, Modeling, Optimization and Information Technology, 2022, vol. 10, no. 3(38) (In Russian). EDN: DIECUK. DOI: https://doi.org/10.26102/2310-6018/2022.38.3.009.
6. Achoh A., Melnikov S., Bondarev D. Electrochemical characteristics of the MF-4SK membrane doped with the hyperbranched phosphorylated dendrimer BOLTORN H20, In: Ion Transport in Organic and Inorganic Membranes, Conference Proceedings (Sochi, 22-27 May 2023). Krasnodar, 2023, pp. 15-17. EDN: RPXPHA.
7. Bondarev D., Eterevskova S., Zabolotsky V., et al. Homogeneous anion-exchange membrane with heterocyclic functional groups, In: Ion Transport in Organic and Inorganic Membranes, Conference Proceedings (Sochi, 22-27 May 2023). Krasnodar, 2023, pp. 31-32. EDN: BPLDVU.
8. Melnikov S. Experimental and theoretical study of ion transport through bilayer ionexchange membranes, In: Ion Transport in Organic and Inorganic Membranes, Conference Proceedings (Sochi, 22-27 May 2023). Krasnodar, pp. 185-187. EDN: WKHJME.
9. Prosviryakov E. Yu. Non-helical exact solutions to the Euler equations for swirling axisymmetric fluid flows, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2019, vol.23, no. 4, pp. 764-770. EDN: WITBIY. DOI: https://doi.org/10.14498/vsgtu1715.
10. The splitting of Navier-Stokes equations for a class of axisymmetric flows, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2020, vol.24, no. 1, pp. 163-173 (In Russian). EDN: RUDNVJ. DOI: https://doi.org/ 10.14498/vsgtu1740.
11. Zabolotskii V. I., Nikonenko V. V. Perenos ionov v membranakh [Ion Transport in Membranes]. Moscow, Nauka, 1996, 392 pp. (In Russian)
12. Uzdenova A. M., Kovalenko A. V., Urtenov M. K., Nikonenko V. V. Theoretical analysis of the effect of ion concentration in solution bulk and at membrane surface on the mass transfer at overlimiting currents, Russ. J. Electrochem., 2017, vol. 53, no. 11, pp. 1254-1265. EDN: XXDDNZ. DOI: https://doi.org/10.1134/S1023193517110179.
