Alexandria Engineering Journal (2016) xxx, xxx-xxx

HOSTED BY

Alexandria University Alexandria Engineering Journal

www.elsevier.com/locate/aej www.sciencedirect.com

ORIGINAL ARTICLE

Steady fully developed natural convection flow in a vertical annular microchannel having temperature dependent viscosity: An exact solution

Basant K. Jha, Babatunde Aina *, Zubairu Rilwanu

Department of Mathematics, Ahmadu Bello University, Zaria, Nigeria Received 9 October 2015; revised 21 February 2016; accepted 25 February 2016

KEYWORDS

Annular micro-channel; Natural convection; Temperature dependent viscosity;

Velocity slip and temperature jump

Abstract This study is devoted to investigate the steady fully developed natural convection flow in a vertical annular micro-channel having temperature dependent viscosity in the presence of velocity slip and temperature jump at the annular micro-channel surfaces. The governing equations of the motion are a set of ordinary differential equations and their analytical solutions in dimensionless form have been obtained for the temperature field and velocity field. The effect of various flow parameters entering into the problem is discussed with the aid of line graphs. During the course of numerical investigation, it is found that increase in viscosity variation parameter enhances the fluid velocity and velocity slip. Furthermore, an increase in viscosity variation parameter leads to increase in the volume flow rate and skin friction.

© 2016 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

1. Introduction

Fluid flow in micro-channel has continued to attract interest because of its practical applications in space systems, manufacturing and material processing operations, and in high-power-density chips in supercomputers and other electronics. Several investigations have been accomplished on forced convection fluid flow in micro-channels and micro-tube. However, only a few studies have been carried out on natural and mixed convection in vertical micro-channels and micro-tube. Chen and Weng [1] studied the flow mechanism in a vertical

* Corresponding author.

E-mail addresses: basant777@yahoo.co.uk (B.K. Jha), ainavicdydx@ gmail.com (B. Aina), zubairurilwanu72@gmail.com (Z. Rilwanu). Peer review under responsibility of Faculty of Engineering, Alexandria University.

micro-channel and obtained an exact solution of the fully developed natural convection in an open-ended vertical parallel-plate micro-channel due to asymmetric heating of micro-channel walls. They found that the rarefaction and fluid-wall interaction have significant effects on the flow and thermal fields. Jha et al. [2] extended this work by taking into account suction/injection on the micro-channel walls. They concluded in their work that skin friction as well as rate of heat transfer is strongly dependent on suction/injection parameter. The transient hydrodynamics and thermal behaviors of fluid flow in an open-ended vertical parallel-plate micro-channel, under the effect of the hyperbolic-heat-conduction model, were investigated semi-analytically in [3]. They concluded that, as Knudsen number increases, velocity slip and temperature jump increase at the boundaries. Haddad et al. [4] numerically investigated the developing hydrodynamical behaviors of free convection gas flow in a vertical open-ended parallel-plate

http://dx.doi.Org/10.1016/j.aej.2016.02.023

1110-0168 © 2016 Faculty of Engineering, Alexandria University. Production and hosting by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

Nomenclature

B viscosity variation parameter w dimensional gap between the cylinders

CPo specific heat at constant pressure thermal and tangential momentum accommoda-

ln fluid-wall interaction parameter, bt/ bm tion coefficients, respectively

g gravitational acceleration

ki radius of the inner cylinder Greek letters

k2 radius of the outer cylinder a thermal diffusivity

Kn Knudsen number, k/w bo coefficient of thermal expansion

q volume flow rate bt; bv dimensionless variables

Q dimensionless volume flow rate c ratio of specific heats

Pr Prandtl number l dynamic viscosity

r dimensional radial coordinate lo dynamic viscosity at T — T0

R dimensionless radial coordinate h dimensionless temperature

R specific gas constant Po density

T temperature of fluid V fluid kinematic viscosity (i0/p0)

To reference temperature g ratio of radii (ki/k2)

Ti temperature at outer surface of the inner cylinder k molecular mean free path

u axial velocity ko thermal conductivity

U dimensionless axial velocity s skin-friction

