Characterization of the Induced Magnetic Field on Third-Grade Micropolar Fluid Flow Across an Exponentially Stretched Sheet

Abstract

The current research article discusses the two-dimensional, laminar, steady, and incompressible third-grade viscoelastic micropolar fluid flow along with thermal radiation caused by an exponentially stretched sheet. The primary goal of this extensive study is to improve thermal transportation. Thermophoresis and Brownian motion are two key causes of nanoparticle migration in nanofluids, and their impacts on the thermophysical properties of nanofluids are significant. Micropolar fluids are investigated due to their micro-motions that are significant in convective thermal and mass transport polymer formation, nanotechnology, and electronics. The consequences of third-grade fluid parameters, thermophoresis and Brownian motion, induced magnetic field, micro-polarity, and micro-inertia density on the stream of an electrically conductive fluid are analyzed. A homogeneous magnetic field is supplied perpendicularly to the surface, and the liquid is believed to be electrically conducting. As the flow has a significant magnetic Reynolds number, the contribution of the evoked magnetic field is properly accounted in the governing equations. A mathematical model in the form of partial differential equations (PDEs) is built under certain assumptions. By invoking the suitable similarity transformation, the non-linear PDEs are modified into dimensionless coupled ordinary differential equations (ODEs). The MATLAB numerical technique bvp4c is employed to settle the subsequent ODEs together with the boundary constraints. The consequences of numerous physical parameters on the non-dimensional concentration, temperature, micropolar, velocity, and induced magnetic field profiles are portrayed in graphs. It is found that the concentration boundary layer, thermal boundary layer, and micropolar boundary layer thickness decelerate with the increment in the micro-polarity of the fluid.

1 Introduction

Magnetohydrodynamics (MHD) is a branch of fluid mechanics that studies the movement of an electrically conductive liquid in the existence of the magnetic field. Alfven [1] was the first who found that the movement of conducting fluid between magnetic field lines generates potential differences, that, in turn induces the flow of electric currents. The magnetic fields coupled with these electric currents alter the magnetic force that generates them. Alternatively, the fluid stream modifies the system’s electromagnetic structure. The propagation of electric current through a magnetic field, on the other hand, is coupled with a body force, known as the Lorentz force, which affects fluid motion. Scientists and researchers in the field of fluid dynamics have identified its use in metal dispersion, mining, fusion reactors, targeted drug delivery, MHD-based laser beam scanning, construction of MHD pumps, MHD generators, and MHD flow meters. Andersson [2] probed the MHD stream of a viscoelastic liquid past over an expanding sheet and showed that external magnetic field affects the flow similar as viscoelasticity. Hayat and Abbas [3] investigated the radiation impacts on the MHD stream of a viscoelastic liquid across a permeable surface. Nadeem and Hussain [4] explored the shrinking solutions in the existence of MHD stream of a viscous liquid toward a non-linear contracting surface. The effects of an exponentially stretched surface on MHD boundary layer stream and heat transmission characteristics in a permeable medium have been investigated by Ahmad et al. [5]. Hayat et al. [6] discussed the MHD boundary layer stream of a viscous fluid due to non-linear extending cylinder with thermal slip stratification and radiation. Sohail et al. [7] probed the boundary layer stream of the steady MHD Carreau liquid with bioconvection across the heated disk. They determined that the intensity of the magnetic field reduces fluid’s particle velocity, although mass transfer and fluid temperature increase with increasing magnetic field values. Riaz et al. [8] probed the MHD and entropy formation modeling of nanoliquid with three-dimensional (3-D) peristaltic cylindrical confinement. Employing molecular dynamical modeling, AbdulHussein et al. [9] investigated the boiling mechanism of several fluids in micro-channels in the existence of an exterior electromagnetic field.

The study of fluid flow across an expanding sheet is crucial for a variety of applications including, extrusion, glass blowing, cord depiction, copper spiraling, strengthening and tinning of copper wires, glass blowing, warm progressing, thermal conductivity of heat sinks, and melts of high molecular weight polymers. Hussain et al. [10] explored the influences on the boundary layer stream of a micropolar liquid flowing on the stretched surface. Abbas et al. [11] inspected the micropolar fluid flow behavior across the stretched sheet due to a variety of significant applications. Awan et al. have investigated different types of flows over the stretching surface and porous media, with various physical implications [12, 13]. Nadeem and Khan [14] studied the rotating Maxwell nanofluid flow between linear and exponential expanding sheets. Riaz et al. [15] exposed the bioconvection mechanism in the stream of magnetically polarized Williamson nanoparticles along with activating energy and heat source/sink.

