# Hall and Ion Slip Effects on Mixed Convection Flow of Eyring-Powell Nanofluid over a Stretching Surface

Authors: Wubshet Ibrahim; Temesgen Anbessa
Journal: Advances in Mathematical Physics
 (2020)

---

## Abstract

The purpose of this research is to inspect the mixed convection flow of Eyring-Powell nanofluid over a linearly stretching sheet through a porous medium with Cattaneo–Christov heat and mass flux model in the presence of Hall and ion slip, permeability, and Joule heating effects. Proper similarity transforms yield coupled nonlinear differential systems, which are solved using the spectral relaxation method (SRM). The story audits show that the present research problem has not been studied until this point. Efficiency of numerous parameters on velocity, temperature, and concentration curves is exposed graphically. Likewise, the numerical values of skin friction coefficients, local Nusselt, and Sherwood numbers are computed and tabulated for some physical parameters. It is manifested that fluid velocities, skin friction coefficients, local Nusselt, and Sherwood numbers promote with the larger values of Eyring-Powell fluid parameterε. It is also noticed that primary velocity promotes with larger values of mixed convection parameter λ, Hall parameter βe, and ion slip parameter βi, while the opposite condition is observed for secondary velocity, temperature, and concentration. Furthermore, comparative surveys between the previously distributed writing and the current information are made for explicit cases, which are examined to be in a marvelous understanding.


---

## Body

## 1. Introduction

