# MHD Nonlinear Mixed Convection Flow of Micropolar Nanofluid over Nonisothermal Sphere

Authors: Wubshet Ibrahim; Chaluma Zemedu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

In this paper, two-dimensional steady laminar boundary layer flow of a nonlinear mixed convection flow of micropolar nanofluid with Soret and magnetic field effect over a nonisothermal sphere is evaluated. The mathematical formulation for the flow problem has been made with appropriate similarity transformation and dimensionless variables, and the main nonlinear boundary value problems were reduced into mixed high-order nonlinear ordinary differential equations. Solution for velocity, microrotation, temperature, and concentration has been obtained numerically. The equations were calculated using method bvp4c from Matlab software for various quantities of main parameters. The effects of various parameters on skin friction coefficientf″0, wall duo stress coefficient -G′0, and convection mass transfer coefficient -Φ′0 are analysed and presented through the graphs and tables. The convergence test has been maintained. For the number of points greater than the suitable mesh number of points, the precision is not influenced but the set time is increased. Moreover, a comparison with a previous paper, obtainable in the literature, has been presented and an excellent agreement is obtained. The findings indicate that an increase in the values of nonisothermal parameters (m, P), magnetic Ma, thermal and solutal nonlinear convection (λ, s) parameter, and Soret number is to enhance the temperature difference between the boundary layer and ambient fluid to diffuse which increases the velocity profile f′ζ and their boundary layer thicknesses near the surface of the sphere.


---

## Body

## 1. Introduction

The investigation of MHD flow, heat and mass transfer has drawn attention of several researchers as a result of its request to numerous technological and different engineering and industrial developments, for instance, magnetic materials deal out, distillation of crude oil, magneto hydrodynamic electrical power creation, glass production, paper manufacturing, heat exchange devices, etc. Prasad et al. [1] have presented changeable thickness on mixed convection and heat transfer past an elastic vertical sheet with viscous electrically conducting fluid flow. However, the massive recirculation in reverse flow structure on combined convection past a heated sphere at small Reynolds number has been discussed by Kotouc et al. [2].Subba et al. [3], Chamkha and Al-Mudhaf [4], and Balla and Naikoti [5] have reported the effects of convection slip on non-Newtonian nanofluid, magnetic field, and thermal radiation on a permeable sphere at uniform heat and mass fluxes with time dependent, correspondingly. The influences of changeable fluid density on the combined convection heat transfer have been analysed by Liping et al. [6]. Results show that the effect of changeable fluid density on velocity and temperature under gravity contains two branches: the free convection and forced convection parts. Gaffar et al. [7] have studied the behaviors of heat and mass on non-Newtonian Powell–Eyring fluid flow past an isothermal sphere. Findings indicate that temperature and concentration distributions exhibit increasing performance for enhance of the values of dimensionless stream-wise coordinate. Furthermore, the influence of changeable viscosity on variation of time in free convection MHD flow toward an isothermal sphere in a rotating system has been demonstrated by Dogbalous et al. [8].Alzgool et al. [9] and Salleh and Nazar [10] have examined the mixed convection boundary layer flow starting from the solid sphere including thermal radiation effect and Newtonian heating in MHD flow micropolar fluid. Wubshet and Shanker [11] have reported the influences of magnetic field and thermal radiation on nanofluid past nonisothermal sphere along stretching sheet including nonisothermal surface temperature as well as nonisothermal heat flux. Recently, Wubshet and Chaluma [12] have investigated the impact of multiple slip conditions and thermal and solutal nonlinear convection on micropolar nanofluid flow over an isothermal sphere.All the abovementioned studies have been on the flow over a plane, isothermal sphere, and in many other fields. Here, we examine numerically the nonlinear mixed convection flow of micropolar nanofluid over a nonisothermal sphere in the presence of magnetic and nonisothermal indexes and dimensionless stream-wise coordinate impacts, using bvp4c from Matlab. The impacts of physical parameters on fluid velocity, microrotation, temperature, and concentration are discussed and shown in graphs and tables


## 2. Mathematical Formulation

