Next Article in Journal
Peri-Implant Mucosa Augmentation with an Acellular Collagen Matrix
Next Article in Special Issue
Fouling Mitigation via Chaotic Advection in a Flat Membrane Module with a Patterned Surface
Previous Article in Journal
Model-Based Analysis of Low Stoichiometry Operation in Proton Exchange Membrane Water Electrolysis
Previous Article in Special Issue
Effect of Modifying the Membrane Surface with Microcapsules on the Flow Field for a Cross-Flow Membrane Setup: A CFD Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Continuous Extraction of Li+ from High Mg2+/Li+ Ratio Brines Based on Free Flow Ion Concentration Polarization Microfluidic System

1
College of Mechanical and Electrical Engineering, Wenzhou University, Wenzhou 325000, China
2
School of Mechanical Engineering, Hebei University of Technology, Tianjin 300401, China
3
National Engineering Research Center for Technological Innovation Method and Tool, Tianjin 300401, China
*
Authors to whom correspondence should be addressed.
Membranes 2021, 11(9), 697; https://doi.org/10.3390/membranes11090697
Submission received: 20 July 2021 / Revised: 8 September 2021 / Accepted: 8 September 2021 / Published: 10 September 2021
(This article belongs to the Special Issue Numerical Modelling in Membrane Processes)

Abstract

:
Ion concentration polarization (ICP) is a promising mechanism for concentrating and/or separating charged molecules. This work simulates the extraction of Li+ ions in a diluted high Mg2+/Li+ ratio salt lake brines based on free flow ICP focusing (FF-ICPF). The model solution of diluted brine continuously flows through the system with Li+ slightly concentrated and Mg2+ significantly removed by ICP driven by external pressure and perpendicular electric field. In a typical case, our results showed that this system could focus Li+ concentration by ~1.28 times while decreasing the Mg2+/Li+ ratio by about 85% (from 40 to 5.85). Although Li+ and Mg2+ ions are not separated as an end product, which is preferably required by the lithium industry, this method is capable of decreasing the Mg2+/Li+ ratio significantly and has great potential as a preprocessing technology for lithium extraction from salt lake brines.

1. Introduction

