Fully Developed Flow of a Nanofluid through a Circular Micropipe in the Presence of Electroosmotic Effects

Niazi, Muhammad D. K.; Xu, Hang

doi:https://doi.org/10.1155/2020/1723256

Mathematical Problems in Engineering

On this page

Abstract Introduction Results Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 1723256 | https://doi.org/10.1155/2020/1723256

Fully Developed Flow of a Nanofluid through a Circular Micropipe in the Presence of Electroosmotic Effects

Muhammad D. K. Niazi¹and Hang Xu¹

Academic Editor: Cedrick A. K. Kwuimy

Received15 Aug 2019

Revised06 Dec 2019

Accepted19 Dec 2019

Published18 Feb 2020

Abstract

Microscale heat sinks based on channels or pipes are designed to restrict the temperatures of microelectromechanical systems, which have a wide range of applications in the modern engineering and mechanics. In this context, this work aims to study heat convection and entropy generation of a fully developed nanofluid flow in a circular micropipe in the presence of an electrical double layer. Buongiorno’s model is employed to exhibit the nanofluid behavior. The governing equations are reduced to a system of nonlinear ordinary differential equations through appropriate similarity transformations. Particularly, we rectify the pressure term as an unknown constant, which makes our flow model compatible with those well-known fluid flow models in macrosize. Highly accurate solutions are obtained and verified. Analysis for physical properties of electric field, velocity field, temperature, and nanoparticle distributions is discussed followed by an investigation of the entropy evolution in the flow. The results show that flow behavior and total entropy of the system depend on the electroosmosis, thermophoresis, and fluid viscosity. However, the influence of the electrical double layer on the flow and system entropy is negligible when the electroosmotic parameter exceeds a maximum value.

1. Introduction

Microelectromechanical systems (MEMS) are of great significance due to their vast applications in modern electronic equipment. In some of these devices, order of heat dissipation is expected to be as high as [1], which could cause breakdowns and failures. In this context, cooling pumps with microsized channels or pipes are deployed. Tuckerman and Pease [2] were the first to suggest the idea of microchannel flow as heat sinks. Their concept was then adopted by several researchers [3, 4] to study the thermal properties and behaviors in such channels/pipes. Particularly, Phillips [5] concluded that the thermal conductivity and the viscosity are two essential properties that affect heat transfer behavior of microchannel flows.

Nanoscale particles, with a diameter of 1 nm–100 nm, are classified as nanoparticles, and a combination of these particles with a base liquid is known as nanofluid. Choi and Eastman [6] were the first to suggest the use of nanoparticles with the base fluid to improve its thermal conductivity. Their approach has been used to enhance the thermal conductivity of fluids in various studies [7–10]. Specifically, a wide range of particles have been used to observe their specific behavior in a base fluid [11, 12]. From a theoretical point of view, several researchers have attempted to mathematically model the nanoparticle behavior. Particularly, Jang and Choi [13] concluded the existence of a link between nanoparticle diffusion and thermal properties of the nanofluid and proposed a thermal conductivity model based on temperature. Various other models were designed to study nanoparticle behavior in nanofluid, where the two-phase models [14–16] and the Tiwari–Das model [17] have been applied widely in the literature.

Buongiorno [18] proposed a better approach in 2006 which incorporated nanoparticle and base fluid properties to obtain a more accurate model. In his study, the fluid layers in the neighborhood of a nanoparticle were considered as a continuum and seven mechanisms of molecular disposition were used. However, thermophoretic and Brownian diffusion are two most important phenomena considered in his model [19]. In contrast to Buongiorno’s model, the Tiwari–Das model uses nanoparticle volumetric fraction, whereas the two-phase models use physical phenomena such as gravity, interphase friction, sedimentation, Brownian diffusion, and dispersion. However, these approaches do not consider sufficient effects such as thermophoresis [20]. In this context, Tsai [21] considers mechanisms such as electrophoresis and convection in addition to Brownian and thermophoretic diffusion as the major primary nanoparticle deposition mechanisms. Therefore, Buongiorno’s model has been used in various studies to model nanoparticle behavior. In particular, Kuznetsov and Nield [22] studied the nanoparticles’ impact on a natural convection flow on a vertical plate, while Hwang et al. [23] studied their impact in a laminar forced-convection flow. Nield and Kuznetsov [24, 25], Malvandi et al. [26], Sheikhzadeh et al. [27], and Saleem et al. [28] performed extensive studies on nanofluid convective transport. All these studies highlight the adequacy of Buongiorno’s model.

Mala et al. [29] experimentally confirmed the significance of the electrokinetic effects and other factors such as channel material and size on the characteristics of microchannel flow. Electroosmosis is an electrokinetic effect defined as the motion of the ionized liquid in relation to stationary charged surfaces, and one of the most important mechanisms due to this motion is the electrical double layer. Several theoretical and experimental studies investigating the electroosmotic effects on flow behavior can be found in the literature [30, 31]. In particular, Hunter [32] observed that electrical double layer (EDL) near the microchannel walls are generated by mutual attraction between free ions in the fluid and the static ions in the channel boundaries. Since most of the MEMS operate in the presence of electrical double layer, it is essential to investigate their influence on flow and heat transfer behavior. Mala et al. [33] modeled this behavior with the assumption that the electrostatic potential vanishes in the middle part of the channel. However, this assumption is only valid when the height between channel walls is very large. Ren et al. [34] made an improvement to this model with the symmetric flow pattern being assumed. From then on, many researchers [35, 36] considered various aspects of microchannel flow behavior.

