Next Article in Journal
Research on Low-Cost Attitude Estimation for MINS/Dual-Antenna GNSS Integrated Navigation Method
Previous Article in Journal
Characterization of Al Incorporation into HfO2 Dielectric by Atomic Layer Deposition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermally Fully Developed Electroosmotic Flow of Power-Law Nanofluid in a Rectangular Microchannel

Institute of Architecture and Civil Engineering, Guangdong University of Petrochemical Technology, Maoming 525000, China
Micromachines 2019, 10(6), 363; https://doi.org/10.3390/mi10060363
Submission received: 6 May 2019 / Revised: 24 May 2019 / Accepted: 28 May 2019 / Published: 30 May 2019
(This article belongs to the Section A:Physics)

Abstract

:
The hydrodynamic and thermal behavior of the electroosmotic flow of power-law nanofluid is studied. A modified Cauchy momentum equation governing the hydrodynamic behavior of power-law nanofluid flow in a rectangular microchannel is firstly developed. To explore the thermal behavior of power-law nanofluid flow, the energy equation is developed, which is coupled to the velocity field. A numerical algorithm based on the Crank–Nicolson method and compact difference schemes is proposed, whereby the velocity, temperature, and Nusselt number are computed for different parameters. A larger nanoparticle volume fraction significantly reduces the velocity and enhances the temperature regardless of the base fluid rheology. The Nusselt number increases with the flow behavior index and with electrokinetic width when considering the surface heating effect, which decreases with the Joule heating parameter. The heat transfer rate of electroosmotic flow is enhanced for shear thickening nanofluids or at a greater nanoparticle volume fraction.

Graphical Abstract

1. Introduction

The development of the micro-electric mechanical system (MEMS) has received great attention because of its innovative application in chemical, medical, and biological-related industries. For instance, Lab-on-a-chip, which conducts experiments by processing bioliquids or chemical solutions on a microchip, has become increasingly popular [1]. Electroosmotic flow (EOF) being a major electrokinetic effect refers to the motion of an ionized liquid inside microchannels with respect to the stationary charged channel walls under an external electric field applied tangentially along the microchannel [2]. An electroosmosis-based micropump has been successfully developed [3] and became one of the most important components in lab-on-a-chip since it has advantages, such as the production of pulse-free and plug-like flow, and the dependence on non-mechanical parts [4] over the conventional pressure-driven flow (PDF). Therefore, a fundamental understanding of the mechanism of EOF is vital for the precise control and optimal design of microfluidic devices.
The hydrodynamic behavior of EOF of different fluids is investigated in different microchannels, such as slit [5], rectangular [6,7], elliptic [8], and circular microchannels [9]. The core of the attention on EOF has shifted to the heat transfer characteristics of EOF [9,10,11], two-layer EOF [12], rotating EOF [13,14], and pressure effects on EOF [15]. In biological and chemical industries, biofluids, such as blood and DNA, solutions manipulated in microfluidic devices show nonlinear rheological behavior, such as the viscosity dependent on shear rate, which cannot be modeled by the linear constitutional relation. Therefore, the theory of non-Newtonian fluids offers researchers in the fields of mathematics and physics challenges in developing solutions to the nonlinear governing equations. The EOF of Jeffrey fluid [16] and Maxwell fluid [17] under an alternating current (AC) electric field are investigated in terms of the influence of nonlinear rheological behavior on flow performance. The corresponding analytical solution for the velocity field has been developed by adopting the method of variable separation. Das and Charkraborty [5] first proposed a non-Newtonian power-law model for the EOF of biofluids and acquired analytical solutions of velocity, temperature, and concentration distribution. The power-law model has received much more attention because of the concise expression and breadth coverage for working liquids in a MEMS. Zhao et al. highlighted the dynamic properties of EOF of power-law fluid under an AC electric field along a rectangular microchannel [18]. My previous works conducted parametric studies on the EOF of power-law fluids in terms of the flow behavior index and electrokinetic width in a cylindrical circular microcapillary [19] and in a rectangular microchannel [20]. It showed that the nonlinear rheological behavior of power-law fluid has a significant influence on the flow pattern and volumetric flow rate of EOF. Particularly, the shear thinning fluid needs more developed temporal length and is more sensitive to the change in parameters than that of shear thickening fluid.
The thermal characteristics have to be carefully considered due to the involvement with fluid flow and the applications in the cooling systems of a MEMS. The Joule heating effect that is an inherent phenomenon in EOF resulting from the ohmic resistance of the electrolyte dominates the heat generation of fluid flow. It results in the change of temperature-dependent electrical/transport properties, such as viscosity and electric conductivity of the working liquids, which, in turn, influences the hydrodynamic behavior of EOF. Concerning the Joule heating-induced heat transfer of EOF, the models in a microcapillary [9] and in a rectangular microchannel [21] without consideration of viscous dissipation have been developed. A model in a slit microchannel with consideration of viscous dissipation and thermal radiation is developed by Shit et al. [10], which has been extended to power-law fluid in a circular microchannel [22]. The entropy generation rate for the EOF of a power-law fluid in a microchannel has been obtained, and it turns out that entropy generation is dominated by the Joule heating effect [23]. In addition, the microchannel flow can be assisted by the mixed effects of electroosmosis and pressure gradient. The hydrodynamic and thermal behavior of combined EOF and PDF of Newtonian fluid in a microchannel was studied by solving the velocity distribution, temperature distribution, and Nusselt number [24]. The model above has been extended to the case of power-law fluid in a slit microchannel [25] and in a rectangular microchannel [26], where the temperature distribution and Nusselt number have been numerically solved by taking into account the Joule heating effect and viscous dissipation. Lately, the effects of magnetic field on the heat generation of EOF has been given importance [27].
Compared with the conventional fluid, the nanofluid containing nanoparticles with a diameter of 1 to 100 nm shows larger overall thermal conductivity. This is because the thermal conductivity of solid particles, such as Cu, is 700 times higher than the thermal conductivity of water, and the thermal conductivity of metal oxide, such as AI2O3, is much greater than that of a pure fluid, and hence, the thermal performance of nanofluid is dramatically enhanced in the presence of these nanoparticles. Furthermore, the nanofluid is stable, which has no additional problems like pressure drop, sedimentation, or blockage in microchannels. Choi firstly proposed the model of nanofluid [28]. In the fields of energy, power, electronics, etc., the trend for using suspensions of nanoparticles in a fluid medium is growing [29], which has led to several researches on the nanofluid flow [30,31]. The comparison among several existing models for thermal conductivity was carried out to assist in the accurate prediction of thermal behavior of nanofluid in [32]. A more comprehensive survey on the turbulent heat transfer of nanofluid flow was conducted [33]. Since the characteristics of nanofluid mentioned above meet the increasingly higher requirements in heat transfer of heat exchange equipment of MEMS, the nanofluid flow in microscale has been paid more and more attention recently. The electrokinetic effects of nanofluids induced by the streaming potential on PDF in microtubes were theoretically studied [34]. The transient thermal characteristics of the combined EOF and PDF of nanofluids in a microchannel under the effect of a magnetic field were investigated [35]. Furthermore, as heat transfer performance of fluid flow is increasingly encountered in micro electrics, the base fluid is not only limited to Newtonian fluid, but also the emphasis has been laid on the non-Newtonian fluid flow. Therefore, in microscale devices with characteristics of compact structure, to achieve the efficient thermal management of the non-Newtonian fluid, the power-law nanofluid was developed, where nanoparticles are added, and the base fluid rheology is described by the power-law model. The EOF of power-law nanofluid in a parallel plate microchannel was studied where one-dimensional momentum equation and energy equation were analytically solved [36]. A mixed convection flow and the corresponding heat transfer of shear thinning power-law nanofluid was respectively investigated by Ellahi et al. [37] and by Si et al. [38]. It turns out that the solid volume fractions exert special effects on both the velocity field and temperature field of fluid flow.
An up-to-date literature review indicates that for non-Newtonian fluids, the studies on the thermal behavior of EOF lack detailed insights into the substance of microstructure, namely the role that nanoparticles play, and the influence of channel geometry, owing to the complexity of governing equations for non-Newtonian nanofluid EOF in complex microchannels. Moreover, the thermal behavior of pure EOF of a power-law fluid in complex microchannels has not been studied yet. Due to the practical and fundamental significance of EOF of power-law nanofluid in a complex microchannel, such as a rectangular channel, and motivated by the works reviewed above, it is necessary to provide a fundamental understanding of the Joule heating induced heat transfer characteristics. This paper develops a model for thermally fully developed EOF of power-law nanofluid in a rectangular microchannel. The viscosity of power-law nanofluid depends both on the flow behavior index and the volume fraction of nanoparticles. By numerically solving the modified momentum equation and the energy equation, the velocity distribution, temperature distribution, and Nusselt number under the influence of the Joule heating effect are evaluated for different nanoparticle volume fraction, electrokinetic width, flow behavior index. The results are useful for the optimal design and active control of microfluidic systems.

