MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2018/9840862 9840862 Research Article MHD Mixed Convection Micropolar Fluid Flow through a Rectangular Duct http://orcid.org/0000-0001-9460-0469 Ayano Mekonnen Shiferaw 1 http://orcid.org/0000-0002-1518-6990 Sikwila Stephen T. 2 http://orcid.org/0000-0002-5829-8634 Shateyi Stanford 3 de Monte Filippo 1 University of Swaziland Private Bag 4 Kwaluseni Swaziland uniswa.sz 2 Sol Plaatje University Private Bag X5008 Kimberley South Africa spu.ac.za 3 University of Venda Private Bag X5050 Thohoyandou South Africa univen.ac.za 2018 912018 2018 05 07 2017 01 12 2017 05 12 2017 912018 2018 Copyright © 2018 Mekonnen Shiferaw Ayano et al. 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.

Mixed convection flow through a rectangular duct with at least one of the sides of the walls of the rectangle being isothermal under the influence of transversely applied magnetic field has been analyzed numerically in this study. The governing differential equations of the problem have been transformed into a system of nondimensional differential equations and then solved numerically. The dimensionless velocity, microrotation components, and temperature profiles are displayed graphically showing the effects of various values of the parameters present in the problem. The results showed that the flow field is notably influenced by the considered parameters. It is found that increasing the aspect ratio increases flow reversal, commencement of the flow reversal is observed after some critical value, and the applied magnetic field increases the flow reversal in addition to flow retardation. The microrotation components flow in opposite direction; also it is found that one component of the microrotation will show no rotational effect around the center of the duct.

1. Introduction

The theory of fluids with microstructures has been the subject of a large number of investigations. These are realistic and important from a technological point of view. The classical theories of continuum mechanics are inadequate to explain the microscopic manifestations of complex hydrodynamic behaviour. Microcontinuum theory or generalized continuum theories incorporate independent deformations of the microstructure inside of a material point. There are a number of microcontinuum theories by Eringen [1, 2], namely, couple stress, micropolar, microstretch, and micromorphic. These theories impose more or less constraints on the motion of microstructure inside of a material point. In microstretch theory, it is assumed that the microstructure of each material point can undergo expansion or contraction independently in addition to translation and rigid rotation. This theory is a generalization of micropolar theory, in which the microstructure can only have translation and rigid rotation. Micromorphic theory constitutes extensions of the classical field theories concerned with the deformations, motions, and electromagnetic interactions of material media, as continua, in microscopic time and space scales. In micromorphic theory, a material body is considered as a continuous collection of deformable particles, each with finite size and inner structure. Axtell et al.  discussed blood circulation in tissues as a multistate, multiphase porous media problem by simulating the fluid phase as a micromorphic continuum. While blood is known to be Newtonian, it also has electrically conducting properties. Using the finite element method, Ansari et al.  studied micromorphic first-order shear deformable plate element and the microstructure effect on the analysis of small-scale structures. Applications of micromorphic theory in micro-/nanoscale are discussed by formulating generalized micromorphic solid and fluid; see Lee and Wang . Within the framework of micromorphic and elasticity theory, microstructure effect is developed to describe the bending behaviour of microplates by Ansari et al. . Micropolar fluids are a subset of the micromorphic fluid theory introduced in a pioneering paper by Eringen . Micropolar fluids constitute an important branch of non-Newtonian fluid dynamics where microrotation effects as well as microinertia are exhibited. Traditional Newtonian fluids cannot precisely describe the characteristic of fluid with suspended particles. Micropolar fluids formulated by Eringen  include certain microscopic effects arising from the local structure and micromotions of the fluid elements and provide a mathematical model for the non-Newtonian fluid flow behaviour such as exotic lubricants, polymers, animal blood, and real fluids with suspensions. Eringen  extended the micropolar fluid theory and developed the theory of thermomicropolar fluids. An excellent review of micropolar fluids and their applications was given by Ariman et al. . Extensive review of foundations of the micropolar continuum mechanics is found in the book by Eremeyev et al. .