Mixed convection streams emerge in numerous transport processes both in nature and in engineering applications. They play a vital role, for instance, in air limit layer streams, heat exchangers, atomic reactors, solar collectors, and in electronic hardware. Such procedures happen when the impacts of buoyancy forces in constrained convection or the impacts of the constrained stream in natural convection become noteworthy. The interaction of natural and constrained convection is particularly articulated in circumstances where the constrained stream velocity is low as well as the temperature contrasts are huge. This stream is likewise an important kind of stream showing up in numerous mechanical procedures, for example, assembling and extraction of polymer and elastic sheets, paper creation, wire drawing and glass-fiber creation, dissolve turning, and consistent throwing. Moreover, mixed convection in permeable media has numerous applications, for example, food handling and storage, metallurgy, geophysical framework, fibrous insulation, and underground removal of atomic waste. Then again, thermal conductivity of the ordinary heat transport liquids, for instance, water, oil, and ethylene glycol are exceptionally low. Thus, expanding thermal conductivity of the traditional liquids prompts improves the heat transport of these liquids. As of late, nanofluids are presented with upgraded thermal conductivity. A nanofluid is a suspension of nanoparticles with normal sizes beneath 100 nm in the base liquids. Due to the upgraded thermal conductivity, nanofluids are proposed for some mechanical applications, for example, transportation, atomic reactors, and nourishment. In this manner, numerous ongoing analysts have generous enthusiasm for the mixed convection stream of nanofluid over an extending sheet in a permeable medium due to their impressive use in the modern and innovative applications. Accordingly, Ameen et al. [1] investigated a 3D turning stream of carbon nanotubes (CNTs) over a permeable stretchable sheet for warmth and mass exchange with thought of kerosene oil as a base fluid and announced that both velocity fields delineate expanding conduct for enormous value of the CNT nanoparticle. Besides, Rasool et al. [2] examined a numerical investigation of the MHD Williamson nanofluid flow maintained to flow through permeable medium bounded by a nonlinearly stretching flat surface and reported that an enhancement in the values of magnetic parameter augments the drag force, whereas a decrement in the drag force can be observed as the values of Weissenberg number enhances. Also, Ibrahim and Anbessa [3] examined the mixed convection stream of a Maxwell nanofluid over a growing surface in a penetrable medium utilizing the spectral relaxation method (SRM) and found that a climb in Deborah number decreases both the stream and transverse velocity profiles, while the backwards design is seen with enlargement in the mixed convection parameter. Furthermore, Waini et al. [4] examined the consistent mixed convection stream along a vertical surface implanted in a permeable medium with hybrid nanoparticles utilizing the bvp4c solver in Matlab programming. They considered both assisting and opposing streams and found that there exist dual solutions for the case of opposing stream. Further, numerous other ongoing works have been done around there as found in references [4–12].Investigation of nonNewtonian liquid has incredible significance because of its numerous industrial and engineering applications. Specifically, these liquids are used in the material preparing, concoction and atomic ventures, bioengineering, oil repository designing, polymeric fluids, and groceries. A few liquids like paints, paper mash, shampoos, ketchup, fruit purée, slurries, certain oils, and polymer arrangements are instances of nonNewtonian liquids. NonNewtonian liquids are logically perplexed when compared with Newtonian liquids due to nonlinear association among the stress and strain rate. Various models have been proposed in the writing for the examination of nonNewtonian fluids; however, not a sole model is built up that shows all properties of nonNewtonian liquids. In general, nonNewtonian fluids are mainly characterized into three kinds: to be explicit differential, rate, and integral kinds. Among the nonNewtonian liquid models, the Eyring-Powell model has accomplished exceptional consideration of the analysts because of its distinctive attributes in current science. This could be derived from kinetic theory of liquids rather than from its empirical relation. Likewise, the Eyring-Powell model decreases to Newtonian liquid qualities for low and high shear rate. Warmth and mass exchange in the Eyring-Powell liquid model assumes a significant role in the procedures which include formation and spread of haze, plotting of concoction handling instrumentation, environmental contamination, drying of permeable slides, raised oil recuperation, warm protection, and underground energy transport. Khan et al. [13] explored blended convection stream of the Eyring-Powell liquid model with variable viscosity and convective limit conditions over a slanted extending sheet utilizing the homotopy analysis method (HAM) and revealed that velocity profile increments by expanding liquid parameter M and Grashof numbers Gm and Gr, while velocity profile diminishes by expanding variable viscosity parameter A, liquid parameter k, and magnetic field parameter Ha. Moreover, Khan et al. [14] examined warm dispersion and dissemination thermo impacts on the wobbly flow of electrically conducting Eyring-Powell liquid over an oscillatory extending sheet with convective limit conditions utilizing HAM and found that the bigger values of Eyring-Powell liquid parameter improve the amplitude of velocity and limit the layer thickness. However, inverse impacts are seen in temperature and concentration profiles. Furthermore, Ishaq et al. [15] explored two dimensional nanofluid film stream of Eyring-Powell liquid with variable warmth transmission in the presence of MHD on a shaky permeable extending sheet and announced that porosity parameter diminishes the movement of the fluid movies, and enlarging the nanoparticle concentration effectively expands the rubbing characteristic of Eyring-Powell nanofluid. In recent times, various researches have been made in the field of Eyring-Powell nanofluids that can be found in references [16–23].Heat transport in the presence of strong magnetic field is significant in different parts of MHD power generation, nanotechnological processing, nuclear energy systems exploiting fluid metals, and blood stream control. Further, Hall and ion slip impacts become noteworthy in strong magnetic fields and can impressively influence the current density in hydromagnetic heat transport. Most recently, Krishna and Chamkha [24] explored Hall and ion slip impacts on the MHD free convective pivoting stream of nanofluids in a permeable medium past a moving vertical semiboundless level plate and announced the velocity increments with Hall and ion slip parameters. Furthermore, Rani et al. [25] explored Hall and ion slip impacts on the MHD natural convective pivoting stream of Ag-water-based nanofluid past a semiboundless penetrable moving plate with consistent warmth source utilizing the perturbation method. They found that in the limit layer area, liquid velocity diminishes with the expanding values of rotation parameter, magnetic field parameter, and suction parameter, while it increments with the expanding values of Hall and ion slip parameters. Further examinations identified with the MHD stream issues with Hall and ion slip impacts can be found in references [26–31].The temperature differences between two distinct bodies cause heat transport device. It assumes a significant task in the creation of energy, cooling of atomic reactors, and biomedical applications such as medication targeting and heat conduction in tissues. Fourier [32] was the spearheading individual who at the outset portrayed the marvels of heat transport, which is the parabolic energy equation for temperature field and has disadvantage that initial interruption is felt right away all through the entire medium. Later on, so as to control this limitation, Cattaneo [33] improved the Fourier law of heat conduction by including the thermal relaxation term which causes heat transportation as thermal waves with limited speed. Furthermore, Christov [34] utilized Oldroyd’s upper convected derivative instead of time derivative so as to achieve the material-invariant detailing. This new model is named as the Cattaneo–Christov heat flux model. After the spearheading work of Christov [34], a number of ongoing investigations in this setting are given in the references [35–45].The aforesaid writing surveys affirmed that no endeavor has been made at this point to investigate the MHD mixed convection flow of Eyring–Powell nanofluid past a vertical linearly extending sheet through a porous medium in the presence of the Cattaneo–Christov heat and mass flux model with the consolidated effects of Hall and ion slip, porosity, and Joule heating. With this perspective, the current correspondence aims to fill this gap inside the current writing. Consequently, the fundamental object of the present work is to examine the impact of Hall and ion slip on the mixed convection stream of Eyring–Powell nanofluid past a vertical linearly extending sheet through a penetrable medium in the presence of Joule heating with the Cattaneo–Christov heat and mass transition model. Numerical clarifications were accomplished employing the spectral relaxation method [3, 46–48]. The effects of embedding parameters on the velocity, temperature, and concentration fields are analyzed graphically. Moreover, numerical estimations of skin friction coefficients, local Nusselt, and Sherwood numbers for some pertinent parameters are computed and tabulated.


