# Heat and Mass Transfer in Stagnation Point Flow of Maxwell Nanofluid Towards a Vertical Stretching Sheet with Effect of Induced Magnetic Field

Authors: Tadesse Walelign; Eshetu Haile; Tesfaye Kebede; Gurju Awgichew
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

This paper presents a mathematical model analysis of heat and mass transfer in a two-dimensional flow of electrically conducting, thermally radiative, and chemically reactive Maxwell nanofluid towards a vertical stretching and permeable sheet embedded in a porous medium. Boundary layer approximation and suitable transformations are used to reduce the governing differential equations convenient for computation. Eventually, the transformed nonlinear differential equations along with the corresponding boundary conditions are solved in the framework of optimal homotopy analysis method. The effects of induced magnetic field, buoyancy force, viscous dissipation, heat source, Joule heating, and convective boundary condition are analyzed in detail. The rates of heat, mass, and momentum transfer with respect to the relevant parameters are also examined in terms of the local Nusselt number, Sherwood number, and skin friction coefficients, respectively. Among the many results of the study, it is shown that the induced magnetic field, flow velocity, and temperature profiles are increasing functions of the Maxwell parameter. The results of the present study are also in a close agreement with previously published results under common assumptions.


---

## Body

## 1. Introduction

Based on the deformation of fluids in response to the applied shear stress, fluids can be classified into two categories, namely, Newtonian and non-Newtonian. Fluids such as honey, corn starch, lubricating oils, and paints are non-Newtonian fluids having unpredictable viscosity and responses to applied forces. In 1867, Maxwell developed a constitutive model that is used to describe stress relaxation of air. Now, this model has been extensively utilized to analyze the viscoelastic flow characteristics of many industrially important fluids such as liquid polymers, paints, paper pulps, shampoos, flour dough, certain engine oils, and other viscoelastic fluids. On the other hand, flows near stagnation point are commonly encountered in flows around tips of aircrafts, rockets, ships, submarines, etc. These flows also occur in the cooling of electronic devices and nuclear reactors and in many other hydrodynamic processes. The first classical investigation of stagnation point flow of viscous fluids over a flat stationary plate was made by Hiemenz in 1911. He developed an exact solution to the Navier–Stokes equations governing the flow phenomena. The study of fluid flow, heat transfer, and mass diffusion near stagnation point has significantly contributed in understanding and solving problems in manufacturing and engineering. In order to enhance the heat transfer capabilities of traditional fluids, Choi and Eastman [1] introduced the concept of nanofluids in 1995, and they showed experimentally that embedding of nanometer-sized particles with the common base fluids such as water, oil, and ethylene glycol mixture radically increased the thermal conductivity of the fluid. Due to their wide range of applications, a number of experimental and theoretical investigations have been conducted to examine the flow properties of nanofluids in a variety of flow situations. For instance, the heat transfer phenomena of Nanofluid moving over a wedge surface are presented by Ibrahim and Tulu [2]. Recently, Zhang et al. [3] analyzed the three-dimensional flow of nanofluid among rotating circular plates. Arain et al. [4] deal with the boundary layer flow characteristics of Carreau fluid in the presence of nanoparticles. Further, Alwawi et al. [5] analyzed the effects of magnetic field and nanoparticles in boundary flow of Casson nanofluids.A number of researchers have been interested in conducting several investigations on stagnation point flow of Maxwell nanofluids. For instance, Nagendramma et al. [6] considered a two-dimensional stagnation point flow of Maxwell nanofluid over a stretching surface in the presence of thermal radiation and viscous dissipation. Ramesh et al. [7] carried out a numerical study to analyze and discuss a stagnation point flow of a Maxwell nanofluid towards a permeable stretching sheet. Bai et al. [8] investigated a two-dimensional stagnation point flow of Maxwell nanofluid over a stretching surface by taking both the convective heat transfer and the Brownian motion into account. Haritha et al. [9] investigated the combined effect of Navier slip and convective boundary conditions on an unsteady two-dimensional flow of a Maxwell fluid over a stretching surface in the presence of magnetic field, thermal radiation, heat source, and chemical reaction. Mushtaq et al. [10] studied Buoyancy effects in stagnation point flow of Maxwell fluid utilizing non-Fourier heat flux approach. Recently, Khan et al. [11] investigated the heat and mass diffusion effects in the stagnation point flow of Maxwell nanofluid over a stretchable sheet with slip conditions and chemical reaction. Aziz and Shams [12] studied the heat transfer processes due to the flow of electrically conducting Maxwell nanofluid in the context of volumetric entropy generation.In the aforementioned papers, and in many other early works, the effect of induced magnetic field is usually neglected by assuming an extremely small magnetic Reynolds number and/or getting the mathematical analysis of the problem simple. However, in various practical situations such as power generation, manufacturing of glass, purification of crude oil, and geophysics, the induced magnetic field effect is significant. To this end, extremely few studies are available in open literature that examine the effect of induced magnetic field in stagnation point flow of Maxwell nanofluids. Ibrahim [13] considered the effect of induced magnetic field in examining a stagnation point flow and heat transfer due to upper-convected Maxwell nanofluid over a stretching sheet. In his numerical study, the skin friction coefficient, local Nusselt number, and Sherwood number are found to be decreasing with the increase in the stretching parameter and magnetic field strength. His study also showed that the increase in the convective parameter causes a decrease in the local Sherwood number but an increase in the local Nusselt number. The review of related literature reveals that even though several investigations have been made in the area under consideration, the field remains an active area of research with a number of study gaps. Thus, motivated by the aforementioned early works, this study attempts to examine the influences of magnetic field, thermal radiation, heat source, viscous dissipation, Joule heating, chemical reaction, permeability, and porosity parameters on stagnation point flow of Maxwell nanofluid towards a vertical stretching sheet. Further, a relatively recent and more reliable optimal homotopy analysis method has been implemented successfully to solve the resulting model equations. The results of the present study are also in close agreement with previous works under some restricted assumptions.


