# Similarity Solution for Free Convection Flow of a Micropolar Fluid under Convective Boundary Condition via Lie Scaling Group Transformations

Authors: Ch. RamReddy; T. Pradeepa; D. Srinivasacharya
Journal: Advances in High Energy Physics
 (2015)

---

## Abstract

The free convective
flow of an incompressible micropolar
fluid along permeable vertical
plate under the convective boundary condition is investigated. The Lie scaling group of
transformations is applied to get the similarity representation for the system of partial differential equations and then the resulting systems of equations are solved using spectral quasi-linearisation method. A quantitative comparison of the numerical results is made
with previously published results for special cases and the results are found to be in good
agreement. The results of the physical parameters on the developments of
flow, temperature,
concentration, skinfriction, wall couple stress, heat transfer, and mass transfer characteristics
along vertical plate are given and the salient features are discussed.


---

## Body

## 1. Introduction

In the past few decades, most of the researchers considered convective heat transfer problems with either constant wall temperature (CWT), constant heat flux (CHF), or Newtonian heating (NH) in a Newtonian and/or non-Newtonian fluid. Recently, a novel mechanism for the heating process has drawn the involvement of many researchers, namely,convective boundary condition (CBC), where the heat is supplied to the convecting fluid through a bounding surface with a finite heat capacity. Further, this results in the heat transfer rate through the surface being proportional to the local difference in temperature with the ambient conditions (Merkin [1]). Besides, it is more general and realistic, particularly in various technologies and industrial operations such as transpiration cooling process, textile drying, and laser pulse heating. Aziz [2] reported similarity solution for thermal boundary layer flow over a flat plate in a uniform stream of fluid with the convective boundary condition and he concluded that a similarity solution is possible if the convective heat transfer related to hot fluid on the lower surface of the plate is proportional to the inverse square root of the axial length. In the presence of an internal heat generation local similarity solution for free convection heat transfer from a moving vertical plate with the convective boundary condition is discussed by Makinde [3]. The laminar natural convection flow over a semi-infinite moving vertical plate under the convective boundary condition is examined by Ibrahim and Bhashar Reddy [4]. RamReddy et al. [5] investigated the influence of the prominent Soret effect on mixed convection in a nanofluid under the convective boundary conditions. The nonsimilar result has been presented for the free convection boundary layer flow along a solid sphere under the convective boundary conditions by Alkasasbeh et al. [6]. More recently, a note on the natural convection along convectively heated vertical plate is given by Pantokratoras [7].One of the best established theories of fluids with microstructure is the theory of micropolar fluids and this theory can be found in the books by Lukaszewicz [8] and Eremeyev et al. [9]. It has gathered a good deal of attention due to the obvious reasons that the Navier Stokes equation for Newtonian fluids cannot successfully explain the attributes of fluids with a substructure. Physically, the micropolar fluids may be treated as non-Newtonian fluids consisting of dumb-bell molecules or rigid cylindrical element, polymer fluids, fluid suspension, animal blood, and so forth. Further, the theory of micropolar fluids includes microrotation as well as microinertia effects. This theory studies viscous fluids in which microconstituents are rigid and spherical or randomly oriented as well. The subject of free convection boundary layer flow in a micropolar fluid has been keyed out by several investigators due to its immense applications in the engineering problems such as solar energy collecting devices, air conditioning of a room, material processing, and passive cooling of nuclear reactors. The boundary layer flow over a semi-infinite flat plate is considered for analyzing theory of micropolar fluid and its application to low concentration suspension flow by Ahmadi [10]. Rees and Pop [11] discussed the free convection boundary layer flow of a micropolar fluid from a vertical flat plate. The nonsimilarity transformations are used to analyze the effects of double stratification on free/mixed convective transport in a micropolar fluid by Srinivasacharya and RamReddy [12–14] (also see the references cited therein). The problems of a steady laminar stagnation point flow towards a stretching/shrinking sheet in an incompressible micropolar fluid under the convective surface boundary condition are discussed by Yacob and Ishak [15] and Zaimi and Ishak [16]. Merely from the literature, it is noted that the majority of the researchers have found the local similarity or nonsimilarity solutions for the problems involvingconvective boundary conditions, since most of the researchers have taken a convective heat transfer coefficient as a function x for getting the similarity solutions in their problems. Nevertheless, the assumption of a heat transfer coefficient varying along the plate as a function of x is not realistic and very difficult to be obtained in practice. For that cause, it could be supposed that the above works have only theoretical value.In the recent past, several researchers are focused on obtaining the similarity solutions of the convective transport phenomena problems arising in fluid dynamics, aerodynamics, plasma physics, meteorology, and some branches of engineering by using different procedures. One such procedure isLie group analysis. The concept of Lie group analysis also called symmetry analysis is developed bySophius Lie to determine transformations which map a given differential equation to itself and it unifies almost all known exact integration techniques (see [17–19]). It provides a potent, sophisticated, and systematic tool for generating the invariant solutions of the system of nonlinear partial differential equations (PDEs) with relevant initial or boundary conditions. A special form of Lie group transformations, known as the scaling group, has been suggested by various researchers to study convection flows of different flow phenomena (see Tapanidis et al. [20], Hassanien and Hamad [21], Kandasamy et al. [22], Aziz et al. [23], Mutlag et al. [24], etc.; they are worth observing).From the literature survey, it seems that the problem of the free convective heat and mass transport along permeable vertical plate in a micropolar fluid under the convective boundary condition has not been investigated so far. Motivated by all these works, this paper attempts to present the new similarity transformations and corresponding similarity solution to investigate the free convection flow of a micropolar fluid under the convective boundary condition using the Lie group transformations. The mathematical model involving the convective boundary conditions becomes slightly more complicated leading to the complex interactions of the flow, heat, and mass transfer mechanism. Further, the analytical solution is out of scope in the present set-up and hence a numerical solution is obtained for the current problem. Also, the influence of important parameters, namely, micropolar, suction/injection, and convective heat transfer parameters, on the physical quantities of the flow, heat, and mass transfer rates is analyzed in different flow situations.


