MHD Powell–Eyring dusty nanofluid flow due to stretching surface with heat flux boundary condition

Introduction Magnetohydrodynamics (MHD) tends to describe the behavior of electrically conducting fluids like liquid metals, plasmas, and astrophysical systems as well as studying the magnetic field and velocity in these fluids. Characteristics of heat transfer and the effects of Hall on magnetohydrodynamic (MHD) flow in a rotating channel are analyzed by Ghosh et al. [1]. With regard to the MHD power generators, Carabineanu [2] introduced mathematical theory in the simplified form, disregarding the effect of Abstract

Lorentz force contrary to the direction of fluid flow which takes the stationary flow of plasma. Some flows of MHD in a porous medium of the second-degree fluid with Darcy law for modeling the flow are investigated by Hayat et al. [3]. Alfv´en [4,5] put the magnetohydrodynamics (MHD) equations requisite basic and he recognized the importance of the generated magnetic field as well as the electric currents loaded by plasma. Thermal radiation impact on the flow of magnetohydrodynamic (MHD) over a vertical plate with considering convective boundary condition into the model of flow has been discussed by Etwire and Seini [6]. Flow problem of mixed convection boundary layer and mass transfer over a vertical plate in the presence of a magnetic field and constant heat flux through a porous medium is illustrated by Makinde [7]. In the presence of the magnetic field, Afrand et al. [8] presented a numerical approach of natural convection flow in three-dimensional inside a vertical annulus loaded by gallium with taking lower and upper parts of the annulus in the account. The flow of MHD of Maxwell fluid caused by a moving flat plate subject to second-order slip is investigated by Liu and Guo [9]. Studying magnetohydrodynamic (MHD) flow of incompressible fluid with heat transfer effects between two circular discs has been presented by Ayaz et al. [10]. Recently, Patel [11] utilized the technique of Laplace transform to examine the flow of MHD of Casson fluid through porous medium over an oscillating vertical plate in the presence of Hall current, thermal radiation, and heat generation.
Nanofluid is a distinctive category of heat transfer fluids and can be gained by uniform dispersion of nanoparticles as well as stabilized suspension of these particles in base fluids like water, ethylene glycol, kerosene, oil, and bio-fluids. Nanofluid which has additional thermal properties achieves the largest and highest potential thermal properties utilizing particles at very small concentrations suitable. Furthermore, nanofluids are allocated by thermal conductivity. Nanoparticles are defined as metal or non-metallic particles with a size below than 100 nm that move a random movement known as Brownian motion and improve energy transfer in nanofluids. Choi [12] was the first to introduce the concept of nanofluid as a mixture of disperse particles (nanoparticles) and liquid. Among the many of nanofluid applications, we mention transportation, biomedical, power generation, and microelectronics. An important application of nanofluid is to improve the coefficient of heat transfer of the heat transfer fluids. Xuan and Li [13] offered the major reasons for improving the instrument of heat transfer of fluids by the concept of nanofluids. The research which deals with the methods of preparation of nanofluids, its applications, and its stability, as well as its thermophysical properties was presented by Yu and Xie [14]. MHD boundary layer flow with natural convection of nanofluid over a moving vertical flat plate in the presence of radiative heat flux and a magnetic field was described by Das and Jana [15]. Papers on nanofluid flow can be found by Mahdy [16] and Hady [17,18]. Ibrahim and Makinde [19] investigated numerically the problem of heat transfer and MHD stagnation point flow with convective heating and velocity slip impacts of power-law nanofluid towards a stretching sheet. Recently, Darcy-Forchheimer model through a porous medium is applied for viscous nanofluid flow because of a curved stretchable surface by Saif et al. [20]. They noted that increasing the parameter of Brownian motion tends to reduce local Nusselt number while enhances from local Sherwood number. The method of lattice Boltzmann has been utilized by Ma et al. [21] to examine MHD natural convection of nanofluids within U-shaped enclosures with a baffle.
Powell-Eyring fluid model is regarded as an important type of non-Newtonian fluid models for its essential uses in diverse geophysical, natural, and industrial processes. Its importance is due to the fact that most fluids used in the industry have the same nature as non-Newtonian fluids. There are many patterns that have been developed to describe the behavior of non-Newtonian fluid flow; this is because of the nature of complex non-Newtonian fluids. Powell-Eyring fluid model is characterized by following the same behavior as Newtonian fluid in the case of rates of low and high shear. It is also derived from the kinetic theory of liquid but not from the empirical relation as it was before. In 1944, Powell and Eyring suggested a full and adequate mathematical model is Eyring-Powell fluid model [22]. Boundary layer flow of Powell-Eyring fluid is described across a stretching surface by many researchers in their work. This due to the importance of a stretching surface in various industrial applications including electronic chips, fiber yarn, polymer industries, and glass blowing. Here we will display a pool of previous works related to Powell-Eyring fluid. Javed et al. [23] presented Powell-Eyring fluid boundary layer flow over a stretching sheet. The numerical approach of MHD boundary layer flow of non-Newtonian Eyring-Powell fluid towards a stretching surface has been presented by Akbar et al. [24]. Hayate et al. [25] examined the problem of steady flow and heat transfer of non-Newtonian Eyring-Powell fluid in the presence of convective boundary conditions over a moving surface. The peristaltic flow with mass and heat transfer of Powell-Eyring fluid over a curved channel has been illustrated by Hina et al. [26]. Ghadikolaei et al. [27] considered the effect of Joule heating, thermal radiation, and Lorentz force when studying the unsteady problem of Eyring-Powell fluid flow by Akbari-Ganji Method (AGM). Powell-Eyring nanofluid flow due to gyrotactic microorganisms and magnetic field effect towards a stretched surface has been examined by Naseem et al. [28]. Panigrahi et al. [29] applied the method of Homotopy analysis to solve the mixed convective problem of Powell-Eyring fluid flow under the effect of Dufour number and Soret number through a nonlinear stretching surface. They found that the profiles of temperature and concentration improve with increasing Dufour number and Soret number, respectively. The problem of MHD boundary layer flow of a non-Newtonian nanofluid enforcing the model of Powell-Eyring past an impermeable nonlinear stretching sheet with variable thickness has explored by Hayat et al. [30]. They observed that an increase in the thermophoresis parameter leads to an increment of both concentration and temperature profiles. Recently, Kumar et al. [31] presented an incompressible laminar flow and the characteristics of heat transfer on non-Newtonian Powell-Eyring fluid at a shrinking wedge in the presence of the effects of the magnetic field, irregular heat source /sink, and radiative heat flow.
Dusty fluid flow phenomenon subsists in the flow of fluid compressible and incompressible such as gas and liquid, respectively, which containing solid particles in the size of the micrometer. Model of compressible gas flow through porous media in the presence of dust particles has been presented by Hamdan and Barron [32]. The flow of the boundary layer over a stretching surface of the dusty fluid under an applied magnetic field is studied numerically by Jalil et al. [33]. Ramesh and Gireesha [34] analyzed the influence of radiation on the flow of boundary layer of an incompressible dusty fluid past a stretching sheet. Dust fluid flow modeling is important because of its application's atmospheric fallout, paint spray, guided missiles, and dust collection. Furthermore, several works related to the flow of dust fluids were presented by a number of authors [35][36][37][38][39][40][41]. The object of our paper is to present an analysis of the problem of MHD flow of non-Newtonian Powell-Eyring fluid in the presence of nanoparticles and dust particles towards a stretching vertical plate with heat flux condition. The Runge-Kutta method is applied in the fourth-order to obtain numerical solutions of ordinary differential equations subject to appropriate boundary conditions. Perform a comparison to check the integrity and accuracy of the method used. The impacts of material fluid parameters, nanofluid parameters, and particle-phase parameters, Prandtl number, Grashof number, magnetic parameter (Hartmann number), and Eckert number on local Nusselt number and skin friction coefficient are estimated in tables.

