# Investigation of Effects of Thermal Radiation, Magnetic Field, Eckert Number, and Thermal Slip on MHD Hiemenz Flow by Optimal Homotopy Asymptotic Method

Authors: Solomon Bati Kejela; Mitiku Daba; Abebe Girum
Journal: Journal of Mathematics
 (2021)

---

## Abstract

Analytical investigation of thermal radiation, Prandtl number, Eckert number, permeability parameter, magnetic field, velocity, and thermal slip effects on magnetohydrodynamic Hiemenz flow over a permeable plate with forced convection has been presented. Similarity variable conversion method has been applied to transmute the fundamental governing equations of the fluid dynamics in flow into a pair of nonlinear third-order ordinary differential equations and is analytically solved by the optimal homotopy asymptotic method (OHAM). The influences of several relevant physical parameters in the model on velocity and temperature of the fluid have been studied and analysed profoundly by use of graphs and tables. It is detected that, with mounting value of suction/blowing parameter and magnetic field parameter, the skin friction coefficient enhances. Likewise, it is seen that the Nusselt number increases with enhancing value of magnetic parameter. It is also witnessed that the velocity increases as the Eckert number, blowing/suction parameter, and permeability parameter increase, but it decays against magnetic field and velocity slip parameter. Moreover, the result reveals that the fluid temperature upsurges along with snowballing the radiant heat, magnetic field parameter, and the Eckert number. However, it descends against thermal slip parameter, Prandtl number, wall temperature exponent, and velocity slip parameter. A comparison with previous studies has been made, and the result shows an excellent agreement.


---

## Body

## 1. Introduction

The study of Hiemenz (stagnation-point) flow plays an important role in the study of plentiful natural and industrial phenomena because of its applications in discovering flows over the tips of submarines, the tip of ships, and aircraft. It is also vital in various engineering disciplines like hydrodynamic processes, refrigeration of electronic materials using fans, freezing of nuclear apparatus, and so forth. The theoretical inquiry of the tide of Newtonian fluid near a two-dimensional stop off-point flow in the route of a motionless semi-infinite barricade has been studied before all by Hiemenz [1]. Hiemenz applied a similarity variable conversion method to reduce the Navier Stokes equations to equivalent solvable nonlinear ODEs. The study of magnetohydrodynamic (MHD) boundary layer ﬂow on a continuous stretching sheet has numerous applications in industrial manufacturing processes such as the aerodynamic extrusion of plastic sheets, wire drawing, glass ﬁber and paper production, liquid ﬁlm, hot rolling, drawing of plastic ﬁlms, metal, polymer extrusion, and metal spinning. The MHD stagnation-point flow problems have theoretical and practical applications in manufacturing processes like boundary layer along material handling conveyers, blood flow problems, extrusion of plastic sheets, geothermal energy extractions, cooling of infinite metallic plate in cooling bath, petroleum industries, and so forth; due to this, they have attracted the attention of many researchers. Ishak et al. [2] discussed the MHD stagnation-point flow towards a stretching surface with variable surface temperature, whereas the MHD effect past a shrinking sheet has been analysed by Mahapatra et al. [3]. Since then, several studies have been conducted on stretching/shrinking problems [4–12].The effects of thermal radiation and heat transfer play an essential role in the fluid flow problem. The modern system of electric power generation, plasma, space vehicles, astrophysical flows, and cooling of nuclear reactors are governed by applications of thermal radiation and heat transfer of fluid flow. Because of its importance, the impact of thermal radiation, thermal slip, and heat transfer on MHD stagnation-point flow for different geometrical configurations has been tested by numerous researchers [13–19]. Makinde [20] investigated the hydromagnetic varied convection Hiemenz flow towards upright plate entrenched in an extremely permeable medium with radiation and internal heat generation. He discussed some aspects of the pertinent parameters on the dimensionless axial velocity, temperature and the concentration profiles, local skin friction, and local Nusselt number. Wubishet [21] examined the effect of induced magnetic field on MHD Hiemenz flow and heat transfer due to upper-convicted Maxwell fluid over a stretching sheet in the company of nanoparticles which is heated convectively. Numerical results are obtained for velocity, temperature, concentration profiles, skin friction coefficient, local Nusselt number, and Sherwood number. His result indicated that skin friction coefficient, the local Nusselt number, and Sherwood number increase with an increase in velocity ratio and Deborah number and decrease as the values of magnetic field increase.An analysis for the magnetohydrodynamic (MHD) axisymmetric Hiemenz flow and heat transmission over a shrinking sheet was presented by Mahapatra and Nandy [22], in the paper, the temperature profiles are obtained for the sheet with prescribed surface temperature and the sheet with prescribed surface heat flux. The effects of various physical parameters on Nusselt number and skin friction coefficient on the flow and heat transfer characteristics are discussed. Khan et al. [23] have investigated the effects of homogeneous-heterogeneous reactions in MHD flow of Casson fluid flow over a stretched surface and they discussed the behavior of velocity, temperature, concentration, drag force, and heat transfer rate. They found that, with larger magnetic parameter, the velocity field decreases, whereas surface drag force enhances. Temperature profile and heat transfer rate decay for higher estimation of Prandtl number. Shateyi and Fazle [24] analysed the problem of MHD mixed-convection stagnation-point flow towards a nonlinearly stretching vertical sheet in the presence of thermal radiation and viscous dissipation. It is noticed that the thermal boundary layer thicknesses and velocity are decreasing against the value of the nonlinearity parameter. The effect of injection on the MHD mixed-slip flow is to enhance the velocity field.Khan and Faris [25] developed Cattaneo-Christov double-diffusion model with nanofluid for Williamson fluid and constructed convergent solution for the model. In this research work, they discussed the physical features of MHD stagnation-point flow of Williamson nanomaterial over a stretched surface. Their results showed velocity boosts up significantly via velocity ratio parameter and variation of Weissenberg number yields improvement of temperature. The effect of the magnetic field, the radiation parameter, the uniform suction/injection parameter, and Hall parameter on MHD flow of a viscous, Newtonian, and electrically conducting fluid past a porous rotating infinite disk in consequence of Hall effect was studied by Devi et al. [26]. A comprehensive study on entropy generation in nanofluid flow towards a curved stretched surface with combined effects of activation energy, Brownian motion, viscous dissipation, nonlinear mixed convection, MHD, Joule heating, and thermophoresis diffusion is presented by Khan and Faris [27]. It is reported that velocity field is an increasing function of curvature parameter, while opposite impact is observed for Deborah number and velocity slip parameter. Nusselt number is increased via larger Eckert number and declines against thermophoretic parameter.Khan and Faris [28] presented mathematical modeling and numerical simulation for the steady, incompressible two-dimensional Darcy-Forchheimer nanofluid flow of viscous material towards a stretched surface and discussed the behavior of pertinent flow parameters on the velocity, temperature, Bejan number, and entropy generation for both nanoparticles for both nanoparticles (silicon dioxide and molybdenum). The effect of entropy generation in flow of viscous fluid of hybrid nanoparticles over a stretchable rotating disk subject to nonlinear thermal radiation and slip conditions is analysed by Khan [29]. The impacts of mixed-convection parameter, porosity parameter, velocities as well as thermal slips parameters, stretching parameter, nonlinear radiation parameter, and Reynolds number on radial velocity profile, tangential velocity profile, and temperature profile are studied.Fluids showing slip are important in the areas of technology and industry such as in the polishing of artificial heart valves and internal cavities. Fang et al. [30] analytically solved the MHD flow under slip conditions over a permeable shrinking surface, and they reported that the velocity slip at the shrinking surface greatly affects the velocity distribution and drag forces on the wall. El-Aziz and Ahmed [31] numerically analysed the influences of slip velocity and induced magnetic field on MHD stagnation-point flow and heat transfer of Casson fluid over a stretching sheet. Their result revealed that, with an increase in slip parameter, the velocity distribution reduces, whereas the temperature profile enhances.However, to the best of the authors’ knowledge, no study has been previously reported on the problem of MHD Hiemenz boundary layer flow at stagnation region against flat plate through a holey medium with thermal radiation, induced magnetic field, and velocity and thermal slip effect. In view of this and importance of MHD stagnation-point flow in engineering and in various technological applications and in numerous production processes such as the aerodynamic squeezing of polymers, rolling at high temperature, cooling control technology, and glass fibber production, the authors in the present paper aim to examine the effect of thermal radiation, magnetic field, Eckert number, Prandtl number, penetrability parameter, velocity, and thermal slip parameters on MHD Hiemenz flow at stagnation region against flat plate over a porous medium with thermal radiation and slip effect. The governing equations and their associated boundary conditions are at first changed into dimensionless form with the help of similarity variables; then the resulting system of nonlinear high-order ordinary differential equations is solved by the optimal homotopy asymptotic method (OHAM). The results obtained are then compared with those from the available literature for some particular values of the physical parameters, and it is found that they are in an excellent agreement. Graphical results for the velocity and temperature of the flow are discussed.