## 2. Mathematical Formulation

Consider the steady, laminar, and free convective flow of an incompressible micropolar fluid with the free stream temperature and concentration,T∞ and C∞, respectively. Choose the coordinate system such that the x¯-axis is along the vertical plate and y¯-axis normal to the plate, as shown in Figure 1. The suction/injection velocity distribution is assumed to be vw. The plate is either heated or cooled from left by convection from a fluid of temperature Tf with Tf>T∞ corresponding to a heated surface (assisting flow) and Tf<T∞ corresponding to a cooled surface (opposing flow), respectively. On the wall concentration is taken to be constant and is given by Cw.Figure 1
Physical model and coordinate system.By employing Boussinesq approximation and making use of the standard boundary layer approximations, the governing equations for the micropolar fluid [10] are given by(1)∂u¯∂x¯+∂v¯∂y¯=0,(2)ρu¯∂u¯∂x¯+v¯∂u¯∂y¯=μ+κ∂2u¯∂y¯2+κ∂ω¯∂y¯+ρg∗βTx¯T-T∞+βCx¯C-C∞,(3)ρju¯∂ω¯∂x¯+v¯∂ω¯∂y¯=γ∂2ω¯∂y¯2-κ2ω¯+∂u¯∂y¯,(4)u¯∂T∂x¯+v¯∂T∂y¯=α∂2T∂y¯2,(5)u¯∂C∂x¯+v¯∂C∂y¯=D∂2C∂y¯2,where u¯ and v¯ are the velocity components in x¯ and y¯ directions, respectively, ω¯ is the component of microrotation whose direction of rotation lies in the x¯y¯-plane, T is the temperature, C is the concentration, g∗ is the acceleration due to gravity, ρ is the density, μ is the dynamic coefficient of viscosity, βT(x¯) is the volumetric coefficient of thermal expansion, βC(x¯) is the volumetric coefficient of solutal expansions, κ is the vortex viscosity, j is the microinertia density, γ is the spin-gradient viscosity, α is the thermal diffusivity, and D is the solutal diffusivity of the medium.The boundary conditions are(6a)u¯=0,v¯=vw,ω¯=-n∂u¯∂y¯,-k∂T∂y¯=hfTf-T,C=Cw,aty¯=0,(6b)u¯=0,ω¯=0,T=T∞,C=C∞,asy¯⟶∞,where subscripts w and ∞ indicate the conditions at the wall and at the outer edge of the boundary layer, respectively, hf is the convective heat transfer coefficient, k is the thermal conductivity of the fluid, and n is a material constant. Further, we follow the work of many recent authors by assuming that γ=μ+κ/2j. This assumption is invoked to allow the field of equations to predict the correct behavior in the limiting case when the microstructure effects become negligible and the total spin ω¯ reduces to the angular velocity [10].


## 3. Nondimensionalization of the Governing Equations

