Impact of asymmetric zeta-potential on natural convection flow in a vertical microannulus with electroosmotic and Joule heating effects

Abstract

Analytical solutions of temperature distributions and flow formation of joint buoyancy and electroosmotic flow in a vertical microtube formed by two concentric microcylinders are presented. The central equations describing flow formations and thermal energy are offered and solved in closed-form in terms of modified Bessel’s function of first and second kinds. These solutions have significant application in predicting and analyzing flow formation and thermal behavior of Newtonian fluids in micropumps and microchips. Based on the exact solutions obtained, the effects of flow parameters are clearly explained with the use of line graphs. Based on the simulation of results using MATLAB, it is found that fluid temperature distributions and fluid velocity in the vertical microannulus could be enhanced by increasing the radius ratio of the concentric microcylinders.

Keywords

Electroosmotic flow zeta-potential microannulus Joule heating

Introduction

The study of electroosmotic flow (EOF) gains its relevance in micro geometries, due to the comparable channel width to the electric double layer length. This phenomenon is generated by the combined application of low electric voltage at the surfaces of the microchannel and the interaction between the walls and the charged fluid. A lot has been discussed over the years about this phenomenon owing to its significance in micromechanical systems (MEMS). Experimental existence of the phenomenon and theoretical background were respectively presented by Reuss¹ and Probstein.² Since then, many applicable extensions for this phenomenon have been carried out in the literature. Rice and Whitehead³ in a narrow cylindrical capillary, discussed the significance of electrokinetic flow in the presence of low zeta potential. Kang et al.⁴ revisited the work of Rice and Whitehead ³ in the presence of with high zeta potentials. The electrokinetic flow in ultrafine capillary slits was carried out by Burgreen and Nakache⁵ while Debye and Hückel⁶ offered a linearized form of the nonlinear Poisson equation governing the electric potential in the EDL for low electric potential. Also, Tsao⁷ examined the EOF in an annulus, whereas, Wang et al.⁸ gave an analytical solution of EOF in a semi-circular microchannel. In all these articles, the influence of buoyancy force induced by unequal heating at the surfaces of the microcylinders was not accounted for. Other related articles can be seen in Refs.^9–17

The natural convection flow (NCF) is a kind of flow generated as a consequence of unequal density induced by unequal wall temperature. This kind of heat transfer is significant in cooling nuclear reactors, automobiles and other devices at microlevel. The investigation of NCF in microgeometries extends its applicability in cooling microchips and micro-electro-mechanical systems (MEMS). Some of the earliest works on this phenomenon include the works of Elenbaas,¹⁸ Ostrach,¹⁹ Sparrow and Gregg²⁰ and Bodoia and Osterle²¹ where they respectively studied the NCF in parallel plates, with and without heat sources, in vertical plate with uniform heat flux and uniform temperature. Later, Aung²² and Aung et al.²³ analyzed respectively the fully developed and developing free convection flow between vertical flat plates. Since then, many articles on different physics of flows have been published. To mention a few, Weng and Chen²⁴ studied the drag reduction of NCF in vertical microchannel while Weng and Chen²⁵ investigated the NCF in a vertical microchannel. Jha and Oni²⁶ studied the time-dependent NCF in a vertical annulus heated/cooled asymmetrically. Jha et al.²⁷ scrutinized the mixed convection flow in a vertical microannulus with heat source/sink. They gave separate solutions for the heat source/sink cases and concluded that heat generating case has higher velocity and fluid temperature compared with the case of heat sink in the microannulus. Other related articles on NCF in an annular geometry include Refs.^28–33

In spite the contributions above, precise theory of electroosmotic flow in vertical microannulus with buoyancy force and asymmetric heat flux is still rare. Such solution has significant application in chemical separation at microlevel, which could be used in drug delivery to a targeted cell in the body. In view of this, the purpose of this article is to discuss the role of asymmetric zeta-potential on NCF in a vertical microannulus with electroosmotic and Joule heating effect. The novelty of this article is the presentation of analytical mathematical models to predict flow formation and skin-friction at the surfaces of the microcylinders of NCF in a vertical microannulus with electroosmotic effect.

Mathematical study

Consider the flow on an incompressible balance $(z : z)$ electrolyte, where $z$ is the valence of the electrolyte in a vertical microannulus induced by the combined influence of buoyancy (caused by density difference), electroosmotic force caused by the EDL at the interface between the microannulus surfaces and the fluid and an externally applied electric field as revealed in Figure 1. The outer and inner cylinders are of radii $k_{2}$ and $k_{1}$ respectively such that $k_{2} > k_{1}$ and the length of the cylinders is strictly greater than the width $(k_{2} - k_{1})$ of the microannulus. Constant isofluxes, $q_{1}$ and $q_{2}$ are prescribed at the inner and outer cylinders respectively.

Figure 1.

Schematic diagram of the problem.

The flow is presumed to be hydrodynamically and thermally fully developed. Other assumptions made include:

(i) The flow is laminar with equivalent Reynolds number less than unity.

(ii) The charge in the EDL agrees with the Boltzmann distribution and the Debye length is strictly smaller than the annular gap.

(iii) The EDL width is not affected by the temperature gradients.

(iv) The temperature jump effect as well as the velocity slip are incorporated at the microannulus surfaces.