Many scientists from all over the world are eager to learn more about non-Newtonian fluid flow. The reason for such a motivation in the study of these fluids is as a result of their use in industries and technologies including suspension fabrication, detergent and paint production, skincare creams, polymer production, spinning of metal etc. The non-Newtonian fluids in comparison to the Newtonian fluids are rheological in structure having non-linear correlation between shearing stress and velocity gradient. Non-Newtonian liquids are classified into three kinds depending upon their characteristics: rate type, differential type, and integral type fluids. The differential type liquid models have found to be well known among them. The second-grade fluid is the most basic subclass of these viscoelastic models, capturing typical stress variations but not anticipating shear thinning/thickening processes. In contrast, the third-grade liquid model can predict both ordinary stress and shear thinning/expanding processes. Abbasbandy et al. [16] investigated the both exact and series approaches for third-grade fluid by using thin film. Hayat et al. [17] probed the rotating stream of third-grade liquid among two permeable sheets by applying MHD effects. Hayat et al. [18] probed the MHD nanofluid stream of second-grade fluid caused by a non-linear stretched surface. In the vicinity of nanoparticles, the non-transient motions of a third-grade fluid driven by a pressure sort die are investigated by Mahanthesh and Joseph [19]. The characteristics of an applied magnetism and entropy formation on Jeffrey nanoliquid in between an annulus region of two small non-concentric pipes was disclosed by Riaz et al. [20]. Mondal et al. [21] discussed nanoliquid stream across a permeable vertical plate along with interior heat production and non-linear thermal radiations. Sangeetha and De [22] analyzed the bioconvection in nanoliquid stream with viscous losses and Ohmic heating. Riaz et al. [23] presented a comparative analysis of entropy assessment on a (3-D) wavy stream of Eyring–Powell nanoliquid.

Micropolar fluids have a microstructure and are coupled to fluids with a non-symmetrical stress tensor. Because of the liquid particles’ localized composition and micro-motions, micropolar liquids display particular microscopic features. They characterize fluids composed of stiff, arbitrarily oriented, or cylindrical particles dispersed in a viscous material, where fluid’s particle distortion is neglected. Eringen [24] was the first to investigate the hypothesis of micropolar fluids. In this hypothesis, the continuum is defined as the set of systematic particles that have not only momentum but also a sub-structure. In other words, each material volume element is made up of micro volume components that can translate and twist independently of each other. Gorla et al. [25] scrutinized the impacts of buoyancy on driven convection in an axially symmetric stagnation stream of micropolar fluids across a vertically placed cylinder. Rehman and Sattar [26] explained the MHD convective stream of a micropolar liquid over the consistently moving permeable sheet. The temperature boundary layer stream caused by a linearly stretched surface submerged in a constant density micropolar liquid together with radiation effects is probed by [27]. Gaffar et al. [28] examined the free convective boundary layer stream of viscoelastic third-grade micropolar liquid passing over a vertically positioned isothermal cone. Ali et al. [29] revealed the significance of MHD on the micropolar nanoliquid stream across an expanding surface along with radiation and heat stratification influences. Jiang et al. [30] conducted the numerical assessment of the passive usage of phase-change elements and the active usage of nanoliquid inside a rectangular channel.

To the extent of the writer’s insights, no research has been conducted to investigate the magnetohydrodynamics (MHD), and heat transmission effects of a viscous, micropolar, and third-grade fluid passed across an exponentially stretched sheet. The current research is presented to analyze the consequences of third-grade fluid parameters, thermophoresis and Brownian motion, induced magnetic field, micro-polarity, and micro-inertia density on the stream of an electrically conductive fluid. The primary purpose of this extensive study is to improve thermal transportation subject to the existence of micro-rotations of tiny nanoparticles. The mathematical model is built under considered flow assumptions and Buongiorno model. The similarity transformation is a technique that is commonly used for solving various flow problems, in which the system of partial differential equations (PDEs) is modified into ordinary differential equations (ODEs) to solve them analytically or numerically. The current model is simplified by applying similarity transformation. The resultant ODEs are numerically tackled in MATLAB using the bvp4c algorithm, and the numerical tables are constructed to ensure the validity of results. The obtained findings have significant engineering and technological applications.

2 Mathematical Analysis and Flow Geometry

2.1 Constitutive Model for Third-Grade Fluid

For an incompressible fluid, the equations of motion and continuity are as follows:

Here, ρ, V, T, b, J, and B are density of the fluid, velocity vector, Cauchy stress tensor, body force, electric current, and external magnetic field, respectively. Following Rajagopal and Fosdick [31], the third-grade fluid’s stress tensor is as follows:Here, p, I, T, and S are the pressure, identity tensor, Cauchy stress tensor, and the extra stress tensor, respectively. Furthermore, are metallic constants and B_i(i = 1, 2, 3) is the kinematic tensor defines as:Here, represents the substantial derivative, and it is stated as:

In the case of third-grade fluid, the material moduli satisfy Clausius–Duhem inequality stated as:

So, the stress tensor takes the following form.

2.2 Constitutive Model for Micropolar Fluid

The field equations for the micropolar fluid following Papautsky et al. [32] are stated as follows:

