Abstract

The aim of this numerical research is to study the stagnation point flow of the electrical magnetohydrodynamic micropolar nanofluid with slip conditions past a stretching sheet. The phenomenon of linear thermal radiation, Ohmic and internal heating, has also been considered in the energy equation. The modelled PDEs are converted into ODEs via similarity transformation, and converted ODEs are tackled via the shooting technique. The features of assorted parameters on the axial and angular velocities and energy and concentration fields are sketched. The numerical values of the Sherwood and Nusselt numbers have been computed numerically and displayed in the form of tables. Our analysis shows that the heat transfer rate is decreased as the thermal slip parameter and the diffusion slip parameter are enhanced. The present study illustrates that the energy and concentration distribution are decreased with each of the mass free convection parameter, stagnation parameter, and thermal free convection parameter.

1. Introduction

The phenomenon of transfer of heat, which has significant application in many engineering and industry disciplines, is affected positively by the implementation of an appropriate magnetic field. The investigation of magnetohydrodynamic flow past a heated surface has gained significant attention due to its vast applications in engineering problems, i.e., magnetohydrodynamic power generators, petroleum industries, and crystal growth. Swedish scientist Alfven [1] was the first to introduce the magnetohydrodynamic fluid flow. He won the Nobel Prize in Physics for his work on MHD in 1970. He described the class of magnetohydrodynamic waves, which are now known as Alfven waves. Zheng et al. [2] reported the magnetohydrodynamic 2-D (dimensional) flow past a porous shrinking surface with slip conditions with conclusion that an acclivity in the shrinking parameter enhances the thermal boundary layer. By considering the slip effects in a porous medium, Ullah et al. [3] explored the magnetohydrodynamic (MHD) Casson fluid and noticed that boosting the unsteadiness parameter enhances the wall shear stress. Rahbari et al. [4] examined the magnetohydrodynamic Maxwell fluid flowing through parallel plates and determined that increase in the Deborah number increases the velocity. By investigating heat transfer in the magnetohydrodynamic flow past a radially shrinking/stretching sheet, Soid et al. [5] concluded that dual solutions exist only in case of suction and for small values of magnetic parameters. Decline in surface drag due to increment in the squeezed flow parameter in Carreau fluid with thermal radiation and magnetohydrodynamic effect past a sensor surface was reported by Atif et al. [6]. Transverse magnetohydrodynamic effect on nonlinear stretching sheet was ascertained by Ramana et al. [7]. They used modified Fourier flux law and concluded that the relaxation and retardation time have opposite effects on the thermal profile. For further studies, see [8, 9].

Stagnation point flow is one of the fields in which scientists and engineers show keen interest. Some of the recent studies include the following: For MHD viscoelastic nanofluid, the dual solution of stagnation point past a porous stretching surface with radiation effect was reported by Juosh [10]. It was found that an acclivity in the Deborah number contributes to upsurge in the drag coefficient. Impact of MHD on stagnation point flow of a nanofluid with nonuniform thermal reservoir was analyzed by Rashid et al. [11]. Pal [12] put light on magnetohydrodynamic stagnation point flow with suction effect and reported that the Sherwood number was decreased as the Lewis number increased. Bioconvective stagnation point of Maxwell nanofluid flow past a convectively heated surface was reported by Abbasi et al. [13]. They observed that energy, concentration, and density profiles were higher for nonconvective surfaces than in convective heated surfaces. Lund et al. [14] performed the stability analysis and reported the dual solution of MHD stagnation point of Casson fluid. Their main observation was that the sign of the smallest eigenvalues shows that the first solution was stable. Effect of solar radiation on MHD stagnation point nanofluid flow was discussed by Ghasemia and Hatami [15] with a key finding that the energy profile is hiked as the Biot number is increased rapidly. The phenomena of thermal radiation have much significance in the transfer of heat and were discussed by many authors in the literature [1619].