## 2. Basic Principles of Optimal Homotopy Asymptotic Method

The OHAM is an amendment of the homotopy asymptotic method (HAM), which is based on reducing the residual error. In this proposed method, the rheostat and adjustment of the confluence region are set properly. To simply illustrate the underlying rules of this method as expanded by Marinca and Herisanu [32] and different scholars, we take the nonlinear DE expressed as(1)Lvx+rx+Nvx=0,x∈Ω,with boundary condition:(2)Bx,dvdx=0,where Ω is problem domain, L and N represent the linear and nonlinear operators, vx is an unknown function, and rx is a known function. The optimal homotopy equation (also called deformation equation) is constructed as(3)1−pLϕx;p+ax=hpLϕx;p+ax+Nϕx;p+ax,where 0≤p≤1 is an inserting parameter.hp=∑k=1mpkck is an auxiliary function that determines convergence of a problem’s solution. The auxiliary function hp helps to adjust the convergence domain in addition to controlling a convergence region. By the suggested method, we will achieve approximate solution in a series form whenever the function ϕx;p,cj is expressed in Taylor’s series about p.(4)ϕx;p,cj=v0x+∑k=1∞vkx,cjpk,j=1,2,3,….As it is identified by many scholars, the convergence of equation (4) depends upon the values of convergence parameters cjj=1,2,3,…,m. If the series in equation (4) is convergent, then(5)v˜=v0x+∑k=1mvkx;cj.Plugging equation (5) in equation (1) results in the following residual:(6)Rx;cj=Lv˜x;cj+rx+Nv˜x;cj.If the function in equation (6) will vanish or equals zero, then v˜ will be an accurate solution of the given problem but this is not usually the case in nonlinear problems. The values of cjj=1,2,3,…,m can be optimally determined via various methods like colocation method, least-square method, Galerkin’s method, Ritz method, the Kantorovich method, and so forth, which can be utilized. Finally, substituting the values of these convergence control parameters in equation (5), one can get the analytical approximate solution of a problem. In our case, we applied the least-square method, because it is a powerful method to find optimal values of convergence control parameter presented in a solution of a problem in using different numerical and analytical methods [33].


## 3. Mathematical Formulation