For fully developed flow, the flow becomes unidirectional ( $\vec{V} = 0 \vec{r} + 0 \vec{ϑ} + u (r) \vec{z}$ ). Following,³⁴ the governing electric potential, continuity, momentum and energy equations in dimensional form respectively become:

\frac{1}{r} \frac{d}{dr} (r \frac{d ψ^{'}}{dr}) = \frac{2 F z^{'} C_{0}}{ε'} \sinh (\frac{z^{'} F ψ^{'}}{R^{'} \hat{T}})

(1)

\frac{\partial u}{\partial z} = 0

(2)

\frac{\partial u}{\partial t} = \frac{μ}{r} \frac{\partial}{\partial r} (r \frac{\partial u}{\partial r}) + β g ρ_{0} (T - T_{0}) + ρ_{e} E'_{z}

(3)

\begin{matrix} \frac{\partial T}{\partial t} + u \frac{\partial T}{\partial z} = \frac{k^{'}}{ρ_{0} C_{p}} \frac{1}{r} \frac{\partial}{\partial r} (r \frac{\partial T}{\partial r}) + \frac{{E^{' 2}}_{z}}{ρ_{0} C_{p} σ} \cosh^{2} \\ (\frac{zF ψ^{'}}{R^{'} \hat{T}}) + \frac{μ}{ρ_{0} C_{p}} {(\frac{\partial u}{\partial r})}^{2} \end{matrix}

(4)

subject to

\begin{matrix} ψ^{'} = ζ_{1}, \frac{\partial T}{\partial r} = - \frac{q_{1}}{k}, u (r) = α' \frac{\partial u}{\partial r} at r = k_{1} \\ ψ^{'} = ζ_{2}, \frac{\partial T}{\partial r} = \frac{q_{2}}{k}, u (r) = - α' \frac{\partial u}{\partial r}, at r = k_{2} \end{matrix}

(5)

where $ψ^{'}$ is the dimensional electric potential due to EDL, $F$ is the Faraday’s constant, “ $z^{'}$ is the valence number of ions in the solution,” $ρ_{0}$ is the density of the fluid to $T_{0}$ , $T_{0}$ is the ambient fluid temperature and $q_{i = 1, 2}$ are the prescribed heat fluxes at the surfaces of the microcylinders. Equation (1) is the Poisson-Boltzmann equation which describes the EDL distributions in the microannulus, the fourth term of equation (3) represents the electric body force, which is obtained from the combination of applied electric gradient and EDL potential. While the last two terms in the energy equation (4) signify the viscous dissipation and Joule heating terms respectively. Equation (2) on the other hand represents the continuity equation. All other physical quantities employed in equations (1)–(5) are well defined in nomenclature.

Introducing the following dimensionless quantities and parameters¹³:

\begin{matrix} R = \frac{r - k_{1}}{w}, Z = \frac{z}{w}, w = k_{2} - k_{1}, T = \frac{θ q_{2} w}{k} + T_{s 1}, \\ U = \frac{u}{u_{0}}, α = δ = \frac{λ}{w}, Gr = \frac{g β ρ_{0} q_{2} w^{4}}{κ μ ν}, \\ ξ_{1}^{*} = \frac{z^{'} F ζ_{1}}{R^{'} T_{0}}, ξ_{2}^{*} = \frac{z^{'} F ζ_{2}}{R^{'} T_{0}}, Γ = \frac{ε^{'} R^{'} \hat{T} ζ_{1}}{z^{'} F ρ_{0} ν^{2}}, \\ E_{z} = \frac{{E^{'}}_{z} w}{ζ_{1}}, η = \frac{k_{1}}{k_{2}}, s = \frac{ζ_{1} w}{σ (T_{1} - T_{0})}, ψ = \frac{z^{'} F ψ'}{R^{'} T_{0}}, \\ κ = \frac{w}{λ_{D}}, λ_{D} = {[\frac{ε^{'} R^{'} T_{0}}{2 F^{2} z^{' 2} C_{0}}]}^{1 / 2}, Ra = \frac{ρ_{0} c_{p_{0}} u_{0} w}{k^{'}}, \\ r_{q} = \frac{q_{1}}{q_{2}}, α = \frac{α^{'}}{w} \end{matrix}

(6)

Considering the steady, fully developed region and a very small EDL size relative to $w$ , so that equations (1), (3), and (4) on applying equation (6) respectively become:

\frac{1}{(η + (1 - η) R)} \frac{d}{dR} ((η + (1 - η) R) \frac{d ψ}{dR}) - k^{2} ψ = 0

(7)

\begin{matrix} \frac{1}{(η + (1 - η) R)} \frac{d}{dR} ((η + (1 - η) R) \frac{dU}{dR}) \\ + Gr θ - E_{z} Γ k^{2} ψ = 0 \end{matrix}

(8)

\begin{matrix} RaU \frac{\partial θ}{\partial Z} = \frac{1}{(η + (1 - η) R)} \frac{\partial}{\partial R} ((η + (1 - η) R) \frac{\partial θ}{\partial R}) \\ + s {E_{z}}^{2} \end{matrix}

(9)

Subject to

