MHD Nonlinear Mixed Convection Flow of Micropolar Nanofluid over Nonisothermal Sphere

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.+e 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. +e equations were calculated using method bvp4c from Matlab software for various quantities of main parameters. +e effects of various parameters on skin friction coefficient f′′(0), wall duo stress coefficient -G′(0), and convection mass transfer coefficient -Φ′(0) are analysed and presented through the graphs and tables. +e 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. +e 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.


Introduction
e 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. e 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. e impacts of physical parameters on fluid velocity, microrotation, temperature, and concentration are discussed and shown in graphs and tables

Mathematical Formulation
Consider a two-dimensional steady laminar boundary layer flow of a viscous micropolar nanofluid over a nonisothermal spherical body of unchanging radius c. Assume a nonisothermal surface temperature T w and the micropolar nanofluid is a constant concentration C w . e temperature and concentration at far, correspondingly, are T ∞ and C ∞ .
e 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.
Following Amanulla et al. [13], the differential equations governing the present problem are given as follows: Subject to the proper velocity, microrotation, and concentration boundary conditions, e proper thermal boundary conditions are (1) e fixed surface temperature (FST): B 0 magnetic field g gravity field T ∞ , C ∞ free stream condition Figure 1: Flow diagram and coordinate system.

Mathematical Problems in Engineering
(2) e fixed heat flux (FHF): 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. e case δ � 1 is meant for modeling the turbulent boundary layer flows. μ stands for the coefficient of fluid viscosity, ρ is the density, ] � μ/ρ is kinematic viscosity, c p stands for the specific heat, j � d 2 / �� � Gr √ is the microinertia per unit mass, φ � (μ + (κ/2))j is spin gradient viscosity, κ is the vortex viscosity, T stands for temperature, K is the thermal conductivity of the fluid, T w � (T ∞ + ΔT) and C w � (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.
where R 0 and R * 0 are constants and R t and R s are the constant coefficients of thermal and volumetric expansion, respectively. is relation will be nonlinear density temperature and concentration (NDTC) variation. c � (ρc) p /(ρc) f is the ratio between the effective heat capacity of the nanoparticle material and the heat capacity of the fluid, D ct represents the Soret-type diffusivity, T m is fluid mean temperature, T k is the thermal diffusion ratio, D T stands for the thermophoretic diffusion coefficient, σ is electrical conductivity, and B 0 represents magnetic field. Let r(x) � c sin(x/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: For thermal boundary conditions, we consider the nondimensional quantities as follows: (a) For fixed surface temperature, m is the wall temperature parameter, T w is the temperature at the wall, and c is the constant. e dimensionless temperature in FST can be selected as (b) For fixed heat flux FHF, the wall heat flux is taken to be changeable with radial distance r(x) from the symmetrical axis to the surface of the sphere, and we introduce the dimensionless temperature as where T w − T ∞ � (x P /c)(c/Gr 1/4 ). e continuity equation can be integrated by introducing a stream function ψ(x, y) as ur � zr ψ/zy and vr � − (zr ψ/zx), which satisfies the continuity equation (1). Hence, equations (2)-(5) can be written in the following nonlinear system of PDEs: (a) For FST case, 1 Pr

Mathematical Problems in Engineering 3
Subject to the boundary conditions, Subject to the boundary conditions, where prime represent differential with respect to ζ, ε � x/c stands for dimensionless stream-wise coordinate, bv � aGr 1/4 /c is velocity slip parameter, bh � δc + δ 2 /ε denotes microrotation slip parameter, bm � eGr 1/4 /c stands for solutal jump parameter, α � κ/μ the is material parameter, Gr � c 3 gR 0 (T w − 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 � c 3 gR * 0 (C w − 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 ) 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.
represents thermophoresis parameter, Ma � (σ/ρ)B 2 0 (c 2 /]Gr 1/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, 1 Pr 4 Mathematical Problems in Engineering Subject to the boundary conditions, e physical quantities of awareness in this problem are the limited skin friction coefficient C f , wall couple stress m w , the Nusselt number Nu, and the Sherwood Sh number, and they can be written as follows: where Using the nondimensional variables (9) and (10) and boundary condition (18), the limited skin friction coefficient C f , wall couple stress m w , the Nusselt number Nu, and the Sherwood number Sh are

Numerical Solution
e 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.
e 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].

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.

Velocity and Microrotation
Profiles. e dimensionless velocity profile graphs of f ′ (ζ), 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       , λ), 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 3 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.

Temperature and Concentration Profiles.
e 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),       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. Figures 22 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.
Figures 28-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. e presentation of grid-independence test is shown in Tables 1 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.
Tables 3 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]. e outcomes show that the numerical used in the information is truthful and highly accurate.
Tables 5 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). is 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).

Conclusions
is 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. e boundary layer equations of the flow problem are reduced into a duo of high-order nonlinear ordinary differential equations with the similarity transformation.
e obtained differential equations are solved numerically with bvp4c from Matlab software. Numerical results are acquired for various main parameters of the flow problem. e results of governing parameters are reported by means of figures and tables. e main findings are (1) An increase in the values of m, 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) e existence of Soret number So and nonisothermal parameter of FST m 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) e 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. is result happens due to improved density of the fluid (6) e upsurge in values of Ma from 1 to 2 results in the improved density of the fluid which is associated