In theoretical and experimental studies on microchannels, the pressure term was commonly prescribed as a constant. This is true in real circumstance. However, this term is closely associated with other physical quantities such as velocity, temperature, and fluid properties. Once it is prescribed, several other parameters are automatically determined. For example, in the well-known Poiseuille flow problem, the constant pressure is closely associated with the Reynolds number; we can determine the flow pattern either by the constant pressure or by the Reynolds number when the fluid properties remain unchanged. However, most researchers were not aware of this issue, which caused their models to be incompatible. Recently, Xu et al. [37] overcame this incompatibility in their model by assuming that the Reynolds number is known while the pressure constant is unknown. Their proposed model was further testified in a series of alternative microchannel configurations by Xu and his collaborators [38, 39].

There has been growing interest in optimization of system entropy in order to achieve better performance. According to the second law of thermodynamics, every real process is irreversible. However, this amount of irreversibility can be evaluated through entropy generation rate of the system. It is known that the quality of energy is lower for a higher entropy system [40]. Bejan [41] suggested that heat transfer and viscous friction are the major factors of entropy generation in a forced-convection viscous flow. Many researchers [42, 43] made efforts to understand the mechanism of entropy generation. Apart from the factors of heat transfer and viscous friction, other parameters such as Joule heating, electromagnetic diffusion, magnetic field, and nanoparticle distribution could be significant to affect energy generation in a microchannel system. The corresponding analysis was performed by Xie and Jian [44] who considered six factors that include heat transfer, viscous friction, Joule heating, electromagnetic diffusion, magnetic field, and nanoparticle distribution to analyze the fluid irreversibility rate in a microchannel with both magnetic and electrical fields being taken into account simultaneously. Xie and Jian considered a two-layer magnetohydrodynamic electroosmotic microchannel flow and analyzed the entropy of the system. A similar research was conducted by Escandón et al. [45], who studied entropy generation of non-Newtonian fluids affected by electroosmotic effects in a microchannel. Earlier in 2010, Zhao and Liu [46] studied evolution of entropy generation in open-end and closed-end microchannel flows under influence of electroosmosis. Some recent studies on entropy analysis are from Hayat et al. [47], Rashid et al. [48], and Zeeshan et al. [49].

Motivated from the aforementioned research, this work is an investigation of entropy generation in a fully developed nanofluid flow in a microradius circular pipe in the presence of electrical effects from an EDL. To the best of the authors’ knowledge, entropy analyses of electroosmotic nanofluid flow in a microradius pipe are rare in the literature. However, this research is necessary due to wide range of applications of MEMS and the importance of microradius structures in the design of comparable heat sinks. Buongiorno’s nanofluid model is used to describe the nanofluid flow behavior. The pressure term is assumed as an unknown constant as proposed by Xu et al. [37]. The governing equations embodying the conservation of electrical potential, total mass, momentum, energy, and nanoparticle volume fraction are reduced to a system of ordinary differential equations through appropriate similarity transformations. Highly accurate solutions are obtained by the improved homotopy analysis method (HAM) [50]. The physical behaviors of the flow are presented and discussed. The dimensionless parameters’ effects on the system entropy are analyzed in order to optimize system performance.

2. Mathematical Modeling

Consider a laminar, fully developed flow of a nanofluid driven in a microradius pipe by a constant pressure gradient in the presence of an electrical double layer. The pipe is placed along the z-axis with its radius being R and its length being L. It is assumed that the length is large enough so that the effects from the openings at the end of the pipe can be eluded. The physical sketch is presented in Figure 1 with the cylindrical coordinate system being used.

Because of the small radius of the pipe, parallel flow hypothesis can be assumed. The flow is assumed to be symmetrical in angular direction, i.e., . Since the zeta potential is constant throughout the pipe surface, the electrical potential depends only on the radial component, i.e., (refer to [51]). It shall be noted that the velocities in r-, ϕ-directions, based on the parallel flow hypothesis [52] (i.e., the length L is far greater than the radius R of the pipe), are negligible, i.e., and . Hence, only the z-direction velocity component is considered. The flow is thermally fully developed with a constant wall heat flux which implies that the temperature gradient in z-direction is constant, i.e., (refer to [53]). Furthermore, the nanoparticle mass flux on the wall is assumed to be constant, i.e., . The induced electric field vector is , while the flow is considered in the absence of external electric and magnetic fields. Since the flow is laminar, inertia effect, diffusiophoresis, Magnus effect, fluid drainage, and gravity effect can be neglected when considering Buongiorno’s model [18].

With those assumptions, the continuity equation is automatically satisfied. Hence, the governing equations embodying the conservation of the electrostatic potential, total mass, momentum, energy, and nanoparticle volume fraction are written aswhere ψ is the electric potential, is the velocity in z-direction, T is the temperature, C is the nanoparticle volume fraction, is the streaming potential in z-direction, is the thermal diffusivity, is the heat capacity ratio, ρ is the fluid density, c is the specific heat, is the thermal conductivity with the subscripts f and p denoting the quantity for the fluid or solid particles, respectively, is the charge density, ε is the fluid medium permittivity, p is the pressure, μ is the viscosity, σ is the electrical conductivity, is the reference temperature, is the coefficient of Brownian diffusion, and is the coefficient of thermophoretic diffusion.

