# Triple Solutions with Stability Analysis of MHD Mixed Convection Flow of Micropolar Nanofluid with Radiation Effect

Authors: Hazoor Bux Lanjwani; Muhammad Saleem Chandio; M. Imran Anwar; Amnah S. Al-Johani; Ilyas Khan; Md. Nur Alam
Journal: Journal of Nanomaterials
 (2022)

---

## Abstract

This paper deals with two-dimensional steady boundary layer flow, heat, and mass transfer characteristics of micropolar nanofluid past on exponentially stretching/shrinking surface. The effect of different physical parameters like magnetic field, buoyancy, thermal radiation, and connective heat transfer are examined. Furthermore, similarity solutions are obtained by similarity transformation on the governing system of partial differential equations. The shooting method with help of the Maple software is used to achieve the numerical solutions of the equations. For the different ranges of the applied parameters, triple solutions are obtained for both cases of the surface. In view of the triple solutions, stability analysis is performed by bvp4c in the MATLAB software, where only first solution is found feasible which is discussed. The main findings of the first solution indicate the skin friction, drag force, heat, and mass transfer rates are increasing for theλ>0 and decreasing for λ<0 as the K is enhanced. The velocity profiles decrease with increase in magnetic, slip velocity, and suction parameters. The temperature profiles increase with increase in magnetic, thermophoresis, thermal radiation, and Brownian motion parameters, whereas concentration profiles reduce with increase in Schmitt number and Brownian motion.


---

## Body

## 1. Introduction

The study of non-Newtonian fluid flows have gained much importance, because of the common Newtonian fluids may not completely satisfy the properties of the fluid flow in many industrial applications, examples of such fluids are biological fluids, polymeric fluids, fluids containing additives, liquid crystals, and paint colloidal solutions. Moreover, the class of the non-Newtonian fluids containing different kinds of complex properties are Casson fluids, Maxwell fluids, and micropolar fluids. The micropolar fluid introduced by Eringen [1] possesses a microscopic effect due to the microstructure and micromotion of particles present in fluid. These microstructure particles are of different shapes which rotate independently to the motion of the fluid perticles (Anwar et al. [2]). The system of the micropolar fluid flow equations contains a microrotating vector besides the classical velocity vector. These fluids contain smaller rigid particles which rotate about the centroid of volume particles that predicts the flow behaviors at rotation and microscale independently that is defined by a microrotation vector. Therefore, micropolar fluids are very important in fluid dynamics, especially in studying some flows around some important surfaces such as stretching surfaces or shrinking surfaces. In this regard, MHD micropolar fluid flow on the inclined plate was investigated by Kasim et al. [3]. The micropolar fluid flow on the inclined surface with different physical parameters was also studied by Das [4]. Srinivasacharya and Bindu [5] examined the entropic generation of the micropolar fluid with parallel plates on the inclined channel. Rahman et al. [6] examined the heat flux process in the micropolar fluid with incorporation different kinds of the fluid properties. Dero et al. [7] examined MHD micropolar nanofluid flow on exponential shrinking/stretching surface with radiation, Brownian motion, and thermophoresis effects. Srinivasacharya et al. [8] examined double diffusion impacts on the micropolar fluid flow on a slanted sheet. In [9], the authors studied peristaltic flow of nanofluids through a tapered channel with a porous medium.Whereas magnetohydrodynamics (MHD) is concerned with physical and mathematical scaffold which shows magnetic dynamics in the electrical conducting fluid, the applications of magnetohydrodynamics are especially used in the modern industrial and engineering areas, such as in drawing of the plastic wires and films, polymer extrusion in melt spinning process, crystal growth, paper production and glass fiber, fluid film in the condensations processes, manufacturing in foods, electronic chips, electrochemical process, flow through the filtering devices, and thermal energy storage. Several investigators have worked on magnetohydrodynamics flow overstretching surfaces by using different physical parameters [10].Moreover, nanotechnology has achieved much attention of the investigators because of having many applications in industrial and biological sciences. In industries, clay nanocomposites are used for the production of the impermeable wrapping films, silver nanoparticles are particularly used for productions of bins storage, while particles of the zinc oxides can be used for protection of the materials from the ultraviolet radiations and the carbon nanotubes pervaded graphite are used for the productions of light tennis rackets. In biological science, the investigators are busy developing the nanocapsules for the injections replacement that can be passed through the stomach and can accessed in the blood stream easily, whereas nanoparticles’ antibacterial properties make them fit to cut and stitch the bandages. The nanofluid is a subclass of such a widely growing research field. Actually, nanofluids are engineered manufacturing fluids that are made by suspensions of the smallest sized particles (10−9 to 10−11 m) in common traditional fluids. Choi [11] called such flu as the fluids because of the addition of the smallest nanosized particles. Numerous engineering problems are studied by many researchers for more than two decades on the progress of nanofluids [12, 13], whereas the extensive convective rate of the heat transfers in nanofluids was studied by Buongiorno [14] by developing a two-phase nanofluid model. Furthermore, the boundary layer nanofluid flow phenomena on the shrinking/stretching surfaces have been discussed by many researchers [15, 16]. The nanofluid stagnation point flow was considered by Nazar et al. [17]. Kuznetsov and Nield [18] considered the boundary layer flow of the nanofluid through a vertically moving plate with buoyancy effect, while Chamkha et al. [19] considered mixed convective magnetohydrodynamics nanofluid flow with buoyancy effect by using Buongiorno’s (2006) model.Motivated from the abovementioned work, two-dimensional laminar magnetohydrodynamic boundary layer flow of micropolar nanofluid on exponentially shrinking/stretching surface with mixed convection, thermal radiation, magnetic, Brownian motion, and thermophoresis parameters has been studied. Furthermore, convective heat transfer and porous medium are also taken into account. There has been extension on the work done by Bidin and Nazar [20]. Their work was on the viscous fluid with unique solution while the present work is much updated. The graphical achieved results are used to study the influence of various physical parameters on the velocity, microrotation, temperature, and concentration profiles. Few other interesting studies in this direction are given in [21–24].To the best of the authors’ knowledge, such an attempt on MHD mixed convection flow of a micropolar nanofluid in the existence of the thermal radiation and magnetic field on an exponentially vertical stretching/shrinking surface by obtaining triple solutions with stability analysis has not been tried before. The results are demonstrated graphically with detailed discussion. It can be expected that the present research work will prove helpful for the researchers which are interested to study nanofluids particularly micropolar-based nanofluids as well as the importance of the stability analysis in case of the occurrence of multiple solutions, where stability analysis decides about real and extraneous solutions.