## 2. Mathematical Formulations

Considering a laminar, steady, viscous, and incompressible mixed convection stream of an electrically conducting Eyring-Powell nanofluid past a vertical linearly extending sheet through a porous medium with velocityUwx=ax (where a is a positive constant) towardx-axis, the Cattaneo–Chirstov heat and mass flux model are utilized to explore the heat and mass transfer qualities. The surface temperature Tw and concentration Cw are assumed consistent at the extending surface. As y extends to infinity, the encompassing values of temperature and concentration are, respectively, T∞ and C∞. Besides, a strong uniform magnetic field B0 is applied in the direction along y−axis as shown in Figure 1. Because of this strong magnetic field force, the electrically conducting liquid has the Hall and ion slip impacts which produces a crossflow in the z-direction, and consequently, the stream becomes three dimensional. The x-axis is along the course of movement of the surface, y-axis is perpendicular to the surface, and z-axis is normal to the xy-plane. The induced magnetic field can be overlooked conversely with the applied magnetic field on the assumption that the stream is steady, and the magnetic Reynolds number is very little.Figure 1
Physical representation of the flow problem.The extra stress tensorτij of the Eyring–Powell liquid model is characterized as [15, 21]:
(1)τij=μ∂ui∂xj+1βsin−11b∂ui∂xj,where μ is the dynamic viscosity, and β and b are the material fluid parameters of the Eyring–Powell fluid model such that
(2)sin−11b∂ui∂xj≅1b∂ui∂xj−161b∂ui∂xj3and1b∂ui∂xj≪1.Generalized Ohm’s law with Hall and ion slip impacts is given by [20, 23]
(3)J=σΕ+V×Β−ωeτeB0J×Β+ωeτeβiΒ02J×Β×Β,where J=Jx,Jy,Jz is the current density vector,Ε is the intensity vector of the electric field,V is the velocity vector, Β is the magnetic field, ωe is the cyclotron frequency, andτe is the electrical collision time.Applying the Boussinesq approximations alongside the above suppositions, the basic governing equations are [30, 36, 43, 44, 49]
(4)∂u∂x+∂v∂y=0(5)u∂u∂x+v∂u∂y=υ+1βbρ−12βb3ρ∂u∂y2∂2u∂y2−σB02ραe2+βe2αeu+βew+gβTT−T∞+βCC−C∞−υku(6)u∂w∂x+v∂w∂y=υ+1βbρ−12βb3ρ∂w∂y2∂2w∂y2+σB02ραe2+βe2βeu−αew−υkw(7)u∂T∂x+v∂T∂y=α∂2T∂y2−γTu∂u∂x+v∂u∂y∂T∂x+u∂v∂x+v∂v∂y∂T∂y+2uv∂2T∂x∂y+u2∂2T∂x2+v2∂2T∂y2+τDB∂T∂y∂C∂y+DTT∞∂T∂y2+σB02ρCpαe2+βe2u2+w2(8)u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2−γCu∂u∂x+v∂u∂y∂C∂x+u∂v∂x+v∂v∂y∂C∂y+2uv∂2C∂x∂y+u2∂2C∂x2+v2∂2C∂y2,with the relating boundary conditions [3, 36]:
(9)u=Uwx=ax,v=w=0,T=Tw,C=Cw,aty=0,u→0,w→0,T→T∞,C→C∞,asy→∞,where u,v,andw are the velocity components along the x−,y−,andz—directions, respectively, and υ is the kinematic viscosity, ρ is the density of the fluid, g is the acceleration due to gravity, βe is the Hall parameter, βi is the ion-slip parameter,αe=1+βeβi is a constant, BT is the coefficient of thermal expansion, BC is the solutal coefficient of expansion, κ is the permeability of the porous medium, T and C are the fluid temperature and concentration, respectively, γT and γC are the heat and mass flux relaxation times correspondingly, Cp is the specific heat capacity, τ is the ratio of the effective heat capacity of nanoparticle to the heat capacity of the fluid, i.e.,τ=ρCp/ρCf,DB is the Brownian diffusion coefficient, and DT is the thermophoresis diffusion coefficient.The overseeing equations can be changed to coupled nonlinear ordinary differential equations utilizing the subsequent similarity transformations [3, 36]:
(10)η=aυy,ψ=aυxfη,w=axgη,θη=T−T∞Tw−T∞,ϕη=C−C∞Cw−C∞,where η is the similarity variable, and ψ is the stream function defined as
(11)u=∂ψ∂yandv=−∂ψ∂x.Equation (4) is identically fulfilled, and equations (5)–(9) become
(12)1+εf′′′−εδf′′2f′′′−f′2+ff′′−Mαe2+βe2αef′+βeg+λθ+Nrϕ−Κf′=0,(13)1+εg′′−εδg′2g′′−f′g+fg′+Mαe2+βe2βef′−αeg−Κg=0,(14)1Prθ′′+fθ′−γ1ff′θ′+f2θ′′+Nbθ′ϕ′+Ntθ′2+M.Ecαe2+βe2f′2+g2=0,(15)ϕ′′+Pr.Lefϕ′−γ2ff′ϕ′+f2ϕ′′+NtNbθ′′=0.The transformed boundary conditions are(16)f0=0,f′0=1,g0=0,θ0=1,ϕ0=1,f′⟶0,g⟶0,θ⟶0,ϕ⟶0,asη⟶∞.The dimensionless parameters showing up in equations (12)–(15) can be characterized asM=σB02/ρa, αe=1+βeβi,λ=Grx/Rex2, Grx=gβTTw−T∞x3/υ2, Rex=ax2/υ, Nr=βCCw−C∞/βTTw−T∞, K=υ/aκ, ε=1/βbμ,δ=a3x2/2b2υPr=υ/α, Nb=τDBCw−C∞/υ, Nt=τDTTw−T∞/υT∞, τ=ρCp/ρCf, Ec=ax2/CpTw−T∞, Le=α/DB,γ1=γTa and γ2=γCa.The skin friction coefficientsCfxandCgz, the local Nusselt number Nux, and the Sherwood number Shx are the physical quantities of interest which can be characterized as follows:
(17)Cfx=τwxρUw2,Cfz=τwzρUw2,Nux=xqwkTW−T∞,Shx=xjmDBCW−C∞,where τwx and τwz are the wall shear stresses in the directions of x and z, respectively, and qw is the surface heat flux, and jm is the surface mass flux, which are given by
(18)τwx=μ+1βb∂u∂y−16βb3∂u∂y3y=0,τwz=μ+1βb∂w∂y−16βb3∂w∂y3y=0,qw=−k∂T∂yy=0,jm=−DB∂C∂yy=0.From equations (10), (11), (17), and (18), we obtain
(19)CfxRex=1+εf′′0−εδ3f′′03,CgzRex=1+εg′0−εδ3g′03,NuxRex−1=−θ′0,ShxRex−1=−ϕ′0.