Consider a two-dimensional steady laminar boundary layer flow of a viscous micropolar nanofluid over a nonisothermal spherical body of unchanging radiusc. Assume a nonisothermal surface temperature Tw and the micropolar nanofluid is a constant concentration Cw. The temperature and concentration at far, correspondingly, are T∞ and C∞. The coordinates x and y are chosen such that x measures the distance toward the surface of the sphere from the lower point and y measures the distance normal to the surface of the sphere, as shown in Figure 1.Figure 1
Flow diagram and coordinate system.Following Amanulla et al. [13], the differential equations governing the present problem are given as follows:(1)∂ur∂x+∂vr∂y=0,(2)u∂u∂x+v∂u∂y=ν+κρ∂2u∂y2+κρ∂H∂y+gR0T−T∞+RtT−T∞2+gR0C−C∞+RsC−C∞2sinxc+σB02ρu,(3)u∂H∂x+v∂H∂y=φρj∂2H∂y2−κρj2H+∂u∂y,(4)u∂T∂x+v∂T∂y=κρcp∂2T∂y2+γDB∂C∂T∂y∂y+γDT1T∞∂T∂y2,(5)u∂C∂x+v∂C∂y=DB∂2C∂y2+DctTkTm∂2T∂y2+DTT∞∂2T∂y2.Subject to the proper velocity, microrotation, and concentration boundary conditions,(6)u=a∂u∂y,u=a∂u∂x,v=0,H=−δ∂u∂y,H=−δ∂u∂x,C=Cw+e∂C∂y,C=Cw+e∂C∂x,aty=0,u⟶u∞,H⟶0,C⟶C∞,asy⟶∞.The proper thermal boundary conditions are(1)
The fixed surface temperature (FST):(7)T=Tw=T∞+xmc,∂T∂x=0,aty=0,T⟶T∞,asy⟶∞.(2)
The fixed heat flux (FHF):(8)−K∂T∂y=qw=xPc,−K∂T∂x=xPc,aty=0,T⟶T∞,asy⟶∞,where u and v are the components of velocity in the x and y axes, correspondingly, a stands for momentum slip factor, e is solutal jump factor, m and P are integer constant, H is the component of microrotation perpendicular to the xy plane. Here, δ is a constant which lies between 0 and 1 stands for different categories of variations for concentration. When δ=0, that is, the movement particle concentration is concentrated, then the microelements near to the wall surface are incapable to rotate; when δ=1/2, the antisymmetric part of the stress tensor disappears, which stands for puny concentration. The case δ=1 is meant for modeling the turbulent boundary layer flows. μ stands for the coefficient of fluid viscosity, ρ is the density, ν=μ/ρ is kinematic viscosity, cp stands for the specific heat, j=d2/Gr is the microinertia per unit mass, φ=μ+κ/2j is spin gradient viscosity, κ is the vortex viscosity, T stands for temperature, K is the thermal conductivity of the fluid, Tw=T∞+ΔT and Cw=C∞+ΔC represent the variables temperature and concentration at surface, where ΔT and ΔC are constants which give the rate of growth of temperature and concentration alongside the surface and T∞ and C∞ stand for the uniform temperature, concentration of the free stream, and g represents gravity. ∇ρ/ρ=R0T−T∞+RtT−T∞2+R0∗C−C∞+RsC−C∞2, where R0 and R0∗ are constants and Rt and Rs are the constant coefficients of thermal and volumetric expansion, respectively. This relation will be nonlinear density temperature and concentration (NDTC) variation. γ=ρcp/ρcf is the ratio between the effective heat capacity of the nanoparticle material and the heat capacity of the fluid, Dct represents the Soret-type diffusivity, Tm is fluid mean temperature, Tk is the thermal diffusion ratio, DT stands for the thermophoretic diffusion coefficient, σ is electrical conductivity, and B0 represents magnetic field. Let rx=csinx/c be the radial distance from the symmetrical axis to the surface of the sphere.Following Amanulla et al. [13], the non dimensional variables are as follows:(9)ε=xc,ζ=Gr1/4cy,r=r∗c,f=ψνεGr1/4,G=c2νεGr3/4H,Φ=C−C∞Cw−C∞.For thermal boundary conditions, we consider the nondimensional quantities as follows:(a)
For fixed surface temperature,m is the wall temperature parameter, Tw is the temperature at the wall, and c is the constant. The dimensionless temperature in FST can be selected as(10)ϑ=T−T∞Tw−T∞.(b)
For fixed heat flux FHF, the wall heat flux is taken to be changeable with radial distancer(x) from the symmetrical axis to the surface of the sphere, and we introduce the dimensionless temperature as(11)θ=T−T∞Tw−T∞,whereTw−T∞=xP/cc/Gr1/4.The continuity equation can be integrated by introducing a stream functionψx,y as ur=∂rψ/∂y and vr=−∂rψ/∂x, which satisfies the continuity equation (1). Hence, equations (2)–(5) can be written in the following nonlinear system of PDEs:(a)
For FST case,(12)1+αf‴+1+εcotεf+ε∂f∂εf″−ε∂f′∂εf′−f′2+Maf′+αG′+ϑ+λϑ2sinεε+GmGrΦ+sΦ2sinεε=0,(13)1+α2G″+1+εcotεfG′−Gf′−α2G+f″+ε∂f∂εG′−∂G∂εf′=0,(14)1Prϑ″+1+εcotεfϑ′+ε∂f∂εϑ′−∂ϑ∂εf′+Ntϑ′2−mϑf′+NbΦ′ϑ′=0,(15)1PrLeΦ″+Soϑ″+NtNbPrLeϑ″+1+εcotεf+ε∂f∂εΦ′−ε∂Φ∂εf′=0.Subject to the boundary conditions,(16)atζ=0:f′0=bvf″0,∂f0∂ε=0,ϑ0=1,∂ϑ0∂ε=0,G0=−δf″0,∂G0∂ε=bhf″0,Φ0=1+bmΦ′0,∂Φ0∂ε=Gr1/4Φ′0,asζ⟶∞:f′=0,ϑ=0,G=0,Φ=0.(b)
For FHF case,(17)1Prθ″+1+εcotεfθ′+ε∂f∂εθ′−∂θ∂εf′+Ntθ′2−Pθf′+NbΦ′θ′=0.Subject to the boundary conditions,(18)atζ=0:f′0=bvf″0,∂f0∂ε=0,θ′0=−1,∂θ0∂ε=−1−Pθε,G0=−δf″0,∂G0∂ε=bhf″0,Φ0=1+bmΦ′0,∂Φ0∂ε=Gr1/4Φ′0,asζ⟶∞:f′=0,θ=0,G=0,Φ=0,where prime represent differential with respect toζ, ε=x/c stands for dimensionless stream-wise coordinate, bv=aGr1/4/c is velocity slip parameter, bh=δc+δ2/ε denotes microrotation slip parameter, bm=eGr1/4/c stands for solutal jump parameter, α=κ/μ the is material parameter, Gr=c3gR0Tw−T∞/ν2 represents the ratio of the buoyancy force which arise from temperature difference to the viscous force times inertia force to viscous force termed as thermal Grashof number, Gm=c3gR0∗Cw−C∞/ν2 represents the ratio of the buoyancy force which arise from concentration difference to the viscous force times inertia force to viscous force termed as solutal Grashof number, and λ=Rt/R0Tw−T∞ and s=Rs/R0∗Cw−C∞ are thermal and solutal nonlinear convection parameters, respectively. We observe that, for λ=s=0, the flow of equation (12) becomes a convective micropolar nanofluid. Also, the relation of buoyancy force (Gm : Gr) can be seen in the equation is the dimensionless parameter indicating the ratio between the buoyancy force caused by concentration change and the buoyancy force owing to temperature change. Furthermore, Gm=0 is for nonbuoyancy influence because of mass diffusion, Gr⟶∞ i.e., Gm/Gr≈0 is for nonbuoyancy effect owing to thermal diffusion and Gm/Gr=1 is for thermal and mass buoyancy forces of the same strength. So=DctTkTw−T∞/TmνCw−C∞ stands for Soret number, Pr=ρνcp/K is Prandtl number, Le=K/DBρcp is Lewis number, Nb=ρcp/ρcfDB/νCw−C∞ stands for Brownian motion parameter, Nt=ρcp/ρcfDT/T∞νTw−T∞ represents thermophoresis parameter, Ma=σ/ρB02c2/νGr1/2 is magnetic parameter, and m stands for power index of terms in the nonisothermal surface.It can be appeared that, in the lower point of the sphereε≈0, equations (12)–(18) reduced to the next nonlinear system of ordinary differential equations.(a)
For FST case,(19)1+αf′″+ff″−f′2+Maf′+αG′+ϑ+λϑ2+GmGrΦ+sΦ2=0,(20)1+α2G″+fG′−Gf′−α2G+f″=0,(21)1Prϑ″+fϑ′+Ntϑ′2−mϑf′+NbΦ′ϑ′=0,(22)1PrLeΦ″+Soϑ″+NtNbPrLeϑ″+fΦ′=0.(b)
For FHF case,(23)1Prθ″+fθ′+Ntθ′2−Pθf′+NbΦ′θ′=0.Subject to the boundary conditions,(24)atζ=0:f′0=bvf″0,ϑ0=1,G0=−δf″0,Φ0=1+bmΦ′0,asζ⟶∞:f′=0,ϑ=0,G=0,Φ=0.The physical quantities of awareness in this problem are the limited skin friction coefficientCf, wall couple stress mw, the Nusselt number Nu, and the Sherwood Sh number, and they can be written as follows:(25)Cf=τwρuw,Nu=qwκTw−T∞,Sh=qmDBCw−C∞,mw=−ν+κ2ρj∂H∂yy=0,whereτw=μ+κ∂u/∂yy=0+κH0, qw=−K∂T/∂yy=0, and qm=−DB∂C/∂yy=0.Using the nondimensional variables (9) and (10) and boundary condition (18), the limited skin friction coefficient Cf, wall couple stress mw, the Nusselt number Nu, and the Sherwood number Sh are(26)f″0=CfaεGr3/4c21+α1−δ,−G′0=mwcν2εa1+α/2Gr,NucGr1/4=−ϑ′0,ShcGm1/4=−Φ′0.