2. Mathematical Modeling

As sketched in Figure 1, a rectangular microchannel of width 2b and height 2a is considered. The channel is filled with power-law nanofluid with a dielectric constant ε and is subject to constant heat flux qs. The channel wall is uniformly charged with zeta potential ξ, and the electric double layers (EDLs) will not overlap. As an electric field is applied tangentially along the channel (z-direction), the liquid inside is set in motion due to the existence of free ions in EDLs. Owing to the symmetry, the following analysis is limited to a quarter of the channel, Ω , and the coordinate system used is shown in Figure 1.

2.1. Electric Potential Distribution

Based on the electrostatics theory, the relationship between the electric potential distribution ψ(x, y) and the local net charge density ρe is expressed by the Poisson equation.
2 ψ x 2 + 2 ψ y 2 = ρ e ε ε 0 .
The local net charge density ρe is subject to Boltzmann distribution.
ρ e = 2 χ e n 0 sinh ( χ e ψ k b T 0 ) ,
where χ denotes the valence of ions, e denotes the elementary charge, n0 is the ion density of the bulk liquid, kb is the Boltzmann constant, and T0 implies the absolute temperature. Substituting Equation (2) to Equation (1) yields the two-dimensional Poisson–Boltzmann (P-B) equation for the electric potential distribution of EDL in the rectangular microchannel,
2 ψ x 2 + 2 ψ y 2 = 2 χ e n 0 ε sinh ( χ e ψ k b T 0 ) .
The corresponding boundary conditions are given as
∂ψ/∂x| x = 0 = 0, ∂ψ/∂y|y = 0 = 0, ψ| x = b = ξ, ψ|y = a = ξ.
To facilitate the discussion, the following scales are introduced x ¯ = x / a , y ¯ = y / a , ξ ¯ = χ e ξ / ( k b T 0 ) , α = b/a, K = κ a , with K defined as the dimensionless electrokinetic width and 1/ κ as the thickness of EDL, where κ = [2χe2n0/(εkbT0)]1/2. Accordingly, the dimensionless P-B equation and the boundary conditions are obtained
2 ψ ¯ x ¯ 2 + 2 ψ ¯ y ¯ 2 = K 2 sinh ( ψ ¯ )
ψ ¯ / x ¯ | x ¯ = 0 = 0 , ψ ¯ / y ¯ | y ¯ = 0 = 0 ,   ψ ¯ | x ¯ = α = ξ ¯ ,   ψ ¯ | y ¯ = 1 = ξ ¯ .
For small zeta potential (ξ ≤ 0.025 V), the known Debye–Hückel linearization principle [5] sinh ( ψ ¯ ) ψ ¯ is used to linearize the P-B equation. In this case, an analytical solution to Equations (5) and (6) can be readily solved by using variable separation method as
ψ ¯ ( x ¯ , y ¯ ) = 2 ξ ¯ M = 1 ( 1 ) M + 1 cos ( μ M y ¯ ) cosh ( μ M 2 + K 2 x ¯ ) μ M cosh ( μ M 2 + K 2 α ) + 2 ξ ¯ N = 1 ( 1 ) N + 1 cos ( σ N x ¯ ) cosh ( σ N 2 + K 2 y ¯ ) α σ N cosh ( σ N 2 + K 2 )
where μ M = ( 2 M 1 ) π / 2 , σ N = ( 2 N 1 ) π / ( 2 α ) , M,N = 1,2,3,…

2.2. Velocity Distribution

The vector formation of Cauchy momentum equation reads
ρ ( V t + V · V ) = p + · τ + E ρ e
where V is the velocity field of EOF, ρ represents the density of mass, t denotes the variable of time, p is the pressure gradient, E is the electric field strength, τ ˙ = 2 μ e f f γ ˙ e denotes the stress tensor which can be expressed by a second invariant of the strain rate tensor eij = [(∂ui/∂xj) + (∂uj/∂xi)]/2 as γ ˙ = ( 2 e k l e k l ) 1 / 2 . μeff denotes the effective viscosity coefficient. For a fully developed, incompressible and laminar EOF of power-law nanofluid, it is assumed that there is only an axial velocity w = w(x, y) [7,20], accordingly, the effective viscosity of the power-law nanofluid is μeff = μf/(1 − ϕ)5/2. The viscosity of the base fluid is μf = η(|∂w/∂x|n − 1, |∂w/∂y|n − 1) which is a function of strain rate, flow consistency index η of dimension (Nm−2sn) [7,20], and volume fraction of nanoparticle ϕ [39]. Finally, one can express the constitutive equation of power-law nanofluid as
τ = μ e f f · γ ˙ = η ( 1 ϕ ) 5 / 2 ( | w x | n 1 w x + | w y | n 1 w y ) ,
where τ is the shear stress, and n is the flow behavior index. n = 1 represents a Newtonian nanofluid, n < 1 and n > 1 represents the shear thinning nanofluids (pseudoplastic nanofluids) and the shear thickening nanofluids (dilatant nanofluids), respectively. Consequently, the modified Cauchy momentum equation for power-law nanofluid EOF is
n η ( 1 ϕ ) 5 / 2 ( | w x | n 1 2 w x 2 + | w y | n 1 2 w y 2 ) + E ρ e = 0 .
The corresponding boundary conditions are
w / x | x = 0 = 0 , w / y | y = 0 = 0 , w | x = b = 0 , w | y = a = 0 .
Introducing the following scales w ¯ = w / W with W = ε k b E T 0 / ( μ 0 χ e ) where μ0 is the viscosity coefficient of Newtonian fluid of dimension (Nm−2s), one can obtain the dimensionless modified Cauchy momentum equation and the boundary conditions:
n η μ 0 ( 1 ϕ ) 5 / 2 ( W a ) n 1 ( | w ¯ x ¯ | n 1 2 w ¯ x ¯ 2 + | w ¯ y ¯ | n 1 2 w ¯ y ¯ 2 ) + K 2 sinh ( ψ ¯ ) = 0 ,
w ¯ / x ¯ | x ¯ = 0 = 0 , w ¯ / y ¯ | y ¯ = 0 = 0 , w ¯ | x ¯ = α = 0 , w ¯ | y ¯ = 1 = 0 .
When considering Newtonian nanofluid, i.e., n = 1, and using the known Debye–Hückel principle, an analytical solution can be obtained by the method of variable separation and the method of constant variation [40], which is expressed as
w ¯ ( x ¯ , y ¯ ) = I = 1 cos ( β I x ¯ ) U I ( y ¯ ) ,
where β I = ( 2 I 1 ) π / ( 2 α ) with I = 1,2,3,…, A = 8 α K 2 ( 1 ϕ ) 5 / 2 ξ ¯ ,
U I ( y ¯ ) = A C I cosh ( C I ) ( M = 1 A M I G M I D M + N = 1 B N I Q N I E N ) cosh ( C I y ¯ ) ,
+ A { M = 1 A M I [ cosh ( C I y ¯ ) cos ( μ M y ) ] D M ( μ M 2 + C I 2 ) + N = 1 B N I [ cosh ( σ N 2 + K 2 y ) cosh ( C I y ¯ ) ] E N ( σ N 2 + K 2 C I 2 ) } ,
A M I = [ μ M 2 + K 2 sinh ( μ M 2 + K 2 α ) cos ( β I α ) + β I cosh ( μ M 2 + K 2 α ) sin ( β I α ) ] μ M 2 + K 2 + β I ,
B N I = 1 2 { sin [ ( σ N + β I ) α ] σ N + β I + sin [ ( σ N β I ) α ] σ N β I } , C I 2 = β I D M = ( 1 ) M + 1 ( 2 M 1 ) π cosh ( μ M 2 + K 2 α ) , E N = ( 1 ) N + 1 ( 2 N 1 ) π cosh ( σ N 2 + K 2 ) , G M I = C I [ cosh ( C I ) cos ( μ M ) ] μ M 2 + C I 2 ,
Q N I = C I [ cosh ( σ N 2 + K 2 ) cosh ( C I ) ] σ N 2 + K 2 C I 2 .
In the case of non-Newtonian power-law nanofluid flow (n ≠ 1), a numerical algorithm is proposed to obtain the velocity field, which is presented in the next section.

2.3. Temperature Distribution