Governing equations
Steady MHD boundary layer flow of an incompressible Powell-Eyring nanofluid over a stretching surface in the presence of dust particles is investigated. The flow is in twodimensional and the stretching plate coincides in Cartesian coordinates ( x, y) . So that one of these coordinates x-axis is in the direction of flow and the other y-axis is in the opposite direction of flow. Figure 1 shows the geometric representation of our problem. The heat flux condition and the influences of thermophoresis and Brownian motion are taken into nanofluid model. As well as the thermophoresis effect, the normal flux of nanoparticles at the boundary is zero. The stretching velocity is denoted by U w and at the plate surface the temperature is given by T w . Away from the stretching plate, the temperature and nanoparticle volume fraction are represented, respectively, by T ∞ and C ∞ . Using a transverse magnetic field and uniform B o , the fluid becomes electrically conductive. The induced magnetic field has very little effect because the magnetic Reynolds number is also small without applying the voltage. After the previous view, the equations of continuity, momentum, nanoparticle concentration, energy, take the following form: For the fluid phase: For the dust phase: In the above equations, there are fluid velocity components û,v , particle-phase velocity components û p ,v p , kinematic viscosity ν f , density of dust particles and fluid ρ p and ρ , material parameters of Powell-Eyring model β f and c , volume thermal expansion coefficient of the fluid β , acceleration of gravity vector g , strength of the magnetic field B 0 , fluid electrical conductivity σ , dust particle thermal relaxation time and dust particle velocity relaxation time τ T , and τ m , dust particles temperature and fluid temperature T p and T , fluid concentration C , thermal diffusivity of the base fluid α = k/(ρc p ) f , nanofluid heat capacity ratio τ = The pertinent boundary conditions as follow: By using the following conversions, an appropriate form of the equations to be processed is obtained: For the fluid phase: For the dust phase: where primes refer to differentiation for η , Prandtl number is denoted by Pr = ν α , parameters of material fluid are given as � = a 3x2 2νc 2 and ω = 1 ρ f β f cν , magnetic parameter (Hartmann number) is represented in M = σ B 2 0 aρ f , Grashof number is presented as , Schmidt number is symbolized by Sc = ν D B , Buoyancy ratio parameter is written as Nr = The dimensionless form of skin friction and Nusselt number is represented as where the local Reynolds number is defined as Re x = U 0 x ν .