Introduce the following dimensionless variables:(7)x=x¯L,y=y¯LGr1/4,u=LνGr1/2u¯,v=LνGr1/4v¯,ω=L2νGr3/4ω¯,θ=T-T∞Tf-T∞,ϕ=C-C∞Cw-C∞,where Gr=g∗βT0(Tf-T∞)L3/ν2 is the Grashof number.In view of the continuity equation (1), we introduce the stream function ψ by(8)u=∂ψ∂y,v=-∂ψ∂x.Using (7) and (8) into (2)–(5), we get the following momentum, angular momentum, energy, and concentration equations: (9)Δ1=∂ψ∂y∂2ψ∂x∂y-∂ψ∂x∂2ψ∂y2-11-N∂3ψ∂y3-N∂ω∂y-g∗βTxTf-T∞L3ν2Grθ-g∗βCxCw-C∞L3ν2Grϕ=0,Δ2=∂ψ∂y∂ω∂x-∂ψ∂x∂ω∂y-2-N2-2N∂2ω∂y2+N1-N2ω+∂2ψ∂y2=0,Δ3=∂ψ∂y∂θ∂x-∂ψ∂x∂θ∂y-1Pr∂2θ∂y2=0,Δ4=∂ψ∂y∂ϕ∂x-∂ψ∂x∂ϕ∂y-1Sc∂2ϕ∂y2=0.In usual definitions, ν is the kinematic viscosity, Pr=ν/α is the Prandtl number, Sc=ν/D is the Schmidt number, N=κ/(μ+κ) (0≤N<1) is the coupling number [25], and the microinertia density is taken to be j=L2/Gr1/2.Now boundary conditions (6a) and (6b) become (10a)∂ψ∂y=0,∂ψ∂x=fw,ω=-n∂2ψ∂y2,∂θ∂y=-Bi1-θ,ϕ=1,aty=0,(10b)∂ψ∂y=0,ω=0,θ=0,ϕ=0,asy⟶∞,where fw=-(L/νGr1/4)vw is the suction/injection parameter. It is worth mentioning that fw determines the transpiration rate at the surface, with fw>0 for suction and fw<0 for injection, and fw=0 corresponds to an impermeable surface. Further, Bi=hfL/kGr1/4 is the Biot number. It is a ratio of the internal thermal resistance of the plate to the boundary layer thermal resistance of the hot fluid at the bottom of the surface.


## 4. Similarity Equations via Lie Scaling Group Transformations

A one-parameter Lie scaling group of transformations, which is a simplified form of Lie group transformation, is selected as (for more, see [26–31])(11)Γ:x∗=xeɛα1,y∗=yeɛα2,ψ∗=ψeɛα3,ω∗=ωeɛα4,θ∗=θeɛα5,ϕ∗=ϕeɛα6,βT∗=βTeɛα7,βC∗=βCeɛα8.Here ɛ≠0 is the parameter of the group and αi (where i=1,2,3,…,8) are arbitrary real numbers whose interrelationship will be determined by our analysis. Transformations in (11) may be treated as a point transformation, transforming the coordinates(12)x,y,ψ,ω,θ,ϕ,βT,βC=x∗,y∗,ψ∗,ω∗,θ∗,ϕ∗,βT∗,βC∗.We now investigate the relationship among the exponents αi (where i=1,2,3,…,8) such that(13)Δjx∗,y∗,u∗,v∗,…,∂3ψ∗∂y∗3=Hjx,y,u,v,…,∂3ψ∂y3;a·Δjx,y,u,v,…,∂3ψ∂y3,j=1,2,3,4.This is the requirement that the differential forms Δ1, Δ2, Δ3, and Δ4 are conformally invariant under transformation (11). Substituting transformations (11) in (9), we have (14a)Δ1=eɛα1+2α2-2α3∂ψ∗∂y∗∂2ψ∗∂x∗∂y∗-∂ψ∗∂x∗∂2ψ∗∂y∗2-11-Neɛ3α2-α3∂3ψ∗∂y∗3-N1-Neɛα2-α4∂ω∗∂y∗-g∗βT∗Tf-T∞L3ν2Gre-ɛα5+α7θ∗-g∗βC∗Cw-C∞L3ν2Gre-ɛα6+α8ϕ∗=0,(14b)Δ2=eɛα1+α2-α3-α4∂ψ∗∂y∗∂ω∗∂x∗-∂ψ∗∂x∗∂ω∗∂y∗-2-N2-2Neɛ2α2-α4∂2ω∗∂y∗2+N1-N2ω∗e-ɛα4+eɛ2α2-α3∂2ψ∗∂y∗2=0,(14c)Δ3=eɛα1+α2-α3-α5∂ψ∗∂y∗∂θ∗∂x∗-∂ψ∗∂x∗∂θ∗∂y∗-1Preɛ2α2-α5∂2θ∗∂y∗2=0,(14d)Δ4=eɛα1+α2-α3-α6∂ψ∗∂y∗∂ϕ∗∂x∗-∂ψ∗∂x∗∂ϕ∗∂y∗-1Sceɛ2α2-α6∂2ϕ∗∂y∗2=0.Now, boundary conditions (10a) and (10b) become (15a)eɛα2-α3∂ψ∗∂y∗=0,eɛα1-α3∂ψ∗∂x∗=fw,e-ɛα4ω∗=-neɛ2α2-α3∂2ψ∗∂y∗2,eɛα2-α5∂θ∗∂y∗=-Bi1-e-ɛα5θ∗,e-ɛα6ϕ∗=1,aty∗=0,(15b)eɛα2-α3∂ψ∗∂y∗=0,e-ɛα4ω∗=0,e-ɛα5θ∗=0,e-ɛα6ϕ∗=0,asy∗⟶∞.The system remains invariant under the group transformation Γ. We then have the following relationships for the parameters:(16)α1+2α2-2α3=3α2-α3=α2-α4=-α5-α7=-α6-α8;α1+α2-α3-α4=2α2-α4=-α4=2α2-α3;α1+α2-α3-α5=2α2-α5;α1+α2-α3-α6=2α2-α6;α1-α3=0;-α4=2α2-α3;α2-α5=0=-α5;α6=0.Solving linear system (16), we have the following relationship among the exponents:(17)α1=α3=α4=α7=α8,α2=α5=α6=0.The set of transformations Γ reduces to(18)x∗=xeɛα1,y∗=y,ψ∗=ψeɛα1,ω∗=ωeɛα1,θ∗=θ,ϕ∗=ϕ,βT∗=βTeɛα1,βC∗=βCeɛα1.Expanding by the Taylor series in power of ɛ, keeping the term up to the first degree (neglecting higher power of ɛ), we get(19)x∗-x=ɛα1x,y∗=y,ψ∗-ψ=ɛα1ψ,ω∗-ω=ɛα1ω,θ∗=θ,ϕ∗=ϕ,βT∗-βT=ɛα1βT,βC∗-βC=ɛα1βC.The characteristic equations are(20)dxα1x=dy0=dψα1ψ=dωα1ω=dθ0=dϕ0=dβTα1βT=dβCα1βC.Solving the above characteristic equations, we have the following similarity transformations:(21)η=y,ψ=xfη,ω=xgη,βT=βT0x,βC=βC0x,θ=θη,ϕ=ϕη,where βT0 and βC0 are constant thermal and mass coefficient of expansion.Using (21) into (9), we get the following similarity equations:(22)11-Nf′′′+ff′′-f′2+N1-Ng′+θ+Bϕ=0,2-N2-2Ng′′+fg′-f′g-N1-N2g+f′′=0,1Prθ′′+fθ′=0,1Scϕ′′+fϕ′=0,where the primes indicate differentiation with respect to η alone and B=βC0(Cw-C∞)/βT0(Tf-T∞) is the buoyancy ratio.Boundary conditions (10a) and (10b) in terms of f, g, θ, and ϕ become (23a)η=0:f0=fw,f′0=0,g0=-nf′′0,θ′0=-Bi1-θ0,ϕ0=1,(23b)η⟶∞:f′∞=0,g∞=0,θ∞=0,ϕ∞=0.