Assuming the thermophysical properties are constant, the heat generation induced by the viscous dissipation is much smaller than that induced by Joule heating [23,41] and thus, is neglected in the present analysis. The energy equation for the flow field in a rectangular microchannel is
( ρ c p ) e f f w T z = k e f f 2 T + σ E 2
with ( ρ c p ) e f f = ϕ ( ρ c p ) s + ( 1 ϕ ) ( ρ c p ) f and k e f f = k s + 2 k f + 2 ( k s k f ) ( 1 + ω ) 3 ϕ k s + 2 k f 2 ( k s k f ) ( 1 + ω ) 3 ϕ k f [32,42]. T denotes temperature field, ω implies the ratio of the nanolayer thickness to the original particle radius. σ, k, and (ρcp) denote electrical conductivity, thermal conductivity, and heat capacity of power-law nanofluid at the reference pressure, respectively. The subscripts s, f, and eff denote the solid nanoparticles, base fluid, and nanofluid, respectively. The second term at the right-hand side of Equation (13) denotes the volumetric heat generation owing to the Joule heating effect. The corresponding boundary conditions are
T / x | x = 0 = 0 , T / y | y = 0 = 0 , T / x | x = b = q s / k e f f , T / y | y = a = q s / k e f f   ( or   T | y = a = T w ( z ) ) ,
where qs = h(Tw − Tm) implies the constant wall heat flux and h is the convective heat transfer coefficient. Subscripts w and m denote the local wall value and mean value, respectively. As assumed in the preceding section, for thermally fully developed EOF, one has T / z = d T m / d z = d T w / d z = 0 and thus 2 T / z 2 = 0 . It is noted that the temperature variation with respect to z-direction cannot be neglected and thus, the temperature on the wall is a function of z. Plus, an overall energy balance at an elemental control volume on a length of duct dz gives:
d T m d z = q s ( a + b ) + σ E 2 a b ( ρ c p ) e f f a b w ¯ m
with the mean velocity w ¯ m = w m / W = 1 α 0 α 0 1 w ¯ d x ¯ d y ¯ .
Introducing T ¯ = k f ( T T w ) / ( q s a ) and the Joule heating parameter S = σ E 2 a / q s , namely, the ratio of Joule heating to the heat flux from the channel wall, the dimensionless energy equation and the corresponding boundary conditions are obtained as
2 T ¯ x ¯ 2 + 2 T ¯ y ¯ 2 = k f k e f f [ ( 1 + 1 α + S ) w ¯ w ¯ m S ] ,
T ¯ / x ¯ | x ¯ = 0 = 0 ,   T ¯ / y ¯ | y ¯ = 0 = 0 ,   T ¯ / x ¯ | x ¯ = α = k f / k e f f , T ¯ / y ¯ | y ¯ = 1 = k f / k e f f , ( or   T ¯ | y ¯ = 1 = 0 ) .
In the case of Newtonian fluid and using the known Debye–Hückel linearization principle, an analytical solution of temperature distribution is obtained by the method of variable separation and the method of constant variation,
T ¯ ( x ¯ , y ¯ ) = I = 1 cos ( β I x ¯ ) Y I ( y ¯ ) ,
where Y I ( y ¯ ) = L I ( 0 ) { k 1 A L I ( 1 ) L I ( 0 ) [ M = 1 A M I G M I D M + N = 1 B N I Q N I E N ]
+ k 1 A M = 1 A M I ( L I ( 1 ) L M I ( 2 ) ) D M ( μ M 2 + C I 2 ) + k 1 A N = 1 B N I ( L N I ( 3 ) L I ( 1 ) ) E N ( σ N 2 + K 2 C I 2 ) + B I C I ( 1 cosh C I ) } cosh ( C I y ¯ )
+ B I C I 2 [ 1 cosh ( C I y ¯ ) ] k 1 A C I L I ( 0 ) ( M = 1 A M I G M I D M + N = 1 B N I Q N I E N ) y ¯ sinh ( C I y ¯ ) 2
+ k 1 A C I M = 1 A M I D M ( μ M 2 + C I 2 ) { y ¯ sinh ( C I y ¯ ) 2 + C I [ cos ( μ M y ¯ ) cosh ( C I y ¯ ) ] C I 2 + μ M 2 }
+ k 1 A C I N = 1 B N I E N ( σ N 2 + K 2 C I 2 ) { C I [ cosh ( σ N 2 + K 2 y ¯ ) cosh ( C I y ¯ ) ] σ N 2 + K 2 C I 2 y ¯ sinh ( C I y ¯ ) 2 } ,
L I ( 0 ) = 1 C I cosh ( C I ) , L I ( 1 ) = 1 2 sinh ( C I ) , L M I ( 2 ) = C I [ cos ( μ M ) cosh ( C I ) ] C I 2 + μ M 2 , and
L N I ( 3 ) = C I [ cosh ( σ N 2 + K 2 ) cosh ( C I ) ] σ N 2 + K 2 C I 2 .
The detailed procedure is presented in Appendix A.
For non-Newtonian power-law nanofluid, the energy equation can be solved using the numerical algorithm developed later.

2.4. Entropy Generation Analysis

With the temperature distribution obtained, Nusselt number indicating heat transfer rate from the channel wall to the liquid inside the channel is deduced as:
N u = h a k e f f = k f k e f f T ¯ m ,
where the mean temperature is
T ¯ m = 0 1 0 α T ¯ w ¯ d x ¯ d y ¯ 0 1 0 α w ¯ d x ¯ d y ¯ = k f ( T m T w ) q s a
and the convective heat transfer coefficient is h = qs/(Tw − Tm).
The entropy generation rate per unit volume is presented, which evaluates the amount of useful energy destroyed during the process, namely, the thermodynamic irreversibility of EOF [35,43,44].
S L = S H + S J = k e f f T 2 ( 2 T x 2 + 2 T y 2 ) + σ E 2 T .
Since the heat generated by viscous dissipation is negligible [23,41], the irreversibility of the system equals the summation of local volumetric entropy generation rate due to heat transfer and Joule heating. The corresponding nondimensional form is
S ¯ L = S ¯ H + S ¯ J = k e f f k f ( T ¯ + Θ ) 2 [ ( T ¯ x ¯ ) 2 + ( T ¯ y ¯ ) 2 ] + S E 2 T ¯ + Θ ,
where Θ = T w k f / q w a [44]. The global total entropy generation rate can be obtained by integrating Equation (22) over the entire domain Ω . Then an important dimensionless number, Bejan number, representing the relative dominance of entropy generation due to heat transfer and Joule heating [45], is presented as
B e = Ω S H d Ω Ω S L d Ω .
It is noted that Be > 0.5 denotes that heat conduction irreversibility dominates and Be < 0.5 denotes the Joule heating irreversibility dominates.

3. Numerical Algorithm