Results and discussion
The systems of differential Eqs. (10)- (14) and (15)     determine the accuracy of the numerical method or the estimated numerical results. This comparison actually exists in Table 1 while Tables 2 and 3 present the values of local Nusselt number and skin friction coefficient under the influence of various values of the emerging physical parameters. In this study, the results obtained are valid enough after the comparison. The effect of Prandtl number, material fluid parameters, nanofluid parameters, and particle-phase parameters, etc., on the field of velocity, temperature, and concentration is examined, and the results of these tests are illustrated through figures. M as shown in Fig. 4. Lorentz force is an anti-flow force produced by the high values of M and therefore takes place decreases in motion. It is also known that this force improves the thermal boundary layer thickness any increase in temperature. Figure 5 (a-b) is plotted to display the influence of dust parameters on the velocity and temperature profiles. It is observed from Fig. 5 (a) that an enhancement takes place in particle phase velocity F ′ due to improvement in values of fluid-particle interaction parameter (α d = 0.01, 0.02, 0.03, 0.04) . The particle phase reduces the velocity of the fluid till it reaches the same fluid velocity; this is why the velocity of the particle phase increases, and the velocity of the fluid phase decreases. The behavior of θ and θ p when studying the mass concentration of the dust particles D ρ = 0, 1, 10, 15 is presented in Fig. 5 (b). It is detected that the temperatures deteriorate with an increment in D ρ . This behavior can be explained by the fact that the thermal conductivity of the particle phase is improved as a result of the increased dust particle's mass concentration, meaning that the rate of heat transfer increases, causing a shrinking in temperature distribution. Utilizing Figs. 6 and 7, the temperature profiles for particle phase and fluid phase are depicted under the effect of Powell-Eyring fluid parameter (ω = 0.0, 0.4, 0.9, 1.2) , Buoyancy ratio parameter (Nr = 0.1, 0.3, 0.5, 0.6) , Prandtl number (Pr = 0.7, 0.8, 0.9, 1) , and Eckert number (Ec = 0.08, 0.2, 0.4, 0.6) . Figure 6 (a) describes that for the rising values of ω , the temperature profiles θ and θ p diminish. Fig. 6 (b) displays the temperature profiles θ and θ p for distinct values of Nr . The temperatures ameliorate with increased Nr as shown in Fig. 6 (b). Characteristics of Pr on the temperatures θ and θ p are represented in Fig. 7 (a). It is clear from this figure that the temperatures diminish. The decrease in temperature distribution is due to that the increase in Pr leads to reduce the thickness of the thermal boundary layer. It is noted in Fig. 7 (b) that temperatures θ and θ p positively affected by the increase in Ec.
The impact of several parameters like Brownian motion parameter . It is clear that the concentration profile increases as Nt increment and it decreases with increasing NB as confirmed in Fig. 8 (a-b). The behavior of Nt can be traced back to increasing the distribution of nanoparticles near the surface due to the transition of these particles from the hot fluid system to the surface, and this occurs when the surface is cold, i.e. Nt > 0 . Progress in the concentration profile appears by Ec in Fig. 9 (a). Figure 9 (b) reveals the outcome of Sc effect on the concentration profile. The concentration profile deteriorates with rising values of Sc . Figure 10 illustrates that the magnitude of the velocity with Powell-Eyring nanofluid is greater when compared with Newtonian nanofluid, while the temperature distribution of Newtonian is higher than those of Powell-Eyring.

Conclusions
In our work entitled numerical investigation of MHD Powell-Eyring dusty nanofluid flow towards a stretching vertical surface with heat flux boundary condition, we computed the values of Nusselt number and skin friction coefficient for various values of the governing parameters. Further inserted the numerical results of velocity, temperature, and concentration distributions for both particle and fluid phases graphically. The significant main findings of the present work are filed as follows: 1. The concentration of nanoparticle improves with thermophoresis parameter and Eckert number but deteriorates with Brownian motion parameter and Schmidt number. 2. The profiles of velocity for particle and fluid phases enhance with the increase in Powell-Eyring fluid parameter and Grashof number, while they reduce with height in magnetic parameter values. The velocity of the particle phase increases with increment fluid-particle interaction parameter. 3. An enhancement in Prandtl number, Powell-Eyring fluid parameter, and mass concentration of the dust particles reduces the temperature of both the particle and fluid phases. 4. Rising the values of both the Buoyancy ratio parameter and Eckert number tends to optimize the temperature distribution in two-particle and fluid phases.