## 3. Numerical Solution

The five-pair high-order ordinary differential equations (19)–(23), subjected to the boundary condition equation (18), are examined numerically using the function bvp4c from Matlab software for different values of physical parameters, that is, slip parameters (bv, bh, bm), Prandtl number Pr, thermal Grashof number Gr, Solutal Grashof number Gm, Brownian motion parameter Nb, thermophoresis parameter Nt, Lewis number Le, thermal nonlinear convection parameter λ, solutal nonlinear convection parameter s, material parameter α, microrotation parameter δ, dimensionless stream-wise coordinate ε, Magnetic parameter Ma, and nonisothermal parameters (m, P).Mathematical results are found using Matlab BVP solver bvp4c from Matlab which is a finite difference code that make the three-stage Lobatto IIIa formulation real. To apply bvp4c from Matlab, first, equations (12)–(18) are changed into a system of first-order equations.Second, establish the boundary value problem and use bvp solver bvp4c from Matlab. In solving the BVP by means of Matlab, bvp4c has only three points of views: a function ODEs for calculation of the residual in the boundary conditions, and a building solint that gives a guess for a mesh. The ODEs are handled exactly as in Matlab IVP solvers. Further clarification on the procedure of bvp4c is found in the book by Shampine et al. [14].Lety1=f, y2=f′, y3=f″, y4=h, y5=h′, y6=ϑ, y7=ϑ′, y8=Φ, y9=Φ′, y10=∂f/∂ε, y11=∂h/∂ε, y12=∂ϑ/∂ε, y13=∂Φ/∂ϵ, and y14=∂f′/∂ϵ, and y=f,f′, f″, h, h′, ϑ, ϑ′, Φ, Φ′, ∂f/∂ε, ∂h/∂ε, ∂ϑ/∂ε, ∂Φ/∂ε, ∂f′/∂εT gives(27)ddηy1y2y3y4y5y6y7y8y9y10y11y12y13y14=y2y311+α∗y2∗y2−y1+εcotε∗y1∗y3+ε∗y14∗y2−y10∗y3−α∗y5−Ma∗y2−y6+λ∗y6∗y6+GmGr∗y8+s∗y8∗y8∗sinεεy511+0.5∗α∗y4∗y2+α∗2∗y4+y3−y1∗y5−εcotε∗y1∗y5−α∗y10∗y5−y11∗y2y7−Pr∗y1∗y7+εcotε∗y1∗y7+ε∗y10∗y7−y12∗y2+Nt∗y7∗y7−m∗y6∗y2+Nb∗y9∗y7y9Le∗Pr∗So+NtNb∗Pr∗y1∗y7+εcotε∗y1∗y7+ε∗y10∗y7−y12∗y2+Pr∗Nt2Nb∗y7∗y7+Pr∗Nt∗y9∗y7−Pr∗m∗y6∗y2−Le∗Pr∗y1∗y9+εcotε∗y1∗y9+y10∗y9+y13∗y2y10y11y12y13y14.


## 4. Results and Discussion

In this part, the effects of different governing physical parameters on dimensionless velocity, temperature, concentration, microrotation, skin friction, wall couple stress coefficients, and mass transfer coefficient have been discussed.


### 4.1. Velocity and Microrotation Profiles