In developing a model that describes the flow under consideration, we consider a two-dimensional steady laminar enforced convection magnetohydrodynamic Hiemenz flow at the stagnation area against a smooth plate through a holey medium with radiant heat, where a plate is placed aligned to thex-axis and the y-axis is ordinary to the sheet. The liquid in flow is assumed to be electrically conducted and incompressible. A continuous magnetic field of strength B0 is applied along the y direction, the wall temperature has been taken as variable, and the surface mass fluctuation is supposed to be uniform. The induced magnetic field, the Hall effects, and the viscus dissipation terms are neglected. The schematic diagram representing the physical phenomena of the problem is illustrated in Figure 1. Following the Raptis and Takhar [34] model designed for the permeable medium and by prefacing the boundary layer approximation, the governing equations describing the flow are as follows:Figure 1
Geometry of the flow problem.Continuity:(7)∂u∂x+∂v∂y=0.Momentum:(8)u∂u∂x+v∂u∂y=−1ρ∂P∂x+υ∂2u∂y2−υKu−σB02ρu.Energy:(9)u∂T∂x+v∂T∂y=α∂2T∂y2−1ρcp∂qr∂y+υ1ρcp∂u∂y2,where u and v are, respectively, the velocity constituents on the way to the x- and y-directions, P is the pressure, ρ represents density, K is porosity parameter, μ is the fluid viscosity coefficient, υ=μ/ρ is the kinematic viscosity, σ is the electrical conductivity, T is a fluid temperature and a permeable medium which are in local thermal equilibrium, α=k/ρcp is the equivalent thermal diffusivity, where k is a heat conduction coefficient, cp is a specific heat at unchanged pressure, and qr represents the radiative heat flux.The boundary conditions of the problems are as follows:(10)v=vw,u=uw+uslip=0+Nυ∂u∂y,T=Tw+Tslip=T∞+Axλ+D∂T∂y,,aty=0,(11)u=ul=ax,T=T∞,aty⟶∞,where vw is an invariable surface mass flux, which is positive for blowing and negative for suction; uslip is velocity slip, which is comparative with a shear stress in the layer of fluid next to wall of a plate and is given by Nυ∂u/∂y, where N is Navier's constant slip length; λ is the exponent of a wall temperature; D is temperature slip factor, and ul=mx is the unrestricted flow velocity where m denotes every positive number.Using the equationul=mx in (8), the equation becomes reduced to(12)ulduldx=−1ρ∂P∂x−υKul−σB02ρul.From equation (12), we have(13)−1ρ∂P∂x=ulduldx+υKul+σB02ρul.Upon using equation (13) in equation (8), we obtain(14)u∂u∂x+v∂u∂y=ulduldx+υKul−u+σB02ρul−u+υ∂2u∂y2.From (14), it is obvious that ul affects the fluid flow. Thus, convection of heat will be affected considerably, and this shows that the flow is a type of forced convection.Using Rosseland approximation, we can define or express the heat energy transmitted by electromagnetic surfs in contrast to heat transmitted by conduction in terms of temperature as follows:(15)qr=−4σ′3k1∂T4∂y,where σ′ represents the Stefan-Boltzmann constant and k1 represents the mean immersion coefficient. Assuming that the temperature variances in the interior of the flow are appropriately small and expanding T4 about the unrestricted stream temperature T∞ and disregarding higher-order expressions, we get(16)T4≈4T∞3−3T∞4.Upon using equation (16) in equation (15), we get(17)∂qr∂y=−16σ′T∞33k1∂2T∂y2.Substituting equation (17) into equation (9), we obtain(18)u∂T∂x+v∂T∂y=α∂2T∂y2+16σ′T∞33ρcpk1∂2T∂y2+υ1ρcp∂u∂y2,where α=k/ρcp is the thermal diffusivity, and this s equation illustrates that radiation enhances thermal diffusivity.Let us define the stream function satisfying equation (7) as follows:(19)u=∂ϕ∂y,v=−∂ϕ∂x.To solve equations (8) and (9), the following similarity variables have been introduced:(20)η=yaα,ϕ=aαxfη,θη=T−T∞Tw−T∞.From equations (19) and (20), we get(21)u=axf′η,v=−aαfη,where f′η=df/dη, Tw denotes a constant surface temperature, and θη is a nondimensional form of the temperature.Using equation (21) in equation (14), we obtain(22)Prf‴+1−f′2+ff″+Ω+M21−f′=0.We have the following boundary conditions:(23)f0=s,f′0=βf″0,aty=0,f′η=1,asη⟶∞,where Pr=υ/α is the Prandtl number, Ω=υ/ka is the permeability parameter, M=σB02/ρa is the charismatic parameter, s=−vw/aα is the mass transfer parameter, and β=Nυa/α is the dimensionless velocity slip parameter.From equations (10), (11), and (20), we have(24)∂T∂y=θ′Tw−T∞aα,∂2T∂y2=aαθ″Tw−T∞,∂T∂x=Aλθxλ−1.Plugging equation (24) into equation (18), we obtain(25)1+4R3θ″+fθ′−λf′θ+PrEcf″2=0.The above equation is subject to the following boundary conditions:(26)θ0=1+γθ′0,θη=0,asη⟶∞,where γ=Da/α is dimensionless thermal slip parameter, R=4σ•T∞3/k1k is radiation parameter, and Ec=ul2/cpTw−T∞ is the Eckert number.


## 4. Solution of the Problem Using the OHAM