\begin{matrix} \frac{d θ}{dR} = - r_{q}, U (R) = α \frac{dU}{dR} at R = 0 \\ \frac{d θ}{dR} = 1, U (R) = - α \frac{dU}{dR}, at R = 1 \end{matrix}

(10)

By using the transformation $χ = η + (1 - η) R$ , it is convenient to examine the influence of annular gap on flow formation and skin-friction at the surfaces of the microannulus. Equations (7)–(10) are transformed as follows²⁴:

\frac{1}{χ} \frac{d}{d χ} (χ \frac{d ψ}{d χ}) - κ^{2} ψ = 0

(11)

\frac{1}{χ} \frac{d}{dR} (χ \frac{dU}{d χ}) + \frac{Gr θ}{{(1 - η)}^{2}} - \frac{E_{z} Γ κ^{2} ψ}{{(1 - η)}^{2}} = 0

(12)

\frac{RaU}{{(1 - η)}^{2}} \frac{\partial θ}{\partial Z} = \frac{1}{χ} \frac{\partial}{\partial χ} (χ \frac{\partial θ}{\partial χ}) + \frac{s {E_{z}}^{2}}{{(1 - η)}^{2}}

(13)

Subject to:

\begin{matrix} ψ (η) = ξ_{1}^{*}, \frac{d θ}{d χ} = \frac{- r_{q}}{(1 - η)}, U (χ) = α (1 - η) \frac{dU}{d χ} at χ = η \\ ψ (1) = ξ_{2}^{*}, \frac{d θ}{d χ} = \frac{1}{(1 - η)}, U (χ) = - α (1 - η) \frac{dU}{dR}, at χ = 1 \end{matrix}

(14)

where $κ = \frac{k}{(1 - η)}$ , $ξ_{1}^{*}, ξ_{2}^{*}$ are the non-dimensional asymmetric zeta-potentials at the microannulus surfaces.

The equations above describe the mathematical equations governing electroosmotic NCF in a vertical microannulus with Joule heating effect. The term $r_{q}$ explains the level of unevenness in heat fluxes whereas $α$ signifies the dimensionless velocity slip-length.

For fully developed flow, the expression $\frac{\partial θ}{\partial Z}$ is a constant and is computed by multiplying through by $χ$ , integrating equation (13) and substituting equation (14) as follows:

\frac{\partial θ}{\partial Z} = \frac{1}{Ra} [\frac{2 (1 + η r_{q}) + s {E_{z}}^{2} (1 + η)}{(1 + η)}]

(15)

So that equation (13) becomes:

\frac{1}{χ} \frac{d}{d χ} (χ \frac{d θ}{\partial χ}) + \frac{s {E_{z}}^{2}}{{(1 - η)}^{2}} = \frac{2 (1 + η r_{q}) + s {E_{z}}^{2} (1 + η) U}{(1 + η) {(1 - η)}^{2}}

(16)

The solution of equation (11) with boundary condition (14) can be conveniently written in form of modified Bessel’s function of first and second kinds of order zero as (Olver et al.³⁵):

ψ (χ) = C_{1} I_{0} (χ κ) + C_{2} K_{0} (χ κ)

(17)

where $K_{0}$ and $I_{0}$ are the modified Bessel’s function of second and first kinds of order zero and $C_{1}$ and $C_{2}$ are constants obtained as:

\begin{matrix} C_{1} = \frac{ξ_{2}^{*} K_{0} (η κ) - ξ_{1}^{*} K_{0} (κ)}{I_{0} (κ) K_{0} (η κ) - I_{0} (η κ) K_{0} (κ)}, \\ C_{2} = \frac{ξ_{1}^{*} I_{0} (η κ) - ξ_{2}^{*} I_{0} (κ)}{I_{0} (κ) K_{0} (η κ) - I_{0} (η κ) K_{0} (κ)} \end{matrix}

(18)

Equations (12) and (16) are the mathematical models governing the velocity profile and temperature distributions in the microannulus respectively. These equations are coupled and linear, so that their solutions are obtained by first making $θ$ subject of formula from equation (12) as follows:

θ = \frac{{(1 - η)}^{2}}{Gr} [\frac{E_{z} Γ κ^{2} ψ}{{(1 - η)}^{2}} - \frac{1}{χ} \frac{d}{dR} (χ \frac{dU}{d χ})]

(19)

So that on applying the operator $\frac{1}{χ} \frac{d}{dR} (χ \frac{d}{d χ})$ to both sides of (19) and comparing with equation (16), we have:

\begin{matrix} \frac{1}{χ} \frac{d}{dR} (χ \frac{d}{d χ}) \frac{1}{χ} \frac{d}{dR} (χ \frac{dU}{d χ}) \\ + \frac{Gr [2 (1 + η r_{q}) + s {E_{z}}^{2} (1 + η)] U}{(1 + η) {(1 - η)}^{4}} \\ = \frac{sGr {E_{z}}^{2}}{{(1 - η)}^{4}} + \frac{E_{z} Γ k^{2} ψ}{{(1 - η)}^{4}} \frac{1}{χ} \frac{d}{dR} (χ \frac{d ψ}{d χ}) \end{matrix}

(20)

The homogeneous part of equation (20) can be written as:

(L^{4} - m^{4}) U = 0

(21)