The flow is assumed to satisfy the no-slip condition, and the zeta potential, the temperature, and the nanoparticle volume fraction on the boundary are assumed to be constant, given as , , and , respectively. So, the appropriate boundary conditions for equations (1)–(4) arealong with the average velocity condition due to the conservation of flow flux:

With a uniform dielectric constant and the assumption that the thermal energy is significantly larger than the electrical potential, the Poisson–Boltzmann equation (1) is reduced to [54]where is the Debye–Hückel parameter, whose reciprocal determines the width of EDL. Here, is the bulk concentration of ions, is the ion valence, e is the fundamental charge, is the Boltzmann constant, and is the absolute temperature.

3. Nondimensional Reduction

The following similarity variables are defined to reduce the governing equations to a nondimensional form:

It shall be noted that the domain [0, R] is transformed to [0, 1] after application of transformations (8), and the respective governing equations becomewhich are subject to the following boundary conditions:and the average velocity condition becomes

Here, represents the dimensionless zeta potential and κ, , , , , , , and denote the nondimensional parameters for the electroosmosis, i.e., the pressure gradient, the streaming potential, the thermophoretic diffusion, the Brownian diffusion, the Brinkman number, the constant temperature, and the nanoparticle concentration gradients in z-direction, respectively, which are defined as

The nondimensional streaming potential parameter in equation (11) is unknown, which is derived by using the equilibrium condition of the streaming current produced from the movement of charged particles in the diffuse layer and the conduction current produced by the streaming potential generated from the streaming current (refer to [51]), which is obtained aswhere G5 = LW_mn₀z^e/(ζλ₀).

3.1. Entropy Analysis

Entropy of the system can be computed after the solutions for the system of equations (9)–(14) are obtained. According to Bejan’s conclusion [42], the entropy generation rate is given bywhere , , , and are the entropy terms for the reversibility due to heat transfer, Joule heating, viscous friction, and the nanoparticle distribution, respectively, which are written as

The similarity variables defined in (8) are used to obtain dimensionless form of equation (18). The characteristic entropy rate is given as

In the above expression, , , and are given as

Finally, the total entropy in the pipe can be computed by taking an integral over the pipe radius:

4. Solution Procedure

Solutions of equations (9)–(14) are obtained using the homotopy analysis method (HAM) with the help of Mathematica computational software. The first step of solving a system of differential equations by the HAM is the formulation of the deformation equation which is given aswhere represents a mapping function of , denotes the chosen linear operator, denotes the nonlinear operator, the embedding parameter , and denotes the convergence control parameter. From the deformation equation, the mth-order HAM deformation equation can be obtained by taking m-times derivative with respect to q and dividing the resulting equation by , then finally letting . So, from equation (23), it becomeswhere

In order to solve the problem in the HAM framework, the unknown solutions for , W, θ, and s are assumed to be of the following form:and the unknown pressure gradient constant is assumed to be

The linear operator is chosen to be

With the assumed solution series expressions (26) and the linear operator (28), the mth-order HAM deformation equations can be written assubject to the boundary conditionswith the following additional condition:

The initial guesses are chosen as

In equation (29), , , , and denote the convergence control parameters, whereas , , , and arewhere . The solutions of equation (29) are expressed as follows:where , , , and are particular solutions, defined byand is the inverse linear operator which holds the following properties:

The constants , , defined in equations (34)–(37) can be obtained using the boundary conditions (30). It shall be noted that the streaming potential parameter and the pressure gradient constant are unknowns in equations (10) and (11) which are computed during the HAM computations. In order to compute , the additional condition (15) is used on each iteration, and as a result, the series approximation for is obtained. is obtained in a similar manner, but since the series approximation is not needed, the solution from the highest order computation is used in equation (17).

5. Results

This section is divided into two parts. In the first part, the HAM solution is verified using the respective errors from the equations and comparing the electric potential obtained through the analytical and HAM approach. The latter part is a discussion on the physical aspects of the obtained solutions.

5.1. Solution Validation

The convergence of the HAM technique has been tested and is presented in Figure 2. The 20th-order HAM solutions are used to determine the convergence domains for the parameters. It shall be noted that, for every convergence control parameter, there exists a convergence domain. Therefore, it can be stated that the applied technique is convergent. For the validation of the obtained solution, the approach used by [54] is used. Accuracy of the solution is checked using the maximum total average squared error which is given aswhere m represents the order of the HAM computation and denote the functions for error evaluation, which are defined, based on equations (9)–(12), given as

Equation (40) is used to compute the mth-order HAM computation error. Table 1 illustrates the errors computed using the following parameter values: , , , , , and . Note in the table that the error decreases as the HAM computational order is increased; this shows that the results are convergent. It should be noted that the homotopy-Padé technique [50] was used to increase the convergence rate and obtain more accuracy, indicating the insignificance of the convergence control parameters which can hence be given any nonzero value. The 50th-order HAM solutions have a good enough accuracy to be used for the foregoing analysis. The HAM approach for solution computation can be further verified using the analytical solution obtained for equation (9). After application of the boundary conditions (15) and (16) and neglecting the imaginary parts, the solution is obtained aswhere is the zeroth-order modified Bessel function of the first kind.

A comparison of the HAM exact solutions and the analytical solutions for the electric potential field is shown in Figure 3(a). Note that the solutions obtained through both techniques are smooth throughout the domain. Furthermore, the exact solution calculated through the HAM technique is in complete agreement with the analytical solution from equation (42). Similarly, Figure 3(b) presents a comparison of the HAM velocity solution with the analytical solutions obtained by Sun et al. [55], which is given as

(a)

(b)

