# Nonlinear Convection Flow of Micropolar Nanofluid due to a Rotating Disk with Multiple Slip Flow

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

---

## Abstract

In this analysis, steady, laminar, and two-dimensional boundary layer flow of nonlinear convection micropolar nanofluid due to a rotating disk is considered. The mathematical formulation for the flow problem has been made. By means of appropriate similarity transformation and dimensionless variables, the governing nonlinear boundary value problems were reduced into coupled high-order nonlinear ordinary differential equations with numerically solved. The equations were calculated using method bvp4c from matlab software for various quantities of main parameters. The influences of different parameters on skin friction coefficientsf″0 and G′0, wall duo stress coefficients H1′0, -H2′0, and -H3′0, the Nusselt number -θ′0, and Sherwood number Ω′0, as well as the velocities, temperature, and concentration are analysed and discussed through tables and plotted graphs. The findings indicate that an increase in the values of thermal and solutal nonlinear convection parameters allow to increase the value of velocities f′η and Gη near surface of the disk and reduce at far away from the disk as well as thermal and solutal Grashof numbers tolerate to increase in the value of radial velocity f′η near surface of the disk.


---

## Body

## 1. Introduction

Eringen [1] was the first researcher who presents the theory of micropolar fluid for which the classical Navier–Stokes theory is has a limitation for its full description. Consequently, Hamzeh [2] evaluated the behaviour of micropolar Casson fluid on natural convective flow past a solid sphere. It was found that material parameter diminishes the values of the narrow skin friction coefficient, but grows both values of the local Nusselt number and angular velocity profiles as its value rises. Wubshet and Chaluma [3] also reported the variation of density on micropolar nanofluid flow about an isothermal sphere. Moreover, wall shear stress and angular velocity gradient at the wall enhance with an enlargement in the material parameter as illustrated by Mandal and Mukhopadhyay [4]. Furthermore, the flow of MHD and heat transfer due to stretching rotating disk were examined by Mustafa [5] and Akhter et al. [6].Mostafa and Shimaa [7] reported the influence of magnetic field on flow and heat transfer of a micropolar fluid about a stretching surface in the presence of heat generation (absorption). Also, the effects of microrotation parameter on flow between a rotating and stationary disk was examined by Anwar and Guram [8]. Moreover, Rashidi and Freideonimehr [9] presented the effects of velocity and temperature slip on the entropy generation past rotating disk. Also, the effects of diffusion-thermo and thermodiffusion on radially stretching disk have been evaluated by Khan et al. [10].The impact of Prandtl number on radiative flow due to stretchable rotating disk with variable thickness was computed by Tasawar et al. [11]. Shamshuddin et al. [12] have evaluated the influence of this parameter on numerical study of heat transfer and viscous flow in a dual rotating extendable disk system employing a non-Fourier heat flux model. They indicated that, with an improvement in the Prandtl number, there is a strong decrease in temperature of fluid as well as radial skin friction. Moreover, the Stefan blowing effect on bioconvective flow and heat transfer of nanofluid over a rotating stretchable disk was reported by Lafiff et al. [13] and Yin et al. [14].Muhammad and Naeem [15] and Noor et al. [16] have examined the influences of velocity slip with magnetic field in micropolar nanofluid flow along rotating disk and in mixed convection lower flow of a micropolar nanofluid along vertically elongating surface accordingly. Also, the flow and heat transfer computations for nanofluid fluid flow and the impact of heat generation/consumption and thermal radiation over a rotating disk were analysed by Narfifah et al. [17], Mushtaq and Mustafa [18], and Anwer et al. [19]. Moreover, Tasawar et al. [20] had discussed the impact of thermal slip condition on MHD flow of Cu-water nanofluid due to a rotating disk. They found that thermal boundary thickness increased for lower thermal slip parameter values, but increasing the values of it reduced heat transfer from the disk to the adjacent fluid. Moreover, the numerical study of nanofluid flow and heat transfer over a rotating disk was examined by Ahmad et al. [21]. MHD mixed convection movement of a nanofluid over nonlinear enlarging sheet including variable Brownian and thermophoretic diffusion coefficient have been evaluated by Sumalatha and Shanker [22], and the results illustrated that the influence of nonlinear stretching parameter drops both the flow of the fluid as well as temperature distribution.All the above research articles have been on the flow over a plane, over a stretching surface, on MHD boundary flow, or in many other areas. In this paper, we evaluate numerically the nonlinear convection flow of micropolar nanofluid due to a rotating disk in the presence of multiple slip conditions, using bvp4c from matlab. The outcomes of physical parameters on fluid velocity, temperature, and concentration were discussed and indicated in graphs and tables as well.


## 2. Mathematical Formulation