Applying OHAM on the nonlinear ODEs equations (22) and (25), we construct the following equations:(27)1−pf‴−H1pPrf‴+1−f′2+ff″+Ω+M21−f′=0=01−pθ″−H2p1+4R3θ″+fθ′−λf′θ+PrEcf″2=0,where p∈0,1 and f,θ,H1p,H2p are expanded as follows:(28)f=f0+pf1+p2f2+p3f3,θ=θ0+pθ1+p2θ2+p3θ3,H1p=pc1+p2c2+p3c3,H2p=pc4+p2c5+p3c6.Using equation (28) in equation (27) and simplifying, rearranging, and collecting terms with common powers of p, we get the following zeroth-, first-, second-, and third-order problems together with their boundary conditions.The zeroth-order problem is(29)f‴η=0,θ″η=0,with boundary conditions:(30)f00=s,f0′0=βf0″0,f0′η=1,asη⟶∞,θ00=1+γθ0′0,θ0η=0,asη⟶∞.The solution for equation (29) with BCs (30) is(31)f0η=1201−βf0″0η2+βf0″0η+s,θ0η=1+γθ0′01−η10.The first-order problem is(32)f1″′η,c1=c1Pr1−f0′2+f0f0″+Ω+M21−f0′,θ1″η,c4=3c43+4Rf0θ0′−λf0′θ0+PrEcf0″2,with boundary conditions:(33)f10=0,f1′0=0,θ10=0,f1′η=0,θ1η=0,asη⟶∞.Its solution is(34)f1η,c1=5c1Pr1+Ω+M21−βf0″0+s10−βf0″0η330−η22+16βf0″0βf0″0−1+Ω+M2f0″0−110η4200−η2+16βf0″01−βf0″02−12η55000−η22,θ1η,c4=−3c41+γθ0′06+8Rs+βf0″0η2−10η+13βf0″0+λ1−2βf0″0η310−10η+1121−2λ1−βf0″0η4100−10η+3c4PrEc6+8R1+βf0″0βf0″0−2η2200−η10.The second-order problem is(35)f2″′η,c1,c2=c1+c2f1″′+c2Pr1−f0′2+f0f0″+Ω+M21−f0′+f0f1″+f1f0″−2f0′f1′−f1,θ2″η,c4,c5=1+c4θ1″+3c43+4Rf1θ0′+f0θ1′−λf1′θ0+f0′θ1+2PrEcf0″f1″+3c53+4Rf0θ0′−λf0′θ0+PrEcf0″2.The above problem is subject to the following boundary conditions:(36)f20=0,f2′0=0,θ20=0,f2′0=0,θ20=0,asη⟶∞.The solution for this second-order problem is(37)f2η,c1,c2=10Pr1+Ω+M21−βf0″0+s10−βf0″0s10+βf0″0c1+c12η36−η24+c12Prη4240−η312+5η212+1−βf0″0−η614400+η5800+25η224−2βf0″0+1η51200−η4480+25η212βf0″0βf0″0−1+Ω+M2βf0″0−110c1+c12η42400−η212+c12Prη51200−250η29+50η219+c12Pr1−βf0″0−η7720000+η52400−65η272−2βf0″0+1η67200−η41440+35η236+10Prβf0″01−βf0″02−12c1+c12+c2η560000−η224+c12Prη6360000−η372+η28+1−βf0″0η51200−η4480+25η212+η8201600000−η514400+310η22016+c12Prβf0″0βf0″0−2βf0″0+12η7520000−η4288+5η29+1−βf0″0η840320−η57200+125η2504+1+Ω+M2η360−η24+10βf0″0βf0″0−1+Ω+M2βf0″0−1η42400−η212+s−βf0″0s+10βf0″0−10βf0″0Ω+M2η3600−η240,θ2η,c4,c5=−3c41+γθ0′3+4R3c4s+100βf0″06+8R1+c43+4R3c4η2100−1+1210−βf0″0η510000−η41200+53+sη3300−η220+53+λ10−βf0″0η41200−η520000−103+βf0″0η4600−η360+λη36−η41200−253c4βf0″0+λ1−2βf0″06+8R1+c43+4R3c4η330−1003+1210−βf0″0η63000−η436+5009+βf0″0η5200−5η39+5009+sη412000−5η23+2503+λ1−βf0″0η436−η69000−5003+λ5η39−η5600−350093c41−βf0″0+2λβf0″0−13+4R1+c43+4R3c4η4600−503+13610−βf0″0η63000−η436+5009+βf0″0η73500−η44+25007+η64500−5η318+s6η5500−5η2+λ181−βf0″0η44−η714000−125007+λβf0″065η33−η63000−40003−3c41+γθ0′3+4Rc1Pr1+Ω+M21−βf0″0+s10−βf0″010s+10βf0″0η31200−η448+125+λη412−5η36−η5400+125+βf0″02−βf0″0+Ω+M2βf0″0−110η67200−η4144+5009+λη51200−5η318+η472−η618000+10009+βf0″0−12βf0″02−12η72520000−η4288+λη636000−5η336+η4144−η750400+3c4PrEc1−2βf0″0+βf0″026+8R1+c13+4R3c4η2100−1+5−βf0″02η5100000−η412000−16+βf0″060η610−η3+sη3300−η220+53+λ1−βf0″0η41200−η520000−103+λβf0″0η360−η41200−253−3c51+γθ0′3+4R12−12βf0″0+λβf0″0−1η41200−253−βf0″0+λ1−2βf0″0η360−503+s+10λβf0″02η210−10+3c5PrEc3+4Rβf0″02−2βf0″0+1η210−5η.The third-order problem is(38)f3″′η,c1,c2,c31+c1f2″′+c2f1″′+c3Prf0f0″−f0′2+Ω+M21−f0′+c2Prf0f1″+f1f0″−2f0′f1′−f1'Ω+M2−c1Pr2f0′f2′+f12+f2′Ω+M2,θ3″η,c4,c5,c6=1+c4θ2″+c5θ1″+3c43+4Rf0θ2′+f1θ1′+f2θ0′−λf0′θ2+f1′θ1+f2′θ0+PrEc2f0″f2″+f1'2+3c53+4Rf1θ0′+f0θ1′−λf1′θ0+f0′θ1+PrEc2f0″f1″−3c6λ3+4Rf0′θ0,with boundary conditions:(39)f30=0,f3′0=0,θ30=0,f3′η=0,θ3η=0,asη⟶∞.The solution for equation (38) subject to the boundary conditions in equation (39) is obtained in a manner similar to that in equations (29), (32), and (35). However, its expression is very large to be included in this paper.Collecting the solution components of equations (31) , (34), and (37) and the solution of equation (38), the four-term solution obtained by OHAM for p = 1 is(40)f˜η,c1,c2,c3=f0η+f1η,c1+f2η,c1,c2+f3η,c1,c2,c3,θ˜η,c4,c5,c6=θ0η+θ1η,c4+θ2η,c4,c5+θ3η,c4,c5,c6.To find the value ofci′s,i=1,2,…,6 in equation (40) for Pr=0.72, β=γ=0.2, R=Ω=M=s=λ=1,and Ec=0.1, we applied the least-square method to decide the values of convergence control parameters and we get c1=−0.0137964;c2=4.181133;c3=−0.006856;c4=−0.04417926;c5=−9.05104899;c6=2.236272.By substituting these values into equation (40), we obtain the solution for the problem.


