Starting Electroosmosis in a Fibrous Porous Medium with Arbitrary Electric Double-Layer Thickness

Abstract

The transient electroosmotic response in a charged porous medium consisting of a uniform array of parallel circular cylindrical fibers with arbitrary electric double layers filled with an electrolyte solution, for the stepwise application of a transverse electric field, is analyzed. The fluid momentum conservation equation is solved for each cell by using a unit cell model, where a single cylinder is surrounded by a coaxial shell of the electrolyte solution. A closed-form expression for the transient electroosmotic velocity of the bulk fluid in the Laplace transform is obtained as a function of the ratio of the cylinder radius to the Debye screening length and the porosity of the fiber matrix. The effect of the fiber matrix porosity on the continuous growth of the electroosmotic velocity over time is substantial and complicated. For a fiber matrix with larger porosity, the bulk fluid velocity takes longer to reach a certain percentage of its final value. Although the final value of the bulk fluid velocity generally increases with increasing porosity, early velocities may decrease with increasing porosity. For a given fiber matrix porosity, the transient electroosmotic velocity is a monotonically increasing function of the ratio of the cylinder radius to the Debye length.

1. Introduction

The electroosmotic flows of electrolyte solutions in charged porous media or microchannels are of great importance in many technological areas, both fundamentally and practically. Although the elementary mathematical analysis of electroosmosis was mainly performed in the steady state [1,2,3,4], its transient behavior is sometimes as important as its steady-state behavior. In various applications, pulsating or alternating electric fields are used for electroosmotic measurements and manipulations [5,6,7,8]. Particularly, the time evolution of the fluid velocity is relevant to electroosmosis in microfluidics in the sub-millisecond range, such as biochip operations in high-speed electrokinetic separation processes, the use of short-duration pulsed electric fields to distinguish between colloidal particle velocities and background electroosmotic flow or to suppress thermal-zone broadening, and micro-pumping by means of AC or modulated DC fields [9]. The time-periodic and transient starting electroosmotic flows of electrolyte solutions in microchannels with various cross-sectional shapes have been analyzed in the past for the electric double layers of small [10,11] and arbitrary [12,13,14] thicknesses.

In the real applications of electroosmotic flows, a porous medium is frequently encountered. To avoid the difficulties of complex geometries, unit cell models are often used to predict the effects of interactions between small solid entities constituting the porous medium on electroosmotic velocity [15,16,17,18,19,20,21,22,23,24]. These models allow a group of entities to be divided into numerous identical cells, each with an entity at its center. Therefore, the boundary value problem of multiple solid bodies is reduced to dealing with a single element consisting of a solid body and a concentric shell of fluid surrounding it. Different shapes of unit cells can be chosen, but a spherical or cylindrical shape is the most convenient.

The transient starting electroosmotic flow of electrolyte solutions in the axial direction of a fibrous porous medium consisting of a uniform array of parallel circular dielectric cylinders, with arbitrary electric double layer thickness and zeta potential after the application of a step-function electric field, was analyzed using a unit cell model [25]. This analysis was also extended approximately to the starting electroosmotic flow in the transverse direction of the fibrous porous medium for the case of constant zeta potential and a thin but finite double layer (say, $\kappa a\ge 10$, where ${\kappa }^{-1}$ is the Debye screening length and $a$ is the radius of cylinders) [26].

In the practical case of electroosmosis in charged porous media, the thickness of the electric double layers can be comparable to the pore size and the dependence of the dynamic response of the electroosmotic flow on $\kappa a$ over a wide range must be considered. In this paper, we present an analytical study of the transient electroosmosis induced by the sudden application of a constant transverse electric field in a matrix of parallel circular cylindrical fibers with arbitrary double-layer thickness using a unit cell model. The zeta potential (or surface charge density) of the cylinders is assumed to be uniform and small. The transient response of the electroosmotic velocity of the electrolyte solution as a function of relevant parameters is explicitly obtained.

2. Analysis

We consider a Newtonian liquid solution of arbitrary electrolytes with the viscosity $\eta$ and density $\rho$, filled in a charged fibrous porous medium consisting of a uniform array of parallel dielectric cylinders of arbitrary radius $a$ and Debye screening length ${\kappa }^{-1}$ (characteristic thickness of the electric double layer). Since the initial time $t=0$, a constant electric field $E_{\infty }$ is applied in the $x$ direction normal to the axes of the cylinders. As a result, the electrolyte solution starts flowing in the same direction with an instantaneous bulk velocity $-U(t)$ (where $U$ is positive and negative if the cylinders are positively and negatively charged, respectively) due to electroosmosis. A unit cell model is used in which each circular cylinder is surrounded by a coaxial cylindrical shell of liquid with an outer radius $b$, as shown in Figure 1; therefore, the liquid-to-cell volume ratio $1-\varphi$ is equal to the porosity of the fiber matrix, where $\varphi ={(a/b)}^2$. The origin of the polar coordinate system $(r,\theta )$ is set on the axis of the cylinder and $\theta =0$ means the $x$ axis. Clearly, this two-dimensional problem for a cell is symmetric about the $x$ axis. Our goal is to determine the transient electroosmotic velocity $U(t)$.