Global demand for lithium has increased significantly over the past decades, driven by the expanding requirement of rechargeable lithium batteries for portable electronic devices, electric vehicles, and grid storage applications. Generally, lithium has been obtained from two sources—hard rock ores and continental brines [1,2]. Extraction of lithium from ores/minerals is mainly through calcination. In this method, lithium ore is calcined or roasted and then leached to dissolve lithium into the liquid phase. In this process, a large amount of hydrogen chloride gas is produced, which has serious pollution consequences [3]. As a result, lithium production from alternative resources, i.e., salt lake brines, is increasing steadily over the last decades [4].
Essentially, the key task of lithium extraction from salt lake brines is the separation of Li+ and Mg2+ ions. This is because magnesium and lithium elements possess similar chemical properties, and magnesium must be removed in most lithium applications [5,6]. So far, several methods have been developed for extracting lithium from salt lake brines, such as precipitation [7,8], adsorption [9], solvent extraction [10,11], calcination [12], and membrane approaches [13,14], etc.
Generally, the precipitation method uses chemical reactions to convert Li+ into insoluble substances and then separate them from the solution. This method is commonly used in low Mg2+/Li+ ratio brines of salt lakes in South America, yielding the largest portion of the world’s lithium productions today [15]. However, this procedure needs a large area and a long processing time [16], and it does not work for high Mg2+/Li+ ratio salt lake brines. The adsorption method uses ion sieve adsorbents to absorb Li+ selectively from the brine and extract lithium in subsequent stages. It works well for high Mg2+/Li+ ratio brines. Still, its processes are complicated, time-consuming, and suffering from the cost and efficiency problems of the absorbents. Solvent extraction uses an organic extractant to extract Li+ from the brine. However, the organic extractant may solve into the brine in this process, causing environmental pollution [17]. Calcination technology uses dried salts obtained from the brines through heating and dissolution. The production procedures and the suffering problems are the same as those in the calcination of mineral ores [18]. Membrane technologies use membranes to separate Li+ and Mg2+, which could be driven by mechanical pressure (nanofiltration) [13,14] or by the electric field (electrodialysis) [19]. In a nanofiltration system, high pressure is used to push Li+ through the nanoporous membranes while Mg2+ cannot pass through. The advantages of nanofiltration are low energy consumption, environmental friendliness, and easy operation. However, the total salinity of the feed solution must be low, and the cost of an ion-selective nanofiltration membrane is high [20]. In an electrodialysis system, a monovalent ion exchange membrane is used to transport monovalent Li+ through the membrane while blocking the permeation of multivalent Mg2+ [21,22]. This method requires high-quality membranes, the cost of which is still hindering its industrial applications. Further research is required to develop new technologies based on completely new fundamental principles. In this aspect, the ion concentration polarization (ICP) based method may play a promising role.
A preferable method for lithium extraction must be able to concentrate Li+ and remove Mg2+ simultaneously. For the first task, i.e., the concentration of a specific kind of ions, or generally a charged species, Wang et al. developed a micro-nanochannel ICP system to preconcentrate proteins and peptides [23]. In their pioneer system, two microchannels (the main channel and a buffer one) are connected by an array of nanochannels. An electric field is applied through nanochannels from the main channel to the buffer channel. Because the walls of nanochannels are negatively charged, only cations and positively charged species can pass through nanochannels. As a result, the concentration of ions in the main channel near micro-nanochannel interfaces is lower than that in the buffer channel, i.e., ICP is induced across two ends of nanochannels [24,25,26,27]. If the electric field across two microchannels is so high that the concentration of the ions is close to zero near the openings of nanochannels in the main channel, an ion depletion zone (IDZ) is formed there. The electric field in IDZ is significantly higher than those at other locations [28]. This locally amplified electric field in IDZ can hinder the fluid-drag motion of the co-charged species (with respect to the charges on the nanochannel walls). In Wang’s work [23], when the negatively charged proteins or peptides migrate from the inlet to the outlet (carried by fluid flow), they are subjected to a sharply increasing resistive electric force at the boundary of the IDZ. For a specific molecule, if the electric force (FE) applied to it becomes equal to the fluid drag force (FD) at a specific location, it will be trapped there. When a large number of molecules arrive at this location, they gather and get concentrated there [29]. Using this method, Wang et al. achieved a million-fold concentration for protein and peptide molecules [23].
If there is more than one type of co-charged species, their force balances (between FD, proportional to the size, and FE, proportional to the charge) will be different. As a result, they will be focused at different locations in the microchannel. If one applies a pressure-driven flow to make the fluid drag force FD of a charged molecule greater than the maximum FE at the boundary of IDZ, this molecule may squeeze through IDZ and flow out of the microchannel. In the meantime, other co-charged molecules remain trapped at the boundary of IDZ because their fluid drag forces are lower than the maximum electric forces at that location. In this way, selective trapping and separation of co-charged species are realized. Based on this understanding, Ouyang et al. selectively focused DNA molecules in the microchannel but let the proteins pass through [30]. Gong et al. studied the applicability of this mechanism for the separation of Li+ and Mg2+ through numerical simulation. They achieved a very low Mg2+/Li+ ratio in their product solution, but Li+ concentration is also very low [31,32]. Later the same group proposed a method to collect Li+ enriched solution selectively through a branch channel, achieving simultaneous concentration of Li+ and removal of Mg2+ [32].
In recent work, Papadimitridou et al. reported a new system to focus (and separate) charged particles [33]. In their system, sample solutions were pumped horizontally from left to right (as in Figures of [33]) through a macroscale chamber, with all the chamber walls formed by large arrays of microchannels. In the meantime, an electric field was applied in the perpendicular direction (vertically from top to bottom as in Figures of [33]). Nafion membranes are embedded across the microchannel arrays outside the bottom boundary of the chamber to facilitate the cation-selective transport and generate the ICP effect there. The vertical electric field also produces an electroosmotic flow from top to bottom, which carries all the species in the solution to move downwards and facilitates focusing of negatively charged particles while migrating from left to right along with the pressure-driven flow. Based on the fact that the solution runs continuously through the macroscale chamber, with charged species focused (through ICP) and collected at the bottom right corner, this device was referred to as a free flow ICP focusing (FF-ICPF) system. Their FF-ICPF device operates in two modes: peak mode (the concentration of the focused analytes is significantly lower than that of buffer ions) and plateau mode (concentrations of the focused analytes are high enough to neutralize the counterions in the buffer). Separation of different species can be performed in the peak mode by collecting the solutions at different groups of microchannels at the outlet (at the right boundary of the chamber). The FF-ICPF system in [28] can process the sample at a flow rate of 10 µL/min, significantly faster than typical microfluidic systems. However, its enrichment factor is low: ~17 times for three dyes in peak mode and 4–5 times in plateau mode, which is orders lower than the single-channel devices proposed by Wang et al. [23]. In fact, concentrations of focused analytes in the peak mode and the plateau mode had been termed as “electrokinetic limit” and “neutrality limit” by Ouyang et al. in their theoretical work, where approximated analytical expressions and scaling laws were given rigorously [34].
Now that the FF-ICP system has the potential to focus and separate the charged species at a macroscopic scale, one may ask if such a system is capable of concentration and separation of simple ions. This paper will answer this question through numerical simulation. More specifically, we will evaluate the performance of the FF-ICP system in separating Li+ and Mg2+ ions and check its potential role in lithium extraction from high Mg2+/Li+ ratio salt lake brines.

2. Method

2.1. Physical Setup