Let us choose the cylindrical systemr,ω,z in the direction component of the flow velocity u,v,w and the angular velocity components H1,H2,H3 correspondingly. The study considers incompressible, laminar, and nonlinear convection flow of micropolar nanofluid over a circular disk at z=0. The disk rotates with uniform microrotation rM0 about the z-axis. As the result of revolving symmetry the end products in the azimuthal direction may be considered. The wall of the gyrating disk has constant temperature Tw and concentration Cw despite the fact the ambient temperature and concentration are represented by T∞ and C∞ correspondingly, as shown in Figure 1. By means of Mandal and Mukhopadhyay [4], Anwer and Guram [8], Tasawar et al. [11], Noor et al. [16], and Sajjad et al. [23] the governing differential equations of the flow are given as follows:(1)∂u∂r+∂w∂z+ur=0,(2)u∂u∂r+w∂u∂z−v2r=μ+κρ∂2u∂r2+1r∂u∂r−ur2+∂2u∂z2−κρ∂N2∂z+gσT−T∞+σtT−T∞2+gσC−C∞+σsC−C∞2sinω,(3)u∂v∂r+w∂v∂z+uvr=μ+κρ∂2v∂r2+1r∂v∂r−vr2+∂2v∂z2+κρ∂N1∂z−∂N3∂r+gσT−T∞+σtT−T∞2+gσC−C∞+σsC−C∞2sinω,(4)u∂w∂r+w∂w∂z=μ+κρ∂2w∂r2+1r∂w∂r+∂2w∂z2+κρN2r,(5)u∂N1∂r+w∂N1∂z−vN2r=a+b+cρj∂∂r∂N1∂r+N1r+∂N3∂z−cρj∂∂z∂N3∂r−∂N1∂z−κρj∂v∂z−2κρjN1,(6)u∂N2∂r+w∂N2∂z+vN1r=cρj∂∂r∂N2∂r+N2r+∂2N2∂z2+κρj∂u∂z−2κN2ρj,(7)u∂N3∂r+w∂N3∂z=a+b+cρj∂∂z∂N1∂r+N1r+∂H3∂z−cρj∂∂rr∂N1∂z−∂N3∂r+κρj2v−2κN3ρj,(8)u∂T∂r+w∂T∂z=Kρcp∂2T∂r2+1r∂T∂r+∂2T∂z2+γDB∂C∂T∂r∂r+∂C∂T∂z∂z+γDTT∞∂T∂r2+∂T∂z2,(9)u∂C∂r+w∂C∂z=DB∂2C∂r2+1r∂C∂r+∂2C∂z2+DTT∞∂2T∂r2+1r∂T∂r+∂2T∂z2,with boundary conditions(10)u=L∂u∂z,u=L∂u∂r,v=rM0+L∂v∂z,v=L∂v∂r,w=0,N1=n∂N1∂z,N1=n∂N1∂r,N2=n∂N2∂z,N2=n∂N2∂r,N3=rM0,T=Tw+h∂T∂z,T=Tw+h∂T∂r,(11)C=Cw+m∂C∂z,C=Cw+m∂C∂r,atz=0,(12)u⟶u∞,v⟶v∞,N⟶0,T⟶T∞,C⟶C∞,asz⟶∞,where N1,N2,N3 are angular velocity in the r, ω, and z axes, correspondingly, κ,a,b,c are material constants (viscosity coefficients), L stands for momentum slip factor, n stands for angular slip factor, h represents thermal factor, m is solutal jump factor, μ stands for the coefficient of fluid viscosity, ρ is the density, cp stands for the specific heat, g represents the gravity, j=χ/M0 is the microinertia per unit mass, κ 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 variable temperature and concentration at the surface, where ΔT and ΔC being constants 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∇ρ/ρ=σT−T∞+σtT−T∞2+σ∗C−C∞+σsC−C∞2, and σ and σ∗ are constants, σt and σs are the constant coefficients of thermal and volumetric expansion, respectively. This relation will be nonlinear density, temperature, and concentration (NDTC) variation. Y=ρcp/ρcf is the ratio between the effective heat capability of the nanoparticle material and the heat capability of the fluid and DB and DT stand for the Brownian and the thermophoretic diffusion coefficient, respectively. Let rx=dsinx/d be the radial distance from the symmetrical axis to the surface.Figure 1
Physical model and coordinate system.By using the nondimensional variables such asω=x/d, η=z2M0/χ, r=rω/d, u=rM0f′r,η, v=rM0Gr,η, w=−2M0χfr,η, χ=μ+κ/ρ, N1=rM02M0/χH1r,η, N2=rM02M0/χH2r,η,(13)N3=rM0H3r,η,θr,η=T−T∞Tw−T∞,Ωr,η=C−C∞Cw−C∞.The continuity equation can be integrated so that it satisfies equation (1). Hence, equations (2)–(9) can be written in the following nonlinear system of ODEs:(14)2f″′+B10.5H2f−2H2′+0.5f2tanω∂f∂ω−2f′−tanωf′∂f′∂ω+G2−f′2+sinωGrRe2θ+λθ2+GmRe2GmΩ+sΩ2=0,(15)2G″+B1H1′+G′f−2f′G−tanωf∂G∂ω+sinωGrRe2Grθ+λθ2+GmRe2GmΩ+sΩ2=0,(16)H1″−0.5B2G′+2H1−0.5B3f′H1+tanω∂H1∂ωf′−2fH1′−GH2=0,(17)2H2″+B2f″−2H2−B3f′H2+tanω∂H2∂ωf′−H2f+GH1=0,(18)2H3″+4B4H1′+2B5G−2H3−B6−f′tanω∂H3∂ω−fH3+2fH3′=0,(19)2Prθ″+NbΩ′θ′+2Ntθ′2+2θ′f−tanωf′∂θ∂ω=0,(20)Ω″−0.5PrNtNbNbΩ′θ′+2Ntθ′2+θ′f−tanωf′∂θ∂ω+PrLefΩ′−tanω∂Ω∂ωf′=0,with the boundary conditions.At(21)η=0:f0=0,f′0=s0f″0,∂f′0∂ω=2Re1/2sinωf″0,∂f0∂r=0,G0=1+s0G′0,∂G0∂ω=2Re1/2sinωG′0,H10=s1H1′0,∂H10∂ω=2Re1/2sinωH1′0,H20=s1H2′0,∂H20∂ω=2Re1/2sinωH2′0,H30=1,∂H30∂r=0,θ0=1+s2θ′0,∂θ0∂ω=2Re1/2sinωθ′0,Ω0=1+s3Ω′0,∂Ω0∂ω=2Re1/2sinωΩ′0,asη⟶∞:f′=G=0,θ=0,H1=H2=H3=0,Ω=0,where prime represents differential with respect to η, s0=L2M0/χ is the velocity slip parameter, s1=n2M0/χ denotes microrotation slip parameter, s2=h2M0/χ stands for thermal parameter, s3=m2M0/χ withstands for solutal parameter, (B1=κ/μ+κ, B2=κχ/cM0, B3=χρj/c, B4=a+b/a+b+c, B5=κχ/M0a+b+c, and B6=ρjχ/a+b+c are angular velocity parameters, Gr/Re2 and Gm/Re2 denote thermal and solutal mixed convection parameters (or Archimedes numbers), Gr=r3gσTw−T∞/χ2 represents thermal Grashof number, Gm=r3gσ∗Cw−C∞/χ2 represents solutal Grashof number, Re=M0r2/χ is Reynolds number, and λ=σt/σTw−T∞ and s=σs/σ∗Cw−C∞ are thermal and solutal nonlinear convection parameters, correspondingly. We observe that, for λ=s=0, the flow of equations (14) and (15) become mixed convective micropolar nanofluid. Also, the relation of buoyancy forces (Gr:Re2 and Gm:Re2) seen in the these equations are the dimensionless parameters, indicating the ratios between the buoyancy forces caused by temperature, concentration change, and the ratio inertial to viscous forces owing to kinematic viscosity change. Furthermore, Gr and Gm are zero for nonbuoyancy influence because of thermal and mass diffusions, as Re⟶∞, i.e., Gr/Re2andGm/Re2≈0 are for nonbuoyancy effect owing to momentum diffusion and when Gr/Re2=Gm/Re2=1 are for thermal, mass buoyancy forces, and inertial to viscous forces of the same strength. Pr=ρχcp/K is Prandtl number, Sc=PrLe=χ/DB is Schmidt number, Nb=ρcp/ρcfDB/χCw−C∞ stands for Brownian motion parameter, and Nt=ρcp/ρcfDT/T∞χTw−T∞ represents thermophoresis parameter.It can appear that, due to the rotational symmetry, the derivative in the radial direction may be neglected and equations (14)–(21) are reduced to the next nonlinear system of ordinary differential equations:(22)2f″′+B10.5H2f−2H2′−f2f′+G2−f′2=0,2G″+B1H1′+G′f−2f′G=0,H1″−0.5B2G′+2H1−0.5B3f′H1−2fH1′−GH2=0,2H2″+B2f″−2H2−B3f′H2−H2f+GH1=0,2H3″+4B4H1′+2B5G−2H3−B6−fH3+2fH3′=0,2Prθ″+NbΩ′θ′+2Ntθ′2+2θ′f=0,Ω″−0.5PrNtNbNbΩ′θ′+2Ntθ′2+θ′f+PrLefΩ′=0,with the boundary conditions.At(23)η=0:f0=0,f′0=s0f″0,∂f′0∂ω=∂f0∂r=0,G0=1+s0G′0,∂G0∂ω=0,H10=s1H1′0,∂H10∂ω=0,H20=s1H2′0,∂H20∂ω=0,H30=1,∂H30∂r=0,θ0=1+s2θ′0,∂θ0∂ω=0,Ω0=1+s3Ω′0,∂Ω0∂ω=0,asη⟶∞:f′=G=0,θ=0,H1=H2=H3=0,Ω=0.The physical measures of awareness in this problem are the narrow skin friction coefficientsCfr,ω, surface couple stresses mwr,ω, the Nusselt number Nu, and the Sherwood Sh number, and they can be written as follows:(24)Cfr=τzrρrM02,Cfω=τzωρrM02,Mr=mwrρM0χjr2,Mω=mwωρM0χjr2,Mr=mrρM0χjr,Nu=rqwκTw−T∞,Sh=rqmDBCw−C∞,where τzr=μ+κ∂u/∂z+∂w/∂r, τzω=μ+κ∂v/∂z+∂w/∂r, qw=−K∂T/∂zz=0, qm=DB∂C/∂zz=0, Mr=−a+b+c∂N2/∂z, Mω=c∂N1/∂z, and mr=a+b+c∂N3/∂z. Using the dimensionless variables (12) and the boundary conditions (20), the narrow skin friction coefficients surface couple stresses, the Nusselt number and the Sherwood number are obtained:(25)f″0=CfrRe2,G′0=CfωRe2,H1′0=B32Remω,−H2′0=B62Remr,−H3′0=B6Mrχ2M0,Nu12Re=−θ′0,Sh12Re=Ω′0.


## 3. Numerical Solution

Pairs of seven harmonized high order ordinary differential equations (14)–(20), subjected to the boundary conditions, equation (21), are answered numerically using the function bvp4c from matlab software for various values of physical parameters and numbers.Statistical results are found using Matlab BVP solver bvp4c from matlab which is a finite difference code that realize the three-stage Lobatto IIIa formulation. To apply bvp4c from matlab, first, equations (14)–(20) are converted into a system of first-order equations.Second, assemble a boundary value problem (bvp) and using the bvp solver in matlab to numerically solve this system, including the above boundary condition and income on a suitable finite value for the far field boundary condition, that is,η⟶∞, say η∞=1 and the step size is taken as △η=0.01, the numerical result is obtained; it is exact to the fifth decimal place as the measure of convergence. In solving the BVP by means of matlab, bvp4c has only two point of views: a function ODEs for calculation of the residual in the boundary conditions and a building solint that provides a guess for a mesh. The ODEs are handled exactly as in the Matlab IVP solvers. Further clarification on the procedure of bvp4c is found in the book by Shampine et al. [24].Lety1=f, y2=f′, y3=f″, y4=G, y5=G′, y6=H1, y7=H1′, y8=H2, y9=H2′, y10=H3, y11=H3′, y12=θ, y13=θ′, y14=Ω, y15=Ω′, y16=∂f/∂ω, y17=∂f′/∂ω, y18=∂G/∂ω, y19=∂H1/∂ω, y20=∂H2/∂ω, y21=∂H3/∂ω, y22=∂θ/∂ω, y23=∂Ω/∂ω, and y = [f, f′, f″, G, G′, H1, H1′, H2, H2′, H3, H3′, θ, θ′, Ω, Ω′, ∂f/∂ω, ∂f′/∂ω, ∂G/∂ω, ∂H1/∂ω, ∂H2/∂ω, ∂H3/∂ω, ∂θ/∂ω, ∂Ω/∂ω]T give(26)ddηy1y2y3y4y5y6y7y8y9y10y11y12y13y14y15y16y17y18y19y20y21y22y23=y2y3−0.5∗B1∗0.5∗y8∗y1−2∗y2−0.25∗y1∗y1∗B∗y16−2∗y2−y4∗y4+y2∗y2+B∗y2∗y17−0.5∗C∗Gr/Re∗Re∗y12+L∗y12∗y12+Gm/Re∗Re∗y14+s∗y14∗y14y5y2∗y4−0.5∗y1∗y5+0.5∗B∗y1∗y18−0.5∗B1∗y7−0.5∗C∗Gr/Re∗Re∗y12+L∗y12∗y12+Gm/Re∗Re∗y14+s∗y14∗y14y70.5∗B2∗y5+2∗y6+B3∗y2∗y6+B∗y19∗y2−2∗y1∗y7−y4∗y8y9−0.5∗B2∗y3−2∗y8−B3∗y2∗y8+B∗y20∗y2−y1∗y8+y4∗y6y11−2∗B4∗y7+B5∗y4−y10+0.5∗B6∗−y1∗y10−B∗y2∗y21+2∗y1∗y11y13−Pr∗0.5∗Nb∗y13∗y15−B∗y1∗y22+Nt∗y13∗y13+y1∗y13y152∗Nt∗Pr∗0.5∗Nb∗y13∗y15−B∗y1∗y22+Nt∗y13∗y13+y1∗y13−Pr∗Le∗2∗y1∗y15−B∗y23y16y17y18y19y20y21y22y23.


## 4. Results and Discussion

In this section, the results of different governing physical parameters on nondimensional velocity, temperature, concentration, skin friction and wall couple stress coefficients and confined Nusselt and Sherwood numbers have been discussed. In this study, the comparison with the literature value is not compared, since there is no related article.


### 4.1. Velocity Profiles

The dimensionless velocity profile graphs off′η,Gη for different values of angular velocity parameter B1, thermal and solutal nonlinear convection parameters (λ, s), and thermal and solutal Grashof numbers (Gr, Gm) are shown in Figures 2–11. Figures 2–9 reveal that an increase in values of Gr, Gm, λ, and s reduces the kinematic viscosity of the fluid flow which falls opposite to the flow of the fluid that results in an increase of the magnitude of radial and azimuthal velocity profiles f′ηandGη near the surface of the disk, but they decrease with an enhancement of Gm, λ, and s far away from the surface. Moreover, boost in the values of angular velocity parameter B1 raises the vortex viscosity of the fluid flow which adds opposition to flow of the fluid that results in declination of the fluid velocity profile f′η, but the reverse effects are observed on the fluid velocity profile Gη as seen in Figures 10 and 11.Figure 2
Graphf′ profile for different values of Gr when Gm=100,B1=5,B2=B5=s2=0.2,s=0,B6=0.3,B3=B4=s1=Nb=0.1,Le==Nt=λ=0.5,s=1,Re=0.01,Pr=0.72,andω=1.Figure 3
GraphGη for different values of Gr when Gm=20.Figure 4
Graphf′η profile for different values of Gm when Gr=70.Figure 5
GraphGη for different values of Gm when B1=2.Figure 6
Graphf′ profile for different values of λ when s=0.1.Figure 7
GraphGη for different values of λ.Figure 8
Graphf′ profile for different values of s when Gm=9,Gr=80,andλ=0.1.Figure 9
GraphGη for different values of s when Gm=9,Gr=80,andλ=0.1.Figure 10
Graphf′ profile for different values of B1 when Gm=0.9andGr=80.Figure 11
GraphGη for different values of B1 when Gm=9andGr=80.


### 4.2. Temperature and Concentration Profiles

The influences of thermal and solutal jump parameterss2 and s3, nonlinear convection parameters (λ, s), and thermal and solutal Grashof numbers (Gr, Gm) on temperature and concentration sketches are prearranged in Figures 12–19. These figures indicate that the values of temperature and concentration disseminations and their boundary layer thickness decrease with increase of thermal and solutal jump parameters s2 and s3, nonlinear convection parameters (λ, s), and thermal and solutal Grashof numbers (Gr, Gm). These effects happen due to the decline in the kinematic viscosity of the fluid which reduces thermal diffusion, that leads to a decrease in thermal and solutal boundary layer thicknesses.Figure 12
Graphθη profile for different values of Gr when Gm=80, B2=B5=s2=0.2, s0=0, B6=3, B3=s1=s=Nb=0.1, B4=5, s3=Le==Nt=λ=0.5, Re=0.01, Pr=0.72, and B1=ω=1.Figure 13
Graphθη for different values of Gm when Gr=10.Figure 14
Graphθη profile for different values of λ when s=0.5.Figure 15
GraphΩη for different values of λ when Gr=20.Figure 16
Graphθη profile for different values of s when Gr=1andλ=0.5.Figure 17
GraphΩη for different values of s when Nb=0.5andNt=λ=0.1.Figure 18
Graphθη profile for different values of s2 when Gm=0.1.Figure 19
GraphΩη for different values of s3 when Gr=1.


### 4.3. Skin Frictions and Wall Couple Stresses

The impacts of thermal and solutal nonlinear convection parameters (λ, s), thermal and solutal Grashof numbers (Gr, Gm), angular velocity parameters (B1–B5) on skin friction coefficients (f″0, G′0), and wall duos coefficients (H′10,−H′20,−H′30) sketches are prearranged in Figures 20–35. Table 1 and Figures 20 and 29 illustrate that upsurge in values of B1 reduces the kinematic viscosity of the fluid flow which falls opposite to the flow of the fluid that results in decline of the skin friction coefficients f″0 and wall duos coefficients -H2′0; Table 1 demonstrates that increase in velocity slip parameter s0 reduces skin friction coefficients f″0,G′0, whereas the opposite effects are observed on the skin friction coefficients G′0 and wall duos coefficients H1′0. Moreover, Figures 22–27 show that boost in values of λ, s, Gr, and Gm raise density of the fluid flow which increases resistance to the flow of the result upsurge skin friction coefficients f″0andG′0. Similarly, Figures 30 and 31 indicate that an increase in values of B2 improves vortex viscosity of the fluid flow which enlarges resistance to rotate the fluid that results in upsurge wall duos coefficients H1′0 and -H2′0. Furthermore, increase in values of B3, improve wall couple coefficients -H2′0 but bring down wall couple coefficients H1′0 as indicated in Figures 32 and 33. Next, wall couple coefficients -H2′0 reduces with raise values of B4 and B5 in the boundary layer, as shown in Figures 34 and 35.Figure 20
Graphf″0 profile for different values of B1 when Gm=Gr=10, B2=B3=B4=s=0.5, s0=0, B6=s2=0.2, B5=s3=λ=0.1, Nb=0.5.s3=Le=0.7, Nt=2, Re=0.01, Pr=0.6, and ω=1.Figure 21
GraphG′0 for different values of B1.Figure 22
Graphf″0 profile for different values of λ.Figure 23
GraphG′0 for different values of λ.Figure 24
Graphf″0 profile for different values of s.Figure 25
GraphG′0 for different values of s.Figure 26
Graphf″0 profile for different values of Gr.Figure 27
Graphf″0 for different values of Gm.Figure 28
GraphH1′0 profile for different values of B1 when Gm=Gr=10, B2=B3=B4=s=0.5, s0=0, B6=s2=0.2, B5=s3=λ=0.1, Nb=0.5.s3=Le=0.7, Nt=2, Re=0.01, Pr=0.6, and ω=1.Figure 29
Graph -H2′0 for different values of B1.Figure 30
GraphH1′0 profile for different values of B2.Figure 31
Graph -H2′0 for different values of B2.Figure 32
GraphH1′0 profile for different values of B3 when Gr = Gm = 1.Figure 33
Graph -H2′0 for different values of B3 when Gr = 1 and Gm = 12.Figure 34
Graph -H3′0 profile for different values of B4.Figure 35
Graph -H3′0 for different values of B5.Table 1
The computed values of skin friction coefficientsf″0 and G′0, when B3=0.2, B4=B5=Gr=Gm=s2=s=λ=Nt=0.1, B6=0.3, Re=0.01, Nb=Le=s3=0.5, Pr=0.72, and ω=1.
B1B2s0f″0G′00.10.10.10.81160.76180.20.10.10.81040.76240.30.10.10.80930.76310.130.10.76420.150.10.76560.170.10.76680.10.10.30.25710.53040.10.10.40.21520.45710.10.10.50.18500.4002


### 4.4. Nusselt and Sherwood Numbers

The impacts of Prandtl Pr, Lewis Le numbers, thermal, solutal (s2, s3), and Brownian motion Nb parameters on Nusselt −θ′0, and Sherwood numbers Ω′0 are presented in Figures 36–41. Figure 36 illustrates that boost in values of Pr bring down thermal diffusion of the temperature that cause upsurge Nusselt number −θ′0 near the surface of the disk. Figure 37 demonstrates that a rise in the values of Le reduces the density of the fluid , as a result the Sherwood number Ω′0 increased, whereas an increment in the values s2 and s3 reduces the Nusselt number −θ′0 and Sherwood number Ω′0 as demonstrated in Figures 38-39. Figures 40 and 41 show the impact of Brownian motion parameter Nb on Nusselt number −θ′0 and Sherwood number Ω′0 . From the figures, it is possible to see that Nb favors the Nusselt number and opposite effect on the Sherwood number.Figure 36
Graph -θ′0 profile for different values of Pr when Gm=10, Gr=B2=B3=5, s0=0, B6=0.2, B5=s2=λ=Gr=0.1, Nb=B4=s=s3=Le=0.5, Nt=2, Re=0.01, and B1=ω=1.Figure 37
GraphΩ′0 for different values of Le.Figure 38
Graph -θ′0 profile for different values of s2.Figure 39
GraphΩ′0 for different values of s3.Figure 40
Graph -θ′0 profile for different values of Nb.Figure 41
GraphΩ′0 for different values of Nb.


## 4.1. Velocity Profiles

The dimensionless velocity profile graphs off′η,Gη for different values of angular velocity parameter B1, thermal and solutal nonlinear convection parameters (λ, s), and thermal and solutal Grashof numbers (Gr, Gm) are shown in Figures 2–11. Figures 2–9 reveal that an increase in values of Gr, Gm, λ, and s reduces the kinematic viscosity of the fluid flow which falls opposite to the flow of the fluid that results in an increase of the magnitude of radial and azimuthal velocity profiles f′ηandGη near the surface of the disk, but they decrease with an enhancement of Gm, λ, and s far away from the surface. Moreover, boost in the values of angular velocity parameter B1 raises the vortex viscosity of the fluid flow which adds opposition to flow of the fluid that results in declination of the fluid velocity profile f′η, but the reverse effects are observed on the fluid velocity profile Gη as seen in Figures 10 and 11.Figure 2
Graphf′ profile for different values of Gr when Gm=100,B1=5,B2=B5=s2=0.2,s=0,B6=0.3,B3=B4=s1=Nb=0.1,Le==Nt=λ=0.5,s=1,Re=0.01,Pr=0.72,andω=1.Figure 3
GraphGη for different values of Gr when Gm=20.Figure 4
Graphf′η profile for different values of Gm when Gr=70.Figure 5
GraphGη for different values of Gm when B1=2.Figure 6
Graphf′ profile for different values of λ when s=0.1.Figure 7
GraphGη for different values of λ.Figure 8
Graphf′ profile for different values of s when Gm=9,Gr=80,andλ=0.1.Figure 9
GraphGη for different values of s when Gm=9,Gr=80,andλ=0.1.Figure 10
Graphf′ profile for different values of B1 when Gm=0.9andGr=80.Figure 11
GraphGη for different values of B1 when Gm=9andGr=80.


## 4.2. Temperature and Concentration Profiles

The influences of thermal and solutal jump parameterss2 and s3, nonlinear convection parameters (λ, s), and thermal and solutal Grashof numbers (Gr, Gm) on temperature and concentration sketches are prearranged in Figures 12–19. These figures indicate that the values of temperature and concentration disseminations and their boundary layer thickness decrease with increase of thermal and solutal jump parameters s2 and s3, nonlinear convection parameters (λ, s), and thermal and solutal Grashof numbers (Gr, Gm). These effects happen due to the decline in the kinematic viscosity of the fluid which reduces thermal diffusion, that leads to a decrease in thermal and solutal boundary layer thicknesses.Figure 12
Graphθη profile for different values of Gr when Gm=80, B2=B5=s2=0.2, s0=0, B6=3, B3=s1=s=Nb=0.1, B4=5, s3=Le==Nt=λ=0.5, Re=0.01, Pr=0.72, and B1=ω=1.Figure 13
Graphθη for different values of Gm when Gr=10.Figure 14
Graphθη profile for different values of λ when s=0.5.Figure 15
GraphΩη for different values of λ when Gr=20.Figure 16
Graphθη profile for different values of s when Gr=1andλ=0.5.Figure 17
GraphΩη for different values of s when Nb=0.5andNt=λ=0.1.Figure 18
Graphθη profile for different values of s2 when Gm=0.1.Figure 19
GraphΩη for different values of s3 when Gr=1.


## 4.3. Skin Frictions and Wall Couple Stresses

The impacts of thermal and solutal nonlinear convection parameters (λ, s), thermal and solutal Grashof numbers (Gr, Gm), angular velocity parameters (B1–B5) on skin friction coefficients (f″0, G′0), and wall duos coefficients (H′10,−H′20,−H′30) sketches are prearranged in Figures 20–35. Table 1 and Figures 20 and 29 illustrate that upsurge in values of B1 reduces the kinematic viscosity of the fluid flow which falls opposite to the flow of the fluid that results in decline of the skin friction coefficients f″0 and wall duos coefficients -H2′0; Table 1 demonstrates that increase in velocity slip parameter s0 reduces skin friction coefficients f″0,G′0, whereas the opposite effects are observed on the skin friction coefficients G′0 and wall duos coefficients H1′0. Moreover, Figures 22–27 show that boost in values of λ, s, Gr, and Gm raise density of the fluid flow which increases resistance to the flow of the result upsurge skin friction coefficients f″0andG′0. Similarly, Figures 30 and 31 indicate that an increase in values of B2 improves vortex viscosity of the fluid flow which enlarges resistance to rotate the fluid that results in upsurge wall duos coefficients H1′0 and -H2′0. Furthermore, increase in values of B3, improve wall couple coefficients -H2′0 but bring down wall couple coefficients H1′0 as indicated in Figures 32 and 33. Next, wall couple coefficients -H2′0 reduces with raise values of B4 and B5 in the boundary layer, as shown in Figures 34 and 35.Figure 20
Graphf″0 profile for different values of B1 when Gm=Gr=10, B2=B3=B4=s=0.5, s0=0, B6=s2=0.2, B5=s3=λ=0.1, Nb=0.5.s3=Le=0.7, Nt=2, Re=0.01, Pr=0.6, and ω=1.Figure 21
GraphG′0 for different values of B1.Figure 22
Graphf″0 profile for different values of λ.Figure 23
GraphG′0 for different values of λ.Figure 24
Graphf″0 profile for different values of s.Figure 25
GraphG′0 for different values of s.Figure 26
Graphf″0 profile for different values of Gr.Figure 27
Graphf″0 for different values of Gm.Figure 28
GraphH1′0 profile for different values of B1 when Gm=Gr=10, B2=B3=B4=s=0.5, s0=0, B6=s2=0.2, B5=s3=λ=0.1, Nb=0.5.s3=Le=0.7, Nt=2, Re=0.01, Pr=0.6, and ω=1.Figure 29
Graph -H2′0 for different values of B1.Figure 30
GraphH1′0 profile for different values of B2.Figure 31
Graph -H2′0 for different values of B2.Figure 32
GraphH1′0 profile for different values of B3 when Gr = Gm = 1.Figure 33
Graph -H2′0 for different values of B3 when Gr = 1 and Gm = 12.Figure 34
Graph -H3′0 profile for different values of B4.Figure 35
Graph -H3′0 for different values of B5.Table 1
The computed values of skin friction coefficientsf″0 and G′0, when B3=0.2, B4=B5=Gr=Gm=s2=s=λ=Nt=0.1, B6=0.3, Re=0.01, Nb=Le=s3=0.5, Pr=0.72, and ω=1.
B1B2s0f″0G′00.10.10.10.81160.76180.20.10.10.81040.76240.30.10.10.80930.76310.130.10.76420.150.10.76560.170.10.76680.10.10.30.25710.53040.10.10.40.21520.45710.10.10.50.18500.4002


## 4.4. Nusselt and Sherwood Numbers

The impacts of Prandtl Pr, Lewis Le numbers, thermal, solutal (s2, s3), and Brownian motion Nb parameters on Nusselt −θ′0, and Sherwood numbers Ω′0 are presented in Figures 36–41. Figure 36 illustrates that boost in values of Pr bring down thermal diffusion of the temperature that cause upsurge Nusselt number −θ′0 near the surface of the disk. Figure 37 demonstrates that a rise in the values of Le reduces the density of the fluid , as a result the Sherwood number Ω′0 increased, whereas an increment in the values s2 and s3 reduces the Nusselt number −θ′0 and Sherwood number Ω′0 as demonstrated in Figures 38-39. Figures 40 and 41 show the impact of Brownian motion parameter Nb on Nusselt number −θ′0 and Sherwood number Ω′0 . From the figures, it is possible to see that Nb favors the Nusselt number and opposite effect on the Sherwood number.Figure 36
Graph -θ′0 profile for different values of Pr when Gm=10, Gr=B2=B3=5, s0=0, B6=0.2, B5=s2=λ=Gr=0.1, Nb=B4=s=s3=Le=0.5, Nt=2, Re=0.01, and B1=ω=1.Figure 37
GraphΩ′0 for different values of Le.Figure 38
Graph -θ′0 profile for different values of s2.Figure 39
GraphΩ′0 for different values of s3.Figure 40
Graph -θ′0 profile for different values of Nb.Figure 41
GraphΩ′0 for different values of Nb.


## 5. Conclusions

This report considers the effects of thermal and solutal nonlinear convection as well as manifold slip conditions on dimensional boundary layer flow of a micropolar nanofluid due to a rotating disk were analysed. The boundary layer equations of the flow problem are reduced into a pairs of high-order nonlinear ordinary differential equations by using the similarity transformation. The obtained differential equations are solved numerically using bvp4c from matlab software. Numerical outcomes are acquired for different main parameters of the flow problem. The effects of governing parameters are presented by means of figures and tables. The main findings are(1)
An increasing in the values of thermal and solutal nonlinear convection parameters allow to increase in the value of velocitiesf′η and Gη near surface of the disk and reduce at far from the disk as well as thermal and solutal Grashof numbers tolerate to increase the value of radial velocity f′η near surface of the disk.(2)
Velocity profileGη is positive throughout its boundary, indicating in that way there is outward flow from the disk which reaches the particular point and it drops to zero when far from the surface of the disk.(3)
The existence of the thermal and solutal jump, nonlinear convection parameters, and Grashof numbers allows to decline the temperature and concentration distributions near the disk surface.(4)
Radial skin friction coefficients exhibit increasing behavior for increasing thermal and solutal nonlinear convection parameters and Grashof numbers; skin friction coefficientsG′0 exhibit increasing behavior for increasing thermal and solutal nonlinear convection parameters.(5)
An increase in quantities of Pr andLe causes enhancement in the Nusselt and Sherwood numbers while they fall as magnitudes of s2 and s3 increase, accordingly.(6)
The value of wall couple coefficients (H1′0,−H2′0) increase with B2. However, they rotate in opposite directions.(7)
An improvement in the values of angular velocity parameters (B1, B3) allows to rotate the wall couple stresses (H1′0 and H2′0) in opposite directions.



---

*Source: 4735650-2020-06-17.xml*