After using the boundary conditions (13) and (14), equation (43) becomes

The values for unknown constants and obtained through the HAM formulation are used with . It can be seen that the 50th-order HAM solution is fully consistent with the analytical solution.

5.2. Analysis and Discussion

The electroosmotic parameter κ depends on the pipe radius and Debye–Hückel parameter. An increase in this parameter implies an increase in the channel radius and a decrease in the width of EDL compact layer. Because of this change, the electrostatic potential in the pipe is expected to decrease. Additionally, EDL has minimal effects in channels with larger radius. Therefore, the electrostatic potential is expected to be immune to a change in κ when this parameter is significantly large. It can be seen in Figure 3(a) that the electrostatic potential decreases in the mid part of the pipe as κ increases. It can also be seen that the effect on electrostatic potential decreases as κ increases and it becomes negligible when κ in increased to a sufficiently large value. This trend is consistent with the study conducted by Mala et al. [33], which analyzed the Boltzmann–Poisson model for EDL.

The two unknown constants, streaming potential and pressure gradient , are computed during the HAM computations. Figure 4 illustrates the trend in these quantities with respect to the electroosmotic parameter κ. Gong and Wu [56] report the streaming potential as the electric resistance number and state that the larger value of this constant means higher resistance due to EDL. Note that, as κ increases, the streaming potential decreases; but, this decrease becomes negligible after a significantly large value of κ as the approaches to zero. The effects of this decrease can be seen in the velocity and temperature profiles hereinafter. It shows that, as the EDL thickness decreases, the EDL effects decrease and become negligible eventually. In Figure 4(b), it can be seen that the pressure gradient constant starts decreasing rapidly from 12; but, as κ increases to a significantly large value, this change becomes negligible and eventually becomes constant. Physically, κ is inversely proportional to the thickness of EDL compact layer; so, after a significantly large value of κ, the radius does not change enough to influence the pressure gradient constant.

(a)

(b)

Malvandi and Ganj [57] noted that the nanoparticle concentration in the mid part of the channel decreased with increase in thermophoretic parameter. They drew a similar conclusion in case of microtubes in another study [59]. The particle acceleration is higher for higher thermophoretic diffusion which causes a drop in the fluid nanoparticle concentration. Figure 5 illustrates the change in the nanoparticle volume fraction s(η) when the thermophoretic diffusion parameter and the Brinkman number are varied. Note in Figure 5(a) that increasing the thermophoretic diffusion in the flow causes a decrease in the nanoparticle concentration throughout the domain. Figure 5(b) shows change in the nanoparticle concentration when the Brinkman number is varied. Note that the nanoparticle concentration decreases with increase in Br due to the increased fluid viscosity.

(a)

(b)

The electrostatic potential is expected to affect the flow velocity in such a way that the resistance due to EDL will initially increase near the pipe boundary. However, this resistance will start to decrease as κ exceeds the max effect value . This behavior can be observed in Figure 6. The dimensionless skin friction coefficient is given as , where is the Reynolds number. It can be seen in Figure 6(a) that the resistance on the pipe boundary initially increases with increase in κ, but this trends changes as κ exceeds . is obtained through maximization of with respect to κ.

(a)

(b)

The velocity is observed against η with variation in κ in Figure 6(b). Note that there is a small change in the velocity profile when κ is changed. This change is maximum at and smaller for larger values of κ, which indicates the decrease of EDL effects with increase in the pipe radius; the velocity near the boundary decreases with increase in κ, but this decrease gets less prominent for higher values of κ and eventually, the velocity profile approaches to the case with no EDL. This is consistent with the results obtained by Gong and Wu [56], according to which the resistance due to EDL is greater near the boundaries as compared to the mid part of the channel. Additionally, according to Zhao et al. [54], the resistance in the velocity due to EDL is reduced with increase in κ, which is also evident in this case.

Similarly, the electrostatic potential profile will affect the temperature distribution due to Joule heating effect in the energy equation as illustrated in Figure 7. The dimensionless Nusselt number, given as , determines the dominance of convective heat transfer through the pipe/channel boundary. Plotted in Figure 7(a), the Nusselt number increases with increase in κ. However, this increase sufficiently slows down as and eventually becomes negligible when κ is sufficiently large. The temperature profile with change in κ is observed in Figure 7(b). Note that the temperature throughout the pipe decreases as κ increases. Since the EDL effects decrease with increase in κ, the effects from this term also decrease as a consequence, and this change becomes negligibly small as κ increases to a sufficiently large value.

(a)

(b)

It is important to analyze the influence of the thermophoretic diffusion, Brownian diffusion, and viscosity on the temperature profile in a microradius flow. An increase in the thermophoretic diffusion in the flow can cause an increase in the fluid-particle acceleration, which as a consequence can cause an increase in the flow temperature profile. This trend was observed by Malvandi and Ganj [57] in their study of water/alumina nanofluid in a vertical microchannel. Their results state that when the thermophoretic parameter is increased, the temperature in the channel is increased. Similarly, there is a higher amount of molecular collision, which causes rise in the flow temperature, when Brownian diffusion in a flow is increased. This was observed by Animasaun et al. [58] in detail. According to their hypothesis, the flow temperature is bound to increase with an increase in the Brownian motion in the flow. The Brinkman number (Br) is the ratio of the heat generation due to viscosity to external heating. So, an increase in the Brinkman number implies an increase in the fluid viscosity, which can raise the temperature in the fluid. It shall be noted in Figure 8(a) that, as the thermophoretic parameter increases, the temperature profile increases throughout the domain. Similarly, the temperature distribution enhances with an increase in the Brownian diffusion parameter, which is shown in Figure 8(b). Similarly, the temperature increases with increase in the Brinkman number, as shown in Figure 8(c). All these results are consistent with the previous studies on heat convection in fluid flow.