The dimensionless velocity profile graphs off′ζ, skin friction f″0, wall couple stresses -G′(0), and angular velocity profile Gζ for various magnitude of the microrotation slip parameter bh, velocity slip parameter bv, solutal and thermal nonlinear convection (s, λ), material parameters α, dimensionless stream-wise coordinate ε, thermal and solutal Grashof numbers (Gr, Gm), nonisothermal parameters (m, P), and magnetic parameter Ma are represented in Figures 2–15. Figures 2 and 6, indicate that enhancing of α, n, and bh reduces the density of the fluid which raise the opposition to rotate the fluid that cause increase the wall couple stresses -G′(0) and their boundary layer thickness.Figure 2
−G′(0) profile for different values ofα when Gr = 1.5, Gm = 10, Nt = 0.7, m = 1, ε = 20, α = s = Ma = 0.1, bm = bv = Le = n = λ = Nb=So = 0.5, and Pr = 0.72.Figure 3
f″(0) graph for different values of Gr when Gm = 6.5 and Nt = m = 1.Figure 4
f″ profile for different values of Gm when Gr = 2 and Nt = m = 1.Figure 5
f″(0) graph for different values of λ when Gr = 2, Gm = 5, and Nt = m = 1.Figure 6
G′(0) profile for different values of n when Gr = 1.5, m = 1, Gm = 9, and Nt = 0.7.Figure 7
f″ profile for different values of s when Gr = 2, Gm = 5, and Nt = m = 1.Figure 8
f′ profile for different values of m when Gr = bv = 0.1, Gm = 10, s = 0.24, Ma = 0.11, Le = 0.54, and Nt = 0.1.Figure 9
f′ζ graph for different values of Ma when ε=0.2,Gr=0.1,m=1,s=0.27,Nt=0.4,Pr=α=0.72−b±b2−4ac/2a.Figure 10
f′ profile for different values of λ when P=1 and Nt = Gr = Gm = 0.1.Figure 11
f′ζ graph for different values of s when PP=Gr=Gm=1.Figure 12
f′ζ profile for different values of So when P=Gr=Gm=bv=1 and Nt = 0.1.Figure 13
f′ζ graph for different values of P when Gr=1.2,Gm=2.1,bv=0.5,s=0.72,Nt=0.5,α=5,Nb=bm=bh=0.1,andλ=0.7.Figure 14
Gζ profile for different values of ε when Gm=20,m=1,s=0.1,Nb=1,λ=0.5,Gr=18,andNt=0.1.Figure 15
f′ζ graph for different values of ε when Gr = 13, Gm = m = 1, bv = 0.35, Nb = 0.1, and Nt = 0.28.Figure3 shows that improving of Gr and bv decreases the kinematic viscosity of the fluid which raise the opposition to flow the fluid that cause decrease the skin friction coefficients f″0 and their boundary layer thickness. Figures 4, 5, and 7 reveal that an increase in Gm, s, λ, and bv improves thermal and solutal volumetric expansion coefficients which improves the particle movement, and as a result, the skin friction coefficent f″0 is enhanced. Also, Figures 8–13 illustrate that increase of m, Ma, λ, s, So, and P boost the temperature difference between the boundary layer and ambient fluid to diffuse, which increases the velocity profile f′ζ and their boundary layer thicknesses near the surface of the sphere and some distance far away from the wall (Figure 9). Moreover, when the values of m = 1 and m = 120, the graph of velocity profile f′ζ is concave downwards near the points (0.3648, 0.1273) and (0.2679, 0.1874), respectively, as shown in Figure 8; as well as, when the values of P=7 and P=35, this profile is concave downwards near the points (0.9555, 0.03508) and (1.008, 0.03008), accordingly, as shown in Figure 13; while, for the value of λ, the velocity profile f′ζ and its boundary layer thickness is decreasing in the interval of ζ [0.84, 3], as shown in Figure 10. Furthermore, Figures 14 and 15 show that the increase of ϵ brings down the microrotation profile and the velocity profile near the wall of the sphere.


### 4.2. Temperature and Concentration Profiles