## 5. Results and Discussion

In this work, we employed OHAM to develop the convergent solutions for nonlinear system equations (22) and (25) subject to the BCs equations (23) and (26). Analyses of the effects of various embedded parameters in the flow model on temperature, momentum of the fluid, skin friction coefficient, and local Nusselt number have been made, and the fallouts of the investigative computations are presented via tables and different figures and are discussed.


### 5.1. Comparative Analysis

To confirm the accuracy and appropriateness of the method, we have compared the value of the skin friction coefficientf″0 and the temperature ramp at the wall −θ′0 obtained by the present method with the results reported by previous studies under limiting conditions [35–37]. The comparisons of the obtained outcome are presented in Tables 1 and 2 , illustrating an excellent agreement, and thus we are confident that the method is suitable for investigation of the problem. Table 1 is prepared to show an association of the discrepancy of the friction coefficient between a fluid and its enclosing surface for some values of suction/blowing and magnetic meadow parameters with previous studies. From the table, it is observed that our result is in excellent agreement with the results reported previously by Bana et al. [37], Bhatti et al. [35], and Kechil and Hashim [36]. The table also reveals that local skin friction coefficient f″0 increases with growing values of s and M.Table 1
Comparison of the values of local skin friction coefficientf″0 obtained by the OHAM and other numerical methods in the existing literature for altered values of M and s whenPr=1,β=Ω=0.
sMPresent resultBana et al. [37]Bhatti et al. [35]Kechil and Hashim [36]−100.7565680.7565750.75650.755979466211.1164301.1164211.11641.116429235021.8776231.8776211.87761.877622149254.6675234.667525—4.6675255187109.5859159.585914—9.5859131993001.2325641.2325861.23251.23065249311.5853121.5853311.58531.585203958522.3466402.3466632.4662.346663511055.1479415.14793—5.14796460251010.0747310.074742—10.074741096101.8893221.8893121.889310.07474109612.2029172.2029332.20292.201433545822.9201132.9201102.920292.920114202855.6768315.676830—5.67683034211010.58836810.588367—10.5883674767Table 2
Comparison of values of the local Nusselt number−θ′0 obtained by the present method and other numerical methods in the existing literature for different s and M values whenPr=1,β=Ω=λ=0.
sMPresent resultBana et al. [37]Bhatti et al. [35]−100.116760.116750.1167710.140040.14000.1400020.173130.173120.17312000.570510.570460.5703510.595370.595350.5953920.634160.634130.63418101.323691.323681.3236811.338051.338041.3380421.364471.364461.36446In Table2, a relationship of the change in local Nusselt number on behalf of numerous values of suction/blowing parameter and magnetic flux parameter with previously reported data in the cited literature was presented. As is illustrated in Table 2, the results are in good agreement with the results reported by Bana et al. [37] and Bhatti et al. [35]. From the table, it is also obvious that the Nusselt number rises with the charismatic field parameter (M) value.


### 5.2. Physical Description

The optimal homotopy analysis method described in Section2 of this paper is employed to find the solution for nonlinear differential system. Here, we will discuss the impact of relevant parameters in the flow model on velocity and temperature distribution.


### 5.3. Velocity

The influence of induced magnetic field, permeability parameter, velocity slip parameter, suction/blowing parameter, and Eckert number on velocity profile is illustrated in Figures2–6 .Figure 2
Velocity profile for numerous values ofM when Pr=0.72, Ω=λ=R=s=1, β=0.2, γ=0.2, and Ec=0.1.Figure 3
Velocity profile for various values ofΩ when Pr=0.72, M=λ=R=s=1, β=0.2, γ=0.2, and Ec=0.1.Figure 4
Velocity profile for altered values ofβ when Pr=0.72, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure 5
Velocity sketch for diverse values ofs when Pr=0.72, Ω=M=λ=R=1, β=0.2, γ=0.2, and Ec=0.1.Figure 6
Velocity graph for altered values of Ec whenPr=0.72, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure2 illustrates the magnetic field effects on the fluid velocity distribution. It is detected from the figure that boundary layer thickness and the fluid velocity profile decline against the charismatic field M. This result is manifested from reality, because the presence of a magnetic field produces a body force, called Lorentz force, which opposes the fluid motion and, hence, it reduces the fluid freedom of movement. Therefore, when the magnetic flux increases, the retardation force also increases, and this resistance offered to the flow is responsible for decreasing the fluid velocity. Figure 3 demonstrates the effect of penetrability parameter on the fluid velocity dispersal. It has been realized that the rising of the penetrability parameter value results in increasing the velocity profile, which causes an increment in the width of the velocity boundary wall layer. Physically, it is clear that porous medium reduces the motion of the fluid. Thus, as the permeability parameter Ω increases resistance to the motion of the fluid along the surface drops, in other words, when the value of porosity parameter K declines, the velocity boosts with the increase of the value of Ω.Figure4 exhibits the graphical illustration of the impact of velocity slip parameter β on momentum (velocity) profiles. From the figure, it is seen that intensification in slip parameter results in deterioration of velocity profile. This is due to the fact that an upsurge in the velocity slip factor produces friction force, which allows more fluid to slip past the sheet, and consequently the motion of the flow decelerates. As it is exhibited by Figure 5, the momentum profile was raised with suction/blowing parameter s. In fact, the action of suction increases the momentum boundary layer thickness, which results in accelerating the motion of the fluid. The alteration of value of the ratio of the advective mass transfer to the heat dissipation potential Ec does not affect velocity of the fluid as is displayed in Figure 6. Likewise, altering radiation parameter, wall temperature exponent, and thermal slip values does not disturb the momentum profile of the fluid in motion; this is evident from the fact that the velocity equation is not being coupled with the energy equation of the problem.