(a)

(b)

(c)

Figures 9 and 10 present an entropy analysis for the flow. Note in equation (20) that the flow entropy is highly dependent on the temperature, velocity, and nanoparticle concentration of the flow. From the former analysis, it is obvious that the electroosmotic parameter κ effects the flow velocity and temperature fields. The effects on the entropy due to change in κ are illustrated in Figure 9. Note that, in Figure 9(a), the total entropy of the flow is maximum when the value of κ is near zero. The total flow entropy decreases with increase in κ. From previous analysis, it is evident that the EDL effect is maximum at . This effect starts to decrease as , which is evident in the flow entropy trend shown in the figure. The observed decrease is due to the decrease in temperature profile with increase in κ, which can be further understood by observing Figure 9(b), where the entropy field of the flow is observed at different values of κ. It can be seen that, for lower values of κ, i.e., , there is more disturbance in the entropy near the pipe wall and the entropy is relatively higher. It reduces throughout the channel when κ is increased, but after a sufficiently large value of κ, this change becomes negligible.

(a)

(b)

The effects of the thermophoretic diffusion parameter () and the Brinkman number () on the total entropy of the flow can be seen in Figure 10. The total entropy of the channel increases with increase in Br, whereas this trend is reversed in case of where it decreases with increase in the thermophoretic parameter. The reason for this trend is illustrated in Table 2. It shall be noted that, with increase in Br, there is a notable increase in the entropy due to heat transfer and viscous friction, whereas a minute decrease can be observed due to the nanoparticle and Joule heating. Hence, the total entropy of the flow is increased. On the other hand, entropy due to nanoparticle has a prominent decrease when is increased; this causes a decrease in the total entropy of the flow. This can also be seen in Figures 5(a) and 8(a) where the decrease in the nanoparticle concentration (due to ) is greater compared to the increase in the temperature field.

6. Conclusions

A fully developed nanofluid flow in a microradius pipe with the presence of electrical effects from an electrical double layer has been studied. The governing equations embodying, the Poisson–Boltzmann equation, momentum equation, energy equation, and Buongiorno’s nanofluid model equation, have been transformed using the appropriate similarity transformations. The unknown pressure gradient constant has been computed using the average velocity condition, and the streaming potential due to EDL has derived for the considered problem. Highly accurate exact solutions have been obtained using the homotopy analysis method. The findings of this study are as follows.

The model used for formulation of the nanofluid flow is consistent with the experimental studies in the past. Change in the electroosmotic parameter has high influence on the flow when it is smaller than the max effect value , which was found to be 2.30675. The EDL effects reduce and eventually become negligible as . The electrostatic potential, velocity, temperature, and nanoparticle volume fraction are affected by the evolution in viscosity, thermophoretic diffusion, Brownian diffusion, and electroosmosis. The flow entropy increases with increase in viscosity and thermophoretic diffusion, while it decreases when the EDL compact layer is reduced. Furthermore, the nanoparticle and temperature concentration are found to be two of the main factors causing the change in flow entropy.

Nomenclature

:	Brinkman number
C:	Nanoparticle volume fraction
:	Reference nanoparticle volume fraction
:	Nanoparticle volume fraction on the pipe surface
:	Local skin friction coefficient
:	Specific heat of fluid and nanoparticles
:	Brownian diffusion coefficient ()
:	Thermophoretic diffusion coefficient ()
e:	Charge of a proton (C)
:	Error for HAM computation order m
:	Dimensionless streaming potential parameter
:	Electric field in z direction ()
k:	Debye–Hückel parameter ()
:	Boltzmann constant (J·)
:	Thermal conductivity of the fluid (W·)
L:	Length of the pipe (m)
:	Dimensionless mass diffusion parameter
:	The bulk ionic concentration ()
:	Dimensionless Brownian motion parameter
:	Dimensionless thermophoresis parameter
:	Local Nusselt number
:	Pressure (Pa)
:	Heat flux on the pipe surface (W·)
R:	Radius of the pipe (m)
:	Universal gas constant
:	Reynolds number
s:	Nondimensional nanoparticle volume fraction
:	Entropy concentration in the pipe
:	Total entropy generated in the pipe
T:	Temperature (K)
:	Reference temperature (K)
:	Temperature on the pipe surface (K)
:	Absolute temperature (K)
W:	Nondimensional velocity of the fluid
:	Average velocity of the fluid ()
:	z-component of the fluid velocity ()
:	The valence of ions

Greek symbols

α:	Thermal diffusivity of the nanofluid ()
ε:	Dielectric constant of the medium
:	Permittivity of vacuum (C·)
η:	Nondimensional spatial variable
:	Nondimensional pressure gradient parameter
κ:	Electroosmotic parameter
:	Max effect electroosmotic parameter value
:	Dimensionless reference nanoparticle parameter
μ:	Dynamic viscosity of the fluid (kg·)
:	Dimensionless reference temperature parameter
θ:	Dimensionless temperature distribution
:	Density of the fluid and nanoparticles (kg·)
:	Charge density (C·)
:	Shear stress on the surface of the pipe (Pa)
ψ:	Electrostatic potential (V)
:	Dimensionless electrostatic potential
ζ:	Dimensionless zeta potential
:	Zeta potential