## 2. Problem Formulation

There is considered two-dimensional laminar magnetohydrodynamic (MHD) mixed convection flow, heat, and the mass transfer of the micropolar nanofluid on the permeable exponentially vertical shrinking/stretching surface with radiation effect. An electrically conducting nanofluid is considered with influence of the magnetic fieldBx which is used perpendicular to the surface. The value of the Reynolds number is negligible that is ignored. The flow is considered along the x-axis perpendicular to the y-axis. The surface velocity is considered as uwx=Uwex/L at y=0. A graphical presentation of flow model along the coordinate system is presented in Figure 1.Figure 1
The flow model and the coordinate system.Applying boundary layer approximations, the continuity, momentum, microrotation, energy, and concentration equations for the micropolar nanofluid flow of present study are written as(1)∂u∂x+∂v∂y=0,(2)u∂u∂x+v∂u∂y=ϑ+K1ρ∂2u∂y2+K1ρ∂N∂y−σBo2xuρ+g∗βTT−T∞+g∗βcC−C∞,(3)u∂N∂x+v∂N∂y=1ρjγ∗∂2N∂y2−K12N+∂u∂y,(4)u∂T∂x+v∂T∂y=α∂2T∂y2+τwDB∂C∂y∂T∂y+DTT∞∂T∂y2−ϑcp∂qr∂y,(5)u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2,where velocity components in the directions of the x and yare symbolized by u and v, respectively.σ is an electrical conduction,cp is the specific heat of the constant pressure, ϑ is the kinematic viscosity, Box is the constant of the magnetic field Bx, βT is the thermal expansion coefficient, βCis concentration expansion coefficient, and ρ is the density of the fluid. N is the microrotation, K1 is the vortex viscosity, γ∗ represents the spin gradient viscosity, j is the ratio of microinertia, α is thermal diffusion, and kthermal conduction of the fluid. m is the constant whose range is 0≤m≤1. Where at m=0, there will be N=0 that represents strong concentration that show no rotation due to microelements nearest to the solid surface, while m=1/2 represents weaker concentration. Furthermore, m=1 shows turbulences in flows of fluid (Ishak et al. [21]). Several researchers have taken γ∗=μ+K1/2j=μ1+K/2j,where K1=μK is the material parameter and =2lve−x/l/Uw. T is the temperature of the boundary layer flowing fluid, Tw is the surface temperature, T∞ is the ambient fluid temperature,C is the nanoparticle volumetric fraction, qr is the radiative thermal flux, and DBandDT are Brownian and thermophoresis diffusions coefficients, respectively. The radiative heat flux qr is written by
(6)qr=4σ∗3k∗∂T4∂y.Here,k∗ is a mean absorptions coefficient and σ∗ is the Stefan Boltzmann constant. At smaller difference in the temperature of the flow, the T4 will be described as the function of T. Using Taylor’s series to T∞ and by ignoring higher terms due to smallest values, we have
(7)T4≅4T∞3T−3T4.By using Equations (5) and (6) in Equation (3), we get
(8)u∂T∂x+v∂T∂y=α1+4.Rd3∂2T∂y2+τwDB∂C∂y∂T∂y+DTT∞∂T∂y2.The considered boundary conditions of this problem are(9)v=vwx;u=λuwx+Adudy,N=−m∂u∂y,T=Twx=T∞+T0ex/2l,,C=Cwx=C∞+C0ex/2laty=0,u⟶0;N⟶0;T⟶T∞;C⟶C∞asy⟶∞.Magnetic fieldBx for exponential surface is written as
(10)B=B0ex/l.Furthermore, the following exponential type of the similarity transformation is applied to achieve the similarity solutions:(11)ψ=2ϑlUwex/2Lfη,(12)N=Uwe3x/2LUw2ϑLgη,(13)θη=T−T∞Tw−T∞,(14)∅η=C−C∞Cw−C∞,(15)η=yUw2ϑLex/2L,here, ψx,y is an stream function that in velocity form is defined as u=∂ψ/∂yand v=−∂ψ/∂x. Thus, we get
(16)u=Uwex/Lf′η,v=−Uwϑ2Lex/2Lfη+ηf′η.Whereas using Equation (15) in Equations (2) to (5), we get following system of equations:
(17)1+Kf′′′+ff′′−2f′2+Kg′−2Mf′+2λTθ+2λCϕ=0,(18)1+K2g′′+fg′−3gf′−2Kg−Kf′′=0,(19)1Pr1+43Rdθ′′+fθ′+Nb∅′θ′+Ntθ′2=0,(20)ϕ′′+Scfϕ′+NtNbθ′′=0.Boundary conditions take the form:(21)f0=S,f′0=λ+δf′′0,g0=−mf′′0,θ′0=1,ϕ0=1,(22)f′η⟶0,gη⟶0,θη⟶0,ϕη⟶0asη⟶∞,here, prime represents derivative along the η, K=K1/μ is the non-Newtonian parameter, Pr=ϑ/α is the Prandtl number,Rd=4σ∗T∞3/kk∗ is the radiation parameter, Sc=ϑ/DB is the Schmidt number, Nt=τwDTTw−T∞/νT∞ is the thermophoresis parameter, Nb=τwDBCw−C∞/ν is the parameter of Brownian motion, λ is the shrinking/stretching surface, here λ<0 denotes shrinking case and λ>0 denotes stretching case and S>0 denotes the suction of mass by the surface due to porosity, and δ=AUw/ϑ is a velocity slip parameter.λT=Gr/Rex2=g∗βTTw−T∞/uw2 where λT is ratio between buoyancy and inertia forces that is used as the criterion to create dominant regions of the forced and free convections. The λT=1 stands for mixed (free and forced) convection and λT=Gr/Rex2=g∗βTTw−T∞/uw2 stands for mass convective parameters. The skin friction coefficient Cf, Nusselt number Nux, and the Sherwood number Sh are written as
(23)Cf=μ+K∂u/∂y+KNy=0ρuw2,Nux=−x∂T/∂yy=0Tw−T∞,Sh=−x∂C/∂yy=0Cw−C∞,using Equation (15), we get
(24)CfRex1/22lx=1+1−mKf′′0,NuRex−1/2=−θ′0,ShRex−1/2=−∅′0,whereRex=uwx/ϑ is the Reynolds number.