The importance of solutal jump parameter bm, Lewis and Prandtl numbers (Le and Pr), Brownian motion Nb, thermophoresis Nt, thermal nonlinear convectionλ, solutal nonlinear convection parameters s, dimensionless stream-wise coordinate ε, Soret number So, nonisothermal parameters (m, P), magnetic parameter Ma on temperature, concentration, and Sherwood numbers graphs are shown in Figures 16–33. Figures 16–21 indicate that increase of magnetic Ma, thermal nonlinear convection λ, and solutal nonlinear convection parameters s, enhances the density of the fluid which result in decline in both temperature distributions θ,ϑ and their thermal boundary layer thickness. It is also observed that the temperature in both cases for FST and FHF reduces when the values of magnetic Ma, thermal nonlinear convection λ, and solutal nonlinear convection parameters s improve at fixed values of ζ. Similarly, the concentration distribution exhibits the reduction behavior for enlarge the values of solutal nonlinear convection parameters s, as shown in Figure 26. Moreover, Figure 24 demonstrates that boost of Soret number So enhances density of the fluid which associate through decrease in temperature and thermal boundary layer thicknesses; however, Figure 25 shows that enlarging of Lewis number Le enhances thermal conductivity and concentration difference between the boundary layer and ambient fluid to diffuse the temperature which is associated with the increase in temperature distribution as well as thermal boundary layer thickness.Figure 16
ϑζ profile for different values of Ma when Gr = 1.25, Gm = s = 0.1, α = 0.3, λ = bv = bm = bh = Le = Nb = So = n = 0.5, m = 1, Nt = 3, Pr = 0.77, and ε=20, for FST.Figure 17
θζ profile for different values of Ma when Gm = Gr = bm = s = bh = α = λ = Nb = 0.1, bv = Le = So = n = 0.5, P=1, Nt = 0.7, and ε=20, for FHF.Figure 18
θζ graph for different values of λ when Gm = Gr = bm = bh = Nt = 0.1 and bv=P=1, for FHF.Figure 19
ϑζ graph for different values of λ when Gr = 1.25, Gm = 0.1, α = 0.3, and Pr = 0.77, for FST.Figure 20
ϑζ profile for different values of s when Gm = Gr = 0.1, Nb = 1, α=0, and Nt = 3, for FST.Figure 21
θζ graph for different values of s when Gm = 0.72, Gr = bm = bh = Nt = 0.1, and bv=P=1, for FHF.Figure 22
ϑζ profile for different values of ε when Gr = 20, Gm = Nb = 0.1, λ = bv = bm = bh = 0.5, and Nt = 3, for FST.Figure 23
θζ graph for different values of ϵ when Gr = bm = bh = Nt = 0.1, λ = bv = Nb = Gm = 0.5, and P=1, for FHF.Figure 24
ϑζ profile for different values of So when Gm = Gr = α = 0.1, λ = bv = bm = bh = 0.5, m = 1, and Nt = 0.4, for FST.Figure 25
ϑ0 graph for different values of Le when Gm = m = NB = 1, Gr = α = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.4, for FST.Figure 26
Φζ graph for different values of s when Gm = 0.1, Gr = 1, and Nt = 0.4.Figure 27
Φζ profile for different values of m when Gm = 1, Gr = 0.2, Pr = Nb = 2, = bh = 0.5, and Nt = 0.1.Figure 28
−Φ′0 profile for different values of Pr when Gr = bh = 0.1, λ = bv = bm = bh = 0.5, Gm = m = 1, and Nt = 0.7.Figure 29
−Φ′0 graph for different values of So when Gm = m = 1, Gr = bh = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.2.Figure 30
−Φ′0 profile for different values of Nb when Gm = 1, Gr = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.2.Figure 31
−Φ′ζ profile for different values of Le when Gm = Gr = m = 0.1 s = α = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.4.Figure 32
Grid-independence test forf′ζ when Gr=13,m=1, λ=bv=bm=Nt=n=So=Le=05, Gm=bh=bv=s=Nb=Ma=α=0.1,n=09,Nb=3,andPr=0.72.Figure 33
Grid-independence test forGζ when Gr=10,Gm=m=1,λ=bv=bm=So=Le=05, bh=s=NtMa=α=0.1,n=09,Nb=3,andPr=0.72.Figures22 and 23 show that improving of ϵ increases the thermal and solutal boundary layer of fluid which is associated with the enhance in both temperature distributions cases, i.e., in FST ϑζ and FHF θζ. Similarly, Figure 27 reveals that boost of power index m enlarge the solutal boundary layer thickness in the case of FST near the surface of the sphere.Figures28–31 demonstrate that upsurge of Prandtl number, Soret number, Lewis number, and Brownian motion parameter causes increase in the density of the fluid that results rise of the local Sherwood number. Figures 32–35 illustrate the grid-independence test. From these figures, it is observed that, for the number of points greater than suitable mesh number of points, the precision is not influenced but the set time is raised.Figure 34
Grid-independence test forϑζ when Gr=10,Gm=m=1,λ=bv=bm=So=Le=05, bh=s=NtMa=α=0.1,n=09,Nb=3,andPr=0.72.Figure 35
Grid-independence test forΦζ when Gr=10,Gm=m=1,λ=bv=bm=So=Le=05, bh=s=NtMa=α=0.1,n=09,Nb=3,andPr=0.72.The presentation of grid-independence test is shown in Tables1 and 2, to maintain the point of exactness known as grid convergence test. It started by the common mesh with 30 numbers of points. By enhancing number of points twice and triple, we get the medium mesh, 90 points of accuracy, and suitable mesh, 180 points of accuracy for velocity, angular, temperature, and concentration gradient values. For the number of points greater than suitable mesh number of points, the precision is not influenced but the set time is raised.Table 1
Grid-independence test for velocity gradientf″ζandangularvelocitygradientG′ζ when n=Nt=bm=So=Le=0.5, Pr=0.7,Nb=bv=Gm=bh=s=0.1,andGr=13.
εCommon mesh, 30Medium mesh, 90Suitable mesh, 180f″ζG′ζf″ζG′ζf″ζG′ζ00.03190.02760.03190.0270.03190.02710.03000.02500.03000.02500.03000.025020.02870.02340.02780.02210.02780.022130.02780.02210.02780.02210.02780.022140.02690.02100.02690.02100.02690.021050.02620.02010.02620.02010.02620.0201Table 2
Grid-independence test for Nusselt number–ϑ′ζandSherwoodnumber−Φ′ζ when n=0.9,Gm=m=1,bv=bm=So=Le=0.5, Pr=0.7,Nt=bh=s=α=Ma=0.1,Nb=3,andGr=10.
εCommon mesh, 30Medium mesh, 90Suitable mesh, 180−ϑ′ζ−Φ′ζ−ϑ′ζ−Φ′ζ−ϑ′ζ−G∗Φ′ζ00.36970.32370.36970.32370.36970.323710.36070.29550.36070.29550.36070.295520.34970.28530.34970.28530.34970.285330.34010.27770.34010.27770.34010.277740.33150.27170.33150.27170.33150.271750.32370.26650.32370.26650.32370.2665Tables3 and 4 are drawn to compare the exactness of the method used in association with earlier reported data feasible in the literature studies which have been carried out. From Tables 1 and 2, it can be seen that the numerical magnitudes of the Nusselt number −ϑ′0, in the present study, for various values of χ and ϵ, n = 0.5, Pr = 0.7 and Pr = 7, are in an excellent agreement with the pervious results of published articles by Nazar et al. [15] and Wubshet and Chaluma [12]. The outcomes show that the numerical used in the information is truthful and highly accurate.Table 3
Comparison of Nusselt number -ϑ′0 when n = 0.5 and Pr = 0.7, for different values of α and ϵ, with previously published result.
Present resultNazar et al. [15]Wubshet and Chaluma [12]αε−ϑ′0−ϑ′0−ϑ′00.00.00.45760.45760.45770.50.43360.43360.43341.00.41650.41660.41661.50.40370.40350.40402.00.39310.39300.39320.0100.45610.45650.4560200.45350.45330.4534300.44810.44800.4486400.44060.44050.4406500.43090.43080.4309600.41890.41890.4188700.40480.4046800.38780.3879900.36870.3684Table 4
Comparison of Nusselt number−ϑ′0 when n = 0.5 and Pr = 7, for different values of α and ε, with previously published result.
Present resultNazar et al. [15]αε−ϑ′0−ϑ′00.00.00.95980.95950.50.89040.89051.00.84430.84431.50.80960.80962.00.78090.7805Tables5 and 6 point out that upsurge in values of Ma causes enhancement in the skin friction coefficient f″0, wall couple stress coefficient -G′0, Nusselt number in both cases −ϑ′0,−θ′0, and Sherwood Number -Φ′0. This result happens due to improved density of the fluid which boosts the viscosity of it to flow, rotate, and diffuse that associate through enhancement of these physical quantities. However, as the magnitudes of m and P increase, reduce the viscosity of the fluid flow which lowers opposition to flow and rotate the fluid that results in reduction of skin friction and wall couple stress coefficients. By similarl reasoning, Sherwood Number -Φ′0 of the fluid exhibits decreasing behavior as values of P increase, but Nusselt number-ϑ′0 upsurge with m. However, increase in values of ϵ improves skin friction coefficient, Sherwood Number -Φ′0, and Nusselt number-ϑ′0.Table 5
The computed values of skin friction coefficientf″(0), wall couple stress -G′0, Nusselt number -ϑ′0, and Sherwood number -Φ′0 when Gr = Nb = Le = λ = 0.1, bv = bh = s = n = So = 0.5, Nt = 0.4, α = 1, Gm = 5, and Pr = 0.72, for different values of Ma, P, and ϵ.
Mamϵf″0−G′0−ϑ′0−Φ′011102.28122.07113.0357−5.56621.51102.31442.11593.0536−5.592821102.34802.16143.0712−5.619310.1102.2911-2.08413.010910.5102.2867-2.07833.022111102.2812-2.07113.035911.5102.2757-2.06403.04971112.42273.0979−5.712611102.28123.0359−5.589611202.14662.8395−5.436011301.97952.6535−5.235511401.80712.5937−4.9782Table 6
The computed values of skin friction coefficientf″0, wall couple stress -G′0, and Sherwood number -Φ′0 when Gr = Nb = Le = λ=0.1,bv=bh=s=n=So=0.5,ε=10Nt=0.4, α=1,Gm=5,andPr=0.72, for different values of Ma and P.
MaPf″0−G′0−Φ′0110.8882−0.4811−0.67931.510.9226−0.5008−0.6812210.9464−0.5211−0.682910.10.9598−0.5366−0.915910.50.9342−0.5129−0.8157110.8992−0.4811−0.679311.50.8601−0.4465−0.5282