Eqs. 8–10, respectively, denote conservation of mass, linear momentum, and angular momentum. Furthermore, ρ represents the micropolar fluid density, G is the (gyration) micro-rotation vector, V depicts velocity vector, j represents micro-inertia, p is the pressure, l denotes body couple per unit mass vector , λ∗ reflects Eringen second-degree viscosity parameter, f depicts the body strength per unit mass vector, κ represents coefficient of vortex viscosity, μ denotes dynamic viscosity, and , γ∗, and denote the swirl gradient viscosity coefficient. Also, (′) denotes the time derivative. For and in the absence of l and f, the micro-rotation vector G approaches zero, and Eq. 9 is simplified to the Navier–Stokes equations. For κ = 0, the gyration vector G and velocity vector V become detached, and the micro-movements have no influence on the entire motion of the fluid.

2.3 Problem Formulation

In this section, the steady, laminar, two-dimensional, and incompressible third-grade viscoelastic micropolar fluid stream along with thermal radiation caused by an exponentially stretched surface is reported. The x-axis is marked parallel to the sheet, while the y-axis is directed orthogonal to the sheet. The uniform magnetic field is supplied to the surface in a normal direction, and the fluid is believed to be electrically conductive. The sheet and the third-grade micropolar liquid are both originally maintained at the identical temperature. The temperature is elevated to T_f, and T_f > T_∞ where T_∞ is the ambient temperature, which remains constant as (see Figure 1). The concentration at the sheet is taken as C_f, while the ambient concentration is C_∞ (C_f > C_∞). This is supposed that the flow’s magnetic Reynolds number is not low in magnitude, and thus the produced magnetic field is not negligible. The evoked magnetic field H₂ is also assumed to be supplied along the x-axis. The parallel component H₁ of the induced magnetic field tends to H_e(x) in the free stream flow, and the slope of the magnetic field approaches zero along y-axis.

FIGURE 1

Using Eq. 7 in Eq. of motion (Eq. 2) along with the boundary layer estimations [33–36] in case of third-grade fluid, notably, inside the boundary layer , and u are O(1), v and y are O(δ), and ν be O(δ²), and being O(δ⁴), and the components of O(δ) are ignored (δ is boundary layer width). The equations for conservation of mass, magnetic field, and momentum are given as Eqs. 11–12. Equations 13–14 rise from the conservation of energy and concentration. Eq. 15 rises from micropolar boundary layer approximation by Papautsky et al. [32], and Eq. 16 resulted from the conservation of induced magnetic field. Hence, the mathematical model overseeing the considered flow problem is as follows:

The concerned boundary conditions are as follows:

Here, u, v are the velocity coefficients in x and y directions, ρ depicts the liquid’s density , k denotes the vertex viscosity, μ represents dynamic viscosity, N denotes the micropolar fluid’s angular velocity, T represents the temperature, U₀ depicts the reference velocity, C depicts the concentration, B₀ denotes the uniform magnetic field, η₀ denotes the magnetic diffusivity, σ denotes the electrical conductance, C_p represents specific heat, and τ* reveals the quotient of the latent heat of the nanoparticles to the latent heat of the base liquid. Furthermore, D_B reflects the coefficient of Brownian motion, D_T denotes the coefficient of thermophoresis diffusion, and the micro-inertial density is depicted by j. The viscosity of the spin gradient γ* is given as follows:where is the micropolar parameter. Temperature-based thermal conductivity is defined as follows:

The Rosseland radiative heat flow (q_r) along y-axis is defined as:

Here, Stefan–Boltzman constant and average absorption coefficient are expressed by σ∗ and k∗, respectively.

To simplify the analysis, the following appropriate similarity transformations are implemented.

Here, ψ and ψ₁ are stream functions and can be represented as:where η denotes the dimensionless variable, and f′(η), θ(η), ϕ(η), h(η), and g(η) are the dimensionless velocity curve, temperature curve, concentration curve, micropolar curve, and induced magnetic field curve, respectively. Implying Eqs. 23, 24, the continuity equation and magnetic flux equation holds true, and the Eqs. 12–16 are modified into the non-dimensional ODEs as stated:

The associated non-dimensional boundary constraints are given as:

The non-dimensional parameters are defined as:

Here, α₁, α₂ are the non-dimensional viscoelastic coefficients, β is the third-grade fluid constant, M depicts the magnetic parameter, δ denotes the thermal slip parameter , Ec represents the Eckert no, Rd is the radiation constant, Nt represents the thermophoresis constant, K denotes the micropolar fluid coefficient, Pr stands for Prandtl number, Nb indicates the Brownian diffusivity coefficient, Le denotes the Lewis no, B is the micro-inertia density coefficient, and Pr_m denotes the magnetic Prandtl number. Physical variables of interest, such as local couple stress (M_s), Sherwood number (Sh_x), Nusselt number (Nu_x), and skin friction coefficient (Cf_x) are defined as follows:where

Utilizing the similarity transformation (Eq. 22), the aforementioned expressions in non-dimensional aspects are stated as follows:where X = x/l and is the local Reynolds number.