The electrical potential distribution $\psi (r,\theta )$ in the fluid around the cylinder in a unit cell can be expressed as the equilibrium potential distribution ${\psi }_{\mathrm{eq}}(r)$ induced by the surface charge of the cylinder and mobile ions in the surrounding electrical double layer (which satisfies the Poisson–Boltzmann equation) added with the perturbed potential distribution ${\psi }_{\mathrm{a}}(r,\theta )$ caused by the external electric field $E_{\infty }$ (which satisfies the Laplace equation),

(1)$\psi ={\psi }_{\mathrm{eq}}+{\psi }_{\mathrm{a}},$

where [22]

(2)${\psi }_{\mathrm{eq}}=\zeta {\displaystyle \frac{I_1(\kappa b)K_0(\kappa r)+I_0(\kappa r)K_1(\kappa b)}{I_1(\kappa b)K_0(\kappa a)+I_0(\kappa a)K_1(\kappa b)}}={\displaystyle \frac{\sigma }{\epsilon \kappa }}{\displaystyle \frac{I_1(\kappa b)K_0(\kappa r)+I_0(\kappa r)K_1(\kappa b)}{I_1(\kappa b)K_1(\kappa a)-I_1(\kappa a)K_1(\kappa b)}},$

(3)${\psi }_{\mathrm{a}}=-{\displaystyle \frac{E_{\infty }}{\chi }}(r+{\displaystyle \frac{a^2}{r}})\cos \theta .$

Here, $\zeta$ and $\sigma$ are the zeta potential and surface charge density, respectively, of the cylinder, and $I_n$ and $K_n$ are the modified Bessel functions of order n of the first and second kinds, respectively, $\chi =1+\varphi$ if the Dirichlet condition is used for the perturbed potential distribution on the outer boundary of the cell ($r=b$), $\chi =1-\varphi$ if the Neumann condition is used there, and $\epsilon$ is the dielectric permittivity of the fluid. Note that Equation (2) is valid when the electrical potential is small (in which the Debye–Hückel approximation applies). In comparison with the Neumann condition at $r=b$, the Dirichlet condition has an advantage that both the radial and angular components of the perturbed potential gradient are specified [16,20].

The transient creeping flow of the electrolyte solution is governed by the modified Stokes equation [27]

(4)${\nabla }^2({\nabla }^2-{\displaystyle \frac{1}{\nu }}{\displaystyle \frac{\mathsf{\partial }}{\mathsf{\partial }t}})\Psi ={\displaystyle \frac{\epsilon {\kappa }^2{\psi }_{\mathrm{a}}\mathrm{d}{\psi }_{\mathrm{eq}}}{\eta r\mathrm{d}r}}\tan \theta ,$

where $\nu$ is the kinematic viscosity ($=\eta /\rho$) of the fluid, $t$ is time, and $\Psi$ is the Lagrange stream function. The boundary and initial conditions for the fluid velocity components $v_r$ and $v_{\theta }$ in a unit cell are [24,26]

(5)$r=a: v_r=-{\displaystyle \frac{1}{r}}{\displaystyle \frac{\mathsf{\partial }\Psi }{\mathsf{\partial }\theta }}=0, v_{\theta }={\displaystyle \frac{\mathsf{\partial }\Psi }{\mathsf{\partial }r}}=0,$

(6a)$r=b: v_r=-U\cos \theta ,$

(6b)${\tau }_{r\theta }=\eta [r{\displaystyle \frac{\mathsf{\partial }}{\mathsf{\partial }r}}({\displaystyle \frac{v_{\theta }}{r}})+{\displaystyle \frac{1}{r}}{\displaystyle \frac{\mathsf{\partial }{v}_r}{\mathsf{\partial }\theta }}]=0\hspace{1em}\mathrm{for} \mathrm{the} \mathrm{Happel} \mathrm{model},$

(6c)${(\nabla \times v)}_z={\displaystyle \frac{1}{r}}{\displaystyle \frac{\mathsf{\partial }}{\mathsf{\partial }r}}(r{v}_{\theta })-{\displaystyle \frac{1}{r}}{\displaystyle \frac{\mathsf{\partial }{v}_r}{\mathsf{\partial }\theta }}=0\hspace{1em}\mathrm{for} \mathrm{the} \mathrm{Kuwabara} \mathrm{model},$

(7)$t=0: v_r=v_{\theta }=0,$

where $-U(t)$ is the transient bulk fluid velocity of electroosmosis to be determined, and $U(0)=0$. The steady-state electroosmotic velocity (as $t\to \infty$) is available in the literature [22]. For suspensions of charged spheres in electrolyte solutions, the Happel (free surface) and Kuwabara (surface with zero vorticity) models result in comparable quantities in the electrophoretic mobility. However, the Happel model has a significant advantage in that it does not require an exchange of mechanical energy between the unit cell and its surroundings [28]. Thus, the Happel model should be more satisfactory than the Kuwabara model.

In terms of the solution form $\Psi =(\epsilon \zeta E_{\infty }a/\eta )g(r,t)\sin \theta$ and the Laplace transform, Equation (4) can be expressed as