There have been numerous studies for the laminar flow and heat transfer through rectangular ducts . Cheng and Ou  have also studied the mixed convection heat transfer in the thermal entrance region of horizontal rectangular channels. Aung and Worku  analyzed mixed convection flow in ducts with asymmetric wall heat fluxes. Mahaney et al.  studied development of laminar mixed convection flow in a horizontal rectangular duct with uniform bottom heating and Abou-Ellail and Morcos  also carried out investigations for similar flow geometry buoyancy effects in the entrance region. Huang et al.  examined laminar mixed convection heat and mass transfer in vertical rectangular ducts with film evaporation and condensation. Yan and Lin [25, 26] numerically investigated the laminar or turbulent mixed convection flow in a vertical channel under simultaneous influence of the combined buoyancy effects of the thermal and mass diffusion for an air-water vapour system.

To the knowledge of the authors and in all of the above-mentioned studies, the combined effect of an applied magnetic field under strong buoyancy force on mixed convection on the micropolar fluid model which represents realistic non-Newtonian flow in a rectangular duct has not been well evaluated. This motivates the present investigation. In this paper, we study the problem of steady laminar, incompressible mixed convection micropolar fluid flow through a rectangular duct subject to applied magnetic field with at least one side of the rectangle being isothermal. Mathematical formulation has been presented for the governing equations and the solution is found numerically using the finite difference method. Results are discussed through velocity, microrotation, and temperature.

2. Mathematical Formulation

Consider a steady, laminar, incompressible, mixed convective heat transfer micropolar fluid through an infinite vertical rectangular duct given in Figure 1, with sides of length a and b in the presence of an external applied magnetic field. In the formulation of the problem, the following assumptions are made: Boussinesq approximations are valid, the flow is parallel and hence the only nonvanishing component is the Z (vertical) direction component, a uniform magnetic field is applied normal to the surfaces, and the induced magnetic field is assumed to be very small. The last assumption is considered to be valid for the case of small magnetic Reynolds number as done in previous works by Hayat et al. . At least one of the duct walls is kept isothermal and the temperature is uniform along the duct vertical direction; hence T/z=0.

The base of the rectangular duct and the coordinate axes.

The equations of motion of a micropolar fluid by Eringen  are as follows:(1)·ρV=0,λ+2μ+κ·V-μ+κ××V+κ×ϕ-P+ρf=ρV·V,α+β+γ·ϕ-γ××ϕ+κ×V-2κϕ+ρI=ρjV·ϕ,α2T=ρV·T,where j is microinertia, f and I are body force and body couple per unit mass, α, β, λ, μ, and κ are material constants (viscosity coefficients), P is the thermodynamic pressure, V=(0,0,w) is velocity, ϕ=(Γ1,Γ2,0) is microrotation, α is the thermal diffusivity, and T is the temperature of the fluid.

