**Reasons for the Formation and Properties of Soliton-Like Charge Waves in Membrane Systems When Using Overlimiting Current Modes**

#### **Makhamet Urtenov 1,\* , Natalia Chubyr <sup>2</sup> and Vitaly Gudza <sup>1</sup>**


Received: 22 July 2020; Accepted: 13 August 2020; Published: 16 August 2020

**Abstract:** The study of ion transport in membrane systems in overlimiting current modes is an important problem of physical chemistry and has an important application value. The influence of the space charge on the transport of salt ions under overlimiting current modes was first studied in the work of Rubinstein and Shtilman and later in the works of many authors. The purpose of this research is to study, using the method of mathematical modeling, the reasons of formation and properties of the local maximum (minimum) space charge in membrane systems under overlimiting current conditions. It is shown that, in the diffusion layer of the cation-exchange membrane (CEM), the local maximum of the space charge appears due to the limited capacity (exchange capacity) of the membrane at a given potential jump, i.e., the local maximum of space charge appears due to the presence of a local minimum of space charge at the surface of the CEM. The local maximum of the space charge moves as a single soliton-like wave into the depth of the solution. Unlike real solitons, this charged wave changes its size and shape, albeit quite slowly. In the section of the desalination channel, the situation is completely different. First, the space charge of the anion-exchange membrane (AEM) has a negative value, so we should be talking about the local minimum (or the maximum of the absolute value of the charge). However, this is an insignificant clarification. Secondly, the space charge waves of different signs begin to interact, which leads to a new effect, namely the effect of the breakdown of the space charge. The dependence of the local maximum on the input parameters—the cation diffusion coefficient, the growth rate of the potential jump, and the initial and boundary concentrations—is studied.

**Keywords:** ion-exchange membrane; mathematical modelling; using overlimiting current modes; membrane systems; cation-exchange membrane; effect of the breakdown of the space charge

#### **1. Introduction**

A thorough development of membrane devices has been done over the decades [1]. Of course, progress in separation through phase boundaries/interfaces has been accompanied by issues. The problem of ion transport across phase boundaries is one of the fundamental problems of physical chemistry and electrochemistry and is also important for membrane technologies. Studies [2–7] have shown the prospects of using intensive current modes.

The influence of space charge on the structure of the diffusion layer was first studied by Rubinstein and Shtilman [8]. Instead of the traditional equations of electroneutrality to the system of equations of the Nernst–Planck diffusion layer, the authors introduced the Poisson equation, and the ion-exchange membrane was taken as selective (the effective transport number of counterions was taken independent of current density). The problem was solved numerically. In the future, various mathematical methods and approaches to solving the Nernst–Planck–Poisson equations for electromembrane systems under extreme current conditions were developed [9–14].

The study of the non-stationary transfer of binary electrolyte in the diffusion layer is interesting because it allows us to determine the structure of the diffusion layer and its change over time, which is necessary, for example, for the asymptotic analysis of problems and the establishment of simple engineering calculation formulas for the dependence of the concentration distribution and electric field strength on the parameters of the problem. Research on non-stationary problems is limited to studies [15–18]. In these articles, the main attention is paid to the overlimiting potential dynamic mode and the analysis of the setting time depending on the parameters of the problem.

Studies [19,20] are devoted to the study of «shock electrodialysis», in which a deionization wave propagates through a microchannel or a porous medium with a sharp boundary between the concentrated and depleted zones. The deionization waves can be compared with the charge wave, since the deionization region actually coincides with the space charge region and a «powerful gradient» occurs in the desalination channel. Shock electrodialysis is a newly developed method for the desalination of water and deionization in micro-scale pores near an ion-selective element. In contrast to «shock electrodialysis», we study electro-membrane systems (the diffusion layer and the section of the desalination channel) with macroscale dimensions of the order of millimeters. In addition, it is shown that, in these systems, the interaction of charge waves is possible, up to their destruction (breakdown).