(8)$({\displaystyle \frac{{\mathrm{d}}^2}{\mathrm{d}{r}^2}}+{\displaystyle \frac{1}{r}}{\displaystyle \frac{\mathrm{d}}{\mathrm{d}r}}-{\displaystyle \frac{1}{r^2}})({\displaystyle \frac{{\mathrm{d}}^2}{\mathrm{d}{r}^2}}+{\displaystyle \frac{1}{r}}{\displaystyle \frac{\mathrm{d}}{\mathrm{d}r}}-{\displaystyle \frac{1}{r^2}}-{\displaystyle \frac{s}{\nu }})\overline{g}(r)={\displaystyle \frac{{\kappa }^3}{sa}}G(r),$

where the transformed variable is represented by an overbar, $s$ is the transform parameter, and

(9)$G(r)={\displaystyle \frac{1}{\chi }}[1+{({\displaystyle \frac{a}{r}})}^2]{\displaystyle \frac{I_1(\kappa b)K_1(\kappa r)-I_1(\kappa r)K_1(\kappa b)}{I_1(\kappa b)K_1(\kappa a)-I_1(\kappa a)K_1(\kappa b)}}.$

The general solution to Equation (4) or (8) satisfying the initial condition in Equation (7) is [27]

(10)$\overline{\Psi }={\displaystyle \frac{\epsilon }{\eta }}\zeta E_{\infty }a[C_1{\displaystyle \frac{r}{a}}+C_2{\displaystyle \frac{a}{r}}+C_3{I}_1(Ar)+C_4{K}_1(Ar)+{\overline{g}}_{\mathrm{p}}(r)]\sin \theta ,$

with the particular solution

(11)${\overline{g}}_{\mathrm{p}}(r)={\displaystyle \frac{-{\kappa }^3}{2s{A}^2}}[{\displaystyle \frac{r}{a}}{f}_1(r)-{\displaystyle \frac{a}{r}}{f}_2(r)+2K_1(Ar)f_3(r)-2I_1(Ar)f_4(r)],$

where

(12a)$f_1(r)={\displaystyle {\int }_a^{r}G(r)}\mathrm{d}r,$

(12b)$f_2(r)={\displaystyle {\int }_a^{r}G(r)}{({\displaystyle \frac{r}{a}})}^2\mathrm{d}r,$

(12c)$f_3(r)={\displaystyle {\int }_a^r{I}_1(Ar)G(r)}({\displaystyle \frac{r}{a}})\mathrm{d}r,$

(12d)$f_4(r)={\displaystyle {\int }_a^r{K}_1(Ar)G(r)}({\displaystyle \frac{r}{a}})\mathrm{d}r,$

and $A=\sqrt{s/\nu }$. The coefficients $C_1$, $C_2$, $C_3$, and $C_4$ in Equation (10) result from the boundary conditions in Equations (5) and (6) as

(13a)$C_1={\displaystyle \frac{\eta \overline{U}{H}_1}{\epsilon \zeta E_{\infty }{H}_0}}+{\displaystyle \frac{a{\kappa }^3{H}_2}{A^{2}s{H}_0}},$

(13b)$C_2={\displaystyle \frac{\eta \overline{U}{H}_3}{\epsilon \zeta E_{\infty }{H}_0}}+{\displaystyle \frac{a{\kappa }^3{H}_4}{A^{2}s{H}_0}},$

(13c)$C_3={\displaystyle \frac{\eta \overline{U}{H}_5}{\epsilon \zeta E_{\infty }{H}_0}}+{\displaystyle \frac{a{\kappa }^3{H}_6}{A^{2}s{H}_0}},$

(13d)$C_4={\displaystyle \frac{\eta \overline{U}{H}_7}{\epsilon \zeta E_{\infty }{H}_0}}+{\displaystyle \frac{a{\kappa }^3{H}_8}{A^{2}s{H}_0}},$

where the nine coefficients $H_0$, $H_1$, …, and $H_8$ are given in Appendix A.

The transformed hydrodynamic force acting on the unit cell at its virtual surface $r=b$ (in the $x$ direction) per unit length is [27]

(14)$\overline{F}(s)={{\displaystyle {\int }_0^{2\mathsf{\pi }}[\eta r^3\frac{\mathsf{\partial }}{\mathsf{\partial }r}(\frac{{\nabla }^2\overline{\Psi }}{r})-\rho r^{2}s\frac{\mathsf{\partial }\overline{\Psi }}{\mathsf{\partial }r}+\epsilon {\kappa }^{2}r{\psi }_{\mathrm{eq}}\frac{\mathsf{\partial }{\psi }_{\mathrm{a}}}{s\mathsf{\partial }\theta }]}}_{r=b}\sin \theta \mathrm{d}\theta .$

The cell is electrically neutral as a whole (the total charge of solid plus liquid). Therefore, both the electrostatic force and the net (i.e., only hydrodynamic) force acting on it must be zero. Applying this constraint to Equation (14) (with the substitution of Equations (2), (3), (10), (11) and (13)), we obtain the normalized transient electroosmotic velocity in the Laplace transform as

${\displaystyle \frac{\eta \overline{U}\left(s\right)}{\epsilon \zeta E_{\infty }}}={\displaystyle \frac{\nu {\kappa }^3}{2s^2}}\{{\displaystyle \frac{2(a^2+b^2)H_0}{{\kappa }^{2}b\chi [I_1(\kappa b)K_0(\kappa a)+I_0(\kappa a)K_1(\kappa b)]}}-2a{b}^2{H}_2+2a^3{H}_4$

(15)$+b^2{f}_1(b)H_0+a^2{f}_2(b)H_0-2ab{K}_1(Ab)[a{H}_8-f_3(b)H_0]-2ab{I}_1(Ab)[a{H}_6$

$+f_4(b)H_0]\}{\{b^2{H}_1-a^2{H}_3+ab[I_1(Ab)H_5+K_1(Ab)H_7]\}}^{-1}.$