## 3. Method of Solution

Equations (12)–(15) subject to the boundary conditions (16) are explained utilizing the spectral relaxation method. The spectral relaxation algorithm uses the thought of the Gauss-Seidel method to decouple the arrangement of overseeing equations (12)–(15). The strategy is created by assessing the linear terms at the current iteration level r+1 and nonlinear terms at the former iteration level r. The Chebyshev pseudo-spectral technique is used to unravel the decoupled equations. For the subtleties of the technique, interested readers can refer [3, 46–48].Hence, to utilize the SRM, we start by decreasing the request for the momentum (equation (12)) from third to second order presenting the change f′=h so that f′′=h′andf′′′=h′′.Along these lines, equations (12)–(15) become
(20)f′=h(21)1+εh′′−εδh′2h′′−h2+fh′−Mαe2+βe2αeh+βeg+λθ+Nrϕ−Κh=0,(22)1+εg′′−εδg′2g′′−hg+fg′+Mαe2+βe2βeh−αeg−Κg=0,(23)1Prθ′′+fθ′−γ1fhθ′+f2θ′′+Nbθ′ϕ′+Ntθ′2+M.Ecαe2+βe2h2+g2=0,(24)ϕ′′+Pr.Lefϕ′−γ2fhϕ′+f2ϕ′′+NtNbθ′′=0.The transformed boundary conditions are(25)f0=0,h0=1,g0=0,θ0=1,ϕ0=1,h⟶0,g⟶0,θ⟶0,ϕ⟶0,asη⟶∞.Executing the SRM to equations (20)–(25), we get the following iteration scheme:
(26)1+εhr+1′′−εδhr′2hr′′−hr2+frhr+1′−Mαe2+βe2αehr+1+βegr+λθr+Nrϕr−Κhr+1=0,(27)fr+1′=hr+1,(28)1+εgr+1′′−εδgr′2gr′′−hr+1gr+1+fr+1gr+1′+Mαe2+βe2βehr+1−αegr+1−Κgr+1=0,(29)1Prθr+1′′+fr+1θr+1′−γ1fr+1hr+1θr+1′+fr2θr+1′′+Nbθr+1′ϕr+1′+Ntθr′2+M.Ecαe2+βe2hr2+gr2=0,(30)ϕr+1′′+Pr.Lefr+1ϕr+1′−γ2fr+1hr+1ϕr+1′+fr2ϕr+1′′+NtNbθr+1′′=0.The transformed boundary conditions are(31)fr+10=0,hr+10=1,gr+10=0,θr+10=1,ϕr+10=1,hr+1⟶0,gr+1⟶0,θr+1⟶0,ϕr+1⟶0,asη⟶∞.Actualizing the Chebyshev spectral collocation method to equations (26)–(30), we get the ensuing matrix equations:
(32)A1hr+1=B1,hr+1ξN=1,hr+1ξ0=0,(33)A2fr+1=B2,fr+1ξN=0,(34)A3gr+1=B3,gr+1ξN=0,gr+1ξ0=0,(35)A4θr+1=B4,θr+1ξN=1,θr+1ξ0=0,(36)A5ϕr+1=B5,ϕr+1ξN=1,ϕr+1ξ0=0,where
(37)A1=1+εD2+diagfrD−Mαeαe2+βe2+KI,B1=εδhr′2hr′′+hr2+Mβeαe2+βe2gr−λθr+Nrϕr,(38)A2=D,B2=hr+1,A3=1+εD2+diagfr+1D−diaghr+1+Mαeαe2+βe2+KIB3=εδgr′2gr′′−Mβeαe2+βe2hr+1,A4=diag1Pr−γ1fr2D2+diagfr+1−γ1fr+1hr+1+Nbϕr′D,B4=−M.Ecαe2+βe2hr2+gr2−Ntθr′2,A5=diag1−PrLeγ2fr2D2+PrLediagfr+1−γ2fr+1hr+1D,B5=−NtNbθr+1′′.Here,D=2D/η∞ where D is the Chebyshev differentiation matrix (Trefethen [50]), η∞ is a limited length which is adequately big so that we can simply incorporate the condition at perpetuity in this point, I and diag [.] are the identity and diagonal matrices of order N+1×N+1, N is the number of grid points, andf,h,g,θand ϕ are, respectively, the values of f, h, g, θ, andϕ, when assessed at the grid points.The decoupled equations (32)–(36) can be solved separately by picking appropriate introductory suppositions:
(39)f0η=1−e−η,g0η=ηe−η,θ0η=e−η=ϕ0η.