### 5.4. Temperature

The impact of Prandtl number, suction/blowing parameter, permeability parameter, velocity slip parameter, radiation, Eckert number, wall temperature parameter, thermal slip parameter, and magnetic field parameter on temperature distribution is displayed in Figures7–14 . Figure 7 shows that the fluid temperature is reducing contrary to Prandtl number Pr. This implies that the thermal boundary layer viscosity diminishes in opposition to the value of Pr. This happens basically because high value of Prandtl number relates to weak thermal diffusivity, which upshots in a thinner thermal boundary layer. As a matter of fact, intensification in the Prandtl number infers an upsurge in fluid viscosity, which consecutively causes deterioration in temperature distribution. Hence, Pr can be used to upturn the rate of cooling in flows. Figure 8 displays the alteration in temperature pertaining to the suction/blowing parameter s. It is noticed from the figure that an upsurge in the suction/blowing parameter value results in decreasing the heat transmission proportion from the wall. For the reason that, with strong blowing, the thermal boundary thickness increases and accordingly the heated fluid moves farther from the wall and forms an insulating layer of nearly the equivalent temperature as that of the wall, suction can be used as a means for freezing a surface.Figure 7
Temperature distribution for numerous values of Pr whenΩ=M=λ=R=s=1, γ=0.2,β=0.2, and Ec=0.1.Figure 8
Variation of temperature against various values ofs when Pr=0.72,Ω=M=λ=R=1, γ=0.2, β=0.2, and Ec=0.1.Figure 9
Variation in temperature against various values ofΩ when Pr=0.72,M=λ=R=s=1, γ=0.2, β=0.2, and Ec=0.1.Figure 10
Temperature alteration for numerous values ofβ when Pr=0.72, γ=0.2, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure 11
Alteration in temperature for various values ofR when Pr=0.72, β=0.2, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure 12
Temperature profile for altered values of Ec whenPr=0.72, β=0.2, Ω=M=λ=R=s=1, and γ=0.2.Figure 13
Temperature sketch for altered values ofλ when Pr=0.72, Ω=M=R=s=1, β=0.2, γ=0.2, and Ec=0.1.Figure 14
Temperature sketch for various values ofγ on when Pr=0.72, Ω=M=λ=R=s=1, β=0.2, and Ec=0.1.Figure9 describes the impact of penetrability of the medium on temperature, and it reveals that the thermal boundary layer thickness decreases on account of an increase in the value of penetrability parameter Ω. Physically, permeability of porous medium increases means porosity of the medium, opposing the fluid motion, which decreases. This reduction in resistance to the flow is responsible for offering a decrease in the temperature. Figure 10 represents the temperature sketch on behalf of numerous input values of velocity slip parameter β, and it reveals that the thermal boundary layer thicknesses declines with increased value of β. Therefore, stronger unrestricted velocity results in reduction of the temperature and the width of thermal boundary layer. In Figure 11, the effect of radiation R on temperature of the fluid was illustrated and the figure displayed that the thermal boundary layer width and the temperature dissemination rise as R increases. This happens since the divergence of the radiative heat change increases with the decreasing value of the Rosseland radiative absorptivity. This, in turn, leads to an upsurge in the rate of radiative heat transmission to the fluid, which in turn enhances the fluid temperature.The distinction in the fluid temperature owing to the alteration of Eckert number Ec has been illustrated in Figure12. The figure shows that, with increased value of Ec, the fluid temperature intensifies. This happens because heat is produced in the fluid as the value of Ec increases due to frictional heating. Physically, Eckert number is the ratio of kinetic energy to the specific enthalpy difference between wall and fluid. Therefore, an increase in Eckert number causes the transformation of kinetic energy into internal energy by work that is done against the viscous fluid stresses. Due to this, increasing Ec enhances the temperature of the fluid. In Figure 13, the impact of wall temperature exponent λ on the layer of a liquid between the free stream and a heat-exchange surface is displayed, and it has been detected that the fluid temperature decreases with the increasing value of λ.The effect of thermal slipγ on temperature of the fluid was displayed in Figure 14. It is detected that the temperature dispersion within the fluid in flow is decreased when the parameter that denotes thermal slip gradually upsurges in value. Physically, this is evident from the point that, following in increasing thermal slip, less heat is shifted from the boundary to the fluid, which is responsible for the temperature distribution reduction. In Figure 15, the influence of magnetic field M on the temperature profile of the fluid has been displayed. From the figure, it is noticeable that the temperature boosts with rise in the magnetic field parameter M. Physically, when M is increased, an opposing force (Lorentz force) that retards the flow motion will be produced in the fluid. In fact, the presence of magnetic field decelerates velocity and consecutively persuades the temperature field that increases the temperature profiles.Figure 15
Temperature sketch for altered values ofM when Pr=0.72, γ=0.2, Ω=λ=R=s=1, β=0.2, and Ec=0.1.


## 5.1. Comparative Analysis