## 2. Mathematical Formulation of the Flow System

In this study, it is assumed that as a steady laminar flow of Maxwell nanofluid is approaching a convectively heated vertical sheet, the flow divides into two streams keeping the stagnation point fixed at the origin. Also, a uniform magnetic fieldB of strength B0 is applied normally to the vertical sheet. Magnetic Reynolds number is taken to be large enough to consider the effect of induced magnetic field. Further, an incompressible, electrically conducting, thermally radiative and chemically reactive fluid is considered, and the permeable sheet is stretching with speed Uw=ax, where a>0 is a constant. It is also assumed that the flow system is embedded in a porous medium as illustrated in Figure 1.Figure 1
Sketch of the flow problem.Accordingly, the equations governing the aforementioned flow system can be extended from Ibrahim [13], Hayat and Nadeem [14] as(1)∂u∂x+∂v∂y=0,(2)∂H1∂x+∂H2∂y=0,(3)u∂u∂x+v∂u∂y−μmρfH1∂H1∂x+H2∂H1∂y−HedHedx=Ue∂Ue∂x+υ∂2u∂y2−λ0u2∂2u∂x2+v2∂2u∂y2+2uv∂2u∂x∂y−σB02ρf+υK0u−Ue±1ρfgβ0ρf1−C∞T−T∞−ρp−ρ∞C−C∞,(4)∂H1∂t+u∂H1∂x+v∂H1∂y=H1∂u∂x+H2∂u∂y+1σμm∂2H1∂y2,(5)u∂T∂x+v∂T∂y=α∂2T∂y2+υCp∂u∂y2+1ρCpσ∂H1∂y2+τDB∂C∂y∂T∂y+DTT∞∂T∂y2+σB02ρCpfu−Ue2+16σ∗T∞33ρCpfk∗∂2T∂y2+Q0ρCpfT−T∞,(6)u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2−KrC−C∞,where u,v and H1,H2 denote the components of fluid velocity vector V and induced magnetic field vector H along x and y directions, ρf, μm, λ0, σ, κ and β0 denote density, magnetic permeability, viscoelasticity, electric conductivity, thermal conductivity, and thermal expansion of the nanofluid, respectively; the quantities g and K0 represent the acceleration due to gravity and permeability of the porous medium, respectively; υ=μ/ρf is kinematic viscosity of the nanofluid with μ denoting coefficient of dynamic viscosity; Ue and He are the free stream velocity and induced magnetic field; T and C denote temperature and nanoparticle volume fraction; T∞ and C∞ are the corresponding ambient values of temperature and concentration; α=κ/ρCpf is thermal diffusivity, and τ=ρCpp/ρCpf is the ratio of effective heat capacities of nanoparticle and the ordinary fluid; DB and DT are the Brownian and thermophoresis diffusion coefficients; k∗ and σ∗ are the mean absorption and the Stefan–Boltzmann constants, respectively; the coefficient Q0 stands for heat generation (when Q0>0) and heat absorption (when Q0<0); the rates Kr>0 and Kr<0 denote destructive and constructive reaction rates, respectively.Taking the boundary conditions,aty = 0,(7)u=Uwx=ax,∂H1∂y=H2=0,v=Vwt=−V0,−k∂T∂y=hfTw−T,C=Cw,and as y⟶∞,(8)u⟶uex=bx,v⟶0,H1⟶Hex=xH0,T⟶T∞,C⟶C∞,where Uw,Tw,Cw denote the surface velocity, temperature, and nanoparticle concentration, respectively; Vw is the mass transmission at the surface of the stretching sheet; V0 is the constant value of velocity; k and hf are thermal diffusivity and convective heat transfer coefficients, respectively; H0 is the initial induced magnetic field strength; a and b are positive constants denoting velocity rate of the stretching sheet and the fluid in the free stream, respectively.In order to further simplify the mathematical analysis, we use the following transformations:(9)η=yaυ,ψx,y,t=aυxfη,H1=xH0g′η,H2=−H0υagη,θη=T−T∞Tw−T∞,φη=C−C∞Cw−C∞,where η is a dimensionless similarity variable; f′η,gη,θη and φη denote the dimensionless functions for velocity, induced magnetic field, temperature, and concentration profiles.Introducing the stream functionψx,y with the property u=∂ψ/∂y and v=−∂ψ/∂x, the continuity equation for velocity in equation (1) is identically satisfied. Also, from the definition of H1 and H2 in equation (9), the continuity equation for induced magnetic field equation (2) was satisfied identically. Computing the required partial derivatives and substituting the values into the governing equations (3)–(6), the following system of nonlinear differential equations can be obtained:(10)f‴−ff″−Mg′2−gg″−1+λ2ff′f″−f2f‴+Ha2+Kpδ−f′+rθ−Nrφ+δ2=0,εg‴+g″f−gf″=0,1Pr1+43Rdθ″+fθ′+Nbθ′φ′+Ntθ′2+Ecf″2+Ha2f′2−2δf′+δ2+Qθ=0,φ″+NtNbθ″+Scfφ′−γφ=0,where the prime indicates differentiation with respect to η; M=μmH02/ρυ2 and Ha2=σB02/aρfu−Ue are parameters representing the induced and the external magnetic field strengths, respectively; λ=aλ0≥0 is the Deborah number representing the Maxwell viscoelastic parameter; Kp=υx/UwK0 is the porosity parameter; δ=b/a is a velocity ratio representing the sheet stretching rate, with δ≥1 and δ<1, respectively, implying that the fluid free stream velocity is greater than the sheet velocity and vice versa; r=±β0gTw−T∞1−C∞/a2x is a gravity dependent parameter for assisting (when r>0) and opposing (when r<0) flows; Nr=ρp−ρ∞Cw−C∞/ρfβ0Tw−T∞1−C∞ is the buoyancy ratio; ε=1/μmσυ>0 is the magnetic diffusivity number; Pr=υ/α and Sc=υ/DB are the Prandtl number and the Schmidt number, respectively; Rd=4σ∗T∞3/ρCpfk∗ is the thermal radiation parameter; Ec=Uw2/CpfTw−T∞ is the Eckert number that represents dissipation parameter; Nb=τDBCw−C∞/υ and Nt=τDTTw−T∞/υT∞ are the nanofluid parameters that denote Brownian motion, which represent mass transfer and thermophoresis effects, which represent heat transfer characteristics, respectively; Q=xQ0/ρCpfUw is the heat generation (for Q>0) or absorption (for Q<0) parameter; and γ=Krx/Uw is the chemical reaction parameter, with γ>0 and γ<0 indicating the destructive and generative chemical reaction rates, respectively. Also employing the transformation in equation (9), the boundary conditions in equations (7) and (8) can be reduced as follows:(11)f0=S,f′0=1,g0=g″0=0,θ′0=−Bi1−θ0,φ0=1,f′η⟶δ,g′η⟶1,θη⟶0,φη⟶0,asη⟶∞,where the parameter S=V0/aυ denotes injection and suction for S>0 and S<0, respectively; Bi=hf/kυ/a is the Biot number representing the convective parameter.From a practical point of view, we exploit the behavior of Skin frictionCf, local Nusselt number Nux, and Sherwood number Shx in the region of stagnation point defined, respectively, as(12)Cfx=2τwρfUw2,Nux=xqwκTw−T∞,Shx=xJwDBCw−C∞,where τw=μ01+λ∂u/∂yy=0 is the shear stress; qw=−κ+16σ∗T∞3/3ρCpfk∗∂T/∂yy=0 is the heat flux; and Jw=−DB∂C/∂yy=0 is the mass flux at the stretching surface.Upon substitution, this gives(13)Rex1/2Cfx=21+λf″0,Nux=−1+43RdRex1/2θ′0,Rex−1/2Shx=−φ′0,where Rex=xUw/υ=ax2/υ is the local Reynolds number.


