Design and Numerical Study of Micropump Based on Induced Electroosmotic Flow

Induced charge electroosmotic flow is a new electric driving mode. Based on the Navier–Stokes equations and the Poisson– Nernst–Planck (PNP) ion transport equations, the finite volume method is adopted to calculate the equations and boundary conditions of the induced charge electroosmotic flow. In this paper, the formula of the induced zeta potential of the polarized solid surface is proposed, and a UDF program suitable for the simulation of the induced charge electroosmotic is prepared according to this theory. At the same time, on the basis of this theory, a cross micropump driven by induced charge electroosmotic flow is designed, and the voltage, electric potential, charge density, and streamline of the induced electroosmotic micropump are obtained. Studies have shown that when the cross-shaped micropump is energized, in the center of the induction electrode near the formation of a dense electric double layer, there exist four symmetrical vortices at the four corners, and they push the solution towards both outlets; it can be found that the average velocity of the solution in the cross-flowmicrofluidic pump is nonlinear with the applied electric field, which maybe helpful for the practical application of induced electroosmotic flow in the field of micropump.


Introduction
e 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.
e 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.
e 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.
e 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][2][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 [4].With the reduction of the characteristic scale of the study of fluid motion, a new flow effect was found.e 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 [7].e phenomenon of induced electroosmotic flow mainly depends on the interaction of polarizable solids with an applied electric field to generate an electromotive phenomenon.e induced potential on the polarizable surface is critical to the induced charge electroosmotic flow.e magnitude of its zeta potential is related to the applied electric field.
e 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 [8] 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.
e induced charge electroosmotic flow (ICEOF) has been studied and applied to the microfluidic systems extensively in the last two decades.e 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][10][11] based on it, which can be applied to the driving of microfluidic chip.e model can successfully predict new phenomena when the applied voltage is too small to disrupt the salt concentration.

Theoretical Analysis
As the flow is considered steady and incompressible, the governing equations are shown below: where λ 2 0 � ε f kT/(2z 2 e 2 a 2 p n ∞ ), Pe � Ua p /D, Re � ua p /v � PeSc.e other variables are the characteristic speed u, the characteristic length a p , the kinematic viscosity v, the dielectric constant ε, the valence of ions z, the absolute temperature T, the ion concentration n ∞ , 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.e 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 where ε is the dielectric constant, ε 0 is the dielectric constant of vacuum, r is the radius, μ is the dynamic viscosity, and E is the applied electric field strength.

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.

Journal of Nanotechnology
Under the two-dimensional uniform electric eld, the analytical formula of zeta potential at ideal polarizable cylindrical surface is shown below: ).In this simulation, the ow eld, the applied electric eld, and the zeta potential control equation of the wall surface of the polarizable obstacle are shown in (4).e water used in the solution medium is related to the physical parameter:

Results and Discussions.
First of all, the electric eld of the cross channel is analyzed.In the simulation, an additional electric eld is added to the two inlets to generate an electric eld 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 eld, 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.
e 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 eld will be generated between the positive and negative electrodes.erefore, when the center of the polarizable solid surface produces an electric double layer under the action  of an applied electric eld, the ions in the solution are attracted by the electric double layer, and nally the liquid is driven to form an induced electroosmotic ow. Figure 6 shows the micropump ow diagram of induced electroosmotic ow in the cross channel under di erent electric eld intensities obtained from simulation.What can be seen from the diagram is that some of the uids will ow 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.is is mainly due to the fact that the ion concentration in the di usion layer in the electrical double layer is smaller in distance from the polarizable solid and less in drag force on the uid driven by the external electric eld, so that the uid can not ow out from the outlet but do swirling movement in volatile solids around.
As the applied electric eld 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 φ 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 φ 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.e above results show that the greater the voltage, the more easily the uid ows into the outlet channel and also can result in a more e cient driving e ect.
In general, the performance of a micropump is mainly measured by its micro uidic driving ability, which can be compared with the uid velocity at the exit.
is paper mainly simulates cross induced electroosmotic micropumps with a two-dimensional structure.erefore, 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.erefore, the speed of one of the outlets will be studied separately in this paper.Figure 7 shows the velocity pro le 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 gure, the velocity at the outlet is parabolic, and the larger the voltage is, the greater the velocity is, and the driving e ect of the micropump is better.When the voltage at the inlet is 100 V, the maximum uid velocity at the outlet of the electrical microchannel reaches 10 mm/s, and as the voltage increases, the drive speed increases faster and faster.is shows that the use of the micropump can produce a good driving e ect.
Figure 8 shows the relation between the average velocity of single outlet and applied electric eld 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 gure 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 t the numerical simulation results to get the cross structure of the micropump single-exit average velocity and applied electric eld curve: y 0.001x 2 − 0.081x + 2.0618.

Conclusions
In summary, the mechanism of induced electroosmotic ow is studied in depth.e analytical solution of induced zeta potential at polarizable solid surface is proposed by analyzing the governing equations of induced electroosmotic ow.Based on this theory, a UDF program suitable for induced charge electroosmotic ow simulation is developed.At the same time, the cross micropump driven by induced charge electroosmotic ow was designed, and the voltage, potential, charge density, and ow eld of the induced micropump were obtained.e results show that the cross induced charge electroosmosis micropump has a nonlinear relationship with the applied electric eld, which is more powerful than that of the traditional electroosmotic pump.

Conflicts of Interest
e authors declare that they have no con icts of interest.

2. 1 .
Governing Equations and Boundary Conditions.In this study, the theoretical model is based on the Navier-Stokes equation[12] of viscous fluid flow, combined with the Poisson-Nernst-Planck (PNP) ion transport equation.

3. 1 .( 2 ) 3 )
Model and Boundary Conditions.As shown in Figure2of the cross-shaped induced electroosmotic micropump, a cylindrical polarizable solid is embedded in the middle of the cross-shaped channel.e distance between the energized electrodes is L 2000 μm, the width of the microchannel is W 200 μm, and the diameter of a circular polarizable solid is ϕ 100 μm.Using the Gambit software to mesh the 2D micropump model and pass the grid independency veri cation, the cross geometry model of micropump is shown in Figure3, and the total number of grids nally con rmed is 20,000.e model's boundary conditions are set as follows: (1) Boundary conditions of the surface potential Inlet and outlet:φ inlet−1 φ a , φ inlet−2 0 V; φ outlet−1 φ outlet−2 0 V;e surface of polarizes solid: zφ/zn Boundary conditions of ion concentration Inlet and outlet: c 2; Side wall: c 2; e surface of polarized solid: zc/zn Boundary conditions of ion density Inlet and outlet: q 0; Side wall: q 0; e surface of polarized solid: zq/zn ⇀ −c(zφ/zn ⇀

Figure 4 :Figure 5 :
Figure 4: Voltage diagram in the cross channel.