#### **2. Materials and Methods**

#### *Mathematical Model of One-Dimensional Non-Stationary Ion Transport in Membrane Systems*

Two electromembrane systems will be considered below: The depleted diffusion layer at the cation-exchange membrane (CEM), then *x* = 0, the beginning of the diffusion layer, and *x* = *H* is the conditional boundary of the solution/CEM. Additionally, the section of the desalination channel, in this case *x* = 0, is the conditional boundary of the anion-exchange membrane (AEM)/solution, and *x* = *H* has the same meaning. The transfer of 1:1 salt ions in both cases is described by the same equations—the difference in the boundary conditions at *x* = 0. System of equations.

The non-stationary transport of salt ions for a 1:1 electrolyte is described by the following system of equations:

$$\frac{\partial \mathbf{C}\_i}{\partial t} = -\frac{\partial j\_i}{\partial \mathbf{x}} \qquad i = 1,2\tag{1}$$

$$j\_i = -\frac{F}{RT} z\_i D\_i \mathbf{C}\_i \frac{\partial \phi}{\partial \mathbf{x}} - D\_i \frac{\partial \mathbf{C}\_i}{\partial \mathbf{x}} \qquad \qquad i = 1, 2 \tag{2}$$

$$\frac{\partial^2 \phi}{\partial \mathbf{x}^2} = -\frac{F}{\varepsilon\_d} (\mathbb{C}\_1 - \mathbb{C}\_2) \tag{3}$$

$$I\_{\mathbb{C}} = F(j\_1 - j\_2) \tag{4}$$

Here, (1) is the equation of material balance, (2) is the equation the Nernst–Planck for fluxes of sodium *<sup>i</sup>* = <sup>1</sup> <sup>↔</sup> *Na*<sup>+</sup> and chloride *<sup>i</sup>* = <sup>2</sup> <sup>↔</sup> *Cl*<sup>−</sup> ions, the charge number of cations *<sup>z</sup>*<sup>1</sup> = <sup>1</sup> and anions *z*<sup>2</sup> = −1, (3) is a Poisson equation for the electric field potential, (4) is the equation of the current flow, which means that the current flowing through the diffusion layer is determined by the flow of ions, ε*<sup>a</sup>* is the dielectric permeability of the solution, *F* is the Faraday number, *R* is the universal gas constant, φ is the potential, *E* = − ∂φ ∂*x* is the electric field intensity, and *C<sup>i</sup>* , *ji* , *D<sup>i</sup> I<sup>c</sup>* are concentration, flux, the diffusion coefficient of the i-th ion, and current density determined by the flux of ions, respectively.

Boundary conditions for the diffusion layer.

The boundary conditions for a diffusion layer consist of the following boundary and initial conditions:

$$\mathbb{C}\_{1}(t,0) = \mathbb{C}\_{0}, \mathbb{C}\_{2}(t,0) = \mathbb{C}\_{0}, \phi(t,0) = 0 \tag{5}$$

$$\mathbb{C}\_{1}(t,H) = \mathbb{C}\_{1m}(t), \left(\frac{\partial \mathbb{C}\_{2}}{\partial \mathbf{x}} + \frac{F}{RT}\mathbb{C}\_{2}\mathbb{E}\right)(t,H) = 0, \phi(t,H) = \Delta\_{\mathbf{r}}(t),$$

$$\mathbb{C}\_{10}(\mathbf{x}) = \mathbb{C}\_{0}, \mathbb{C}\_{20}(\mathbf{x}) = \mathbb{C}\_{0}, \phi\_{0}(\mathbf{x}) = \mathbf{0},$$

where ∆*r*(*t*) is the potential jump. As a rule, ∆*r*(*t*) is either constant, or ∆*r*(*t*) = −*d*·*t*, where *d* is the growth rate of the potential jump and has the dimension V/s.