## 3. Method of Solution

The homotopy analysis method (HAM) is one of the relatively recent mathematical methods showing greater efficiency over the past few years. The method is also known for its flexibility and high accuracy to give analytic approximations to the solutions of many important problems in science and engineering. Details of the method can be found in Liao [15, 16]. In order to implement the homotopy analysis method, we first choose the following auxiliary linear operators:(14)Lff=d3fdη3+d2fdη2,Lgg=d3gdη3+d2gdη2,Lθθ=d2θdη2+dθdη,Lφφ=d2φdη2+dφdη,satisfying the conditions(15)LfC1+C2η+C3e−η=0,LgC4+C5η+C6e−η=0,LθC7+C8e−η=0,LφC9+C10e−η=0,where Cii=1−10 are integral constants to be determined by the boundary conditions. Next, based on the linear operators and the boundary conditions, we take the following initial approximations:(16)f0η=S+δη+1−δ1−e−η,g0η=η,θ0η=Bi1+Bie−η,φ0η=e−η.The nonzero auxiliary functions can also be selected as(17)Hfη=Hgη=Hθη=Hφη=e−η.Finally, the parameter valuesδ=ε=1.5,Pr=Sc=1,Ha=0.04,Nb=Nr=Ec=r=0.1,Nc=Bi=0.5,M=λ=Kp=Nt=Rd=Q=S=0.2,andγ=0.1 are used throughout this study unless otherwise stated. In order to do the required computations, we employ a powerful computational software, namely, the BVPh2.0, which is a HAM-based Mathematica package, developed by Zhao and Liao [17].Now, before presenting the results of the study, we give the following convergence analysis of the method. Since convergence of the HAM solutions depends strongly on the proper choice of the so-called convergence-control parameters, we determine the admissible values ofℏf,ℏg,ℏθ and ℏφ as follows. According to Liao [15], one of the mechanisms to identify the range of such values is to sketch the ℏ-curves as shown in Figure 2. This result illustrates that the ℏ-curves are nearly horizontal in the intervals:(18)−1.7<ℏf<−0.2,−1.6<ℏg<0.2,−1.7<ℏθ<−0.1,−1.6<ℏφ<−0.3.Figure 2
Theℏ-curves at 15th order HAM approximations.So, any value in these intervals may be taken to get convergent solutions to our flow problem.More systematically, the optimal values for the convergence-control parameters are determined by reducing the squared residual error:(19)εkℏ≈11+N∑j=0Nℵ∑n=0kunηj2.Upon using the BVPh2.0 Mathematica package, it has been found that(20)ℏf≈−1.4565,ℏg≈−1.2499,ℏθ≈−0.8458,ℏφ≈−1.1501.Further, the respective residual errors and the required values of derivatives are presented against the orders of HAM approximations in Table1.Table 1
Convergence of some HAM solutions forδ=γ=r=Nr=Ec=0.1,Pr=Sc=1,ε=1.5, Bi=0.5,λ=Nb=Nt=Kp=M=Rd=Q=S=0.2,Ha=0.04.
Order of approx.f″0−θ′0−φ′0Squared residual errorsεfεgεθεφ20.19880.28320.74577.06×10−65.38×10−69.40×10−52.82×10−360.19660.28230.78122.19×10−86.59×10−72.12×10−53.08×10−5100.19650.28580.78004.94×10−98.65×10−94.09×10−61.23×10−5140.19640.28760.77923.65×10−91.06×10−89.81×10−74.90×10−6180.19640.28830.77971.79×10−92.25×10−94.44×10−72.08×10−6220.19640.28860.78038.24×10−109.91×10−102.90×10−71.15×10−6260.19640.28860.78083.97×10−104.62×10−101.93×10−77.92×10−7300.19640.28850.78112.07×10−101.96×10−101.21×10−75.79×10−7It is clear from Table1 and Figure 3 that increasing the order of HAM approximation decreases the squared residual errors and leads the series to converge.Figure 3
Total squared residual error for the 20th HAM approximation.The sketch in Figure3 clearly indicated that the total error is decaying extremely fast for increasing orders of the first few iterations. This guarantees that the method is convergent for the selected optimal values of the convergence parameters. To further ensure the validity of our work, we made a comparison of our results with respect to some other previous reports in the absence of extended physical effects as presented in Table 2.Table 2
Comparisons of the present study with previously published work on the coefficientsf″0, −θ0, and −φ0 when ε=1.5,Bi=0.5,δ=1.1,Pr=Sc=1,Ec=r=Nr=Nb=Nt=0.1, λ=Ha=Kp=Rd=Q=γ=S=0 against some values of M.
MCfRe1/2NUxRe−1/2ShxRe−1/2[14]Present study[14]Present study[14]Present study0.00.1683760.16832880.1588170.15880100.7443610.74489270.20.1593330.15945070.1587950.15877880.7439060.74394170.50.1437760.14365470.1587300.15873520.7411250.7421857Table2 verifies that, for the selected values of the parameters, the values of the physical quantities of practical interest determined in this study are in excellent agreement with those of the results reported by Hayat and Nadeem [14].