## 3. Stability Analysis

To check the stability, there are taken unsteady case for Equations (2)–(5) with addition of the new time dependent variableτ as mentioned by Roşca and Pop [22].
(25)∂u∂t+u∂u∂x+v∂u∂y=ϑ+K1ρ∂2u∂y2+K1ρ∂N∂y−σB2uρ+g∗βCT−T∞+g∗βCC−C∞,(26)∂N∂t+u∂N∂x+v∂N∂y=1ρjγ∗∂2N∂y2−κ2N+∂u∂y,(27)∂T∂t+u∂T∂x+v∂T∂y=α+16σ∗T∞33K∗ρcp∂2T∂y2+τwDB∂C∂y∂T∂y+DTT∞∂T∂y2,(28)∂C∂t+u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2.To get the similarity equations, the similarity transformation (15) takes the form as follows:
(29)ψ=2ϑlUwex/2lfη,τ;N=Uwe3x/2lUw2ϑlgη,τ,(30)θη,τ=T−T∞Tw−T∞,(31)∅η,τ=C−C∞Cw−C∞,(32)η=yUw2ϑlex/2l,(33)τ=Uw2lex/l.t.Using Equation (33) in Equations (25)–(28),the following is obtained:
(34)1+K∂3f∂η3−2∂f∂η2−2τ∂f∂η∂2f∂η∂τ+f∂2f∂η2+K∂g∂η−M∂f∂η+2λTθη+2λC∅η−∂2f∂η∂τ=0,(35)1+K2∂2g∂η2+f∂g∂η−3∂f∂η.g−2τ∂f∂η.∂g∂τ−2Kg−K∂2f∂η2−∂g∂τ=0,(36)1Pr1+4Rd3∂2θ∂η2−2τ∂f∂η∂θ∂τ+f∂θ∂η+Nb∂∅∂η∂θ∂η+Nt∂θ∂η2−∂θ∂τ=0,(37)∂2∅∂η2+Scf∂∅∂η−2τSc∂f∂η∂∅∂τ+NtNb∂2θ∂η2−∂∅∂τ=0.The boundary conditions given in Equation (21) are written as
(38)f0,τ=S,(39)∂f0,τ∂η=λ+δ∂2f0,τ∂η2,(40)g0,τ=−m∂2f0,τ∂η2,(41)∂θ0,τ∂η=1,∅0,τ=1,(42)∂fη,τ∂η⟶0,gη,τ⟶0,θη,τ⟶0,∅0,τ⟶0asη⟶∞.For the stability of steady flow solutionsfη=f0η, gη=g0η, θη=θ0η, and ∅η=∅0ηwhich placate the problem given in Equations (17)–(21), we further define
(43)fη,τ−e−γτFη,τ=f0η,gη,τ−e−γτGη,τ=g0η,θη,τ−e−γτHη,τ=θ0η,∅η,τ−e−γτSη,τ=∅0η,whereγ denotes unknown eigenvalues and Fη,τ, Gη,τ, Hη,τ, and Sη,τ are smaller functions related to the f0η, g0η, θ0η, and ∅0η, respectively. By substitution of Equation (43) in Equations (34)–(42), we have
(44)1+KF0′′′+f0F0′′+f0′′F0−4f0′F0′+KG0′−MF0′+2λTH0+2λCS0+2γτf0′G0+γF0′=0,1+K2G0′′+f0G0′+F0g0′−3f0′G0−3g0F0′−2KG0−KF0′′+2γτf0′G0+γG0=0,1Pr1+43RdH0′′+f0H0′+F0θ0′+Nb∅0′H0′+NbS0′θ0′+Ntθ0′H0′+2γτf′0H0+γH0=0,S0′′+Scf0S0′+F0∅0′+NtNbH0′′+2Scγτf′0S0+γS0=0.For to be achieved steady solution, there will be takenτ=0, so, we get
(45)1+KF0′′′+f0F0′′+f0′′F0−4f0′F0′+KG0′−MF0′+2λTH0+2λCS0+γF0′=0,(46)1+K2G0′′+f0G0′+F0g0′−3f0′G0−3g0F0′−2KG0−KF0′′+γG0=0,(47)1Pr1+43RdH0′′+f0H0′+F0θ0′+Nb∅0′H0′+NbS0′θ0′+Ntθ0′H0′+γH0=0,(48)S0′′+Scf0S0′+F0∅0′+NtNbH0′′+γS0=0,Boundary conditions are(49)F00=0;F0′0=δF0′′0;F0′η⟶0asη⟶∞,(50)G00=−mF0′′0;G0η⟶0asη⟶∞,(51)H00=0;H0η⟶0asη⟶∞,(52)S00=0;S0η⟶0asη⟶∞.Solutions of the eigenvalue problem of Equations (45)–(48) that satisfy Equations (49) are achieved by the solver functions of bvp4c through MATLAB software. These solutions provide a set of eigenvalues γ1<γ2<γ3. According to Harris et al. [23], the required smallest eigenvalue can be obtained by relaxing one of boundary condition to initial condition. In this regard, in the present problem, we have relaxed F0′η⟶0 as η⟶∞ in the form of the F0′′0=1. Where the smallest positive values of the γ1 demonstrates initial deterioration of the disturbance, the solutions concerned to the flow possessing smallest positive eigenvalues are considered physically feasible and stable. Converse to it, the negative γ1 demonstrates initial growth of the disorder; therefore, solutions of the fluid flow are considered not practically feasible and unstable.