The constitutive equations, giving the stress tensor τkl and the couple stress tensor Mkl, are given by(2a)τkl=-P+λvr;rgkl+μvk;l+vl;k+κvl;k-ϵklrvr,(2b)Mkl=αvr;rgk;l+γvl;k,where gkl and ϵklr are the metric tensor and the covariant ϵ-symbol, respectively. The semicolon denotes the covariant partial differentiation with respect to a space coordinate and repeated indices are summed. The stress vector τ and couple stress vector M at a point on a surface with normal n are given by ταknα=τk and Mαknα=Mk. The material constants must satisfy the following inequalities, derived from the Clausius-Duhem inequality:(3a)3λ+2μ+κ0,2μ+κ0,κ0,(3b)3α+β+γ0,γβ.The momentum, angular momentum, and energy equations may be written within boundary layer approximation as(4)κΓ2x-Γ1y-P0+ρg1-βTT-T0+μ+κ2wx2+2wy2-σB02w=0,-2κΓ1+κwy-γ2Γ2xy-2Γ1y2+α+β+γ2Γ1x2+2Γ2xy=0,-2κΓ2-κwx+γ2Γ2x2-2Γ1xy+α+β+γ2Γ1xy+2Γ2x2=0,2Tx2+2Ty2=0,T0=0adx0bdyT,where w is flow velocity, Γ1 and Γ2 are components of microrotation, T is the temperature, ρ is fluid density, μ is dynamic coefficient of viscosity, κ is vortex viscosity, α, β, γ are material constants (the spin, gradient, and viscosity), g is the acceleration due to gravity, βT is the coefficient of thermal expansion, T0 is the reference (mean) temperature, B0 is magnetic induction, and σ is electric conductivity of the fluid. We introduce the nondimensional variables(5a)x=x~a,y=y~a,w=w~U0,Γ1=υ1U0a,(5b)Γ2=υ2U0a,T-T0=T1-T2θ,σ=baand the mean velocity(6)U0=1ab0adx0bwdy.Substituting the nondimensional variables in (4) and dropping the tildes, we get(7)N1-Nυ2x-υ1y-λ+GrReθ+11-N2wx2+2wy2-Ha2w=0,(8)-2υ1+wy-2-Nm22υ2xy-2υ1y2+1l22υ1x2+2υ2xy=0,(9)-2υ2-wx+2-Nm22υ2x2-2υ1xy+1l22υ1xy+2υ2y2=0,(10)2θx2+2θy2=0with the additional conditions(11)01dx0y0dyθ=00adx0y0wdy=y0,where N=κ/(μ+κ) is the coupling number, l2=2a2κ/α+β+γ is a nondimensional parameter, Gr=gβT(T2-T1)d3/ν2 is the thermal Grashof number, Re=Uod/ν is the Reynolds number, Ha2=σB02a2/ρν is the Hartman number, ν is the kinematic viscosity, m2=d2κ(2μ+κ)/γ(μ+κ) is micropolar parameter, d=2ab/a+b is the hydraulic diameter, and λ=a2Po/μUo is a dimensionless parameter.

The boundary conditions for velocity and microrotations in dimensionless form are as follows:(12a)wx,0=wx,σ=w0,y=w1,y=0,(12b)υ1x,0=υ1x,σ=υ10,y=υ11,y=0,(12c)υ2x,0=υ2x,σ=υ20,y=υ21,y=0.

The thermal boundary conditions for two cases are set as follows.

Case 1 (thermal boundary conditions).

(13a) T x x = 0 = 0 , T x , 0 , Z = T 1 , (13b) T x x = a = 0 , T x , b , Z = T 2 .

Case 2 (thermal boundary conditions).

(14a) k T y y = 0 = - q w , T 0 , y , Z = T 1 , (14b) k T y y = b = q w , T a , y , Z = T 1 . Applying the dimensionless parameters, taking the reference temperature difference ΔT=T1-T2, and assuming T1>T2 for the thermal boundary conditions (Case 1) and ΔT=aqw/k for the thermal boundary conditions (Case 2), we get the dimensionless thermal boundary conditions as follows.

Two Isothermal and Two Adiabatic Walls. Surfaces parallel to the y-axis are assumed adiabatic and temperature is assumed to be greater than the other parallel surface(15a)θxx=0=0,θxx=1=0,(15b)θx,0=s,θx,σ=s-1.

Two Isoflux and Two Isothermal Walls. Surfaces parallel to the y-axis are assumed isothermal(16a)θyy=0=-1,θyy=σ=1,(16b)θ0,y=s,θ1,y=s,where s=T1-T0/T1-T2.

3. Results and Discussion