where $L^{2} = \frac{1}{χ} \frac{d}{dR} (χ \frac{d}{d χ})$ , and $m = \frac{i^{1 / 2} {Gr [2 (1 + η r_{q}) + s {E_{z}}^{2} (1 + η)]}^{1 / 4}}{{(1 + η)}^{1 / 4} (1 - η)}$ .

So that equation (21) becomes:

(\frac{1}{χ} \frac{d}{dR} (χ \frac{dU}{d χ}) + m^{2} U) (\frac{1}{χ} \frac{d}{dR} (χ \frac{dU}{d χ}) - m^{2} U) = 0

(22)

Equation (22) is the combination of standard Bessel’s function and modified Bessel’s function. On solving equation (20) with boundary conditions (14), we obtained the velocity profile in closed form as:

\begin{matrix} U = C_{3} I_{0} (m χ) + C_{4} K_{0} (m χ) + C_{5} J_{0} (m χ) + C_{6} Y_{0} (m χ) \\ + \frac{E_{z} Γ κ^{4} [C_{1} I_{0} (χ κ) + C_{2} K_{0} (χ κ)]}{κ^{4} - m^{4}} \\ - \frac{2 s {E_{z}}^{2} (1 + η)}{2 (1 + η r_{q}) + s {E_{z}}^{2} (1 + η)} \end{matrix}

(23)

By using equation (19), the temperature distributions in the microannulus is obtained as:

\begin{matrix} θ = \frac{E_{z} Γ}{Gr} [κ^{2} - \frac{κ^{6} {(1 - η)}^{2}}{(κ^{4} - m^{4})}] [C_{1} I_{0} (χ κ) + C_{2} K_{0} (χ κ)] \\ - \frac{{(1 - η)}^{2} m^{2}}{Gr} [C_{3} I_{0} (m χ) + C_{4} K_{0} (m χ) \\ - C_{5} J_{0} (m χ) - C_{6} Y_{0} (m χ)] \end{matrix}

(24)

By putting the boundary conditions in equations (14) on equations (23) and (24), the constants $C_{3}, C_{4}, C_{5},$ and $C_{6}$ are obtained as follow; “

\begin{matrix} C_{3} = \frac{det (\begin{matrix} X_{5} & X_{2} & X_{3} & X_{4} \\ X_{10} & X_{7} & X_{8} & X_{9} \\ X_{16} & X_{13} & X_{14} & X_{15} \\ X_{21} & X_{18} & X_{15} & X_{20} \end{matrix})}{det (\begin{matrix} X_{1} & X_{2} & X_{3} & X_{4} \\ X_{6} & X_{7} & X_{8} & X_{9} \\ X_{12} & X_{13} & X_{14} & X_{15} \\ X_{17} & X_{18} & X_{19} & X_{20} \end{matrix})} C_{4} = \frac{det (\begin{matrix} X_{1} & X_{5} & X_{3} & X_{4} \\ X_{6} & X_{10} & X_{8} & X_{9} \\ X_{12} & X_{16} & X_{14} & X_{15} \\ X_{17} & X_{21} & X_{15} & X_{20} \end{matrix})}{det (\begin{matrix} X_{1} & X_{2} & X_{3} & X_{4} \\ X_{6} & X_{7} & X_{8} & X_{9} \\ X_{12} & X_{13} & X_{14} & X_{15} \\ X_{17} & X_{18} & X_{19} & X_{20} \end{matrix})} \\ C_{5} = \frac{det (\begin{matrix} X_{1} & X_{2} & X_{5} & X_{4} \\ X_{6} & X_{7} & X_{10} & X_{9} \\ X_{12} & X_{13} & X_{16} & X_{15} \\ X_{17} & X_{18} & X_{21} & X_{20} \end{matrix})}{det (\begin{matrix} X_{1} & X_{2} & X_{3} & X_{4} \\ X_{6} & X_{7} & X_{8} & X_{9} \\ X_{12} & X_{13} & X_{14} & X_{15} \\ X_{17} & X_{18} & X_{19} & X_{20} \end{matrix})} C_{6} = \frac{det (\begin{matrix} X_{1} & X_{2} & X_{3} & X_{5} \\ X_{6} & X_{7} & X_{8} & X_{10} \\ X_{12} & X_{13} & X_{14} & X_{16} \\ X_{17} & X_{18} & X_{19} & X_{21} \end{matrix})}{det (\begin{matrix} X_{1} & X_{2} & X_{3} & X_{4} \\ X_{6} & X_{7} & X_{8} & X_{9} \\ X_{12} & X_{13} & X_{14} & X_{15} \\ X_{17} & X_{18} & X_{19} & X_{20} \end{matrix})} \end{matrix}

(25)