Due to the nonlinearity of the governing equations above, the following numerical algorithm is developed to solve the velocity and temperature. Let t ¯ l = l Δ t ¯ , x ¯ i = i Δ x ¯ , y ¯ j = j Δ y ¯ , where Δ t ¯ is the temporal step, Δ x ¯ and Δ y ¯ are the spatial steps with l = 1, 2,…, L, and i = 1, 2, …, I, j = 1, 2,…, J. f l = [ f i , j l ] I × J represents the matrix obtained by discretizing a variable f over the points above. A numerical algorithm is proposed using the time splitting method in terms of time and the compact difference method in space. The complete P-B equation is numerically carried out based on the compact difference schemes, and the corresponding discretized form can be found in [20] and not presented here for conciseness.
To obtain the discretized velocity field w ¯ l = [ w ¯ i , j l ] I × J with w ¯ i , j l = w ¯ ( x ¯ i , y ¯ j , t ¯ l ) and the discretized temperature field T l = [ T ¯ i , j l ] I × J by using the iterative method, one can introduce / t ¯ on the left-hand side of Equations (11) and (16) and set a nonzero initial value for them,
w ¯ t ¯ = n η μ 0 ( 1 ϕ ) 2.5 ( W a ) n 1 ( | w ¯ x ¯ | n 1 2 w ¯ x ¯ 2 + | w ¯ y ¯ | n 1 2 w ¯ y ¯ 2 ) + K 2 sinh ψ ¯ ,
T ¯ t ¯ = 2 T ¯ x ¯ 2 + 2 T ¯ y ¯ 2 k f k e f f [ ( 1 + 1 α + S ) w ¯ w ¯ m S ] .
More specifically, each time step is split into two steps. In the first step, from t ¯ l to t ¯ * , one can focus on the third term on the right-hand side of Equations (19) and (20), and rewrite them as
w ¯ t ¯ = K 2 sinh ψ ¯ ,
T ¯ t ¯ = k f k e f f [ ( 1 + 1 α + S ) w ¯ w ¯ m S ] .
Accordingly, the velocity field at t ¯ = t ¯ * is given as w ¯ * = [ w ¯ i , j * ] I × J = [ K 2 sinh ψ i , j ] I × J . To obtain temperature field T * at t ¯ = t ¯ * , Equation (22) is discretized as
[ T ¯ t ¯ ] i , j l = c 1 c l + c 2 ,
where c 1 = k f k e f f ( 1 + 1 α + S ) , c 2 = k f k e f f S , and c l = [ w ¯ i , j l / ( w i , j l ) m ] I × J . The subscript m implies the mean value of w i , j l . Setting d l = c 1 c l + c 2 and the discretized temperature at t ¯ = t ¯ l as T l , the Runge–Kutta method is applied to Equation (23):
T ( 1 ) = α 1 T l + β 1 Δ t ¯ d l , T ( 2 ) = α 2 T l + β 2 [ T ( 1 ) + Δ t ¯ d l ] , T * = α 3 T l + β 3 [ T ( 2 ) + Δ t ¯ d l ] ,
where α1 = β1 = 1, α2 = 3/4, β2 = 1/4, α3 = 1/3, β2 = 2/3.
In the second step, from t ¯ * to t ¯ l + 1 , one can focus on the second-order derivatives in Equations (19) and (20), which are generally expressed as
f t ¯ = C 1 2 f x ¯ 2 + C 2 2 f y ¯ 2 .
To note, as solving w ¯ l i , j , C 1 = n η / [ μ 0 ( 1 ϕ ) 2.5 ] ( W / a ) n 1 | w ¯ / x ¯ | n 1 and C 2 = n η / [ μ 0 ( 1 ϕ ) 2.5 ] ( W / a ) n 1 | w ¯ / y ¯ | n 1 , the corresponding discretized forms have been given in Appendix B. When solving T ¯ l i , j , C1 = 1 and C2 = 1. According to the Crank–Nicolson technique, the discretized form of Equation (25) is
f l + 1 f * Δ t 1 2 { [ C 1 2 f x ¯ 2 ] l + 1 + [ C 1 2 f x ¯ 2 ] * } 1 2 { [ C 2 2 f y ¯ 2 ] l + 1 + [ C 2 2 f y ¯ 2 ] * } = O ( Δ t 2 )
Using the compact difference schemes and following the procedure presented in Appendix B, fi,jl+1 can be obtained as
A 1 P j 1 l + 1 + A 0 P j l + 1 + A 1 P j + 1 l + 1 = Β 1 P j 1 * + Β 0 P j * + Β 1 P j + 1 * , j = 2 , 3 , , J 1 ,
where P j l = ( f 1 , j l , f 2 , j l , , f I , j l ) T , A 0 = t r i d i a g ( a 1 2 a 3 , 1 2 a 1 2 a 2 + 4 a 3 , a 1 2 a 3 ) , A 1 = t r i d i a g ( a 3 , a 2 2 a 3 , a 3 ) , Β 0 = t r i d i a g ( b 1 2 b 3 , 1 2 b 1 2 b 2 + 4 b 3 , b 1 2 b 3 ) , Β 1 = t r i d i a g ( b 3 , b 2 2 b 3 , b 3 ) , and h x ¯ = Δ t ¯ Δ x ¯ 2 , h y ¯ = Δ t ¯ Δ y ¯ 2 , a 1 = 1 12 C 1 h x ¯ 2 , a 2 = 1 12 C 2 h y ¯ 2 , a 3 = 1 144 C 1 h x ¯ 24 C 2 h y ¯ 24 , b 1 = 1 12 + C 1 h x ¯ 2 , b 2 = 1 12 + C 2 h y ¯ 2 , b 3 = 1 144 + C 1 h x ¯ 24 + C 2 h y ¯ 24 .
As a result, replacing f with w ¯ and T ¯ , the discretized solutions w ¯ l i , j and T ¯ l i , j can be computed. Finally, a specified criterion Δ L is given to identify that if the velocity is fully developed, i.e., w ¯ l w ¯ l + 1 < Δ L . Then, the fully developed velocity w ¯ i , j and temperature T ¯ i , j are obtained based on the iterative method above.

4. Results

A parametric study for velocity, temperature, and Nusselt number is conducted in the case of different flow behavior index n, electrokinetic width K, volume fraction of nanoparticles ϕ, and Joule heating parameter S. The nanoparticle is regarded as aluminum oxide [38]. The other physical parameters of the base fluid come from [9] and [20]. The typical values are presented in Table 1 below.
For validating the numerical algorithm proposed above, a test of grid dependence is conducted. As a result, the volumetric domain Ω is discretized to a grid system of 121 × 121. The numeric algorithm shows that the numerical solutions can achieve the accuracy of O ( Δ t ¯ 2 + Δ x ¯ 4 + Δ y ¯ 4 ) , which is demonstrated in numerical computation by obtaining the numerical error of the order 10−8 between the numerical velocity and analytical velocity of Newtonian nanofluid EOF. Furthermore, the iterative criterion is given as ΔL < 10−10 to make sure the EOF is fully developed. For the EOF of Newtonian nanofluid, in Figure 2a, the numerical velocity profile at y ¯ = 0 is compared with the analytical velocity obtained from Equation (13) and the numerical temperature profile at y ¯ = 0 is compared with the analytical temperature obtained from Equation (19) when K = 10, ϕ = 0.06, and S = 3. To clearly show the comparison, 25 grid points of the numerical solution are plotted in Figure 2. It shows that the numerical solutions show good agreement with the analytical solutions and thus, the numerical algorithm is valid, which can be used to numerically solve velocity, temperature, Nusselt number, and Bejan number of power-law nanofluid EOF.
Figure 3 illustrates velocity distributions of shear thickening nanofluid, Newtonian nanofluid, and shear thinning nanofluid at different volume fractions when K = 10. Figure 3a,c,e show the velocity distributions of shear thinning fluid, Newtonian fluid, and shear thickening fluid without consideration of nanoparticles in fluid, namely at ϕ = 0. Figure 3b,d,f show the velocity distributions of shear thinning nanofluid, Newtonian nanofluid, and shear thickening nanofluid at ϕ = 0.06. It is noted from these two groups of figures that the velocity shows homogeneous abatement across the channel as the volume fraction of nanoparticles increases no matter what fluid rheology is considered. This is because the nanoparticles in fluid enhance the effective viscosity of power-law nanofluid in response to the sear rate and cause greater dispersion of the velocity distribution. As observed in [7], although considering the EOF of power-law nanofluid, the shear thinning feature of base fluid leads to greater velocity compared to the shear thickening feature of the base fluid.
Without consideration of nanoparticles in fluid (ϕ = 0) and when K = 10, the temperature distributions of shear thinning fluid, Newtonian fluid, and shear thickening fluid at S = 0 are displayed in Figure 4a,c,e. For the same prescribed conditions, the temperature distributions at S = 3 are displayed in Figure 4b,d,f. The comparison between these two groups of figures shows that the variation of temperature along the channel becomes greater with the increase of the Joule heating parameter S irrespective of the fluid type. Furthermore, the dimensionless temperature increases with the flow behavior index n, indicating that the difference of heat transfer characters caused by different shear rates of fluid type is remarkable, especially for shear thinning base fluid.
To observe the influence of volume fraction of nanoparticles on the heat transfer characteristics, the temperature distributions are evaluated in Figure 5 when the nanoparticle volume fraction ϕ is increased to 0.06, and other parameters keep unchanged. The corresponding temperature profiles at y ¯ = 0 is illustrated in Figure 6 for different volume fractions of the nanoparticles. From Figure 6 and the comparison between Figure 4 and Figure 5, the dimensionless temperature distribution shows an increment when the volume fraction of nanoparticles ϕ increases, or the flow behavior index n increases. It means that the combined effects of larger thermal conductivity of nanoparticles and the change in shear stress leads to the increase in thermal diffusion effect and in further the temperature.
To observe the influence of the Joule heating parameter on the temperature, in Figure 7, the temperature profiles at y ¯ = 0 for different Joule heating parameters S are illustrated in the case of shear thinning nanofluid, Newtonian nanofluid, and shear thickening nanofluid when ϕ = 0.06 and K = 10. As seen in Figure 5, no matter what nanofluid is considered, as the Joule heating parameter increases the variation of temperature increases, namely, the temperature difference between the centerline and the channel wall enlarges. It implies that although the energy from Joule heating effect is uniformly enhanced, the heat transferred through convection near the channel wall is much less than near the centerline. And this effect is more notable when considering shear thinning base fluid.
The variations of the Nusselt number with the electrokinetic width are shown in Figure 8 for different fluid types and Joule heating parameters when ϕ = 0.06. The cases S < 0 and S > 0 correspond to the surface cooling effect and surface heating effect, respectively. Importantly, in the case of S < 0, the Nusselt number decreases as electrokinetic width K increases, which increases when S ≥ 0 and other parameters remain unchanged. Surface heating (cooling) has a uniform influence on the temperature distribution, which leads to the increasing (decreasing) tendency of Nusselt number with electrokinetic width. Furthermore, the decrease in flow behavior index results in the pronounced abatement of Nusselt number, meaning that the heat transfer rate of shear thinning nanofluid is much smaller than the heat transfer rate of Newtonian fluid and shear thickening fluid. This is consistent with the prediction in [22].
Figure 9 depicts the variations of the Nusselt number with the Joule heating parameter for different fluid types and different nanoparticle volume fraction when K = 10. It shows that the Nusselt number is inversely proportional to the Joule heating parameter, and the decreasing rate becomes smaller as nanofluid changes from shear thickening to shear thinning. Since it is considered that Joule heating has a spatially uniform effect, the greater Joule heating parameter corresponds to relatively uniform heat distribution. Consequently, the convective heat transfer rate is reduced as presented in Figure 7, and one can observe the decreasing trend of the Nusselt number, namely, heat transfer rate. Being consistent with the prediction in Figure 7, the smaller convective heat transfer rate of shear thinning nanofluid leads to a smaller Nusselt number. Furthermore, the comparison between Figure 9a,b show that no matter what base fluid and what Joule heating parameter is taken, the increase in the volume fraction of nanoparticles ϕ significantly improves the heat transfer rate.
The variation of the Nusselt number Nu and the variation of Bejan number Be with flow behavior index n for different volume fractions of nanoparticles are depicted in Figure 10 when K = 20 and S = 3. The combined effects of nonlinear rheological behavior and nanoparticle volume fraction on Nusselt number are investigated. To note, the dependence of Nusselt number on the flow behavior index shows nonlinear. For shear thinning base fluid, the Nusselt number is more amenable to the change in flow behavior index n. The increase in viscosity for shear thinning fluid resulting from the decrease in flow behavior index reduces the Nusselt number, and thus, the heat transfer rate of shear thinning nanofluid is obviously less than that of shear thickening nanofluid. Furthermore, the variation of the Nusselt number with the volume fraction of nanoparticles is almost linear. It shows that the heat transfer performance of nanofluid is clearly enhanced compared to the conventional fluid regardless of the base fluid rheology. These results are also witnessed in Figure 10b. The increase in nanoparticle volume fraction and flow behavior index leads to the enhancement of heat transfer, which further reduces the “loss” of useful energy, i.e., the total irreversibility of the system. Moreover, the Be > 0.5 implies that the primary irreversibility is contributed to the heat conduction when S = 3.

