Mixed convective-radiative dissipative magnetized micropolar nanofluid flow over a stretching surface in porous media with double stratification and chemical reaction effects : ADM-Padé computation

Title Mixed convective-radiative dissipative magnetized micropolar nanofluid flow over a stretching surface in porous media with double stratification and chemical reaction effects : ADM-Padé computation Authors Pattnaik, PK, Bhatti, MM, Mishra, SR, Abbas, MA and Beg, OA Publication title Journal of Mathematics Publisher Hindawi Type Article USIR URL This version is available at: http://usir.salford.ac.uk/id/eprint/62843/ Published Date 2022


Introduction
e study of boundary layer flow and heat transfer in a mixture of fluids flowing across a continuous extending surface has numerous industrial applications, including thermal polymeric coating flows, surface deposition and finishing, enrobing, and mechanical handling equipment such as conveyors [1,2]. Crane [3] was among the first researchers to study the viscous Newtonian boundary layer equations for two-dimensional continuous flow across a flat, incompressible elastic sheet traveling in a straight line. He established a good foundation for subsequent studies involving multiple thermal and electromagnetic effects. Mishra et al. [4] recently investigated heat and mass transfer in magnetized viscoelastic fluid with conducting of a viscoelastic (Walters B') fluid on a stretching surface. is study also included chemical reaction and non-uniform heat source effects, which are commonly encountered in materials processing.
Mass transfer is a common occurrence in a wide variety of industrial technologies including metallurgical processing, food manufacture, and electromagnetic (smart) coating dynamics. It may involve multiple species diffusing through a fluid and is often accompanied with heat transfer. Whereas Fick's law is usually deployed for species diffusion, for thermal conduction the Fourier law is required. Due to the diverse applications in chemical, nuclear, and mechanical process and materials fabrication systems (extrusion, coating, enrobing, etc.), mathematical models of the combined effect of mass transfer and chemical reaction both with and without electromagnetic fields have stimulated considerable interest among engineers and scientists [5,6]. In high temperature applications, radiative heat transfer is also encountered (in addition to thermal conduction and thermal convection) and usually an algebraic flux model is employed in simulations to compute the effect of radiative flux on transport phenomena. Numerous authors [7][8][9] investigated the radiative convective heat and mass transport in magnetohydrodynamic (MHD) flows. Pattnaik et al. [10] investigated the unsteady magnetohydrodynamics free convection flow, heat, and mass transfer with variable temperature and concentration from an inclined surface embedded in a saturated constant permeability porous medium. ey observed that the transport characteristics and substantially modified with heat source, inclination, porous drag force, and magnetic body force effects. Pattnaik et al. [11] also computed the influence of heat flux caused by a concentration and temperature gradient on an unsteady hydromagnetic flow from an infinite vertical porous plat, noting that the temperature and velocity are significantly altered by the magnetic field strength. In [12], it was noted that chemical reaction, heat source, and suction/injection at the wedge wall all have a significant effect on MHD boundary layer flow.
In many materials processing operations, non-Newtonian (rheological) fluids are utilized. ese exhibit characteristics which deviate significantly from viscous Newtonian fluid models. Many non-Newtonian models are available including viscoplastic (Bingham), viscoelastic (Maxwell, Oldroyd-B), shear thinning pseudoplastic (Carreau, Cross, etc.). However, these models merely modify the shear stress tensor and do not include microstructural effects. e particles suspended in liquids deployed in many materials processing operations can spin; i.e., these fluids have a microstructure. Such effects cannot be simulated with standard non-Newtonian models. erefore, a new branch of fluid dynamics, micromorphic fluid mechanics, was introduced by Eringen in the 1960s [13] which successfully simulates the angular velocity or microrotation of particles and molecules in non-Newtonian fluids. A simplification of the general micromorphic fluid model is the micropolar model [14] which restricts the degrees of freedom of the particles and also neglects axial contraction/ extension of microelements. is model and its modification to include thermal effects (thermo-micropolar model [15]) has proved immensely popular and has been deployed in many excellent studies since many different colloidal suspensions, smart fluids, and biological fluent media (e.g., blood, synovial fluid) can be accurately simulated by micropolar theory. A further advantage of micropolar fluid theory is that, by neglecting the spin (gyratory motion of microelements), the classical Navier-Stokes Newtonian model can be retrieved. Motivated by coating applications, many researchers have analyzed the external boundary layer flows of micropolar fluids from different geometries. ese investigations have usually deployed advanced numerical methods since the boundary value problems investigated are highly nonlinear. Bhargava and Takhar [16] used a variational finite element method to compute the stagnation micropolar heat transfer characteristics on a translating surface. A number of magnetohydrodynamics coating flows of micropolar fluids have also been studied in recent years. Ezzat et al. [17] studied MHD micropolar coating heat transfer with thermal relaxation effects using Laplace transform and numerical approaches. Many other interesting applications of magnetic coating flows have been reviewed by Douglas [18] including microstructural effects, for which the micropolar fluid model is very appropriate. Mishra et al. [19] have further examined the effects of double stratification, i.e., thermal and solutal stratification in magneto-convective flows of micropolar fluids with heat generation and buoyancy effects. All these studies have confirmed the substantial modification in thermal, species, and flow characteristics computed with micropolar models compared with either Newtonian or simpler non-Newtonian models.
Nowadays, conventional fluids are being progressively superseded in favor of nanofluids due to their superior thermal properties. Choi and Eastman [20] introduced nanofluids as colloidal suspensions of base fluids containing nanoparticles. Base fluids include oils, polymers, water, etc. Nanoparticles include metallics and non-metal oxides, carbides, and nitrates. Popular nanoscale models used for simulating nanofluid transport characteristics are the Tiwari-Das model and the Buongiorno two-component nanofluid model. e latter features a species concentration equation for nanoparticles and also emphasizes thermophoretic body force and Brownian motion dynamics effects. Wakif et al. [21] developed a generalized Buongiorno nanofluid model with the purpose of analyzing the effect of magnetic field on hydrodynamic stability. Nanofluids improve the thermal conductivity of lowconducting fluids such as ethylene glycol, engine oil, or water. Nanofluids therefore achieve enhanced heat conductivities which makes them more favourable for engineering and technological processes than conventional base fluids. Nanofluids have been employed extensively in a wide variety of technical applications, including heat exchangers, thermal engineering, cancer therapy, propulsion, tribology, and biomedicine. Additionally, there are several applications of magnetic fields in physics, medicine, and engineering. In magnetic materials processing, as noted earlier, external magnetic fields can be deployed to regulate flow characteristics and also manipulate the thermo-solutal dynamics in smart electrically conducting functional liquids such as electroconductive polymers. Flow patterns are significantly influenced by the intensity and organization of the magnetic field in the fluid. e magnetic field exerts an effect on the colloidal fluid particles, altering the fluid's concentration and heat transmission characteristics. A magnetic nanofluid is a nanofluid with electrically conducting properties imparted by nanoparticles which permit the liquid to be manipulated by an external magnetic field. Flows of magnetic nanofluids can be simulated by combining magnetohydrodynamics of viscous fluids with nanofluid models such as the Buongiorno model. Hayat et al. [22] computed the influence of magnetic field on the flow of Casson magnetic nanofluid over a stretching surface. Bhatti et al. [23] evaluated the impact of thermodiffusion and diffusion thermal parameters on a hyperbolic tangent magnetized nanofluid under the effect of Hall current. Khan et al. [24] investigated computationally the nonlinear thermal radiation flux effects on hydromagnetic Sisko nanofluid flow with convective boundary conditions. e above studies did not simultaneously consider microstructural and nanofluid characteristics in magnetohydrodynamic flow from a stretching surface in a porous medium. Motivated by developing a more comprehensive model for actual magnetic non-Newtonian nanofluids, the current work investigates the radiative-convective chemically reacting flow of magnetized nanofluid over a vertical stretching surface to a porous medium with thermal and solutal stratification, heat source, and viscous heating effects. e Rosseland flux model is used for thermal radiation, Eringen's micropolar fluid model is deployed for microstructural rheological features, and the Buongiorno nanoscale model is utilized to simulate nanoscale effects. e novelty of the present article is therefore the simultaneous consideration of multiple effectschemical reaction, thermal radiation, double stratification, micropolar non-Newtonian, and nanofluid characteristics in magnetohydrodynamic stretching sheet flow. To solve the novel boundary value problem, a semi-analytical/numerical approach is employed, i.e., Adomian decomposition method (ADM), is used in collaboration with Padé approximants. Graphs and tables are used to illustrate the impact of key parameters on heat, mass, and momentum characteristics. Extensive interpretation of the results is included from a physical viewpoint. Comparison of the ADM-Padé approach with numerical Runge-Kutta quadrature is also included. e simulations are relevant to magnetic nano-polymer coating fabrication processes.

Mathematical Model
We consider two-dimensional steady incompressible MHD chemically reacting micropolar nanofluid boundary layer coating flow over a vertical stretching surface to a Darcian isotropic porous medium under thermal and species buoyancy effects in an (x, y) coordinate system. ermal and solutal stratifications, as well as nonlinear thermal radiation, are all taken into account. e x-axis is plotted along the vertical surface and the y-axis is normal to this ( Figure 1). A uniform magnetic field of strength B 0 is applied transverse to the vertical surface, invoking a Lorentzian magnetic body force along the surface plane. e vertical sheet stretches with a linear velocity U w � ax where a > 0 represents the stretching rate constant. Hall current and Ohmic dissipation are ignored. However, a heat source is present and viscous dissipation is included. A first-order homogenous chemical reaction (destructive) is considered. Buongiorno's model is used for the nanofluid and Eringen's micropolar model for microstructural rheology characteristics.
Under the above assumptions, the governing boundary layer conservation equations for mass, momentum, angular momentum (microrotation), energy, and nanoparticle species (conservation) for the regime may be shown to take the following form: Journal of Mathematics 3 e relevant boundary conditions at the surface (wall) and in the free stream are prescribed as Here, equation (1) indicates the mass conservation equation i.e. the equation of continuity. Equation (2) is the linear momentum conservation equation. Equation (3) is the angular momentum conservation equation. Equations (4) and (5) are the energy and mass (solutal) conservations equations, respectively. e left-hand sides of each of equations (2)-(5) represent the convective derivative and the first term on the right-hand side terms denotes diffusion term for the corresponding profiles. In the Eringen micropolar model, there are two different viscosities, i.e., fluid dynamic viscosity and vortex viscosity. e third term of equation (2) is due to the resistivity offered by the magnetic field and the Darcian drag due to the porous medium fiber impedance. e last two terms represent the mixed convection effects dictated by both thermal and solutal buoyancy. Equation (3), i.e., the microrotation/angular momentum conservation, includes spin effects of the microelements. e third term on the right-hand side of the energy conservation equation (4) is the thermal radiation flux term (Rosseland), the fourth term indicates the coupling effect of viscous dissipation, and the last term is the additional heat source/sink. Further, the 2 nd term on the right-hand side of the species (concentration) conservation equation (5) is due to thermal diffusion and the last term represents a first-order chemical reaction effect. Following Rosseland's approximation [25] for nonlinear thermal radiation, the net radiation heat flux can be defined as follows: Implementing in the energy boundary layer equation (4) leads to Darcian isotropic porous medium saturated with reactive magnetic micropolar nanofluid Nanoparticles Vertical stretching surface

Journal of Mathematics
We introduce the following similarity transformations: Substitution of equations (10) in (2)-(9) generates the transformed dimensionless boundary layer equations for linear momentum, angular momentum, energy, and species, with associated boundary conditions, as follows: 1 In equations (11)-(14) the nondimensional parameters arising are Journal of Mathematics 5 All parameters are defined in the nomenclature.

Physical Quantities of Engineering Interest
e skin friction (C f ), wall couple stress coefficient (C s ), Nusselt number (Nu x ), and Sherwood number (Sh x ) are important quantities for materials processing operations. ese take the following form: Here,

Adomian Decomposition Method.
Adomian [26] developed a robust semi-analytical method to solve linear and nonlinear differential/algebraic equations named the Adomian decomposition method (ADM). is method also uses the concept of series solutions. ADM has been shown to be generally excellent for solving highly strong nonlinear scientific problems [27]. e principal objective of ADM is to achieve a unified series solution of linear/nonlinear ordinary differential equations (ODEs) or partial differential equations (PDEs) in initial and boundary value problems (IVPs and BVPs). In this method, the nonlinear operator is decomposed into a series of functions. e nonlinear coupled transformed ordinary differential equations (11)- (14) are solved semi-analytically using ADM. e details of the procedure of ADM for the present magnetic micropolar nanofluid stretching sheet problem are as follows: To rewrite equations (11)-(14), we introduced the op- (11)-(14) can be expressed as 6 Journal of Mathematics Hence, taking the inverse of equations (20)-(23) we get e corresponding initial conditions are e initial guess solutions as well as recursive expressions are as follows: Journal of Mathematics 7 Equations (35) and (36) can be simplified by truncating the binomial expansion of the denominator of equation (35). Assuming all the linear as well as nonlinear terms of the right-hand side of equations (33)-(36) in series form and taking m � 0, 1, 2, . . ., the recursive relations (33)-(36) can be solved. Evaluating the unknowns numerically using the MATLAB BVP solver and taking some specific values of pertinent parameters as , and s � 0.5, we get the approximate analytical solutions for stream function, microrotation, temperature, and nanoparticle concentration profiles as follows: All the functions and coefficients are used in truncated forms. To increase the accuracy in the solution of the present boundary value problem, we have adopted Padé approximants after applying ADM. [a, b] can be expressed as a rational function that contains two polynomials, say P N (η) as the numerator and Q M (η) a denominator of degrees N and M, respectively [28]. e notation [N/M] is used to denote the quotient. To obtain better results, we have combined the differential transform method and the diagonal Padé approximants [N/N] and used them in the process. Padé approximants are widely used (Baker and Morris, [28]) to accelerate the convergence of a given series. Arain et al. [29] also used Padé approximants in conjunction with a differential transform method (DTM) to compute the viscoelastic magnetic nanofluid flow between rotating disks. We have considered the function f(η) in power series form as

Padé Approximant. A Padé approximation of f(η) on
e notation c i , i � 0, 1, 2, . . ., represents the coefficients of the power series#, and f(η) is the associated function. e [L/M] Padé approximant is a rational fraction defined by Equating the coefficients of η L+1 , η L+2 , . . . , η L+M , we get Taking i < 0, c i � 0 and b 0 � 1, equation (44) can be written in matrix form as e values of b i can be found from the above expression and from equation (43) we obtain a 0 � c 0 , Combining equations (43)

Result and Discussion
e graphical results produced for variation of several emerging parameters will be discussed in this section. Earlier, Ramzan et al. [30] validated the present computer programs by including certain additional parameters, such as the porosity parameter in the momentum equation, Eckert number, and heat source parameter in the energy equation. e solution process is a strong advantage of the present ADM approach. We use approximate analytical ADM solutions to refine each profile and thereafter deploy the Padé approximant of various orders to refine each profile. Both results are, however, compared to the traditional numerical technique (Runge-Kutta shooting method). Figure 2 shows comparison plot for all the profiles with solution obtained by ADM-Padé, ADM, and numerical shooting, with some default assigned values of the physical parameters K � 0.1, M � 1, Kp � 0.2, λ � 5, δ � 3, Tr � 2, Pr � 2, S � 0.01, Ec � 0.5, Nb � 0.5, Nt � 0.5, Le � 3, t � 0.5, R d � 2, Kc � 1, m � 0.5, and s � 0.5. is indicates a good correlation between the convergence procedures of the methodologies employed in this study. e comparison of the stream function, velocity, microrotation, temperature, and concentration profiles is illustrated in Figures 3-7 for the numerical method, Adomian decomposition method (ADM), and the novel semi-Journal of Mathematics analytical method, ADM-Padé-approximant. ese graphs show the importance of the Padé-Approximant method as it converges more rapidly than ADM on its own when compared with numerical methods such as the 4 th -order Runge-Kutta method (shooting technique).
is encouraged us to study the behavior of different profiles of a micropolar-nanofluid flow by deploying Padé-approximant after using ADM. e study shows since ADM requires some computer programs to calculate the values of f ″ (0), ω ′ (0), θ ′ (0), and ϕ ′ (0) to find the initial solution, recurrence functions using these initial solutions give the final solutions as equations (37)-(40). After getting the recurrence solutions, the diagonal Padé Approximant of an order [3/3] is used for the graphical and tabular presentations. e study shows the diagonal Padé approximant gives accelerated convergent solutions. Figures 3-7 confirm these observations.
In Figure 8, the impacts of magnetic parameter (M), microgyration parameter (m), mixed convection parameter (λ), and temperature difference parameter (δ) on velocity parameter have been studied. Figure 8

(η), θ(η), ϕ(η)
Previous work, dotted profiles-Ramzan et al. [30] Previous work present  e Darcian body force is inversely proportional to Kp as observed in Figure 9(a), via the term. − (M + Kp − 1 )f ′ in the momentum equation (11). Higher Kp values imply greater permeability and therefore a decrease in solid matrix fibers in the porous medium. is reduces the Darcian drag force and accelerates the boundary layer flow leading to a decrease in momentum boundary layer thickness. Figure 8(b) shows that velocity is strongly damped with an increment in micropolar parameter, m. is parameter appears in the wall boundary condition (14) and relates to the nature of the microelements spin effect at the boundary. For the Newtonian case m � 0, velocity is a maximum. Figure 8(c) shows that for increasing values of mixed convection parameter (λ) that enhances the fluid velocity near the boundary for η ≤ 1 and afterward it decelerates the fluid. ermal buoyancy as simulated in the parameter λ � (Gr/Re 2 ) therefore has a complex relationship with the linear velocity field. Figure 8(d) presents the influence of the buoyancy ratio is parameter expresses the relative contribution of solutal (nanoparticle species) buoyancy force to thermal buoyancy force. For negative values of δ, these buoyancy forces oppose each other and for positive values they assist each other. It is apparent that for δ � − 3 there is a strong acceleration in the flow (linear profiles) whereas the opposite trend is computed for δ � +3 (parabolic profiles). ese trends are sustained at all values of transverse coordinate, δ.
Figures 9(a)-9(c) visualize the velocity evolution with different values of the micropolar material parameter (K), Prandtl number (Pr), and thermal stratification parameter (t). Here, non-zero values of material parameter K indicate the non-Newtonian characteristics of the fluid. K � (k/μ) expresses the ratio of micropolar vortex viscosity to Newtonian dynamic viscosity. Momentum boundary layer thickness is clearly reduced with greater micropolar parameter since the velocity is elevated strongly at all values of transverse coordinate. With elevation in Prandtl number, there is a weak enhancement in velocity. At Prandtl number of 1, both momentum and thermal diffusion rates are equal for the magnetic micropolar nanofluid. However, as Pr is increased, the thermal diffusion rate (and thermal conductivity) is decreased, and the momentum diffusion rate is much higher. is results in acceleration of the boundary layer flow and decrease in momentum boundary layer thickness. Higher values of Prandtl number (Pr � 10) correspond to actual magnetic coating materials [18]. Figure 9(c) shows that as a thermal stratification parameter (t � (d/b)) increases, the sheet temperature reduces which induces a boost in the velocity profiles.
is parameter features in the wall boundary condition, θ � 1 − t, in equation (15) and clearly exerts a strong influence on velocity distribution. Greater thermal stratification produces a decrease in momentum boundary layer thickness in the stretching flow regime. Figure 10 visualizes the impact of (a) M (b) m (c) λ (d) δ on microrotation profile, i.e., angular velocity, in the presence/absence of Kp. Maximum angular velocity is generally observed at the wall (sheet). With increment in magnetic parameter, M, there is a strong increase in microrotation, i.e., microelements spin faster in the magnetic micropolar nanofluid.
is is maintained at all values of transverse coordinate. Greater permeability Kp however generates a reduction in angular velocity, i.e., suppresses the gyratory motion of microelements. Figure 10(b) also demonstrates that the microrotation is elevated significantly with an increase in micropolar surface parameter (m). It is noteworthy that when the value of the microgyration parameter (m)  increases from m � 0, i.e., when microelements are unable to rotate (strong concentration) to m � 1 (turbulent boundary layer flows), the microrotation profile is enhanced. Figure 10(c) shows that angular velocity is an increasing function of mixed convection parameter (λ) in conjunction with the permeability parameter, Kp. Microrotation is therefore clearly elevated with greater mixed convection parameter values. Figure 10(d) shows that positive increment in buoyancy ratio parameter (δ) also produces an acceleration in the flow; i.e., the velocity is enhanced for the assisting thermal and species buoyancy case (δ > 0) whereas it is reduced for the opposing buoyancy case (δ < 0). It is also noteworthy that all profiles converge asymptotically to the free stream, confirming the imposition for an adequate infinity boundary condition in the ADM-Padé computations. Figures 11(a)-11(c) illustrate the impact of micropolar vortex viscosity material parameter (K), Prandtl number (Pr), and thermal stratification parameter (t) on angular velocity. Increasing values of K are observed in Figure 11(a) to strongly damp the angular velocity profile, indicating that with increased vortex viscosity the microelements are constrained in their gyratory motion. e opposite effect is computed with increasing Prandtl number, for which in Figure 11(b), there is a boost in the angular velocity (microrotation). e greater momentum diffusivity relative to thermal diffusivity associated with higher Prandtl numbers assists the spin of microelements with reverse trend as expected in case of increasing higher values of Prandtl number (Pr). Figure 11(c) shows that the microrotation is depleted with greater thermal stratification parameter (t � (d/b)). As with other plots, the effect is most profound near the sheet surface (wall). e decrease in wall sheet temperature associated with higher thermal stratification reduces linear momentum and angular momentum in the boundary layer regime. Angular velocity is therefore supressed as microelement rotation is subdued. Figure 12 portrays the impact of (a) thermal radiation, Rd (b) temperature ratio parameter, Tr, (c) Prandtl number, and (d) heat source, S on the temperature distribution. Figure 12(a) shows that the radiation from electromagnetic waves energizes the flow and elevates temperatures further from the sheet towards the free stream. However closer to the wall there is a decrement in temperatures with greater thermal radiation effect. With larger temperature ratio, as seen in Figure 12(b), temperature is suppressed near the wall whereas it is elevated significantly further from the wall. Increasing Prandtl number results in a decrease in temperatures computed further from the wall (Figure 12(c)). Increasing heat source parameter, S, clearly consistently boosts the temperatures at all locations from the wall to the free stream. is parameter represents hot spots in materials processing which can be used to provide a boost in thermal energy during fabrication operations [18]. ermal boundary layer thickness is elevated with greater heat source parameter. In all profiles, a temperature overshoot is also computed near the wall. Figure 13 illustrates the influence of (a) K and t (b) Ec (c) Nb (d) Nt on temperature profiles. Figure 13(a) shows that with increment in stratification parameter (t) temperatures are reduced and thermal boundary layer thickness is decreased. However, with increase in micropolar vortex parameter, K, temperatures are strongly enhanced, and thermal boundary layer thickness is increased. Increasing Eckert number, Ec, in Figure 13(b) produces a strong elevation in temperatures. e viscous dissipation is boosted owing to a greater conversion of kinetic energy in the flow to frictional heat with greater Eckert number. is energizes the flow and increases thermal boundary layer thickness. e temperature peak is also progressively displaced further from the wall (sheet) with greater Eckert number. Figure 13(c) shows that, with greater Brownian motion parameter (Nb), temperatures are initially suppressed near the wall whereas they are strongly elevated shortly thereafter, and this temperature enhancement is sustained towards the free stream. ermal boundary layer thickness is therefore enhanced generally with greater random motion of the nanoparticles (Brownian dynamics). Figure 13(d) shows that with increasing thermophoresis parameter (Nt) there is a sustained elevation in temperature at all locations from the wall to the free stream. e migration of nanoparticles under a temperature gradient is known as thermophoresis. is energizes the flow and increases thermal boundary layer thickness. Figure 14 illustrates the influence of (a) Lewis number, Le and thermal stratification parameter t (b) chemical reaction parameter, Kc (c) Brownian motion parameter Nb, and thermophoresis parameter Nt on nanoparticle concentration profiles. Figure 14(a) shows that there is a strong decrement in concentration magnitudes with greater Lewis number and also stratification parameter. e Lewis number is the ratio of thermal diffusivity with respect to mass diffusivity. With the increase in Lewis number, a slower diffusion process of the solutal concentration is encountered, resulting in the concentration boundary layer decreasing and a suppression in species (nanoparticle) diffusion rate. Figure 14(b) indicates that with an increase in chemical reaction parameters again there is a depletion in nanoparticle concentration values. e reaction is first-order homogenous and destructive, as simulated in the equation (14), via the term − PrLeKcϕ. is implies that nanoparticles react and are converted to another species. erefore, physically the concentration of original nanoparticles, i.e., ϕ(η), must decrease with more intense chemical reaction effect (greater Kc values) as verified in Figure 14  However increasing thermophoretic body force effect, Nt, encourages the migration of nanoparticles under the temperature gradient towards colder regions in the boundary layer and this elevates nanoparticle concentration and therefore produces a greater solutal (nanoparticle) concentration boundary layer thickness.
Finally, Tables 1-3 provide the results of the skin friction C f , wall couple stress (dimensionless microrotation gradient at the wall) C s , local Nusselt number Nu x , and local Sherwood number, Sh x ADM-Padé computations with respect to various contributing parameters. Table 1 shows the influence of micropolar material parameter, K, magnetic parameter, M, permeability parameter, Kp, mixed convection parameter, λ, and buoyancy ratio parameter, δ on skin friction and couple stress coefficients. e skin friction (linear shear stress function) and the couple stress coefficient both increase as the material parameter increases. With increasing M values (greater Lorentzian drag effect), the skin friction and couple stress coefficient are however suppressed. With increasing permeability effect (Kp), there is a decrease in Darcian porous impedance and the skin friction and wall couple stress are both enhanced. Elevation in mixed convection parameter λ depletes skin friction and also couple stress coefficient, whereas an increment in the buoyancy ratio parameter δ produces the opposite effect; i.e., it increases the skin friction at the sheet and also boosts the couple stress coefficient. e local Nusselt number Nu x quantifies the relative contribution of thermal convection to thermal conduction at the sheet surface and is tabulated for various parameters in Table 2. It is evident that increasing thermal radiation (Rd) and temperature ratio (Tr) parameters strongly enhances local Nusselt number; i.e., they produce a greater transfer of heat from the boundary layer to the wall. An increment in Eckert number (Ec) and the heat source parameter (S) however suppresses the local Nusselt number Nu x since these parameters heat the boundary layer and elevate temperatures.
is will suppress the heat transferred to the wall and will decrease Nu x magnitudes. Table 3 shows the variation of local Sherwood number Sh x for different values of Brownian motion (Nb), thermophoresis parameter (Nt), Lewis number (Le), and chemical reaction parameter (Kc). Increment in Brownian motion

Conclusions
Inspired by more robust simulations of magnetic non-Newtonian nanofluid materials processing operations, a mathematical model has been developed for electrically conducting micropolar nanofluid flow from a vertical stretching surface adjacent to a porous medium under a transverse magnetic field. Eringen's micropolar model has been deployed for non-Newtonian characteristics and the Buongiorno nanofluid model employed for nanoscale effects (thermophoresis and Brownian motion). e model also includes double stratification (thermal and solutal) and also chemical reaction effects, heat source, and viscous dissipation. Darcy's model is employed for the porous medium and a Rosseland diffusion flux approximation for nonlinear thermal radiation. e nonlinear governing partial differential conservation equations are rendered into nonlinear ordinary differential equations via relevant transformations. An innovative semi-numerical methodology combining the Adomian decomposition method (ADM) with Padé approximants and known as ADM-Padé is deployed to solve the emerging non-linear ordinary differential boundary value problem with appropriate wall and freestream conditions in MATLAB software. A detailed parametric study of the influence of key parameters on stream function, velocity, microrotation (angular velocity), temperature, and nanoparticle concentration profiles is conducted. Furthermore, skin friction coefficient, wall couple stress coefficient, Nusselt number, and Sherwood number are displayed in tables. e validation of both numerical techniques used, i.e., ADM and ADM-Padé, against a conventional numerical 4 th -order Runge Kutta method is also included and significant acceleration in convergence of solutions achieved with the ADM-Padé approach. e main findings of the current study may be summarized as follows: (i) e flow is decelerated with greater buoyancy ratio parameter whereas microrotation (angular velocity) is enhanced (ii) Increasing thermal and solutal stratification suppresses microrotation (iii) Concentration magnitudes are boosted with greater chemical reaction parameter and Lewis number (iv) Temperatures are significantly enhanced with radiative parameter (v) Increasing Brownian motion parameter depletes concentration values whereas increasing thermophoresis parameter enhances concentrations (vi) Increases in Lewis number and chemical reaction parameter both enhance the local Sherwood number (vii) An increment in Eckert number and the heat source parameter reduces local Nusselt number whereas they elevate temperatures and increase thermal boundary layer thickness (viii) Greater magnetic parameter (stronger Lorentzian drag effect) reduces the skin friction and couple stress coefficient (ix) With increasing permeability effect, there is a decrease in Darcian porous impedance and the skin friction and wall couple stress are both elevated (x) Elevation in mixed convection parameter λ reduces skin friction and also couple stress coefficient, whereas greater buoyancy ratio parameter increases the skin friction at the sheet and also boosts the couple stress coefficient (xi) ADM-Padé method achieves faster convergence and improved accuracy relative to conventional ADM or Runge-Kutta numerical approaches e present study finds applications in thermomagnetic coating processes involving nanomaterials with microstructural characteristics. Future studies may consider entropy generation minimization via the second law of thermodynamics.  Journal of Mathematics