## 4. Results and Discussions

The nonlinear differential equations (12)–(15) with the limit conditions (16) have been explained numerically via the spectral relaxation method, and consequences of the numerical calculations for the velocity, temperature, and concentration profiles as well as skin friction coefficients, local Nusselt number, and Sherwood number have been obtained for different input parameters and presented through graphs and tables. Also, we contrasted our outcomes and those of the current writing as appeared in Table 1, and the outcomes give a generally excellent understanding which gives a certainty of our current numerical outcomes. In this study, for our numerical computations, we fixed the numerical values of the physical parameters as =0.1, Le=1.0, M=10.0, βe=βi=5.0, Pr=0.72, K=0.5, Nr=Nb=Nt=0.1, Ec=0.02, δ=0.2, and ε=γ1=γ2=0.3 unless otherwise stated.Table 1
Comparison of results for the Nusselt number−θ′0 and Sherwood number −ϕ′0 for various values of Nb and Nt whenLe=1,Pr=10, and M=λ=βe=βi=Nr=ε=δ=γ1=γ2=K=Ec=0 are fixed.
NbNtNoghrehabadi et al. [51]Goyal and Bhargava [52]Ibrahim and Gadisa [53]Present results−θ′0−ϕ′0−θ′0−ϕ′0−θ′0−ϕ′0−θ′0−ϕ′00.10.10.95237682.12939380.952442.129490.952392.129380.95237682.12939380.20.69317432.27402150.693182.274010.693182.274010.69317432.27402160.30.52007902.52863820.520252.528550.520192.528550.52007902.52863870.40.40258082.79517010.402602.795200.402582.795190.40258232.79516280.50.32105433.03514250.321053.035110.321093.035110.32104913.03516860.20.10.50558142.38187060.505612.381860.505572.381900.50558152.38187060.30.25215602.41001880.252182.410090.252172.410020.25215622.41001880.40.11940592.39965020.119402.399700.119402.399660.11940582.39965020.50.05425342.3835712————0.05425282.3835713Figures2–5 reveal the impact of the Eyring-Powell fluid parameter ε on the primary velocity, secondary velocity, temperature, and concentration distributions, respectively. By expanding the liquid parameter ε, the primary velocity profile rises and on the contrary, the temperature and concentration profiles diminish. Because the fluid parameter ε=1/μβb has an inverse relation with viscosity of fluid μ; thus, fluid turns out to be less viscous for bigger values of ε which upsurges fluid velocity. Moreover, the secondary velocity profile reduces near the surface and enhances far off from the surface. This is possibly due to the fact that the viscosity can be varied as the Eyring-Powell fluid parameter varies resulting in an increasing viscous force as a result decreasing its velocity up to some maximum value, and then the velocity starts to increase due to decrease in viscosity. The optimum value is obtained around the point ≈1.Figure 2
Effect ofε on the primary velocity f′η.Figure 3
Effect ofε on the secondary velocity gη.Figure 4
Effect ofε on the temperature θη.Figure 5
Effect ofε on the concentration ϕη.Figures6–9 divulge the impact of magnetic field parameter M on primary velocity, secondary velocity, temperature, and concentration profiles, respectively. The transverse magnetite field which is applied vertical to the direction of flow provides a resistive force known as the Lorentz force. Physically, the Lorentz force opposes the flow of nanofluid; hence, the velocity diminishes. Moreover, this force tends to hinder the movement of the liquid in the boundary layer. This assists with declining the primary velocity profile and upgrades the secondary velocity, thermal, and concentration profiles as M increments as appeared in Figures 6–9. It is additionally seen that an augmentation in M prompts a lower in primary velocity boundary layer thickness and superior in thermal and concentration boundary layer thicknesses, though the secondary velocity boundary layer thickness diminishes near the surface and upsurges far off from the surface.Figure 6
Effect ofM on the primary velocity f′η.Figure 7
Effect ofM on the secondary velocity gη.Figure 8
Effect ofM on the temperature θη.Figure 9
Effect ofM on the concentration ϕη.Figures10–13 illustrate the impact of the Hall parameter βe on primary velocity, secondary velocity, temperature, and concentration profiles, respectively. It is seen that the primary velocity profile increments with an expansion in βe, though the contrary condition is watched for secondary velocity, thermal, and concentration profiles. Besides, as the Hall parameter βe enlarges, the primary velocity boundary layer thickness upgrades, and thermal and concentration boundary layer thicknesses lessen, while the secondary velocity boundary layer thickness improves near the surface and diminishes distant from the surface. Physically, the consideration of Hall parameter decreases the efficient conductivity and subsequently drops the magnetic resistive force. Thus, the Joule heating impact is condensed, and quantity of heat owing to ohmic dissipation can be abridged by utilizing partially ionized fluid presented to magnetic field. The comparable conduct is seen with expanding ion slip parameter βi as revealed in Figures 14–17. Physically, the expansion in βi improves the efficient conductivity and therefore diminishes the damping force on the dimensionless velocity, and hence the dimensionless velocity expands.Figure 10
Effect ofβe on the primary velocity f′η.Figure 11
Effect ofβe on the secondary velocity gη.Figure 12
Effect ofβe on the temperature θη.Figure 13
Effect ofβe on the concentration ϕη.Figure 14
Effect ofβi on the primary velocity f′η.Figure 15
Effect ofβi on the secondary velocity gη.Figure 16
Effect ofβi on the temperature θη.Figure 17
Effect ofβi on the concentration ϕη.Figures18–21 portray the impact of blended convection parameter λ on primary velocity, secondary velocity, temperature, and concentration profiles, respectively. It is deduced that as λ upgrades, the primary velocity and secondary velocity profiles ascend. In contrast, the temperature and concentration profiles reduce. These outcomes physically hold since upgrade in blended convection parameter causes upsurge in buoyancy forces which speed up fluid motion.Figure 18
Effect ofλ on the primary velocity f′η.Figure 19
Effect ofλ on the secondary velocity gη.Figure 20
Effect ofλ on the temperature θη.Figure 21
Effect ofλ on the concentration ϕη.Figures22–25 outline the impact of penetrability parameter K on the primary velocity, secondary velocity, thermal, and concentration distributions, respectively. It is obvious from the figures that as K augments, both primary and secondary velocity profiles decline while the opposite condition is watched for thermal and concentration profiles. Physically, the porousness builds obstruction of the permeable medium which will in general, diminish the fluid velocity. It is normal that an expansion in the porousness of the permeable medium prompts the ascent in the flow of fluid through it. When the outlets of the permeable medium become big, the obstruction of the medium might be ignored.Figure 22
Effect ofK on the primary velocity f′η.Figure 23
Effect ofK on the secondary velocity gη.Figure 24
Effect ofK on the temperature θη.Figure 25
Effect ofK on the concentration ϕη.Figure26 proves the impact of the thermal relaxation parameter γ1 on temperature distribution. It uncovers that the fluid temperature and energy boundary layer decrease with expanding value of γ1. Physically, thermal relaxation time is the time required by the fluid particles to move heat energy to its neighboring particles. Subsequently, as γ1ascends, the material particles need additional time to transfer heat to its adjoining particles, and this prompts less transport of heat from the surface to the fluid. Similarly, the concentration relaxation parameter γ2emerges from the relaxation time of mass flux. Less relaxation time of mass flux prompts superior concentration. and the fluid with increasingly mass flux relaxation time brings about smaller concentration. Therefore, increasingly mass flux relaxation time happens owing to the bigger value of γ2 which diminishes the concentration profile as revealed in Figure 27.Figure 26
Effect ofγ1 on the temperature θη.Figure 27
Effect ofγ2 on the concentration ϕη.Figure28 depicts the impact of Eckert number Ec on the temperature profile. It is seen that the temperature profile upsurges with an expansion in Ec. Physically, by expanding the Eckert number, the heat energy is hoarded in the fluid as a result of the frictional or drag forces. Accordingly, the fluid temperature field enhances.Figure 28
Effect ofEc on the temperature θη.Figure29 demonstrates the effect of Lewis number Le on the concentration profile ϕη. By expanding Le, the concentration and the boundary layer thicknesses decline. This is because of the way that the mass transfer rate enhances by expanding Le which improves the volume fraction of nanoparticles.Figure 29
Effect ofLe on the concentration ϕη.Table2 specifies the impact of physical parameters on the skin fraction coefficients, local Nusselt, and Sherwood numbers. The primary skin friction coefficient −ReCfx augments as magnetic parameter M, Eyring-Powell fluid parameter ε, thermal relaxation parameter γ1, concentration relaxation parameter γ2, and permeability parameter K rise and decreases for increasing values of mixed convection parameter λ, Hall parameter βe, ion slip parameter βi, and Eckert number Ec. Also, it is observed that the secondary skin friction coefficient ReCgz amplifies as magnetic parameter M, mixed convection parameter λ, and Eyring-Powell fluid parameter ε increase and declines as Hall parameter βe, ion slip parameter βi, thermal relaxation parameter γ1, concentration relaxation parameter γ2, and permeability parameter K increase, while it remains constant as the Eckert number Ec enhances. Similarly, both the local Nusselt number –θ′0 and Sherwood number –ϕ′0 enhance for larger values of mixed convection parameter λ, Hall parameter βe, ion slip parameter βi, and Eyring-Powell fluid parameter ε and decrease for larger values of magnetic parameter M and permeability parameter K. Further, opposite behaviors of the local Nusselt number –θ′0 and the local Sherwood number –ϕ′0 are observed for increasing values of thermal relaxation parameter γ1, concentration relaxation parameter γ2, and Eckert number Ec.Table 2
Numerical values of the skin friction coefficientsReCfx=1+εf′′0−εδ/3f′′03 and ReCgz=1+εg′0−εδ/3g′03, local Nusselt number –θ′0, and local Sherwood number –ϕ′0 for various values of M, λ, βe, βi, ε, γ1, γ2, K, and Ec when Nr=1.0, Pr=0.72, Nb=0.2, Nt=0.1, Le=1, and δ=10 are fixed.
Mλβeβiεγ1γ2KEc−ReCfxReCgz−θ′0−ϕ′05.00.55.05.00.30.50.52.00.021.476040.014320.437510.3305010.01.534470.027570.430110.3238915.01.590390.039870.422960.3175510.00.21.774360.024700.378740.274710.41.615850.026710.415400.309920.61.452920.028360.443280.336310.55.01.534470.027570.430110.3238910.01.477130.014850.437400.3303715.01.456990.010160.439930.332655.05.01.534470.027570.430110.3238910.01.478830.007630.437190.3301815.01.458390.003500.439760.332495.00.01.382790.024960.421300.316140.51.635960.029140.435810.329011.01.885550.032700.449400.341560.30.01.531320.027690.407680.335890.51.534470.027570.430110.323891.01.538130.027430.456280.310080.50.01.530500.027710.436010.303510.51.534470.027570.430110.323891.01.539120.027400.423550.347720.51.01.187780.034740.473140.363592.01.534470.027570.430110.323893.01.798300.023030.394310.292622.00.01.534470.027570.430190.323850.11.534460.027570.429800.324040.21.534440.027570.429410.32422