The dimensionless temperature equation (10) and the coupled velocity and microrotation equations ((7)–(9)) together with the boundary conditions ((12a), (12b), (12c), (15a), (15b), (16a), and (16b)) are solved numerically using the finite difference method with 40 mesh points in both directions. Throughout the study unless specified otherwise the values N=0.5, l=1, m=3, λ=5, Ha=3, s=-1.6, and σ=1 are used. For validation, we use the following values: N=0, Ha=0, λ=-1, and Gr/Re0. Table 1 shows a comparison of the present values and those obtained by Barletta  and it can be seen that they are in in good agreement.

Values of w(x,y) for different positions in the duct.

(x,y) Barletta  Present
(0.5, 0.5) 0.073671 0.073560
(0.6, 0.6) 0.068744 0.068639
(0.75, 0.75) 0.045286 0.045206
(0.85, 0.85) 0.032666 0.023547
(0.9, 0.9) 0.013071 0.013026

The coupling number N characterizes the coupling of linear and rotational motion arising from the micromotion of the fluid molecules. N signifies the coupling between the Newtonian (μ) and rotational viscosities (κ) and hence 0N<1. For large values of N, the effect of microstructure becomes significant, whereas for small values of N the individuality of the substructure is much less pronounced. As κ0 (i.e., N0), the micropolarity is lost and the fluid behaves as a nonpolar fluid. Hence, N0 corresponds to a viscous fluid. More on this relation can be found. Detailed study comparing the classical and micropolar theory can be found in the work of Ansari et al. .

The results are discussed first for two of the sides (sides parallel to the y-axis) that are adiabatic and the others that are isothermal, that is, thermal boundary condition (Case 1). The profiles for velocity and microrotations v1 and v2 are presented in Figures 2(a)2(e). From Figures 2(a)2(c), it is observed that as the aspect ratio increases the flow reversal near the hot plate (plate with temperature T1) increases. The profile of microrotation v1 is antisymmetric, at y=0.4, whereas the microrotation component v2 shows antisymmetric rotation about the center for σ=1. It is also noted that microrotation components flow in opposite direction as it can be seen from Figures 2(d) and 2(e).

Case 1: different aspect ratios, with Gr/Re=50 and S=0.5.

Velocity profile

Velocity profile

Velocity profile

Microrotation V1

Microrotation V2

Figures 3(a)3(c) show the effect of Gr/Re on velocity and microrotations for the aspect ratio σ=1. For Gr/Re>0, upward flow, when mean flow U0>0, the velocity decreases with increase in Gr/Re whereas the microrotations increase. It is also noted that the microrotation component V2 shows a flow reversal near the plate as the value of Gr/Re increases. The influence of the applied magnetic field is shown in Figure 3(d) for Gr/Re=200 and aspect ratio σ=1. It is seen that the presence of the magnetic field influences the flow behaviour at the center in addition to the decreasing nature when the value of Ha increases.

Case 1: effects for thermal boundary condition.

The effect of Gr/Re on velocity

The effect of Gr/Re on microrotation V1

The effect of Gr/Re on microrotation V2

The effect of Ha effect on velocity for S=0.5 for Gr/Re=200

In the following discussion, we consider the boundary conditions where two of the sides are isothermal (sides parallel to the y-axis) and the other two sides are isoflux. The dimensionless temperature distribution (10) is solved independently of the velocity and microrotation from (10) and (16a) and (16b). The solution is shown in Figure 4(a). The profile of velocity and microrotations for different values of Gr/Re for two isoflux and two isothermal walls is shown in Figures 4(b)4(e). It is interesting to see that the flow reversals for velocity are dependent on the sign of Gr/Re with respect to the axes. It is also noted that microrotation V1 is almost constant (no rotation) around the middle of the duct as seen in Figure 4(d).

Case 2: thermal boundary condition with different values of Gr/Re.

Temperature profile

Velocity profile

Velocity profile

Microrotation

Microrotation

