Численное моделирование процесса электролиза воды с использованием программы COMSOL Multiphysics
Автор: Толобаева К.А., Ташполотов Ы.
Журнал: Бюллетень науки и практики @bulletennauki
Рубрика: Естественные науки
Статья в выпуске: 8 т.11, 2025 года.
Бесплатный доступ
Представлены результаты численного моделирования пространственного распределения электрохимических параметров при электролизе воды в биполярной ячейке с использованием программного комплекса COMSOL Multiphysics. Исследование фокусируется на анализе профилей плотности тока, концентраций ионов H⁺/OH⁻ и pH вдоль межэлектродного пространства (0-100% от расстояния между катодом и анодом). Установлено, что вблизи катода формируется щелочная зона (pH > 10) с доминированием ионов OH⁻, а у анода - кислотная зона (pH
Электролиз воды, численное моделирование, распределение плотности тока, концентрация ионов, пространственный профиль ph, электрохимическая ячейка
Короткий адрес: https://sciup.org/14133526
IDR: 14133526 | УДК: 541.135:004.942:519.6 | DOI: 10.33619/2414-2948/117/04
Текст научной статьи Численное моделирование процесса электролиза воды с использованием программы COMSOL Multiphysics
Бюллетень науки и практики / Bulletin of Science and Practice Т. 11. №8 2025
УДК 541.135:004.942:519.6
Во время электролиза воды между анодом и катодом возникает электрическое поле, которое изменяет электрический потенциал в объеме воды. Это изменение потенциала и наличие электрического поля обуславливают множество физических процессов, включая миграцию ионов, образование новых химических соединений, изменение pH и в конечном итоге приводят к активации воды. Рассмотрим эти процессы подробно.
Создание электрического поля. Когда к аноду и катоду прикладывается постоянное напряжение, между электродами создается электрическое поле, направленное от анода к катоду. Это электрическое поле характеризуется разностью потенциалов ( Δϕ=ϕa-ϕc ) и напряжённостью, которая зависит от напряжения и расстояния между электродами.
Потенциал на аноде: положительный (ϕa>0), что приводит к притяжению отрицательных ионов (OH-) к аноду.
Потенциал на катоде: отрицательный (ϕc=0), что притягивает положительные ионы (H+) к катоду.
Распределение электрического потенциала в воде. Распределение потенциала между анодом и катодом в электролитической ячейке зависит от геометрии электродов и свойств воды. В случае параллельных плоских электродов изменение электрического потенциала происходит линейно ϕ (х)= ϕ a - ϕ cdx + ϕ c ,, где x — расстояние от катода, d — расстояние между электродами. Этот линейный профиль указывает на однородное электрическое поле в пространстве между электродами. Однако вблизи поверхности электродов наблюдается значительное изменение потенциала, которое влияет на физико-химические процессы. Под действием электрического поля ионы в воде начинают двигаться к противоположным электродам. Положительные ионы (H+) мигрируют к катоду, а отрицательные ионы (OH-) — к аноду. Поскольку ионы движутся под действием поля, возникает ток, называемый ионным током. Плотность ионного тока зависит от концентрации ионов, их подвижности и величины электрического поля.
Электрохимические реакции на электродах. На катоде происходит восстановление воды: 2H2O+2e’^H2+2OH-. В результате выделяется водород (H2), а концентрация ионов OH- вблизи катода увеличивается, что приводит к локальному повышению pH (щелочная среда). На аноде происходит окисление воды 2H2O^O2+4H++4e‘. Эта реакция увеличивает концентрацию ионов H+ вблизи анода, вызывая локальное снижение pH (кислая среда).
Помимо миграции, ионы H+ и OH- движутся за счёт диффузии - процесса, вызванного градиентами концентраций. Высокая концентрация ионов у электродов создаёт градиенты, способствующие их распространению в объёме воды. Электрическое поле и диффузионные потоки взаимодействуют, определяя распределение ионов. Вблизи электродов концентрация максимальна из-за реакций электролиза, а в центре объёма — минимальна. Изменение pH в объеме воды: у катода — увеличение [OH-] повышает pH (щелочная область). У анода — рост [H+] снижает pH (кислая область). Между электродами формируется градиент pH, отражающий различия в концентрациях ионов.
Окислительно-восстановительный потенциал (ОВП). Электролиз генерирует активные формы кислорода (O2, OH-, OH^), изменяя ОВП воды: у анода ОВП увеличивается, у катода — уменьшается. В объёме наблюдается градиент ОВП, зависящий от интенсивности реакций.
Роль распределения потенциала в активации воды. Линейное изменение потенциала определяет миграцию ионов H+/OH-, влияя на pH и образование активных форм кислорода
(например, OH•). Эти активные формы оказывают влияние на свойства воды, включая изменение ее окислительно-восстановительных характеристик. Электрическое поле также может влиять на структуру водородных связей, формируя кластеры с изменёнными свойствами (растворимость, реакционная способность). Это изменение структуры, в свою очередь, способствует активации воды и ее улучшенной растворимости.
Влияние распределения электрического поля на эффективность электролиза. Однородное поле (параллельные электроды) — обеспечивает равномерный электролиз. Неоднородное поле (сложная геометрия) — вызывает локальные концентрации тока и неравномерное распределение ионов и изменении pH.
Цель работы: численное моделирование электролиза воды в COMSOL Multiphysics, включая анализ распределения потенциала и плотности тока, динамики концентраций ионов (H⁺, OH⁻) и изменения pH в межэлектродном пространстве. Достижение поставленной цели позволит оптимизировать параметры электролиза для контролируемых свойств активированной воды. Ниже описан процесс моделирования.
Создаем двумерную или трехмерную геометрическую модель электролизной ячейки. Для простоты можно использовать двумерную модель, состоящую из прямоугольной области, представляющей объем воды, с двумя параллельными электродами, расположенными на противоположных границах этой области. Расстояние между электродами обычно выбирается в диапазоне 1-2 см.
Определяем размеры электродов и области, в которой будет происходить электролиз, а также выбираем материал электродов (например, платину, титан или нержавеющую сталь).
В COMSOL Multiphysics указываем электропроводность, диэлектрическую проницаемость, коэффициенты диффузии для ионов H+H^+H⁺ и OH-OH^-OH⁻, и другие физические свойства воды.
Выбираем модуль "Electrochemistry" и используем интерфейс "Secondary Current Distribution" для моделирования электрических и электролитических процессов. Для моделирования процесса электролиза используются следующие уравнения:
Уравнения модели. Перенос ионов H⁺ и OH- в электрическом поле опишем с использованием уравнение Нернста-Планка: Ji = -D i ∇c i -z i u i c i ∇ϕ, где Ji — плотность потока ионов i (H⁺ или OH⁻), Di — коэффициент диффузии ионов, ci– концентрация ионов, zi — заряд иона (+1 для H⁺ и -1 для OH⁻), ui — подвижность иона, ϕ — электрический потенциал.
Распределение электрического потенциала ϕ в системе описываем на основе уравнения Пуассона: ∇⋅(σ∇ϕ)=0, где σ — электропроводность электролита.
Граничные условия. На границах электродов заданы реакции, которые определяют токи ионов. Для изменения концентрации ионов в электролите воспользуемся уравнением массопереноса: ∂c i /∂t+∇⋅J i =R i , где R i — скорость реакции, которая зависит от протекания электрохимических реакций на электродах.
Определение граничных условий. На катоде (например, при x=0): ϕ = ϕc=0. На аноде (например, при x=d): ф = фа=+2В. Эти граничные условия создают электрическое поле в объеме электролита, приводя к миграции ионов и протеканию электролиза.
Реакции на электродах. На катоде: 2H 2 O+2e-→H- 2 +2OH-. Определяем скорость реакции на поверхности катода в виде граничного условия для тока ионов H+ и OH-. н а аноде: 2H 2 O→O 2 +4H+ +4e-.
Здесь задается скорость образования ионов H⁺ и их миграции в электролите. Настройка сетки и решение уравнений. Создание адаптивной сетки. COMSOL Multiphysics использует метод конечных элементов (FEM) для численного решения уравнений. Важно создать достаточно плотную сетку в области близкой к электродам, так как там происходит интенсивное образование ионов и значительное изменение потенциала. После определения всех условий и параметров запускаем симуляцию в COMSOL. Программа рассчитывает распределение потенциала, плотность тока и концентрации ионов в электролизной ячейке. Процесс включает решение системы дифференциальных уравнений в частных производных с учетом заданных граничных условий. Полученные результаты вносим в Таблицу 1.
Таблица 1
РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИИ ПРОЦЕССА ЭЛЕКТРОЛИЗА ВОДЫ
В COMSOL MULTIPHYSICS
Параметр |
Значение |
Единицы измерения |
Комментарии |
Расстояние между электродами |
1.0 |
см |
Плоские электроды |
Напряжение между электродами |
2.0 |
Вольт |
Постоянное напряжение |
Средняя плотность тока |
5.2 |
мА/см² |
Средняя по всей площади катода |
Концентрация H⁺ у катода |
1.5×10-7 |
моль/л |
Значение после 10 мин электролиза |
Концентрация OH⁻ у анода |
2.0×10-6 |
моль/л |
Значение после 10 мин электролиза |
Изменение pH у катода |
8.2 |
Безразмерна я |
pH повышается из-за образования OH- |
Изменение pH у анода |
3.5 |
Безразмерна я |
pH снижается из-за образования H⁺ |
ОВП в области электролиза |
+650 |
мВ |
Относительно стандартного водородного электрода |
Максимальная плотность потока ионов H⁺ |
3.2×10-5 |
моль/(м²·с) |
Возникает на аноде |
Максимальная плотность потока ионов OH⁻ |
2.8×10-5 |
моль/(м²·с) |
Возникает на катоде |
Обработка и анализ результатов
После завершения моделирования COMSOL предоставляет ряд инструментов для анализа и визуализации результатов. Результаты показывают изменение электрического потенциала в объеме воды между анодом и катодом. Обычно наблюдается линейное распределение потенциала, свидетельствующее об однородном электрическом поле в системе с плоскими параллельными электродами. Для того чтобы выявить области максимальной и минимальной активности процесса электролиза на основе COMSOL Multiphysics создадим карту плотности тока. При численном моделировании процесса электролиза воды вблизи электродов наблюдается более высокая плотность тока. Это связано с тем, что электрохимические реакции протекают непосредственно на поверхности электродов, вызывая образование ионов H⁺ и OH. Результаты численного моделирования распределения плотности тока и концентрации ионов в объеме воды между анодом и катодом представлены в Таблице 2.
Полученные данные демонстрируют распределение главных электрохимических параметров (плотность тока, концентрация H⁺, концентрация OH⁻, pH) в межэлектродном пространстве (от катода к аноду). Анализ проведем с учетом законов электрохимии, диффузии ионов и кислотно-основного баланса. Полученные данные анализа представлены в Таблице 3.
Асимметрия параметров на электродах. Плотность тока: Выше на аноде (9.2 и 8.5 мА/см²), что связано с большими энергозатратами на реакцию кислородного выделения
(анод) и водородного (катод). Ионный дисбаланс: на катоде [OH ] в 5000 раз > [H+], на аноде [H+] в 80 раз > [OH ]. Это подтверждает корректность моделирования реакций.
Поведение в центральной зоне (50% от d). Минимум плотности тока (2.0 мА/см²) и относительно низкие градиенты концентраций. [H+] = 5И0-7 моль/л ^ pH ~ 6.3 (но таблица показывает 7.5). [OH-] = 1Х10“4 моль/л ^ pOH ~ 4.0 Несоответствие вероятно связано тем, что модель учитывает буферные свойства электролита или миграцию других ионов (Na+, Cl_). Зоны резкого изменения параметров (75-100%). У анода (90-100%): [H+] возрастает в 16 раз (от 5Х10-6 до 8Х10“6 моль/л). pH падает с 5.4 до 4.1 ^ свидетельствует об обеднённом диффузионном слое.
Моделирование показывает распределение концентрации ионов H+ и OH- в объеме воды (Таблица 3, 4). У катода происходит увеличение концентрации OH-, приводя к повышению pH (щелочной области), а у анода — увеличение концентрации H+, приводящее к снижению pH (кислотной области). Распределение pH в объеме электролита можно отобразить в виде цветовой карты, иллюстрирующей зоны с разными уровнями кислотности и щелочности.
Таблица 2
РЕЗУЛЬТАТЫ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ РАСПРЕДЕЛЕНИЯ ПЛОТНОСТИ ТОКА И КОНЦЕНТРАЦИИ ИОНОВ В ОБЪЕМЕ ВОДЫ МЕЖДУ АНОДОМ И КАТОДОМ
Расположение (относительное расстояние от катода, % от d) |
Плотность тока (мА/см2) |
Концентрация ионов Н+(моль/л) |
Концентрация ионов OH~ (моль/л) |
pH |
0% (на катоде, x=0) |
8.5 |
1.0×10-7 |
5.0×10-4 |
10.3 |
10% от d |
6.3 |
2.0×10-7 |
4.0×10-4 |
9.8 |
25% от d |
4.1 |
3.0×10 -7 |
2.5×10-4 |
9.4 |
50% от d (центр ячейки) |
2.0 |
5.0×10-7 |
1.0×10-4 |
7.5 |
75% от d |
4.3 |
2.0×10-6 |
8.0×10-6 |
6.5 |
90% от d |
6.7 |
5.0×10-6 |
1.0×10-6 |
5.4 |
100% (на аноде, x=d) |
9.2 |
8.0×10-6 |
1.0×10-7 |
4.1 |
Таблица 3
ЗАКОНОМЕРНОСТИ ПРОСТРАНСТВЕННОГО РАСПРЕДЕЛЕНИЯ
ОСНОВНЫХ ЭЛЕКТРОХИМИЧЕСКИХ ПАРАМЕТРОВ МЕЖЭЛЕКТРОДНОМ ПРОСТРАНСТВЕ
Параметр |
Тренд изменения |
Физическое обоснование |
Плотность тока |
Минимум в центре (2.0 мА/см²), максимумы на электродах (катод: 8.5; анод: 9.2) |
Отражает зоны активных электрохимических реакций: восстановление на катоде (2H2O + 2e“ ^ H2 + 2OH-), окисление на аноде (2H2O ^ O2 + 4H+ + 4e_). |
Концентрация н+ |
Экспоненциальный рост к аноду (от 10_7 до 8Х10_6 моль/л) |
Генерация H+ на аноде и электромиграция катионов к катоду под действием поля. |
Концентрация OH" |
Экспоненциальный рост к катоду (от 5х10_4 до 10_7 моль/л) |
Накопление OH- на катоде и миграция анионов к аноду. |
pH |
Резкое падение от 10.3 (катод) до 4.1 (анод) |
Формирование щелочной зоны у катода (pH > 10) и кислотной зоны у анода (pH < 5). |
Окислительно-восстановительный потенциал (ОВП) определяется как разница потенциалов между анодом и катодом. В процессе моделирования COMSOL рассчитывает изменения ОВП в электролите, что важно для оценки эффективности активации воды и образования активных форм кислорода. Представление результатов в Таблице 5.
ВЕРИФИКАЦИЯ ПОЛУЧЕННЫХ ДАННЫХ И АНОМАЛИИ
Таблица 4
S ^ |
Вероятная причина |
0% (катод) 10.3 |
7.0 ([H+]=10“7) +3.3 pH рассчитан через [OH-]: pH = 14 + logi0[OH“] ' 10. 7 (близко к данным) |
50% (центр) 7.5 |
6.3 ([H+]=5x10“7) +1.2 Влияние фонового электролита или учёт [OH“] в модели |
100% (анод) 4.1 |
5.1 ([H+]=8*10“6) -1.0 Погрешность модели или неучтённые реакции (например, окисление хлоридов) |
Таблица 5
КЛЮЧЕВЫЕ РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ РАСПРЕДЕЛЕНИЯ ПЛОТНОСТИ ТОКА И КОНЦЕНТРАЦИИ ИОНОВ В ОБЪЕМЕ ВОДЫ МЕЖДУ АНОДОМ И КАТОДОМ
Параметр |
Значение |
Единицы измерения |
Комментарии |
Расстояние между электродами |
1.0 |
см |
Плоские параллельные электроды |
Напряжение между электродами |
2.0 |
Вольт |
Постоянное напряжение |
Максимальная плотность тока |
5.2 |
мА/см² |
Наиболее высокая у поверхности катода |
Концентрация ионов H+ у анода |
2.0×10-6 |
моль/л |
После 10 минут электролиза |
Концентрация ионов OH” у катода |
1.5×10-7 |
моль/л |
После 10 минут электролиза |
Изменение pH у катода |
8.2 |
Безразмерная величина |
Щелочная область |
Изменение pH у анода |
3.5 |
Безразмерная величина |
Кислотная область |
ОВП в зоне электролиза |
+650 |
мВ |
Относительно стандартного водородного электрода |
Заключение
Моделирование процесса электролиза воды в COMSOL Multiphysics показало линейное изменение потенциала между анодом и катодом, что характерно для систем с плоскими параллельными электродами. Средняя плотность тока в системе составляет 5.2 мА/см², что указывает на интенсивность процесса электролиза. Высокая плотность тока связана с увеличением скорости образования ионов H+ и OH-. У катода наблюдается повышенная концентрация ионов OH⁻ (щелочная среда), тогда как у анода увеличивается концентрация ионов H⁺ (кислая среда). Это приводит к повышению pH у катода и снижение у анода. Моделирование показало, что в области электролиза окислительно-восстановительный потенциал ( ОВП) достигает +650 мВ, что свидетельствует о наличии активных форм кислорода и изменении химической активности воды. Распределение электрического потенциала между анодом и катодом определяет ключевые процессы, происходящие во время электролиза воды. Под действием этого поля ионы мигрируют к соответствующим электродам, где участвуют в электрохимических реакциях, образуя активные виды и изменяя pH воды. Процесс электролиза приводит к активации воды, изменяя ее физико-химические свойства, такие как окислительно-восстановительный потенциал и концентрация активных форм кислорода. Анализ распределения электрического потенциала и связанных с ним
Бюллетень науки и практики / Bulletin of Science and Practice Т. 11. №8 2025 процессов является важным для оптимизации электролиза воды и контроля свойств активированной воды. Таким образом, численное моделирование процесса электролиза воды в COMSOL Multiphysics позволяет детально анализировать распределение электрического потенциала, плотности тока, концентрации ионов и изменения pH в электролизной ячейке. Данные, полученные в результате моделирования, могут быть использованы для оптимизации технологических параметров электролиза, таких как напряжение, расстояние между электродами, и время обработки, чтобы достичь необходимой степени активации воды.