5. Conclusions

In the present paper, a numerical investigation for hydrodynamically and thermally fully developed EOF of power-law nanofluid in a rectangular microchannel has been conducted. The interactive influences of electrokinetic width K, flow behavior index n, nanoparticle volume fraction ϕ, and Joule heating parameter S on the velocity field, temperature field, Nusselt number, and Bejan number have been studied. The following conclusions can be drawn:
  • The velocity is reduced due to the increase in viscosity of nanofluid resulting from the addition of nanoparticles.
  • The volume fraction of nanoparticles ϕ truly enhances the dimensionless temperature and Nusselt number by changing the thermal conductivity and shear stress of velocity, no matter what type of base fluid is considered. Accordingly, the Bejan number decreases with increasing nanoparticle volume fraction and flow behavior index.
  • Nusselt number is an increasing function of the electrokinetic width K for S ≥ 0 and a decreasing function for S < 0.
  • The Nusselt number shows abatement, and the temperature exhibits greater variation along the channel as Joule heating parameter S increases. The heat transfer rate of shear thickening nanofluid is much larger than that of shear thinning fluid.
To sum up, nanoparticle volume fraction and Joule heating effect show significant influence on the thermal behavior of EOF of power-law nanofluid. It should be noted that the dependence of the Nusselt number on the flow behavior index is nonlinear. The Nusselt number becomes a weak function of electrokinetic width as the electrokinetic width gets larger. The heat transfer rate of shear thinning nanofluid is more amenable to the change in nanofluid rheology. Therefore, the heat transfer of power-law nanofluid has to be carefully treated, and one should make a distinction when manipulating different types of nanofluid in practical problems. This understanding of the mechanism of thermally fully developed EOF can be referred to in more efficient control of fluids and optimization of microfluidic devices.

Funding

This work was supported by the Talent Introduction Foundation of Guangdong University of Petrochemical Technology, grant number 2018rc16, the Scientific Research Foundation of Universities in Guangdong Province for Young Talents, grant number 2018KQNCX165, and the Scientific Research Foundation of Guangdong University of Petrochemical Technology, grant number 513040.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

The velocity distribution of the Newtonian nanofluid is solved following the procedure in [40] and is presented as Equation (13). As the nanoparticle volume fraction vanishes, the velocity of Equation (13) becomes Equation (14) in [21].
The temperature distribution of Newtonian nanofluid EOF is solved as follows: in the first step, based on the method of variable separation [40], the temperature can be expressed as
T ¯ ( x ¯ , y ¯ ) = I = 1 cos ( β I x ¯ ) Y I ( y ¯ )
with β I = ( 2 I 1 ) π / ( 2 α ) , I = 1,2,3,….
In the second step, Y I ( y ¯ ) needs to be determined. Substituting Equation (A1) into Equation (17) yields
I = 1 cos ( β I x ¯ ) χ I ( y ¯ ) = k 1 w ¯ k 2 S ,
[ d 2 Y I ( y ¯ ) d y ¯ 2 β I Y I ( y ¯ ) ] = χ I ( y ¯ ) .
χ I ( y ¯ ) can be expressed from Equation (A1a) using the orthogonality of Fourier series, namely,
χ I ( y ¯ ) = 2 α 0 α [ k 1 w ¯ k 2 S ] cos ( β I x ¯ ) d x ¯ ,
where k 1 = k f k e f f ( 1 + 1 α + S ) and k 2 = k f k e f f . Replacing w ¯ with Equation (13), further rearrangements produce χ I ( y ¯ ) = k 1 [ U I ( y ¯ ) ] B I , where B I = 2 k 2 S ( 1 ) I + 1 / ( α β I ) .
In the third step, with the expression of χ I ( y ¯ ) , the ODE of Equation (A1b) is solved adopting the method of constant variation, the solution form is given as
Y I ( y ¯ ) = Y I G ( y ¯ ) + Y I P ( y ¯ ) ,
Y I G ( y ¯ ) = H 1 e C I y ¯ + H 2 e C I y ¯ ,
Y I P ( y ¯ ) = 1 C I 0 y ¯ χ I ( s ) sinh [ C I ( y ¯ s ) ] d s .
Here Equation (A3a) indicates the general solution of the homogeneous version of Equation (A1b). Equation (A3b) indicates the particular solution of Equation (A1b), which is determined by the method of Laplace transformation. Combining boundary condition, the coefficients are deduced as
H 1 = H 2 = 1 2 cosh ( C I ) 1 C I 0 1 χ I ( s ) sinh [ C I ( 1 s ) ] d s .
As a result, Equation (A3) becomes
Y I ( y ¯ ) = cosh ( C I y ¯ ) cosh ( C I ) 1 C I [ 0 1 k 1 [ U I ( s ) ] sinh [ C I ( 1 s ) ] d s + B I C I ( 1 cosh C I ) ] , + B I C I 2 [ 1 cosh ( C I y ¯ ) ] + k 1 C I 0 y ¯ ( U I ( s ) ) sinh [ C I ( y ¯ s ) ] d s .
Substituting the expression of U I ( y ) into Equation (A5) yields
Y I ( y ¯ ) = L I ( 0 ) { k 1 A L I ( 1 ) L I ( 0 ) ( M = 1 A M I G M I D M + N = 1 B N I Q N I E N ) + k 1 A M = 1 A M I ( L I ( 1 ) L M I ( 2 ) ) D M ( μ M 2 + C I 2 ) + k 1 A N = 1 B N I ( L N I ( 3 ) L I ( 1 ) ) E N ( σ N 2 + K 2 C I 2 ) + B I C I ( 1 cosh C I ) } cosh ( C I y ¯ ) + B I C I 2 [ 1 cosh ( C I y ¯ ) ] k 1 A C I L I ( 0 ) ( M = 1 A M I G M I D M + N = 1 B N I Q N I E N ) y ¯ sinh ( C I y ¯ ) 2 + k 1 A C I M = 1 A M I D M ( μ M 2 + C I 2 ) { y ¯ sinh ( C I y ¯ ) 2 + C I [ cos ( μ M y ¯ ) cosh ( C I y ¯ ) ] C I 2 + μ M 2 } + k 1 A C I N = 1 B N I E N ( σ N 2 + K 2 C I 2 ) { C I [ cosh ( σ N 2 + K 2 y ¯ ) cosh ( C I y ¯ ) ] σ N 2 + K 2 C I 2 y ¯ sinh ( C I y ¯ ) 2 }
where L I ( 0 ) = 1 C I cosh ( C I ) , L I ( 1 ) = 0 1 cosh ( C I s ) sinh [ C I ( 1 s ) ] d s = 1 2 sinh ( C I ) ,
L M I ( 2 ) = 0 1 cos ( μ M y ) sinh [ C I ( 1 s ) ] d s = C I [ cos ( μ M ) cosh ( C I ) ] C I 2 + μ M 2 ,
L N I ( 3 ) = 0 1 cosh ( σ N 2 + K 2 s ) sinh [ C I ( 1 s ) ] d s = C I [ cosh ( σ N 2 + K 2 ) cosh ( C I ) ] σ N 2 + K 2 C I 2 .
Consequently, the analytical solution of temperature for Newtonian nanofluid flow is carried out.