## 4. Numerical Method

There are presented a short view of the shooting technique (Meade et al. [24]) that is used in Equations (17)–(20) that follow the boundary conditions given in Equation (21). Actually, this method converts the BVPs into the IVPs. So, we have
(53)f′=FP,f′′=FPP,1+KFPP′+FFPP−2FP2+KGP−2Mf′+2λTθ+2λCϕ=0,g′=GP,1+K2GP′+FGP−3GFP−2KG−KFPP=0,θ′=θP,ϕ′=ϕP,1Pr1+43RdθP′+FθP+NbϕPθP+NtθP2=0ϕ′=ϕP,ϕP′+ScFϕP+NtNbθP′=0.The boundary conditions are(54)F0=S;FP0=λ+δh1,FPP0=h1,G0=−mh1,GP0=h2,θP=1,θP0=h3,ϕ0=1;ϕP0=h4,where h1, h2, h3, and h4 are the unknown applied initial conditions. Therefore, shooting values of h1, h2, h3, and h4 are carefully taken which satisfy the specified boundary conditions of the problem properly. Furthermore, it is noticed that computations of this problem are done by the help of Maple software where shootlib function is added.


## 5. Result and Discussion

Equations (17)–(20) that subject to the boundary conditions given in Equations (21) are solved numerically by a shooting method in Maple software for the distinct values of used physical parameters. Triple solutions are obtained for different initial guesses at various ranges of specified parameters. Furthermore, for authentication of results, our obtained results have been compared with previously obtained results by the Bidin and Nazar [20], where the comparison indicates a good agreement with our obtained results which are shown in Table 1.Table 1
Comparison for the values of–θ0 at different values of Pr, where Nt, Rd, and Nb are ignored.
Pr−θ0Bidin and Nazar [20]Present results10.95470.955121.47141.471331.86911.8692The numerical computations demonstrate the effects of different used physical parameters on skin friction coefficientf′′0,couple stress coefficientg′0, local Nusselt number −θ′0, and Sherwood number −ϕ′0 as well as on the velocity, microrotation, temperature, and concentration profiles of the micropolar nanofluid which are graphically presented in Figures 2–22. Stability analysis is performed due to occurrence of triple solutions. Therefore, stability analysis is performed by bvp4c in MATLAB programming. The achieved eigenvalues are presented in Table 2, while result of the stability shows that the first solution is the stable and the physically feasible which possesses smallest positive eigenvalues, while negative eigenvalues are found for two other solutions which are declared as unstable and not physically feasible. Therefore, in this article, the result of only first solutions which is physically feasible solutions are presented.Figure 2
The effects of theK and λ on skin friction f′′0.Figure 3
The effects of theK and λ on the couple stress g′0.Figure 4
The effects of theK and λ on heat transfer rate −θ′0.Figure 5
The effects of theK and λ on mass transfer rate −ϕ′0.Figure 6
The effects of theM on velocity profiles f′η.Figure 7
The effects of theλT on the velocity profiles f′η.Figure 8
The effects of theS on the velocity profiles f′η.Figure 9
The effects of theδ on the velocity profiles f′η.Figure 10
The effects of theλ on the velocity profiles f′η.Figure 11
The effects of theλ on the microrotation profiles gη.Figure 12
The effects of theS on the microrotation profilesgη.Figure 13
The effects of theδ on the microrotation profilesgη.Figure 14
The effects of theM on the temperature profiles θη.Figure 15
The effects of theS on the temperature profiles θη.Figure 16
The effects of thePr on the temperature profilesθη.Figure 17
The effects of theRd on the temperature profilesθη.Figure 18
The effects of theNb on the temperature profilesθη.Figure 19
The effects of theNt on the temperature profilesθη.Figure 20
The effects of theS on the concentration profilesϕη.Figure 21
The effects of theNb on the concentration profilesϕη.Figure 22
The effects of theNt on the concentration profilesϕη.Table 2
The smallest obtained eigenvaluesγ1 at several values of the M and K, where the values of other parameters are fixed as S=6, Pr=Sc=λT=λ=1, λC=0.5, Nb=0.4, Nt=0.6, Rd=0.7, m=0.1, and δ=0.2.
MKγ1First solutionSecond solutionThird solution0.10.11.54340-1.98503-2.79110.20.10.09248-0.43835-0.967500.30.12.10262-0.71688-1.598090.10.10.05824-0.40107-2.105210.10.150.63291-0.20561-1.903680.10.20.80136-0.57696-0.93832The influence of the material parameterK on the skin friction coefficient, couple stress coefficient, Nusselt number, and Sherwood number with variation of stretching/shrinking parameter λ are shown in Figures 2–5, where triple solutions are found for λ>λc and unique solution for <λc, while second and third solutions are merged into one another at the point λc and only one (first) solution remained continuous for variation of the λ. Furthermore, Figure 2 and Table 3 show that an increase in K with variation in λ increase the skin friction coefficient for the λ>0 and decrease forλ<0 in the first (stable) solution. Table 3 and Figure 3 indicate that an increase in Kalong the variation in λ, the couple stress rises for the λ>0 and reduces for the λ<0. Furthermore, Table 3 defines the numerical results related to the Figures 2 and 3.Table 3
Computation of skin friction and couple stress coefficient for several values ofλ and K.
f′′0g′0Kλ1st solution2nd solution3rd solution1st solution2nd solution3rd solution0.1-25.463761-3.74033-8.49748.856247172.091-15.61500.247399-8.04981-13.1350.433814197.318-20.92-5.1687-12.8786-17.877-9.33088211.336-26.9790.2-25.2882-3.85674-7.0388.10026756.6338-8.456600.239959-7.82914-11.6160.40804369.3817-14.7912-5.02881-12.4521-16.276-8.6607973.5664-21.0410.3-25.121937.42156800.232959-7.85944-10.1290.38477226.4923-8.22972-4.89749-12.2182-14.789-8.0594928.5981-15.258Figures4 and 5 along Table 4 indicate that an increment in the values of K with variation of λ, the Nusselt number, and the Sherwood number increase for λ>0 and reduce for λ<0in the first solution throughout the flow of the micropolar nanofluid.Table 4
Computation of local Nusselt number and the Sherwood number for several values ofλ and K.
−θ′0−ϕ′0Kλ1st solution2nd solution3rd solution1st solution2nd solution3rd solution0.1-22.400891.959991.775512.29511.539851.3001902.446032.031531.818382.350431.615161.3166522.48682.090111.869532.400121.678631.350930.2-22.397191.952671.822152.289991.522791.3475102.445992.036951.869282.350431.61521.3706722.489682.102951.925752.40421.68971.414740.3-22.393512.2848202.445932.031211.925882.350421.59861.4392622.514072.107771.983662.434911.688811.48869While the effect of magnetic parameterM on velocity profile f′η is shown in Figure 6, it can be seen that velocity decreases with increase in the parameterM. That occurs due to rising Lorentz forces transverse magnetic field that becomes a reason to decelerate the velocity of the fluid. Figure 7 indicates the influence of thermal convective parameter λT on velocity profilef′η. It is seen that the velocity of micropolar nanofluid declines as λT is increased. Figure 8 shows the influence of the suction parameter S on velocity profilef′η. It is examined that the increasing rate of the suction decreases the velocity and momentum boundary layer thickness in flow of the micropolar nanofluid. It is because of the flow approaches nearer to solid surface, so the suction takes place due to porosity of the medium which becomes a reason to decrease the velocity and its boundary layer thickness. Figure 9 illustrates the influence of the velocity slip parameter δ on the velocity profile f′η. It is observed that by rising value of the parameter δ, the velocity profile along the velocity boundary thickness decrease throughout the flow. Actually, the increasing value of δ develops a resistance forces at the surface that becomes the reason to decline the velocity profile, while the effect of parameter λ is found opposite to the influence of the parameter δ on velocity profile of micropolar nanofluid that is shown in Figure 10. Likewise, influence of parameter λ on angular velocity profile gη is shown in Figure 11. The figure shows that any increment in λ rises the velocity of the bulk particles that possesses microrotations during the flow. Figure 12 indicates the impact of the parameter S on angular velocity profile gη. It is seen that an increment in S decelerate the angular velocity related to bulk particles. Figure 13 specifies the effect of the parameter δ on the angular velocity profile gη. The result of this figure shows that the angular velocity of bulk particles reduces during the flow. The influence of various physical parameters on temperature profiles θη of micropolar nanofluid are shown by Figures 14–19. Figure 14 indicates that any increment in the parameter M rises the temperature of the micropolar nanofluid during flow. Against it, Figure 15 shows that the increasing rate of suction decreasing the thermal boundary thickness along the temperature profile θη of the micropolar nanofluid. Actually, this occurs due to rise in suction rate which raises the resistance force between the fluid layers and solid surface. Such sorts of results are observed in Figure 16, when the value of the Prandtl number (Pr) is increased. Actually, higher values of the Prandtl number decreases thermal diffusion of fluid particles which result in the decrease temperature of the fluid. Therefore, increasing the value of the Pr decreases the temperature of the all types of the fluids during the flow. An increase in thermal radiative parameter Rd raises temperature and boundary layer thickness during the flow of the micropolar nanofluid that is shown in Figure 17. Actually, the existence of the thermal radiation leads to enhancing the thermal boundary layer thicknesses. Practically, enlarging radiations yields a major enhancement in temperature of fluid. Figures 18 and 19 show the influence of rising the Brownian motion parameter Nb and thermophoresis parameter Nt on the temperature profilesθη,respectively. The results show that any increment in the mentioned parameters increases the temperature of the particular nanofluid and its boundary layer thicknesses. This takes place because of Nb and Nt increases the temperature gradient at the solid adjacent surface.Figures20–22 show the influence of the Brownian motion parameter Nb, Schmitt number, and thermophoresis parameter Nt on nanoparticle concentration profiles ϕη. Figure 20 shows that an increase in Nb deceases the concentration of the nanoparticles in the micropolar nanofluid throughout the flow. Figure 21 shows the influence of the Schmitt number (Sc) on the concentration profile ϕη. This figure shows that the concentration of the nanoparticles decreases with increase in the Sc. It is due to increasing the value of Sc which decreases the mass diffusion; the larger diffusion becomes a reason to decrease in the concentration field. Moreover, it also decreases the thickness of the concentration boundary layers, while influence of the Nt looks opposite to the influence of Sc on the concentration profile, where increasing rate of the Nt raises concentration of nanoparticles throughout the flow of micropolar nanofluid.