## 4.1. Velocity and Microrotation Profiles

The dimensionless velocity profile graphs off′ζ, skin friction f″0, wall couple stresses -G′(0), and angular velocity profile Gζ for various magnitude of the microrotation slip parameter bh, velocity slip parameter bv, solutal and thermal nonlinear convection (s, λ), material parameters α, dimensionless stream-wise coordinate ε, thermal and solutal Grashof numbers (Gr, Gm), nonisothermal parameters (m, P), and magnetic parameter Ma are represented in Figures 2–15. Figures 2 and 6, indicate that enhancing of α, n, and bh reduces the density of the fluid which raise the opposition to rotate the fluid that cause increase the wall couple stresses -G′(0) and their boundary layer thickness.Figure 2
−G′(0) profile for different values ofα when Gr = 1.5, Gm = 10, Nt = 0.7, m = 1, ε = 20, α = s = Ma = 0.1, bm = bv = Le = n = λ = Nb=So = 0.5, and Pr = 0.72.Figure 3
f″(0) graph for different values of Gr when Gm = 6.5 and Nt = m = 1.Figure 4
f″ profile for different values of Gm when Gr = 2 and Nt = m = 1.Figure 5
f″(0) graph for different values of λ when Gr = 2, Gm = 5, and Nt = m = 1.Figure 6
G′(0) profile for different values of n when Gr = 1.5, m = 1, Gm = 9, and Nt = 0.7.Figure 7
f″ profile for different values of s when Gr = 2, Gm = 5, and Nt = m = 1.Figure 8
f′ profile for different values of m when Gr = bv = 0.1, Gm = 10, s = 0.24, Ma = 0.11, Le = 0.54, and Nt = 0.1.Figure 9
f′ζ graph for different values of Ma when ε=0.2,Gr=0.1,m=1,s=0.27,Nt=0.4,Pr=α=0.72−b±b2−4ac/2a.Figure 10
f′ profile for different values of λ when P=1 and Nt = Gr = Gm = 0.1.Figure 11
f′ζ graph for different values of s when PP=Gr=Gm=1.Figure 12
f′ζ profile for different values of So when P=Gr=Gm=bv=1 and Nt = 0.1.Figure 13
f′ζ graph for different values of P when Gr=1.2,Gm=2.1,bv=0.5,s=0.72,Nt=0.5,α=5,Nb=bm=bh=0.1,andλ=0.7.Figure 14
Gζ profile for different values of ε when Gm=20,m=1,s=0.1,Nb=1,λ=0.5,Gr=18,andNt=0.1.Figure 15
f′ζ graph for different values of ε when Gr = 13, Gm = m = 1, bv = 0.35, Nb = 0.1, and Nt = 0.28.Figure3 shows that improving of Gr and bv decreases the kinematic viscosity of the fluid which raise the opposition to flow the fluid that cause decrease the skin friction coefficients f″0 and their boundary layer thickness. Figures 4, 5, and 7 reveal that an increase in Gm, s, λ, and bv improves thermal and solutal volumetric expansion coefficients which improves the particle movement, and as a result, the skin friction coefficent f″0 is enhanced. Also, Figures 8–13 illustrate that increase of m, Ma, λ, s, So, and P boost the temperature difference between the boundary layer and ambient fluid to diffuse, which increases the velocity profile f′ζ and their boundary layer thicknesses near the surface of the sphere and some distance far away from the wall (Figure 9). Moreover, when the values of m = 1 and m = 120, the graph of velocity profile f′ζ is concave downwards near the points (0.3648, 0.1273) and (0.2679, 0.1874), respectively, as shown in Figure 8; as well as, when the values of P=7 and P=35, this profile is concave downwards near the points (0.9555, 0.03508) and (1.008, 0.03008), accordingly, as shown in Figure 13; while, for the value of λ, the velocity profile f′ζ and its boundary layer thickness is decreasing in the interval of ζ [0.84, 3], as shown in Figure 10. Furthermore, Figures 14 and 15 show that the increase of ϵ brings down the microrotation profile and the velocity profile near the wall of the sphere.


## 4.2. Temperature and Concentration Profiles