Micropolar fluids can be characterized as fluids which exhibit the micro-rotational effects and micro-rotational inertia. Analysis of the micropolar fluids has been an active field of interest for many researchers. This class of fluids possesses certain simplicity and elegance in their mathematical formulation which should appeal to mathematicians. The micropolar fluids can support couple stress and body couples only. Physically, they may represent adequately the fluids consisting of dipole elements. Certain anisotropic fluids e.g., liquid crystals which are made up of dumbbell molecules are of this type. In fact, animal blood happens to fall in this category. Other polymeric fluids and fluids containing minute amounts of additives may be represented by the mathematical model underlying micropolar fluids. Eringen [20, 21], through his pioneering work, invited the attention of the researchers’ community in this interesting area of fluid dynamics. Sui et al. [22] investigated the nonlinear constitutive diffusion model in the micropolar fluid with the main finding that both the velocity and energy profiles are increased as the power exponent is decreased from 1. Heat transfer of the free convective micropolar fluid with heat source past a shrinking sheet was noticed by Mishra et al. [23]. They observed that the fluid motion is declined as the heat generation coefficient is upsurged. Atif et al. [24] analyzed the bioconvective magnetohydrodynamic micropolar nanofluid with stratification and reported that the density distribution decreases as the density stratification and mixed number parameter are hiked. Micropolar nanofluid flow with nonlinear convection and multiple slip effects was examined by Zemedu and Ibrahim [25] with concluding remarks that boosting the solutal nonlinear convection parameter causes an increase in the velocity.

The heat transfer in base fluids like mineral oils, water, and ethylene glycol is not as much effective as in nanofluids [2631]. Nanofluids have the ability to improve the heat transfer properties. Their ability to move through capillaries and microchannels without making any blockage in flow makes them unique. By considering the induced magnetic field, Atif et al. [32] investigated the magnetohydrodynamic micropolar Carreau nanofluid and found that the angular velocity is increased rapidly as the magnetic Prandtl number increases. Three-dimensional Eyring–Powell nanofluid with Arrhenius energy was reported by Taseer et al. [33]. Khan [34] reported that nanoparticle dispersion reduces the Nusselt number in a partially heated vertical annulus. For a solutal-dominated regime, both Nusselt and Sherwood numbers declined for micropolar nanofluid as reported by Manaa et al. [35].

In recent years, researchers have paid serious attention to electrical magnetohydrodynamics. Electrical magnetohydrodynamic stagnation point nanofluid with mixed convection and slip boundary over a stretching surface was scrutinized by Hsiao [36]. A major conclusion was that an acclivity in either the electrical or the magnetic parameter led to an upsurge in the temperature profile. Literature review indicates that the EMHD stagnation point micropolar nanofluid with mixed convection and slip boundary has not been investigated yet. In the present article, linear thermal radiation, Joule’s heating, and heat source have also been incorporated in the energy equation. In this study, four aspects have been focused. First, the heat and mass transfer of micropolar nanofluid are addressed. Second, the impact of thermal radiation and electrical magnetohydrodynamics on different profiles is examined. Third, the stagnation point flow is analyzed. Fourth, the analysis of the mixed convection and slip boundary conditions is performed. The arising ordinary differential equations for the problem are tackled through the shooting method. The influence of all the prominent parameters is examined numerically and displayed graphically.

2. Mathematical Model

An incompressible, 2D, mixed convection micropolar nanofluid flow over a stretching sheet with slip effects has been analyzed. By using Ohm’s law and Maxwell’s equations, the continuity equation, linear and angular momentum, and fluid energy and concentration equations have been formulated. Joule’s heating, thermal radiation, and heat source effects have also been considered in the energy equation. , and denotes the surface concentration, ambient concentration surface temperature, and ambient temperature, respectively. The flow is assumed along the which is considered to be in the upward direction, whereas the is perpendicular to the sheet. A uniform magnetic field has been implemented towards the as illustrated in Figure 1. It is also assumed that the magnetic Reynolds number is very small due to which induced magnetic number is ignored.

In the light of the above assumption, the governing equations of the modelled problem are as follows [36, 37]:

The related BCs are as follows:

In equation (3), the Rosseland radiative heat flux is given by . The spin gradient viscosity is given by , where and represents the microinertia density and vortex viscosity, respectively. For nondimensionlization, the following transformation [36] has been considered:

Continuity equation (1) is satisfied automatically, and equations (2)–(5) yield the following:

The associated boundary conditions in the dimensionless form are as follows:

Here, denotes the Prandtl number, represents the stagnation parameter, denotes the Eckert number, denotes the Brownian motion parameter, represents the micropolar parameter, represents the thermal free convection parameter, represents the magnetic number, represents the electric field parameter, denotes the thermophoresis parameter, denotes the slip parameter, denotes the Schmidt number, denotes the mass free convection parameter, denotes the thermal radiation parameter, denotes the heat generation coefficient, represents the shear stress parameter, represents the temperature slip parameter, and represents the diffusion slip parameter, where and are the slip parameters associated with the reference temperature and concentration, respectively.

3. Quantities of Interest

The dimensionless Nusselt and dimensionless Sherwood numbers are the most concerning quantities in engineering and industries.

In the nondimensional form, Nusselt and Sherwood numbers are given bywhere .

4. Solution Methodology

The system of ODEs (8)–(11) along with BCs (12) is tackled numerically by the shooting technique. Now, we introduce , , , , , , , , and as follows:

The corresponding boundary conditions are

4.1. Code Validation

For the verification of the correctness of the code, the results of the Nusselt and Sherwood numbers which were presented by Khan and Pop [38] and Hsiao [36] are successfully reproduced. Our simulations have a satisfactory agreement with the already published results of Khan and Pop [38] and Hsiao [36] in the literature which can be seen in Table 1.

5. Results and Discussion

Table 2 is displayed to view the effect of the sundry parameters on the dimensionless Nusselt number and Sherwood number. It is observed that a boost in each of the electric parameter , slip parameter , micropolar parameter , thermal radiation parameter , the stagnation parameter , thermal free convection parameter , diffusion slip parameter , and mass free convection parameter , causes an increase in Nusselt number, whereas it decreases for a boost in each of the heat generation coefficient , shear stress parameter , magnetic parameter , and the thermal slip parameter . The Sherwood number is hiked as each of the material parameter , magnetic parameter , thermal free convection parameter , velocity slip parameter , electric parameter , b mass free convection parameter , stagnation parameter , and heat generation coefficient , and thermal radiation parameter is boosted. However, it diminishes as the diffusion slip parameter , thermal slip parameter , and shear stress parameter are increased.

Figures 214 are sketched to study the variations occurring due to dimensionless parameters in temperature distribution . For all graphs of the temperature dynamics, the values of , , , and . Figure 2 is prepared to visualize the fluctuation in the temperature distribution in response to the variation in the magnetic effect . The temperature of the fluid is increased as increases. This supports the general behaviour of the implementation of . Resistance to the flow of the fluid is increased as increases due to which is enhanced. The temperature distribution is diminished as the stagnation parameter is enhanced. This effect is evident from Figures 3 and 4 demonstrated to view the effect of on temperature distribution . These graphs indicate that upsurge in the Prandtl number causes depreciation in the thermal profile. There is decline in the thermal conductivity of the fluid due to which is reduced. Fluctuation due to electrical parameter in the temperature is divulged in Figure 5. An increment in results in an enhancement in thermal profile. As the Lorentz force is associated with the magnetic and electric field, it leads to an increment in the resistance causing the energy distribution to enhance. The influence of the thermal radiation parameter on temperature distribution is chalked out in Figure 6. These graphs reflect that an enhancement in increases the energy profile. To visualize the behaviour of the thermophoresis parameter on temperature, Figure 7 is sketched, which shows that there is an increment in as is hiked. Physically, in thermophoresis, the particles apply force on the other particles due to which particles from the hotter region move towards the colder region. Larger values of denotes more application of the force on the other particles and as a result, more fluid moves from the higher temperature region to the colder region. Figure 8 shows that is increased as the Brownian motion parameter increases. Physically, the Brownian motion heats up the fluid and also aggravates the particles away from the fluid regime and therefore a decrement is seen in concentration profile. Figure 9 is illustrated to view the effect of heat generation coefficient on , which shows that is hiked for escalating values of . The viscous dissipation effect which is represented by the Eckert number on energy field is analyzed in Figure 10. It is a number that represents the relation between the kinetic energy and the change in enthalpy. It is noticed that gradually boosting leads to an increase in . Influence of the dimensionless slip parameter on is presented in Figure 11, and it shows that an upsurge in encourages the energy distribution to decline. The influence of the slip parameter which is associated with temperature on the temperature field is chalked out in Figure 12. The energy profile is found to be increasing as slip parameter goes up. The variation in the thermal profile due to the thermal free convection parameter is shown in Figure 13. The energy distribution declined as the thermal free convection parameter is hiked. Figure 14 depicts that the energy profile is decreases as is boosted.