The bulk fluid velocity can be obtained by performing the inverse Laplace transform of Equation (15) using numerical methods [29,30]. Note that the zeta potential $\zeta$ of the cylinder in Equations (10), (13) and (15) can be replaced by its surface charge density $\sigma$ through their relationship given in Equation (2).

3. Results and Discussion

The normalized fluid velocity components $-\eta v_r/2\epsilon \zeta E_{\infty }\cos \theta$ and $\eta v_{\theta }/2\epsilon \zeta E_{\infty }\sin \theta$ of the transient transverse electroosmosis of an electrolyte solution around an unconfined charged circular cylinder (namely $\phi ={(a/b)}^2=0$) in the radial and angular directions, respectively, are plotted against the scaled radial coordinate $r/a$ in Figure 2 for several values of the scaled time $\nu t/a^2$ after the application of the electric field and the ratio of the cylinder radius to the Debye length $\kappa a$. These normalized fluid velocity components around the circular cylinder in a unit cell with the solid volume fraction $\phi =0.04$ are also plotted against $r/a$ in Figure 3 with the same values of $\nu t/a^2$ and $\kappa a$ for both the Happel and Kuwabara cell models with the Dirichlet condition of the applied electric potential distribution ${\psi }_{\mathrm{a}}(r,\theta )$ at the outer boundary $r=b$ of the cell ($\chi =1+\varphi$). For the given values of $r/a$, $\kappa a$, and $\phi$, the radial component $-\eta v_r/2\epsilon \zeta E_{\infty }\cos \theta$ is always a positive increasing function of $\nu t/a^2$ to a finite constant in the steady limit $\nu t/a^2\to \infty$ as expected (as shown in Figure 2a and Figure 3a), while the angular component $\eta v_{\theta }/2\epsilon \zeta E_{\infty }\sin \theta$ is a positive increasing function of $\nu t/a^2$ only at small values of $r/a$ (locations near the cylinder surface) and may become negative and decrease with an increase in $\nu t/a^2$ at greater values of $r/a$ (as shown in Figure 2b and Figure 3b). For the fixed values of $\nu t/a^2$, $\kappa a$, and $\phi$, as expected, both velocity components in general increase with an increase in $r/a$ from zero at the cylinder surface ($r/a=1$), reach their maxima (whose locations shift toward larger $r/a$ at a greater $\nu t/a^2$ or a smaller $\kappa a$), and then decrease with a further increase in $r/a$ (interestingly, $\eta v_{\theta }/2\epsilon \zeta E_{\infty }\sin \theta$ may thereafter become negative and reach a minimum, whose location also shifts toward larger $r/a$ at a greater $\nu t/a^2$ or a smaller $\kappa a$, and then increase with a further increase in $r/a$). In the steady limit $\nu t/a^2\to \infty$, the radial component $-\eta v_r/2\epsilon \zeta E_{\infty }\cos \theta$ is a monotonically increasing function of $r/a$ and the angular component $\eta v_{\theta }/2\epsilon \zeta E_{\infty }\sin \theta$ is always positive. The 2D map of the electroosmotic flow field around a circular cylinder is similar to the corresponding map of the creeping flow field around a sphere, as shown in Happel and Brenner [28]. For specified values of $\nu t/a^2$, $r/a$, and $\phi$, both velocity components generally increase with an increase in $\kappa a$, and there are some exceptions for $\eta v_{\theta }/2\epsilon \zeta E_{\infty }\sin \theta$. The Happel and Kuwabara cell models lead to qualitatively alike and quantitatively equivalent results of the normalized fluid velocity components (as shown in Figure 3).

The normalized fluid velocity $\eta U/\epsilon \zeta E_{\infty }$ of the bulk electrolyte solution undergoing transient transverse electroosmosis in a fibrous porous medium is plotted against the scaled time $\nu t/a^2$ and fiber volume fraction $\varphi$ in Figure 4 and Figure 5, respectively, for the different values of the ratio of the cylinder radius to the Debye length $\kappa a$ for both the Happel and Kuwabara cell models with the Dirichlet condition ($\chi =1+\varphi$). It is understood that $\eta U/\epsilon \zeta E_{\infty }$ equals unity in the limiting case of $\nu t/a^2\to \infty$, $\phi =0$, and $\kappa a\to \infty$ [26]. The effect of $\varphi$ on the continuous growth of the electroosmotic velocity with time is substantial and complicated, and a greater $\nu t/a^2$ is needed for $\eta U/\epsilon \zeta E_{\infty }$ to attain a definite percentage of its terminal value (as $\nu t/a^2\to \infty$) for a fiber matrix with a greater porosity (smaller $\varphi$). For a given value of $\kappa a$, this terminal value generally increases with a decrease in $\varphi$ (an increase in the porosity) as expected [22], while $\eta U/\epsilon \zeta E_{\infty }$ at an early time may decrease with a decrease in $\varphi$, since the characteristic time of the transient electroosmosis is greater and the fluid velocity of the bulk electrolyte solution at the outer edge of the unit cell is hardly influenced by the transient electroosmosis starting from near the charged surface of the fiber cylinder for the porous medium with a smaller $\varphi$. The Happel model in general results in a greater value of the normalized bulk fluid velocity than the Kuwabara model does, but there are exceptions (qualitatively different terminal velocities may be obtained from them, as shown for a thicker double layer and a small solid volume fraction in Figure 5a). Since the Happel model has a significant advantage in that it does not require an exchange of mechanical energy between the unit cell and the environment, the Happel model should be more accurate than the Kuwabara model.