micro-channel filled with porous media. Their results showed that the slip in velocity and jump in temperature decreased in the axial direction of the flow. Numerical solutions were obtained by Buonomo and Manca [5] for natural convection in parallel-plate vertical microchannels due to asymmetric heating by imposing constant heat flux on the boundaries. Buonomo and Manca [6] further performed a numerical study on transient mixed convection in a vertical micro-channel due to asymmetric as well as symmetric heat fluxes on the microchannel surfaces. Chen and Weng [7] numerically studied the creep effect on the flow and heat transfer characteristics, for developing natural convective microflow in the same geometry.

In another related article, Weng and Chen [8] examined the influence of wall-surface curvature on the flow and thermal fields as well as the corresponding characteristics over the heated wall. The results of this work show the nonlinear behavior in temperature. It was also concluded that under certain rarefaction and wall interaction condition, by decreasing the curvature radius ratio, skin-friction decreases while rate of heat transfer increases. Jha et al. [9] extended the work of Weng and Chen [8] to the case in which the cylindrical surfaces forming the annulus are permeable, i.e. when there is a suction or injection through the annulus surfaces. They concluded that as suction/injection on the cylinder walls increases, the fluid velocity and temperature are enhanced. In another work, Jha and Aina [10] analyzed the work of Weng and Chen [8] by incorporating the pressure gradient in the vertical direction. It is observed that the probability of reverse flow formation increases at inner surface of outer cylinder of the micro-annulus by increasing curvature radius, while it decreases with increase in Knudsen number and fluid wall interaction parameter. Avci and Aydin [11,12] presented exact solutions for fully developed mixed convection in a vertical parallel-plate microchannel with constant plate temperature and constant heat flux on the plates respectively. Also, Avci and Aydin [13] studied the fully developed mixed convective heat transfer of a Newtonian

fluid in a vertical micro-annulus formed by two concentric micro-tubes. It is found that increasing mixed convection parameter enhances heat transfer while rarefaction effects considered by the velocity slip and the temperature jump in the slip flow regime decrease it. Jha and Aina [14] further extended the work of Avci and Aydin [13] to the case when suction/injection is imposed on the annulus surfaces. They concluded in their study that as suction/injection on the micro-porous-annulus (MPA) increases, the fluid velocity and temperature increase. Recently, Jha and Aina [15] investigated steady fully developed mixed convection flow in a vertical micro-annulus in the presence of transverse magnetic field. Das et al. [16] presented a theoretical analysis to investigate the effect of buoyancy force on mixed convective Couette flow of a reactive viscous incompressible nanofluid between two concentric cylindrical pipes under bimolecular, Arrhenius and sensitized reaction rates. In another article, Das et al. [17] studied the fully developed mixed convection flow in a vertical channel filled with nanoflu-ids in the presence of a uniform transverse magnetic field. They reported that the magnetic field tends to enhance the nanofluid velocity in the channel. Makinde [18] investigated the thermal analysis of a reactive generalized Couette flow of power law fluids between concentric cylindrical pipes.

On the other hand, most of the existing analytical studies for such problems are based on the constant physical properties of the fluid. However, accurate prediction for the flow formation and heat transfer can be achieved by considering variation of physical properties with temperature [19], especially for fluid viscosity. Makinde and Chinyoka [20] presented numerical solution of unsteady flow of a variable viscosity reactive fluid in a slit with wall suction/injection. Tshehla et al. [21] investigated the entropy generation rate in a variable viscosity liquid flowing steadily through a cylindrical pipe with convective cooling at the pipe surface. They reported that a decrease in the fluid viscosity and an increase in viscous heating enhance total entropy generation in the flow fluid. Klemp

et al. [22] studied numerically the effect of temperature dependent viscosity on the entrance flow in a channel in the hydro-dynamic case. The effects of variable viscosity on hydromagnetic flow and heat transfer have been studied by Seddeek [23]. Grasset and Parmentier [24] investigated thermal convection in a volumetrically heated fluid with strongly temperature dependent viscosity cooled from above. In another article, Saravanan and Kandaswamy [25] analyzed the hydro-magnetic stability of convective flow of variable viscosity fluids generated by internal heat sources. The effect of variable temperature dependent viscosity on the mixed convection flow from vertical plate is investigated by many authors, e.g. Hady et al. [26] and Mahmud [27]. The free-or-mixed convection boundary layer flow from a horizontal surface in a saturated porous medium taking into account the effect of variable viscosity has been studied by Kumari [28]. Recently, Umavathi and Ojjela [29] studied the effect of variable viscosity on free convection flow in a vertical rectangular duct. They observed that the negative values of viscosity variation parameter show intense velocity contour in the lower half region of the duct whereas positive values of viscosity variation parameter show the intense velocity contours in the upper half region of the duct. In another related article, Umavathi [30] studied analytically as well as numerically the combined effects of variable viscosity and variable thermal conductivity on double-diffusive convection flow of fluid in a vertical channel filled with porous material. However, derivation of any exact solution for steady fully developed natural convection flow in a vertical annular micro-channel with temperature dependent viscosity in the presence of velocity slip and temperature jump at the annular micro-channel surfaces with interfacial slip has not been attempted. It is well known that exact solutions have their own theoretical meaning, and many exact solutions played key roles in the early development of fluid mechanics and heat conduction [31,32]. Besides their theoretical importance, exact solutions can also be applied to checking the accuracy, convergence and effectiveness of various numerical computation methods and improving differencing schemes, grid generation ways and so on. Exact solutions are therefore very useful even for the newly rapidly developing computational fluid dynamics and heat transfer.

The objective of this work was to present exact solution of steady fully developed natural convection flow of fluid having temperature dependent viscosity in a vertical annular microchannel in the presence of velocity slip and temperature jump at the annular micro-channel surfaces. The mathematical model employed herein represents a generalization of the work discussed in [8] to include temperature dependent viscosity.

2. Mathematical analysis

Consider a steady fully developed natural convection flow in a vertical annular micro-channel with temperature dependent viscosity in the presence of velocity slip and temperature jump at the annular micro-channel surfaces. A schematic geometry of the problem under investigation is shown in Fig. 1, where X-axis is parallel to the gravitational acceleration g but in the opposite direction while the r-axis is in the radial direction. The radius of the inner and outer cylinder walls is k and k2, respectively. The outer surface of the inner cylinder is heated

Figure 1 Flow configuration and coordinate system.

to a temperature (Ti) greater than that of the surrounding fluid having temperature (T0) and the inner surface of the outer cylinder is maintained at temperature (T0). Due to this temperature difference, natural convection occurs in the vertical annular micro-channel. Since the flow is fully developed and cylinders are of infinite length, the flow depends only on radial coordinate (r). Using Boussinesq's approximation and considering temperature dependent viscosity, the governing momentum and energy equations, describing the present physical situation can be written in dimensional form as follows: Conservation of momentum

1 d du

-rdr^l dt) + gp0^T - T0 Conservation of energy

1 d (dT\ _ r dry dr I

The boundary conditions for the velocity and temperature field in dimensional form are as follows:

i(r — k1) — -

av . du ~kdr jr—ki

, j , 2 - av ,du,

U(r — k2) — -^~ k dr jr—k2

T(r — ki) — Ti + T(r — k2) — To -

2 — at 2y k dT

a, y + 1 Pr dr lr=ki

2 — a, 2y k dT ! a, y + 1 Pr dr

The mathematical model used in the present work to capture the viscosity variation with temperature is [33]

l — 1oexp[—b(T — To)]

where i0 is the viscosity when temperature is T0 while the coefficient b(0C-1) determine the strength of dependency between i and T.

Introducing the following dimensionless quantities in Eqs. (1)-(7),

r - k1 , w — w ■ k2 - k1 , U — h uc T -T1 - T0 T0

P0gp0(T1 - T)w2 Pr — CP01 Kn — II ln

l0 ? k0 '

2 — r.

ki « g — T ' Pm —

B = b(Ti - To),

_ 2 - rt 2y 1 \/pRTo/21o -, Pt —--" —, k —-

rt y + 1 Pr

Eqs. (1) and (2) can be written in dimensionless form as

[g + (1 - g)R] dR

[g + (1 - g)R] dR

[g + (1 - g)R] exp (-Bh) dR

[g + (1 - g)R]

h — 0 (9)

subject to the following dimensionless boundary conditions [8]

U(0)—P'KndU

' U(1)—-PKdU

h(0) — 1 + PvKn ln—

, h(1) —-PvKn ln —

(11) (12)

The physical quantities used in the above equations are defined in the nomenclature.

By using the transformation Z — g + (1 - g)R, the Eqs. (9)-(12) can be written as follows:

1 d ZdZ

1 d Z~dZ

Z exp(-Bh) f

(1 - g)2

subject to the boundary conditions

U(g)— PmKn(1 - g)dZ

U(1) — -PmKn(1 - g)dZ

h(g) — 1 + PvKn ln(1 - g)dZ

— -PmKn ln(1 - g)—

Integrating Eq. (14) and applying the boundary conditions (16) give the following:

h(Z)— A0 + A1 In(Z) (17)

In(g)-PmKnF(1 - g)( 1 + g

A0 — -PmKnF(1 - g)A1

Substituting Eq. (17) into the momentum Eq. (13) and solving it using the boundary condition (15) give

U(Z) — B0 +

B1C1(Z)

A0C1(Z)E

(1 - g)2 ,2 (E1 ln(Z) - 1) -

(Z)E1 4E,

Two important parameters for convective micro-flow are the volume flow rate, and skin-friction. The temperature solution as well as rate of heat transfer is exactly the same as discussed by Chen and Weng [8].

The dimensionless volume flow rate is as follows:

2%w2ur

(1 - g)2

ZU(Z)dZ

By substituting Eqs. (19) into (20) and integrating Eq. (20) give the following:

(1 - g)2 1

B1I1C1

(1 - g)4

A0C1I2 , r 1

J2_ 4E1

The skin-frictions (s) at the cylinder walls are as follows:

dU dR 1

— i(t^|r—0

S0 — e-Bh(1 - g)ddU |

dZ|Z—g

S0 — exp [-BfA0 -

- A1 ln(g)g] [B1C1(g)(C2-1)- F2] (1

l^dU. S1 — l(T) dR jR—1

S1 — ^ 9(1 - g)1

S1 — exp[-BA0][B1C1 - F6](1 - g)

where C1; C2, B0, B1; E1, E2,11 ,... ,/4 , F1 ,... ,F10 are all constants given in the Appendix A.

In order to verify the accuracy of the present work, we have computed the numerical value for the velocity for small value of B. Table gives a comparison of the numerical values of the velocity obtained in the present work when B ! 0 with those obtained by Weng and Chen [8] for ln = 1.667. As can be seen from Table 1, the solutions of the present work agree with those of Weng and Chen [8] for small values of B.

3. Results and discussion

MATLAB program is written to compute and generate line graphs for velocity, volume flow rate, and skin friction at both cylinders for different values of the dimensionless parameters, such as rarefaction parameter (pvKn), fluid-wall interaction parameter (ln), radius ratio (g) and viscosity variation parameter (B) so as to comment on their relative significance in the flow formation. The present parametric study has been performed in the continuum and slip flow regimes (Kn 6 0.1).

Table 1 Comparison of the values of velocity obtained in the present work with those obtained by Weng and Chen [8].

g R Velocity (ln = 1.64)

Weng and Chen [8] Present work (B ! 0)

0.8 0.2 0.0713 0.0712

0.4 0.0838 0.0837

0.6 0.0746 0.0744

0.8 0.0505 0.0503

0.5 0.2 0.0651 0.0650

0.4 0.0730 0.0725

0.6 0.0624 0.0621

0.8 0.0409 0.0405

0.2 0.2 0.0527 0.0526

0.4 0.0538 0.0536

0.6 0.0431 0.0430

0.8 0.0270 0.0269

Also, for air and various surfaces, the values of bv and bt range from near 1 to 1.667 and from near 1.64 to more than 10, respectively. So, this study has been performed over the reasonable ranges of 0 6 bvKn 6 0.1 and 0 6 ln 6 10. The selected reference values of pvKn, and ln for the present analysis are 0.05 and 1.64 respectively as given in Weng and Chen [8].

Temperature solution as well as rate of heat transfer is exactly the same as discussed by Weng and Chen [8]; therefore, we have only discussed the variation of velocity, volume flow rate and skin friction for different controlling parameters.

Fig. 2 shows the velocity distribution for different values of rarefaction parameter (bvKn) for fixed values of B — 0.5 and ln — 1.64. It is observed that, as rarefaction parameter (pvKn) increases, the velocity slip at the cylindrical surfaces increases which reduces the retarding effect of the boundaries. This yields an observable increase in the fluid velocity. Furthermore, the slip induced by rarefaction effect increases as radius ratio (g) increases.

Figure 2 Velocity profile for different values of Kn with ln = 1.64, B = 0.5.

Fig. 3 depicts the effect of fluid-wall interaction parameter (ln) on the velocity distribution for fixed values of pvKn — 0.05 and B — 0.5. It is evident from Fig. 3 that, the increase in fluid-wall interaction parameter (ln) leads to the decrease in fluid velocity and increase in slip velocity near the outer surface of the inner cylinder while reverse trend is observed at inner surface of the outer cylinder. In addition, the slip induced by fluid-wall interaction parameter (ln) increases as radius ratio (g) decreases while the impact of fluid-wall interaction parameter (ln) on the slip is more visible for smaller radius ratio (g).

Fig. 4 exhibits the effect of viscosity variation parameter (B) on velocity distribution for fixed values of pvKn — 0.05 and ln — 1.64 It is observed that increase in viscosity variation parameter (B) enhances the fluid velocity and slip velocity while the profiles for constant viscosity (B — 0) lie above B < 0 and below B > 0 on velocity field. This is due to the fact that as viscosity variation parameter (B > 0) increases, the viscosity of the working fluid decreases due to exponentially decaying nature of viscosity which causes higher velocity and velocity slip while as the viscosity variation parameter (B < 0) increases, the viscosity of the working fluid increases which causes lower velocity and velocity slip. It is also evident in Fig. 4 that velocity slip increases for B > 0 and decreases for B < 0 with increase in radius ratio (g).

Fig. 5 illustrates the effects of rarefaction parameter (bvKn) and fluid-wall interaction parameter (ln) on volume flow rate (Q). It is interesting to note that the volume flow rate (Q) is a decreasing function of fluid-wall interaction parameter (ln). Furthermore, it is found that increase in radius ratio (g) and rarefaction parameter (bvKn) leads to increase in the volume flow rate.

Fig. 6 reveals the influences of rarefaction parameter (pvKn) and viscosity variation parameter (B) on volume flow rate (Q) It is observed from Fig. 6 that, the volume flow rate (Q) increases, as viscosity variation parameter (B) increases. In addition, it is interesting to note that increase in rarefaction parameter (bvKn) leads to increase in the volume flow rate.

Fig. 7 exhibits the effects of rarefaction parameter (pvKn) and fluid-wall interaction parameter (ln) on skin-friction at

Figure 3 Velocity profile for different values of ln with Kn = 0.05, B = 0.5.

Figure 4 Velocity profile for different values of B with Kn = 0.05, ln = 1.64.

Figure 5 Variation of volume flow rate with Kn at different values of ln.

Figure 6 Variation of volume flow rate with Kn at different values of B.

Figure 7 Variation of skin friction on the inner cylinder with Kn at different values of ln.

outer surface of inner cylinder (R = 0). It is clear that the skin-friction decreases with the increase in fluid-wall interaction parameter (ln) and rarefaction parameter (bvKn). It is interesting to note that the impact of these parameters is significant for small value of radius ratio (g).

Fig. 8 shows the effects of rarefaction parameter (bvKn) and fluid-wall interaction parameter (ln) on skin-friction at inner surface of outer cylinder (R — 1). It is evident from Fig. 8 that increases in radius ratio (g) and rarefaction parameter (bvKn) lead to the increase in the skin-friction at inner surface of outer cylinder.

Figs. 9 and 10 illustrate the effects of rarefaction parameter (bvKn) and viscosity variation parameter (B) on skin-friction at outer surface of inner cylinder (R — 0) and inner surface of outer cylinder (R — 1), respectively. It is observed from the these Figures that skin friction increases, as viscosity parameter (B) increases.

Figure 8 Variation of skin friction on the outer cylinder with Kn at different values of ln.

<= 0.8,;-.-.

■ 1 ■

=-1.5,-1.0,-0.5,0.5,1.0,1.5

0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 Kn

Figure 9 Variation of skin friction on the inner cylinder with Kn at different values of B.

' 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 Kn

Figure 10 Variation of skin friction on the outer cylinder with Kn at different values of B.

4. Conclusions

Steady fully developed natural convection flow formation in a vertical annular micro-channel having temperature dependent viscosity in the presence of velocity slip and temperature jump at the annular micro-channel surfaces is considered. The role of radius ratio (g), viscosity variation parameter (B), rarefaction parameter (bvKn), and fluid-wall interaction parameter (ln) on the fluid velocity, volume flow rate and skin-friction is investigated. This study agrees with the findings of Weng and Chen [8] for vanishing viscosity variation parameter (B). The main findings are as follows:

I. Increasing the values of viscosity variation parameter (B) enhances the velocity and velocity slip.

II. The volume flow rate increases with the increase in the values of both viscosity variation parameter (B) and Knudsen number (Kn).

III. Skin friction increases as viscosity variation parameter (B) increases.

Appendix A

Constants used in the present work.

C1 = exp(B^), C2 = BAU E = C2 + 2, E2 = E + 1,

B _ F8F9-F7F10 B _ F8-F10

-L>0 F9-F7 > 1 F7-F9 ,

F1 —

-g)2 h ^ L

A0C1(g)E1 2E1 +

A0C1(g)(C2+1'

A^j^ (Eln(g) —1) —<4^}], 1 }].

-A,C,{ ^ln( g)

= C1(g)C C2

^ lAf — A1C>

{2E2 + '

F4 — C1(g)C2—1; [A

F - C1 [Aa - A1I

F6 — (1— g)2 L 2 4 J '

F7 — F - bvKn(1 — g)F4,F8 — F — ßvKn(1 — g)F2, F — (§) + ßvKn(1 — g)C1,

F10 — ßvKn(1 — g)Fs + F5,

[1—(g)E1 ] , [1—(g)E;'

L, I2 —J

_ _(#2 £2

^ ln(g)—£2 [1 — (g)E2 ],

14 — 2 [1 — (g)'],

References

[1] C.K. Chen, H.C. Weng, Natural convection in a vertical microchannel, J. Heat Transf. 127 (2005) 1053-1056.

[2] B.K. Jha, B. Aina, S.B. Joseph, Natural convection flow in a vertical micro-channel with suction/injection, J. Process Mech. Eng. 228 (3) (2014) 171-180.

[3] A.F. Khadrawi, A. Othman, M.A. Al-Nimr, Transient free convection fluid flow in a vertical micro-channel as described by the hyperbolic heat conduction model, Int. J. Thermophys. 26 (3) (2005) 905-918.

[4] O.M. Haddad, M.A. Al-Nimr, M.M. Abuzaid, The effect of frequency of fluctuating driving force on basic slip micro-flows, Acta Mech. 179 (2005) 249-259.

[5] B. Buonomo, O. Manca, Natural convection flow in a vertical micro-channel with heated at uniform heat flux, Int. J. Therm. Sci. 49 (2012) 1333-1344.

[6] B. Buonomo, O. Manca, Transient natural convection flow in a vertical micro-channel with heated at uniform heat flux, Int. J. Therm. Sci. 56 (2012) 35-47.

[7] C.K. Chen, H.C. Weng, Developing natural convection with thermal creep in a vertical microchannel, J. Phys. D 39 (2006) 3107-3118.

[8] H.C. Weng, C.K. Chen, Drag reduction and heat transfer enhancement over a heated wall of a vertical annular microchannel, Int. J. Heat Mass Transf. 52 (2009) 1075-1079.

[9] Basant Kumar Jha, Babatunde Aina, S.A. Muhammad, Combined effects of suction/injection and wall surface curvature on natural convection flow in a vertical annular micro-channel, J. Thermophys. Aeromech. 22 (2) (2015) 217228.

[10] Basant K. Jha, Babatunde Aina, Steady fully developed mixed convection flow in a vertical annular microchannel: an exact solution, ASIAN J. Math. Comput. Res. 10 (4) (2016) 332-347.

[11] M. Avci, O. Aydin, Mixed convection in a vertical parallel plate microchannel, ASME J. Heat Transf. 129 (2007) 162-166.

[12] M. Avci, O. Aydin, Mixed convection in a vertical parallel plate microchannel with asymmetric wall heat fluxes, ASME J. Heat Transf. 129 (2007) 1091-1095.

[13] M. Avci, O. Aydin, Mixed convection in a vertical microannulus between two concentric microtubes, ASME J. Heat Transf. 131 (2009) 014502-014504.

[14] B.K. Jha, Babatunde Aina, Mathematical modelling and exact solution of steady fully developed mixed convection flow in a vertical micro-porous-annulus, J. Afr. Matematika 26 (2015) 1199-1213.

[15] Basant K. Jha, Babatunde Aina, MHD mixed convection flow in a vertical micro-annulus: an exact solution, Int. J. Fluid Mech. Res. 42 (6) (2015) 537-552.

[16] S. Das, S. Chakraborty, R.N. Jana, O.D. Makinde, Mixed convective Couette flow of reactive nanofluids between concentric vertical cylindrical pipes, J. Nanofluids 4 (4) (2015) 485-493.

[17] S. Das, R.N. Jana, O.D. Makinde, Mixed convective magnetohydrodynamic flow in a vertical channel filled with nanofluids, Eng. Sci. Technol. Int. J. 18 (2) (2015) 244-255.

[18] O.D. Makinde, Thermal analysis of a reactive generalized Couette flow of power law fluids between concentric cylindrical pipes, European Phys. J. Plus 129 (2014) 1-9.

[19] H. Herwig, G. Wicken, The effect of variable properties on laminar boundary layer flow, Warme-und Stof-Fubertrag 20 (1986) 47-57.

[20] O.D. Makinde, T. Chinyoka, Analysis of unsteady flow of a variable viscosity reactive fluid in a slit with wall suction or injection, J. Petrol. Sci. Eng. 94 (95) (2012) 1-11.

[21] M.S. Tshehla, O.D. Makinde, G.E. Okecha, Heat transfer and entropy generation in a pipe flow with temperature dependent viscosity and convective cooling, Sci. Res. Essays 5 (23) (2010) 3730-3741.

[22] K. Klemp, H. Herwig, M. Selmann, Entrance flow in channel with temperature dependent viscosity including viscous dissipation effects, Int. Proceedings of the Third International Congress of Fluid Mechanics, Cairo, Egypt 3 (1990) 1257-1266.

[23] M.A. Seddeek, Effects of radiation and variable viscosity on a MHD free convection flow past a semi-infinite flat plate with an

aligned magnetic field in the case of unsteady flow, Int. J. Heat Mass Transf. 45 (2002) 931-935.

[24] O. Grasset, E.M. Parmentier, Thermal convection in a volumetrically heated, infinite Prandtl number fluid with strongly temperature-dependent viscosity: implications for planetary thermal evolution, J. Geophys. Res. 103 (1998) 18171-18181.

[25] S. Saravanan, P. Kandaswamy, Hydromagnetic stability of convective flow of variable viscosity fluids generated by internal heat sources, Z. Angew. Math. Phys. 55 (2004) 451-467.

[26] F.M. Hady, A.Y. Bakier, R.S.R. Gorla, Mixed convection boundary layer flow on a continuous flat plate with variable viscosity, Heat Mass Transf. 31 (1996) 169-172.

[27] M.A.A. Mahmud, A note on variable viscosity and chemical reaction effects on mixed convection heat and mass transfer along a semi-infinite vertical plate, Math. Probab. Eng. (2007), doi:1155/2007/41323 Article ID 413237.

[28] M. Kumari, Variable viscosity effects on free and mixed convection boundary-layer flow from a horizontal surface in a saturated porous medium - variable heat flux, Mech. Res. Commun. 28 (2001) 339-348.

[29] J.C. Umavathi, Odelu Ojjela, Effect of variable viscosity on free convection in a vertical rectangular duct, Int. J. Heat Mass Transf. 84 (2015) 1-15.

[30] J.C. Umavathi, Combined effect of variable viscosity and variable thermal conductivity on double-diffusive convection flow of a permeable fluid in a vertical channel, Transp. Porous Med. 108 (2015) 659-678.

[31] H. Lamb, Hydrodynamics, sixth ed., Cambridge Univ. Press, London, 1932.

[32] M. Jacob, Heat Transfer, John Wiley and Sons Inc, New York, 1949.

[33] H.T. Syeda, M. Shohel, Entropy generation in a vertical concentric channel with temperature dependent viscosity, Int. Commun. Heat Mass Transf. 29 (2002) 907-918.