Similar to the device in [33], the microfluidic device we study is composed of a chamber surrounded by four microchannels arrays. Surfaces of the microchannels at the upper and lower sides of the chamber are positively charged. Under the microchannels at the upper side of the chamber, an anion exchange membrane (AEM) is embedded. Mixed Li+, Na+, Mg2+, K+, and Cl ions flow from left to right, driven by external pressure.
In the meantime, microchannels at the upper and lower sides of the chamber are connected to reservoirs of NaCl solution. A strong electric field is applied from top to bottom. This electric field will apply an upward FE to anions and a downward FE to cations and produce IDZ near the AEM at the interfaces between microchannels and the upper side of the chamber. This electric field also generates an upward electroosmotic flow, which applies an additional component to FD of all the species in the chamber, on top of the fluid drag imposed by the rightward pressure-driven flow (PDF). As a result, when a cation, e.g., Li+ or Mg2+, enters the chamber from the left side, it is subjected to an up-right directed FD and a downward FE. If the drag force in the y-direction FD,y is greater than the y-component of the electric force FE,y, the ion will move upward when moving rightward with the fluid.
On the other hand, if FD,y is smaller than FD,y, the ion moves downward. If the y-components of two forces are equal, the ion will be balanced in that location. Because of the ICP effect, the electric field at the boundary of IDZ is significantly higher than in the other regions. These uneven distributions of electric force define different equilibrium positions in y-direction for different cation species. As shown in Figure 1a, as Li+ ion enters the chamber at the center of the inlet, it is subject to a larger FD,y than FE,y (because the downward electric field there is low). Hence, it will move upward as it is moving with fluid in the right direction. When it is closer to the IDZ region, the electric field will become stronger, and the electric force will increase. This upward motion will stop when the two forces (FD,y and FE,y) are equal. It will move rightward and run out through microchannels at the right side of the chamber. In the meantime, Mg2+ is subject to a 2-times higher electric field and a 1.46-times (calculated from the diffusion coefficients shown in Table 1) higher fluid drag force as compared to those of Li+, the balance position for Mg2+ will be lower than that of Li+. If we optimize the fluid flow speeds and the electric field such that the focused peaks of Li+ and Mg2+ run out of the system through a different group of microchannels, these two ion species are separated.

2.2. Simulation Model

We use a simplified model for our simulation study, as shown in Figure 1b. In this model, the chamber is of length L and width H. At the left and right sides of the chambers, there are a number of microchannels of length Lc and width Hc. On the upper and lower boundaries of the chamber, there are a series of barriers of length Lb (representing the barrier between microchannels in Figure 1a), the distance between, which is Lm (representing the microchannels). We apply a constant inlet flow of u1 and a constant outlet flow of speed u2 at the location of the microchannels on the lower and upper boundaries, respectively. This treatment is based on the fact that the electroosmotic flow has the plug shape in a microchannel, and the concentration of the solution in the lower and upper sides is different. At the inlets of microchannels at the left side of the chamber, constant pressure of P1 is applied. At the outlet of microchannels on the right side, a constant flow speed uout is set.
To simplify our computation model, we move the AEM embedded to the microchannels at the upper side of the chamber to the chamber boundary (the green sections). For the electric field, a positive electric potential V is applied on these membranes’ segments. At the lower boundary of the chamber, the segments corresponding to the microchannels are grounded.

2.3. Governing Equations

For the fluid flow, the Navier–Stokes equation for an incompressible Newton fluid is employed.
ρ ( U t + ( U ) U ) = P + η 2 U ρ e Φ
ρ U = 0
Here ρ is the fluid mass density, U is the fluid velocity, P is the pressure, and Φ is electric potential.
Transport of charged species is modeled via the Nernst–Planck equations and the mass conservation equation.
J i = D i C i Z i ( D i F / R T ) C i Φ + U C i
J i = C i t
where J i , C i and Z i are flux, concentration, and valence, all for species i . For convenience, we set i = 1 for Li+, i = 2 for Na+, and i = 3, 4, 5 for Mg2+, K+, and Cl, respectively. Symbols F , R , and the parameter T are the Faraday constant, gas constant, and temperature, respectively. The Poisson equation governed the distribution of the electric field.
( ε Φ ) = ρ e
where ρ e = e Z i C i is the charge density, with e representing the elementary charge and ε denoting the dielectric permittivity of the solution.

2.4. Boundary Conditions

As shown in Figure 1b, the boundary conditions illustrate as follows:
On the inlet boundaries of microchannels at the left side, (i) the pressure is P1; (ii) the concentration of all ionic species is equal to those in the reservoir:
P = P 1 ,   C i = C i , 0 ,   i = 1 ,   2 ,   ,   5
On the outlet boundaries of microchannels at the right side, (i) the free boundary condition is applied for mass transport; (ii) the fluid flow speed is the constant u out = ( u out , 0 ) :
U = u out ,   C i n = 0 ,   i = 1 ,   2 ,   ,   5
where n represents the outer-pointing normal vector perpendicular to the boundary of the fluid domain.
On the microchannel portion of the lower boundary, (i) the electric potential of the lower boundary is 0; (ii) the velocity of the fluid is constant u 1 = ( 0 , u 1 ) ; (iii) the concentration of Na+ is C2,1, the concentration of Cl is the opposite number of the sum of other four cations:
Φ = 0 ,   U = u 1 ,   C 2 = C 2 , 1 ,   C 5 = C i Z i ,   i = 1 ,   2 ,   3 ,   4
where c i and Z i represent the variable concentration and the valence of ion i.
On the red dashed boundary, (i) a voltage V is applied on the membrane to generate electric field; (ii) the velocity of the fluid is the constant velocity ( u 2 = ( 0 , u 2 ) ); (iii) the concentration of Cl at the upper boundary is C m ; (iv) fluxes of cations across the membrane are zero. The corresponding equations can be expressed as:
Φ = V ,   U = u 2 ,   C 5 = C m ,   J i n = 0 ,   i = 1 ,   2 ,   3 ,   4
The other black lines and black dashed lines are the microchannel walls, (i) no-slip condition for fluid velocity; (ii) impermeability to all anions and cations:
U = 0 ,   J i n = 0 ,   i = 1 ,   2 ,   ,   5