Figures 1523 have been outlined to study the fluctuations in the concentration field due to variation in the governing parameters. For all the graphical presentations of , we have considered . Figure 15 is given to study the impact of stagnation parameter on . The concentration distribution is diminished as the stagnation point parameter is enhanced. Figure 16 represents the graph of for growing values of . From these curves, it is clear that increasing electric parameter diminishes . Figure 17 presents the role of Schmidt number in the variation of . The concentration field is diminished as is rapidly increased. The effect of on dimensionless concentration field is reported in Figure 18. These graphs present that concentration field is enhanced as is gradually increased. Figure 19 is presented to view the fluctuation in dimensionless caused by the increase in . An increase in dimensionless parameter causes a reduction in . The impact of the variation of the dimensionless velocity slip parameter on the dimensionless is shown in Figure 20. These curves indicate that with an increment in the velocity slip parameter , the dimensionless declined. Figure 21 depicts the graphs of the concentration profile for various values of the dimensionless diffusion slip parameter . From these curves, it is noticed that an enhancement in the dimensionless diffusion slip parameter causes a decrease in . The fluctuation in the dimensionless concentration distribution due to the thermal free convection parameter is shown in Figure 22. The concentration field is reduced as thermal free convection parameter is increased. Figure 23 depicts that is decreased as the mass free convection parameter is increased.

6. Concluding Remarks

In this study, two-dimensional free convection electrical magnetohydrodynamic micropolar nanofluid is analyzed. Some of the key observations are as follows:(i)The energy field declined with an acclivity in the stagnation parameter , slip parameter , thermal free convection parameter , and thermal slip parameter (ii)The Nusselt number is escalated for the increasing values of slip parameter , electric parameter , stagnation parameter , and thermal radiation parameter (iii)The concentration field is diminished with an increase in stagnation parameter , electric parameter , diffusion thermal slip parameter , and slip parameter (iv)The Sherwood number is increased as slip parameter , electric parameter , heat generation coefficient , stagnation parameter , and mass free convection parameter are increased

Nomenclature

:Applied magnetic field
:Fluid concentration inside the boundary layer
:Fluid concentration outside the boundary layer
:Skin friction coefficient
:Specific heat
:Concentration at wall surface
:Coefficient of mass diffusion
:Brownian diffusion coefficient
:Thermophoresis diffusion parameter
:Eckert number
:Electrical parameter
:Reduced streamfunction
:Local surface heat flux transfer coefficient
:Microinertia density
:Local mass flux
:Thermal conductivity
:Material parameter
:Magnetic number
:Nusselt number
:Thermophoresis parameter
:Brownian motion parameter
:Angular velocity
:Prandtl number
:Radiative heat flux
:Thermal free convection parameter
:Heat transfer rate
:Thermal radiation parameter
:Local Reynolds number
:Schmidt number
:Boundary layer temperature
:Surface temperature
:Ambient temperature
:Time
:Velocity in direction
:Characteristics velocity
:Velocity in direction
:Stretching rate
:Stagnation parameter
:Slip parameter
:Kinematic viscosity
:Fluid density
:Dynamic viscosity
:Electric charge density
:Dimensionless temperature
:Dimensionless concentration
:Spin gradient viscosity
:Dimensionless boundary layer thickness
:Dimensionless shear stress parameter
:Thermal slip parameter
:Diffusion mass slip parameter
:Heat capacity of the nanoparticles
:Heat capacity of the fluid
:Mass free convection parameter.

Data Availability

The experimental data used to support the findings of this study are included within this paper.

Conflicts of Interest

The authors declare that they have no conflicts of interest regarding the study.

Authors’ Contributions

All authors equally contributed to this work and read and approved the final manuscript.

Acknowledgments

The authors thank Dr. Muhammad Kashif Iqbal for his assistance in proofreading of the manuscript.