## 5. Skin Friction, Wall Couple Stress, and Heat and Mass Transfer Coefficients

The wall shear stress and the wall couple stress are(24a)τw=μ+κ∂u¯∂y¯+κω¯y¯=0,mw=γ∂ω¯∂y¯y¯=0,and the heat and mass transfers from the plate, respectively, are given by (24b)qw=-k∂T∂y¯y¯=0,qm=-D∂C∂y¯y¯=0.The nondimensional skin frictionCf=2τw/ρu¯∗2, wall couple stress Mw=mw/ρu¯∗2x¯, the local Nusselt number Nux¯=qwx¯/k(Tf-T∞), and local Sherwood number Shx¯=qmx¯/D(Cw-C∞) are given by(25)CfGrx¯1/4=21-nN1-Nf′′0,MwGrx¯1/2=2-N2-2Ng′0,Nux¯Grx¯1/4=-θ′0,Shx¯Grx¯1/4=-ϕ′0,where u¯∗2 is the characteristic velocity and Grx¯=g∗βT0(Tf-T∞)x¯3/ν2 is the local Grashof number.


## 6. Numerical Solution Using the Spectral Quasi-Linearization Method (SQLM)

In this section, we describe the quasi-linearization method (QLM) for solving the governing system of (22) along with boundary conditions (23a) and (23b). This QLM is a generalization of the Newton-Raphson method and was proposed by Bellman and Kalaba [32] for solving nonlinear boundary value problems.Assume that the solutionsfr, gr, θr, and ϕr of (22) at the (r+1)th iteration are fr+1, gr+1, θr+1, and ϕr+1. If the solutions at the previous iteration are sufficiently close to the solutions at the present iteration, the nonlinear components of (22) can be linearised using one-term Taylor series of multiple variables so that (22) give the following iterative sequence of linear differential equations:(26)11-Nfr+1′′′+a1,rfr+1′′+a2,rfr+1′+a3,rfr+1+N1-Ngr+1′+θr+1+Bϕr+1=R1,r,2-N2-2Ngr+1′′+b3,rgr+1′+b4,rgr+1+b1,rfr+1′+b2,rfr+1-N1-Nfr+1′′=R2,r,c1,rfr+1+1Prθr+1′′+c2,rθr+1′=R3,r,d1,rfr+1+1Scϕr+1′′+d2,rϕr+1′=R4,r,where the coefficients as1,r(s1=1,2,3), bs2,r(s2=1,2,…,4), cs3,r(s3=1,2), ds4,r(s4=1,2), and Rs5,r(s5=1,2,…,4) are known functions (from previous iterations) and are defined as(27)a1,r=fr,a2,r=-2fr′,a3,r=fr′′,R1,r=frfr′′-fr′2,b1,r=-gr,b2,r=gr′,b3,r=fr,b4,r=-fr′-2N1-N,R2,r=frgr′-fr′gr,c1,r=θr′,c2,r=fr,R3,r=frθr′,d1,r=ϕr′,d2,r=fr,R4,r=frϕr′,subject to boundary conditions(28)fr+10=fw,fr+1′=0,fr+1′∞=0,gr+1=-nfr+1′′0,gr+1∞=0,θr+1′0=-Bi1-θ0,θr+1∞=0,ϕr+10=1,ϕr+1∞=0.System (26) constitutes a linear system of coupled differential equations with variable coefficients and can be solved iteratively using any numerical method for r=1,2,3,…. In this work, as will be discussed below, the Chebyshev pseudospectral method was used to solve the QLM scheme (26) (for more details, refer to the works of Motsa et al. [33, 34]):(29)f0η=fw+1-e-η,g0η=-ne-η,θ0=e-ηBiBi+1,ϕ0=e-η,and starting from these sets of initial approximations f0, g0, θ0, and  ϕ0, the iteration schemes (26) can be solved iteratively for fr+1(η),gr+1(η),θr+1(η),  and  ϕr+1(η) when r=0,1,2,…. For this, we discretise the equation using the Chebyshev spectral collocation method. The unknown functions are approximated by the Chebyshev interpolating polynomials in such way that they are collocated at the Gauss-Lobatto collocation points defined as(30)τj=cos⁡πjN¯,j=0,1,2,…,N¯,where N¯ is the number of collocation points. The physical region [0,∞) is transformed into the region [-1,1] using the domain truncation technique in which the problem is solved on the interval [0,η∞] instead of [0,∞). This leads to the mapping(31)ηη∞=τ+12,-1≤τ≤1,where η∞ is the scaling parameter used to invoke the boundary condition at infinity. The functions f, g, θ, and ϕ are approximated at the collocation points by(32)fτ=∑k=0N¯fτkTkτj,gτ=∑k=0N¯gτkTkτj,θτ=∑k=0N¯θτkTkτj,ϕτ=∑k=0N¯ϕτkTkτj,j=0,1,2,…,N¯,where Tk is the kth Chebyshev polynomial defined as(33)Tkτ=cos⁡kcos-1⁡τ.The derivatives of the variables at the collocation points are represented as(34)dpfdηp=∑k=0N¯Dlkpfτk,dpgdηp=∑k=0N¯Dlkpgτk,dpθdηp=∑k=0N¯Dlkpθτk,dpϕdηp=∑k=0N¯Dlkpϕτk,l=0,1,…,N¯,where p is the order of the derivative and D=2D/η∞ is the Chebyshev spectral differentiation matrix and its entries are clearly defined in Canuto et al. [35].Substituting (31)–(34) into (26) leads to the matrix equation(35)AX=R,subject to the boundary conditions(36)fr+1τN¯=fw,∑k=0N¯DN¯kfτk=0,∑k=0N¯D0kfτk=0,gr+1τN¯=-n∑k=0N¯DN¯k2fτk,gr+1τ0=0,∑k=0N¯DN¯kθr+1τk-Biθr+1τN¯=-Bi,θr+1τ0=0,ϕr+1τN¯=1,ϕr+1τ0=0.In (35) A is a (4N¯+4)×(4N¯+4) square matrix and X and R are (4N¯+1)×1 column vectors defined by(37)A=A11A12A13A14A21A22A23A24A31A32A33A34A41A42A43A44,X=Fr+1Gr+1Θr+1Φr+1,R=R1R2R3R4,where(38)F=fr+1τ0,fr+1τ1,…,fr+1τN¯T,G=gr+1τ0,gr+1τ1,…,gr+1τN¯T,Θ=θr+1τ0,θr+1τ1,…,θr+1τN¯T,Φ=ϕr+1τ0,ϕr+1τ1,…,ϕr+1τN¯T,A11=11-ND3+diag⁡a1,rD2+diag⁡a2,rD+diag⁡a3,r,A12=N1-ND,A13=I,A14=BI,A21=-N1-ND2+diag⁡b1,rD+diag⁡b2,r,A22=2-N2-2ND2+diag⁡b3,rD+diag⁡b4,r,A23=0,A24=0,A31=diag⁡c1,r,A32=0,A33=1PrD2+diag⁡c2,rD,A34=0,A41=diag⁡d1,r,A42=0,A43=0,A44=1ScD2+diag⁡d2,rD,R1=R1,r,R2=R2,r,R3=R3,r,R4=R4,r,and here I is an identity matrix, the size of the matrix 0 is (N¯+1)×1, and diag is a diagonal matrix of size (N¯+1)×(N¯+1). Subscript r denotes the iteration number. After modifying matrix system (35) to incorporate boundary condition (36), the solution is obtained as(39)X=A-1R.