To confirm the accuracy and appropriateness of the method, we have compared the value of the skin friction coefficientf″0 and the temperature ramp at the wall −θ′0 obtained by the present method with the results reported by previous studies under limiting conditions [35–37]. The comparisons of the obtained outcome are presented in Tables 1 and 2 , illustrating an excellent agreement, and thus we are confident that the method is suitable for investigation of the problem. Table 1 is prepared to show an association of the discrepancy of the friction coefficient between a fluid and its enclosing surface for some values of suction/blowing and magnetic meadow parameters with previous studies. From the table, it is observed that our result is in excellent agreement with the results reported previously by Bana et al. [37], Bhatti et al. [35], and Kechil and Hashim [36]. The table also reveals that local skin friction coefficient f″0 increases with growing values of s and M.Table 1
Comparison of the values of local skin friction coefficientf″0 obtained by the OHAM and other numerical methods in the existing literature for altered values of M and s whenPr=1,β=Ω=0.
sMPresent resultBana et al. [37]Bhatti et al. [35]Kechil and Hashim [36]−100.7565680.7565750.75650.755979466211.1164301.1164211.11641.116429235021.8776231.8776211.87761.877622149254.6675234.667525—4.6675255187109.5859159.585914—9.5859131993001.2325641.2325861.23251.23065249311.5853121.5853311.58531.585203958522.3466402.3466632.4662.346663511055.1479415.14793—5.14796460251010.0747310.074742—10.074741096101.8893221.8893121.889310.07474109612.2029172.2029332.20292.201433545822.9201132.9201102.920292.920114202855.6768315.676830—5.67683034211010.58836810.588367—10.5883674767Table 2
Comparison of values of the local Nusselt number−θ′0 obtained by the present method and other numerical methods in the existing literature for different s and M values whenPr=1,β=Ω=λ=0.
sMPresent resultBana et al. [37]Bhatti et al. [35]−100.116760.116750.1167710.140040.14000.1400020.173130.173120.17312000.570510.570460.5703510.595370.595350.5953920.634160.634130.63418101.323691.323681.3236811.338051.338041.3380421.364471.364461.36446In Table2, a relationship of the change in local Nusselt number on behalf of numerous values of suction/blowing parameter and magnetic flux parameter with previously reported data in the cited literature was presented. As is illustrated in Table 2, the results are in good agreement with the results reported by Bana et al. [37] and Bhatti et al. [35]. From the table, it is also obvious that the Nusselt number rises with the charismatic field parameter (M) value.


## 5.2. Physical Description

The optimal homotopy analysis method described in Section2 of this paper is employed to find the solution for nonlinear differential system. Here, we will discuss the impact of relevant parameters in the flow model on velocity and temperature distribution.


## 5.3. Velocity

The influence of induced magnetic field, permeability parameter, velocity slip parameter, suction/blowing parameter, and Eckert number on velocity profile is illustrated in Figures2–6 .Figure 2
Velocity profile for numerous values ofM when Pr=0.72, Ω=λ=R=s=1, β=0.2, γ=0.2, and Ec=0.1.Figure 3
Velocity profile for various values ofΩ when Pr=0.72, M=λ=R=s=1, β=0.2, γ=0.2, and Ec=0.1.Figure 4
Velocity profile for altered values ofβ when Pr=0.72, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure 5
Velocity sketch for diverse values ofs when Pr=0.72, Ω=M=λ=R=1, β=0.2, γ=0.2, and Ec=0.1.Figure 6
Velocity graph for altered values of Ec whenPr=0.72, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure2 illustrates the magnetic field effects on the fluid velocity distribution. It is detected from the figure that boundary layer thickness and the fluid velocity profile decline against the charismatic field M. This result is manifested from reality, because the presence of a magnetic field produces a body force, called Lorentz force, which opposes the fluid motion and, hence, it reduces the fluid freedom of movement. Therefore, when the magnetic flux increases, the retardation force also increases, and this resistance offered to the flow is responsible for decreasing the fluid velocity. Figure 3 demonstrates the effect of penetrability parameter on the fluid velocity dispersal. It has been realized that the rising of the penetrability parameter value results in increasing the velocity profile, which causes an increment in the width of the velocity boundary wall layer. Physically, it is clear that porous medium reduces the motion of the fluid. Thus, as the permeability parameter Ω increases resistance to the motion of the fluid along the surface drops, in other words, when the value of porosity parameter K declines, the velocity boosts with the increase of the value of Ω.Figure4 exhibits the graphical illustration of the impact of velocity slip parameter β on momentum (velocity) profiles. From the figure, it is seen that intensification in slip parameter results in deterioration of velocity profile. This is due to the fact that an upsurge in the velocity slip factor produces friction force, which allows more fluid to slip past the sheet, and consequently the motion of the flow decelerates. As it is exhibited by Figure 5, the momentum profile was raised with suction/blowing parameter s. In fact, the action of suction increases the momentum boundary layer thickness, which results in accelerating the motion of the fluid. The alteration of value of the ratio of the advective mass transfer to the heat dissipation potential Ec does not affect velocity of the fluid as is displayed in Figure 6. Likewise, altering radiation parameter, wall temperature exponent, and thermal slip values does not disturb the momentum profile of the fluid in motion; this is evident from the fact that the velocity equation is not being coupled with the energy equation of the problem.


## 5.4. Temperature