Appendix B

In the case of f = w ¯ , the compact difference schemes for computing C1 and C2 are given as
1 6 F x ¯ , i + 1 , j ( w ¯ ) + 2 3 F x ¯ , i , j ( w ¯ ) + 1 6 F x ¯ , i 1 , j ( w ¯ ) = ( w ¯ i + 1 , j w ¯ i 1 , j ) / 2 Δ x ¯ ,
1 6 F y ¯ , i , j + 1 ( w ¯ ) + 2 3 F y ¯ , i , j ( w ¯ ) + 1 6 F y ¯ , i , j 1 ( w ¯ ) = ( w ¯ i , j + 1 w ¯ i , j 1 ) / 2 Δ y ¯ ,
in which F x ¯ , i , j ( w ¯ ) and F y ¯ , i , j ( w ¯ ) are the difference approximations of w ¯ / x ¯ and w ¯ / y ¯ , respectively.
To derive the further discretized form of Equation (32), setting δ x ¯ 2 f i , j = ( f i + 1 , j 2 f i , j + f i 1 , j ) / Δ x ¯ 2 , δ y ¯ 2 f i , j = ( f i , j + 1 2 f i , j + f i , j 1 ) / Δ y ¯ 2 , the following finite difference operators are introduced:
δ ¯ x ¯ 2 f i , j = ( 1 + Δ x ¯ 2 12 δ x ¯ 2 ) f i , j ,   δ ¯ y ¯ 2 f i , j = ( 1 + Δ y ¯ 2 12 δ y ¯ 2 ) f i , j .
Note that,
δ x ¯ 2 f i , j = δ ¯ x ¯ 2 [ 2 f x ¯ 2 ] i , j + O ( Δ x ¯ 4 ) , δ y ¯ 2 f i , j = δ ¯ y ¯ 2 [ 2 f y ¯ 2 ] i , j + O ( Δ y ¯ 4 )
or symbolically,
δ ¯ x ¯ 2 δ x ¯ 2 f i , j = [ 2 f x ¯ 2 ] i , j + O ( Δ x ¯ 4 ) ,   δ ¯ y ¯ 2 δ y ¯ 2 f i , j = [ 2 f x ¯ 2 ] i , j + O ( Δ y ¯ 4 ) ,
where δ ¯ x ¯ 2 = ( δ ¯ x ¯ 2 ) 1 and δ ¯ y ¯ 2 = ( δ ¯ y ¯ 2 ) 1 are the inverse of δ ¯ x ¯ 2 and δ ¯ y ¯ 2 , respectively.
The fourth-order compact approximations (B5) is applied to the derivatives in Equation (32), and it yields:
f l + 1 f * Δ t ¯ C 1 2 ( δ ¯ x ¯ 2 δ ¯ x ¯ 2 f i , j l + 1 + δ ¯ x ¯ 2 δ ¯ x ¯ 2 f i , j * ) C 2 2 ( δ ¯ y ¯ 2 δ ¯ y ¯ 2 f i , j l + 1 + δ ¯ y ¯ 2 δ ¯ y ¯ 2 f i , j * ) = O ( Δ t ¯ 2 + Δ x ¯ 4 + Δ y ¯ 4 ) .
Multiplying the above equation by the finite difference operator δ ¯ x ¯ 2 δ ¯ y ¯ 2 , one has
( δ ¯ x ¯ 2 δ ¯ y ¯ 2 Δ t ¯ C 1 2 δ ¯ y ¯ 2 δ x ¯ 2 Δ t ¯ C 1 2 δ ¯ x ¯ 2 δ y ¯ 2 ) f i , j l + 1 = ( δ ¯ x ¯ 2 δ ¯ y ¯ 2 + Δ t ¯ C 1 2 δ ¯ y ¯ 2 δ x ¯ 2 + Δ t ¯ C 1 2 δ ¯ x ¯ 2 δ y ¯ 2 ) f i , j * + O ( Δ t ¯ 3 + Δ t ¯ Δ x ¯ 4 + Δ t ¯ Δ y ¯ 4 ) .
Rearranging Equation (B7) and dropping the term O ( Δ t ¯ 3 + Δ t ¯ Δ x ¯ 4 + Δ t ¯ Δ y ¯ 4 ) , the left-hand side and the right-hand side can be approximated, respectively, as
( δ ¯ x ¯ 2 δ ¯ y ¯ 2 Δ t ¯ C 1 2 δ ¯ y ¯ 2 δ x ¯ 2 Δ t ¯ C 1 2 δ ¯ x ¯ 2 δ y ¯ 2 ) f i , j l + 1 = ( 1 2 a 1 2 a 2 + 4 a 3 ) f i , j l + 1 + ( a 1 2 a 3 ) ( f i + 1 , j l + 1 + f i 1 , j l + 1 ) + ( a 2 2 a 3 ) ( f i , j + 1 l + 1 + f i , j 1 l + 1 ) + a 3 ( f i + 1 , j + 1 l + 1 + f i + 1 , j 1 l + 1 + f i 1 , j + 1 l + 1 + f i 1 , j 1 l + 1 )
( δ ¯ x ¯ 2 δ ¯ y ¯ 2 + Δ t ¯ C 1 2 δ ¯ y ¯ 2 δ x ¯ 2 + Δ t ¯ C 1 2 δ ¯ x ¯ 2 δ y ¯ 2 ) f i , j * = ( 1 2 b 1 2 b 2 + 4 b 3 ) f i , j * + ( b 1 2 b 3 ) ( f i + 1 , j * + f i 1 , j * ) + ( b 2 2 b 3 ) ( f i , j + 1 * + f i , j 1 * ) + b 3 ( f i + 1 , j + 1 * + f i + 1 , j 1 * + f i 1 , j + 1 * + f i 1 , j 1 * )
with h x ¯ = Δ t ¯ / Δ x ¯ 2 , h y ¯ = Δ t ¯ / Δ y ¯ 2 , a 1 = 1 / 12 C 1 h x ¯ / 2 , a 2 = 1 / 12 C 2 h y ¯ / 2 , b 1 = 1 / 12 + C 1 h x ¯ / 2 , b 2 = 1 / 12 + C 2 h y ¯ / 2 , a 3 = 1 / 144 C 1 h x ¯ / 24 C 2 h y ¯ / 24 , b 3 = 1 / 144 + C 1 h x ¯ / 24 + C 2 h y ¯ / 24 .
The following second-order one-side difference approximations are applied to the boundary conditions expressed as Equations (12) and (18):
[ 2 f x ¯ 2 ] 1 , j = ( 2 f 1 , j 5 f 2 , j + 4 f 3 , j f 4 , j ) / Δ x ¯ 2 ,   [ 2 f x ¯ 2 ] I , j = ( 2 f I , j 5 f I 1 , j + 4 f I 2 , j f I 3 , j ) / Δ x ¯ 2 , [ 2 f y ¯ 2 ] i , 1 = ( 2 f i , 1 5 f i , 2 + 4 f i , 3 f i , 4 ) / Δ y ¯ 2 ,   [ 2 f y ¯ 2 ] i , J = ( 2 f i , J 5 f i , J 1 + 4 f i , J 2 f i , J 3 ) / Δ y ¯ 2 .
Consequently, the discretized form of f can be expressed as a column vector P j l = ( f 1 , j l , f 2 , j l , , f I , j l ) T , and the corresponding coefficients can be presented as the symmetric tridiagonal matrices given in Equation (33).