The importance of solutal jump parameter bm, Lewis and Prandtl numbers (Le and Pr), Brownian motion Nb, thermophoresis Nt, thermal nonlinear convectionλ, solutal nonlinear convection parameters s, dimensionless stream-wise coordinate ε, Soret number So, nonisothermal parameters (m, P), magnetic parameter Ma on temperature, concentration, and Sherwood numbers graphs are shown in Figures 16–33. Figures 16–21 indicate that increase of magnetic Ma, thermal nonlinear convection λ, and solutal nonlinear convection parameters s, enhances the density of the fluid which result in decline in both temperature distributions θ,ϑ and their thermal boundary layer thickness. It is also observed that the temperature in both cases for FST and FHF reduces when the values of magnetic Ma, thermal nonlinear convection λ, and solutal nonlinear convection parameters s improve at fixed values of ζ. Similarly, the concentration distribution exhibits the reduction behavior for enlarge the values of solutal nonlinear convection parameters s, as shown in Figure 26. Moreover, Figure 24 demonstrates that boost of Soret number So enhances density of the fluid which associate through decrease in temperature and thermal boundary layer thicknesses; however, Figure 25 shows that enlarging of Lewis number Le enhances thermal conductivity and concentration difference between the boundary layer and ambient fluid to diffuse the temperature which is associated with the increase in temperature distribution as well as thermal boundary layer thickness.Figure 16
ϑζ profile for different values of Ma when Gr = 1.25, Gm = s = 0.1, α = 0.3, λ = bv = bm = bh = Le = Nb = So = n = 0.5, m = 1, Nt = 3, Pr = 0.77, and ε=20, for FST.Figure 17
θζ profile for different values of Ma when Gm = Gr = bm = s = bh = α = λ = Nb = 0.1, bv = Le = So = n = 0.5, P=1, Nt = 0.7, and ε=20, for FHF.Figure 18
θζ graph for different values of λ when Gm = Gr = bm = bh = Nt = 0.1 and bv=P=1, for FHF.Figure 19
ϑζ graph for different values of λ when Gr = 1.25, Gm = 0.1, α = 0.3, and Pr = 0.77, for FST.Figure 20
ϑζ profile for different values of s when Gm = Gr = 0.1, Nb = 1, α=0, and Nt = 3, for FST.Figure 21
θζ graph for different values of s when Gm = 0.72, Gr = bm = bh = Nt = 0.1, and bv=P=1, for FHF.Figure 22
ϑζ profile for different values of ε when Gr = 20, Gm = Nb = 0.1, λ = bv = bm = bh = 0.5, and Nt = 3, for FST.Figure 23
θζ graph for different values of ϵ when Gr = bm = bh = Nt = 0.1, λ = bv = Nb = Gm = 0.5, and P=1, for FHF.Figure 24
ϑζ profile for different values of So when Gm = Gr = α = 0.1, λ = bv = bm = bh = 0.5, m = 1, and Nt = 0.4, for FST.Figure 25
ϑ0 graph for different values of Le when Gm = m = NB = 1, Gr = α = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.4, for FST.Figure 26
Φζ graph for different values of s when Gm = 0.1, Gr = 1, and Nt = 0.4.Figure 27
Φζ profile for different values of m when Gm = 1, Gr = 0.2, Pr = Nb = 2, = bh = 0.5, and Nt = 0.1.Figure 28
−Φ′0 profile for different values of Pr when Gr = bh = 0.1, λ = bv = bm = bh = 0.5, Gm = m = 1, and Nt = 0.7.Figure 29
−Φ′0 graph for different values of So when Gm = m = 1, Gr = bh = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.2.Figure 30
−Φ′0 profile for different values of Nb when Gm = 1, Gr = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.2.Figure 31
−Φ′ζ profile for different values of Le when Gm = Gr = m = 0.1 s = α = 0.1, λ = bv = bm = bh = 0.5, and Nt = 0.4.Figure 32
Grid-independence test forf′ζ when Gr=13,m=1, λ=bv=bm=Nt=n=So=Le=05, Gm=bh=bv=s=Nb=Ma=α=0.1,n=09,Nb=3,andPr=0.72.Figure 33
Grid-independence test forGζ when Gr=10,Gm=m=1,λ=bv=bm=So=Le=05, bh=s=NtMa=α=0.1,n=09,Nb=3,andPr=0.72.Figures22 and 23 show that improving of ϵ increases the thermal and solutal boundary layer of fluid which is associated with the enhance in both temperature distributions cases, i.e., in FST ϑζ and FHF θζ. Similarly, Figure 27 reveals that boost of power index m enlarge the solutal boundary layer thickness in the case of FST near the surface of the sphere.Figures28–31 demonstrate that upsurge of Prandtl number, Soret number, Lewis number, and Brownian motion parameter causes increase in the density of the fluid that results rise of the local Sherwood number. Figures 32–35 illustrate the grid-independence test. From these figures, it is observed that, for the number of points greater than suitable mesh number of points, the precision is not influenced but the set time is raised.Figure 34
Grid-independence test forϑζ when Gr=10,Gm=m=1,λ=bv=bm=So=Le=05, bh=s=NtMa=α=0.1,n=09,Nb=3,andPr=0.72.Figure 35
Grid-independence test forΦζ when Gr=10,Gm=m=1,λ=bv=bm=So=Le=05, bh=s=NtMa=α=0.1,n=09,Nb=3,andPr=0.72.The presentation of grid-independence test is shown in Tables1 and 2, to maintain the point of exactness known as grid convergence test. It started by the common mesh with 30 numbers of points. By enhancing number of points twice and triple, we get the medium mesh, 90 points of accuracy, and suitable mesh, 180 points of accuracy for velocity, angular, temperature, and concentration gradient values. For the number of points greater than suitable mesh number of points, the precision is not influenced but the set time is raised.Table 1
Grid-independence test for velocity gradientf″ζandangularvelocitygradientG′ζ when n=Nt=bm=So=Le=0.5, Pr=0.7,Nb=bv=Gm=bh=s=0.1,andGr=13.
εCommon mesh, 30Medium mesh, 90Suitable mesh, 180f″ζG′ζf″ζG′ζf″ζG′ζ00.03190.02760.03190.0270.03190.02710.03000.02500.03000.02500.03000.025020.02870.02340.02780.02210.02780.022130.02780.02210.02780.02210.02780.022140.02690.02100.02690.02100.02690.021050.02620.02010.02620.02010.02620.0201Table 2
Grid-independence test for Nusselt number–ϑ′ζandSherwoodnumber−Φ′ζ when n=0.9,Gm=m=1,bv=bm=So=Le=0.5, Pr=0.7,Nt=bh=s=α=Ma=0.1,Nb=3,andGr=10.
εCommon mesh, 30Medium mesh, 90Suitable mesh, 180−ϑ′ζ−Φ′ζ−ϑ′ζ−Φ′ζ−ϑ′ζ−G∗Φ′ζ00.36970.32370.36970.32370.36970.323710.36070.29550.36070.29550.36070.295520.34970.28530.34970.28530.34970.285330.34010.27770.34010.27770.34010.277740.33150.27170.33150.27170.33150.271750.32370.26650.32370.26650.32370.2665Tables3 and 4 are drawn to compare the exactness of the method used in association with earlier reported data feasible in the literature studies which have been carried out. From Tables 1 and 2, it can be seen that the numerical magnitudes of the Nusselt number −ϑ′0, in the present study, for various values of χ and ϵ, n = 0.5, Pr = 0.7 and Pr = 7, are in an excellent agreement with the pervious results of published articles by Nazar et al. [15] and Wubshet and Chaluma [12]. The outcomes show that the numerical used in the information is truthful and highly accurate.Table 3
Comparison of Nusselt number -ϑ′0 when n = 0.5 and Pr = 0.7, for different values of α and ϵ, with previously published result.
Present resultNazar et al. [15]Wubshet and Chaluma [12]αε−ϑ′0−ϑ′0−ϑ′00.00.00.45760.45760.45770.50.43360.43360.43341.00.41650.41660.41661.50.40370.40350.40402.00.39310.39300.39320.0100.45610.45650.4560200.45350.45330.4534300.44810.44800.4486400.44060.44050.4406500.43090.43080.4309600.41890.41890.4188700.40480.4046800.38780.3879900.36870.3684Table 4
Comparison of Nusselt number−ϑ′0 when n = 0.5 and Pr = 7, for different values of α and ε, with previously published result.
Present resultNazar et al. [15]αε−ϑ′0−ϑ′00.00.00.95980.95950.50.89040.89051.00.84430.84431.50.80960.80962.00.78090.7805Tables5 and 6 point out that upsurge in values of Ma causes enhancement in the skin friction coefficient f″0, wall couple stress coefficient -G′0, Nusselt number in both cases −ϑ′0,−θ′0, and Sherwood Number -Φ′0. This result happens due to improved density of the fluid which boosts the viscosity of it to flow, rotate, and diffuse that associate through enhancement of these physical quantities. However, as the magnitudes of m and P increase, reduce the viscosity of the fluid flow which lowers opposition to flow and rotate the fluid that results in reduction of skin friction and wall couple stress coefficients. By similarl reasoning, Sherwood Number -Φ′0 of the fluid exhibits decreasing behavior as values of P increase, but Nusselt number-ϑ′0 upsurge with m. However, increase in values of ϵ improves skin friction coefficient, Sherwood Number -Φ′0, and Nusselt number-ϑ′0.Table 5
The computed values of skin friction coefficientf″(0), wall couple stress -G′0, Nusselt number -ϑ′0, and Sherwood number -Φ′0 when Gr = Nb = Le = λ = 0.1, bv = bh = s = n = So = 0.5, Nt = 0.4, α = 1, Gm = 5, and Pr = 0.72, for different values of Ma, P, and ϵ.
Mamϵf″0−G′0−ϑ′0−Φ′011102.28122.07113.0357−5.56621.51102.31442.11593.0536−5.592821102.34802.16143.0712−5.619310.1102.2911-2.08413.010910.5102.2867-2.07833.022111102.2812-2.07113.035911.5102.2757-2.06403.04971112.42273.0979−5.712611102.28123.0359−5.589611202.14662.8395−5.436011301.97952.6535−5.235511401.80712.5937−4.9782Table 6
The computed values of skin friction coefficientf″0, wall couple stress -G′0, and Sherwood number -Φ′0 when Gr = Nb = Le = λ=0.1,bv=bh=s=n=So=0.5,ε=10Nt=0.4, α=1,Gm=5,andPr=0.72, for different values of Ma and P.
MaPf″0−G′0−Φ′0110.8882−0.4811−0.67931.510.9226−0.5008−0.6812210.9464−0.5211−0.682910.10.9598−0.5366−0.915910.50.9342−0.5129−0.8157110.8992−0.4811−0.679311.50.8601−0.4465−0.5282