2.5. Numerical Methods

Simulations were carried out using COMSOL Multiphysics software (version 5.6) on a Dell workstation (Precision 7920) equipped with an Intel Xeon processor (Gold 6128) and 112 GB of RAM. Steady-state simulations were used for all studies. Solution convection was modeled with the “Creeping Flow” interface. Moreover, the “Transport of Diluted Species” and the “Electrostatics” interfaces were coupled to solve the “Nernst–Planck–Poisson” equation; 3,833,759 quadrilateral elements were utilized for meshing. Near the membrane region, extremely fine meshes are used to ensure sufficient solution accuracy. To obtain the highly nonlinear solution under high electric potential, we needed to start from low voltage and sweep the high voltage parameter. Initially, the down boundary should be set as no flux for Li+, Mg2+, and K+ to converge the simulation. After the system is stable, we set it out to facilitate ion exchange between the chamber and the reservoir below.
As shown in Figure 1b, in the simulation model, we kept the horizontal microfins to analyze the average ion fluxes of the outflow because the value taken on the boundary will greatly affect the calculation result. However, we ignore the vertical ones as they do not affect the results.
To analyze the behavior of the separation system, we started with the setting of a particular parameter. Then, we studied the effects of two critical operational parameters, the voltage V and velocity u, and clarified how these parameters affected system performance to prove the feasibility of the proposed ion separation method.

3. Results and Discussions

In this simulation, the geometric parameters of the chamber are L = 90 μm, H = 31.5 μm (see Figure 1b). The length of the microfins channel is Lc = 3 μm, and the width is Hc = 1.5 μm, while the distance between the channels is Hb = 1.5 μm. The length of membrane segments is Lm = 1.5 μm, Lb = 1.5 μm.
Using a simplified brine consisting of only five ions (Li+, Na+, K+, Mg2+, Cl) [31]. After diluting the raw brine, we selected the following typical concentrations. The concentrations of ions (the left inlet boundary) are: C 1 , 0 = 0.001   mM , C 2 , 0 = 0.125   mM , C 3 , 0 = 0.04   mM , C 4 , 0 = 0.04   mM , C 5 , 0 = 0.211   mM . The ion concentrations of the lower edge are C 2 , 1 = 0.375 mM, C 5 , 1 = C i Z i ,   i = 1 ,   2 ,   3 ,   4 , where NaCl enters the chamber from the lower boundary as a supplement buffer to maintain electrical neutrality. In the ideal simplified model of the ion-selective membrane, the results of fixed voltage and fixed counter ion concentration are accurate in most cases, especially in the case of high voltage and/or high charge density [36]. At the membrane boundary, the voltage V is from V = 0 to V = 30VT, with the thermal voltage VT equal to 25.8 mV [31]. The assumption of charge selectivity is a fixed counterion boundary. We set the fixed concentration of Cl at C m = 10 C 5 , 0 [37]. As for other parameters, we set P1 = 100 Pa, T = 300 K, ρ = 1000 kg/m3, and η = 0.001 Pa⋅s, respectively.
The concentration of Na+ is not an essential research item. Although the content of this ion in the brine is high, the Li+ can be collected in the subsequent steps [31]. Therefore, it is easy to separate Li+ ions from it. Instead, we have to focus on the extraction quantity of Li+ ions, especially the Mg2+/Li+ flux ratio in the microchannel and use this parameter as a separate indicator. Steady-state solutions are used to study the performance of the system. The governing Equations (1)–(5) are solved with the above boundary conditions (6)–(10).
We simulate the steady-state behavior of this system at V = 30VT, and the PDF constrains the fluid velocity u out = 0.5 mm/s in the x-direction, and the velocity rate of y-direction is u 1 = 0.63 mm/s and u 2 = 1 mm/s.
Comparing the two pictures (Figure 2a,b), we can find that in the black dashed box range (18 μm ≤ y ≤ 31.5 μm), the Li+ ions concentration was much higher than the feed. We put the collector there to get a high-resolution separation effect. The box indicates the location of the collector.
Since the initial concentrations were different, we paid more attention to the flux of the two ions. Most Li+ flow through the upper right of the chamber where the collector box is located (Figure 3). The flux of Mg2+ at the right outlet of the chamber is much lower than the feed (Figure 4). These two figures demonstrated that the proposed system could continuously extract Li+ ions with a concentration of 1.28 times in the raw brine while simultaneously removing Mg2+ ions.

3.1. Effect of the Voltage

To reduce the vortices, we control the voltage within 30VT. It is noticed that only a small fraction of current across the horizontal direction. Hence, the electric field distribution is zero in that direction. When the voltage increases, the IDZ area becomes more extensive, causing the cations to move downwards. Since the FE of Mg2+ ions is significantly greater than other ions, most move to the downside reservoir. The effect of the Mg2+/Li+ separation is significant before reaching the unstable fluid motion caused by the electric instability (Figure 5) [38,39,40,41,42].

3.2. Effect of the Velocity uout

As shown in Figure 6, when the horizontal speed uout increases, the Mg2+/Li+ flux ratio increased. Although the flux of Li+ ions increases, there is not enough time for Mg2+ to reach the downside reservoir, resulting in a significant rise in the average flux of Mg2+ and poor separation effect, specifically reflected in the increase in the Mg2+/Li+ flux ratio.