## 7. Results and Discussions

It is noticed that the present problem reduces to free convection heat transfer along an impermeable vertical plate in a micropolar fluid without the convective boundary condition whenfw=0, Bi→∞, and B=0. Also in the limit as N→0, governing equations (2)–(5) reduce to the corresponding equations for a free convection heat and mass transfer in a viscous fluid. In order to validate the code generated the results of the present problem have been compared with the results obtained by Merkin [36], Nazar et al. [37], and Molla et al. [38] as a special case by taking N=0, n=0, B=0, Pr=1, Bi→∞, and fw=0 and it was found that they are in good agreement, as presented in Table 1. Also, the comparison of heat transfer coefficient has been made with the results obtained by Nazar et al. [37] as shown in Table 2 when n=0.5, B=0, Pr=1, Bi→∞, and fw=0. To study the effects of coupling number N, suction/injection parameter fw, Biot number Bi, and material parameter n, computations were carried out in the cases of B=1.0, Pr=0.71, and Sc=0.22.Table 1
Comparison of-θ′(0) for free convection along a vertical flat plate in Newtonian fluid when N=0, n=0, B=0, Pr⁡=1, Bi→∞, and fw=0.
Merkin [36]
Nazar  et al. [37]
Molla et al. [38]
Present
0.4214
0.4214
0.4214
0.4214313Table 2
Comparison of-θ′(0) for free convection flow in a micropolar fluid obtained by Nazar  et al.  [37] when n=0.5, B=0, Pr⁡=1, Bi→∞, and fw=0.
N
Nazar  et al.  [37]
Present
0.00
0.4214
0.4214
0.33
0.3991
0.3990
0.50
0.3834
0.3834
0.60
0.3709
0.3709
0.66
0.3608
0.3608
0.71
0.3522
0.3522
0.75
0.3447
0.3447The effects of the coupling numberN on the dimensionless velocity, microrotation, temperature, and concentration are illustrated in Figures 2(a)–2(d) with fixed values of other parameters. The coupling number N characterizes the coupling of linear and rotational motion arising from the fluid particles. In the case of N=0 (i.e., as κ tends to zero) the micropolarity is absent and fluid becomes nonpolar fluid. With a large value of N effect of microstructure becomes significant, whereas with a diminished value of N the individuality of the substructure is much less articulated. As N increases, it is found from Figure 2(a) that the maximum velocity decreases in amplitude and the location of the maximum velocity moves farther away from the wall. Since N→0 corresponds to viscous fluid, the velocity in case of a micropolar fluid has been less compared to that of viscous fluid case. It can be observed from Figure 2(b) that, as N increases, initially microrotation profiles tend to become flatter and then approach their free stream values far away from the wall. This happens due to the vanishing of antisymmetric part of the stress on the boundary that corresponds to a weak concentration of microelements. This is because an increment in the value of N implies a higher vortex viscosity of fluid which promotes the microrotation of micropolar fluids. It is seen from Figures 2(c) and 2(d) that thermal and concentration boundary layers of the fluid increase with increase in coupling number N. Hence, temperature and concentration in case of micropolar fluids are more than those of the viscous fluid case.Figure 2
Effect ofN on (a) velocity, (b) microrotation, (c) temperature, and (d) concentration profiles.
(a)
(b)
(c)
(d)The Biot numberBi is the ratio of internal thermal resistance of a solid to boundary layer thermal resistance. When Bi=0 the plate is totally insulated, internal thermal resistance of the plate is extremely high, and no convective heat transfer to the cold fluid on the upper part of the plate takes place. Figure 3(a) depicts fluid velocity profiles for different values of the Biot number with N=0.5,fw=0.5, and n=0.5. Generally, fluid velocity is zero at plate surface and increases gradually away from plate to free stream value satisfying boundary conditions. It is interesting to note that an increase in the intensity of convective surface heat transfer Bi produces significant enhancement in fluid velocity within the momentum boundary layer. In Figure 3(b), we bring out the behavior of microrotation with different values of Biot number Bi for fixed values of other parameters. As the parameter value Bi increases microrotation showing reverse rotation near the two boundaries. Hence, the condition of vanishing of antisymmetric part of the stress on the boundary results in a drastic change of the microrotation profiles. Given that convective heating increases with Biot number, Bi→∞ simulates the isothermal surface which is clearly seen from Figure 3(c), where θ(0)=1 as Bi→∞. In fact, a high Biot number indicates higher internal thermal resistance of the plate than boundary layer thermal resistance. In this fluid temperature is maximum at the plate surface and decreases exponentially to zero value far out from the plate satisfying boundary conditions. As a consequence, an increment in the Biot number leads to increase of fluid temperature efficiency. Figure 3(d) illustrates the variation of dimensionless concentration for different values of Bi. It is clear that the concentration of fluid decreases with increase of Bi.Figure 3
Effect ofBi on (a) velocity, (b) microrotation, (c) temperature, and (d) concentration profiles.
(a)
(b)
(c)
(d)The effect offw on velocity profile is depicted in Figure 4(a). Here, fw>0 represents suction and fw<0 denotes injection. The lower velocity is noticed in case of suction when compared to case of injection. From Figure 4(b), we note that microrotation is showing reverse rotation near two boundaries with both suction and injection parameter. The dimensionless temperature for different values of suction/injection parameters is drawn in Figure 4(c). It is readable that the temperature of the fluid is more in case of injection, whereas it is less in case of suction in comparison with the impermeable surface case (fw=0). Figure 4(d) demonstrates dimensionless concentration for different values of suction/injection parameters. It is determined that the concentration of fluid is more with injection, whereas it is less with suction when compared to the impermeable surface case (fw=0). As a finale, the thermal and solutal boundary layer thicknesses increase in case of injection compared to case of suction as displayed in Figures 4(c) and 4(d).Figure 4
Effect offw on (a) velocity, (b) microrotation, (c) temperature, and (d) concentration profiles.
(a)
(b)
(c)
(d)In Figures5(a)–5(d), the effect of material parameter n on dimensionless velocity, microrotation, temperature, and concentration is presented for fixed values of other parameters, since the material parameter n signifies the microrotation effects (i.e., for n=0, particles are not free to revolve near the surface whereas, as n increases from 0 to 1, the microrotation term gets augmented and induces flow enhancement). As n increases, it is found from Figure 5(a) that the minimum velocity increases in amplitude and the location of the minimum velocity moves farther away from the wall. From Figure 5(b), we observe that the microrotation is decreasing with increasing value of material parameter n within the boundary layer. It is clear from Figures 5(c)-5(d) that the thermal and solutal boundary layer thicknesses decrease with increase of material parameter n.Figure 5
Effect of material parametern on (a) velocity, (b) microrotation, (c) temperature, and (d) concentration profiles.
(a)
(b)
(c)
(d)The variations of-θ′(0) and -ϕ′(0) versus coupling number N are shown in Figures 6–8. It can be noticed from these figures that the heat and mass transfer coefficients are less in case of micropolar fluids when compared to the viscous fluids. This is because as N increases, the thermal and solutal boundary layer thicknesses become larger, thus giving rise to a small value of local heat and mass transfer rates. The effect of the Biot number Bi with fixed fw=0.5 and n=0.5 on local heat transfer coefficient is exhibited in Figure 6. It is found from Figure 6 that local heat transfer rate increases nonlinearly with the increase in Biot number Bi. The influence of the Biot number Bi on local mass transfer coefficient is shown in Figure 6. Figure 6 reveals that the local mass transfer coefficient is enhanced by the growth in the Biot number Bi. In Figure 7, the effect of the suction/injection parameter fw with fixed Bi=0.1 and n=0.5 on local heat and mass transfer coefficients is displayed. It is found from Figure 7 that the local heat and mass transfer coefficients are less in the case of injection fw<0 in comparison with the case of suction fw>0. Figure 8 is prepared to analyze the effect of the material parameter n with fixed Bi=0.1 and fw=0.5 on local heat and mass transfer coefficients. Figure 8 reveals that the local heat and mass transfer coefficients are enhanced by the increase in material parameter n. This is because when n increases from 0 to 1, the microrotation term gets augmented and induces flow enhancement.Figure 6
Effect ofBi on (a) heat transfer rate and (b) mass transfer rate.
(a)
(b)Figure 7
Effect offw on (a) heat transfer rate and (b) mass transfer rate.
(a)
(b)Figure 8
Effect of material parametern on (a) heat transfer rate and (b) mass transfer rate.
(a)
(b)The variations ofCfGrx¯1/4 and MwGrx¯1/2, which are proportional to the coefficients of skin friction and wall couple stress, are shown in Table 3 with different values of the coupling number N for fixed n=0.5, fw=0.5, and Bi=0.1. It indicates that the skin friction factor is higher for micropolar fluid than the viscous fluids (N=0). Since micropolar fluids offer a heavy resistance (resulting from vortex viscosity) to fluid movement and cause larger skin friction factor compared to viscous fluid, the results as well suggest larger values of coupling number N and lower wall couple stresses. Since the skin friction coefficient f′′(0) and wall couple stress coefficient as well as high temperature and mass transport rates are more depressed in the micropolar fluid comparing to the viscous fluid, which may be beneficial in flow, temperature, and concentration control of polymer processing, thus, the presence of microscopic effects arising from the local structure and of the fluid elements reduces the high temperature and mass transfer coefficients. The effect of Bi on CfGrx¯1/4 and MwGrx¯1/2 for fw=0.5, n=0.5, and N=0.5 is illustrated in Table 3. It can be noticed that the skin friction and wall couple stress coefficients are increasing with increase of Bi for fixed values of other parameters. This notice is consistent with physical profiles presented in Figure 3. The effects of suction/injection parameter on the skin friction and wall couple stress coefficients are also shown in Table 3. It is noted that the skin friction and wall couple stress coefficients are less with injection case, whereas they are more with suction case when compared to the case of impermeable surface. Finally, the detailed behavior of the material parameter n is given in Table 3. The skin friction decreases and wall couple stress increases with increase of material parameter n.Table 3
Effects of skin friction and wall couple stress for varying values of Biot numbers Bi, micropolar parameterN, material parameter n, and suction/injection parameter fw.
N
Bi
f
w
n
C
f
G
r
x
¯
1
/
4
M
w
G
r
x
¯
1
/
2
0.1
0.1
0.5
0.5
2.424227
0.63279
0.3
0.1
0.5
0.5
2.502716
0.639966
0.5
0.1
0.5
0.5
2.647987
0.630863
0.7
0.1
0.5
0.5
2.943561
0.603879
0.9
0.1
0.5
0.5
3.874631
0.550891
0.5
0.1
0.5
0.5
2.647987
0.630863
0.5
1.0
0.5
0.5
3.211755
0.805879
0.5
5.0
0.5
0.5
3.533141
0.908078
0.5
20.0
0.5
0.5
3.629791
0.939137
0.5
0.1
−0.5
0.5
2.426265
0.336219
0.5
0.1
0.0
0.5
2.530014
0.465429
0.5
0.1
1.0
0.5
2.737553
0.820899
0.5
0.1
2.0
0.5
2.721244
1.193596
0.5
0.1
0.5
0.0
2.916655
−0.289339
0.5
0.1
0.5
0.5
2.647987
0.630863
0.5
0.1
0.5
1.0
2.232556
2.066191