3 Numerical Procedure

Due to the extreme non-linearity, the coupled ODEs Eqs. 25–29 along with the boundary constraints (Eq. 30) are unable to solve analytically. To address non-linear boundary value problems, different numerical techniques are applied in MATLAB. The bvp4c approach is a useful technique to solve such problems numerically. The underlying partial differential equations (PDEs) are transfigured into ODEs utilizing the similarity analysis. The resulting system of ODEs is resolved numerically employing the built-in bvp4c technique in MATLAB. The solution strategy is described as follows:

Utilizing the aforementioned notations, the coupled ODEs are transformed into the following first-order ODEs:

The concerned boundary constraints in MATLAB script are stated as:

The MATLAB code is processed to get the numerical solutions. In order to determine the validity of numerical technique, we have analyzed our outcomes for limiting cases with already existing research work. The outcomes are reported to be in excellent concordance as shown in Tables 1, 2. It is also probed that for β = α₂ = 0, the mathematical model reduces to the second-grade fluid model, and for β = α₁ = α₂ = 0, we get the equations of motion for the classical viscous liquid.

4 Computational Results

In this article, the micropolar, third-grade, nanofluid flow past across an exponentially expanding surface is considered. The mathematical model is built by analyzing fluid flow assumptions. The system of dimensionless ODEs (Eqs. 25–29) are solved numerically, utilizing the bvp4c strategy in MATLAB in addition to the boundary constraints (Eq. 30), and the impacts of numerous physical variables are explored. The significance of these coefficients on the temperature curve θ(η), micro-rotation curve h(η), velocity curve f′(η), produced magnetic field profile g(η), and concentration profile ϕ(η) is highlighted via tables and figures.

4.1 Velocity Profile f′(η)

Figure 2 depicts the tendency of velocity distribution with the rising values of material constant α₁. It is clear that velocity and momentum boundary layer thickness increases for the higher values of viscoelastic constant α₁. Physically, the viscosity of the liquid is inversely proportional to the material constant α₁ because of that when the stress is applied, it reduces the strain and strengthens the elastic effects between the adjacent layers, and hence the velocity profile enhances. Figure 3 denotes the impact of second material constant α₂ on the velocity profile. The velocity curve is the decreasing function for larger values of α₂. This parameter causes shear thickening of the fluid and an increase in resistance that reduces the boundary layer flow, and originates a decrement in the size of the momentum boundary layer width. Figure 4 exhibits the ascending behavior of velocity with the augmentation of the third-grade fluid coefficient β. This coefficient is inversely proportionate to the square of the liquid viscosity. Alternatively, higher β readings indicate superior third-grade material characteristics (higher liquid elasticity) and smaller fluid viscosity. This causes the boundary layer stream to accelerate, resulting in higher f′(η) values. As this value is raised, the fluid needs less stress to flow, encouraging flow acceleration. In Figure 5, the influence of the magnetic field coefficient on velocity profile is examined. It is analyzed that the velocity curve declines monotonically with the augmented values of M ,and the velocity diminishes far away from the sheet. This process assists in controlling the size of the boundary layer. This is because of the fact that the existence of magnetism in an electrically conductive liquid generates a force known as the Lorentz force that operates opposite to the flow direction and forces velocity profile to decline. The velocity diminishes as the retardation to the flow enhances. In Figure 6 the influence of micropolar fluid parameter on the structure of velocity distribution is noticed. The increase in vortex viscosity strongly accelerates the fluid flow. The results show that the momentum transfer layer-by-layer is significantly affected by the rise in viscosity caused by the micro-rotation of the molecules. The micro-elements rotate more strongly, which helps to accelerate the liquid motion in the boundary layer. As a result, linear momentum diffusion is enhanced by micro-polarity, which explains why suspension fluids have thinner boundary layers than regular fluids.

FIGURE 2

FIGURE 3

FIGURE 4

FIGURE 5

FIGURE 6

4.2 Temperature Profile θ(η)