The impact of Prandtl number, suction/blowing parameter, permeability parameter, velocity slip parameter, radiation, Eckert number, wall temperature parameter, thermal slip parameter, and magnetic field parameter on temperature distribution is displayed in Figures7–14 . Figure 7 shows that the fluid temperature is reducing contrary to Prandtl number Pr. This implies that the thermal boundary layer viscosity diminishes in opposition to the value of Pr. This happens basically because high value of Prandtl number relates to weak thermal diffusivity, which upshots in a thinner thermal boundary layer. As a matter of fact, intensification in the Prandtl number infers an upsurge in fluid viscosity, which consecutively causes deterioration in temperature distribution. Hence, Pr can be used to upturn the rate of cooling in flows. Figure 8 displays the alteration in temperature pertaining to the suction/blowing parameter s. It is noticed from the figure that an upsurge in the suction/blowing parameter value results in decreasing the heat transmission proportion from the wall. For the reason that, with strong blowing, the thermal boundary thickness increases and accordingly the heated fluid moves farther from the wall and forms an insulating layer of nearly the equivalent temperature as that of the wall, suction can be used as a means for freezing a surface.Figure 7
Temperature distribution for numerous values of Pr whenΩ=M=λ=R=s=1, γ=0.2,β=0.2, and Ec=0.1.Figure 8
Variation of temperature against various values ofs when Pr=0.72,Ω=M=λ=R=1, γ=0.2, β=0.2, and Ec=0.1.Figure 9
Variation in temperature against various values ofΩ when Pr=0.72,M=λ=R=s=1, γ=0.2, β=0.2, and Ec=0.1.Figure 10
Temperature alteration for numerous values ofβ when Pr=0.72, γ=0.2, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure 11
Alteration in temperature for various values ofR when Pr=0.72, β=0.2, Ω=M=λ=R=s=1, γ=0.2, and Ec=0.1.Figure 12
Temperature profile for altered values of Ec whenPr=0.72, β=0.2, Ω=M=λ=R=s=1, and γ=0.2.Figure 13
Temperature sketch for altered values ofλ when Pr=0.72, Ω=M=R=s=1, β=0.2, γ=0.2, and Ec=0.1.Figure 14
Temperature sketch for various values ofγ on when Pr=0.72, Ω=M=λ=R=s=1, β=0.2, and Ec=0.1.Figure9 describes the impact of penetrability of the medium on temperature, and it reveals that the thermal boundary layer thickness decreases on account of an increase in the value of penetrability parameter Ω. Physically, permeability of porous medium increases means porosity of the medium, opposing the fluid motion, which decreases. This reduction in resistance to the flow is responsible for offering a decrease in the temperature. Figure 10 represents the temperature sketch on behalf of numerous input values of velocity slip parameter β, and it reveals that the thermal boundary layer thicknesses declines with increased value of β. Therefore, stronger unrestricted velocity results in reduction of the temperature and the width of thermal boundary layer. In Figure 11, the effect of radiation R on temperature of the fluid was illustrated and the figure displayed that the thermal boundary layer width and the temperature dissemination rise as R increases. This happens since the divergence of the radiative heat change increases with the decreasing value of the Rosseland radiative absorptivity. This, in turn, leads to an upsurge in the rate of radiative heat transmission to the fluid, which in turn enhances the fluid temperature.The distinction in the fluid temperature owing to the alteration of Eckert number Ec has been illustrated in Figure12. The figure shows that, with increased value of Ec, the fluid temperature intensifies. This happens because heat is produced in the fluid as the value of Ec increases due to frictional heating. Physically, Eckert number is the ratio of kinetic energy to the specific enthalpy difference between wall and fluid. Therefore, an increase in Eckert number causes the transformation of kinetic energy into internal energy by work that is done against the viscous fluid stresses. Due to this, increasing Ec enhances the temperature of the fluid. In Figure 13, the impact of wall temperature exponent λ on the layer of a liquid between the free stream and a heat-exchange surface is displayed, and it has been detected that the fluid temperature decreases with the increasing value of λ.The effect of thermal slipγ on temperature of the fluid was displayed in Figure 14. It is detected that the temperature dispersion within the fluid in flow is decreased when the parameter that denotes thermal slip gradually upsurges in value. Physically, this is evident from the point that, following in increasing thermal slip, less heat is shifted from the boundary to the fluid, which is responsible for the temperature distribution reduction. In Figure 15, the influence of magnetic field M on the temperature profile of the fluid has been displayed. From the figure, it is noticeable that the temperature boosts with rise in the magnetic field parameter M. Physically, when M is increased, an opposing force (Lorentz force) that retards the flow motion will be produced in the fluid. In fact, the presence of magnetic field decelerates velocity and consecutively persuades the temperature field that increases the temperature profiles.Figure 15
Temperature sketch for altered values ofM when Pr=0.72, γ=0.2, Ω=λ=R=s=1, β=0.2, and Ec=0.1.


## 6. Conclusion

In this work, magnetohydrodynamic Hiemenz boundary layer flow over a plate with the stagnation area considered against a smooth plate through a holey medium with thermal and velocity slip effects is considered. The liquid in flow is assumed to be electrically conducted and incompressible. A continuous magnetic field is applied along the vertical direction, the wall temperature is taken as variable, and the surface mass fluctuation is supposed to be uniform. The Hall effects and the internal friction terms are disused. The nonlinear equation describing the flow problem is converted to third-order nonlinear ODEs via similarity variable transformation system and then solved by the use of the OHAM. The effects of various parameters in the model of the problem on temperature and velocity are analysed. In order to confirm the applicability and accuracy of the method, the obtained outcomes are verified with some successful and available theoretical data in the literature [35–37] and the comparative outcomes showed an excellent agreement with the current computational analysis.From the findings of the study, the following conclusions were drawn:(1)
The thickness of the velocity boundary layer and the profile of the fluid velocity decreases asM (the magnetic field) increases; however, the temperature increased as the charismatic field parameter (M) increased in value.(2)
The fluid velocity rises with the increased values of both the penetrability parameter (Ω) and suction/blowing parameter (s).(3)
The flow movement, the thermal boundary layer, and gradient of the temperature reduce against the velocity slip parameter (β).(4)
The temperature of the fluid rises in conjunction with the Eckert number (Ec) but the velocity of the fluid is not influenced by varying the value of Ec.(5)
An upsurge in thermal slip parameter results in a discount in the thermal boundary layer thickness and so reduces the temperature profile.(6)
The thermal boundary layer thickness upturns following an increased value of both the magnetic field parameter and radiation parameter, but it declines in conjunction with an increased value of Prandtl number (Pr) and the exponent of wall temperature parameter (λ).(7)
The skin friction coefficient and the Nusselt number increase with rising value of suction/blowing and magnetic field parameter.



---

*Source: 5590657-2021-04-12.xml*