Design and Numerical Study of Micropump Based on Induced Electroosmotic Flow.
The huge technological advances have made people's demand for technology products continue to move toward the direction of portability, miniaturization, and intelligence. With the improvement of living standards, more and more attention is being paid to health problems and the accuracy of the results, and the monitoring methods have also been set at higher standards. At this moment, the microfluidic chip is a good choice for further development and popularization of real-time diagnostic technology.
The so-called microfluidics chip refers to a chemical or biological lab built on a chip that is only a few square centimeters. It integrates the processes of biological and chemical reactions and separation and detection into microchannels, while using the design of microchannel networks for microfluidic control and transport and ultimately enables various functions in chemical or biological laboratories, that is, lab on a chip. Much research has been made under the leadership of a large number of scholars such as Lin Bingcheng, Qin Jianhua, and so on, all of which contributed to the development of microfluidic chip technology in China and even in the world. The internationally renowned magazine Lab on a Chip even published a special album titled "Focus on China" on the 10th anniversary of its founding, which is to affirm the important contribution made by Chinese scholars to the research of microfluidic technology.
The delivery mixing, reaction, separation, and control of microfluidics are key components of microfluidics system. Because of its small characteristic scale, most of the fluid flowing in the microchannels is laminar. Particles, droplets, or bubbles generally within the microchannels belong to the field of low Reynolds number flow theory [1-3]. Due to the sharp decrease of the volume to surface area ratio, the study found that the laws and phenomena of fluid movement at the microscale are different from the macro environment; the continuity equation in the three equations of hydrodynamics may no longer be suitable for use in microfluidics . With the reduction of the characteristic scale of the study of fluid motion, a new flow effect was found. The coupling of the electric field force, the flow field, the temperature field and the ion motion within the microchannel, the electroosmotic flow, electrophoresis, induced electroosmotic flow, and other electrical phenomena can be used to achieve microfluidic (microparticle) transport and control [5, 6]. Electrokinetic phenomena in microfluidics are caused by the interaction between the applied electric field and the diffusion layer in the electric double layer. More electrokinetic phenomena studied at present include electrophoresis, dielectric electrophoresis, electroosmotic flow, and induced electroosmotic flow. Electrokinetic phenomena can be divided into linear (electrophoresis and electroosmotic flow) and nonlinear (dielectric electrophoresis, electrophoresis, and induced electroosmotic flow) electrokinetic phenomena according to whether the zeta potential in the electrokinetic phenomenon changes with an applied electric field. In this paper, the phenomenon of induced electroosmotic flow is used.
Induced electroosmosis (ICEO) is a phenomenon driven by electrostatic forces under applied electric field and is a variant of the electroosmotic phenomenon . The phenomenon of induced electroosmotic flow mainly depends on the interaction of polarizable solids with an applied electric field to generate an electromotive phenomenon. The induced potential on the polarizable surface is critical to the induced charge electroosmotic flow. The magnitude of its zeta potential is related to the applied electric field. The earliest induced electroosmotic flow was discovered by Romans et al. at the end of the 20th century. Subsequently, in 2004, Bazant and Squires perfected the relevant theory and formally proposed the concept of inducing electroosmotic flow. And the study of the mixing  and transporting of the fluid in the simple microchannel is accomplished by using this theory. By 2005, Levitan used experimental methods to confirm the correctness of the basic model of induced electroosmotic flow.
The induced charge electroosmotic flow (ICEOF) has been studied and applied to the microfluidic systems extensively in the last two decades. The phenomenon is used by Wu and Li to realize the function of fluid mixing and flow regulation in microfluidic chips; Zhao and Bau used induced electroosmotic flow to enhance chaotic flow to improve the mixing efficiency of microfluidics; Yariv, Bau, and Li et al. gave attention and conducted preliminary studies on inducing particle-wall effect in electroosmotic flow; Peng then experimentally found that the higher the zeta potential of the electrical double layer around the surface of the polarizable solid, the more particles agglomerated; demonstrating the feasibility of using micronanoparticle manipulation to induce electroosmosis. Harbin Institute of Technology, Peng and Jia innovated the use of ITO conductive glass as the electrode, based on the principle of induced electroosmotic flow and implementation of micronanoparticle manipulation.
Compared with the classical electroosmotic flow, the induced electroosmotic flow can obtain a higher driving speed under the same voltage, so we design a micropump [9-11] based on it, which can be applied to the driving of microfluidic chip.
2. Theoretical Analysis
2.1. Governing Equations and Boundary Conditions. In this study, the theoretical model is based on the NavierStokes equation  of viscous fluid flow, combined with the Poisson-Nernst-Planck (PNP) ion transport equation. The model can successfully predict new phenomena when the applied voltage is too small to disrupt the salt concentration.
As the flow is considered steady and incompressible, the governing equations are shown below:
[Pe.sup.-1][nabla] x ([- or +][n.sub.+] [nabla][phi] - [nabla][n.sub.[+ or -]]) + [??] x [nabla][n.sub.[+ or -]] = 0, (1)
[nabla] x (Pe[??]q) = [nabla] x (cV[phi] + [nabla]q), (2)
[nabla] x (Pe[??]c) = [nabla] x (q[nabla][phi] + Vq), (3)
2[[lambda].sup.2.sub.0][[nabla].sup.2][phi] = -q, (4)
[nabla] x [??] = 0, (5)
[mathematical expression not reproducible], (6)
where [[lambda].sup.2.sub.0] = [[epsilon].sub.f]kT/(2[z.sup.2][e.sup.2][a.sup.2.sub.p][n.sub.[infinity]]]), Pe = U[a.sub.p]/D, Re = u[a.sub.p]/v = PeSc. The other variables are the characteristic speed u, the characteristic length [a.sub.p], the kinematic viscosity v, the dielectric constant [epsilon], the valence of ions z, the absolute temperature T, the ion concentration [n.sub.[infinity]], the diffusion coefficient D, and Boltzmann's constant k.
Equations (2)-(4) are solved to obtain ion concentration and density distribution, and then (5) and (6) are solved to get the information of flow field. The zeta potential in the electroosmotic flow is induced by an applied electric field, and the magnitude of the potential depends on the applied electric field. According to the relevant theoretical study, it is found that the induced tangent slip velocity of the electric double layer on the polarizable solid surface in the electroosmotic flow is
[u.sub.ICEO] = [[[epsilon][[epsilon].sub.0]r/[mu]] [E.sup.2], (7)
where [epsilon] is the dielectric constant, [[epsilon].sub.0] is the dielectric constant of vacuum, r is the radius, [mu] is the dynamic viscosity, and E is the applied electric field strength.
2.2. Zeta Potential Verification. Induced charge electroosmosis flow (ICEOF) phenomenon, which is caused by the interaction between the applied electric field and the electric double layer formed on the polarizable surface, and zeta potential changes on the polarizable solid are shown in Figure 1.
Under the two-dimensional uniform electric field, the analytical formula of zeta potential at ideal polarizable cylindrical surface is shown below:
[zeta] = 2[E.sub.0]a cos[beta]. (8)
3. Numerical Simulation
3.1. Model and Boundary Conditions. As shown in Figure 2 of the cross-shaped induced electroosmotic micropump, a cylindrical polarizable solid is embedded in the middle of the cross-shaped channel. The distance between the energized electrodes is L = 2000 [micro]m, the width of the microchannel is W = 200 [micro]m, and the diameter of a circular polarizable solid is [phi] = 100 [micro]m. Using the Gambit software to mesh the 2D micropump model and pass the grid independency verification, the cross geometry model of micropump is shown in Figure 3, and the total number of grids finally confirmed is 20,000.
The model's boundary conditions are set as follows:
(1) Boundary conditions of the surface potential Inlet and outlet: [[phi].sub.inlet-1] = [[phi].sub.a], [[phi].sub.inlet-2] = 0 V; [[phi].sub.outlet-1] = [[phi].sub.outlet-2] = 0 V;
The surface of polarizes solid: [partial derivative][phi]/[partial derivative][??] = -([partial derivative]c/[partial derivative][??])/q.
(2) Boundary conditions of ion concentration
Inlet and outlet: c = 2;
Side wall: c = 2;
The surface of polarized solid: [partial derivative]c/[partial derivative][??] = -q([partial derivative][phi]/[partial derivative][??]).
(3) Boundary conditions of ion density
Inlet and outlet: q = 0;
Side wall: q = 0;
The surface of polarized solid: [partial derivative]q/[partial derivative][??] = -c([partial derivative][phi]/[partial derivative][??]).
In this simulation, the flow field, the applied electric field, and the zeta potential control equation of the wall surface of the polarizable obstacle are shown in (4). The water used in the solution medium is related to the physical parameter: [[epsilon].sub.r] = 80, [[epsilon].sub.0 = 8.85e - 12 C x [V.sup.-1][m.sup.-1], [mu] = 1.003[e.sup.-3 ]kg x [m.sup.-1][s.sup.-1], [rho] = 998.2 kg/[m.sup.3].
3.2. Results and Discussions. First of all, the electric field of the cross channel is analyzed. In the simulation, an additional electric field is added to the two inlets to generate an electric field from the positive electrode to the negative electrode in the solution medium in the microchannel, as shown in Figure 4. At the same time, under the action of an applied electric field, the centrally located polarizable electrode is polarized, and the opposite ion in the adsorption solution forms a close-packed charge layer on the surface, eventually producing an electric double layer near the surface. The potential is the zeta potential, and the charge density around the polarizable solid is shown in Figure 5. In the program, the negative terminal defaults to zero, so the potential and charge in the positive direction will be more dense, but after the power is applied, an electric field will be generated between the positive and negative electrodes.
Therefore, when the center of the polarizable solid surface produces an electric double layer under the action of an applied electric field, the ions in the solution are attracted by the electric double layer, and finally the liquid is driven to form an induced electroosmotic flow. Figure 6 shows the micropump flow diagram of induced electroosmotic flow in the cross channel under different electric field intensities obtained from simulation. What can be seen from the diagram is that some of the fluids will flow along the polarizable solid surface from the left and right inlet to the outlet. Fluid at a distance farther away from the polarizable solids does not enter the exit channel but instead creates vortices around the polarizable solids. This is mainly due to the fact that the ion concentration in the diffusion layer in the electrical double layer is smaller in distance from the polarizable solid and less in drag force on the fluid driven by the external electric field, so that the fluid can not flow out from the outlet but do swirling movement in volatile solids around.
As the applied electric field increases, the shape of the vortex around the polarizable solid can also be found from Figure 6 above. When the voltage at the inlet is [[phi].sub.a] = 10 V, the four vortices are basically at four corners and distributed evenly. With the increase of voltage, the four vortices around the polarizable solid gradually move toward the exit channel. When the voltage at the inlet is [[phi].sub.a] = 300 V, it can be clearly seen that the four vortices basically entered the interior of the exit channel. At the same time, the distance between the two vortices of polarizable solids increases with increasing voltage. The above results show that the greater the voltage, the more easily the fluid flows into the outlet channel and also can result in a more efficient driving effect.
In general, the performance of a micropump is mainly measured by its microfluidic driving ability, which can be compared with the fluid velocity at the exit. This paper mainly simulates cross induced electroosmotic micropumps with a two-dimensional structure. Therefore, it is necessary to study the speed of its outlet. According to the simulation results, under the ideal conditions, micropump at the upper and lower exit has the same speed. Therefore, the speed of one of the outlets will be studied separately in this paper. Figure 7 shows the velocity profile at one outlet, where the vertical axis is the exit speed v (mm/s) and the abscissa is the distance between the solution and the exit distance l (mm). As can be seen from the figure, the velocity at the outlet is parabolic, and the larger the voltage is, the greater the velocity is, and the driving effect of the micropump is better. When the voltage at the inlet is 100 V, the maximum fluid velocity at the outlet of the electrical microchannel reaches 10 mm/s, and as the voltage increases, the drive speed increases faster and faster. This shows that the use of the micropump can produce a good driving effect.
Figure 8 shows the relation between the average velocity of single outlet and applied electric field strength, where the ordinate is the average speed at a single exit v (mm/s) and the abscissa is the voltage at the power source U (V). It can be seen from the figure that the average speed at a single outlet is a quadratic nonlinear relationship with the power supply voltage, and when the power supply voltage is higher, the average speed of the micropump increases faster. When the voltage is greater than 100 V, the average speed of the microchannel outlet at this time is already close to 10 mm/s. At this point, we linearly fit the numerical simulation results to get the cross structure of the micropump single-exit average velocity and applied electric field curve: y = 0.001[x.sup.2] - 0.08Th + 2.0618.
In summary, the mechanism of induced electroosmotic flow is studied in depth. The analytical solution of induced zeta potential at polarizable solid surface is proposed by analyzing the governing equations of induced electroosmotic flow. Based on this theory, a UDF program suitable for induced charge electroosmotic flow simulation is developed. At the same time, the cross micropump driven by induced charge electroosmotic flow was designed, and the voltage, potential, charge density, and flow field of the induced micropump were obtained. The results show that the cross induced charge electroosmosis micropump has a nonlinear relationship with the applied electric field, which is more powerful than that of the traditional electroosmotic pump.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This study was supported by the National Natural Science Foundation of China (11472260).
 J. Z. Lin, Y. L. Wang, P. J. Zhang, and X. K. Ku, "Mixing and orientation behaviors of cylindrical particles in a mixing layer of an Oldroyd-B fluid," Chemical Engineering Science, vol. 176, pp. 270-284, 2018.
 J. Lin, X. Pan, Z. Yin, and X. Ku, "Solution of general dynamic equation for nanoparticles in turbulent flow considering fluctuating coagulation," Applied Mathematics and Mechanics, vol. 37, no. 10, pp. 1275-1288, 2016.
 M. Z. Yu, J. Z. Lin, and T. Chan, "Effect of precursor loading on non-spherical Ti[O.sub.2] nanoparticle synthesis in a diffusion flame reactor," Chemical Engineering Science, vol. 63, no. 9, pp. 2317-2329, 2008.
 M. Z. Yu, Y. Y. Liu, J. Z. Lin, and M. Seipenbusch, "Generalized TEMOM scheme for solving the population balance equation," Aerosol Science and Technology, vol. 49, no. 11, pp. 1021-1036, 2015.
 K. Ward and Z. Hugh Fan, "Mixing in microfluidic devices and enhancement methods," Journal of Micromechanics and Microengineering, vol. 25, no. 9, p. 094001, 2015.
 K. Masilamani, S. Ganguly, C. Feichtinger, D. Bartuschat, and U. Rude, "Effects of surface roughness and electrokinetic heterogeneity on electroosmotic flow in microchannel," Fluid Dynamics Research, vol. 47, no. 3, p. 035505, 2015.
 C. Canpolat, "Induced-charge electro-osmotic flow around cylinders with various orientations," Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 231, no. 21, pp. 4507-4066, 2017.
 C. Wang, Y. Song, X. Pan, and R. Dongqing Li, "A novel microfluidic valve controlled by induced charge electroosmotic flow," Journal of Micromechanics and Microengineering, vol. 26, no. 7, p. 075002, 2016.
 K. Bengtsson and N. D. Robinson, "A large-area, all-plastic, flexible electroosmotic pump," Microfluidics and Nanofluidics, vol. 21, p. 178, 2017.
 A. K. R. Lai, C. C. Chang, and C. Y. Wang, "Optimizing electroosmotic pumping rates in a rectangular channel with vertical gratings," Physics of Fluids, vol. 29, no. 8, p. 082002, 2017.
 M. Gao and L. Gui, "A handy liquid metal based electroosmotic flow pump," Lab on a Chip, vol. 14, no. 11, pp. 18661872, 2014.
 M. Z. Yu, J. Z. Lin, and T. Chan, "Numerical simulation of nanoparticle synthesis in diffusion flame reactor," Powder Technology, vol. 181, no. 1, pp. 9-20, 2008.
Kai Zhang (iD), Lengjun Jiang, Zhihan Gao, Changxiu Zhai, Weiwei Yan, and Shuxing Wu
Institute of Fluid Engineering of China Jiliang University, Hangzhou 310018, China
Correspondence should be addressed to Kai Zhang; firstname.lastname@example.org
Received 1 February 2018; Accepted 5 April 2018; Published 9 May 2018
Academic Editor: Xiaoke Ku
Caption: Figure 1: Zeta potential on the polarizable solid surface under the electric filed of [E.sub.0].
Caption: Figure 2: The geometry of micropump.
Caption: Figure 3: Mesh division of cross micropump.
Caption: Figure 4: Voltage diagram in the cross channel.
Caption: Figure 5: Charge density diagram the in cross channel.
Caption: Figure 6: Streamline diagram of cross channel under different electric field intensities. (a) [[phi].sub.a] = 10 V, (b) [[phi].sub.a]= 100 V, (c) [[phi].sub.a]= 200 V, and (d) [[phi].sub.a]= 300 V.
Caption: Figure 7: Velocity profiles at individual outlets at different voltages.
Caption: Figure 8: Relation between average velocity of single outlet and applied electric field strength.
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Zhang, Kai; Jiang, Lengjun; Gao, Zhihan; Zhai, Changxiu; Yan, Weiwei; Wu, Shuxing|
|Publication:||Journal of Nanotechnology|
|Date:||Jan 1, 2018|
|Previous Article:||Simple Modeling of the Ratio of Fields at a Tip and a Contacting Surface with External Illumination.|
|Next Article:||Green Synthesis of Iron Nanoparticles: Application on the Removal of Petroleum Oil from Contaminated Water and Soils.|