Figure 7 explains that temperature in the boundary layer increases with progressing micropolar coefficient, that is, higher K values. The regime is greatly heated, and the size of the temperature boundary layer enhances. The vortex’s enhanced viscosity promotes thermal diffusion and serves as a rotator. This boosts up the capacity of thermal diffusion within the fluid’s regime from the micro to the macro level and rapidly carries heat from the sheet boundary into the liquid body with greater intensity. Figure 8 depicts the relationship between the temperature profile and the Prandtl number Pr. The fraction of the momentum diffusion coefficient to heat diffusion coefficient is called the Prandtl number. This is examined that the temperature curve declines with increasing the Prandtl number. The higher values of Prandtl number affect the thermal diffusion. Prandtl number is inversely proportionate to the thermal diffusion; hence, higher Pr values indicate lower thermal diffusion, leading in lower temperature and a weaker temperature boundary layer. When the Prandtl number increases, rate of thermal conductivity gets lower. Consequently, heat is dissipated more quickly, and hence, the temperature boundary layer width and temperature of the liquid both diminish. As a result, the Prandtl number is employed to enhance the cooling tendency of fluids. The temperature-dependent thermal conductivity parameter’s effect over the temperature curve is explored in Figure 9. This is evident that the temperature curve is increasing with the augmentation of ϵ. The consequence of various Eckert number values over the temperature curve is explored in Figure 10. Ec is known as the quotient of the kinetic energy and enthalpy of heat transfer. As the Eckert number (Ec) increases, the liquid’s temperature rises. The fluctuation in the thermal curve θ(η) is examined in Figure 11 in respect to distinct amounts of the radiation parameter Rd. It is identified that a higher thermal radiation constant is related to a greater temperature and a wider temperature boundary layer. As the temperature field grows, the enhanced radiation transfers a significant amount of heat to the fluid. Physically, the involved radiation produces more heat which raises the fluid’s temperature. The influence of thermophoresis coefficient over the thermal distribution is probed in Figure 12. Thermophoresis is a phenomenon that occurs in mixture of sub-micron sized particles, in which the different particles respond differently to the force of a thermal gradient. The particle’s velocity is known as the thermophoretic velocity, and the stress exerted on the dispersed particles because of the temperature difference is known as the thermophoretic force. The thermophoretic force increases as the value of Nt grows, causing the temperature to increase.

FIGURE 7

FIGURE 8

FIGURE 9

FIGURE 10

FIGURE 11

FIGURE 12

4.3 Concentration Profile ϕ(η)

The influence of Brownian diffusivity and thermophoresis coefficients over the concentration distribution ϕ(η) is portrayed in Figures 13, 14. According to these graphs, nanofluid constants have opposite impacts on concentration distributions. Particularly, as the thermophoresis constant rises, the width of the concentration boundary layer enhances; nevertheless, as the Brownian motion constant increases, the ϕ(η) values decrease. The thermophoresis constant amplifies the thermophoretic force, resulting in the transfer of nanoparticles from warm to cool locations and an increment in nanoparticle volume. Thermophoresis has numerous uses, including radioactive particle deposition in nuclear reactors, silicon thin film deposition, and aerosol technologies. Furthermore, the progressing amount of the Brownian motion coefficient reduces the micro-mixing of nanoparticles into the fluid’s zone, which diminishes the boundary layer thickness of concentration distribution.

FIGURE 13

FIGURE 14

The behavior of ϕ(η) for the numerous values of Lewis number (Le) is examined in Figure 15. The Lewis number is stated as the rate of heat-to-mass diffusion coefficient. It is utilized to express the flow of liquid, in which heat and momentum transfer occur simultaneously. The concentration distribution becomes steeper when Lewis number is increased. The higher values of ‘Le’ imply the lesser values of mass diffusivity ‘D_B,’ which causes a weaker penetration depth for the concentration boundary layer. In Figure 16, the concentration curve for the various amounts of the micropolar coefficient is explored. However, when K increases, the micro-rotation velocity ϕ(η) drops down within a shorter range in the concentration boundary layer.

FIGURE 15

FIGURE 16

4.4 Micropolar Profile h(η)

Figures 17, 18 are portrayed for the different values of micro-inertia density coefficient B and micropolar coefficient K. This is clear that the micropolar (angular) speed of the sub-micron sized particles decreases with the increase in the micro-polarity of the fluid. Similarly, the micro-inertia density coefficient B serves as a restricting force for the micro-rotation profile h(η). So, for the progressing amounts of B, the micropolar boundary layer width declines. This is possible because micropolar fluids provide a high barrier to fluid motion.

FIGURE 17

FIGURE 18

4.5 Induced Magnetic Field Profile h(η)

The influence of magnetic Prandtl number Pr_m over induced magnetic field curve g(η) is exhibited in Figure 19 , where a rise in Pr_m leads to an improvement in the produced magnetic field distribution. This is because the magnetic diffusivity over the boundary layer’s surface decreased while the fluid viscous dispersion rate enhanced. Magnetic Prandtl no.(Pr_m) is a non-dimensional quantity in Magnetohydrodynamics that estimates the ratio of momentum diffusion coefficient (ν) to magnetic diffusion coefficient (η).

FIGURE 19

4.6 Numerical Results

The influence of all the non-dimensional coefficients used in the considered article on the skin friction parameter , Nusselt number. , Sherwood number. , and coupled wall stress (Re_xM_s) are presented in Tables 3, 4. Table 3 probes the significance of involved coefficients on the skin friction. The values of the coefficients other than used in Table 3 are taken to be fixed Pr = 0.7, ϵ = 0.3, Ec = 0.1, K = 0.2, Rd = 0.3, Nt = 0.3, Nb = 1.0, and Le = 1.0. By increasing the value of material parameters α₁ and α₂, the skin friction coefficient reduces. It is because of the reason that the viscosity near the surface of sheet is enhanced with the progression in material constant, which causes the decline in values. As the rise in shear thickening coefficient β, improves the boundary layer thickness that give rise to the skin friction coefficient. The increase in micropolar parameter K, give rise to the value of , but the opposite behavior is seen for couple stress Re_xM_s. A reduction is caused by a rise in the magnetic field constant M in both and Re_xM_s values. By increasing the micro-inertia density B, the values of get enhanced but reverse behavior is seen for couple stress Re_xM_s.