where det represents the determinant of the matrix and $X_{i' s}; i = 1, 2, 3, \dots$ are constants defined in Appendix.

In order to incorporate the temperature jump effect, following (Avci and Aydin³⁶), the temperature distribution in the vertical microchannel is obtained as:

\begin{matrix} θ^{*} (Y) = \frac{(T - T_{1})}{q_{2} D_{h} / k} = \frac{E_{z} Γ}{Gr} [κ^{2} - \frac{κ^{6} {(1 - η)}^{2}}{(κ^{4} - m^{4})}] \\ [C_{1} I_{0} (χ κ) + C_{2} K_{0} (χ κ)] \\ - \frac{{(1 - η)}^{2} m^{2}}{Gr} [C_{3} I_{0} (m χ) + C_{4} K_{0} (m χ) \\ - C_{5} J_{0} (m χ) - C_{6} Y_{0} (m χ)] - β_{t} r_{q} \end{matrix}

(26)

From the velocity solution obtained, the skin frictions at the micro-annular surfaces are given as:

\begin{matrix} {\frac{dU}{d χ} |}_{χ = η} = (1 - η) m [C_{3} I_{1} (m η) - C_{4} K_{1} (m η) \\ - C_{5} J_{1} (m η) - C_{6} Y_{1} (m η)] \\ + \frac{(1 - η) E_{z} Γ κ^{5}}{(κ^{4} - m^{4})} {C_{1} I_{1} (κ η) - C_{2} K_{1} (κ η)} \end{matrix}

(27)

\begin{matrix} {\frac{dU}{d χ} |}_{χ = 1} = (1 - η) m [C_{3} I_{1} (m) - C_{4} K_{1} (m) - C_{5} J_{1} (m) \\ - C_{6} Y_{1} (m)] + \frac{(1 - η) E_{z} Γ κ^{5}}{(κ^{4} - m^{4})} {C_{1} I_{1} (κ) - C_{2} K_{1} (κ)} \end{matrix}

(28)

Also, the volumetric flow-rate in the microannulus is obtained as:

\begin{array}{l} Q = \frac{1}{{(1 - η)}^{2}} \int_{η}^{1} χ U (χ) d χ \\ = \frac{1}{{(1 - η)}^{2}} {\begin{cases} \frac{C_{3}}{m} [I_{1} (m) - η I_{1} (m η)] - \frac{C_{4}}{m} [K_{1} (m) - η K_{1} (m η)] - \\ \frac{C_{5}}{m} [J_{1} (m) - η J_{1} (m η)] - \frac{C_{6}}{m} [Y_{1} (m) - η Y_{1} (m η)] + \frac{E_{z} Γ κ^{3}}{(κ^{4} - m^{4})} \\ {C_{1} [I_{1} (κ) - η I_{1} (κ η)] - C_{2} [K_{1} (κ) - η K_{1} (κ η)]} - \frac{s E_{z}^{2} (1 - η^{2}) (1 + η)}{2 (1 + η r_{q}) + s E_{z}^{2} (1 + η)} \end{cases}} \end{array}

(29)

Results and discussion

The discussion of results for electroosmotic NCF in a vertical microannulus with asymmetric zeta-potentials at the surfaces of the microcylinders in the existence of Joule heating is carried out in the section. The solutions gotten showed that the sundry parameters influencing flow formations are ratio of radii $(η)$ which is inversely proportional to annular gap, slip-parameter $(α)$ , Debye Hückel parameter $(κ)$ which is inversely proportional to the EDL size, Joule heating parameter $(s)$ which is proportioanl to the thermal energy, ratio of heat fluxes $(r_{q})$ , magnitude of electroosmotic force $(Γ)$ , and the Grashof number $(Gr)$ which defines the strength of the bouyant force. These parameters have significant applications in engineering designs and understanding the physics of flow formation. Therefore, it is expedient to understand how these parameters affect flow formation and minimization of skin-friction in the system. To ascertain this objectives, line graphs are generated for the analytical solutions obtained using Matlab software. Except otherwise stated, we have use the values of $0 < η < 1$ , with fixed value of $η = 0.5$ , $0 \leq α \leq 0.1$ to represent the no-slip $(α = 0.0)$ and slip $(α = 0.1)$ domain, $0 < κ \leq 30$ , $0 < s \leq 5$ , $0 \leq r_{q} \leq 1$ , $0 < Gr \leq 50$ .

Figures 2 and 3 respectively depict the fluid temperature and velocity profile in the vertical microannulus as functions of radius ratio and asymmetry in wall zeta-potential for constant value of other leading parameters. It is noticeable that fluid temperature as well as velocity respond noticeably to change in $η$ and $ξ_{2}^{*}$ . In fact, fluid temperature and velocity both increase with increase in $η$ regardless the value of $ξ_{2}^{*}$ throughout the microannulus. From both figures, the maximum temperature as well as velocity are accomplished for the instance of symmetry zeta-potential at the surfaces of the microannulus regardless the values of $η$ . This is expected, because the annular gap varies inversely with $η$ , which sequentially increases the thermal energy within the small annular gap and thereby increasing the overall fluid temperature and velocity throughout the microannulus.

Figure 2.

Temperature distributions for different values of $η$ and $ξ_{2}^{*}$ .

Figure 3.

Velocity profile for different values of $η$ and $ξ_{2}^{*}$ .

In addition, while the parabolic nature of fluid velocity is preserved, we found a nearly constant temperature throughout the microannulus.

It is pertinent to also understand the role of buoyancy term on temperature distributions as well as velocity profile in the microannulus. To this end, Figures 4 and 5 respectively represent the temperature distributions and fluid velocity as functions of Grashof number $(Gr)$ and $ξ_{2}^{*}$ . It is established from Figure 4 that the impact of buoyancy term is to reduce fluid temperature throughout the microannulus. On the other hand (Figure 5), velocity decreases as $Gr$ increases, around the center of the annulus but increases along the inner surface of the outer cylinder. From both figures, it is obvious that the maximum temperature and velocity can be obtained by subjecting the surfaces of the microannulus to equal zeta-potentials. In addition, for the situation of unequal zeta-potentials, two points of inflection (turning points) are observed of which at these points, fluid velocity is independent of $Gr$ .

Figure 4.

Temperature distributions for different values of $Gr$ and $ξ_{2}^{*}$ .

Figure 5.

Velocity profile for different values of $Gr$ and $ξ_{2}^{*}$ .

For proper investigation of the impact of sundry parameters on volumetric flow-rate, Figures 6 to 8 depict the effect of various leading parameters on the amount of fluids passing through the microannulus (volumetric flow-rate). This phenomenon has significant importance in optimization of flow rate in design of micropumps and drug targeting. Figure 6 revealed that the maximum volumetric flow-rate can be achieved by exposing the microannulus surfaces to equal wall zeta-potentials. Also, flow-rate is found to first increase for small values of $κ$ (corresponding to large EDL) and then gradually decreases until flow formation transforms from combined buoyancy-electroosmotic flow to purely buoyancy driven flow. Throughout the microannulus, the volumetric flow-rate is found to decline with increase in $η$ regardless the values of $ξ_{2}^{*}$ and $κ$ . This could be attributed to the fact that increasing $η$ decreases the annular gap, which then distorted the free motion of flow formation. Furthermore, a relook at this figure shows that as $κ \to \infty$ , the volumetric flow-rate is independent of $η, ξ_{2}^{*}$ and $κ$ . This is physically true as increase in $κ$ corresponds to a smaller EDL, which gradually diminishes the electroosmotic flow effect as $κ \to \infty$ .

Figure 6.

Volumetric flow-rate for different values of $η, κ$ , and $ξ_{2}^{*}$ .

Figure 7.

Volumetric flow-rate for different values of $κ, Γ$ , and $s$ .

Figure 8.

Volumetric flow-rate for different values of $κ, Gr$ , and $r_{q}$ .

Figure 7 also discussed the volumetric flow-rate as functions of $Γ$ and $s$ in the microannulus. It is found that volumetric flow-rate rises with surge in Joule heating parameter. In addition, the maximum volumetric flow-rate is achieved for the case of electroosmotic flow. It is worthy to state that $Γ = 0$ signifies the absence of electroosmotic effect while $Γ = 1$ represents the presence of electroosmotic phenomenon with unity magnitude. This result shows that the volumetric flow-rate could have been wrongly computed should the EDL phenomenon at the surfaces of the microcylinders (which is often present for flows in microgeometries in the presence of electric field) had not been considered.

Figure 8 portrays the influence of $Gr$ and $r_{q}$ “on the volumetric flow-rate. This figure demonstrated that volumetric flow rate decreases with increase in $Gr$ and $r_{q}$ for $κ < 10$ while a uniform flow rate is found for higher values of $κ$ .

Figure 9 shows the effect of unequal wall zeta-potentials, EDL size and ratio of radiuses on skin-frictions at the outer surface of inner cylinder. This figure discloses that the skin-friction for the symmetric wall zeta potential $(ξ_{2}^{*} = 1.0)$ are significantly different from those of asymmetric wall zeta potential, regardless of the annular gap or EDL size. Also, the skin-frictions at this surface is seen to rise with growth in radius ratio and EDL size. This is because increase radius ratio results to a decrease in annular gap, which increases the shear stress at the surfaces of the cylinders therefore increasing skin-friction.

Figure 9.

Skin-friction at the wall $R = 0$ for different values of $ξ_{2}^{*}, κ,$ and $η$ .

Figure 10 depicts the combined impact of slip parameter, magnitude of electroosmotic force and radius ratio on skin-friction at the outer surface of inner cylinder. It is observed from this figure that the role of electroosmotic force is to reduce skin-friction at this surface regardless of the flow regime. Also, the slip flow (which is a closed to real life physics of flow in microannulus) is found to further has a lower skin-friction compared to the no-slip case. This outcome clearly demonstrated that the skin-friction at this surface would have been overestimated, if the surface slip-parameter was not put into consideration.

Figure 10.

Skin-friction at the wall $R = 0$ for different values of $α, Γ$ , and $η$ .

Figure 11 displays the role of $Gr$ and $r_{q}$ on skin-friction at the wall $R = 0$ . It is noticed from this figure that both the buoyancy parameter as well as ratio of heat fluxes (which defines the magnitude of asymmetric wall heat fluxes) respond sparingly to skin-friction. It is also found that the skin-friction is a decreasing function $Gr$ regardless the value of $r_{q}$ .

Figure 11.

Skin-friction at the wall $R = 0$ for different values of $Gr, r_{q}$ , and $η$ .

It is noteworthy to mention that in the absence of electroosmotic force and joule heating effect, this work corresponds to the work of Avci and Aydin³³ when the radius ratio $η$ approaches 1. Also, the finding contained here reduces to the work of Oni and Jha¹³ once the Joule heating effect is relaxed and the geometry is transformed from microannulus to microchannel.

Conclusion

Analytical solution for electroosmotic NCF in a vertical microannulus with asymmetric zeta-potential and Joule heating is presented in this section. The equation of motion, electric potential, temperature distributions are derived and solved exactly. The influences of relevant parameters are presented using line graphs generated by MATLAB. The major conclusions are:

Fluid velocity, temperature, skin-friction and volumetric flow rate respond sharply to asymmetry in zeta-potential.

The volumetric flow rate as well as the skin-friction would have been erroneously calculated should the EDL as well as the slip parameters had not been put into consideration.

Fluid temperature and velocity profiles in the vertical microannulus increase with increase in radius ratio.

The role of electroosmotic force is to increase volumetric flow-rate and decrease skin-friction.

Footnotes

Appendix

\begin{matrix} X_{1} = I_{0} (m η) - m (1 - η) α I_{1} (m η), \\ X_{2} = K_{0} (m η) + m (1 - η) α K_{1} (m η) \\ X_{3} = J_{0} (m η) + m (1 - η) α J_{1} (m η), \\ X_{4} = Y_{0} (m η) + m (1 - η) α Y_{1} (m η) \\ X_{5} = \frac{{sE}_{Z}^{2}}{m^{4} {(1 - η)}^{2}} - \frac{E_{Z} κ^{4}}{(κ^{4} - m^{4})} \\ {C_{1} [I_{0} (κ η) - κ (1 - η) α I_{1} (κ η)] \\ + C_{2} [K_{0} (κ η) + κ (1 - η) α K_{1} (κ η)]}, \\ X_{6} = I_{0} (m) + m (1 - η) α I_{1} (m), \\ X_{7} = K_{0} (m) - m (1 - η) α K_{1} (m), \\ X_{8} = J_{0} (m) - m (1 - η) α J_{1} (m), \\ X_{9} = Y_{0} (m) - m (1 - η) α Y_{1} (m), \\ X_{10} = \frac{{sE}_{Z}^{2}}{m^{4} {(1 - η)}^{2}} - \frac{E_{Z} κ^{4}}{(κ^{4} - m^{4})} \\ {C_{1} [I_{0} (κ) + κ (1 - η) α I_{1} (κ)] \\ + C_{2} [K_{0} (κ) + κ (1 - η) α K_{1} (κ)]}, \end{matrix}

\begin{matrix} X_{11} = \frac{E_{Z} Γ κ^{3}}{Gr} [1 - \frac{κ^{4} {(1 - η)}^{2}}{(κ^{4} - m^{4})}], X_{12} = \frac{{(1 - η)}^{2} m^{3}}{Gr} \\ I_{1} (m η), X_{13} = - \frac{{(1 - η)}^{2} m^{3}}{Gr} K_{1} (m η), \\ X_{14} = \frac{{(1 - η)}^{2} m^{3}}{Gr} J_{1} (m η), X_{15} = \frac{{(1 - η)}^{2} m^{3}}{Gr} \\ Y_{1} (m η), X_{16} = \frac{r_{q}}{(1 - η)} + X_{11} \\ [C_{1} I_{1} (κ η) - C_{2} K_{1} (κ η)], \\ X_{17} = \frac{{(1 - η)}^{2} m^{3}}{Gr} I_{1} (m), X_{18} = - \frac{{(1 - η)}^{2} m^{3}}{Gr} K_{1} (m), \\ X_{19} = \frac{{(1 - η)}^{2} m^{3}}{Gr} J_{1} (m) \\ X_{20} = \frac{{(1 - η)}^{2} m^{3}}{Gr} Y_{1} (m), X_{21} = X_{11} \\ [C_{1} I_{1} (κ) - C_{2} K_{1} (κ)] - \frac{1}{(1 - η)} \end{matrix}

Acknowledgements

The authors are thankful to the University as well as the Department.

Author contributions

Author Roles: M.O.O. solve the problem, prepared the writeups, write the codes and generates all figures. B.K.J. formulated the problem and proofread.

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iD

Michael O Oni

References

Reuss

Sur un nouvel effet de l’électricité galvanique. Mémoires de la Société Impériale des Naturalistes de Moscou 1809; 2: 327–337.

Probstein

Physicochemical hydrodynamics: an Introduction. New York: John Wiley & Sons, Inc, 1994.

Rice

Whitehead

Electrokinetic flow in a narrow cylindrical capillary. J Chem Phys 1965; 69(11): 4017–4024.

Kang

Yang

Huang

Electroosmotic flow in a capillary annulus with high zeta potentials. J Colloid Interface Sci 2002; 253(2): 285–294.

Burgreen

Nakache

FR.

Electrokinetic flow in ultrafine capillary slits. J Chem Phys 1964; 68: 1084–1091.

Debye

Hückel

The theory of electrolytes. I. Lowering of freezing point and related phenomena. Phys Z 1923; 24: 185–206.

Tsao

H-K.

Electroosmotic flow through an annulus. J Colloid Interface Sci 2000; 225: 247–250.

Wang

Liu

Chang

CC.

Analytical solution of electro-osmotic flow in a semicircular microchannel. Phys Fluids 2008; 20: 063105.

Chakraborty

Ray

Chakraborty

Role of streaming potential on pulsating mass flow rate control in combined electroosmotic and pressure-driven microfluidic devices. Electrophoresis 2012; 33: 419–425.

10.

Jha

Oni

MO.

Mathematical modelling of combined pressure driven and electrokinetic effect in a channel with induced magnetic field: an exact solution. J King Saud Univ - Sci 2019; 31(4): 575–585.

11.

Liechty

Webb

Maynes

RD.

Convective heat transfer characteristics of electro-osmotically generated flow in microtubes at high wall potential. Int J Heat Mass Transf 2005; 48: 2360–2371.

12.

Oni

Jha

BK.

Theoretical analysis of transient natural convection flow in a vertical microchannel with electrokinetic effect. J Taibah Univ Sci 2019; 13(1): 1087–1099.

13.

Oni

Jha

. Electroosmotic natural convection flow in a vertical microchannel with asymmetric heat fluxes. SN Appl Sci 2020; 2(10): 1–14.

14.

Oni

Jha

BK.

Joule heating and viscous dissipation effect on electroosmotic mixed convection flow in a vertical microchannel subjected to asymmetric heat fluxes. Propulsion Power Res 2021; 10: 83–94.

15.

Soong

Wang

SH.

Theoretical analysis of electrokinetic flow and heat transfer in a microchannel under asymmetric boundary conditions. J Colloid Interface Sci 2003; 265: 202–213.

16.

Levine

Marriott

Neale

, et al. Theory of electrokinetic flow in fine cylindrical capillaries at high zeta-potentials. J Colloid Interface Sci 1975; 52: 136–149.

17.

Jha

Oni

MO.

Electromagnetic natural convection flow in a vertical microchannel with Joule heating: exact solution. J Taibah Univ Sci 2018; 12(5): 661–668.

18.

Elenbaas

Heat dissipation of parallel plates by free convection. Physica 1942; 9: 1–28.

19.

Ostrach

Laminar natural convection flow and heat transfer of fluids with and without heat sources in channels with constant wall temperatures, Technical Report, Report no. NACA TN 2863, 1 December 1952. Boston, MA: NACA.

20.

Sparrow

Gregg

JL.

Laminar free convection from a vertical plate with uniform surface heat flux. Trans Am Soc Mech Eng 1956; 78: 435–440.

21.

Bodoia

Osterle

JF.

The development of free convection between heated vertical plates. J Heat Transf 1962; 84: 40–43.

22.

Aung

Fully developed laminar free convection between vertical plates heated asymmetrically. Int J Heat Mass Transf 1972; 15: 1577–1580.

23.

Aung

Fletcher

Sernas

Developing laminar free convection between vertical flat plates with asymmetric heating. Int J Heat Mass Transf 1972; 15: 2293–2308.

24.

Weng

Chen

CK.

Drag reduction and heat transfer enhancement over a heated wall of a vertical annular microchannel. Int J Heat Mass Transf 2009; 52: 1075–1079.

25.

Chen

Weng

HC.

Natural convection in a vertical microchannel. J Heat Transf 2005; 127: 1053–1056.

26.

Jha

Oni

MO.

Transient natural convection flow between vertical concentric cylinders heated/cooled asymmetrically. Proc IMechE, Part A: J Power and Energy 2018; 232: 926–939.

27.

Jha

Oni

Aina

Steady fully developed mixed convection flow in a vertical micro-concentric-annulus with heat generating/absorbing fluid: an exact solution. Ain Shams Eng J 2018; 9(4): 1289–1301.

28.

Oni

MO.

Combined effect of heat source, porosity and thermal radiation on mixed convection flow in a vertical annulus: an exact solution. Eng Sci Technol Int J 2017; 20: 518–527.

29.

Ostrach

Combined natural- and forced-convection laminar flow heat transfer of fluids with and without heat sources in channels with linearly varying wall temperatures. Technical Report. Report no. NACA-TN 3141, 01 April 1954. Boston, MA: NACA.

30.

Schaaf

Chambre

Flow of rarefied gases. Princeton, NJ: Princeton University Press, 1961.

31.

Acharya

On the flow patterns and thermal control of radiative natural convective hybrid nanofluid flow inside a square enclosure having various shaped multiple heated obstacles. Eur Phys J Plus 2021; 136(8): 89.

32.

Acharya

Spectral quasi linearization simulation of radiative nanofluidic transport over a bended surface considering the effects of multiple convective conditions. Eur J Mech B/Fluids 2020; 84: 139–154.

33.

Acharya

Spectral quasi linearization simulation on the radiative nanofluid spraying over a permeable inclined spinning disk considering the existence of heat source/sink. Appl Math Comput 2021; 411: 126547.

34.

Oni

Jha

Ajibade

AO.

Interplay of dual streaming potentials on electroosmotic mixed convection flow in a vertical microannulus with Joule heating effect. Int Commun Heat Mass Transf 2022; 131: 105839.

35.

Olver

FWJ

Lozier

Boisvert

, et al. NIST handbook of mathematical functions. New York, NY: Cambridge University Press, 2010, p.233.

36.

Avcı

Aydın

Mixed convection in a vertical parallel plate microchannel with asymmetric wall heat fluxes. ASME J of Heat Transf 2007; 129: 1091–1095.