## 4. Results and Discussions

In this section, we present the most significant results of the study in graphical or tabular forms followed by brief discussions. The impacts of various thermophysical parameters on the dimensionless velocityf′η, induced magnetic field gη, temperature θη, and concentration φη profiles in the stagnation point region are presented as follows:The effects of external magnetic field in the flow field profiles can be expressed in terms of the Hartman numberHa as presented in Figure 4.Figure 4
Impacts of Hartman number on velocity, induced magnetic field, temperature, and concentration profiles.Here, it is worth mentioning that the termσB02/ρu−Ue can be viewed as the combination of the imposed pressure term σB02/ρUe and the Lorentz force σB02/ρu. Unlike free convection, for forced convection regime, the external flow velocity Ue is greater than the velocity u of the fluid in the boundary layer region. That is, the imposed pressure overcomes the Lorentz force. This corresponds to an assisting force in the direction of main flow. Thus, the increase in external magnetic field advances the flow velocity and the induced magnetic field. However, the opposite behavior is observed for temperature profile.The influence of stretching parameterδ on the dimensionless profiles in the region of stagnation point has been investigated and illustrated graphically in Figures 5(a) and 5(b). It is revealed that as the parameter δ increases, the velocity profile also increases for both δ<1 and δ≥1 but with inverted boundary layer structures. This is due to the fact that as the velocity ratio parameter increases, the velocity of the straining motion near the stretching surface also increases from the smooth intersection of the external fluid and the sheet. Also, the induced magnetic field profile was found to be increasing for δ<1 and decreasing for δ≥1. On the other hand, temperature and concentration profiles have remained decreasing as we move away from the stretching sheet.Figure 5
(a) Impacts ofδ on f′η,g′η,θη, and φη for δ<1. (b) Impacts of δ on f′η,g′η,θη, and φη for δ≥1 (right) of the sheet.
(a)(b)The influences of buoyancy ratio parameterNr on dimensionless velocity, induced magnetic field, temperature, and concentration profiles for both assistive and resistive flows are analyzed and given in Figures 6(a) and 6(b).Figure 6
(a) Effects ofNr for resistive (r<0) flow. (b) Effects of Nr for resistive (r<0) assistive (r>0) flow.
(a)(b)For resistive flowr<0, it can be observed from Figure 6(a) that both the velocity and induced magnetic field profiles are increasing functions of Nr. On the other hand, temperature and concentration profiles slowly decrease with Nr. The opposite behavior of the profiles can be noticed for the assistive flow r>0 as shown in Figure 6(b). This is true because, in the assistive flow, the buoyancy force acts against the flow field.The effect of Maxwell viscoelastic parameter expressed in terms of the Deborah numberλ has been investigated and presented in Figures 7(a) and 7(b).Figure 7
(a) Effects ofλ on velocity and induced magnetic field profiles. (b) Effects of λ on temperature and concentration profiles.
(a)(b)As pointed out in Figures7(a) and 7(b), the velocity profile is increasing function of λ. Physically, at higher Deborah numbers, elasticity of the material increases, while viscosity decreases, which enhances the fluid velocity. Similarly, the induced magnetic field also increases with the increase in the values of λ. On the other hand, the temperature gradually decreases, and no significant fluctuation of concentration is observed with increase in the values of λ.The impact of the convective parameterBi has been studied and illustrated in Figure 8. It can be noticed that the velocity, temperature, and concentration profiles can be lowered by increasing the Biot number.Figure 8
Impacts of Biot number on velocity, induced magnetic field, temperature, and concentration profiles.It can be deduced from Figure8 that the influence of Bi is more prominent on temperature profile near the vertical sheet. It is shown that the temperature profile is declining for increasing values of Bi.The reciprocal of magnetic Prandtl numberε is the ratio of the magnetic diffusivity to viscous diffusivity. Its impact in the flow field profiles has been examined and shown in Figures 9(a) and 9(b).Figure 9
(a) Effects ofε on velocity profile. (b) Effects of ε on induced magnetic field profile.
(a)(b)As shown in Figures9(a) and 9(b), the increase in the reciprocal of magnetic Prandtl number ε causes the increase in momentum and induced magnetic field profiles near the stretching surface. This is evident from the relation that the reciprocal of magnetic Prandtl number is inversely proportional to momentum diffusivity, which is low near the stretching sheet.Similarly, the effect of porosity in the flow problem has been investigated and presented as in Figures10(a) and 10(b).Figure 10
(a) Impacts of porosity parameterKp on velocity profile. (b) Impacts of porosity parameter Kp on induced magnetic field profile.
(a)(b)It can be observed in Figures10(a) and 10(b) that the velocity and induced magnetic field profiles increase as porosity of the medium increases in the flow field.The influence of radiationRd in the region of stagnation point flow has been determined as depicted in Figure 11(a) and 11(b).Figure 11
(a) Effects of thermal radiation parameterRd on velocity and induced magnetic field profiles. (b) Effects of thermal radiation parameter Rd on temperature and concentration profiles.
(a)(b)The results in Figures11(a) and 11(b) show that the increment in thermal radiation leads to the growth of thermal boundary layer. Physically, this is due to the fact that larger thermal radiation develops higher heat flux at the surface. It is also noted that as thermal radiation increases, the velocity and induced magnetic field profiles also increase, but the nanoparticle concentration profile reduces in the stagnation point region.Chemical reactions occur in many industrial and engineering activities such as hydrolysis, electroplating, and combustion processes.According to the results in Figures12(a) and 12(b), as the destructive chemical reaction γ>0 increases, the induced magnetic field profile has been increased near the stretching sheet, and the concentration profile decreases in the region of stagnation point. On the other hand, the concentration profile increases with the constructive chemical reaction γ<0 causing the raise of nanoparticle concentration at the sheet. This holds as constructive chemical reactions γ<0 result in consumption, while a destructive chemical reaction γ>0 causes production of the chemical species.Figure 12
(a) Effects of chemical reaction on induced magnetic field profile. (b) Effects of chemical reaction on concentration profile.
(a)(b)In many situations, heat is generated within a body by a chemical, electrical, or nuclear process. The heat source parameterQ has been used to analyze the effects of heat generation and heat absorption as shown in Figures 13(a)–13(d).Figure 13
(a) Effects of heat generation on velocity profile. (b) Effects of heat generation on induced magnetic field profile. (c) Effects of heat generation parameter on temperature profile. (d) Effects of heat absorption parameter on temperature profile.
(a)(b)(c)(d)As presented in Figures13(a) and 13(b), the increase in heat generation results in the increasing of velocity and induced magnetic field profiles. Figures 13(c) and 13(d) show that the temperature of the nanofluid increases with heat generation and decreases with heat absorption. Here, the effect of heat generation yields thermal energy; thereby, the temperature rises, while heat absorption causes the removal of thermal energy in which the temperature profile declines.Viscous dissipation is the conversion of kinetic energy into thermal energy due to the work done against the viscous fluid. The effects of viscous dissipation can be demonstrated by using the Eckert numberEc.It is shown in Figure14(a) that the increases in the dissipation parameter cause the increase in temperature profile. Physically, this is true because the increase in the dissipation effect implies a rise in the frictional heating, which leads to storing of more heat energy in the flow system. Also, Figures 14(b) and 14(c) display that the velocity and induced magnetic field profiles are found to increase with the increase in the dissipation parameter.Figure 14
(a) Impacts of Eckert number on temperature profile. (b) Impacts of Eckert number on velocity profile. (c) Impacts of Eckert number on induced magnetic field profile.
(a)(b)(c)This study also examines the permeability effects in terms of the suction/injection parameterS as depicted in Figures 15(a) and 15(b).Figure 15
(a) Effects ofS on velocity and induced magnetic field. (b) Effects of S on temperature and concentration profiles.
(a)(b)It can be seen from Figures15(a) and 15(b) that both the velocity and induced magnetic field profiles are increasing functions of the injection parameter. On the other hand, the temperature and concentration profiles both decrease with the increase in the injection parameter.The impacts of some parameters on the rates of momentum, heat, and mass transfer near the stretching sheet are examined and expressed in terms of the skin friction, local Nusselt number, and Sherwood number. To this end, the values of the corresponding boundary derivativesf″0, −θ′0, and −φ′0 are computed and provided in Table 3.Table 3
Coefficients of skin friction, Nusselt number, and Sherwood number.
δNrEcεKpRdSγQf″0−θ′0−φ′01.00.015470.287881.026241.51.100460.296961.116042.00.12.400580.282071.239351.52.341570.282911.234633.00.12.292090.283481.231620.22.294540.251201.290540.31.52.297000.218741.349763.02.299980.218681.349964.50.22.302740.218531.350280.42.349230.216931.354730.60.22.394780.215361.359070.52.396430.217551.338320.80.32.397790.217231.325240.42.397790.217231.325240.50.32.478380.222081.385160.42.559190.227011.486930.50.21.058180.296241.019020.41.036380.288530.966440.61.036380.288530.96644-0.22.557510.249851.49336-0.42.556440.259351.47937-0.62.555230.267841.46680-0.12.548620.268391.21123-0.32.546070.268651.11602-0.52.543260.268951.01462Table3 shows that the values of skin friction will grow if any values of δ,Ec,ε,Kp and Rd are increased, or any values of Nr,Q− and γ− are decreased. The local Nusselt number can be increased by increasing Nr,γ and Q− or decreasing Ec,ε and Kp. Also, the local Sherwood number can be maximized by increasing B,Ec,ε and Kp or reducing Nr,Rd,Q− and γ−.