Table 4 depicts the influence of different parameters on the Nusselt number and Sherwood number. The other parameters are taken to be fixed β = 0.5, α₁ = 0.1, B = 0.5, α₂ = 0.5, K = 0.2, and M = 0.3. The increment in the values of Prandtl number causes a decline in the thermal diffusivity, and hence resists the rise in the heat transmission rate at the boundary. Boosting the Prandtl number raises the average Nusselt number at the heated surface, while a reverse effect is examined for the . As ϵ is enhanced, the Nusslet number decreases while the Sherwood number increases. The increase in the value of Ec causes the Nusselt number to decrease while the Sherwood number increases. The boosting values of micropolar parameter K enhances the Nusselt number but declines the Sherwood number. Higher values of the Rd parameter improve convective heat and matter transmission, which ultimately increases the average Nusselt and Sherwood numbers values.

The higher values of thermophoresis coefficient Nt causes a decline in both the and values. The increment in the Nb values results in the enhancement of Sherwood number because the mass transfer rate from higher to lower concentration area get increased. The growth in the values of Lewis number Le improves the thermal diffusion and declines the mass conductivity. Thus, the values of decrease and that of increase.

5 Conclusion

References

• 1.

  AlfvénH. Existence of Electromagnetic-Hydrodynamic Waves. Nature (1942) 150(3805):405–6. 10.1038/150405d0

  • CrossRef
  • Google Scholar

• 2.

  AnderssonHI. MHD Flow of a Viscoelastic Fluid Past a Stretching Surface. Acta Mech (1992) 95(1):227–30. 10.1007/bf01170814

  • CrossRef
  • Google Scholar

• 3.

  AbbasZHayatT. Radiation Effects on MHD Flow in a Porous Space. Int J Heat Mass Transf (2008) 51(5-6):1024–33. 10.1016/j.ijheatmasstransfer.2007.05.031

  • CrossRef
  • Google Scholar

• 4.

  NadeemSHussainA. MHD Flow of a Viscous Fluid on a Nonlinear Porous Shrinking Sheet with Homotopy Analysis Method. Appl Math Mech.-Engl Ed (2009) 30(12):1569–78. 10.1007/s10483-009-1208-6

  • CrossRef
  • Google Scholar

• 5.

  AhmadISajidMAwanWRafiqueMAzizWAhmedMet alMHD Flow of a Viscous Fluid over an Exponentially Stretching Sheet in a Porous Medium. J Appl Maths (2014) 2014:1–8. 10.1155/2014/256761

  • CrossRef
  • Google Scholar

• 6.

  KhanMITamoorMHayatTAlsaediA. MHD Boundary Layer thermal Slip Flow by Nonlinearly Stretching cylinder with Suction/blowing and Radiation. Results Phys (2017) 7:1207–11. 10.1016/j.rinp.2017.03.009

  • CrossRef
  • Google Scholar

• 7.

  SohailMNazirUChuYMAl-KouzWThounthongP. Bioconvection Phenomenon for the Boundary Layer Flow of Magnetohydrodynamic Carreau Liquid over a Heated Disk. Sci Iran (2021) 28(3):1896–907. 10.24200/SCI.2021.53970.3518

  • CrossRef
  • Google Scholar

• 8.

  RiazAAbbasTZeeshanADoranehgardMH. Entropy Generation and MHD Analysis of a Nanofluid with Peristaltic Three Dimensional Cylindrical Enclosures. Hff (2021) 31:2698–714. 10.1108/HFF-11-2020-0704

  • CrossRef
  • Google Scholar

• 9.

  AbdulHusseinWAAbedAMMohammedDBSmaisimGFBaghaeiS. Investigation of Boiling Process of Different Fluids in Microchannels and Nanochannels in the Presence of External Electric Field and External Magnetic Field Using Molecular Dynamics Simulation. Case Stud Therm Eng (2022) 35:102105. 10.1016/j.csite.2022.102105

  • CrossRef
  • Google Scholar

• 10.

  HussainMAshrafMNadeemSKhanM. Radiation Effects on the thermal Boundary Layer Flow of a Micropolar Fluid towards a Permeable Stretching Sheet. J Franklin Inst (2013) 350(1):194–210. 10.1016/j.jfranklin.2012.07.005

  • CrossRef
  • Google Scholar