## 5. Conclusions

This paper takes into account the effects of Soret, thermal and solutal nonlinear convection, nonisothermal, and magnetic parameters as well as manifold slip conditions on a micropolar nanofluid past a nonisothermal sphere, and the effects were analysed. The boundary layer equations of the flow problem are reduced into a duo of high-order nonlinear ordinary differential equations with the similarity transformation. The obtained differential equations are solved numerically with bvp4c from Matlab software. Numerical results are acquired for various main parameters of the flow problem. The results of governing parameters are reported by means of figures and tables. The main findings are(1)
An increase in the values ofm, Ma, λ, s, So, and P allows to enhance the temperature difference between boundary layer and ambient fluid to diffuse, which increases the velocity profile f′ζ and their boundary layer thicknesses near the surface of the sphere.(2)
The existence of Soret number So and nonisothermal parameter of FSTm agrees to reduce momentum kinematic viscosity which is associated with decrease in temperature distributions and their boundary layer thicknesses.(3)
An enhancement in the values of Nt andϵ results in the boost of the temperature difference between the boundary layer and ambient fluid to diffuse the temperature which is associated with an increase in both temperature ϑζ,θζ and concentration distributions.(4)
The tempereture distributionθ within the boundary layer can be enhanced by raising the values of nonisothermal parameter P which increase the density of the fluid; however, an increase in the value of m reduces the surface temprature θ.(5)
An improvement in magnitude ofϵ is made available to improve the microrotation profile, but it reduces the velocity profile within the boundary. This result happens due to improved density of the fluid(6)
The upsurge in values of Ma from 1 to 2 results in the improved density of the fluid which is associated with an enhancement in the skin friction coefficientf″0 by 2.9 percent, wall couple stress coefficient -G′0 by 4.36 percent, Nusselt number −ϑ′0by1.16percent, and Sherwood Number -Φ′0 by 0.94 percent.(7)
Skin friction and wall couple stress coefficients exhibit decreasing behavior by 0.67 and 10 percent and by 0.96 and 16 percent for increasing the magnitudes ofm and P from 0.1 to 1.5. These effects are caused due to a decrease of viscosity of the fluid flow which lowers opposition to flow and to rotate the fluid in the boundary layer.



---

*Source: 3596368-2020-06-09.xml*