For Gr/Re>0, the fluid is heating (qw>0) with upward mean flow; flow reversal is observed in the neighbourhood of the isothermal walls. Figures 5(a)5(c) show that increasing the values of Gr/Re decreases the velocity. In Figures 5(d)5(f), flow reversal is observed for negative values of the parameter Gr/Re next to the isoflux walls. Microrotation components show symmetric effects in general; it is also observed that the microrotation component V2 shows reversal flow near the same plates.

Simulation of velocity and microrotation.

Velocity

Microrotation

Microrotation

Velocity

Microrotation

Microrotation

The effect of the magnetic field is presented in Figure 6(b) which shows the magnetic field slows down the flow. It is also interesting to see that the higher the applied magnetic field in addition to the flow retardation, the more pronounced the flow reversal. In Figure 6(b), the effect of s is presented. The numerical results show that when the value of s is approximately s>-0.1, the commencement of flow reversal is observed for the aspect ratio σ=1, on a square duct.

Simulation of velocity.

4. Conclusion

A model has been developed for mixed convection micropolar fluid through a vertical rectangular duct with transverse magnetic field. Effects of different parameters governing the flow have been studied numerically. Two cases for the boundary conditions have been considered: the first case with two faces with constant different temperature and the others kept adiabatic and the other case with two faces isothermal and the others isoflux. Generally, increasing the aspect ratio shows more flow reversal. It is also observed that as s-0.1, the onset of flow reversal is observed. As the magnetic field increases, flow reversal increases and velocity decreases. For microrotation components, opposite direction flow is observed. It is also seen that, for the microrotation component V1, the rotational effect is suction around the center of the duct for higher values of Gr/Re.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Eringen A. C. Microcontinuum Field Theories I: Foundations and Solids 1999 New York, NY, USA Springer-Verlag Eringen A. C. Microcontinuum Field Theories-II: Fluent Media 2001 New York, NY, USA Springer 10.1007/978-1-4612-0555-5 MR1720520 Axtell N. K. Park M. Cushman J. H. Micromorphic fluid in an elastic porous body: blood flow in tissues with microcirculation International Journal for Multiscale Computational Engineering 2005 3 1 71 84 10.1615/IntJMultCompEng.v3.i1.60 Ansari R. Bazdid-Vahdati M. Shakouri A. Norouzzadeh A. Rouhi H. Micromorphic first-order shear deformable plate element Meccanica 2016 51 8 1797 1809 2-s2.0-84948958996 10.1007/s11012-015-0325-7 Lee J. D. Wang X. Generalized micromorphic solids and fluids International Journal of Engineering Science 2011 49 12 1378 1387 MR2893957 10.1016/j.ijengsci.2011.04.001 2-s2.0-80053983329 Ansari R. Bazdid-Vahdati M. Shakouri A. H. Norouzzadeh A. Rouhi H. Micromorphic prism element Mathematics and Mechanics of Solids 2017 22 6 1438 1461 MR3659624 10.1177/1081286516637115 Eringen A. C. Simple microfluids International Journal of Engineering Science 1964 2 205 217 MR0169468 10.1016/0020-7225(64)90005-9 2-s2.0-49749223867 Eringen A. C. Theory of micropolar fluids Journal of Applied Mathematics and Mechanics 1966 16 1 18 MR0204005 Eringen A. C. Theory of thermomicrofluids Journal of Mathematical Analysis and Applications 1972 38 2 480 496 10.1016/0022-247X(72)90106-0 2-s2.0-33847495482 Ariman T. Sylvester N. D. Turk M. A. Applications of microcontinuum fluid mechanics International Journal of Engineering Science 1974 12 4 273 293 10.1016/0020-7225(74)90059-7 2-s2.0-0016046496 Eremeyev V. A. Lebedev L. P. Altenbach H. Foundations of Micropolar Mechanics 2013 Springer Sci and Business Media MR2963524 Barletta A. Analysis of flow reversal for laminar mixed convection in a vertical rectangular duct with one or more isothermal walls International Journal of Heat and Mass Transfer 2001 44 18 3481 3497 2-s2.0-0035452123 10.1016/S0017-9310(01)00016-3 Ahmed S. E. Numerical study of MHD natural convection in an inclined rectangular cavity with internal heat generation filled with a porous medium under the influence of joule heating Latin American Applied Research 2013 43 67 71 Spiga M. Morino G. L. A symmetric solution for velocity profile in laminar flow through rectangular ducts International Communications in Heat and Mass Transfer 1994 21 4 469 475 2-s2.0-0028468927 10.1016/0735-1933(94)90046-9 Srnivasacharya D. Shiferaw M. MHD flow of micropolar fluid in a rectangular duct with hall and ion slip effects Journal of the Brazilian Society of Mechanical Sciences and Engineering 2008 4 313 10.1590/S1678-58782008000400007 Chu H. H. S. Churchill S. W. Patterson C. V. S. The effect of heater size, location, aspect ratio, and boundary conditions on two-dimensional, laminar, natural convection in rectangular channels Journal of Heat Transfer 1976 98 Sayed-Ahmed M. E. Attia H. A. The effect of Hall current on magnetohydrodynamic flow and heat transfer for Bingham fluids in a rectangular duct Canadian Journal of Physics 2005 83 6 637 651 2-s2.0-27544434542 10.1139/p05-006 Chou F. C. Hwang G. J. Vorticity-velocity method for Graetz problem with the effect of natural convection in a horizontal rectangular channel with uniform wall heat flux Journal of Heat Transfer 1987 109 704 710 Ramakrishna K. Rubin S. G. Khosla P. K. Laminar natural convection along vertical square ducts, Numer Heat Transfer 1982 5 59 79 Cheng K. C. Ou J. W. Convection instability and finite amplitude convection in the entrance region of horizontal rectangular channels heated from below Proceedings of 7th International Heat Transfer Conference 1982 2 189 194 2-s2.0-0020279891 Aung W. Worku G. Mixed convection in ducts with asymmetric wall heat fluxes ASME Journal of Heat Transfer 1987 109 947 951 Mahaney H. V. Incropera F. P. Ramadhyani S. Development of laminar mixed convection flow in a horizontal rectangular duct with uniform bottom heating Numerical Heat Transfer, Part B: Fundamentals 1987 12 2 137 155 2-s2.0-0023175601 10.1080/10407788708913578 Abou-Ellail M. M. Morcos S. M. Buoyancy effects in the entrance region of horizontal rectangular channels Journal of Heat Transfer 1983 105 924 928 Huang C.-C. Yan W.-M. Jang J.-H. Laminar mixed convection heat and mass transfer in vertical rectangular ducts with film evaporation and condensation International Journal of Heat and Mass Transfer 2005 48 9 1772 1784 2-s2.0-17644381539 10.1016/j.ijheatmasstransfer.2004.11.023 Yan W. M. Lin T. F. Effects of wetted wall on laminar mixed convection in a vertical channel Journal of Thermophysics and Heat Transfer 1989 3 1 94 96 2-s2.0-0024482246 10.2514/3.56231 Yan W.-M. Turbulent mixed convection heat and mass transfer in a wetted channel Journal of Heat Transfer 1995 117 1 229 233 2-s2.0-0029135625 10.1115/1.2822311 Hayat T. Nadeem S. Asghar S. Hydromagnetic Couette flow of an Oldroyd-B fluid in a rotating system International Journal of Engineering Science 2004 42 1 65 78 MR2008507 10.1016/S0020-7225(03)00277-5 2-s2.0-0141860819 Ansari R. Shakouri A. H. Bazdid-Vahdati M. Norouzzadeh A. Rouhi H. A nonclassical finite element approach for the nonlinear analysis of micropolar plates Journal of Computational and Nonlinear Dynamics 2017 12 1 2-s2.0-84999751834 10.1115/1.4034678 011019