• 11.

  AbbasNSaleemSNadeemSAlderremyAAKhanAU. On Stagnation point Flow of a Micro Polar Nanofluid Past a Circular cylinder with Velocity and thermal Slip. Results Phys (2018) 9:1224–32. 10.1016/j.rinp.2018.04.017

  • CrossRef
  • Google Scholar

• 12.

  AwanAUAbidSUllahNNadeemS. Magnetohydrodynamic Oblique Stagnation point Flow of Second Grade Fluid over an Oscillatory Stretching Surface. Results Phys (2020) 18:103233. 10.1016/j.rinp.2020.103233

  • CrossRef
  • Google Scholar

• 13.

  AwanAUAbidSAbbasN. Theoretical Study of Unsteady Oblique Stagnation point Based Jeffrey Nanofluid Flow over an Oscillatory Stretching Sheet. Adv Mech Eng (2020) 12(11):1687814020971881. 10.1177/1687814020971881

  • CrossRef
  • Google Scholar

• 14.

  KhanMNNadeemS. A Comparative Study between Linear and Exponential Stretching Sheet with Double Stratification of a Rotating Maxwell Nanofluid Flow. Surf Inter (2021) 22:100886. 10.1016/j.surfin.2020.100886

  • CrossRef
  • Google Scholar

• 15.

  AbdelmalekZKhanSUWaqasHRiazAKhanIATliliI. A Mathematical Model for Bioconvection Flow of Williamson Nanofluid over a Stretching cylinder Featuring Variable thermal Conductivity, Activation Energy and Second-Order Slip. J Therm Anal Calorim (2021) 144(1):205–17. 10.1007/s10973-020-09450-z

  • CrossRef
  • Google Scholar

• 16.

  AbbasbandySHayatTMahomedFMEllahiR. On Comparison of Exact and Series Solutions for Thin Film Flow of a Third-Grade Fluid. Int J Numer Meth Fluids (2009) 61(9):987–94. 10.1002/fld.1994

  • CrossRef
  • Google Scholar

• 17.

  HayatTNazRAlsaediARashidiMM. Hydromagnetic Rotating Flow of Third Grade Fluid. Appl Math Mech.-Engl Ed (2013) 34(12):1481–94. 10.1007/s10483-013-1761-7

  • CrossRef
  • Google Scholar

• 18.

  HayatTAzizAMuhammadTAhmadB. On Magnetohydrodynamic Flow of Second Grade Nanofluid over a Nonlinear Stretching Sheet. J Magnetism Magn Mater (2016) 408:99–106. 10.1016/j.jmmm.2016.02.017

  • CrossRef
  • Google Scholar

• 19.

  MahantheshBJosephTV. Dynamics of magneto-nano third-grade fluid with Brownian motion and thermophoresis effects in the pressure type die. J Nanofluids (2019) 8(4):870–5. 10.1166/jon.2019.1642

  • CrossRef
  • Google Scholar

• 20.

  RiazAEllahiRSaitSMMuhammadT. Magnetized Jeffrey Nanofluid with Energy Loss in between an Annular Part of Two Micro Non-concentric Pipes. Energy Sourc A: Recovery, Utilization, Environ Effects (2020) 2020:1–20. 10.1080/15567036.2020.1798568

  • CrossRef
  • Google Scholar

• 21.

  MondalHDePGoqoSSibandaP. A Numerical Study of Nanofluid Flow over a Porous Vertical Plate with Internal Heat Generation and Nonlinear thermal Radiation. J Por Media (2020) 23(6):517–29. 10.1615/JPorMedia.2020026624

  • CrossRef
  • Google Scholar

• 22.

  SangeethaEDeP. Bioconvection in Nanofluid Flow Embedded in Non-darcy Porous Medium with Viscous Dissipation and Ohmic Heating. J Por Media (2021) 24(1):15–23. 10.1615/JPorMedia.2020036165

  • CrossRef
  • Google Scholar

• 23.

  RiazAZeeshanABhattiMM. Entropy Analysis on a Three-Dimensional Wavy Flow of Eyring-Powell Nanofluid: A Comparative Study. Math Probl Eng (2021) 2021:1–14. 10.1155/2021/6672158

  • CrossRef
  • Google Scholar

• 24.

  EringenA. Theory of Micropolar Fluids. Indiana Univ Math J (1966) 16(1):1–18. 10.1512/iumj.1967.16.16001

  • CrossRef
  • Google Scholar

• 25.

  Subba Reddy GorlaRMansourMAMohammedienAA. Combined Convection in an Axisymmetric Stagnation Flow of Micropolar Fluid. Int J Numer Methods Heat Fluid Flow (1996) 6(4):47–55. 10.1108/09615539610123441

  • CrossRef
  • Google Scholar

• 26.

  RahmanMMSattarMA. Magnetohydrodynamic Convective Flow of a Micropolar Fluid Past a Continuously Moving Vertical Porous Plate in the Presence of Heat Generation/absorption. J Heat Transfer (2006) 128(2):142–52. 10.1115/1.2136918

  • CrossRef
  • Google Scholar

