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. H. Urtenov

Journal of Samara State Technical University, Ser. Physical and Mathematical Sciences, Tome 28 (2024) no. 1, pp. 130-151.

Résumé

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); 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.

Détail

Keywords: asymptotic solution, hybrid numerical-analytical method, rotating membrane disk, axial symmetry, cylindrical coordinate system

Bibliographie

[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, 26:2 (2020), 79–87 (In Russian) | DOI

[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., 42:12 (2006), 1345-1351 | DOI

[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), 58–68 (In Russian)

[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, 31:3 (2021), 505–516 (In Russian) | DOI

[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, 10:3(38) (2022) (In Russian) | DOI

[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, 15–17

[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, 31–32

[8] Melnikov S., “Experimental and theoretical study of ion transport through bilayer ion-exchange membranes”, Ion Transport in Organic and Inorganic Membranes, Conference Proceedings (Sochi, 22–27 May 2023), Krasnodar, 185–187

[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.], 23:4 (2019), 764–770 | DOI

[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.], 24:1 (2020), 163–173 (In Russian) | DOI

[11] Zabolotskii V. I., Nikonenko V. V., Perenos ionov v membranakh [Ion Transport in Membranes], Nauka, Moscow, 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., 53:11 (2017), 1254–1265 | DOI | DOI