Kovalenko Anna, Chubyr Natalia, Uzdenova Aminat, Urtenov Makhamet
Department of Data Analysis and Artificial Intelligence, Kuban State University, Krasnodar 350040, Russia.
Department of Applied Mathematics, Kuban State University, Krasnodar 350040, Russia.
Membranes (Basel). 2022 Oct 26;12(11):1047. doi: 10.3390/membranes12111047.
At present, it is customary to consider the overlimit operating modes of electromembrane systems to be effective, and electroconvection as the main mechanism of overlimiting transfer. The breakdown of the space charge is a negative, "destructive" phenomenon, since after the breakdown the size and number of electroconvective vortices are significantly reduced, which leads to a decrease in mass transfer. Therefore, electromembrane desalination processes must be carried out before space charge breakdown occurs. Thus, the actual problem arises of determining at which potential jumps a breakdown of the space charge occurs at a given concentration of the solution. Electromembrane systems are used for desalination at electrolyte solution concentrations ranging from 1 to 100 mol/m3. In a theoretical study of increasing the efficiency of the desalination process, mathematical modeling is used in the form of a boundary value problem for the system of Nernst-Planck and Poisson (NPP) equations, which refers to "hard" problems that are difficult to solve numerically. This is caused by the appearance of a small parameter at the derivative in the Poisson equation in a dimensionless form, and, correspondingly, a boundary layer at ion-exchange membranes, where concentrations and other characteristics of the desalination process change exponentially. It is for this reason that the numerical study of the boundary value problem is currently obtained for initial concentrations of the order of 0.01 mol/m3. The paper proposes a new numerical-analytical method for solving boundary value problems for the system of Nernst-Planck and Poisson equations for real initial concentrations, using which the phenomenon of space charge breakdown (SCB) in the cross section of the desalination channel in potentiostatic and potentiodynamic modes is studied. The main regularities of the appearance and interaction of charge waves, up to their destruction (breakdown), are established. A simple formula is proposed for engineering calculations of the potential jump depending on the concentration of the solution, at which the breakdown of the space charge begins.
目前,人们通常认为电膜系统的超限运行模式是有效的,并将电对流视为超限传递的主要机制。空间电荷的击穿是一种负面的“破坏性”现象,因为击穿后电对流涡旋的大小和数量会显著减少,从而导致传质减少。因此,电膜脱盐过程必须在空间电荷击穿发生之前进行。这样一来,就出现了一个实际问题,即确定在给定溶液浓度下,空间电荷在哪个电位跃变时会发生击穿。电膜系统用于电解质溶液浓度范围为1至100 mol/m³的脱盐。在提高脱盐过程效率的理论研究中,采用了以能斯特 - 普朗克方程和泊松方程(NPP)系统的边值问题形式进行的数学建模,这属于难以用数值方法求解的“硬”问题。这是由于泊松方程以无量纲形式出现时导数处出现小参数,相应地在离子交换膜处出现边界层,脱盐过程的浓度和其他特性在该边界层呈指数变化。正是由于这个原因,目前边值问题的数值研究是针对0.01 mol/m³量级的初始浓度进行的。本文提出了一种新的数值分析方法来求解能斯特 - 普朗克方程和泊松方程系统对于实际初始浓度的边值问题,利用该方法研究了恒电位和动电位模式下脱盐通道横截面中的空间电荷击穿(SCB)现象。确定了电荷波出现和相互作用直至其破坏(击穿)的主要规律。提出了一个简单公式,用于根据溶液浓度对电位跃变进行工程计算,在该电位跃变时空间电荷开始击穿。