• 27.

  IshakA. Thermal Boundary Layer Flow over a Stretching Sheet in a Micropolar Fluid with Radiation Effect. Meccanica (2010) 45(3):367–73. 10.1007/s11012-009-9257-4

  • CrossRef
  • Google Scholar

• 28.

  GaffarSABégOAPrasadVR. Mathematical Modeling of Natural Convection in a Third-Grade Viscoelastic Micropolar Fluid from an Isothermal Inverted Cone. Iran J Sci Technol Trans Mech Eng (2020) 44(2):383–402. 10.1007/s40997-018-0262-x

  • CrossRef
  • Google Scholar

• 29.

  AliLLiuXAliBMujeedSAbdalSMutahirA. The Impact of Nanoparticles Due to Applied Magnetic Dipole in Micropolar Fluid Flow Using the Finite Element Method. Symmetry (2020) 12(4):520. 10.3390/sym12040520

  • CrossRef
  • Google Scholar

• 30.

  JiangYSmaisimGFMahmoudMZLiZAybarHŞAbedAM. Simultaneous Numerical Investigation of the Passive Use of Phase-Change Materials and the Active Use of a Nanofluid inside a Rectangular Duct in the thermal Management of Lithium-Ion Batteries. J Power Sourc (2022) 541:231610. 10.1016/j.jpowsour.2022.231610

  • CrossRef
  • Google Scholar

• 31.

  FosdickRLRajagopalKR. Thermodynamics and Stability of Fluids of Third Grade. Proc R Soc Lond A (1980) 369(1738):351–77. 10.1098/rspa.1980.0005

  • CrossRef
  • Google Scholar

• 32.

  PapautskyIBrazzleJAmeelTFrazierAB. Laminar Fluid Behavior in Microchannels Using Micropolar Fluid Theory. Sens Actuator A Phys (1999) 73(1-2):101–8. 10.1016/S0924-4247(98)00261-1

  • CrossRef
  • Google Scholar

• 33.

  SajidMHayatTAsgharS. Non-similar Analytic Solution for MHD Flow and Heat Transfer in a Third-Order Fluid over a Stretching Sheet. Int J Heat Mass Transfer (2007) 50(9-10):1723–36. 10.1016/j.ijheatmasstransfer.2006.10.011

  • CrossRef
  • Google Scholar

• 34.

  PakdemirliM. The Boundary Layer Equations of Third-Grade Fluids. Int J Non-Linear Mech (1992) 27(5):785–93. 10.1016/0020-7462(92)90034-5

  • CrossRef
  • Google Scholar

• 35.

  LoganathanKMohanaKMohanrajMSakthivelPRajanS. Impact of Third-Grade Nanofluid Flow across a Convective Surface in the Presence of Inclined Lorentz Force: an Approach to Entropy Optimization. J Therm Anal Calorim (2021) 144(5):1935–47. 10.1007/s10973-020-09751-3

  • CrossRef
  • Google Scholar

• 36.

  BégOABakierAYPrasadVRZuecoJGhoshSK. Nonsimilar, Laminar, Steady, Electrically-Conducting Forced Convection Liquid Metal Boundary Layer Flow with Induced Magnetic Field Effects. Int J Therm Sci (2009) 48(8):1596–606. 10.1016/j.ijthermalsci.2008.12.007

  • CrossRef
  • Google Scholar

• 37.

  ChaudharySSinghSChaudharyS. Thermal Radiation Effects on MHD Boundary Layer Flow over an Exponentially Stretching Surface. Appl Maths (2015) 06(02):295–303. 10.4236/am.2015.62027

  • CrossRef
  • Google Scholar

• 38.

  HayatTKhanMIWaqasMAlsaediAYasmeenT. Diffusion of Chemically Reactive Species in Third Grade Fluid Flow over an Exponentially Stretching Sheet Considering Magnetic Field Effects. Chin J Chem Eng (2017) 25(3):257–63. 10.1016/j.cjche.2016.06.008

  • CrossRef
  • Google Scholar

• 39.

  ElbashbeshyE. Heat Transfer over an Exponentially Stretching Continuous Surface with Suction. Arch Mech (2001) 53(6):643–51.

  • Google Scholar

• 40.

  AliBNaqviRAHussainDAldossaryOMHussainS. Magnetic Rotating Flow of a Hybrid Nano-Materials Ag-MoS2 and Go-MoS2 in C2H6O2-H2o Hybrid Base Fluid over an Extending Surface Involving Activation Energy: FE Simulation. Mathematics (2020) 8(10):1730. 10.3390/math8101730

  • CrossRef
  • Google Scholar

• 41.

  KumariMNathG. Transient Rotating Flow over a Moving Surface with a Magnetic Field. Int J Heat Mass Transfer (2005) 48(14):2878–85. 10.1016/j.ijheatmasstransfer.2005.02.009

  • CrossRef
  • Google Scholar