References

  1. Stone, H.A.; Stroock, A.D.; Ajdari, A. Engineering flows in small devices. Ann. Rev. Fluid Mech. 2004, 36, 381–411. [Google Scholar] [CrossRef]
  2. Bruus, H. Theoretical Microfluidics; Oxford University Press: New York, NY, USA, 2008. [Google Scholar]
  3. Zeng, S.; Chen, C.H.; Mikkelsen, J.C., Jr.; Santiago, J.G. Fabrication and characterization of electroosmotic micropumps. Sens. Actuators B 2001, 79, 107–114. [Google Scholar] [CrossRef]
  4. Horiuchi, K.; Dutta, P. Joule heating effects in electroosmotically driven microchannel flows. Int. J. Heat Mass Transf. 2004, 47, 3085–3095. [Google Scholar] [CrossRef]
  5. Das, S.; Chakraborty, S. Analytical solutions for velocity, temperature and concentration distribution in electroosmotic microchannel flows of a non-newtonian bio-fluid. Anal. Chim. Acta 2006, 559, 15–24. [Google Scholar] [CrossRef]
  6. Datta, S.; Ghosal, S.; Patankar, N.A. Electroosmotic flow in a rectangular channel with variable wall zeta-potential: Comparison of numerical simulation with asymptotic theory. Electrophoresis 2006, 27, 611–619. [Google Scholar] [CrossRef]
  7. Deng, S.Y.; Jian, Y.J.; Bi, Y.H.; Chang, L.; Wang, H.J.; Liu, Q.S. Unsteady electroosmotic flow of power-law fluid in a rectangular microchannel. Mech. Res. Commun. 2012, 39, 9–14. [Google Scholar] [CrossRef]
  8. Srinivas, B. Electroosmotic flow of a power law fluid in an elliptic microchannel. Colloids Surf. A 2016, 492, 144–151. [Google Scholar] [CrossRef]
  9. Zhao, C.L.; Yang, C. Joule heating induced heat transfer for electroosmotic flow of power-law fluids in a microcapillary. Int. J. Heat Mass Transf. 2012, 55, 2044–2051. [Google Scholar] [CrossRef]
  10. Shit, G.C.; Mondal, A.; Sinha, A.; Kundu, P.K. Electro-osmotic flow of power-law fluid and heat transfer in a micro-channel with effects of joule heating and thermal radiation. Physica A 2016, 462, 1040–1057. [Google Scholar] [CrossRef]
  11. Kwak, H.S.; Kim, H.; Hyun, J.M.; Song, T.H. Thermal control of electroosmotic flow in a microchannel through temperature-dependent properties. J. Colloid Interface Sci. 2009, 335, 123–129. [Google Scholar] [CrossRef]
  12. Afonso, A.M.; Alves, M.A.; Pinho, F.T. Analytical solution of two-fluid electro-osmotic flows of viscoelastic fluids. J. Colloid Interface Sci. 2013, 395, 277–286. [Google Scholar] [CrossRef]
  13. Si, D.Q.; Jian, Y.J.; Chang, L.; Liu, Q.S. Unsteady rotating electroosmotic flow through a slit microchannel. J. Mech. 2016, 32, 603–611. [Google Scholar] [CrossRef]
  14. Qi, C.; Ng, C.O. Rotating electroosmotic flow of viscoplastic material between two parallel plates. Colloids Surf. A 2017, 513, 355–366. [Google Scholar] [CrossRef] [Green Version]
  15. Vakili, M.A.; Sadeghi, A.; Saidi, M.H. Pressure effects on electroosmotic flow of power-law fluids in rectangular microchannels. Theor. Comp. Fluid Dyn. 2014, 28, 409–426. [Google Scholar] [CrossRef]
  16. Li, F.Q.; Jian, Y.J.; Xie, Z.Y.; Liu, Y.B.; Liu, Q.S. Transient alternating current electroosmotic flow of a jeffrey fluid through a polyelectrolyte-grafted nanochannel. Res. Adv. 2017, 7, 782–790. [Google Scholar] [CrossRef]
  17. Li, H.C.; Jian, Y.J. Dispersion for periodic electro-osmotic flow of maxwell fluid through a microtube. Int. J. Heat Mass Transf. 2017, 115, 703–713. [Google Scholar] [CrossRef]
  18. Zhao, C.L.; Zhang, W.Y.; Yang, C. Dynamic electroosmotic flows of power-law fluids in rectangular microchannels. Micromachines 2017, 8, 34. [Google Scholar] [CrossRef]
  19. Deng, S.Y. The parametric study of electroosmotically driven flow of power-law fluid in a cylindrical microcapillary at high zeta potential. Micromachines 2017, 8, 344. [Google Scholar] [CrossRef]
  20. Zhu, Q.Y.; Deng, S.Y.; Chen, Y.Q. Periodical pressure-driven electrokinetic flow of power-law fluids through a rectangular microchannel. J. Non-Newton. Fluid Mech. 2014, 203, 38–50. [Google Scholar] [CrossRef]
  21. Su, J.; Jian, Y.J.; Chang, L. Thermally fully developed electroosmotic flow through a rectangular microchannel. Int. J. Heat Mass Transf. 2012, 55, 6285–6290. [Google Scholar] [CrossRef]
  22. Sun, Y.J.; Jian, Y.J.; Chang, L.; Liu, Q.S. Thermally fully developed electroosmotic flow of power-law fluids in a circular microchannel. J. Mech. 2013, 29, 609–616. [Google Scholar] [CrossRef]
  23. Escandón, J.; Bautista, O.; Méndez, F. Entropy generation in purely electroosmotic flows of non-newtonian fluids in a microchannel. Energy 2013, 55, 486–496. [Google Scholar] [CrossRef]
  24. Zade, A.Q.; Manzari, M.T.; Hannani, S.K. An analytical solution for thermally fully developed combined pressure—Electroosmotically driven flow in microchannels. Int. J. Heat Mass Transf. 2007, 50, 1087–1096. [Google Scholar] [CrossRef]
  25. Babaie, A.; Saidi, M.H.; Sadeghi, A. Heat transfer characteristics of mixed electroosmotic and pressure driven flow of power-law fluids in a slit microchannel. Int. J. Therm. Sci. 2012, 53, 71–79. [Google Scholar] [CrossRef]
  26. Vakili, M.A.; Saidi, M.H.; Sadeghi, A. Thermal transport characteristics pertinent to electrokinetic flow of power-law fluids in rectangular microchannels. Int. J. Therm. Sci. 2014, 79, 76–89. [Google Scholar] [CrossRef]
  27. Xie, Z.Y.; Jian, Y.J.; Li, F.Q. Thermal transport of magnetohydrodynamic electroosmotic flow in circular cylindrical microchannels. Int. J. Heat Mass Transf. 2018, 119, 355–364. [Google Scholar] [CrossRef]
  28. Choi, U.S. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab.: Lemont, IL, USA, 1995; Volume 231, pp. 99–103.
  29. Fan, J.; Wang, L. Review of heat conduction in nanofluids. J. Heat Transf. 2011, 133, 40801–40814. [Google Scholar] [CrossRef]
  30. Avsec, J.; Oblak, M. The calculation of thermal conductivity, viscosity and thermodynamic properties for nanofluids on the basis of statistical nanomechanics. Int. J. Heat Mass Transf. 2007, 50, 4331–4341. [Google Scholar] [CrossRef]
  31. Saidur, R.; Leong, K.Y.; Mohammad, H.A. A review on applications and challenges of nanofluids. Renew. Sustain. Energy Rev. 2011, 15, 1646–1668. [Google Scholar] [CrossRef]
  32. Hemmat Esfe, M.; Saedodin, S.; Mahian, O.; Wongwises, S. Thermal conductivity of AL2O3/water nanofluids. J. Therm. Anal. Calorim. 2014, 117, 675–681. [Google Scholar] [CrossRef]
  33. Yu, W.; France, D.M.; Timofeeva, E.V.; Singh, D.; Routbort, J.L. Comparative review of turbulent heat transfer of nanofluids. Int. J. Heat Mass Transf. 2012, 55, 5380–5396. [Google Scholar] [CrossRef]
  34. Zhao, G.P.; Jian, Y.J.; Li, F.Q. Heat transfer of nanofluids in microtubes under the effects of streaming potential. Appl. Therm. Eng. 2016, 100, 1299–1307. [Google Scholar] [CrossRef]
  35. Ganguly, S.; Sarkar, S.; Kumar Hota, T.; Mishra, M. Thermally developing combined electroosmotic and pressure-driven flow of nanofluids in a microchannel under the effect of magnetic field. Chem. Eng. Sci. 2015, 126, 10–21. [Google Scholar] [CrossRef]
  36. Shehzad, N.; Zeeshan, A.; Ellahi, R. Electroosmotic flow of mhd power law AL2O3-PVC nanouid in a horizontal channel: Couette-poiseuille flow model. Commun. Theor. Phys. 2018, 69, 655. [Google Scholar] [CrossRef]
  37. Ellahi, R.; Hassan, M.; Zeeshan, A. A study of heat transfer in power law nanofluid. Therm. Sci. 2016, 20, 2015–2026. [Google Scholar] [CrossRef] [Green Version]
  38. Si, X.; Li, H.; Zheng, L.; Shen, Y.; Zhang, X. A mixed convection flow and heat transfer of pseudo-plastic power law nanofluids past a stretching vertical plate. Int. J. Heat Mass Transf. 2017, 105, 350–358. [Google Scholar] [CrossRef]
  39. Brinkman, H.C. The viscosity of concentrated suspensions and solutions. J. Chem. Phys. 1952, 20, 571. [Google Scholar] [CrossRef]
  40. Jian, Y.J.; Liu, Q.S.; Yang, L.G. Ac electroosmotic flow of generalized maxwell fluids in a rectangular microchannel. J. Non-Newton. Fluid Mech. 2011, 166, 1304–1314. [Google Scholar] [CrossRef]
  41. Gleeson, J.P. Electroosmotic flows with random zeta potential. J. Colloid Interface Sci. 2002, 249, 217–226. [Google Scholar] [CrossRef]
  42. Yu, W.; Choi, S.U.S. The role of interfacial layer in the enhanced thermal conductivity of nanofluids a renovated maxwell model. J. Nanopart. Res. 2003, 6, 355–361. [Google Scholar] [CrossRef]
  43. Jian, Y.J. Transient mhd heat transfer and entropy generation in a microparallel channel combined with pressure and electroosmotic effects. Int. J. Heat Mass Transf. 2015, 89, 193–205. [Google Scholar] [CrossRef]
  44. Zhao, G.; Jian, Y.; Li, F. Streaming potential and heat transfer of nanofluids in parallel plate microchannels. Colloids Surf. A 2016, 498, 239–247. [Google Scholar] [CrossRef]
  45. Bejan, A. The thermodynamic design of heat and mass transfer processes and devices. Heat Fluid Flow 1987, 8, 258–276. [Google Scholar] [CrossRef]