## 5. Conclusions

This paper looks at the mixed convection stream of Eyring-Powell nanofluid with the Cattaneo-Christov heat and mass flux model over linearly stretching sheet through a permeable medium. The impacts of Hall and ion slip, permeability, and Joule heating are also considered. The governing nonlinear partial differential equations are changed into joined nonlinear ordinary differential equations utilizing similitude transformations. These nonlinear ordinary differential equations subject to boundary conditions are handled numerically by means of SRM. Numerical results are presented for different parameters of interest, and results are discussed with the assistance of graphs and tables. The study has a momentous application in the material processing technology. From the introduced investigation, the accompanying ends were drawn:(i)
Primary velocity profile upsurges by increasing the values of Eyring-Powell fluid parameterε, Hall parameter βe, ion slip parameter βi, and mixed convection parameter λ, while opposite condition is noticed for mounting the values of magnetic field parameter M and porosity parameter K(ii)
Secondary velocity profile increases by increasing the values of Eyring-Powell fluid parameterε and magnetic field parameter M, while reverse behavior is observed for increasing values of Hall parameter βe, ion slip parameter βi, mixed convection parameter λ, and porosity parameter K(iii)
Both temperature and concentration profiles enhance with increasing values of magnetic field parameterM and porosity parameter K, whereas reverse trend is noticed for increasing values of Eyring-Powell fluid parameter ε, Hall parameter βe, ion slip parameter βi, mixed convection parameter λ, thermal relaxation parameter γ1, and concentration relaxation parameter γ2. Moreover, temperature profile increases as Eckert numberEc increases, while concentration profile decreases as Lewis numberLe increases(iv)
Both skin friction coefficients increase with increasing values of magnetic field parameterM and Eyring-Powell fluid parameter ε and decrease with increasing values of Hall parameter βe and ion slip parameter βi. On the contrary, opposite behavior is observed for the primary and secondary skin friction coefficients as mixed convection parameter λ, porosity parameter K, thermal relaxation parameter γ1, and concentration relaxation parameter γ2 upsurges. An increase in Eckert number Ec shows insignificant change in the primary skin friction coefficient, while the secondary skin friction coefficient remains constant(v)
Both local Nusselt and Sherwood numbers increase by enhancing the values of mixed convection parameterλ, Hall parameter βe, ion slip parameter βi, and Eyring-Powell fluid parameter ε, while opposite behavior is noticed for increasing values of magnetic field parameter M and permeability parameter K. On the contrary, opposite behavior is observed for local Nusselt and Sherwood numbers as Eckert numberEc, thermal relaxation parameter γ1, and concentration relaxation parameter γ2 increase



---

*Source: 4354860-2020-09-09.xml*