Abstract
To overcome a world-wide water shortage problem, numerous desalination methods have been developed with state-of-the-art power efficiency. Here we propose a spontaneous desalting mechanism referred to as the capillarity ion concentration polarization. An ion-depletion zone is spontaneously formed near a nanoporous material by the permselective ion transportation driven by the capillarity of the material, in contrast to electrokinetic ion concentration polarization which achieves the same ion-depletion zone by an external d.c. bias. This capillarity ion concentration polarization device is shown to be capable of desalting an ambient electrolyte more than 90% without any external electrical power sources. Theoretical analysis for both static and transient conditions are conducted to characterize this phenomenon. These results indicate that the capillarity ion concentration polarization system can offer unique and economical approaches for a power-free water purification system.
Similar content being viewed by others
Introduction
Recent noteworthy investigations in micro/nanofluidic devices enable to develop not only various biomedical but also significant energy and environmental applications1. Among these groundbreaking technologies, micro/nanofluidic platforms have successfully served as the differential shells of macro-scale water purification/desalination systems2,3,4. One can microscopically investigate the salt separation mechanism occurring near nanoporous membranes or porous electrodes by flow field tracking or voltage–current measurements. Although state-of-the-art desalination methods, such as reverse-osmosis and flash distillation, can provide the gigantic amount of fresh water with high energy efficiency5, low energy efficiency of small-scale platforms and the requirement of high-cost infrastructures have become acute nuisances, especially for remote/rural areas of the world. On the other hand, the micro/nanofluidic differential shells can be scalably integrated to a portable-scale system for individual usage or even to a mid-scale system for small village supply.
For such small-sized desalination/purification applications, ion concentration polarization (ICP) has been recently introduced2,6,7,8. ICP indicates the polarization of electrolytic concentrations at the anodic and cathodic side of a permselective membrane under d.c. bias. Typical behaviour is that the concentration becomes largely depleted at the anodic side (also known as an ion-depletion zone) and enriched at the cathodic side (also known as an ion-enrichment zone) in the case with a cation-selective membrane and vice versa9. In addition to nanoporous membranes10,11,12, bipolar electrodes13,14 and nanoporous particles15,16 have been reported for a new practical platform to generate ICP as well. Intensive researches have been conducted to explore its fundamental aspects such as the source of over-limiting current17,18,19, the vortex generation20,21,22 and related instability issues23 and its novel engineering applicability such as biological sample preparations24,25, energy saving mechanisms26 and electro-desalination applications2,14. Although ICP desalination method is unprecedentedly capable of extracting fresh water from highly saline, biologically infected or even heavy metal contaminated water, the necessity of a high external electrical power source4,27 is still the most critical hurdle to be overcome. Therefore, a new desalination mechanism that consumes the minimum level of energy is highly demanded.
In nature, mangrove trees that grow in saline water along the tropical or subtropical coasts are well known for their unique physiological responses to survive in highly saline environments without any external supply of mechanical, chemical or electrical powers28. The phenomenological observations of mangrove28 so far have found that the root is able to provide sap of 90% lower salinity than the external saline water by a reverse osmosis mechanism. Therefore, we are particularly interested in the ability of mangrove to exclude salt without any external energy. Motivated by this fact, we suggest a new spontaneous desalination mechanism by introducing the nanoscale electrokinetic and the hydrodynamic transportation of ions at the interfaces between the nanoporous structure and the saline environment. We visualize the ion movements that illustrate the new class of the ion-depletion phenomenon. By utilizing dry ionic hydrogel, the effective ionic flux through the hydrogel is induced by the capillarity of the hydrogel itself, instead of an external electrical bias, so that we name the currently exploited process as the capillarity ion concentration polarization (CICP). Osmo-poro-elastic process due to the swelling of the hydrogel29,30 is excluded in this work by a mechanical confinement31,32,33 within a microchannel to focus on the role of capillarity. Both analytical and numerical models are obtained to analyse the statics and dynamics of the spontaneous desalting mechanism. Since the mechanism is independent from any electrical power sources, it offers significant merits in terms of a power consumption and a stable operation compared with the conventional ICP mechanism. Therefore, CICP mechanism can be applied to power-free and small-scale (possibly portable) desalination/purification devices that would be particularly useful in remote/rural and disaster-stricken areas.
Results
Imbibition rates through the ionic hydrogel
The ionic hydrogel is capable of absorbing ambient aqueous solution by the capillarity. This phenomenon is critical to the presented CICP system, which utilizes the capillarity as a driving force for a liquid transportation. A hydrogel usually has an osmo-poro-elastic effect due to a significant swelling. However, the swelling has been reported to be minimized by a mechanical confinement31,32,33. In this work, the hydrogel swells within a microchannel and, thus, we are able to analyse the CICP phenomenon only by the capillarity effect. See Supplementary Note 1 for the swelling result either with or without mechanical confinements. As shown in Fig. 1a, there are two significant evidences of the imbibition by the capillarity. Micro-particles are injected near the hydrogel pad and their motions are captured at distinct time intervals. As indicated, their paths converge toward the hydrogel pad from both reservoirs in the centre-connection device, evidencing a mass sink at the centre. See Supplementary Video 1. In addition, the particle convergence in the end-connection device is also shown in Supplementary Note 2. The descriptions of the centre- and the end-connection device are given in the ‘Methods’ section. Since there are no external driving forces (such as gravity or external electric/magnetic/acoustic field) other than the capillary force, the total volume that has disappeared (Vwater) is calculated by the speed of the particles. Since we confirm that there is no leakage at the bonding or hydrogel interfaces, the entire Vwater should be transported into the hydrogel. Because the hydrogel has auto-fluorescence and it disappears when in contact with water molecules, the wetted volume of hydrogel (Vwet) is able to be measured by the progression of an imbibition front. Therefore, the fraction of these two volumes (Vwater/Vwet) is the same as the porosity of the hydrogel and it is calculated to be around 0.1, which would be used as one of the important parameters in the theoretical calculations.
Furthermore, the dynamic measurements of the imbibition length reveal that it grows as the square root of time as shown in Fig. 1b. This result is consistent with the classical Washburn’s equation that describes the temporal evolution of the imbibition length in a porous medium34. Thus, the imbibition rate shown in the inset of Fig. 1b indicates that the capability of absorbing solution decreases as a function of time. This fact significantly affects CICP phenomenon as a restoration phase and would be discussed later on as well.
The formation of the ion-depletion zone by CICP phenomenon
The conventional ICP phenomenon generated under an external d.c. bias results in significant concentration gradients at the anodic and cathodic sides of a permselective nanoporous membrane. Typical behaviour is that the ion concentration becomes extremely low at an anodic side (referred to as an ion-depletion zone) and enriches at a cathodic side (referred to as an ion-enrichment zone) in the case of a cation-selective membrane. Although the mechanism involves complex interactions of electric fields and electrokinetic flows35,36,37, the eye-catching fact in terms of engineering applications is that ICP can play as an electrical filter. This virtual barrier rejects the entrance of charged species into the ion-depletion zone, so that it can be utilized as a water desalination/purification mechanism2,38 and a biomolecular preconcentration mechanism24,25. Instead of an electrical bias, one can achieve the same ICP phenomena if another driving force pushes or pulls an electrolyte solution into the permselective membrane. Therefore, the capillarity is used in this work to initiate ICP phenomena. In this sense, the ion-depletion zone should also be formed in CICP system, which we confirm using a confocal microscope as shown in Fig. 1c. Since the driving force for a liquid flow due to the capillarity is much weaker than that due to an electrical bias in general, the depletion zone initially stays on the ionic hydrogel for a few hours (while the zone quickly expands in a conventional ICP system within a few seconds). A noticeable observation is that the fluorescent trackers are trapped at the top corners of a microchannel as shown in Fig. 1c. If there were no such ion-depletion zone, the trackers should have been accumulated at the bottom of the microchannel. In the meantime, the centre of the fluorescent cloud exhibits the brightest signal, since the flow converges into the hydrogel. Therefore, the confocal images with most of the fluorescent signals gathered at the corners of the microchannel provide strong evidences of CICP formation on the hydrogel pad.
CICP as the spontaneous desalting mechanism
Compared with the conventional ICP system, the formation of the ion-depletion zone in CICP system takes several hours since it relies only on the capillary force. The fluorescent microscopic images are captured as shown in Fig. 2a–c. Each snapshot is captured at different instances as indicated in the caption when the depletion zone expands in maximum. As shown in Fig. 2a, the fluorescent signal linearly decreases from the bulk solution as indicated at the line profile along A–A′. This is because the CICP is driven only by the capillarity, while the conventional ICP has a sudden drop and a flat concentration profile with fast electrokinetic vortices18,19,36. Before the measurement, the reference signals of various concentrations of fluorescent dye are obtained. See Supplementary Note 3. By comparing with the reference signal intensities, the dye molecules that infer the electrolyte concentration at the interface of hydrogel pad is expected to be removed with approximately 80%. Here, the ions are assumed to behave similar to fluorescent dyes. At least, the length of the depletion zone for ions and dye molecules are the same. This is confirmed by additional numerical simulations and experimental demonstrations. See Supplementary Note 4 for detailed explanation. Although Fig. 2a is conducted on the centre-connection device, the depletion zone asymmetrically expanded (that is, the depletion zone is rarely developed on the right side of the hydrogel). This is caused by an unavoidable residual flow from the level difference between the reservoirs. Detailed discussion would be given later on. To avoid the unwanted residual flow, the end-connection device is tested as shown in Fig. 2b,c. The dye concentration profiles have a similar linear drop and the removal efficiency is expected to be around 90%, which is similar to that of absorbed liquid by mangrove28. By comparing Fig. 2b with Fig. 2c, the duration of the depletion behaviour largely depends on the composition of the hydrogel. Increasing the volume fraction of charged group (acrylic acid (AA) in this work) provides the longer duration of CICP zone.
These behaviours, especially for a static case, are solved analytically. Under the assumption that the CICP system has an ideally cation-selective membrane with electroneutrality, a one-dimensional domain ignoring a surface effect and a quasi-equilibrium, the analytical solutions of the non-dimensional concentration () and electric potential () are given by
and
where is the dimensionless spatial coordinate whose origin coincided with the bulk reservoir and is the dimensionless ionic flux induced by the imbibition through the membrane.
Although CICP experiments are conducted without any electric power source, the electric interactions between ionic species should appear in terms of electro-migration. Pe is the Péclet number, which is the ratio of a convective transportation to a diffusive transportation (uL/D in this system). See the ‘Methods’ section for full governing equations, boundary conditions and detailed derivations. These simple analytical solutions imply the most important parameters for enhancing the desalination efficiency near the membrane. The experimental results show that the electrolyte concentration is depleted near the membrane (that is, the concentration gradient is negative, dc/dx<0) so that is greater than Pe, since is equal to (Pe−) exp[Pe ]. Figure 2d,e show the concentration profile plots at the different physicochemical parameters, while Pe and are fixed at 0.1 for each plot because the microfluidic device has Pe of ∼0.1 with D∼O(10−9) m2 s−1, L (the length of the microchannel)∼O(1) mm and u (a flow velocity due to an imbibition)∼O(10) μm min−1. The physical intuitions from the plots are such that higher and lower Pe give a severer depletion. Since is proportional to the charge imbalance (, See equations (13) and (16) in the ‘Methods’ section), the higher permselectivity would enhance the formation of the ion-depletion zone. This is also confirmed in Fig. 2b,c. By comparing the plots, the higher composition of the charged group (AA) leads to a stronger depletion behaviour which implies a longer duration. However, an excess amount of the charged group in the hydrogel results in a brittle structure so that one should tweak the composition as one’s discretion. For the latter case, a diffusion limited system (that is, lower Pe) drives the depletion phenomenon more efficiently.
In Fig. 2a, the depletion zone propagates asymmetrically. An inevitable liquid level difference between each reservoir in the centre-connection device results in a residual flow over hydrogel pad. In the case of the residual flow from right to left (with the imbibition flow still converging to the centre from both reservoirs), the velocity of combined flow would increase at the right interface of the hydrogel and decrease at the left interface of the hydrogel. On the basis of the above theoretical analysis, low Pe would give a better depletion efficiency than high Pe. Since Pe is defined as a flow rate over a diffusion rate, there is a higher chance to have a stronger depletion zone at the left interface in Fig. 2a than the one at the right interface. See Supplementary Note 5 for a detailed description.
The restoration phase by a diminished imbibition
The aforementioned depletion zone is formed and sustained for a few hours, depending on the composition of the ionic hydrogel and the geometry of the microchannel. However, the ion-depletion zone collapses and the concentration restores to the bulk concentration as shown in Fig. 3. See Supplementary Videos 4, 5, 6. This has never been observed in a conventional ICP phenomenon, since it continuously applies an external power. Note that the snapshots are captured every 1, 7 and 8 h in Fig. 3a–c, respectively. One interesting feature other than the appearance of the restoration phase is that the rate of the restoration phase is much slower than the rate of the depletion phase.
The transient analysis of CICP provides a clue for understanding the restoration phase. In CICP, the imbibition rate gradually decreases as a function of time as discussed above using the Washburn’s equation34. Without loss of generality, the imbibition rate is expressed as
where uimb is the imbibition velocity, S is the imbibition parameter and t is the time. As shown above, uimb as the source of CICP is proportional to t−0.5, implying a high initial imbibition rate, but it should decay with time. Detailed descriptions are given in the ‘Methods’ section.
Through the transient model, the dynamic changes of CICP layer are investigated. With chosen parameters in the ‘Methods’ section, the depletion phase as depicted in Fig. 3d lasts for when the non-dimensional imbibition parameter is 100 (normalized by D2τD/L2). However, beyond this point, the phase transits to the restoration phase as shown in Fig. 3e because of the decreased imbibition rate. See Supplementary Video 7. These simulation results largely correspond to the experimental results in Fig. 3a–c. Although the imbibition rate decreases gradually in both the depletion and the restoration phase, the contribution of the imbibition in the depletion phase is still larger than the dissipation of ionic species induced by a diffusion. Briefly, the depletion phase initially appears (that is, the imbibition rate >the dissipation) and changes into the restoration phase (that is, the imbibition rate <the dissipation).
In Fig. 3f, the minimum cation concentrations inside the depletion zone are plotted with varying imbibition parameters, , as a function of time. The parameter could be tuned by the absorbing capability of the ionic hydrogel. The imbibition parameter is defined by (γ dcap)/(4Dη) where γ is the surface tension, dcap is the diameter of a capillary and η is the dynamic viscosity of the fluid. The increase of γ and dcap or the decrease of η are the solutions to enhancing the imbibition rate. Although the Washburn’s equation has been known to describe the imbibition through a single capillary tube, recent reports39,40,41 suggested that the functional form of the imbibition in the case of a capillary network also follows the Washburn’s equation. Nevertheless, the transient analysis with arbitrary provides the insight for understanding the transition from the depletion phase to the restoration phase. Larger generates a higher imbibition rate, which is enough to overcome the dissipation so that the depletion phase sustains longer. Moreover, the ionic hydrogel of a higher charges has more sustainable depletion phase by comparing Fig. 3b,c. The hydrogel used in Fig. 3c has the higher composition of the charged group, AA. In this theoretical model, the non-dimensional Donnan equilibrium concentration () is related to the intrinsic charge of the membrane. The membrane possessing a higher charge generates a higher ionic flux through the membrane, so that the depletion phase can last longer as shown in Fig. 3g with fixed at 100.
Discussion
In this study, we propose a new ICP phenomenon driven by capillarity. Instead of an external electric field, the inherent capillarity of nanoporous network drives a permselective ion transportation through the nanoporous hydrogel network, leading to the ion-depletion zone near the ionic hydrogel. Although the state-of-the-art desalination method utilizing a size exclusion, a thermal or electrical energy have remarkably advanced, the high level of energy requirements has hampered the development of a small-scale desalination/purification system. The fact that the CICP mechanism is working without any electrical power source offers significant merits in terms of the power consumption and the stable operation compared with the conventional electrokinetic ICP mechanism and other state-of-the-art large-scale desalination methods. Therefore, the present spontaneous desalting mechanism can lead to a power-free desalination system that is significantly useful in remote/rural and disaster-stricken areas. To achieve such a system, maintaining a quiescent fluid near the nanoporous membrane (that is, low Péclet number), enhancing the intrinsic surface charge of the membrane (that is, high ionic flux) and sustaining an imbibition at a certain level (that is, preventing restoration phase) are strongly required.
Methods
CICP device fabrication
The main building block of the CICP system is polydimethyl-siloxane (PDMS) material as schematically shown in Fig. 4. The final device consists of two layers of PDMS blocks. The bottom block is for a microchannel filled with the ionic hydrogel and the top block is for a microchannel to inject the saline water sample. Silicon masters are prepared using SU8 photoresist with the microchannel dimension of 400 μm (width) × 50 μm (depth) × 16 mm (length) for bottom PDMS block and 100 μm (width) × 15 μm (depth) × 10 mm (length) for top PDMS block. The PDMS blocks are fabricated as follows. The PDMS material (Sylgard 184 Silicone elastomer kit, Dow Corning, USA) mixed with a curing agent at 10:1 ratio is poured on the silicon masters that have desirable microchannel patterns. Then the polymer solution is cured in an oven for 4 h at 75 °C.
Detailed fabrication process from the first schematic is as follows. First, only the top PDMS for hydrogel mould is treated with oxygen plasma (CUTE-MP, Femto Science, Korea) for adhesion of the hydrogel and attached to a blank bottom PDMS piece. Then the microchannel is filled with a hydrogel precursor material consisting of 2-hydroxyethyl methacrylate (HEMA), acrylic acid (AA), a crosslinker (ethylene glycol dimethacrylate (EGDMA)) and a photoinitiator, (2,2-dimethoxy-2 phenyl-acetophenone (DMPA)) at the weight ratio of 32.51: 6: 0.41: 1.25 (named 5:1 hydrogel) or 32.51: 18: 0.41: 1.25 (named 5:3 hydrogel)12. The ratio of HEMA to AA would have a significant effect on CICP phenomenon as discussed. All the materials are purchased from Sigma-Aldrich. Ultraviolet irradiation for 8 min at the dose of 225 mW cm−2 (WUV-L50, Daihan Scientific, Korea) polymerizes the hydrogel inside the microchannel and then the top PDMS block is detached from the bottom blank PDMS piece. The hydrogel inside the microchannel is initially brittle and stiff (dry state) but the hydrogel absorbs a solution so that the portion of hydrogel becomes flexible and tough as the CICP proceeded. See Supplementary Video 8 for demonstrating only the water-absorbed portion is flexible and tough. Finally the flipped PDMS block for hydrogel and the PDMS for saline water samples are irreversibly bonded with oxygen plasma treatment. Assembled devices placed on glass slide are shown in Fig. 4b. There are two distinct designs, a ‘centre-connection device’ and an ‘end-connection device’. The centre-connection device is designed for easiness of facilitating the experiment such as sample filling and flushing, while a residual flow is able to be completely eliminated in the end-connection device.
The microchannel inside top PDMS block is filled with KCl electrolyte solution at the concentration of 300 mM mixed with a fluorescent dye (Alexa 488, Sigma Aldrich, USA) as a tracker of the ion movements and negatively charged micro-particles (d=5 μm, Invitrogen, USA) for flow field tracking. KCl is used for a simple theoretical interpretation by setting the same diffusivity of cation and anion and the concentration of 300 mM is the average salinity of brackish water.
Experimental setup
Compared with the conventional ICP process, any external voltage source is no longer required for the present CICP device. Instead, the total experimental time takes from a few hours to a few days depending on the device (for example, a conventional ICP occurs within a few seconds) so that the images of flow and concentration changes are captured every 2 min for a few days using a commercial computer software (CellSens, Olympus, Japan) and an inverted fluorescent microscope (IX 53, Olympus, Japan). The location where the hydrogel channel and the sample channel crossed is intensively analysed. A confocal microscope (FV1200, Olympus, Japan) is used for the three-dimensional concentration profiling.
Theoretical domain descriptions with governing equations
The numerical and analytical analyses are conducted for the following domain where CICP occurs. A microchannel with length L is filled with 1:1 electrolyte solution. One side of the microchannel opens to a well-mixed reservoir with bulk concentration, c0 (at x=0), while the other side contacts an ideal cation-selective surface (at x=L). Assuming that the ideal permselectivity allows us to avoid complexity in describing the transport phenomena inside the membrane17,22,23,42,43, the domain from the reservoir to the membrane surface becomes the region of interest to be theoretically investigated. Since the streaming potential generated by the imbibition through the membrane is negligible, the effects of a surface conduction and an electro-osmotic flow on the CICP are excluded. Therefore, the domain is simplified as a one-dimensional configuration.
In general, the transport of ionic species along the x axis is governed by the Nernst–Planck equation,
where t is the time, ci is the concentration of i-th species and Ji is the ionic molar flux of the i-th species along the x axis. In accordance with transport mechanisms44, the ionic flux, Ji is given by
where Di is the diffusivity of i-th species, F is the Faraday constant, R is the gas constant, T is the absolute temperature, ψ is the electric potential and u is the flow velocity of x direction. In equation (5), each term on the right-hand side corresponds to the transport mechanism due to a diffusion, an electro-migration and a convection, respectively. Although CICP experiments are conducted without any electric power source, the term of the electro-migration should be considered to describe the electric interactions between ionic species. Thus, the interactions are explained by the Poisson equation,
where ɛf is the electrical permittivity of the electrolyte solution, and c+ and c− are the concentrations of cation and anion, respectively. The convective component of the ionic flux is usually obtained by solving the Stokes equations and the continuity equation with coupled partial differential equations (4)–(6),. However, in this work, we adopt the simplified approach suggested by Dhopeshwarkar et al.45 to avoid mathematical complication. In their work, the driving force to generate ICP phenomena was the difference of applied electric potentials through the microchannel–nanochannel–microchannel system. However, our system is a microchannel–nanochannel connected system so that the theoretical domain is different and the imbibition through the membrane drives the ICP phenomena. Thus, the domain is reduced to a microchannel part rather than a whole microchannel–nanochannel connected system, while the domain of Dhopeshwarkar et al. analysed a whole microchannel–nanochannel–microchannel system. Since the effective convective transport is induced by the imbibition into the membrane, the flow velocity in the microchannel (u) is related to the imbibition rate (uimb) as
where ϕp is the porosity of the membrane, which is assumed to consist of a bundle of cylindrical capillaries. Equation (7) is derived from the fact that the flow rate in the whole system should be uniform. This imbibition is the driving force for CICP, while the driving force of conventional ICP platforms is usually the difference of applied electric potentials.
Detailed boundary conditions at the reservoir bulk (x=0) and the cation-selective surface (x=L) to describe the CICP would be discussed in the following sections.
The derivation of analytical solutions at the steady state
Since each ionic flux of the i-th species is constant over the domain in the steady state, the Nernst–Planck equation is reduced to simplified balance equations as
and
for 1:1 electrolyte solution. For simplicity, thin double layers and electroneutrality (c+=c−=c) are assumed in the microchannel as the domain. In addition, the diffusivity of each ionic species is fixed to the same value, (D+=D−=D). For the ideal cation-selective surface, no anion flux exists over the microchannel. As a result, the above equations are rewritten as
and
Combining equations (10) and (11) yields
while subtracting equation (11) from equation (10) yields
The equations (12) and (13) are further simplified by non-dimensionalization with
where Pe is known as the Péclet number, the ratio of a convective transport rate and a diffusive transport rate. Consequently, the non-dimensional forms of equations (12) and (13) are,
and
respectively. Because and at , equations (15) and (16) can be integrated under constant which is induced by the imbibition through the membrane. Therefore, the following analytical solutions are obtained with respect to the concentration and the electric potential inside CICP layer as
and
The formulation of transient model
The imbibition into a capillary tube is described by Washburn’s equation34 that is
where Limb is imbibition length, γ is the surface tension of water/capillary wall/air interface and dcap is the diameter of capillary. By differentiating equation (19), the imbibition rate is obtained as
where S is defined by γ dcap/4η. Because uimb as the source of CICP is proportional to t−0.5, the thickness of CICP layer induced by imbibition into the membrane would be decreased as time passed.
The spatiotemporal concentration profile inside the CICP layer is obtained by the following non-dimensional forms of equations (5) and (6) with appropriate boundary conditions,
and
where the characteristic scales of ci, ψ and x are denoted in equation (14), the characteristic time is chosen as the diffusion time scale, τD defined as L2/D, and the characteristic velocity is set to be D/L. In equation (21), is the non-dimensional Debye length defined as
for z:z electrolyte solution. At the reservoir bulk (), the following boundary conditions should be satisfied,
At the ideally cation-selective surface (), the cation flux through the ionic hydrogel pad is generated only by the imbibition,
and the permselectivity is described by
and
Equation (25) implies that the gradient of electrochemical potential () would be zero, so that only a convective transportation into the membrane contributes to form the CICP layer. The Donnan equilibrium concentration17,22 inside the nanoporous membrane is described by equations (26) and (27), which impedes that anion flux into the membrane should be zero due to the ideal cation selectivity. Equations (21) and (22) with boundary conditions of equations (24)–(27), are solved numerically.
Numerical method
Numerical simulations are conducted by COMSOL Mutiphysics 4.4, a commercial finite element method tool, under the one-dimensional domain and the transient analysis mode. The domain is discretized into 1,000 elements, with finer meshes in the proximity to the membrane surface. To carry out a qualitative analysis for CICP system, we assume that the diffusion time scale (τD) and the porosity (ϕp) are 500 s and 0.1, respectively. In addition, the non-dimensional Debye length and Donnan equilibrium concentration are set to be 0.01 and 5. From Washburn’s equation34, the imbibition rate is infinite at t=0 s so that the initial conditions cannot be defined. To avoid this problem, a transient analysis is conducted from t=1 s rather than from t=0 s. The time step and the termination time are set to be 1 and 3,600 s, respectively.
Additional information
How to cite this article: Park, S. et al. Capillarity ion concentration polarization as spontaneous desalting mechanism. Nat. Commun. 7:11223 doi: 10.1038/ncomms11223 (2016).
References
van den Berg, A., Craighead, H. G. & Yang, P. From microfluidic applications to nanofluidic phenomena. Chem. Soc. Rev. 39, 899–900 (2010).
Kim, S. J., Ko, S. H., Kang, K. H. & Han, J. Direct seawater desalination by ion concentration polarization. Nat. Nanotechnol. 5, 297–301 (2010).
Kwak, R., Guan, G. F., Peng, W. K. & Han, J. Y. Microscale electrodialysis: concentration profiling and vortex visualization. Desalination 308, 138–146 (2013).
Nikonenko, V. V. et al. Desalination at overlimiting currents: state-of-the-art and perspectives. Desalination 342, 85–106 (2014).
Shannon, M. A. et al. Science and technology for water purification in the coming decades. Nature 452, 301–310 (2008).
Kwak, R., Pham, V. S., Lim, K. M. & Han, J. Y. Shear flow of an electrically charged fluid by ion concentration polarization: scaling laws for electroconvective vortices. Phys. Rev. Lett. 110, 114501 (2013).
Deng, D. et al. Overlimiting current and shock electrodialysis in porous media. Langmuir 29, 16167–16177 (2013).
Deng, D. et al. Water purification by shock electrodialysis: deionization, filtration, separation, and disinfection. Desalination 357, 77–83 (2015).
Pu, Q., Yun, J., Temkin, H. & Liu, S. Ion-enrichment and ion-depletion effect of nanochannel structures. Nano. Lett. 4, 1099–1103 (2004).
Kim, S. J. & Han, J. Self-sealed vertical polymeric nanoporous junctions for high-throughput nanofluidic applications. Anal. Chem. 80, 3507–3511 (2008).
Heo, J. et al. Ultra-high-aspect-orthogonal and tunable three dimensional polymeric nanochannel stack array for BioMEMS applications. Nanoscale 6, 9681–9688 (2014).
Kim, P., Kim, S. J., Suh, K.-Y. & Han, J. Stabilization of ion concentration polarization using a heterogeneous nanoporous junction. Nano Lett. 10, 16–23 (2010).
Fosdick, S. E., Knust, K. N., Scida, K. & Crooks, R. M. Bipolar electrochemistry. Angew. Chem. Int. Ed. 52, 10438–10456 (2013).
Knust, K. N., Hlushkou, D., Anand, R. K., Tallarek, U. & Crooks, R. M. Electrochemically mediated seawater desalination. Angew. Chem. Int. Ed. 52, 8107–8110 (2013).
Leinweber, F. C. & Tallarek, U. Nonequilibrium electrokinetic effects in beds of ion-permselective particles. Langmuir. 20, 11637–11648 (2004).
Holtzel, A. & Tallarek, U. Ionic conductance of nanopores in microscale analysis systems: where microfluidics meets nanofluidics. J. Sep. Sci. 30, 1398–1419 (2007).
Dydek, E. V. et al. Overlimiting current in a microchannel. Phys. Rev. Lett. 107, 118301 (2011).
Nam, S. et al. Experimental verification of overlimiting current by surface conduction and electro-osmotic flow in microchannels. Phys. Rev. Lett. 114, 114501 (2015).
Kim, S. J., Wang, Y.-C., Lee, J. H., Jang, H. & Han, J. Concentration polarization and nonlinear electrokinetic flow near nanofluidic channel. Phys. Rev. Lett. 99, 044501 (2007).
Rubinstein, S. M. et al. Direct observation of a nonequilibrium electro-osmotic instability. Phys. Rev. Lett. 101, 236101 (2008).
Pundik, T., Rubinstein, I. & Zaltzman, B. Bulk electroconvection in electrolyte. Phys. Rev. E 72, (2005).
Green, Y. & Yossifon, G. Dynamical trapping of colloids at the stagnation points of electro-osmotic vortices of the second kind. Phys. Rev. E 87, (2013).
Rubinstein, I. & Zaltzman, B. Electro-osmotically induced convection at a permselective membrane. Phys. Rev. E 62, 2238–2251 (2000).
Kim, S. J., Song, Y.-A. & Han, J. Nanofluidic concentration devices for biomolecules utilizing ion concentration polarization: theory, fabrication, and application. Chem. Soc. Rev. 39, 912–922 (2010).
Ko, S. H. et al. Massively-parallel concentration device for multiplexed immunoassays. Lab Chip 11, 1351–1358 (2011).
Cho, I., Sung, G. & Kim, S. J. Overlimiting current through ion concentration polarization layer: hydrodynamic convection effects. Nanoscale 6, 4620–4626 (2014).
Kim, S. J., Ko, S. H., Kang, K. H. & Han, J. Direct seawater desalination by ion concentration polarization. Nat. Nanotechnol. 8, 609–609 (2013).
Parida, A. K. & Jha, B. Salt tolerance mechanisms in mangroves: a review. Trees 24, 199–217 (2010).
Beebe, D. J. et al. Functional hydrogel structures for autonomous flow control inside microfluidic channels. Nature 404, 588–590 (2000).
Lee, K. Y. & Mooney, D. J. Hydrogels for tissue engineering. Chem. Rev. 101, 1869–1880 (2001).
Shengqiang, C. & Zhigang, S. Equations of state for ideal elastomeric gels. Europhys. Lett. 97, 34009 (2012).
Marcombe, R. et al. A theory of constrained swelling of a pH-sensitive hydrogel. Soft Matter 6, 784–793 (2010).
Suzuki, A. & Hara, T. Kinetics of one-dimensional swelling and shrinking of polymer gels under mechanical constraint. J. Chem. Phys. 114, 5012–5015 (2001).
Dullien, F. A. L. & Brenner, H. Porous Media: Fluid Transport and Pore Structure Elsevier Science (2012).
Kim, S. J., Li, L. & Han, J. Amplified electrokinetic response by concentration polarization near nanofluidic channel. Langmuir 25, 7759–7765 (2009).
Kim, S. J. et al. Multi-vortical flow inducing electrokinetic instability in ion concentration polarization layer. Nanoscale 4, 7406–7410 (2012).
Lee, S.-H. et al. Sub-10 nm transparent all-around-gated ambipolar ionic field effect transistor. Nanoscale 7, 936–946 (2015).
Kwak, R., Kim, S. J. & Han, J. Continuous-flow biomolecule and cell concentrator by ion concentration polarization. Anal. Chem. 83, 7348–7355 (2011).
Elizalde, E., Urteaga, R. & Koropecki, R. R. Berli CLA. Inverse Problem of Capillary Filling. Phys. Rev. Lett. 112, 134502 (2014).
Huang, Y.-R. et al. Polymer nanocomposite films with extremely high nanoparticle loadings via capillary rise infiltration (CaRI). Nanoscale 7, 798–805 (2015).
Gruener, S., Hofmann, T., Wallacher, D., Kityk, A. V. & Huber, P. Capillary rise of water in hydrophilic nanopores. Phys. Rev. E 79, 067301 (2009).
Rubinstein, I. & Shtilman, L. Voltage against current curves of cation exchange membranes. J. Chem. Soc. Faraday Trans. 75, 231–246 (1979).
Druzgalski, C. L., Andersen, M. B. & Mani, A. Direct numerical simulation of electroconvective instability and hydrodynamic chaos near an ion-selective surface. Phys. Fluids 25, 110804 (2013).
Masliyah, J. H. & Bhattacharjee, S. Electrokinetic and Colloid Transport Phenomena Wiley (2006).
Dhopeshwarkar, R., Crooks, R. M., Hlushkou, D. & Tallarek, U. Transient effects on microchannel electrokinetic filtering with an ion-permselective membrane. Anal. Chem. 80, 1039–1048 (2008).
Acknowledgements
This work is mainly supported by Samsung Research Funding Center of Samsung Electronics under Project Number SRFC-MA1301-02. We also acknowledge the partial financial support to S.P. by the Center for Integrated Smart Sensor (Grant No. CISS-2011-0031870), H.L. by NRF (Grant Nos. 2012-0009563 and 2014-048162), I.C. by the Basic Science Research Program (Grant No. 2013R1A1A1008125) and Y.C., Y.J. by Korean Health Technology RND Project, Grants No. HI13C1468 and No. HI14C0559, respectively. All the authors gratefully acknowledge the financial support from BK21 plus program of Creative Research Engineer Development IT, Seoul National University.
Author information
Authors and Affiliations
Contributions
S.P. performed most of the desalting experiments; Y.J. analysed the steady-state solution; S.Y.S. helped in designing the experiments; I.C. imaged the ion-depletion zone using the confocal microscope; Y.C. performed the fluorescent tracking; H.L. designed and analysed the transient model; H.L., H.-Y.K. and S.J.K. supervised the project, designed the experiments and wrote the paper.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Figures 1-7, Supplementary Tables 1-2, Supplementary Notes 1-5 and Supplementary References (PDF 787 kb)
Supplementary Movie 1
Particle converging in the center connection device. (AVI 794 kb)
Supplementary Movie 2
Particle converging in the end connection device. (AVI 2489 kb)
Supplementary Movie 3
CICP operation with DI and 5:1 hydrogel in the end connection device. (AVI 497 kb)
Supplementary Movie 4
CICP operation with 300 mM sample and 5:1 hydrogel in the center connection device. (AVI 1762 kb)
Supplementary Movie 5
CICP operation with 300 mM sample and 5:1 hydrogel in the end connection device. (AVI 1450 kb)
Supplementary Movie 6
CICP operation with 300 mM sample and 5:3 hydrogel in the end connection device. (AVI 3372 kb)
Supplementary Movie 7
Transient analysis of CICP operation. (AVI 755 kb)
Supplementary Movie 8
The mechanical stiffness of swollen hydrogel. (AVI 7904 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Park, S., Jung, Y., Son, S. et al. Capillarity ion concentration polarization as spontaneous desalting mechanism. Nat Commun 7, 11223 (2016). https://doi.org/10.1038/ncomms11223
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms11223
This article is cited by
-
Eco friendly nanofluidic platforms using biodegradable nanoporous materials
Scientific Reports (2021)
-
Dynamic analysis of the extended space charge layer using chronopotentiometric measurements
Micro and Nano Systems Letters (2020)
-
Effect of evaporation through nanoporous medium on diffusiophoresis
Micro and Nano Systems Letters (2020)
-
Spontaneous diffusiophoretic separation in paper-based microfluidic device
Micro and Nano Systems Letters (2020)
-
Direct Visualization of Perm-Selective Ion Transportation
Scientific Reports (2020)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.