Boundary conditions for the desalination channel.

To model the ion transport, replace the boundary condition (5) with the following:

$$\mathbf{C}\_{2}(t,0) = \mathbf{C}\_{2m} \left( \frac{\partial \mathbf{C}\_{1}}{\partial \mathbf{x}} - \frac{F}{RT} \mathbf{C}\_{1} E \right)(t, H) \\ = \mathbf{0}, \\ \phi(t,0) = \mathbf{0} \tag{6}$$

#### **3. Results**

*3.1. Reasons for the Formation of a Local Maximum (Minimum) Space Charge in the Extended Space Charge Region (Extended SCR)*

When the current passes through the membrane, the concentration of cations decreases and reaches its minimum at the left border of the border layer. When using overlimiting current densities, this minimum is preserved, but the minimum value is greatly reduced, and there are practically no anions in the area of the minimum point. The size of the space charge ρ = *F*(*C*<sup>1</sup> − *C*2) is almost completely determined by the concentration of cations ρ ≈ *FC*1(Figure 1b). Therefore, the space charge also has a local minimum other than zero at this point. The minimum value depends on the applied potential jump (or on the current that is passed through the system). On the other hand, in the depth of the solution, where the condition of electroneutrality is met, the value of the space charge is almost zero. Therefore, the value of the space charge must have a local maximum between them (Figure 1a). In [9,12], it is shown that *x* ∗ = *Inp <sup>I</sup> <sup>H</sup>*; therefore, [0, *<sup>I</sup>np <sup>I</sup> H*) is the region of electroneutrality, and ( *Inp <sup>I</sup> H*, *H*) is the extended space charge region (SCR). In addition, in the field of electroneutrality, a balance is observed between the processes of diffusion and electromigration, and the currents of diffusion and electromigration (ohmic) are equal. At the same time, in the SCR, the electromigration is an order of magnitude greater than the diffusion one. Thus, the current pattern changes at the local maximum point. *Membranes* **2020**, *10*, x FOR PEER REVIEW 4 of 11

**Figure 1.** (**а**) Graph of the space charge in the diffusion layer; (**b**) Comparison with the graph of the concentration of *С*<sup>1</sup> cations near the cation-exchange membrane (CEM) when using overlimiting **Figure 1.** (**a**) Graph of the space charge in the diffusion layer; (**b**) Comparison with the graph of the concentration of *C*<sup>1</sup> cations near the cation-exchange membrane (CEM) when using overlimiting current modes in the dimensional form.

<sup>0</sup> <sup>1</sup> <sup>=</sup> <sup>∂</sup> ∂ *x*

*RT <sup>F</sup> <sup>j</sup>* <sup>∂</sup>

> *x <sup>a</sup>* ε −= <sup>∂</sup>

> > *Fj*<sup>1</sup> *Ic* =

2 2

<sup>ε</sup> <sup>=</sup> ∂

2

*x <sup>a</sup>*

ϕ∂ *xj <sup>D</sup> RT*

∂ ϕ∂

2 1 *<sup>x</sup> <sup>x</sup> <sup>a</sup>* <sup>ε</sup> <sup>=</sup> <sup>∂</sup> ϕ∂

2

*CD*

*<sup>C</sup> <sup>F</sup>*

2 1