3.3. Effect of the Velocity u1

When the speed u1 increases, as shown in Figure 7, in this case, the separation effect becomes better due to the competition of the co-ions [32]. The Na+ supplemented from below increase, the electrophoretic mobility of Na+ is between the Li+ and Mg2+, it would behave as an electrophoretic spacer.

3.4. Effect of the Velocity u2

The velocity u2 plays a crucial role in this system. As shown in Figure 8, when the u2 increases, the enrichment effect of Li+ ions becomes better. As Li+ ions have the lowest electrophoretic mobility, they have the largest upward FD than the downward FE. Therefore, they tend to focus near the AEM. Ions with higher valence, i.e., Mg2+, are driven to the downside reservoir by the FE.

3.5. Advantages and Limitations

The advantage of this system is that the analytes are continuously flowed through the system by a PDF while vertically separated and concentrated by ICP. Compared to other focusing systems [43], it has a simple structure. More significantly, this method does not involve the chemical reaction of ions in the brine. Compared with other methods such as nanofiltration membrane [14]. Li et al. increased the Li+ yield and improved the separation effect of Mg2+ and Li+ by increasing the operating pressure. The Mg2+ rejection rate of 92% was observed at a high-pressure level. Nevertheless, membrane fouling often occurs in nanofiltration separation, and after prolonging the operation time, the separation efficiency is reduced [31]. The typical channel size of our system is the order of tens of microns, so it will not suffer from clogging problems. Zhang et al. added new material to the extraction system of polymer inclusion membrane to extract Li+ from the Mg2+ rich solution (the initial Mg2+/Li+ molar ratio is 15) [44]. They have increased the extraction rate of Li+ by 20%. However, it is necessary to reduce the loss of chemicals additive further to enhance the stability of the membrane. Our system can operate stably without worrying about material loss.
However, some limitations of this work should be noted. First of all, this work is based on first-principles numerical simulation (no experience/parameter approximation). Therefore, the proposed system should be finally verified through experiments. Secondly, in the conventional electro-membrane process, high current density usually produces various nonlinear effects (permeation selectivity loss, heating, over-limit behavior, water splitting effect) [38,45,46,47]. When the electric convection operates in the overcurrent mode, it will form near the ion-selective membrane. The size of these vortices increases with electric potential, which eventually leads to fluid instability. The system can only operate under a certain electric field strength to maintain stability.

4. Conclusions

We have proved the feasibility of continuous separation of Mg2+ and Li+ ions through numerical simulation based on FF-ICPF. This method could perform enrichment perpendicular to the flow direction of the analyte. Furthermore, the maximum Mg2+/Li+ ratio can be decreased by ~85% under the electroneutrality limit. However, after trying various channel structures and parameters, the enrichment factor cannot be increased, which is the inherent disadvantage of this system. Numerical simulation shows that this limiting behavior is confined by the accumulation of the charged particles, which affects the shape of the electric field gradient. In the future, if the Joule heating is added, the simulation would be more accurate.

Author Contributions

