Массоперенос в электромембранных системах в режимах интенсивного тока сопровождается возникновением дополнительных механизмов переноса, которые существенно влияют на эффективность их функционирования. Согласно современным представлениям для разбавленных растворов электролитов среди таких механизмов особенно важными являются электроконвекция и реакции диссоциации/рекомбинации молекул воды. Эти процессы оказывают противоположное действие на эффективность электромембранных технологий.
В исследованиях мембранных систем активно применяются математические модели, учитывающие влияние указанных механизмов, однако они обычно описывают только потенциодинамический режим, при котором устанавливается скачок потенциала в системе. Для интерпретации обширной базы экспериментальных данных по гальванодинамическому режиму (при фиксированной плотности тока) также необходимы инструменты теоретического анализа.
Цель данной работы заключается в разработке математической модели массопереноса в слое раствора электролита у ионообменной мембраны с учетом электроконвекции и диссоциации воды в гальванодинамическом режиме. Модель основана на системе связанных уравнений Нернста—Планка—Пуассона—Навье—Стокса, дополненной новым гальванодинамическим граничным условием для потенциала.
С использованием разработанной модели впервые рассчитаны хронопотенциограммы мембранной системы с учетом влияния как электроконвекции, так и реакции диссоциации/рекомбинации молекул воды. Результаты показали, что отношение концентрации продуктов диссоциации воды к концентрации ионов соли определяет баланс эффектов электроконвекции и диссоциации.
Рассмотрены следующие варианты соотношения эффектов электроконвекции и диссоциации молекул воды:
значимое влияние на массоперенос оказывает электроконвекция, в то время как влияние диссоциации воды минимально;
электроконвекция и диссоциация существенно влияют на процессы переноса:
образование дополнительных носителей заряда в результате диссоциации молекул воды снижает скачок потенциала в слое раствора электролита, что уменьшает интенсивность электроконвекции, в то время как развитие электроконвекции, в свою очередь, замедляет процесс диссоциации;
продукты интенсивной диссоциации молекул воды тормозят развитие электроконвекции.
Идентификаторы и классификаторы
Электромембранные системы (электродиализные аппараты, нано- и микрофлюидные устройства) применяются для обессоливания и деионизации воды, производства кислот и оснований, при преобразовании и производстве энергии и в других областях [1–3]. В основе принципов функционирования электромембранных систем лежит селективный перенос ионообменных мембран, то есть проницаемость для ионов одного знака (противоионов) и ограничение движения ионов другого знака (коионов). В структуре вольт-амперной кривой электромембранной системы выделяются следующие характерные участки [4, 5]:
Список литературы
- Ran J., Wu L., He Y., et al. Ion exchange membranes: New developments and applications // J. Membr. Sci., 2017. vol. 522. pp. 267–291. DOI: https://doi.org/10.1016/j.memsci.2016.09.033.
- Slouka Z., Senapati S., Chang H. C. Microfluidic systems with ion-selective membranes // Annu. Rev. Anal. Chem., 2014. vol. 7. pp. 317–335. DOI: https://doi.org/10.1146/annurev-anchem-071213-020155.
- Gurreri L., Tamburini A., Cipollina A., Micale G. Electrodialysis applications in wastewater treatment for environmental protection and resources recovery: A systematic review on progress and perspectives // Membranes, 2020. vol. 10, no. 7, 146. DOI: https://doi.org/10.3390/membranes10070146.
- Рубинштейн И., Зальцман Б., Прец И., Линдер К. Экспериментальная проверка электроосмотического механизма формирования “запредельного” тока в системе с катионообменной электродиализной мембраной // Электрохимия, 2002. Т. 38, №8. С. 956–967. EDN: RSNSYP.
- Письменская Н. Д., Никоненко В. В., Белова Е. И. [и др.] Сопряженная конвекция раствора у поверхности ионообменных мембран в режимах интенсивного тока // Электрохимия, 2007. Т. 43, №3. С. 325–345. EDN: IACEHN.
- Никоненко В. В., Мареев С. А., Письменская Н. Д. [и др.] Эффект электроконвекции и его использование для интенсификации массопереноса в электродиализе (обзор) // Электрохимия, 2014. Т. 53, №10. С. 1266–1289. EDN: ZNAASD DOI: https://doi.org/10.7868/S0424857017100061.
- Mani A., Wang K. M. Electroconvection near electrochemical interfaces: Experiments, modeling, and computation // Annu. Rev. Fluid Mech., 2020. vol. 52. pp. 509–529. DOI: https://doi.org/10.1146/annurev-fluid-010719-060358.
- Simons R. Strong electric field effects on proton transfer between membranebound amines and water // Nature, 1979. vol. 280. pp. 824–826. DOI: https://doi.org/10.1038/280824a0.
- Frilette V. J. Electrogravitational transport at synthetic ion exchange membrane surfaces // J. Phys. Chem., 1957. vol. 61, no. 2. pp. 168–174. DOI: https://doi.org/10.1021/J150548A010.
- Заболоцкий В. И., Никоненко В. В., Корженко Н. М. [и др.] Влияние гетеролитической диссоциации воды на массоперенос ионов соли в электромембранной системе при нарушении электронейтральности в области диффузионного слоя // Электрохимия, 2002. Т. 38, №8. С. 911–920. EDN: GIGGKH.
- Узденова А. М. Математическое моделирование нестационарного переноса ионов в электро-мембранных системах с учетом реакции диссоциации (рекомбинации) молекул воды в гальванодинамическом режиме // Перспективы науки, 2023. Т. 11, №170. С. 104–112. EDN: USIXRO.
- Mishchuk N. A. Concentration polarization of interface and non-linear electrokinetic phenomena // Adv. Colloid Interface Sci., 2010. vol. 160, no. 1–2. pp. 16–39. EDN: MYAMRJ. DOI: https://doi.org/10.1016/j.cis.2010.07.001.
- Porozhnyy M. V., Shkirskaya S. A., Butylskii D. Y., et al. Physicochemical and electrochemical characterization of nafion-type membranes with embedded silica nanoparticles: effect of functionalization // Electrochim. Acta, 2021. vol. 370, 137689. EDN: INLSZS. DOI: https://doi.org/10.1016/j.electacta.2020.137689.
- Grossman G. Water dissociation effects in ion transport through composite membrane // J. Phys. Chem., 1976. vol. 80, no. 14. pp. 1616–1625. DOI: https://doi.org/10.1021/j100555a020.
- Rubinstein I. A diffusional model of “water splitting” in electrodialysis // J. Phys. Chem., 1977. vol. 81, no. 14. pp. 1431–1436. DOI: https://doi.org/10.1021/j100529a018.
- Rubinstein I., Shtilman L. Voltage against current curves of cation exchange membranes // J. Chem. Soc., Faraday Trans. 2, 1979. vol. 75. pp. 231–246. DOI: https://doi.org/10.1039/F29797500231.
- Коваленко А. В., Уртенов М. Х., Сеидова Н. М., Письменский А. В. Влияние реакции диссоциации/рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 1. Математическая модель // Научный журнал КубГАУ, 2016. Т. 121, 122. EDN: WWSKMP. DOI: https://doi.org/10.21515/1990-4665-121-122.
- Коваленко А. В., Уртенов М. Х., Сеидова Н. М., Письменский А. В. Влияние реакции диссоциации/рекомбинации молекул воды на перенос 1:1 электролита в мембранных системах в диффузионном слое. Часть 2. Асимптотический анализ // Научный журнал КубГАУ, 2016. №122, 017. EDN: XBDYXZ. DOI: https://doi.org/10.21515/1990-4665-122-017.
- Уртенов М. Х., Письменский А. В., Никоненко В. В., Коваленко А. В. Математическое моделирование переноса ионов и диссоциации воды у границы ионообменная мембрана/раствор в интенсивных токовых режимах // Мембр. и мембр. техн., 2018. Т. 8, №1. С. 24–33. EDN: YNIJEE. DOI: https://doi.org/10.1134/S2218117218010054.
- Urtenov M., Gudza V., Shkorkina I., Chubyr N. Theoretical analysis of the stationary transport of 1:1 salt ions in a cross-section of a desalination channel, taking into account the non-catalytic dissociation/recombination reaction of water molecules // Membranes, 2020. vol. 10, no. 11, 342. DOI: https://doi.org/10.3390/membranes10110342.
- Kovalenko A. V., Nikonenko V. V., Chubyr N. O., Urtenov M. Kh. Mathematical modeling of electrodialysis of a dilute solution with accounting for water dissociation-recombination reactions // Desalination, 2023. vol. 550, 116398. DOI: https://doi.org/10.1016/j.desal.2023.116398.
- Коваленко А. В., Уртенов М. Х., Чубырь Н. О. [и др.] Математическое моделирование влияния основных температурных эффектов на стационарный перенос ионов соли в диффузионном слое // Экол. вестн. научн. центров ЧЭС, 2018. Т. 15, №3. С. 78–86. EDN: YABSYX. DOI: https://doi.org/10.31429/vestnik-15-3-78-86.
- Коваленко А. В., Уртенов М. Х., Чубырь Н. О. [и др.] Влияние температурных эффектов, связанных с реакцией диссоциации/рекомбинации молекул воды и джоулевым нагревом раствора на стационарный перенос ионов соли в диффузионном слое // Экол. вестн. научн. центров ЧЭС, 2018. Т. 15, №4. С. 67–84. EDN: YRMLET. DOI: https://doi.org/10.31429/vestnik-15-4-67-84.
- Nikonenko V., Urtenov M., Mareev S., Pourcelly G. Mathematical modeling of the effect of water splitting on ion transfer in the depleted diffusion layer near an ionexchange membrane // Membranes, 2020. vol. 10, no. 2, 22. DOI: https://doi.org/10.3390/membranes10020022.
- Коваленко А. В. Влияние диссоциации воды на развитие электроконвекции в мембранных системах // Конденсированные среды и межфазные границы, 2014. Т. 16, №3. С. 288–293. EDN: SQBSBF.
- Kovalenko A., Urtenov M., Chekanov V. Kandaurova N. Theoretical analysis of the influence of spacers on salt ion transport in electromembrane systems considering the main coupled effects // Membranes, 2024. vol. 14, no. 1, 20. DOI: https://doi.org/10.3390/membranes14010020.
- Newman J., Thomas-Alyea K. E. Electrochemical Systems. NJ, USA: John Wiley and Sons, 2004. xx+647 pp.
- Uzdenova A. M. 2D mathematical modelling of overlimiting transfer enhanced by electroconvection in flow-through electrodialysis membrane cells in galvanodynamic mode // Membranes, 2019. vol. 9, no. 3, 39. DOI: https://doi.org/10.3390/membranes9030039.
- Uzdenova A. M. Time-dependent two-dimensional model of overlimiting mass transfer in electromembrane systems based on the Nernst–Planck, displacement current and Navier–Stokes equations // Computation, 2023. vol. 11, no. 10, 205. DOI: https://doi.org/10.3390/computation11100205.
- Cohen H., Cooley J. W. The numerical solution of the time-dependent Nernst–Planck equations // Biophys. J., 1965. vol. 5, no. 2. pp. 145–162. DOI: https://doi.org/10.1016/s0006-3495(65)86707-8.
- Brumleve T. R., Buck R. P. Numerical solution of the Nernst–Planck and Poisson equation system with applications to membrane electrochemistry and solid state physics // J. Electroanal. Chem., 1978. vol. 90, no. 1. pp. 1–31. DOI: https://doi.org/10.1016/s0022-0728(78)80137-5.
- COMSOL Multiphysics Reference Manual. https://doc.comsol.com/6.1/doc/com.comsol.help.comsol/COMSOL_ReferenceManual.pdf.
- Nikonenko V. V., Vasil’eva V. I., Akberova E. M., et al. Competition between diffusion and electroconvection at an ion-selective surface in intensive current regimes // Adv. Colloid Interface Sci., 2016. vol. 235. pp. 233–246. DOI: https://doi.org/10.1016/j.cis.2016.06.014.
- de Valenc˛a J. C., Wagterveld R. M., Lammertink R. G. H., Tsai P. A. Dynamics of microvortices induced by ion concentration polarization // Phys. Rev. E, 2015. vol. 92, no. 3. pp. 031003. DOI: https://doi.org/10.1103/physreve.92.031003.
- Filippov A. N., Akberova E. M., Vasil’eva V. I. Study of the thermochemical effect on the transport and structural characteristics of heterogeneous ion-exchange membranes by combining the cell model and the fine-porous membrane model // Polymers, 2023. vol. 15, no. 16, 3390. DOI: https://doi.org/10.3390/polym15163390.
- Urtenov M. A. K., Kirillova E. V., Seidova N. M., Nikonenko V. V. Decoupling of the Nernst–Planck and Poisson equations. Application to a membrane system at overlimiting currents // J. Phys. Chem. B, 2007. vol. 111, no. 51. pp. 14208–14222. DOI: https://doi.org/10.1021/jp073103d.
- Krol J. J., Wessling M., Strathmann H. Chronopotentiometry and overlimiting ion transport through monopolar ion exchange membranes // J. Membr. Sci., 1999. vol. 162, no. 1–2. pp. 155–164. DOI: https://doi.org/10.1016/S0376-7388(99)00134-9.
Выпуск
Другие статьи выпуска
Рассматриваются математические модели конвекции–диффузии–реакции, которые относятся к моделям тепломассопереноса и применяются при исследовании природных и техногенных процессов. Для данного класса моделей актуальной является задача идентификации как параметров самой модели, так и входящих в нее граничных условий по результатам измерений значений искомой функции в отдельных точках рассматриваемой области. Задачу усложняет наличие неполных измерений, искаженных случайными помехами.
Решение заключается в разработке комбинированного двухэтапного метода идентификации, основанного на последовательном применении метода минимизации критерия идентификации безградиентного типа и рекуррентного метода оценивания неизвестных входных сигналов. Для применения указанных методов выполняется переход от исходной модели, описываемой уравнениями в частных производных, к дискретной линейной стохастической модели в пространстве состояний, в которой неизвестные граничные условия рассматриваются как неизвестные входные сигналы.
В результате построены новые дискретные линейные стохастические модели конвекции–диффузии–реакции для трех разных типов граничных условий. Предложена общая схема процесса параметрической идентификации, включающая двухэтапную идентификацию неизвестных параметров математической модели и идентификацию неизвестных граничных условий.
Для проверки работоспособности предложенного метода построены компьютерные модели конвекции–диффузии–реакции и выполнена реализация всех алгоритмов на языке MATLAB. Проведена серия вычислительных экспериментов, результаты которых показали, что разработанная двухэтапная комбинированная схема позволяет идентифицировать параметры исходной модели, значения функций, входящих в граничные условия, а также вычислить по неполным зашумленным измерениям оценки функции, описывающей процесс конвекции–диффузии–реакции.
Полученные результаты могут быть использованы не только при исследовании процессов тепломассопереноса, но также при решении задач идентификации параметров моделей дискретных стохастических систем с неизвестными входными сигналами и при наличии случайных помех.
Вопрос исследования гравитационного поля тел сложной формы (не относящийся к шарообразной) представляет большой интерес для геофизики, астрофизики, математической физики и других областей. Статья состоит из двух частей. В первой части представлен краткий литературный обзор различных методов расчета потенциала гравитационного поля однородного куба в рамках классической механики: получение аналитического решения; как частный случай задачи нахождения гравитационного поля полиэдра; методом конечных элементов; методом мультипольного разложения. Более подробно проанализирован метод расчета потенциала гравитационного поля однородного куба с помощью аналитического решения и мультипольного разложения. Во второй части статьи описан релятивистский случай гравитационного поля однородного куба в рамках постньютоновского формализма в первом и втором приближении. Данный метод расчета выбран по причине чрезвычайной сложности получение решения с помощью уравнений Эйнштейна. Ранее подобные задачи для тел с формой куба не рассматривались. Для решения задачи выбрана физическая модель — координатный равновесный куб, заполненный несжимаемой жидкостью с нулевой скоростью и постоянной плотностью. Получены релятивистские поправки для временной и пространственной координаты. Получен точный аналитический вид этих поправок для области вне куба, а также компоненты метрического тензора. Дано краткое сравнение полученных результатов для релятивистского случая с результатами классического ньютоновского случая. Для области внутри куба решение получено с помощью численных методов. Полученные результаты с достаточной точностью определяют параметры гравитационного поля для однородного куба, рассмотренного в рамках релятивистского подхода. Основное приложение этой задачи в рамках релятивистской физики относится к области математической физики (или, шире, математики).
Работа посвящена дальнейшему исследованию и построению конструктивных методов последовательной параметрической оптимизации неизвестных характеристик нестационарных процессов технологической теплофизики на компактном множестве непрерывных и непрерывно-дифференцируемых функций. Предложенная методика распространяет разработанный алгоритмически точный метод решения на многомерную постановку обратных задач технологической теплофизики и позволяет отыскать физически обоснованную идентифицируемую характеристику на последовательно сходящихся компактных множествах.
В качестве объекта исследования рассматривается двумерное осесимметричное тело канонической формы, где искомой функцией является сосредоточенная величина мощности внутренних теплоисточников. Задача сформулирована в равномерной метрике оценивания температурного отклонения расчетного состояния от экспериментального. В качестве математической модели рассматриваемого объекта используется его модальное описание, на основе которого проведена редукция исходной обратной задачи теплопроводности, сформулированной в экстремальной постановке, к задаче оптимального управления.
Использование предварительной параметризации искомой характеристики процесса приводит к ее представлению в форме кусочно-параболических функций, конкретизируемых в рамках выбранной структуры с помощью вектора параметров. Количество учитываемых параметров определяет точное представление идентифицируемой величины, а их значения отыскиваются в результате решения полученной задачи параметрической оптимизации. Для решения полученной задачи математического программирования относительно оптимальных значений вектора параметров используются, аналогично одномерному случаю, альтернансные свойства искомых экстремалей, в результате чего задача сводится к замкнутой системе соотношений.
Полученные результаты демонстрируют эффективность распространения конструктивного метода последовательной параметрической оптимизации, опробованного на одномерных обратных задачах теплопроводности, на решение двумерных задач с использованием их модального представления. Увеличение числа параметров решений, образующих кусочно-параболическую форму искомой зависимости, приводит к уменьшению погрешности восстановления как искомой сосредоточенной функции, так и пространственно-временного температурного поля во всей двумерной области определения пространственных переменных.
Представлен численный метод расчета полей остаточных напряжений в поверхностно упрочненном призматическом образце с несквозной V-образной трещиной, базирующийся на упругопластическом решении задачи. По полученным результатам проведен подробный анализ распределений остаточных напряжений вблизи дефекта по нескольким контурам. Определено, что при глубине трещины 0.3 мм практически все изучаемые компоненты остаточных напряжений сжатия имеют бóльшие (по модулю) значения, чем при глубине 0.1 мм, либо равные значения.
Предлагается новый алгоритм нахождения приближенных симметрий для дробно дифференциальных уравнений с производными типа Римана–Лиувилля и Герасимова–Капуто, порядок которых близок к целому. Алгоритм основан на разложении дробной производной в ряд по малому параметру, выделяемому из порядка дробного дифференцирования. В линейном приближении такое разложение содержит нелокальный интегро-дифференциальный оператор с логарифмическим ядром.
В результате исходное дробно-дифференциальное уравнение приближается интегро-дифференциальным уравнением с малым параметром, для которого могут быть найдены приближенные симметрии. Доказывается теорема о виде продолжения однопараметрической группы точечных преобразований на новую переменную, порождаемую нелокальным оператором, входящим в разложение дробной производной. Знание такого продолжения позволяет применить к рассматриваемому уравнению приближенный критерий инвариантности.
Предлагаемый алгоритм иллюстрируется на задаче нахождения приближенных симметрий для нелинейного дробно-дифференциального уравнения фильтрации субдиффузионного типа. Показано, что размерность алгебры приближенных симметрий такого уравнения оказывается существенно больше размерности алгебры точных симметрий, что открывает возможность построения большого числа приближенно инвариантных решений. Также на примере линейного дробно-дифференциального уравнения субдиффузии показывается, что алгоритм дает принципиальную возможность находить нелокальные приближенные симметрии определенного вида.
Финансовая система является важной составляющей в регулировании глобальных экономических процессов, поскольку обеспечение безопасности или контроль финансовой системы или рынка является ключом к стабилизации экономики.
Целью данного исследования является выяснение, насколько приближенные аналитические решения, полученные с помощью метода остаточного степенного ряда и метода разложения Эльзаки для дробной нелинейной финансовой модели, соответствуют экономической теории. Здесь понятие дробной производной используется в смысле производной Капуто.
Полученные численные результаты показывают, как приближенные решения реагируют на изменения процентной ставки, инвестиционного спроса и индекса цен. Оба метода показали результаты, согласующиеся с экономической теорией. Это означает, что исследователи могут использовать эти два метода для решения различных задач, связанных с дробными нелинейными моделями в финансовых системах.
Предлагается новый гибридный численный метод с использованием производной Капуто для решения нелинейного дробного уравнения Льенара — метод дифференциального преобразования Халуты. Доказана теорема сходимости данного метода при определенных условиях.
Метод дифференциального преобразования Халуты представляет собой полуаналитическую технику, объединяющую два мощных подхода: метод преобразования Халуты и метод дифференциального преобразования. Основное преимущество этого метода заключается в том, что он позволяет очень быстро находить решения и не требует линеаризации, возмущения или каких-либо других предположений. Предложенный метод подробно описан, а его эффективность продемонстрирована на двух числовых примерах. Результаты вычислений хорошо согласуются с точными решениями, что подтверждает надежность и эффективность предложенного подхода.
Издательство
- Издательство
- СамГТУ
- Регион
- Россия, Самара
- Почтовый адрес
- 443100, Самарская обл, г Самара, Октябрьский р-н, ул Молодогвардейская, д 244
- Юр. адрес
- 443100, Самарская обл, г Самара, Октябрьский р-н, ул Молодогвардейская, д 244
- ФИО
- Быков Дмитрий Евгеньевич (РЕКТОР)
- E-mail адрес
- rector@samgtu.ru
- Контактный телефон
- +7 (846) 2784300
- Сайт
- https://samgtu.ru/