In Figure 6, the normalized transient electroosmotic velocity $\eta U/\epsilon \zeta E_{\infty }$ of the bulk electrolyte solution in a fibrous porous medium with $\varphi =0.25$ is plotted against the scaled time $\nu t/a^2$ and ratio of the cylinder radius to the Debye length $\kappa a$. Again, $\eta U/\epsilon \zeta E_{\infty }$ grows incessantly with $\nu t/a^2$ for a given value of $\kappa a$. For fixed values of $\nu t/a^2$ and $\varphi$, this velocity is a monotonically increasing function of $\kappa a$. In comparison with our results in Figure 4, Figure 5 and Figure 6 for an arbitrary value of $\kappa a$, the calculations in the thin-double-layer approximation to the starting electroosmosis in the fibrous porous medium with $\kappa a\ge 10$ [26] overestimates the electroosmotic velocity of the fluid with the same value of $\kappa a$.

4. Concluding Remarks

In this paper, the start-up of the electroosmosis of electrolyte solutions upon the stepwise application of an electric field transversely to a charged fibrous porous medium consisting of parallel circular cylinders with arbitrary electric double layers is analytically investigated. The fluid momentum equation is solved by using both the Happel and Kuwabara unit cell models. A closed-form expression for the transient electroosmotic velocity of the bulk fluid $\eta U/\epsilon \zeta E_{\infty }$ in the Laplace transform is obtained as a function of the scaled elapsed time $\nu t/a^2$, ratio of the cylinder radius to the Debye screening length $\kappa a$, and the porosity of the fiber matrix $1-\varphi$. The effect of $\phi$ on the continuous growth of $\eta U/\epsilon \zeta E_{\infty }$ over time is substantial and complex. For a fiber matrix with smaller $\phi$ (larger porosity), $\eta U/\epsilon \zeta E_{\infty }$ takes longer to reach a certain percentage of its final value. Although the final value of $\eta U/\epsilon \zeta E_{\infty }$ generally decreases with increasing $\phi$, early velocities may increase with increasing $\phi$. For a given $\phi$, the transient electroosmotic velocity is a monotonically increasing function of the ratio of $\kappa a$.

The starting electroosmosis of electrolyte solutions in the axial direction of the fibrous porous medium composed of parallel circular cylinders was examined by using the unit cell model and an analytical solution of the electroosmotic velocity was obtained for arbitrary $\kappa a$, $\varphi$, and $\nu t/a^2$ [25]. For the general problem of transient electroosmosis in the fibrous medium caused by an external electric field in an arbitrary direction with respect to the axes of its constitutive cylinders, the solution of the electroosmotic velocity can be obtained by the vectorial addition of the previous result for the axial flow and the present result for the transverse flow due to the linearity.

In a real system, there should be a slab of fibrous material between two reservoirs and a voltage should be applied. Some entrance effect may be considered in our calculations for such a realistic case to show how large the electroosmotic flow may be and to provide some ideas on how the entrance effect may enter in the problem, like in the case of other electrokinetic phenomena [31,32,33].

Author Contributions

Conceptualization, H.J.K.; methodology, H.J.K. and W.Z.C.; investigation, H.J.K. and W.Z.C.; writing—original draft preparation, H.J.K. and W.Z.C.; writing—review and editing, H.J.K.; supervision, H.J.K.; funding acquisition, H.J.K. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures

Figure 1. Geometric illustration of the unit cell model for the transient electroosmosis in an array of parallel circular cylinders.

Figure 2. Normalized fluid velocity components of the transient electroosmosis around an unconfined circular cylinder ([Forumla omitted. See PDF.]), versus the scaled radial coordinate [Forumla omitted. See PDF.] for various values of the scaled time [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.]; (b) [Forumla omitted. See PDF.]. The solid and dashed curves denote [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.], respectively.

Figure 3. Normalized fluid velocity components of the transient electroosmosis around the circular cylinder in a unit cell with [Forumla omitted. See PDF.] and [Forumla omitted. See PDF.] versus the scaled radial coordinate [Forumla omitted. See PDF.] for the various values of the scaled time [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.]; (b) [Forumla omitted. See PDF.]. The solid and dashed curves denote the Happel and Kuwabara cell models, respectively.

Figure 4. Normalized bulk fluid velocity [Forumla omitted. See PDF.] of the transient electroosmosis in a fibrous porous medium versus the scaled time [Forumla omitted. See PDF.] for various values of the solid volume fraction [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.]; (b) [Forumla omitted. See PDF.]. The solid and dashed curves denote the Happel and Kuwabara cell models, respectively.