Subscripts

:	The fluid and solid particles
:	Physical quantities on the pipe surface

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was partially supported by the National Natural Science Foundation of China (grant no. 11872241).

References

A. Majumdar, J. P. Carrejo, and J. Lai, “Thermal imaging using the atomic force microscope,” Applied Physics Letters, vol. 62, no. 20, pp. 2501–2503, 1993.
View at: Publisher Site | Google Scholar
D. B. Tuckerman and R. F. W. Pease, “High-performance heat sinking for VLSI,” IEEE Electron Device Letters, vol. 2, no. 5, pp. 126–129, 1981.
View at: Publisher Site | Google Scholar
G. M. Mala and D. Li, “Flow characteristics of water in microtubes,” International Journal of Heat and Fluid Flow, vol. 20, no. 2, pp. 142–148, 1999.
View at: Publisher Site | Google Scholar
J. Darabi and K. Ekula, “Development of a chip-integrated micro cooling device,” Microelectronics Journal, vol. 34, no. 11, pp. 1067–1074, 2003.
View at: Publisher Site | Google Scholar
R. J. Phillips, “Microchannel heat sinks,” The Lincoln Laboratory Journal, vol. 1, no. 1, pp. 31–48, 1988.
View at: Google Scholar
S. U. Choi and J. A. Eastman, “Enhancing thermal conductivity of fluids with nanoparticles,” in Proceedings of the 1995 International Mechanical Engineering Congress and Exhibition, San Francisco, CA, USA, October 1995.
View at: Google Scholar
A. Bibi, H. Xu, Q. Sun, I. Pop, and Q. Zhao, “Free convection of a hybrid nanofluid past a vertical plate embedded in a porous medium with anisotropic permeability,” International Journal of Numerical Methods for Heat & Fluid Flow.
View at: Publisher Site | Google Scholar
M. Waqas, S. A. Shehzad, T. Hayat, M. I. Khan, and A. Alsaedi, “Simulation of magnetohydrodynamics and radiative heat transport in convectively heated stratified flow of Jeffrey nanofluid,” Journal of Physics and Chemistry of Solids, vol. 133, pp. 45–51, 2019.
View at: Publisher Site | Google Scholar
M. I. Khan, M. Waqas, T. Hayat, A. Alsaedi, and M. I. Khan, “Significance of nonlinear radiation in mixed convection flow of magneto Walter-B nanoliquid,” International Journal of Hydrogen Energy, vol. 42, no. 42, pp. 26408–26416, 2017.
View at: Publisher Site | Google Scholar
D. Yadav and J. Wang, “Convective heat transport in a heat generating porous layer saturated by a non-Newtonian nanofluid,” Heat Transfer Engineering, vol. 40, no. 16, pp. 1363–1382, 2019.
View at: Publisher Site | Google Scholar
B. Mahanthesh, B. J. Gireesha, I. L. Animasaun, T. Muhammad, and N. S. Shashikumar, “MHD flow of SWCNT and MWCNT nanoliquids past a rotating stretchable disk with thermal and exponential space dependent heat source,” Physica Scripta, vol. 94, no. 8, Article ID 085214, 2019.
View at: Publisher Site | Google Scholar
I. L. Animasaun, B. Mahanthesh, and O. K. Koriko, “On the motion of non-Newtonian eyring–powell fluid conveying tiny gold particles due to generalized surface slip velocity and buoyancy,” International Journal of Applied and Computational Mathematics, vol. 4, no. 6, p. 137, 2018.
View at: Publisher Site | Google Scholar
S. P. Jang and S. U. S. Choi, “Role of Brownian motion in the enhanced thermal conductivity of nanofluids,” Applied Physics Letters, vol. 84, no. 21, pp. 4316–4318, 2004.
View at: Publisher Site | Google Scholar
A. Behzadmehr, M. Saffar-Avval, and N. Galanis, “Prediction of turbulent forced convection of a nanofluid in a tube with uniform heat flux using a two phase approach,” International Journal of Heat and Fluid Flow, vol. 28, no. 2, pp. 211–219, 2007.
View at: Publisher Site | Google Scholar
H. Aminfar, M. Mohammadpourfard, and F. Mohseni, “Two-phase mixture model simulation of the hydro-thermal behavior of an electrical conductive ferrofluid in the presence of magnetic fields,” Journal of Magnetism and Magnetic Materials, vol. 324, no. 5, pp. 830–842, 2012.
View at: Publisher Site | Google Scholar
A. Zeeshan, M. M. Bhatti, N. Ijaz, O. A. Beg, and A. Kadir, “Biologically-inspired transport of solid spherical nanoparticles in an electrically-conducting viscoelastic fluid with heat transfer,” Thermal Science, p. 324, 2018.
View at: Publisher Site | Google Scholar
R. K. Tiwari and M. K. Das, “Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids,” International Journal of Heat and Mass Transfer, vol. 50, no. 9-10, pp. 2002–2018, 2007.
View at: Publisher Site | Google Scholar
J. Buongiorno, “Convective transport in nanofluids,” Journal of Heat Transfer, vol. 128, no. 3, pp. 240–250, 2006.
View at: Publisher Site | Google Scholar
O. D. Makinde and I. L. Animasaun, “Thermophoresis and Brownian motion effects on MHD bioconvection of nanofluid with nonlinear thermal radiation and quartic chemical reaction past an upper horizontal surface of a paraboloid of revolution,” Journal of Molecular Liquids, vol. 221, pp. 733–743, 2016.
View at: Publisher Site | Google Scholar
S. Mansur and A. Ishak, “The flow and heat transfer of a nanofluid past a stretching/shrinking sheet with a convective boundary condition,” Abstract and Applied Analysis, vol. 2013, Article ID 350647, 9 pages, 2013.
View at: Publisher Site | Google Scholar
R. Tsai, “A simple approach for evaluating the effect of wall suction and thermophoresis on aerosol particle deposition from a laminar flow over a flat plate,” International Communications in Heat and Mass Transfer, vol. 26, no. 2, pp. 249–257, 1999.
View at: Publisher Site | Google Scholar
A. V. Kuznetsov and D. A. Nield, “Natural convective boundary-layer flow of a nanofluid past a vertical plate,” International Journal of Thermal Sciences, vol. 49, no. 2, pp. 243–247, 2010.
View at: Publisher Site | Google Scholar
K. S. Hwang, S. P. Jang, and S. U. Choi, “Flow and convective heat transfer characteristics of water-based Al₂O₃ nanofluids in fully developed laminar flow regime,” International Journal of Heat and Mass Transfer, vol. 52, no. 1-2, pp. 193–199, 2009.
View at: Publisher Site | Google Scholar
D. Nield and A. Kuznetsov, “The Cheng–Minkowycz problem for natural convective boundary-layer flow in a porous medium saturated by a nanofluid,” International Journal of Heat and Mass Transfer, vol. 52, no. 25-26, pp. 5792–5795, 2009.
View at: Publisher Site | Google Scholar
D. Nield and A. Kuznetsov, “The Cheng–Minkowycz problem for the double-diffusive natural convective boundary layer flow in a porous medium saturated by a nanofluid,” International Journal of Heat and Mass Transfer, vol. 54, no. 1–3, pp. 374–378, 2011.
View at: Publisher Site | Google Scholar
A. Malvandi, D. D. Ganji, F. Hedayati, and E. Yousefi Rad, “An analytical study on entropy generation of nanofluids over a flat plate,” Alexandria Engineering Journal, vol. 52, no. 4, pp. 595–604, 2013.
View at: Publisher Site | Google Scholar
G. A. Sheikhzadeh, M. Dastmalchi, and H. Khorasanizadeh, “Effects of nanoparticles transport mechanisms on Al₂O₃-water nanofluid natural convection in a square enclosure,” International Journal of Thermal Sciences, vol. 66, pp. 51–62, 2013.
View at: Publisher Site | Google Scholar
S. Saleem, H. Rafiq, A. Al-Qahtani, M. A. El-Aziz, M. Y. Malik, and I. L. Animasaun, “Magneto Jeffrey nanofluid bioconvection over a rotating vertical cone due to gyrotactic microorganism,” Mathematical Problems in Engineering, vol. 2019, Article ID 3478037, 11 pages, 2019.
View at: Publisher Site | Google Scholar
G. Mala, D. Li, C. Werner, H.-J. Jacobasch, and Y. Ning, “Flow characteristics of water through a microchannel between two parallel plates with electrokinetic effects,” International Journal of Heat and Fluid Flow, vol. 18, no. 5, pp. 489–496, 1997.
View at: Publisher Site | Google Scholar
D. Yadav, “Electrohydrodynamic instability in a heat generating porous layer saturated by a dielectric nanofluid,” Journal of Applied Fluid Mechanics, vol. 10, no. 3, pp. 763–776, 2017.
View at: Publisher Site | Google Scholar
D. Yadav, “The influence of pulsating throughflow on the onset of electro-thermo-convection in a horizontal porous medium saturated by a dielectric nanofluid,” Journal of Applied Fluid Mechanics, vol. 11, no. 6, pp. 1679–1689, 2018.
View at: Publisher Site | Google Scholar
R. J. Hunter, Zeta Potential in Colloid Science: Principles and Applications, Academic Press, Cambridge, MA, USA, 2013.
G. M. Mala, D. Li, and J. Dale, “Heat transfer and fluid flow in microchannels,” International Journal of Heat and Mass Transfer, vol. 40, no. 13, pp. 3079–3088, 1997.
View at: Google Scholar
L. Ren, W. Qu, and D. Li, “Interfacial electrokinetic effects on liquid flow in microchannels,” International Journal of Heat and Mass Transfer, vol. 44, no. 16, pp. 3125–3134, 2001.
View at: Publisher Site | Google Scholar
C. Qi and C.-O. Ng, “Electroosmotic flow of a two-layer fluid in a slit channel with gradually varying wall shape and zeta potential,” International Journal of Heat and Mass Transfer, vol. 119, pp. 52–64, 2018.
View at: Publisher Site | Google Scholar
J. Zheng and Y. Jian, “Rotating electroosmotic flow of two-layer fluids through a microparallel channel,” International Journal of Mechanical Sciences, vol. 136, pp. 293–302, 2018.
View at: Publisher Site | Google Scholar
H. Xu, I. Pop, and Q. Sun, “Fluid flow driven along microchannel by its upper stretching wall with electrokinetic effects,” Applied Mathematics and Mechanics, vol. 39, no. 3, pp. 395–408, 2018.
View at: Publisher Site | Google Scholar
H. Xu and H. Huang, “Liquid flow in a porous channel with electrokinetic effects,” Communications in Theoretical Physics, vol. 70, no. 4, p. 391, 2018.
View at: Publisher Site | Google Scholar
H. Xu, A. Raees, and X.-H. Xu, “Entropy generation of nanofluid flow and heat transfer driven through a paralleled microchannel,” Canadian Journal of Physics, vol. 97, no. 6, pp. 678–691, 2018.
View at: Google Scholar
Y. A. Cengel and M. A. Boles, “Thermodynamics: an engineering approach,” Sea, vol. 1000, p. 8862, 2002.
View at: Google Scholar
A. Bejan, “A study of entropy generation in fundamental convective heat transfer,” Journal of Heat Transfer, vol. 101, no. 4, pp. 718–725, 1979.
View at: Publisher Site | Google Scholar
A. Bejan, “Second law analysis in heat transfer,” Energy, vol. 5, no. 8-9, pp. 720–732, 1980.
View at: Publisher Site | Google Scholar
O. D. Makinde, “Entropy-generation analysis for variable-viscosity channel flow with non-uniform wall temperature,” Applied Energy, vol. 85, no. 5, pp. 384–393, 2008.
View at: Publisher Site | Google Scholar
Z.-Y. Xie and Y.-J. Jian, “Entropy generation of two-layer magnetohydrodynamic electroosmotic flow through microparallel channels,” Energy, vol. 139, pp. 1080–1093, 2017.
View at: Publisher Site | Google Scholar
J. Escandón, O. Bautista, and F. Méndez, “Entropy generation in purely electroosmotic flows of non-Newtonian fluids in a microchannel,” Energy, vol. 55, pp. 486–496, 2013.
View at: Publisher Site | Google Scholar
L. Zhao and L. H. Liu, “Entropy generation analysis of electro-osmotic flow in open-end and closed-end micro-channels,” International Journal of Thermal Sciences, vol. 49, no. 2, pp. 418–427, 2010.
View at: Publisher Site | Google Scholar
T. Hayat, M. I. Khan, S. Qayyum, and A. Alsaedi, “Entropy generation in flow with silver and copper nanoparticles,” Colloids and Surfaces A: Physicochemical and Engineering Aspects, vol. 539, pp. 335–346, 2018.
View at: Publisher Site | Google Scholar
M. Rashid, M. I. Khan, T. Hayat, M. I. Khan, and A. Alsaedi, “Entropy generation in flow of ferromagnetic liquid with nonlinear radiation and slip condition,” Journal of Molecular Liquids, vol. 276, pp. 441–452, 2019.
View at: Publisher Site | Google Scholar
A. Zeeshan, N. Shehzad, T. Abbas, and R. Ellahi, “Effects of radiative electro-magnetohydrodynamics diminishing internal energy of pressure-driven flow of titanium dioxide-water nanofluid due to entropy generation,” Entropy, vol. 21, no. 3, p. 236, 2019.
View at: Publisher Site | Google Scholar
S. Liao, Homotopy Analysis Method in Nonlinear Differential Equations, Higher Education Press, Beijing, China, 2012.
R. J. Hunter, “Zeta potential,” in Colloid Science, Academic Press, London, UK, 1981.
View at: Google Scholar
R. Chakraborty, R. Dey, and S. Chakraborty, “Thermal characteristics of electromagnetohydrodynamic flows in narrow channels with viscous dissipation and Joule heating under constant wall heat flux,” International Journal of Heat and Mass Transfer, vol. 67, pp. 1151–1162, 2013.
View at: Publisher Site | Google Scholar
Z.-Y. Xie, Y.-J. Jian, and F.-Q. Li, “Thermal transport of magnetohydrodynamic electroosmotic flow in circular cylindrical microchannels,” International Journal of Heat and Mass Transfer, vol. 119, pp. 355–364, 2018.
View at: Publisher Site | Google Scholar
Q. Zhao, H. Xu, and L. Tao, “Nanofluid flow and heat transfer in a microchannel with interfacial electrokinetic effects,” International Journal of Heat and Mass Transfer, vol. 124, pp. 158–167, 2018.
View at: Publisher Site | Google Scholar
L.-X. Sun, Y.-J. Jian, L. Chang, H.-Y. Zhang, and Q.-S. Liu, “Alternating current electro-osmotic flow of the Maxwell fluids through a circular micro-pipe,” Journal of Mechanics, vol. 29, no. 2, pp. 233–240, 2013.
View at: Publisher Site | Google Scholar
L. Gong and J.-K. Wu, “Resistance effect of electric double layer on liquid flow in microchannel,” Applied Mathematics and Mechanics, vol. 27, no. 10, pp. 1391–1398, 2006.
View at: Publisher Site | Google Scholar
A. Malvandi and D. D. Ganji, “Mixed convective heat transfer of water/alumina nanofluid inside a vertical microchannel,” Powder Technology, vol. 263, pp. 37–44, 2014.
View at: Publisher Site | Google Scholar
I. L. Animasaun, R. O. Ibraheem, B. Mahanthesh, and H. A. Babatunde, “A meta-analysis on the effects of haphazard motion of tiny/nano-sized particles on the dynamics and other physical properties of some fluids,” Chinese Journal of Physics, vol. 60, pp. 676–687, 2019.
View at: Publisher Site | Google Scholar
A. Malvandi and D. D. Ganji, “Magnetohydrodynamic mixed convective flow of Al₂O₃-water nanofluid inside a vertical microtube,” Journal of Magnetism and Magnetic Materials, vol. 369, pp. 132–141, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Muhammad D. K. Niazi and Hang Xu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

446

Downloads

741

Citations