D.Z.: investigation, writing–original draft, writing–review and editing. X.Z.: writing—review and editing. L.X.: data curation. Z.L.: methodology, funding acquisition, resources. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (Grant Nos. 12072100 and 21576130), and the Graduate Scientific Research Foundation of Wenzhou University, China.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Averill, W.A.; Olson, D.L. A Review of extractive processes for lithium from ores and brines. In Lithium Needs and Resources; Penner, S.S., Ed.; Pergamon: Oxford, UK, 1978; pp. 305–313. ISBN 978-0-08-022733-7. [Google Scholar]
  2. Kesler, S.E.; Gruber, P.W.; Medina, P.A.; Keoleian, G.A.; Everson, M.P.; Wallington, T.J. Global Lithium resources: Relative Importance of pegmatite, brine and other deposits. Ore Geol. Rev. 2012, 48, 55. [Google Scholar] [CrossRef]
  3. Swain, B. Recovery and recycling of lithium: A review. Sep. Purif. Technol. 2017, 172, 388–403. [Google Scholar] [CrossRef]
  4. Yu, J.; Zheng, M.; Wu, Q.; Nie, Z.; Bu, L. Extracting lithium from Tibetan Dangxiong Tso Salt Lake of carbonate type by using geothermal salinity-gradient solar pond. Sol. Energy 2015, 115, 133–144. [Google Scholar] [CrossRef]
  5. An, J.W.; Kang, D.J.; Tran, K.T.; Kim, M.J.; Lim, T.; Tran, T. Recovery of lithium from Uyuni salar brine. Hydrometallurgy 2012, 117, 64–70. [Google Scholar] [CrossRef]
  6. Xu, Z.; Zhang, H.; Wang, R.; Gui, W.; Liu, G.; Yang, Y. Systemic and direct production of battery-grade lithium carbonate from a saline lake. Ind. Eng. Chem. Res. 2014, 53, 16502–16507. [Google Scholar] [CrossRef]
  7. Zhang, Y.; Hu, Y.; Wang, L.; Sun, W. Systematic review of lithium extraction from salt-lake brines via precipitation approaches. Miner. Eng. 2019, 139, 105868. [Google Scholar] [CrossRef]
  8. Li, X.; Mo, Y.; Qing, W.; Shao, S.; Tang, C.Y.; Li, J. Membrane-based technologies for lithium recovery from water lithium resources: A review. J. Memb. Sci. 2019, 591, 117317. [Google Scholar] [CrossRef]
  9. Xiao, G.; Tong, K.; Zhou, L.; Xiao, J.; Sun, S.; Li, P.; Yu, J. Adsorption and Desorption behavior of lithium ion in spherical PVC–MnO2 Ion sieve. Ind. Eng. Chem. Res. 2012, 51, 10921–10929. [Google Scholar] [CrossRef]
  10. Xu, W.; He, L.; Zhao, Z. Lithium extraction from high Mg/Li brine via electrochemical intercalation/de-intercalation system using LiMn2O4 materials. Desalination 2021, 503, 114935. [Google Scholar] [CrossRef]
  11. Li, X.; Chao, Y.; Chen, L.; Chen, W.; Luo, J.; Wang, C.; Wu, P.; Li, H.; Zhu, W. Taming wettability of lithium ion sieve via different TiO2 precursors for effective Li recovery from aqueous lithium resources. Chem. Eng. J. 2020, 392, 123731. [Google Scholar] [CrossRef]
  12. Nie, X.Y.; Sun, S.Y.; Sun, Z.; Song, X.; Yu, J.G. Ion-fractionation of lithium ions from magnesium ions by electrodialysis using monovalent selective ion-exchange membranes. Desalination 2017, 403, 128–135. [Google Scholar] [CrossRef]
  13. Yang, G.; Shi, H.; Liu, W.; Xing, W.; Xu, N. Investigation of Mg2+/Li+ separation by nanofiltration. Chin. J. Chem. Eng. 2011, 19, 586–591. [Google Scholar] [CrossRef]
  14. Li, Y.; Zhao, Y.J.; Wang, H.; Wang, M. The application of nanofiltration membrane for recovering lithium from salt lake brine. Desalination 2019, 468, 114081. [Google Scholar] [CrossRef]
  15. Flexer, V.; Baspineiro, C.F.; Galli, C.I. Lithium recovery from brines: A vital raw material for green energies with a potential environmental impact in its mining and processing. Sci. Total Environ. 2018, 639, 1188–1204. [Google Scholar] [CrossRef]
  16. Gao, F.; Zheng, M.-P.; Nie, Z.; Liu, J.-H.; Song, P.-S. Brine lithium resource in the salt lake and advances in its exploitation. Acta Geosci. Sin. 2011, 32, 483–492. [Google Scholar]
  17. Wang, J.; Yang, S.; Bai, R.; Chen, Y.; Zhang, S. Lithium Recovery from the Mother Liquor Obtained in the Process of Li2CO3 Production. Ind. Eng. Chem. Res. 2019, 58, 1363–1372. [Google Scholar] [CrossRef]
  18. Gu, D.; Sun, W.; Han, G.; Cui, Q.; Wang, H. Lithium ion sieve synthesized via an improved solid state method and adsorption performance for West Taijinar Salt Lake brine. Chem. Eng. J. 2018, 350, 474–483. [Google Scholar] [CrossRef]
  19. Ji, P.-Y.; Ji, Z.-Y.; Chen, Q.-B.; Liu, J.; Zhao, Y.-Y.; Wang, S.-Z.; Li, F.; Yuan, J.-S. Effect of coexisting ions on recovering lithium from high Mg2+/Li+ ratio brines by selective-electrodialysis. Sep. Purif. Technol. 2018, 207, 1–11. [Google Scholar] [CrossRef]
  20. Ji, Z.Y.; Zhao, M.Y.; Zhao, Y.Y.; Liu, J.; Peng, J.L.; Yuan, J.S. Lithium extraction process on spinel-type LiMn2O4 and characterization based on the hydrolysis of sodium persulfate. Solid State Ion. 2017, 301, 116–124. [Google Scholar] [CrossRef]
  21. Guo, X.; Hu, S.; Wang, C.; Duan, H.; Xiang, X. Highly Efficient Separation of Magnesium and Lithium and High-Valued Utilization of Magnesium from salt lake brine by a reaction-coupled separation technology. Ind. Eng. Chem. Res. 2018, 57, 6618–6626. [Google Scholar] [CrossRef]
  22. Ji, Z.Y.; Chen, Q.B.; Yuan, J.S.; Liu, J.; Zhao, Y.Y.; Feng, W.X. Preliminary study on recovering lithium from high Mg2+/Li+ ratio brines by electrodialysis. Sep. Purif. Technol. 2017, 172, 168–177. [Google Scholar] [CrossRef]
  23. Wang, Y.C.; Stevens, A.L.; Han, J. Million-fold preconcentration of proteins and peptides by nanofluidic filter. Anal. Chem. 2005, 77, 4293–4299. [Google Scholar] [CrossRef]
  24. Kim, S.J.; Li, L.D.; Han, J. Amplified electrokinetic response by concentration polarization near nanofluidic channel. Langmuir 2009, 25, 7759–7765. [Google Scholar] [CrossRef] [Green Version]
  25. Mani, A.; Zangle, T.A.; Santiago, J.G. On the propagation of concentration polarization from microchannel-nanochannel interfaces Part I: Analytical model and characteristic analysis. Langmuir 2009, 25, 3898–3908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kim, S.J.; Song, Y.A.; Han, J. Nanofluidic concentration devices for biomolecules utilizing ion concentration polarization: Theory, fabrication, and applications. Chem. Soc. Rev. 2010, 39, 912–922. [Google Scholar] [CrossRef] [Green Version]
  27. Kim, S.J.; Ko, S.H.; Kang, K.H.; Han, J. Direct seawater desalination by ion concentration polarization. Nat. Nanotechnol. 2010, 5, 297–301. [Google Scholar] [CrossRef]
  28. Hlushkou, D.; Knust, K.N.; Crooks, R.M.; Tallarek, U. Numerical simulation of electrochemical desalination. J. Phys. Condens. Matter 2016, 28, 194001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Gong, L.; Ouyang, W.; Li, Z.; Han, J. Force fields of charged particles in micro-nanofluidic preconcentration systems. AIP Adv. 2017, 7, 125020. [Google Scholar] [CrossRef] [Green Version]
  30. Ouyang, W.; Li, Z.; Han, J. Pressure-Modulated selective electrokinetic trapping for direct enrichment, purification, and detection of nucleic acids in human serum. Anal. Chem. 2018, 90, 11366–11375. [Google Scholar] [CrossRef]
  31. Gong, L.; Ouyang, W.; Li, Z.; Han, J. Direct numerical simulation of continuous lithium extraction from high Mg2+/Li+ ratio brines using microfluidic channels with ion concentration polarization. J. Membr. Sci. 2018, 556, 34–41. [Google Scholar] [CrossRef] [PubMed]
  32. Gong, L.; Li, Z.; Han, J. Numerical simulation of continuous extraction of highly concentrated Li+ from high Mg2+/Li+ ratio brines in an ion concentration polarization-based microfluidic system. Sep. Purif. Technol. 2019, 217, 174–182. [Google Scholar] [CrossRef]
  33. Papadimitriou, V.A.; Segerink, L.I.; Eijkel, J.C.T. Free Flow ion concentration polarization focusing (FF-ICPF). Anal. Chem. 2020, 92, 4866–4874. [Google Scholar] [CrossRef] [Green Version]
  34. Ouyang, W.; Ye, X.; Li, Z.; Han, J. Deciphering ion concentration polarization-based electrokinetic molecular concentration at the micro-nanofluidic interface: Theoretical limits and scaling laws. Nanoscale 2018, 10, 15187–15194. [Google Scholar] [CrossRef] [Green Version]
  35. Vanýsek, P. Ionic Conductivity and Diffusion at Infinite Dilution, 1992/93 ed.; CRC Press: Boca Raton, FL, USA, 1992. [Google Scholar]
  36. Jiang, J.; Tang, J.; Al-Anzi, B.; Han, J.; Li, Z. On the validity of ion selective membrane simplification in concentration polarization. AIP Adv. 2021, 11, 035116. [Google Scholar] [CrossRef]
  37. Rubinstein, I.; Zaltzman, B. Equilibrium electroconvective instability. Phys. Rev. Lett. 2015, 114, 114502. [Google Scholar] [CrossRef] [Green Version]
  38. Nanthasurasak, P.; Cabot, J.M.; See, H.H.; Guijt, R.M.; Breadmore, M.C. Electrophoretic separations on paper: Past, present, and future—A review. Anal. Chim. Acta 2017, 985, 7–23. [Google Scholar] [CrossRef]
  39. Kim, S.J.; Ko, S.H.; Kwak, R.; Posner, J.D.; Kang, K.H.; Han, J. Multi-vortical flow inducing electrokinetic instability in ion concentration polarization layer. Nanoscale 2012, 4, 7406–7410. [Google Scholar] [CrossRef]
  40. Druzgalski, C.L.; Andersen, M.B.; Mani, A. Direct numerical simulation of electroconvective instability and hydrodynamic chaos near an ion-selective surface. Phys. Fluids 2013, 25, 110804. [Google Scholar] [CrossRef]
  41. Karatay, E.; Druzgalski, C.L.; Mani, A. Simulation of chaotic electrokinetic transport: Performance of commercial software versus custom-built direct numerical simulation codes. J. Colloid Interface Sci. 2015, 446, 67–76. [Google Scholar] [CrossRef] [PubMed]
  42. Park, S.; Kwak, R. Microscale electrodeionization: In situ concentration profiling and flow visualization. Water Res. 2020, 170, 115310. [Google Scholar] [CrossRef] [PubMed]
  43. Phan, D.T.; Jin, L.; Wustoni, S.; Chen, C.H. Buffer-free integrative nanofluidic device for real-time continuous flow bioassays by ion concentration polarization. Lab Chip 2018, 18, 574–584. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, C.; Mu, Y.; Zhao, S.; Zhang, W.; Wang, Y. Lithium extraction from synthetic brine with high Mg2+/Li+ ratio using the polymer inclusion membrane. Desalination 2020, 496, 114710. [Google Scholar] [CrossRef]
  45. Yoon, J.; Do, V.Q.; Pham, V.S.; Han, J. Return flow ion concentration polarization desalination: A new way to enhance electromembrane desalination. Water Res. 2019, 159, 501–510. [Google Scholar] [CrossRef] [PubMed]
  46. De Valença, J.; Jõgi, M.; Wagterveld, R.M.; Karatay, E.; Wood, J.A.; Lammertink, R.G.H. Confined Electroconvective vortices at structured ion exchange membranes. Langmuir 2018, 34, 2455–2463. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Liu, W.; Gong, L.; Zhu, Y.; Li, Z. Augmented electroosmotic flow and simultaneous desalination in microchannels embedded with permselective membranes. Sci. Sin. Technol. 2018, 48, 17–24. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of FF-ICP system for ion concentration and separation. (b) The simulation model.