Figure 5. Normalized bulk fluid velocity [Forumla omitted. See PDF.] of the transient electroosmosis in a fibrous porous medium versus the solid volume fraction [Forumla omitted. See PDF.] for various values of the scaled time [Forumla omitted. See PDF.]: (a) [Forumla omitted. See PDF.]; (b) [Forumla omitted. See PDF.]. The solid and dashed curves denote the Happel and Kuwabara cell models, respectively.

Figure 6. Normalized bulk fluid velocity [Forumla omitted. See PDF.] of the transient electroosmosis in a fibrous porous medium with the solid volume fraction [Forumla omitted. See PDF.]: (a) versus [Forumla omitted. See PDF.] for various values of [Forumla omitted. See PDF.]; (b) versus [Forumla omitted. See PDF.] for various values of [Forumla omitted. See PDF.]. The solid and dashed curves denote the Happel and Kuwabara cell models, respectively.

Appendix A

The coefficients $H_0$, $H_1$, …, and $H_8$ in Equation (13) are given as follows:$$H_0=8-2A{I}_0(Ab)[A(a^2-b^2)K_0(Aa)+2a{K}_1(Aa)]+Ab{I}_1(Ab)[d{K}_0(Aa)+2Aa{K}_1(Aa)]$$$$-A\{2a{I}_1(Aa)R-I_0(Aa)[2A(a^2-b^2)K_0(Ab)+bd{K}_1(Ab)]\},$$$$H_1=4+Ab\{M{K}_0(Aa)-I_0(Aa)[2Ab{K}_0(Ab)+N{K}_1(Ab)]\},$$$$H_2={(2a)}^{-1}\{f_1(b)H_1+Aa[f_2(b)Aa\{Ab{I}_1(Ab)K_0(Aa)-2I_0(Ab)K_0(Aa)+I_0(Aa)R\}+2f_3(b)Q+2f_4(b)P]\},$$$$H_3=b\{Aa[bM{K}_1(Aa)+bN{I}_1(Aa)-2a]K_2(Aa)-N(b-A{a}^{2}T\left)\right\}{\{a[a{K}_1(Ab)-b{K}_1(Aa)]\}}^{-1},$$$$H_4=[16\{b{K}_1(Aa)-a{K}_1(Ab)\left\}\right(b-A{a}^{2}T\left)\right\{2a{f}_2(b)-b{f}_3(b)N-b{f}_4(b)M\}$$$$+8N(b-A{a}^{2}T)\{a[a{f}_2(b)-2b{f}_3(b)K_1(Ab)+2b{f}_4(b)I_1(Ab)]-b^2{f}_1(b)\}$$$$+8Aa{K}_2(Aa)\{bM{K}_1(Aa)+bN{I}_1(Aa)-2a\}\{b^2{f}_1(b)-a[a{f}_2(b)-2b{f}_3(b)K_1(Ab)$$$$+2b{f}_4(b)I_1(Ab)]\}]{\{16a^{2}b[a{K}_1(Ab)-b{K}_1(Aa)]\}}^{-1},$$$$H_5=2a^{-1}[bN-2A{a}^2{K}_2(Aa)],$$$$H_6=a^{-1}[f_1(b)H_5/2+Aa{f}_2(b)Q+A{f}_3(b)\{2a{K}_1(Aa)R+K_0(Aa)[2A(a^2-b^2)K_0(Ab)+bd{K}_1(Ab)]\}$$$$+f_4(b)\{2A{I}_0(Ab)[A(a^2-b^2)K_0(Aa)+2a{K}_1(Aa)]+Ab{I}_1(Ab)[d{K}_0(Aa)-2Aa{K}_1(Aa)]-4\}],$$$$H_7=2a^{-1}[bM-2A{a}^2{I}_2(Aa)],$$$$H_8=a^{-1}[f_1(b)H_7/2+f_2(b)AaP+f_3(b)\{4-2aA{I}_1(Aa)R+A{I}_0(Aa)[2A(a^2-b^2)K_0(Ab)-bd{K}_1(Ab)]\}$$$$+A{f}_4(b)\{2a{I}_1(Aa)[-2I_0(Ab)+Ab{I}_1(Ab)]+I_0(Aa)[2A(a^2-b^2)I_0(Ab)+bd{I}_1(Ab)]\}],$$ for the Happel model, and$$H_0=a^{-1}\{I_1(Ab)[A(a^2-b^2)K_0(Aa)+2a{K}_1(Aa)]+K_1(Ab)[A(a^2-b^2)I_0(Aa)-2a{I}_1(Aa)]\},$$$$H_1=-A{a}^{-1}{b}^2[I_1(Ab)K_0(Aa)+I_0(Aa)K_1(Ab)],$$$$H_2={(2a^2)}^{-1}\{A{I}_0(Aa)K_1(Ab)[-b^2{f}_1(b)+a^2{f}_2(b)]-2a{K}_1(Ab)f_3(b)$$$$+I_1(Ab)[A{K}_0(Aa)\{a^2{f}_2(b)-b^2{f}_1(b)\}+2a{f}_4(b)]\},$$$$H_3=A{a}^{-1}{b}^{2}T,$$$$H_4={(2a^3)}^{-1}\{Aa[b^2{f}_1(b)-a^2{f}_2(b)]T+2b^2{f}_3(b)K_1(Ab)-2b^2{f}_4(b)I_1(Ab)\},$$$$H_5=2a^{-2}{b}^2{K}_1(Ab),$$$$H_6=a^{-3}\{b^2{f}_1(b)K_1(Ab)+a[\{A(a^2-b^2)K_0(Aa)+2a{K}_1(Aa)\}\{f_3(b)K_1(Ab)-f_4(b)I_1(Ab)\}-a{f}_2(b)K_1(Ab)]\},$$$$H_7=-2a^{-2}{b}^2{I}_1(Ab),$$$$H_8=a^{-3}\{a[a{f}_2(b)I_1(Ab)+\{A(a^2-b^2)I_0(Aa)-2a{I}_1(Aa)\}\{f_3(b)K_1(Ab)-f_4(b)I_1(Ab)\}]-b^2{f}_1(b)I_1(Ab)\},$$ for the Kuwabara model, where$$M=2Ab{I}_0(Ab)+(4+A^2{b}^2)I_1(Ab),$$$$N=2Ab{K}_0(Ab)+(4+A^2{b}^2)K_1(Ab),$$$$P=2I_0(Aa)-2I_0(Ab)+Ab{I}_1(Ab),$$$$Q=2K_0(Aa)-2K_0(Ab)-Ab{K}_1(Ab),$$$$R=2K_0(Ab)+Ab{K}_1(Ab),$$$$T=I_2(Aa)K_1(Ab)+I_1(Ab)K_2(Aa),$$$$d=A^2(a^2-b^2)-4$$