( ) *<sup>F</sup> )t(I tjj <sup>c</sup>* <sup>11</sup> ==

ϕ∂ , then, taking into account equation (8), we have:

*j D RT*

> 1 1

1 1

β+

*x*

ϕ∂

*<sup>j</sup>* (7)

<sup>1</sup> −= <sup>11</sup> (8)

ϕ∂ (9)

Using the above assumptions, you can analytically determine the solution in the extended SCR. Indeed, we put in the system of equations <sup>2</sup> *)x,t(C* = 0 , <sup>2</sup> *)x,t(j* = 0 ,

current modes in the dimensional form.

∂ >>

*x C*

, then we get:

∂*x*

*x*

From Equation (7), we get:

Multiply Equation (9) by

∂ ϕ∂

*F <sup>i</sup> ii* ∂

*Cz*

*RT*

0

Using the above assumptions, you can analytically determine the solution in the extended SCR. Indeed, we put in the system of equations *C*2(*t*, *x*) = 0, *j*2(*t*, *x*) = 0, *F RT*0 *ziC<sup>i</sup>* ∂φ ∂*x* >> ∂*C<sup>i</sup>* ∂*x* , then we get:

$$\frac{\partial j\_1}{\partial \mathbf{x}} = \mathbf{0} \tag{7}$$

$$j\_1 = -\frac{F}{RT}D\_1\mathcal{C}\_1\frac{\partial\phi}{\partial\mathbf{x}}\tag{8}$$

$$\frac{\partial^2 \phi}{\partial \mathbf{x}^2} = -\frac{F}{\varepsilon\_d} \mathbf{C}\_1 \tag{9}$$

$$I\_c = Fj\_1$$

From Equation (7), we get:

$$j\_1 = \dot{j}\_1(t) = \frac{I\_c(t)}{F}$$

Multiply Equation (9) by ∂φ ∂*x* , then, taking into account equation (8), we have:

$$\frac{\partial \phi}{\partial \mathbf{x}} \frac{\partial^2 \phi}{\partial \mathbf{x}^2} = \frac{RT}{\varepsilon\_d D\_1} j\_1$$

$$\frac{1}{2} \left[ \frac{\partial \phi}{\partial \mathbf{x}} \right]^2 = \frac{RT}{\varepsilon\_d D\_1} j\_1 \mathbf{x} + \beta$$

$$\frac{\partial \phi}{\partial \mathbf{x}} = -\sqrt{2 \frac{RT}{\varepsilon\_d D\_1} j\_1 \mathbf{x} + \beta} \tag{10}$$

or

where β > 0 is the integration constant. From (8) and (10), we have:

$$\mathcal{C}\_1(t, \mathbf{x}) = \frac{j\_1 R T}{F D\_1 \sqrt{2 \frac{R T}{\varepsilon\_d D\_1} j\_1 \mathbf{x} + \boldsymbol{\beta}}}$$

or

$$\mathcal{C}\_1(t, \mathbf{x}) = \frac{RT}{F^2 D\_1 \sqrt{\frac{2RT}{\varepsilon\_d F D\_1} I\_c(t) \mathbf{x} + \beta}} I\_c(t)$$

From this formula, it can be seen that *C*1(*t*, *x*) monotonously decreases in *x* in the extended SCR and reaches a minimum on its right border. In addition, as the current density increases, the minimum value decreases, which coincides with the numerical solution.

In order to understand the role and significance of the occurrence and growth of a local maximum, it is necessary to study its dependence on the input parameters of the problem, such as the initial and boundary concentrations, the rate of growth of the potential jump, the composition of the electrolyte solution, etc.

#### *3.2. Dependence of the Local Maximum in the Di*ff*usion Layer on the Boundary Concentration of C*1*<sup>m</sup>*

Consider the dependence of the local maximum on the parameter *C*1*m*, which characterizes the exchange capacity of the membrane.

The value in Figure 2 increases by the formula *C*1*<sup>m</sup>* = 0.5·*t*, with the local maximum appearing and disappearing. The point of the local maximum is shifted to the right, but the value does not change much (Figure 2a,b). At the same time, the local minimum is gradually filled with increasing *C*1*<sup>m</sup>* (the value of the local minimum increases) and it is first compared and then surpassed (Figure 2c,d).

3.4.1. Diffusion Layer

*Membranes* **2020**, *10*, x FOR PEER REVIEW 6 of 11

**Figure 2.** Dependence of the local maximum on the exchange capacity of the CEM (NaCl, 0.5 ∆ = *<sup>r</sup>* ϕ ): (**a**) 14 *C*1*<sup>m</sup>* = , (**b**) 350 *C*1*<sup>m</sup>* = , (**c**) 920 *C*1*<sup>m</sup>* = , (**d**) 1335 *C*1*<sup>m</sup>* = **Figure 2.** Dependence of the local maximum on the exchange capacity of the CEM (NaCl, ∆*r*φ = 0.5): (**a**) *C*1*<sup>m</sup>* = 14, (**b**) *C*1*<sup>m</sup>* = 350, (**c**) *C*1*<sup>m</sup>* = 920, (**d**) *C*1*<sup>m</sup>* = 1335.

*3.3. Soliton-like Charge Wave. Dependence of the Local Maximum (Minimum) Charge on the Growth Rate of the Potential Jump d* Thus, it can be concluded that the local maximum of the space charge appears due to the limited exchange capacity of the membrane at a given potential jump, i.e., the local maximum of the space charge appears due to the presence of a local minimum of the space charge at the surface of the CEM.

#### Consider a non-stationary model of transport in the diffusion layer of the CEM (1–5) with a linear growth of the potential jump *(t,H ) d t <sup>r</sup>* ∆ ϕ = − ⋅ . *3.3. Soliton-like Charge Wave. Dependence of the Local Maximum (Minimum) Charge on the Growth Rate of the Potential Jump d*

The calculations are at times greater than a certain critical *kt* , at which point the potential jump *<sup>r</sup> <sup>k</sup> <sup>k</sup>* ∆ ϕ*(t ,H )* = −*d* ⋅*t* corresponds to the limiting current density and a local maximum of the Consider a non-stationary model of transport in the diffusion layer of the CEM (1–5) with a linear growth of the potential jump ∆*r*φ(*t*, *H*) = −*d*·*t*.

space charge appears as a single soliton-like wave begins to move into the depth of the solution. Unlike true solitons, this wave has a charge. In addition, it changes its size and shape. With a decrease in the growth rate of the potential jump *d* , the space charge wave becomes The calculations are at times greater than a certain critical *t<sup>k</sup>* , at which point the potential jump ∆*r*φ(*t<sup>k</sup>* , *H*) = −*d*·*t<sup>k</sup>* corresponds to the limiting current density and a local maximum of the space charge appears as a single soliton-like wave begins to move into the depth of the solution. Unlike true solitons, this wave has a charge. In addition, it changes its size and shape.

less pronounced, but remains. *3.4. Dependence of Charge Waves on the Cation Diffusion Coefficient* With a decrease in the growth rate of the potential jump *d*, the space charge wave becomes less pronounced, but remains.

and maximum and their value in the diffusion layer of CEM differ slightly (Figure 3).

In order to study the dependence of the local maximum on the diffusion coefficient of the cation, solutions of NaCl and KCl are considered. Calculations show that, despite the difference in the diffusion coefficients of sodium and potassium by 1.5 times, the positions of the local minimum

#### *3.4. Dependence of Charge Waves on the Cation Di*ff*usion Coe*ffi*cient*

#### 3.4.1. Diffusion Layer

In order to study the dependence of the local maximum on the diffusion coefficient of the cation, solutions of NaCl and KCl are considered. Calculations show that, despite the difference in the diffusion coefficients of sodium and potassium by 1.5 times, the positions of the local minimum and maximum *Membranes*  and their value in the diffusion layer of CEM differ slightly (Figure 3). **2020**, *10*, x FOR PEER REVIEW 7 of 11

**Figure 3.** Wave charge when *d* = 0*.*1 when *t* ≈ 9*с* (**a**) and *t* ≈100*с* (**b**). The blue and green lines (they merge with great accuracy) correspond to solutions of NaCl and KCl. **Figure 3.** Wave charge when *d* = 0.1 when *t* ≈ 9*c* (**a**) and *t* ≈ 100*c* (**b**). The blue and green lines (they merge with great accuracy) correspond to solutions of NaCl and KCl.

#### 3.4.2. Section of the Desalination Channel

greater than in the second.

3.4.2. Section of the Desalination Channel Let us analyze the issue (1–4), (6). In the section of the desalination channel, the situation is completely different. First, the space charge of AEM has a negative value, so we must talk about the local minimum (or the local maximum of the absolute value of the charge). However, this is an insignificant clarification. Secondly, the space charge waves of different signs begin to interact, Let us analyze the issue (1–4), (6). In the section of the desalination channel, the situation is completely different. First, the space charge of AEM has a negative value, so we must talk about the local minimum (or the local maximum of the absolute value of the charge). However, this is an insignificant clarification. Secondly, the space charge waves of different signs begin to interact, which leads to a new effect, namely the effect of the breakdown of the space charge.

which leads to a new effect, namely the effect of the breakdown of the space charge. Let us first consider the behavior of the local maximum in the KCl solution at d = 0.1. As can be seen from Figure 4a–f, two symmetrical soliton-like lone waves are formed, which move towards each other. Unlike real solitons, these waves have, as noted above, charges, the left wave negative charges and the right positive charges. At first, they practically do not interact, but as they approach each other they begin to attract and their speed of convergence increases and at the moment before contact there is a practical instantaneous breakdown, and they are discharged. A further increase in the potential jump does not lead to the formation of a new wave of local maximum, since the concentration of the solution practically becomes zero (Figure 4f–h), with the exception of narrow border layers in AEM and CEM, due to the fact that the concentration of anions and cations is maintained constant at the boundaries of AEM/solution and solution/CEM. As the potential jump increases, the width of the border layers decreases very slowly. A similar scenario is realized for the NaCl solution, except that the local minimum of negative space charge is generated much later Let us first consider the behavior of the local maximum in the KCl solution at d = 0.1. As can be seen from Figure 4a–f, two symmetrical soliton-like lone waves are formed, which move towards each other. Unlike real solitons, these waves have, as noted above, charges, the left wave negative charges and the right positive charges. At first, they practically do not interact, but as they approach each other they begin to attract and their speed of convergence increases and at the moment before contact there is a practical instantaneous breakdown, and they are discharged. A further increase in the potential jump does not lead to the formation of a new wave of local maximum, since the concentration of the solution practically becomes zero (Figure 4f–h), with the exception of narrow border layers in AEM and CEM, due to the fact that the concentration of anions and cations is maintained constant at the boundaries of AEM/solution and solution/CEM. As the potential jump increases, the width of the border layers decreases very slowly. A similar scenario is realized for the NaCl solution, except that the local minimum of negative space charge is generated much later (Figure 4f) than for the KCl solution. This is due to the fact that the current density in the first case is greater than in the second.

(Figure 4f) than for the KCl solution. This is due to the fact that the current density in the first case is

(**a**) (**b**)

3.4.2. Section of the Desalination Channel

(**a**)(**b**)

**Figure 3.** Wave charge when *d* = 0*.*1 when *t* ≈ 9*с* (**a**) and *t* ≈100*с* (**b**). The blue and

Let us analyze the issue (1–4), (6). In the section of the desalination channel, the situation is completely different. First, the space charge of AEM has a negative value, so we must talk about the local minimum (or the local maximum of the absolute value of the charge). However, this is an insignificant clarification. Secondly, the space charge waves of different signs begin to interact,

Let us first consider the behavior of the local maximum in the KCl solution at d = 0.1. As can be seen from Figure 4a–f, two symmetrical soliton-like lone waves are formed, which move towards each other. Unlike real solitons, these waves have, as noted above, charges, the left wave negative charges and the right positive charges. At first, they practically do not interact, but as they approach each other they begin to attract and their speed of convergence increases and at the moment before contact there is a practical instantaneous breakdown, and they are discharged. A further increase in the potential jump does not lead to the formation of a new wave of local maximum, since the concentration of the solution practically becomes zero (Figure 4f–h), with the exception of narrow border layers in AEM and CEM, due to the fact that the concentration of anions and cations is maintained constant at the boundaries of AEM/solution and solution/CEM. As the potential jump increases, the width of the border layers decreases very slowly. A similar scenario is realized for the

green lines (they merge with great accuracy) correspond to solutions of NaCl and KCl.

which leads to a new effect, namely the effect of the breakdown of the space charge.

**Figure 4.** *Cont*.

(**g**) (**h**)

**Figure 4.** Graphs of the space charge at *d* = 0*.*1 *,* <sup>3</sup> *C*<sup>0</sup> = 0*.*1 *mol / m* . *<sup>t</sup>* <sup>≈</sup> <sup>29</sup>*<sup>с</sup>* (**a**); *t* ≈ 68*с* (**b**); *t* ≈ 82*с* (**c**); *t* ≈ 86*с* (**d**); *t* ≈ 97*с* (**e**); *t* ≈108*с* (**f**); *t* ≈112*с* (**g**);

At *d* = 0*.*005 in NaCl, a single wave occurs in the CEM, which reaches the region of negative space charge in the AEM, where there is no local minimum and maximum, and is discharged at 436 *s*. In KCl, two almost symmetrical waves are formed, one in the CEM and the

*t* ≈116*с* (**h**).

(**e**) (**f**)

(**c**) (**d**)

**Figure 4.** Graphs of the space charge at *d* = 0*.*1 *,* <sup>3</sup> *C*<sup>0</sup> = 0*.*1 *mol / m* . *<sup>t</sup>* <sup>≈</sup> <sup>29</sup>*<sup>с</sup>* (**a**); *t* ≈ 68*с* (**b**); *t* ≈ 82*с* (**c**); *t* ≈ 86*с* (**d**); *t* ≈ 97*с* (**e**); *t* ≈108*с* (**f**); *t* ≈112*с* (**g**); **Figure 4.** Graphs of the space charge at *d* = 0.1, *C*<sup>0</sup> = 0.1 *mol*/*m*<sup>3</sup> . *t* ≈ 29*c* (**a**); *t* ≈ 68*c* (**b**); *t* ≈ 82*c* (**c**); *t* ≈ 86*c* (**d**); *t* ≈ 97*c* (**e**); *t* ≈ 108*c* (**f**); *t* ≈ 112*c* (**g**); *t* ≈ 116*c* (**h**). *Membranes* **2020**, *10*, x FOR PEER REVIEW 9 of 11

*t* ≈116*с* (**h**). At *d* = 0*.*005 in NaCl, a single wave occurs in the CEM, which reaches the region of negative space charge in the AEM, where there is no local minimum and maximum, and is discharged at 436 *s*. In KCl, two almost symmetrical waves are formed, one in the CEM and the At *d* = 0.005 in NaCl, a single wave occurs in the CEM, which reaches the region of negative space charge in the AEM, where there is no local minimum and maximum, and is discharged at 436 *s*. In KCl, two almost symmetrical waves are formed, one in the CEM and the other in the AEM, which meet almost in the middle and discharge at 282 *c*. Thus, the behavior of charge waves at *d* = 0.005 in NaCl and KCl solutions is completely different. other in the AEM, which meet almost in the middle and discharge at 282 *c* . Thus, the behavior of charge waves at *d* = 0*.*005 in NaCl and KCl solutions is completely different. *3.5. Dependence of Charge Waves on the Initial Concentration of C*<sup>0</sup>

#### *3.5. Dependence of Charge Waves on the Initial Concentration of C*<sup>0</sup> The comparison of space charge graphs for <sup>3</sup> *C*<sup>0</sup> = 0*.*01 *mol / m* and for

The comparison of space charge graphs for *C*<sup>0</sup> = 0.01 mol/m<sup>3</sup> and for *C*<sup>0</sup> = 0.1 mol/m<sup>3</sup> (see Figures 4a and 5a) at *t* = 29 *s* shows that the less concentrated solution is desalinated faster and the charge waves are already very different, and, at *t* = 54.3 s (Figure 5b), the charge waves for KCl began to discharge. <sup>3</sup> *C*<sup>0</sup> = 0*.*1 *mol / m* (see Figures 4a and 5a) at *t* = 29 *s* shows that the less concentrated solution is desalinated faster and the charge waves are already very different, and, at *t* = 54*.*3 *s* (Figure 5b), the charge waves for KCl began to discharge.

**Figure 5.** Space charge graph at <sup>3</sup> *C*<sup>0</sup> = 0*.*01 *mol / m* . *<sup>t</sup>* <sup>≈</sup> <sup>29</sup>*<sup>с</sup>* (**a**); *<sup>t</sup>* <sup>≈</sup> <sup>54</sup>*<sup>с</sup>* (**b**) **Figure 5.** Space charge graph at *C*<sup>0</sup> = 0.01 mol/m<sup>3</sup> . *t* ≈ 29*c* (**a**); *t* ≈ 54*c* (**b**).

#### **4. Conclusions**

this phenomenon are studied.

V.G..; Visualization, M.U., N.C., V.G.

desalination in intensive current regimes.

**4. Conclusions** In this paper, the reasons for the formation and properties of the local maximum (minimum) space charge in membrane systems when using overlimiting current conditions are investigated. In this paper, the reasons for the formation and properties of the local maximum (minimum) space charge in membrane systems when using overlimiting current conditions are investigated. The depleted diffusion layer at the cation-exchange membrane and the section of the desalination

jump. It is shown that the local maximum of the space charge in the diffusion layer moves as a single soliton-like wave into the depth of the solution, slowly changing its size and shape. In the section of the desalination channel, the space charge waves of different signs begin to interact, which leads to a new effect, namely the effect of discharge (breakdown) of the space charge. The fundamental laws of

**Author Contributions:** Conceptualization, M.U., N.C., V.G.; Investigation, M.U., N.C.; Formal Analysis, M.U.,

**Funding:** This research was funded by Russian Foundation for Basic Research project number 20-58-12018 NNIO\_a: The influence of electroconvection, water dissociation, and geometry of spacers on electrodialysis

**Acknowledgments:** The study was carried out with the financial support of the Russian Foundation for Basic **Research:** grant number 20-58-12018 NNIO\_a: The influence of electroconvection, water dissociation, and

geometry of spacers on electrodialysis desalination in intensive current regimes.

**Conflicts of Interest:** The authors declare no conflict of interest.

channel are considered as membrane systems. It is shown that the local maximum and minimum of space charge appear due to the limited capacity of ion-exchange membranes at a given potential jump. It is shown that the local maximum of the space charge in the diffusion layer moves as a single soliton-like wave into the depth of the solution, slowly changing its size and shape. In the section of the desalination channel, the space charge waves of different signs begin to interact, which leads to a new effect, namely the effect of discharge (breakdown) of the space charge. The fundamental laws of this phenomenon are studied.

**Author Contributions:** Conceptualization, M.U., N.C., V.G.; Investigation, M.U., N.C.; Formal Analysis, M.U., V.G.; Visualization, M.U., N.C., V.G. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Russian Foundation for Basic Research project number 20-58-12018 NNIO\_a: The influence of electroconvection, water dissociation, and geometry of spacers on electrodialysis desalination in intensive current regimes.

**Acknowledgments:** The study was carried out with the financial support of the Russian Foundation for Basic Research grant number 20-58-12018 NNIO\_a: The influence of electroconvection, water dissociation, and geometry of spacers on electrodialysis desalination in intensive current regimes.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

*Article*