Figure 1. (a) Schematic of FF-ICP system for ion concentration and separation. (b) The simulation model.
Membranes 11 00697 g001
Figure 2. (a) The concentration distribution of Li+ in the chamber (b) The concentration distribution of Mg2+ in the chamber. Values of other parameters: V = 30VT, uout = 0.5 mm/s, u1 = 0.63 mm/s, and u2 = 1 mm/s.
Figure 2. (a) The concentration distribution of Li+ in the chamber (b) The concentration distribution of Mg2+ in the chamber. Values of other parameters: V = 30VT, uout = 0.5 mm/s, u1 = 0.63 mm/s, and u2 = 1 mm/s.
Membranes 11 00697 g002
Figure 3. The flux distribution of Li+. Values of other parameters are: V = 30 VT, uout = 1 mm/s, u1 = 0.63 mm/s, and u2 = 1 mm/s.
Figure 3. The flux distribution of Li+. Values of other parameters are: V = 30 VT, uout = 1 mm/s, u1 = 0.63 mm/s, and u2 = 1 mm/s.
Membranes 11 00697 g003
Figure 4. The flux distribution of Mg2+. Values of other parameters are: V = 30VT, uout = 1 mm/s, u1 = 0.63 mm/s, and u2 = 1 mm/s.
Figure 4. The flux distribution of Mg2+. Values of other parameters are: V = 30VT, uout = 1 mm/s, u1 = 0.63 mm/s, and u2 = 1 mm/s.
Membranes 11 00697 g004
Figure 5. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on V. Values of other parameters are: V = 30VT, uout = 1 mm/s, u1 = 0.63 mm/s, u2 =1 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Figure 5. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on V. Values of other parameters are: V = 30VT, uout = 1 mm/s, u1 = 0.63 mm/s, u2 =1 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Membranes 11 00697 g005
Figure 6. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on uout. Values of other parameters are: V = 30VT, u1 = 0.63 mm/s, u2 =1 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Figure 6. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on uout. Values of other parameters are: V = 30VT, u1 = 0.63 mm/s, u2 =1 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Membranes 11 00697 g006
Figure 7. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on u1. Values of other parameters are: V = 30VT, uout = 0.5 mm/s, u2 = 1 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Figure 7. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on u1. Values of other parameters are: V = 30VT, uout = 0.5 mm/s, u2 = 1 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Membranes 11 00697 g007
Figure 8. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on u2. Values of other parameters are: V = 30VT, uout = 0.5 mm/s, u1 = 0.63 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Figure 8. The dependence of the average flux of Li+, the average flux of Mg2+ and the Mg2+/Li+ flux ratio on u2. Values of other parameters are: V = 30VT, uout = 0.5 mm/s, u1 = 0.63 mm/s, 18 μm ≤ y ≤ 31.5 μm, x = 96 μm.
Membranes 11 00697 g008
Table 1. The transport parameters of the ions in the raw brine [35].
Table 1. The transport parameters of the ions in the raw brine [35].
Index iSpeciesDiffusion Coefficient Di (×10−9 m2/s)Electrophoretic Mobility μi (×10−8 m2/V⋅s)
1Li+1.0293.98
2Na+1.3345.152
3Mg2+0.7067.563
4K+1.9575.457
5Cl2.0327.853
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, D.; Zhang, X.; Xing, L.; Li, Z. Numerical Simulation of Continuous Extraction of Li+ from High Mg2+/Li+ Ratio Brines Based on Free Flow Ion Concentration Polarization Microfluidic System. Membranes 2021, 11, 697. https://doi.org/10.3390/membranes11090697

AMA Style

Zhang D, Zhang X, Xing L, Li Z. Numerical Simulation of Continuous Extraction of Li+ from High Mg2+/Li+ Ratio Brines Based on Free Flow Ion Concentration Polarization Microfluidic System. Membranes. 2021; 11(9):697. https://doi.org/10.3390/membranes11090697

Chicago/Turabian Style

Zhang, Dongxiang, Xianglei Zhang, Leilei Xing, and Zirui Li. 2021. "Numerical Simulation of Continuous Extraction of Li+ from High Mg2+/Li+ Ratio Brines Based on Free Flow Ion Concentration Polarization Microfluidic System" Membranes 11, no. 9: 697. https://doi.org/10.3390/membranes11090697

APA Style

Zhang, D., Zhang, X., Xing, L., & Li, Z. (2021). Numerical Simulation of Continuous Extraction of Li+ from High Mg2+/Li+ Ratio Brines Based on Free Flow Ion Concentration Polarization Microfluidic System. Membranes, 11(9), 697. https://doi.org/10.3390/membranes11090697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop