Scholarly article on topic 'Fully developed MHD natural convection flow in a vertical annular microchannel: An exact solution'

Fully developed MHD natural convection flow in a vertical annular microchannel: An exact solution Academic research paper on "Chemical engineering"

CC BY-NC-ND
0
0
Share paper
OECD Field of science
Keywords
{"Annular micro-channel" / "Natural convection" / "Hartmann number" / "Velocity slip and temperature jump"}

Abstract of research paper on Chemical engineering, author of scientific article — Basant K. Jha, Babatunde Aina, Sani Isa

Abstract An exact solution of steady fully developed natural convection flow of viscous, incompressible, electrically conducting fluid in a vertical annular micro-channel with the effect of transverse magnetic field in the presence of velocity slip and temperature jump at the annular micro-channel surfaces is obtained. Exact solution is expressed in terms of modified Bessel function of the first and second kind. The solution obtained is graphically represented and the effects of radius ratio (η), Hartmann number (M), rarefaction parameter (βvKn), and fluid–wall interaction parameter (F) on the flow are investigated. During the course of numerical computations, it is found that an increase in Hartmann number leads to a decrease in the fluid velocity, volume flow rate and skin friction. Furthermore, it is found that an increase in curvature radius ratio leads to an increase in the volume flow rate.

Academic research paper on topic "Fully developed MHD natural convection flow in a vertical annular microchannel: An exact solution"

Journal of King Saud University - Science (2014) xxx, xxx-xxx

King Saud University Journal of King Saud University Science

www.ksu.edu.sa www.sciencedirect.com

Journal of

King Saud University

-Science

ORIGINAL ARTICLE

Fully developed MHD natural convection flow in a vertical annular microchannel: An exact solution

Basant K. Jha, Babatunde Aina *, Sani Isa

Department of Mathematics, Ahmadu Bello University, Zaria, Nigeria Received 30 September 2014; accepted 2 December 2014

KEYWORDS

Annular micro-channel; Natural convection; Hartmann number; Velocity slip and temperature jump

Abstract An exact solution of steady fully developed natural convection flow of viscous, incompressible, electrically conducting fluid in a vertical annular micro-channel with the effect of transverse magnetic field in the presence of velocity slip and temperature jump at the annular micro-channel surfaces is obtained. Exact solution is expressed in terms of modified Bessel function of the first and second kind. The solution obtained is graphically represented and the effects of radius ratio (g), Hartmann number (M), rarefaction parameter (bvKn), and fluid-wall interaction parameter (F) on the flow are investigated. During the course of numerical computations, it is found that an increase in Hartmann number leads to a decrease in the fluid velocity, volume flow rate and skin friction. Furthermore, it is found that an increase in curvature radius ratio leads to an increase in the volume flow rate.

© 2014 The Authors. Production and hosting by Elsevier B.V. on behalf of King Saud University. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).

1. Introduction

Recently, a growing interest in micro-channel fluid mechanics and heat transfer has emerged because of possible cooling applications in space systems, manufacturing and material processing operations, and in high-power-density chips in supercomputers and other electronics (Al-Nimr and

* Corresponding author. Tel.: +234 8063178775. E-mail addresses: basant777@yahoo.co.uk (B.K. Jha), ainavicdydx @gmail.com (B. Aina), sanimath@yahoo.com (S. Isa). Peer review under responsibility of King Saud University.

Khadrawi, 2004) . As this area continues to grow, it becomes increasingly important to understand the mechanisms and fundamental differences involved with fluid mechanics and heat transfer mechanisms in macro-channel and micro-channel.

A series of investigations have been conducted recently in the field of micro geometry flow. However, to cite a few works in this direction, Chen and Weng (2005) analytically studied the fully developed natural convection in open-ended vertical parallel plate micro-channel with asymmetric wall temperature distribution in which the effect of rarefaction and fluid wall interaction was shown to increase the volume flow rate and decrease the heat transfer. This result is further extended by taking into account suction/injection on the micro-channel walls by Jha et al. (2014). They concluded that skin frictions as well as the rate of heat transfer are strongly dependent on suction/injection parameter. In another work, Weng and

http://dx.doi.Org/10.1016/j.jksus.2014.12.002

1018-3647 © 2014 The Authors. Production and hosting by Elsevier B.V. on behalf of King Saud University. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/3.0/).

Nomenclature

B0 constant magnetic flux density w dimensional gap between the cylinders

Cp0 specific heat at constant pressure thermal and tangential momentum accommoda

F fluid-wall interaction parameter, fit/bv tion coefficients, respectively

g gravitational acceleration

K radius of the inner cylinder Greek letters

k2 radius of the outer cylinder a thermal diffusivity

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

M Hartmann number bt, bv dimensionless variables

q volume flow rate c ratio of specific heats

Q dimensionless volume flow rate l0 dynamic viscosity

Pr Prandtl number h dimensionless temperature

r dimensional radial coordinate P0 density

R dimensionless radial coordinate V fluid kinematic viscosity (i0/p0)

R specific gas constant g ratio of radii (kj/k2)

T temperature of fluid k molecular mean free path

T0 reference temperature k0 thermal conductivity

T1 temperature at outer surface of the inner cylinder r electrical conductivity of the fluid

u axial velocity s skin-friction

U dimensionless axial velocity

Chen (2009) studied the impact of wall surface curvature on steady fully developed natural convection flow in an open-ended vertical micro-annulus with an asymmetric heating of the annulus surface. Recently, Jha et al. (in press) further extended the work of Weng and Chen (2009) by taking into account suction/injection on a vertical annular micro-channel. They discovered that skin-friction decreases at the outer surface of the inner porous cylinder with an increase of fluid-wall interaction parameter in the case of injection at the outer surface of the inner porous cylinder and simultaneous suction at inner surface of the outer porous cylinder while the result is just reverse at the inner surface of outer porous cylinder. In a related article, Avci and Aydin (2009) studied the fully developed mixed convective heat transfer of a Newtonian fluid in a vertical micro-annulus between two concentric micro-tubes. Recently, Jha and Aina (2014) further extended the work of Avci and Aydin (2009) to the case when the vertical micro-annulus formed by two concentric micro-tubes is porous, i.e. where there is suction or injection through the annulus surfaces.

On the other hand, the MHD phenomenon has received considerable attention during the last two decades due to its importance from the energy generation point of view, and one may envisage MHD generators for power generation. MHD pumps are already in use in chemical energy technology for pumping electrically conducting fluids at some of the atomic energy centres. Besides these applications, when the fluid is electrically conducting, the free convection flow is appreciably influenced by an imposed magnetic field. Therefore, to refer to few works in this direction, Sheikholeslamia et al. (2014a) investigated the magnetic field effect on nanofluid flow and heat transfer in a semi-annulus enclosure via control volume based finite element method. Khan and Ellahi (2008) observed the effects of magnetic field and porous medium on some unidirectional flows of a second grade fluid. Farhad et al. (2012) examines the slip effect on hydromagnetic rotating flow of viscous fluid through a

porous space. In another work, Ali et al. (2012) investigated the effects of slip condition on the unsteady magnetohydro-dynamics (MHD) flow of incompressible viscoelastic fluids in a porous channel under the influence of transverse magnetic field and Hall current with heat and mass transfer. An analysis to investigate the combined effects of heat and mass transfer on free convection unsteady magnetohydrody-namics (MHD) flow of the viscous fluid embedded in a porous medium is presented by Ali et al. (2013a)

Some recent works related to the present investigation are found in the literature (Ali et al., 2013b, 2014; Sheikholeslamia et al., 2014b; Sheikholeslami et al., 2012a,b; Ashorynejad et al., 2013). In Ali et al. (2013b), Farhad et al. presented an exact analysis of combined effects of radiation and chemical reaction on the magnetohydrodynamics (MHD) free convection flow of an electrically conducting incompressible viscous fluid over an inclined plate embedded in a porous medium. Ali et al. (2014) studied the unsteady free convection flow of a second grade fluid past an isothermal vertical plate oscillating in its plane with constant viscosity. Also, Sheikholeslamia et al. (2014b) studied the magnetohydrody-namic effect on natural convection heat transfer of Cu-water nanofluid in an enclosure with a hot elliptic cylinder. Sheikholeslami et al. (2012a) numerically examined the natural convection of nanofluids in a concentric annulus between a cold outer square cylinder and a heated inner circular. Flow and heat transfer of a nanofluid over a stretching cylinder in the presence of magnetic field has been investigated by Ashorynejad et al. (2013). Recently, Sheikholeslami et al. (2012b) presented the numerical solution for natural convection of nanofluids in a cold outer circular enclosure containing a hot inner sinusoidal cylinder. Sheikholeslamia et al. (2013) carried out a numerical investigation on natural convection nanofluid flow in a half annulus enclosure with one wall under constant heat flux in the presence of a magnetic field.

However, to the best knowledge of the authors, no studies have been carried out on the fully developed MHD natural

convection flow in a vertical annular micro-channel. This prompts the present work.

, I4 2 — r .du

u(r — kl) —-к-

r—ki

2. Mathematical analysis

The flow considered is a fully developed steady natural convection flow of viscous, incompressible, electrically conducting fluid in vertical annular micro-channel under the effect of transverse magnetic field as presented in Fig. 1. The 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 are k1 and k2 respectively. A magnetic field of uniform strength B0 is assumed to be applied in the direction perpendicular to the direction of flow. It is assumed that the magnetic Reynolds number is very small, which corresponds to negligibly induced magnetic field compared to the externally applied one (Sarris et al., 2006; Jha and Apere, 2013; Crammer and Pai, 1973; Pai, 1962). The outer surface of the inner cylinder is heated to a temperature (T1) 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). In the recent past, Weng and Chen (2009) analysed the drag reduction and heat transfer enhancement of viscous incompressible fluid in vertical micro-annulus. Following Weng and Chen and considering the viscous incompressible electrically conducting fluid in the presence of a transverse magnetic field to the flow direction, the mathematical model for the present physical situation under Boussinesq's approximation in dimensional form is

v d / du\ ^ ь (t r dry dr I 0

— 0;

l d ( dT

r dr dr

— 0.

The boundary conditions for the velocity and temperature field in the presence of velocity slip and temperature jump are (Weng and Chen, 2009; Avci and Aydin, 2009):

Figure 1 Flow configuration and coordinate system.

, , -, 2 — ffv .du

u(r — k2) —--к-

T(r — ki) — Ti +

2 — rt 2y k dT

T(r — k2) — To

rt y + l Pr dr 2 — rt 2y к dT

rt y + l Pr dr

By introducing the following non-dimensional quantities

r — ki u T — T0

R = -, w = k2 — ki, U = —, h = —-—,

w uc Ti — T0

. Pogßo(Ti — To)w lo

CPo 1 ko ' 2 — rt 2y l

к ßt Kn — - , F — ßt ;

kl a 2 — rv „

g —T; ßv —-; ßt — - 1 и '

k2 Гт rt y + l Pr

nRTo/2io

rB0 w2 Pov ;

Using Eq. (7) in Eqs. (1)-(6), the mathematical model representing the present physical situation in dimensionless form in the presence of velocity slip and temperature jump on annular microchannel surfaces is:

[g + (l — g)R] dR

[g + (l — g)R] dR

[g+(l—g)R] dU

[g +(l — g)R] dR

— M2U + h — 0;

The boundary conditions which describe velocity slip and temperature jump conditions at the fluid-wall interface in dimen-sionless form are (Weng and Chen, 2009)

U(0)— ^

0(0) — l + ßvKnF-

U(l) — —

0(l) — —ßvKnF-

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

By using the transformation Z = g + (1 — g)R, the Eqs. (8)-(11) can be written as:

= 0; (12)

l d ZdZ

l d ZdZ

(l — g)2 (l — g)2 0.

Subject to the boundary conditions

U(g)—ßvKn(l — g)dZ

U(l) ——ßvKn(l — g)dZ

0(g) — l + ßvKnF(l — g)—

0(l) — —ßvKnF(l — g)—

Integrating equation (13) and applying the boundary conditions (15) gives:

d(Z)=A0 + A1 In(Z) (16)

where:

A1 =-1--v, A0 = -&KnF(1 - g)A1.

In(g)-bvKnF(1 - g) 1

The general solution of Eq. (12) after substitution equation (16) in Eq. (12) is obtained by rearranging the homogeneous part to take the general form of the Bessel equation (Arpaci, 1966). The particular solution of Eq. (12) is sought by assuming an appropriate form (according to the forcing function i.e., the right hand side of Eq. (12))).The exact solution of Eq. (12) under the appropriate velocity slip condition defined in Eq. (13) is:

U(Z) = QIq(ZE) + C2K0(ZE)+M [Aq + AlIn(Z)}. (18)

Two important parameters for convective micro-flow are the volume flow rate, and skin-friction. The dimensionless volume flow rate is:

ZU(Z)dZ.

2nw2uc (1 - g)2 j g

Substituting equation (18) into Eq. (19), one obtains

C1I1(E)-C2K1(E) g(C2K1(Eg)-C1I1(Eg))

(1 - g)2

+(1 - g)2

2M2 4M2 \ 2M2 "

g2In(g)

The skin-frictions (s) at the cylinder walls are obtained by differentiating the velocity as follows:

.. JU = (1 - g)dZ

{C1I1(gE)-C2K1(gE)}E + -A-M2g

(1 - g);

n \dU = (1 - g)dZ

{C1I1(E)-C2K!(E)E}+ -2

(1 - g);

(21) (22)

where Ci, C2, E, Ei,..., E6 are all constants given in the Appendix A.

3. Results and discussion

Numerical values of the velocity, volume flow rate, and skin friction are computed and the graphs are presented in Figs. 210 so that the influences of radius ratio (g), Hartmann number (M), rarefaction parameter (pvKn), and fluid-wall interaction parameter (F) can be easily seen at a glance. The present parametric study has been performed in the continuum and slip flow regimes (Kn 6 0.1). 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

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Figure 2 Velocity distribution for different values of bvKn with F = 1.64 and M = 2.0.

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R

Figure 3 Velocity distribution for different values of F with bvKn = 0.05 and M = 2.0.

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 R

Figure 4 Velocity distribution for different values of M with bvKn = 0.05 and F = 1.64.

Figure 5

with M =

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 PvKn

Volume flow rate versus bvKn for different values of F 2.0.

0.13 0.12 0.11 ) 0.1 § 0.09 m 0.08 0.07 0.06 0.05

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

Figure 8 Variation of skin friction for different values of F at (R = 1).

Figure 6 Volume flow rate versus bvKn for different values of M with F = 1.64.

Figure 9 Variation of skin friction for different values of M at (R = 0).

Figure 10 Variation of skin friction for different values of M at Figure 7 Variation of skin friction for different values of F at (R =1).

(R = 0).

performed over the reasonable ranges of 0 6 ftvKn 6 0.1 and

0 6 F 6 10. The selected reference values of pvKn, and F for the present analysis are 0.05 and 1.64 respectively as given in Weng and Chen (2009).

The present work extends the work of Weng and Chen (2009) by considering the influence of externally applied transverse magnetic field on a steady fully developed natural convection flow of viscous, incompressible, electrically conducting fluid in a vertical annular micro-channel. The values of Hartmann number (M) are taken over the range of

1 6 M 6 2 with a reference value of M =2.

The expression for the temperature in Eq. (10), the effect of the rarefaction parameter (pvKn), and fluid-wall interaction parameter (F) on the temperature profile and rate of heat transfer which is expressed as the Nusselt number are exactly the same as those given by Weng and Chen (2009).

Fig. 2 shows the velocity distribution for different values of rarefaction parameter (pvKn) with fixed values of M = 2.0 and F = 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 boundaries. This yields an observable increase in the fluid velocity near both boundaries. This is due to the fact that, as (pvKn) increases, the temperature jump increases and this reduces the amount of heat transfer from the cylindrical surfaces to the fluid (Weng and Chen, 2009). This reduction in heat transfer reduces the buoyancy effect, which derives the flow and hence reduces the fluid velocity far from the outer surface of inner cylinder (near the inner surface of outer cylinder). The reduction in velocity due to the reduction in heat transfer is offset by the increase in velocity due to the reduction in the frictional retarding forces near the cylindrical surfaces. Furthermore, the slip induced by rarefaction effect increases as radius ratio (g) increases.

Fig. 3 shows the velocity distribution for different values of fluid-wall interaction parameter (F) with fixed values of ftv-Kn = 0.05 and M = 2.0. It is evident from Fig. 3 that, the increase in fluid-wall interaction parameter (F) leads to the decrease in fluid velocity and increase in slip velocity near the outer surface of the inner cylinder while the reverse trend is observed at the inner surface of the outer cylinder. In addition, the slip induced by fluid-wall interaction parameter (F) increases as the radius ratio (g) decreases while, the impact of fluid-wall interaction parameter (F) on the slip is more visible for smaller radius ratio (g).

Fig. 4 depicts the velocity distribution for different values of the Hartmann number (M) with fixed values of pvKn = 0.05 and F = 1.64. It is observed that as the Hartmann number (M) increases, there is a decrease in the fluid velocity and an increase in slip velocity. It is also observed that there is higher slip at the inner surface of outer cylinder compared to the outer surface of the inner cylinder. In addition, for the fixed Hartmann number (M), as radius ratio (g) increases there is an increase in the slip. Further, it is seen that the particular case of Eq. (18) corresponding M fi 0 reduces to the one in Weng and Chen (2009).

Fig. 5 shows the variation of volume flow rate (Q) against the rarefaction parameter (pvKn) for different values of fluid-wall interaction (F). It is interesting to note that the volume flow rate (Q) is a decreasing function of fluid-wall interaction (F). Furthermore, it is found that an increase in radius ratio (g)

and rarefaction parameter (fivKn) leads to an increase in the volume flow rate.

Fig. 6 reveals the volume flow rate (Q) plotted against the rarefaction parameter (pvKn) for different values of Hartmann number (M). It is observed from Fig. 6 that, the volume flow rate (Q) decreases as Hartmann number (M) increases. Also, it is interesting to note that an increase in the rarefaction parameter (pvKn) leads to an increase in the volume flow rate.

Fig. 7 depicts the variation of skin-friction at the outer surface of the inner cylinder (R = 0). It is clear from Fig. 7 that, the skin - friction decreases with the increase of fluid-wall interaction parameter (F) and rarefaction parameter (bvKn). It is interesting to note that the impact of these parameters is significant for a small value of radius ratio (g).

Fig. 8 shows a variation of skin-friction at the inner surface of the outer cylinder (R = 1). It is obvious from Fig. 8 that increases in radius ratio (g) and rarefaction parameter (pvKn) lead to the increase in the skin-friction at the inner surface of the outer cylinder.

Figs. 9 and 10 presents a variation of skin-friction at the outer surface of the inner cylinder (R = 0) and the inner surface of outer cylinder (R = 1), respectively for different values of the Hartmann number (M). It is observed from these figures that, skin friction decreases as the Hartmann number (M) increases.

4. Conclusions

This study considered the fully developed steady natural convection flow of a viscous, incompressible, electrically conducting fluid in vertical annular micro-channel under the effect of a transverse magnetic field. The influence of the radius ratio (g), Hartmann number (M), rarefaction parameter (pvKn), and fluid-wall interaction (F) on the fluid velocity, volume flow rate and skin-friction is analysed. This study exactly agrees with the finding of Weng and Chen (2009) in the absence of a magnetic field and further concludes as follows:

I. It is found that as the Hartmann number (M) increases, there is a decrease in the fluid velocity and an increase in slip velocity.

II. The slip induced by rarefaction effect and the Hartmann number increases as the radius ratio (g) increases while the slip induced by fluid-wall interaction parameter (F) increases as the radius ratio (g) decreases.

III. The increase in radius ratio (g) and rarefaction parameter (bvKn) leads to an increase in the volume flow rate.

IV. The increase in the Hartmann number (M) leads to a decrease in the volume flow rate.

V. Finally, the skin friction decreases as the Hartmann number (M) increases.

Appendix A

Constants used in the present work.

[E6E2 - E3E5] [E1E6 - E3E4] M

1 [E2E4 - E1E5]; 2 [E1E5 - E2E4]; (1 - g);

E1 = [/0(Eg)-(b„Kn(1 - g)E)h(gE)], E2 = [K>(Eg) + (ß„Kn(1 - g)E)^(gE)],

E3 = -¡4"ß„Kn(1 - g) - -L [A0 + A1/n(g)] M g M

E4 = [/0(E) + ß„Kn(1 - g)EI1(E)], E5 = [K0(E)- ß„Kn(1 - g)EK1(E)],

E6 — —

A0 , ß„Kn(1 - g)A1

References

Ali, Farhad, Khan, Ilyas, Samiulhaq, Mustapha, Norzieha, Shafie, Sharidan, 2012. Unsteady magnetohydrodynamic oscillatory flow of viscoelastic fluids in a porous channel with heat and mass transfer. J. Phys. Soc. Jpn. 81, 064402.

Ali, Farhad, Khan, Ilyas, Shafie, Sharidan, Musthapa, Norzieha, 2013a. Heat and mass transfer with free convection MHD flow past a vertical plate embedded in a porous medium (Special Issue Mathematical Modeling of Heat and Mass Transfer in Energy Science and Engineering (MMTP). Math. Prob. Eng. 2013, http:// dx.doi.org/10.1155/2013/346281.

Ali, Farhad, Khan, Ilyas, Samiulhaq, Shafie, Sharidan, 2013b. Conjugate effects of heat and mass transfer on MHD free convection flow over an inclined plate embedded in a porous medium. PLoS ONE 8 (6), e65223.

Ali, Farhad, Khan, Ilyas, Shafie, Sharidan, 2014. Closed form solutions for unsteady free convection flow of a second grade fluid over an oscillating vertical plate. PLoS ONE 9 (2), e85099.

Al-Nimr, M.A., Khadrawi, A.F., 2004. Thermal behavior of a stagnant gas confined in a horizontal microchannel as described by the dual-phase-lag heat conduction model. Int. J. Thermophys. 25 (6), 1953-1964.

Arpaci, V.S., 1966. Conduction Heat Transfer. AddisonWesley, Reading, MA, pp. 135-136.

Ashorynejad, H.R., Sheikholeslami, Pop, I., Ganji, D.D., 2013. Nanofluid flow and heat transfer due to a stretching cylinder in the presence of magnetic field. Heat Mass Transfer 49, 427-436.

Avci, M., Aydin, O., 2009. Mixed convection in a vertical microan-nulus between two concentric Microtubes. ASME J. Heat Transfer 131, 014502-014504.

Chen, C.K., Weng, H.C., 2005. Natural convection in a vertical microchannel. J. Heat Transfer 127, 1053-1056.

Crammer, K.R., Pai, S., 1973. Magnetohydrodynamics for Engineers and Applied Physicist. McGraw Hill Book Company, New York.

Farhad, A., Norzieha, M., SSharidan, S., Khan, I., Samiuthaq, 2012. On hydromagnetic rotating flow in a porous medium with slip condition and hall current. Int. J. Phys. Sci. 7 (10), 1540-1548.

Jha, B.K., Aina, Babatunde, 2014. Mathematical modelling and exact solution of steady fully developed mixed convection flow in a vertical micro-porous-annulus. J. Afrika Matematika. http:// dx.doi.org/10.1007/s13370-014-0277-4.

Jha, B.K., Apere, C.A., 2013. Time-dependent MHD Couette flow in a porous annulus. Commun. Nonlinear Sci. Numer. Simul. 18, 19591969.

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

Jha, B.K., Aina, Babatunde, Shehu, A.M., in press. Combined effects of suction/injection and wall surface curvature on natural convection flow in a vertical micro-porous-annulus. J. Thermophys. Aerodyn., in press.

Khan, I., Ellahi, R., 2008. Some unsteady MHD flows of a second grade fluid through porous medium. J. Porous Media 11 (4), 389400 (ISI Indexed, I.F 0.684).

Pai, S., 1962. Magnetohydrodynamics and Plasma Dynamics. Springer, Berlin.

Sarris, I.E., Zikos, G.K., Grecos, A.P., Vlachos, N.S., 2006. On the limits of validity of the low magnetic Reynolds number approximation in MHD natural convection heat transfer. Numer. Heat Transfer: Part B - Fundam. 50, 157-180.

Sheikholeslami, M., Gorji-Bandpay, M., Ganji, D.D., 2012a. Magnetic field effects on natural convection around a horizontal circular cylinder inside a square enclosure filled with nanofluid. Int. Commun. Heat Mass Transfer 39, 978-986.

Sheikholeslami, M., Gorji-Bandpy, M., Ganji, D.D., Soleimani, Soheil, Seyyedi, S.M., 2012b. Natural convection of nanofluids in an enclosure between a circular and a sinusoidal cylinder in the presence of magnetic field. Int. Commun. Heat Mass Transfer 39, 1435-1443.

Sheikholeslamia, M., Gorji-Bandpy, M., Ganji, D.D., Soleiman, Soheil, 2013. Effect of a magnetic field on natural convection in an inclined half-annulus enclosure filled with Cu-water nanofluid using CVFEM. Adv. Powder Technol. 24, 980-991.

Sheikholeslamia, M., Gorji-Bandpy, M., Ganji, D.D., Rana, P., Soleimani, Soheil, 2014a. Magnetohydrodynamic free convection of Al2O3-water nanofluid considering thermophoresis and Brown-ian motion effects. Comput. Fluid 94, 147-160.

Sheikholeslamia, M., Gorji Bandpy, M., Ellahi, R., Hassan, Mohsan, Soleimani, Soheil, 2014b. Effects of MHD on Cu-water nanofluid flow and heat transfer by means of CVFEM. J. Magn. Magn. Mater. 349, 188-200.

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