## 5. Concluding Remarks

In this study, the heat and mass transfer phenomena in MHD stagnation point flow of Maxwell nanofluid towards a convectively heated, vertical, and stretching sheet embedded in a porous medium have been analyzed in the presence of induced magnetic field and various thermophysical parameters. The homotopy analysis method (HAM) has been applied to obtain series solutions for the steady two-dimensional laminar flow near the stagnation point. Convergence of the method has been verified. Further, the results obtained in this study are also found in excellent agreement with previous study works under some restricted assumptions. The influences of pertinent parameters on velocity, induced magnetic field, temperature, and concentration profiles are examined as summarized below:(i)
The flow velocity can be accelerated by increasing rates of stretching, external magnetic field, Maxwell parameter, convective parameter, reciprocal of Prandtl number, radiation, porosity, dissipation, heat source, suction, or constructive chemical reaction. It can also be increased by reducing the effects of buoyancy ratio or destructive chemical reaction.(ii)
The induced magnetic field profile can be enhanced by increasing the values of external magnetic field, Maxwell parameter, reciprocal of Prandtl number, radiation, porosity, dissipation parameter, heat source, suction, or constructive chemical reaction. It can also be increased by reducing the effects of buoyancy ratio or stretching parameter.(iii)
The temperature profile can be strengthened near the stagnation point by increasing the values of Maxwell parameter, thermal radiation, buoyancy force, heat source, dissipation parameter, or destructive chemical reaction. The temperature can also be maximized by minimizing the effects of external magnetic field, convective parameter, reciprocal of Prandtl number, porosity of the medium, permeability of the sheet, heat sink, constructive chemical reaction, or stretching parameter.(iv)
The concentration profile of nanoparticles in the region of stagnation point can rise by increasing the buoyancy ratio, heat sink, or constructive chemical reaction. This profile can also be enhanced by decreasing the stretching parameter, radiation, convective parameter, suction, or destructive chemical reaction.(v)
The rate of momentum transfer can be facilitated mainly by increasing the rate of stretching, dissipation parameter, reciprocal of Prandtl number, porosity of the medium, or radiation parameter; it can also be sped up by reducing the effects of buoyancy ratio, heat sink or destructive chemical reaction.(vi)
Heat transfer rate can be enhanced by increasing the buoyancy ratio, chemical reaction, or heat sink. This transfer can also be facilitated by decreasing the dissipation parameter, reciprocal of Prandtl number, or porosity of the medium.(vii)
Mass transfer rate can be sped up by increasing dissipation parameter, reciprocal of Prandtl number, or porosity of the medium. It can also be facilitated by decreasing the effect of buoyancy ratio, radiation, heat sink, or destructive chemical reaction.



---

*Source: 6610099-2021-03-17.xml*