## 8. Conclusions

In this composition, the similarity solution of the free convection flow on a permeable vertical plate of a micropolar fluid under the convective boundary condition is obtained using Lie group transformations. Using the similarity variables, the governing equations are transformed into a set of nondimensional parabolic equations. These equations are solved numerically using spectral quasi-linearisation method. The numerical computation is carried out for various values of nondimensional physical parameters. The main findings are summarized as follows:(i)
The numerical results indicate that velocity distribution is less near the plate but it is more far away from the plate; the wall couple stress coefficient and rate of heat and mass transfers are lower but the temperature and concentration distributions and the skin friction coefficient are higher for the micropolar fluids in comparison with those of viscous fluids. Also, the reverse rotation of microrotation near two boundaries is found with the increasing value ofN.(ii)
An increase in Biot numberBi decreases concentration distribution, whereas it causes an increase in temperature distribution, skin friction and wall couple stress coefficients, and heat and mass transfer rates within the boundary layer. Further, enhancing in Biot number Bi enhances velocity distribution near the plate but shows the reverse behavior far away from the plate. We observe reverse rotation of microrotation near two boundaries within the boundary layer in the presence of Biot number.(iii)
Less velocity, temperature, and concentration distributions are observed, more skin friction and wall couple stress coefficients and heat and mass transfer rates in the case of suction compared to the case of injection. Further, microrotation decreases near the wall and depicts the opposite trend far away from the wall.(iv)
It is found that microrotation, temperature, and concentration distributions and skin friction coefficients are more in the case of a micropolar fluid with strong concentration (i.e.,n=0) when compared to the case of a micropolar fluid with weak concentration (i.e., n=1/2). Further, velocity is less in the case of n=0 when compared to the case of n=1/2 near the wall and shows the opposite trend far away from the wall.



---

*Source: 650813-2015-05-27.xml*