Figure 1. Sketch of the rectangular microchannel with height 2a, width 2b, and constant heat flux qs at the walls.
Figure 1. Sketch of the rectangular microchannel with height 2a, width 2b, and constant heat flux qs at the walls.
Micromachines 10 00363 g001
Figure 2. The comparison between the numerical solution and analytical solution from for (a) velocity distribution at y ¯ = 0 and (b) temperature distribution at y ¯ = 0 .
Figure 2. The comparison between the numerical solution and analytical solution from for (a) velocity distribution at y ¯ = 0 and (b) temperature distribution at y ¯ = 0 .
Micromachines 10 00363 g002
Figure 3. The comparison among the velocity distributions of shear thinning nanofluid, Newtonian nanofluid and shear thickening nanofluid at different nanoparticle volume fraction: (a) ϕ = 0, n = 0.8; (b) ϕ = 0.06, n = 0.8; (c) ϕ = 0, n = 1.0; (d) ϕ = 0.06, n = 1.0; (e) ϕ = 0, n = 1.2; (f) ϕ = 0.06, n = 1.2.
Figure 3. The comparison among the velocity distributions of shear thinning nanofluid, Newtonian nanofluid and shear thickening nanofluid at different nanoparticle volume fraction: (a) ϕ = 0, n = 0.8; (b) ϕ = 0.06, n = 0.8; (c) ϕ = 0, n = 1.0; (d) ϕ = 0.06, n = 1.0; (e) ϕ = 0, n = 1.2; (f) ϕ = 0.06, n = 1.2.
Micromachines 10 00363 g003
Figure 4. The dimensionless temperature distributions of shear thinning fluid, Newtonian fluid, and shear thickening fluid for different Joule heating parameter when nanoparticle volume fraction is ϕ = 0: (a) n = 0.8, S = 0; (b) n = 0.8, S = 3; (c) n = 1.0, S = 0; (d) n = 1.0, S = 3; (e) n = 1.2, S = 0; (f) n = 1.2, S = 3.
Figure 4. The dimensionless temperature distributions of shear thinning fluid, Newtonian fluid, and shear thickening fluid for different Joule heating parameter when nanoparticle volume fraction is ϕ = 0: (a) n = 0.8, S = 0; (b) n = 0.8, S = 3; (c) n = 1.0, S = 0; (d) n = 1.0, S = 3; (e) n = 1.2, S = 0; (f) n = 1.2, S = 3.
Micromachines 10 00363 g004
Figure 5. The dimensionless temperature distributions of shear thinning nanofluid, Newtonian nanofluid, and shear thickening nanofluid for different Joule heating parameters S when nanoparticle volume fraction is ϕ = 0.06: (a) n = 0.8, S = 0; (b) n = 0.8, S = 3; (c) n = 1.0, S = 0; (d) n = 1.0, S = 3; (e) n = 1.2, S = 0; (f) n = 1.2, S = 3.
Figure 5. The dimensionless temperature distributions of shear thinning nanofluid, Newtonian nanofluid, and shear thickening nanofluid for different Joule heating parameters S when nanoparticle volume fraction is ϕ = 0.06: (a) n = 0.8, S = 0; (b) n = 0.8, S = 3; (c) n = 1.0, S = 0; (d) n = 1.0, S = 3; (e) n = 1.2, S = 0; (f) n = 1.2, S = 3.
Micromachines 10 00363 g005
Figure 6. The comparison among temperature profiles at y ¯ = 0 for different nanoparticle volume fractions ϕ in the case of (a) n = 0.8, (b) n = 1.0 and (c) n = 1.2.
Figure 6. The comparison among temperature profiles at y ¯ = 0 for different nanoparticle volume fractions ϕ in the case of (a) n = 0.8, (b) n = 1.0 and (c) n = 1.2.
Micromachines 10 00363 g006
Figure 7. The comparison among temperature profiles at y ¯ = 0 for different Joule heating parameters in the case of (a) n = 0.8, (b) n = 1.0 and (c) n = 1.2.
Figure 7. The comparison among temperature profiles at y ¯ = 0 for different Joule heating parameters in the case of (a) n = 0.8, (b) n = 1.0 and (c) n = 1.2.
Micromachines 10 00363 g007
Figure 8. The variation of Nusselt number with the electrokinetic width K for shear thinning nanofluid, Newtonian nanofluid, shear thickening nanofluid when the Joule heating parameters are given as (a) S = −5; (b) S = 0; (c) S = 5.
Figure 8. The variation of Nusselt number with the electrokinetic width K for shear thinning nanofluid, Newtonian nanofluid, shear thickening nanofluid when the Joule heating parameters are given as (a) S = −5; (b) S = 0; (c) S = 5.
Micromachines 10 00363 g008
Figure 9. The variation of Nusselt number with Joule heating parameter for shear thinning nanofluid, Newtonian nanofluid, and shear thickening nanofluid in the case of (a) ϕ = 0 and (b) ϕ = 0.06.
Figure 9. The variation of Nusselt number with Joule heating parameter for shear thinning nanofluid, Newtonian nanofluid, and shear thickening nanofluid in the case of (a) ϕ = 0 and (b) ϕ = 0.06.
Micromachines 10 00363 g009
Figure 10. The variation of (a) Nusselt number and (b) Bejan number with flow behavior index for different nanoparticle volume fraction.
Figure 10. The variation of (a) Nusselt number and (b) Bejan number with flow behavior index for different nanoparticle volume fraction.
Micromachines 10 00363 g010
Table 1. Typical values of the parameters.
Table 1. Typical values of the parameters.
Parameters (Notation) Value (unit)
The fluid permittivity ε 7.08 × 10−10 C·V−1·m−1
Boltzmann constant kb 1.38 × 10−23 J·K−1
Absolute temperature T0 293 K
Elementary charge e 1.6 × 10−19 C
Half channel height and width, a and b 10 × 10−6 m
Fluid electrical conductivity σ1.2639 × 10−7 S·m−1
Flow consistency index of power-law fluid η 9 × 10−4 N·m2·sn
Viscosity of Newtonian fluid μ0 9 × 10−4 N·m−2·s
Zeta potential ξ 0.025 V
Thermal conductivity of the solid nanoparticle ks 40 W·m−1·K−1
Thermal conductivity of the base fluid kf 0.618 W·m−1·K−1
Valence of ions χ1
Electrokinetic width K 5–50
Flow behavior index n 0.7–1.4
Nanoparticle volume fraction ϕ 0.0–0.06
Joule heating Parameter S −5–5
Ratio of nanolayer thickness to original particle radius ω 1.1

Share and Cite

MDPI and ACS Style

Deng, S. Thermally Fully Developed Electroosmotic Flow of Power-Law Nanofluid in a Rectangular Microchannel. Micromachines 2019, 10, 363. https://doi.org/10.3390/mi10060363

AMA Style

Deng S. Thermally Fully Developed Electroosmotic Flow of Power-Law Nanofluid in a Rectangular Microchannel. Micromachines. 2019; 10(6):363. https://doi.org/10.3390/mi10060363

Chicago/Turabian Style

Deng, Shuyan. 2019. "Thermally Fully Developed Electroosmotic Flow of Power-Law Nanofluid in a Rectangular Microchannel" Micromachines 10, no. 6: 363. https://doi.org/10.3390/mi10060363

APA Style

Deng, S. (2019). Thermally Fully Developed Electroosmotic Flow of Power-Law Nanofluid in a Rectangular Microchannel. Micromachines, 10(6), 363. https://doi.org/10.3390/mi10060363

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop