Analysis of Transient Pulse Electroosmotic Flow of Maxwell Fluid through a Circular Micro-Channel Using Laplace Transform Method ()
1. Introduction
In recent years, microfluidic devices have been vigorously developed and applied in micro-electronic mechanical systems and microbiological sensors, and the electroosmotic flow (EOF) formed in these devices has become more and more attractive [1] [2]. Different from the flow in the conventional macro-sized channel, the existence of the electric double layer (EDL) must be considered when the flow is analyzed in the micro-channel [3]. Generally, when most substances come into contact with polar solutions, they tend to generate negative charges on the surface. The distribution of ions near the wall in the solution will be affected by this phenomenon. The opposite ions with the opposite polarity to the wall will be attracted to the wall, while the same ions will be repelled away from the wall. In this way, an electric double layer (EDL) will be formed [4]. Further, when an external electric field is applied to both ends of the channel, the ions in the EDL will move under the force of the electric field. This is due to the viscosity of the fluid itself, the moving free ions will drive the movement of nearby fluid clusters, and eventually form an electroosmotic flow (EOF).
Among previous research, various theoretical and experimental studies on steady EOF of Newtonian fluids in micro-channels under different geometric regions and physical conditions have been carried out [5] - [10]. Recently, time-dependent EOF as an alternative mechanism of microfluidic transport has attracted more and more attention [11] - [17].
All of the above-mentioned studies are related to Newtonian fluids. However, most of solutions of industry and biopharmaceutical are fluid that has the structural characteristic of non-Newtonian fluids, for example biological fluid and other solutions of long-chain molecules, whose structural characteristics include strain force, normal shear stress, hysteresis effect, variable viscosity, memory effect and so on [18]. Therefore, these fluids cannot be regarded as Newtonian fluids. From many literatures, we know that the research on non-Newtonian fluids at home and abroad mainly focuses on the rheological properties of non-Newtonian fluids, non-Newtonian fluids gas-liquid two-phase flow, non-Newtonian fluids mass and heat transfer, and numerical simulation methods of non-Newtonian fluids. Furthermore, the theoretical research of EOF of non-Newtonian fluids is mainly limited to simple fluid models because of the inherent analytical difficulties introduced by more complex constitutive equations [19]. So far, some work has been done on the simple non-Newtonian fluids models, such as Oldroyd-B fluid model [20] [21] [22], Power-law fluid model [23] [24] [25], Maxwell fluid model [26] - [32] and Jeffrey fluid model [33] [34] [35].
Although we have learned some basic characteristics of EOF of non-Newtonian fluids in the above research, its rich characteristics still need to be studied. Recent study have shown that Maxwell fluid model simulation of blood in narrow conical vessels has achieved an ideal effect, and it is completely possible to analyze the blood-based microfluidics and other microbial fluid transmission systems by means of electric mechanism [26]. Meanwhile, the literature [36] found that low-frequency pulses can stimulate nerve cell tissues, relieve pain and promote local blood circulation, which has been confirmed in clinical medical research. Additionally, the circular channel is very common in practical applications [34]. So we select Maxwell fluid model as a constitutive relation to describe pulse EOF in the circular micro-channel in present work.
Besides, Liu et al. [37] studied the alternating current (AC) EOF of generalized Maxwell fluid through a circular micro-channel and an analytical solution of EOF velocity distribution is derived. Taking into account the limitations of alternating current (such as continuous and time-varying), as well as the wide application of pulse current in engineering in recent years [38] [39] [40], thus another important purpose of this paper is also extending the AC EOF to pulse current (PC) EOF of the generalized Maxwell fluid model in the circular micro-channel. With the help of Laplace transform method, a semi-analytical solution for pulse EOF velocity is obtained by solving the general Cauchy momentum equation. Moreover, we analyzed the influence of several parameters such as the relaxation time, the pulse width and the electrokinetic width on the pulse EOF of Maxwell fluid. The rest of this paper is organized as follows: The physical description of the problem and the analytical solution to the equations governing the pulse EOF of the Maxwell fluid model in a circular micro-channel are given in Section 2. Section 3 discusses the numerical results of the study and the conclusions are presented in Section 4.
2. Mathematical Formulation
2.1. Cauchy Momentum Equation and Constitutive Relation
The transient pulse EOF of an incompressible Maxwell fluid through a circular micro-channel with radius R, the length of the channel is L, assumed to be much larger than the diameteri.e.,
is sketched in Figure 1(a). The pulse EOF is pumped through an ideal rectangle pulse electric field of strength E0 with pulse amplitude of 1, pulse repetition period of
and pulse width of
in Figure 1(b). The ideal rectangle pulse can be expressed as the following form [41]
(1)
Due to the symmetry of the geometry, we only study the semi-section of the micro-channel. Considering the pure pulse EOF and ignoring the pressure gradient, the one-dimensional Cauchy momentum equation can be given as
![]()
Figure 1. Schematic of the transient pulse EOF of Maxwell fluid in a circular micro-channel. (a) Circular micro-channel; (b) Time periodic pulse wave (ideal rectangle pulse wave).
(2)
where
is the velocity along z axial direction,
is the fluid density, t is the time,
is the stress tensor and
is the volume charge density,
is the ideal rectangle pulse electric field of strength E0.
The boundary condition of Equation (2) is supposed no slip and can be written as
,
(3)
For the Maxwell fluid, the constitutive equation satisfies [42]
(4)
where
and
are the relaxation time and the zero shear rate viscosity, respectively.
2.2. Electric Potential Field Solution
The chemical interaction of electrolyte liquid and solid wall generates an electric double layer (EDL), a very thin charged liquid layer at the solid-liquid interface. A cylindrical coordinate system
is introduced. In this theoretical model, it is assumed that the channel wall is uniformly charged, so that the electrical potential in the EDL only varies in this r direction and does not depend on
[43]. For a symmetric binary electrolyte solution, we assume that the electrical potential
of the EDL is stable, and its distribution and the local volumetric net charge density
are given by the Poisson-Boltzmann (P-B) equation
(5)
(6)
where
is the dielectric constant of the electrolyte liquid,
is the electrical potential of the EDL, n0 is the ion density of the bulk liquid,
is the valence, e0 is the electron charge, kb is the Boltzmann constant, T is the absolute temperature and sinh is a sine function.
Combining Equations (5) and (6) gives
(7)
which is subject to the following boundary conditions
,
(8)
where
is wall zeta potential, r is radial coordinate and R is radius of the circular micro-channel.
Provided that the electrical potential is small enough, the Debye-Hückel linearization approximation can be applied, which means physically that the electrical potential is small compared to the thermal energy of the charged species [33]. Then Equation (7) can be simplified as
, where
(9)
where
is the Debye-Hückel parameter and
usually denotes the thickness of the EDL in physical.
The net charge density can be obtained by solving Equation (9) with boundary condition of Equation (8)
(10)
where I0 is first kind modified Bessel function of order zero.
2.3. The Analytical Solutions of the Cauchy Momentum Equation
In order to obtain the solution of velocity field, some dimensionless parameters are given as
(11)
where Ueo denotes dimensionless steady Helmholtz-Smoluchowshi EOF velocity of Newtonian fluids, K is the ratio of the characteristic width of the micro-channel to Debye length.
Using Equation (11), the Equations of (2) and (4) and the corresponding boundary condition (3) can be written as
(12)
(13)
,
(14)
Eliminating
from Equation (12) and Equation (13) yields
(15)
Let us employ the method of Laplace transform defined by
(16)
Obviously
, and the Laplace transform of
is given by the
Appendix. If initial condition satisfies
, the Laplace transforms of Equation (15) and boundary condition (14) can be given as
(17)
(18)
here tanhis a hyperbolic tangent function.
Equation (17) can be simplified as
(19)
where
.
Equation (19) is a linear and inhomogeneous ordinary differential equation, and its solution can be written as the sum of a general solution
corresponding to homogeneous equation and a special solution
.
(20)
Due to the finite of
at
, the homogeneous solution of Equation (19) is written as
(21)
here A is constant, which can be determined from boundary conditions of Equation (18).
Observing the formation of the right hand side of Equation (19), the special solution can be expressed as
(22)
Substituting Equation (22) into Equation (19) yields
(23)
From Equation (9) and Equation (11), we have
(24)
Inserting Equation (24) into Equation (23) and equalizing the coefficients in front of the modified Bessel functions
at the two sides of the equation yields
(25)
Therefore, the solution of velocity
can be given as
(26)
The coefficient A with boundary condition of Equation (18) can be determined as
(27)
Substituting Equation (27) into Equation (26), we can get
(28)
The inverse Laplace transform is defined by
(29)
where
is a vertical line to the right of all singularities of
in the complex s plane. Because of the complexity of the express of
, the numerical computation must be performed by numerical inverse Laplace transform [44]. The method of Laplace inverse transform in the above-mentioned literature needs to be supplemented by the fact that the Laplace inverse transform is integrated using the trapezoidal rule, and the Laplace numerical inverse transform is based on the Fourier series that accelerates the convergence.
3. Results and Discussion
In the second section, we have obtained the semi-analytical solution of the transient pulse EOF velocity of Maxwell fluid through a circular micro-channel, which mainly relies on relevant dimensionless parameters, such as the relaxation time
, the pulse width
and the electrokinetic width K. In this section, we will discuss their influence on the normalized pulse EOF velocity in detail.
Figure 2 shows the variations of normalized pulse EOF velocity with radius in one cycle under different relaxation time
(0.01, 0.05, 0.1 and 0.2). It can be seen from Figure 2 that for fixed electrokinetic width K and pulse width
, the amplitude of the normalized pulse EOF velocity profiles varies significantly with the increase of relaxation time
, and the velocity profiles varies as the time
increases. For fixed relaxation time
, the fluid velocity will gradually stabilize with the increase of time
within a half period. The larger the relaxation time
, the larger the amplitude of the velocity profiles, and the longer the time it takes for the fluid to reach a steady status. The main reason of this fact is that for larger relaxation time, the distorted polymer molecules will have no time to change their structure during the time scale of the fluid flow, and the behavior of
the fluid gradually tends to become a Hookean elastic solid. Thus, the flow takes longer time to attain the steady status [45].
Figure 3 depicts the variations of normalized pulse EOF velocity with time for different relaxation time
(0.1, 0.3, 0.5 and 0.7). As expected, we can find that for fixed pulse width
, the amplitude of the normalized pulse EOF velocity profiles at the center of the micro-channel (that is,
) increases significantly with the increase of the relaxation time
, especially for the smaller pulse width
(see Figures 3(a)-(c)). However, for larger pulse width
, this change will be less obvious. This implies that increasing the pulse width will weaken the effect of relaxation time on velocity. The reason may be a larger pulse width
, which means that the longer the pulse force lasts, the stronger its stability. For fixed relaxation time
, the normalized pulse EOF velocity profiles of frequency are different due to the different pulse width
. With the pulse width
increases, the different change frequency of velocity profiles will slow down, which means a long cycle time [41], and the time required for the velocity profiles to reach a steady state becomes longer.
The variations of normalized pulse EOF velocity with radius for different relaxation time
(0.01, 0.1, 0.3 and 0.7) are presented in Figure 4. For a smaller
![]()
Figure 3. Variation of normalized pulse EOF velocity at different pulse width with time for different relaxation time
when
and
. (a)
; (b)
; (c)
; (d)
.
![]()
Figure 4. Variation of normalized pulse EOF velocity at different electrokinetic width with radius for different relaxation time
when
and
. (a)
; (b)
; (c)
; (d)
.
relaxation time
, it is clearly seen from Figure 4(a) that as the electrokinetic width K value increases, the normalized pulse EOF velocity variations are restricted to a very narrow area close to the EDL. However, with the increase of relaxation time
, the elasticity of the fluid becomes significant. And due to the elasticity of the physical property of the fluid as a whole, so the velocity variation can extend to the whole flow field (see Figures 4(b)-(d)) [46]. In addition, the amplitude of velocity increases with the gradual increase of relaxation time
.
4. Conclusions
A semi-analytical solution of the transient pulse EOF of Maxwell fluid through a circular micro-channel under the Debye-Hückel approximation is presented in this work. The computational results show that the velocity profiles depend mainly on the relaxation time
, the pulse width
and the electrokinetic width K. With the aid of inverse Laplace transform, the following conclusions can be drawn:
• The profiles of normalized pulse EOF velocity vary rapidly and gradually stabilize as the increase of time
within a half period (that is, pulse width
).
• Increasing relaxation time
will lead to larger velocity amplitude. In addition, the time needed to attain the steady status becomes longer with the increase of relaxation time
.
• The velocity profiles at the center of the micro-channel increase significantly with relaxation time
, especially for the smaller pulse width
. However, as the increase of pulse width
, this change will be less obvious. At the same time, the different change frequency of velocity profiles becomes slower, which means a long cycle time, and the time required for the fluid to reach a steady state becomes longer.
• With the increase of electrokinetic width K, the velocity variations are confined to a very narrow area close to the EDL for small relaxation time
. However, as the relaxation time
increases, the elasticity of the fluid becomes conspicuous and the velocity variations can be extended to the entire region of flow.
Acknowledgements
This work was supported by the Scientific Research Project of Inner Mongolia University of Technology (Grant No. ZZ201813).
Appendix
The Laplace transform of
is expressed as follows
(A.1)
here
denotes the period of the ideal rectangle pulse wave.
By shifting the term of Equation (A.1) and using the Equation (1), we have
(A.2)
where
(A.3)