References

1. Rice, C.L.; Whitehead, R. Electrokinetic flow in a narrow cylindrical capillary. J. Phys. Chem.; 1965; 69, pp. 4017-4024. [DOI: https://dx.doi.org/10.1021/j100895a062]

2. Squires, T.M.; Bazant, M.Z. Induced-charge electro-osmosis. J. Fluid Mech.; 2004; 509, pp. 217-252. [DOI: https://dx.doi.org/10.1017/S0022112004009309]

3. Mao, M.; Sherwood, J.D.; Ghosal, S. Electro-osmotic flow through a nanopore. J. Fluid Mech.; 2014; 749, pp. 167-183. [DOI: https://dx.doi.org/10.1017/jfm.2014.214]

4. Di Muccio, G.; della Rocca, B.M.; Chinappi, M. Geometrically induced selectivity and unidirectional electroosmosis in uncharged nanopores. ACS Nano; 2022; 16, pp. 8716-8728. [DOI: https://dx.doi.org/10.1021/acsnano.1c03017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35587777]

5. Hu, Z.; Zhao, T.; Wang, H.; Zhao, W.; Wang, K.; Bai, J.; Wang, G. Asymmetric temporal variation of oscillating AC electroosmosis with a steady pressure-driven flow. Exp. Fluids; 2020; 61, 233. [DOI: https://dx.doi.org/10.1007/s00348-020-03060-z]

6. Ding, H.; Zhong, X.; Liu, B.; Shi, L.; Zhou, T.; Zhu, Y. Mixing mechanism of a straight channel micromixer based on light-actuated oscillating electroosmosis in low-frequency sinusoidal AC electric field. Microfluid. Nanofluid.; 2021; 25, 26. [DOI: https://dx.doi.org/10.1007/s10404-021-02430-1]

7. Yu, Y.; Li, J.-C.; Lin, H.; Li, K.; Yi, F. Study on ethanol driven by alternating current electroosmosis in microchannels. Sens. Actuators A Phys.; 2023; 351, 114174. [DOI: https://dx.doi.org/10.1016/j.sna.2023.114174]

8. Dutta, D.; Smith, K.; Palmer, X. Long-range ACEO phenomena in microfluidic channel. Surfaces; 2023; 6, pp. 145-163. [DOI: https://dx.doi.org/10.3390/surfaces6020011]

9. Yossifon, G.; Frankel, I.; Miloh, T. Macro-scale description of transient electro-kinetic phenomena over polarizable dielectric solids. J. Fluid Mech.; 2009; 620, pp. 241-262. [DOI: https://dx.doi.org/10.1017/S002211200800459X]

10. Hanna, W.T.; Osterle, J.F. Transient electro-osmosis in capillary tubes. J. Chem. Phys.; 1968; 49, pp. 4062-4068. [DOI: https://dx.doi.org/10.1063/1.1670718]

11. Deng, S.; Xiao, T. Transient two-layer electroosmotic flow and heat transfer of power-law nanofluids in a microchannel. Micromachines; 2022; 13, 405. [DOI: https://dx.doi.org/10.3390/mi13030405]

12. Keh, H.J.; Tseng, H.C. Transient electrokinetic flow in fine capillaries. J. Colloid Interface Sci.; 2001; 24, pp. 450-459. [DOI: https://dx.doi.org/10.1006/jcis.2001.7797]

13. Jian, Y.; Yang, L.; Liu, Q. Time periodic electro-osmotic flows through a microannulus. Phys. Fluids; 2010; 22, 042001. [DOI: https://dx.doi.org/10.1063/1.3358473]

14. Sin, J.-S.; Ri, N.-I.; Kim, H.-C.; Hyon, S.-H. Impact of charged soft layers on electroosmosis of Maxwell fluids in soft nanochannels. Phys. Fluids; 2023; 35, 112002. [DOI: https://dx.doi.org/10.1063/5.0170703]

15. Levine, S.; Neale, G.H. The prediction of electrokinetic phenomena within multiparticle systems I. Electrophoresis and electroosmosis. J. Colloid Interface Sci.; 1974; 47, pp. 520-529. [DOI: https://dx.doi.org/10.1016/0021-9797(74)90284-7]

16. Zharkikh, N.I.; Shilov, V.N. Theory of collective electrophoresis of spherical particles in the Henry approximation. Colloid J. USSR Engl. Transl.; 1982; 43, pp. 865-870.

17. Kozak, M.W.; Davis, E.J. Electrokinetic phenomena in fibrous porous media. J. Colloid Interface Sci.; 1986; 112, pp. 403-411. [DOI: https://dx.doi.org/10.1016/0021-9797(86)90108-6]

18. Ohshima, H. Electroosmotic velocity in fibrous porous media. J. Colloid Interface Sci.; 1999; 210, pp. 397-399. [DOI: https://dx.doi.org/10.1006/jcis.1998.5983]

19. Keh, H.J.; Wei, Y.K. Diffusioosmosis and electroosmosis of electrolyte solutions in fibrous porous media. J. Colloid Interface Sci.; 2002; 252, pp. 354-364. [DOI: https://dx.doi.org/10.1006/jcis.2002.8483]

20. Zholkovskiy, E.K.; Masliyah, J.H.; Shilov, V.N.; Bhattacharjee, S. Electrokinetic phenomena in concentrated disperse systems: General problem formulation and spherical cell approach. Adv. Colloid Interface Sci.; 2007; 134–135, pp. 279-321. [DOI: https://dx.doi.org/10.1016/j.cis.2007.04.025] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17599797]

21. Ahualli, S.; Jimenez, M.L.; Carrique, F.; Delgado, A.V. AC electrokinetics of concentrated suspensions of soft particles. Langmuir; 2009; 25, pp. 1986-1997. [DOI: https://dx.doi.org/10.1021/la803171f] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19199730]

22. Keh, H.J.; Wu, Y.Y. Electroosmotic velocity and electric conductivity in a fibrous porous medium in the transverse direction. J. Phys. Chem. B; 2011; 115, pp. 9168-9178. [DOI: https://dx.doi.org/10.1021/jp203890n]

23. Su, Y.W.; Keh, H.J. Electrokinetic flow of salt-free solutions in a fibrous porous medium. J. Phys. Chem. B; 2019; 123, pp. 9724-9730. [DOI: https://dx.doi.org/10.1021/acs.jpcb.9b09008]

24. Chen, W.L.; Keh, H.J. Electroosmosis and electric conduction of electrolyte solutions in charge-regulating fibrous porous media. Colloids Interfaces; 2021; 5, 19. [DOI: https://dx.doi.org/10.3390/colloids5010019]

25. Chen, G.Y.; Keh, H.J. Start-up of electrokinetic flow in a fibrous porous medium. J. Phys. Chem. C; 2014; 118, pp. 2826-2833. [DOI: https://dx.doi.org/10.1021/jp4110568]

26. Chiang, C.C.; Keh, H.J. Transient electroosmosis in the transverse direction of a fibrous porous medium. Colloids Surf. A; 2015; 481, pp. 577-582. [DOI: https://dx.doi.org/10.1016/j.colsurfa.2015.06.021]

27. Li, M.X.; Keh, H.J. Start-up electrophoresis of a cylindrical particle with arbitrary double layer thickness. J. Phys. Chem. B; 2020; 124, pp. 9967-9973. [DOI: https://dx.doi.org/10.1021/acs.jpcb.0c07436] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33085892]

28. Happel, J.; Brenner, H. Low Reynolds Number Hydrodynamics; Nijhoff: Dordrecht, The Netherlands, 1983.

29. Talbot, A. The accurate numerical inversion of Laplace transforms. IMA J. Appl. Math.; 1979; 23, pp. 97-120. [DOI: https://dx.doi.org/10.1093/imamat/23.1.97]

30. Murli, A.; Rizzardi, M. Talbot’s method for the Laplace inversion problem. ACM Trans. Math. Softw.; 1990; 16, pp. 158-168. [DOI: https://dx.doi.org/10.1145/78928.78932]

31. Gadaleta, A.; Sempere, C.; Gravelle, S.; Siria, A.; Fulcrand, R.; Ybert, C.; Bocquet, L. Sub-additive ionic transport across arrays of solid-state nanopores. Phys. Fluids; 2014; 26, 012005. [DOI: https://dx.doi.org/10.1063/1.4863206]

32. Yaroshchuk, A.; Bondarenko, M.P. Interaction of potential sources in infinite 2D arrays: Diffusion through composite membranes, micro-electrochemistry, entrance resistance, and other examples. Adv. Theory Simul.; 2021; 4, 2100128. [DOI: https://dx.doi.org/10.1002/adts.202100128]

33. Baldelli, M.; Di Muccio, G.; Viola, F.; Giacomello, A.; Cecconi, F.; Balme, S.; Chinappi, M. Performance of single nanopore and multi-pore membranes for blue energy. ChemPhysChem; 2024; 25, e202400395. [DOI: https://dx.doi.org/10.1002/cphc.202400395]