## 6. Conclusions

The effects of the various types of physical flow parameters on MHD mixed convection boundary layer flow, heat, and mass transfer of the micropolar nanofluid on exponentially vertical stretching/shrinking surface are studied here. The similarity solutions are achieved by applying similarity transformation from partial differential equations into the ordinary differential equation. The shooting method is implemented to obtain the numerical solution of the equations with variation of the used physical parameters. At various ranges of used parameters, triple solutions are obtained, so the stability is performed by help of the MATLAB through bvp4c. Effects of various used physical parameters such as the magnetic parameter, thermophoresis parameter, microrotation parameter, suction parameter, radiation parameter, Prandtl number, Brownian motion parameter, thermal slip parameter, and Schmidt number are examined, discussed, and also presented graphically. The main results of this study are described below:(1)
Skin friction coefficient, couple stress, Nusselt, and Sherwood numbers are increasing forλ>0 and decreasing for λ<0 with increase in material parameterK(2)
Triple solutions are found at the different ranges of the flow parameters, and from stability analysis, it is found that only first solution is stable as well as physically feasible(3)
The buoyancy parameter and high suction are the leading causes for existence of triple solutions in this problem(4)
The velocity field of the micropolar nanofluid decelerates with increase in the suction parameterS, magnetic parameterM, and velocity slip parameterδ, while increases for increase in thermal convection parameter λT and stretching parameter λ(5)
An increase in parametersM, Rd, Nt, and Nb raises the temperature profiles of the micropolar nanofluid, while it decrease with increase in Pr andS(6)
Any increment inNt rises the concentration boundary layer thickness and its profile, while an increase in Sc and Nb decreases the concentration boundary layer thickness and its profiles



---

*Source: 3147696-2022-04-11.xml*