# Magnetohydrodynamic (MHD) Boundary Layer Flow Past a Wedge with Heat Transfer and Viscous Effects of Nanofluid Embedded in Porous Media

Authors: Wubshet Ibrahim; Ayele Tulu
Journal: Mathematical Problems in Engineering
 (2019)

---

## Abstract

The problem of two-dimensional steady laminar MHD boundary layer flow past a wedge with heat and mass transfer of nanofluid embedded in porous media with viscous dissipation, Brownian motion, and thermophoresis effect is considered. Using suitable similarity transformations, the governing partial differential equations have been transformed to nonlinear higher-order ordinary differential equations. The transmuted model is shown to be controlled by a number of thermophysical parameters, viz. the pressure gradient, magnetic, permeability, Prandtl number, Lewis number, Brownian motion, thermophoresis, and Eckert number. The problem is then solved numerically using spectral quasilinearization method (SQLM). The accuracy of the method is checked against the previously published results and an excellent agreement has been obtained. The velocity boundary layer thickness reduces with an increase in pressure gradient, permeability, and magnetic parameters, whereas thermal boundary layer thickness increases with an increase in Eckert number, Brownian motion, and thermophoresis parameters. Greater values of Prandtl number, Lewis number, Brownian motion, and magnetic parameter reduce the nanoparticles concentration boundary layer.


---

## Body

## 1. Introduction

Fluid flows with convective heat and mass transfer over a wedge shaped bodies is ensured in many thermal engineering applications like crude oil extraction, geothermal systems, thermal insulation, heat exchangers, and the storage of nuclear waste, Nagendramma et al. [1]. A model of steady laminar fluid flow over a wedge has developed for the first time by Falkner and Skan [2] to illustrate the application of Prandtl’s boundary layer theory. Late, Hartree [3] investigated the same problem with similarity transformation and gave numerical results for wall shear stress for different values of the wedge angle. Eckert [4] also solved Falkner-Skan flow along an isothermal wedge and presented the first wall heat transfer values. Afterward, the variety of applications and understanding of the physical features of laminar boundary layer flow past a wedge have motivated many researchers (to mention a few, Martin and Boyd [5], Sattar [6], Kandasamy et al [7], and Turkyilmazoglu [8]).Magnetohydrodynamic (MHD) is the study of fluid flow in electrically conducting fluids with magnetic properties that affect fluid flow characteristics. When a magnetic field is incident in an electrically conducting fluid, a current is induced. This effect polarizes the fluid and as a result the magnetic field is changed (Makanda et al. [9]). Due to extensive practical applications of MHD in technological processes such as plasma studies, petroleum industries, MHD power generator designs, design for cooling of nuclear reactors, and construction of heat exchangers and on the performance of many other systems, there are many studies that considered MHD fluid flow past a wedge. These include the work of Abbasbandy et al. [10] who examined the effects of MHD in the Falken-Skan flow of Maxwell fluid, and El-Dabe et al. [11] considered the MHD boundary layer flow of non-Newtonian Casson fluid on a moving wedge with heat and mass transfer. Khan et al. [12] also analyzed MHD laminar boundary layer flow past a wedge with the influence of thermal radiation, heat generation, and chemical reaction.In the past few years, MHD boundary layer flow with heat and mass transfer of nanofluids has become a major topic of modern-day interest. Nanofluids have a significant role in enhancing the heat transfer properties of fluids. The most important properties of nanofluids are enhanced effective fluid thermal conductivity and heat transfer coefficient. Some recent studies on MHD boundary layer flow of nanofluid include the work of the Srinivasacharya et al. [13] who analyzed the steady laminar MHD flow in a nanofluid over a wedge in the presence of a variable magnetic field. Rasheed et al. [14] also investigated MHD boundary layer flow of nanofluid over a continuously moving stretching surface. The effects of thermal radiation on mixed convection flow of nanofluid over a stretching sheet in the presence of a magnetic field are studied by Nageeb et al. [15].Viscous dissipation effect changes the temperature distribution by playing a role like an energy source, which leads to affected heat transfer rate and hence needs to be considered in heat transfer problems. The analysis of MHD boundary layer flows in porous media with and without the effect of viscous dissipation has been a subject of several recent researchers. Accordingly, Arthur et al. [16] analyzed hydromagnetic stagnation point flow over a porous stretching surface in the presence of thermal radiation and viscous dissipation effects. Also, Ramesh et al. [17] studied the MHD boundary layer flow past a constant wedge within porous media. Heat and mass transfer of MHD flow of nanofluids in the presence of viscous dissipation effects are numerically analyzed by Haile and Shankar [18]. Majety et al. [19] studied the effect of viscous dissipation on MHD boundary layer flow past a wedge through porous medium. They concluded that viscous dissipation produces heat due to drag between the fluid particles, which cause an increase in fluid temperature.Most of the standard methods of solving the boundary layer problems are the numerical approach based on the shooting algorithm with the Runge-Kutta scheme, finite difference method, spectral homotopy analysis method, and Newton-Raphson based methods such as the quasilinearization method and the successive linearization method. Recently, spectral based numerical techniques such as the Spectral Quasilinearization Method and Spectral Relaxation Method have been developed (see Motsa et al. [23], Motsa [24], and Magagula et al. [25]). As indicated by Motsa et al. [26] and Zhou [27] Chebyshev spectral collocation methods are easy to implement and adaptable to various problems and provide more accurate approximations with a relatively small number of unknowns. Gottlieb and Hesthaven [28] also added that the wide use of spectral methods has motivated the researcher by their accuracy and efficiency in solving incompressible Navier-Stokes equations. Furthermore, Motsa et al. [29] stated that the interest in using Chebyshev spectral collocation methods in solving nonlinear PDEs stems from the fact that these methods require less grid points to achieve accurate results and efficient compared to traditional methods like finite difference and finite element methods.The present study deals with MHD boundary layer flow past a wedge with heat and mass transfer of nanofluid embedded in a porous media with viscous dissipation, Brownian motion, and thermophoresis parameter effects. Using appropriate similarity transformation, the governing nonlinear partial differential equations are transformed to nonlinear higher-order ordinary differential equations. Then, these equations are numerically solved using Spectral Quasilinearization Method. The effects of different parameters on velocity, temperature, and concentration fields are investigated and analyzed with the help of their graphical and tabular representations.


## 2. Mathematical Formulation

Consider steady two-dimensional, laminar boundary layer flow past a wedge with heat transfer of incompressible electrically conducting nanofluid embedded in a porous media with viscous dissipation effects. The coordinate system is chosen withx coordinate pointing parallel to the plate in the direction of the flow and y coordinate pointing towards the free stream, as shown in the Figure 1. The wall of the wedge is maintained with uniform and constant temperature Tw and nanoparticle concentration Cw. Tw and Cw are, respectively, greater than the ambient temperature T∞ and ambient nanoparticle concentration C∞. The subscripts w and ∞ denote conditions at wall and ambient, respectively. The fluid is assumed to have constant physical properties.Figure 1
Physical model of flow configuration and coordinate system.It is also assumed that a constant magnetic fieldB0 is applied in the positive y-direction, normal to the walls of the wedge. The induced magnetic field caused by the motion of electrically conducting fluid is neglected, as it is very small compared to applied magnetic field (Ullah et al. [22]). The Naiver-Stokes equation of motion including electromagnetic body force or Lorentz force within conductive media is given by F=Jc×B, where Jc is the conduction current defined as Jc=σ(E+u×B), B is the magnetic field, E is the electric field, u is the velocity vector of the fluid, and σ is the electrical conductivity of the fluid (Nicholas [30]). Since the induced magnetic field caused by the motion of electrically conducting fluid E is negligible as compared to the applied magnetic field B, then Jc=σ(u×B)=σ[(u,v,0)×0,By,0)]=σuByk^. Thus, the Lorentz force F=Jc×B=(0,0,σuBy)×(0,By,0)=-σBy2ui^=-σB02u.We use the general model of the conservation equation for a general scalar variableφ which can be expressed as ( see Versteeg and Malalasekera [31])(1)∂ρφ∂t+∇.ρuφ=∇.Γ∇φ+Sφwhere the velocity vector given by u=u,v,w, ρ is density of the fluid, Γ is the diffusion coefficient, and Sφ is the source term. Equation (1) is the so-called transport equation for property φ, and it clearly highlights the various transport processes.By settingφ equal to 1, u, temperature of fluid T, and nanoparticle concentration C in (1) and selecting appropriate values for diffusion coefficient Γ and source terms Sφ, we obtain special forms of PDEs (Navier-Stokes equations) for the continuity, momentum, energy, and nanoparticle concentration. With the above assumptions and the boundary layer approximations, the continuity, momentum, energy, and nanoparticle concentration equations governing steady two-dimensional MHD boundary layer flow past a wedge embedded in a porous media with viscous dissipation effects and constant fluid properties are given as (see Srinivasacharya et al. [13], Haile and Shankar [18], and Alam et al. [32]).(2)∂u∂x+∂v∂y=0(3)u∂u∂x+v∂u∂y=-1ρf∂p∂x+νf∂2u∂y2-σB02ρf+νfKu(4)u∂T∂x+v∂T∂y=αf∂2T∂y2+νfCp∂u∂y2+τDB∂T∂y∂C∂y+DTT∞∂T∂y2(5)u∂C∂x+v∂C∂y=DB∂2C∂y2+DTT∞∂2T∂y2The appropriate boundary conditions are given as(6)u=0;v=0;T=Tw;C=Cw,aty=0(7)u=Ux=U∞xm;T→T∞;C→C∞,asy→∞(8)u=U∞;T=T∞,atx=0where u and v are, respectively, the x and y velocity components. ρf, νf, αf, and Cp, are respectively, density, kinematic viscosity, thermal diffusivity, and specific heat capacity of the base fluid. K is the permeability of porous medium, DB is the Brownian diffusion coefficient, DT is the thermophoresis diffusion coefficient, and τ is the ratio of the effective heat capacity of the nanoparticle material and the heat capacity of base fluid.The y-momentum equation implies that the pressurep in the boundary layer must be equal to that of the free stream for any given x coordinate. Because the velocity profile is uniform in the free stream, there is no vorticity involved; therefore, simple Bernoulli’s equation can be applied in this high Reynolds number (Falkner and Skan [2]). It is assumed that U(x)=U∞xm is the fluid velocity at the wedge outside the boundary layer, where U∞ is the free stream velocity. For a uniform stream, the momentum equation (3) becomes (see Falkner and Skan [2], and Nageeb et al. [15])(9)-1ρf∂p∂x=UdUdx+σB02ρf+νfKUSubstituting (9) into (3), the momentum equation is written as(10)u∂u∂x+v∂u∂y=UdUdx+νf∂2u∂y2+σB02ρf+νfKU-uHere, x is measured from the tip of the wedge, m is the Falkner-Skan power-law parameter, and β=2m/1+m is the Hartree pressure gradient parameter corresponding to β=Ω/π for the total angle Ω of the wedge (see Figure 1). Physically, m<0 corresponding to an adverse pressure gradient (often resulting in boundary layer separation) while m>0 represents favorable pressure gradient (Nagendramma et al. [1]). In the Blasius solution m=0 corresponding to an angle of attack of zero radians, where m=1 corresponding to stagnation point flow.In order to transform the governing equations (2)-(7) to a set of ordinary differential equations, introduce the stream function ψ(x,y) such that(11)u=∂ψ∂yv=-∂ψ∂xand we use the following transformation variables (see Ullah et al. [22], and Alam et al. [32]):(12)η=y1+m2U∞νfxm-1/2;ψx,η=21+mνfU∞xm+1/2fηf′η=uU;θη=T-T∞Tw-T∞;ϕη=C-C∞Cw-C∞where η is a dimensionless similarity variable, f(η) nondimensional stream function, f′(η) is non-dimensional velocity, θ(η) is nondimensional temperature, and ϕ(η) is nondimensional nanoparticle concentration.Upon substituting similarity variables into (3)-(7), the continuity equation (2) is identically satisfied and the following system of ordinary differential equations is obtained:(13)f′′′+ff′′+β1-f′2+11+mM+κ1-f′=0(14)θ′′+Prfθ′+Ecf′′2+Nbθ′ϕ′+Ntθ′2=0(15)ϕ′′+LePrfϕ′+NtNbθ′′=0The transformed boundary conditions are(16)f=0;f′=0;θ=1;ϕ=1,atη=0;(17)f′→1;θ→0;ϕ→0,asη→∞with(18)M=2σB02x1-mρfU∞;κ=2νfx1-mKU∞;Ec=U2CpTw-T∞;Pr=νfαfLe=αfDB;Nb=τDBCw-C∞νf;Nt=τDTTw-T∞νfT∞;Rex=Uwxνfwhere M is magnetic parameter, Rex is local Reynolds number, Pr is Prandtl number, Ec is Eckert number, κ is the permeability parameter, Le is Lewis number, Nb is the Brownian motion parameter, Nt is the thermophoretic parameter, and prime (′) denotes derivative with respect to η.The physical quantities of engineering interest in the present study are the skin friction coefficientCf, local Nusselt number Nux, and local Sherwood number Shx, respectively, and defined as(19)Cf=2τwρU2x;Nux=xqwkTw-T∞;Shx=xMwDBCw-C∞;where τw, qw, and Mw are the surface shear stress, the surface heat flux, and surface mass flux; and, respectively, they are given as(20)τw=μf∂u∂yy=0;qw=-kf∂T∂yy=0;Mw=-DB∂C∂yy=0The nondimensional skin friction coefficient, local Nusselt number, and local Sherwood number are, respectively, given as(21)CfRex=2m+12f′′0;NuxRex=-m+12θ′0;ShxRex=-m+12ϕ′0


## 3. Numerical Method

The system of nonlinear ODE (13)-(15) subjected to the boundary conditions (17) has been solved numerically using spectral quasilinearization method (SQLM). The main idea behind this method is identifying univariate and multivariate nonlinear terms of function and its derivative in each of the equations of the systems (13)-(15), linearizing the terms and applying Chebychev pseudospectral collocation method (see Motsa [33]).Applying spectral quasilinearization method, (13)-(15) give the following iterative sequence of linear differential equations:(22)fr+1′′′+a1,rfr+1′′-a2,rfr+1′+a3,rfr+1=a4,r(23)θr+1′′+b1,rθr+1′+b2,rfr+1′′+b3,rfr+1+b4,rϕr+1′=b5,r(24)ϕr+1′′+c1,rϕr+1′+c2,rfr+1′′+c3,rfr+1+c4,rθr+1′=c5,rwhere the terms containing r+1 subscripts denote current approximations and terms containing r subscripts denote previous approximations. The corresponding boundary conditions are(25)fr+10=0;fr+1′0=0,fr+1′∞→1(26)θr+10=ϕr+10=1,θr+1∞=ϕr+1∞→0where(27)a1,r=fr;a2,r=2βfr′+11+mM+κ;a3,r=fr′′;a4,r=frfr′′-β1+fr′2-11+mM+κb1,r=Pr⁡fr+Nbϕr′+2Ntθr′;b2,r=2PrEcfr′′;b3,r=Prθr′;b4,r=PrNbθr′b5,r=PrEcfr′′2+PrNbθr′ϕr′+PrNtθr′2+Prfrθr′;c1,r=PrLefr-Ntθr′c3,r=Pr⁡Leϕr′-NtNbθr′;c4,r=-Pr⁡NtNbfr+Nbϕr′+2Ntθr′c2,r=-2PrEcNtNbfr′′;c5,r=PrLefrϕr′-NtNbEcfr′′2+Nbθr′ϕr′+Ntθr′2+frθr′The physical domain on which the system of governing equations (13)-(15) defined in [0,∞) is moved to [-1,1] using the transformation x=2η/L∞-1, where L∞ is a scaling parameter assumed to be large and the interval [0,∞) is replaced by [0,L∞]. Spectral collocation method is applied to the system of (22)-(24); and the differentiation matrix D=2D/L∞ is used to approximate derivatives of unknown variables, where D is (N+1)×(N+1) Chebyshev differentiation matrix (see Trefethen [34]). The system of (22)-(24) is solved as a coupled matrix:(28)Λ11Λ12Λ13Λ21Λ22Λ23Λ31Λ32Λ33Fr+1Θr+1Φr+1=R1R2R3with transformed boundary condition(29)Fr+1xNx=0,Fr+1xNx-1=0,Fr+1x0=1(30)Θr+1xNx=1,Θr+1x0=0;Φr+1xNx=1,Φr+1x0=0where (31)R1=a4,r;R2=b5,r;R3=c5,rΛ11=D3+diaga1,rD2-diaga2,rD+diaga3,r;Λ12=O;Λ13=OΛ21=diagb2,rD2+diagb3,r;Λ22=D2+diagb1,rD;Λ23=diagb4,rDΛ31=diagc2,rD2+diagc3,r;Λ32=diagc4,rD;Λ33=D2+diagc1,rDFr+1=fr+1,0,fr+1,1,---,fr+1,NT, Θr+1=θr+1,0,θr+1,1,---,θr+1,NT, and Φr+1=ϕr+1,0,ϕr+1,1,---,ϕr+1,N,T are vectors of sizes (Nx+1)×1, diag(⋯) represents a diagonal matrix of vectors, and O is a zero vectors of size (Nx+1)×(Nx+1).The suitable initial approximations that satisfy the governing boundary conditions of the boundary layer equations (13)-(17) are(32)f0η=η-1+e-η,f0′η=1-e-ηθ0η=ϕ0η=e-η


## 4. Results and Discussion

Numerical solutions are obtained using SQLM for the velocity, temperature, and concentration profiles across the boundary layer for different values of the governing parameters. The number of collocation points in the spacex variable used to generate the results is Nx=40 in all cases. To ensure the numerical accuracy of the numerical method used, the skin friction coefficient -f′′(0) and local Nusselt number -θ′(0) have been calculated for different values of Falkner-Skan power-law parameter m. From Table 1, it is observed that the data produced by the SQLM code and those reported by Ashwini and Eswara [20], Watanaba [21], and Ullah et al. [22] are excellent agreement. Thus, we are very much confident that the present results are accurate.Table 1
Comparison of the SQLM results of skin friction coefficient-f′′(0) and local Nusselt number -θ′(0) for various values of m for the case of κ=0, M=0, Ec=0, Pr=0.73, Nb=10-5, Nt=0, and Le=0.
- f ′ ′ ( 0 ) - θ ′ ( 0 ) m Present Ashwini [20] Watanaba [21] Ullah [22] Present Watanaba [21] 0.0000 0.46960 0.4696 0.46960 0.4696 0.42016 0.42015 0.0141 0.50461 0.5046 0.50461 0.5046 0.42578 0.42578 0.0435 0.56898 0.5690 0.56898 0.5690 0.43548 0.43548 0.0909 0.65498 0.6550 0.65498 0.6550 0.44730 0.44730 0.1429 0.73200 0.7320 0.73200 0.7320 0.45694 0.45693 0.2000 0.80213 0.8021 0.80213 0.8021 0.46503 0.46503 0.3333 0.92765 0.9277 0.92765 0.9277 0.47814 0.47814 1.0000 1.23258 1.2326  1.2326Table2 illustrates the influence of nondimensional governing parameters on the skin friction coefficient, local Nusselt, and local Sherwood numbers. The skin friction coefficient enhances with increase in pressure gradient, permeability, and magnetic parameter. The local Nusselt number is a decreasing function and a local Sherwood number is an increasing function of the pressure gradient parameter, magnetic parameter, Permeability parameter, Prandtl number, Lewis number, thermophoresis parameter, and Eckert number.Table 2
Computations of the skin friction coefficient-f′′(0), local Nusselt number -θ′(0) and local Sherwood number -ϕ′(0) for various parameters.
β M κ Pr Ec Nb Nt Le -f′′(0) -θ′(0) -ϕ′(0) 0.25 2.0 0.5 0.5 0.5 0.4 0.2 1.5 1.64283184 0.20409760 0.56485047 0.5        1.64925659 0.20245352 0.56792266 1.0        1.66202020 0.19918041 0.57392554  0.25 1.0       1.35213817 0.22650827 0.52246831 3.0       1.88981423 0.18295134 0.60015609 5.0       2.30622396 0.14414670 0.65848609  2.0 0.2      1.56122953 0.21069258 0.55304965 0.5      1.64283184 0.20409760 0.56485047 0.8      1.72061427 0.19761936 0.57603379  0.5 0.1     1.64283184 0.24670503 0.31761537 0.5     1.64283184 0.20409760 0.56485047 1.0     1.64283184 0.09226748 0.81573784  0.5 0.1    1.64283184 0.35278900 0.49744039 0.5    1.64283184 0.20409760 0.56485047 1.0    1.64283184 0.01794788 0.64924507  0.5 0.2   1.64283184 0.22555655 0.58070313 0.4   1.64283184 0.20409760 0.56485047 0.8   1.64283184 0.16405080 0.55624299  0.4 0.1  1.64283184 0.21129537 0.54685694 0.2  1.64283184 0.20409760 0.56485047 0.4  1.64283184 0.19012517 0.60871601  0.2 1.0 1.64283184 0.20846325 0.47924548 1.5 1.64283184 0.20409760 0.56485047 2.0 1.64283184 0.20105976 0.63575168Figure2(a) shows the variation of velocity profiles for different values of pressure gradient parameter β. It clearly demonstrates that the velocity profile increases with an increase in pressure gradient parameter. Because of the increment of wedge angle, the fluid moves much slower and decreases velocity boundary layer thickness. Figure 2(b) shows the effect of permeability parameter κ on the velocity profile. It is observed that increase in κ leads to increase the velocity of the nanofluid on the porous surface and decrease its boundary layer thickness. It is also noticed that both pressure gradient parameter β and permeability parameter κ have no significant effect on both nanofluid temperature and concentration.Figure 2
(a) Velocity profiles for various values ofβ. (b) Velocity profiles for various values of κ.Figures3(a), 3(b), and 4(a) illustrate the influences of magnetic parameter M on the velocity, temperature, and concentration profiles, respectively. Figure 3(a) reveals that the velocity boundary layer thickness decreases with an increase in magnetic parameter. This is due to the fact that the presence of transverse magnetic field sets in Lorentz force, which results in retarding force on the velocity field. Consequently, as the values of magnetic parameter increase, so does the retarding force and hence the velocity profile increase. Figure 3(b) shows that the thermal boundary layer thickness decreases with an increase in magnetic parameter. This is due to additional work expended in dragging the fluid in the boundary layer against the action of the Lorentz force and energy is dissipated as thermal energy which heats the fluid. This reduces the temperature. It is also observed that the concentration profile and its boundary layer thickness decrease with an increase in magnetic parameter as shown in Figure 4(a). Figure 4(b) descries the concentration profile for different values of Lewis number Le. It is clearly observed that the concentration profile and its boundary layer thickness reduce considerably as the Lewis number increases.Figure 3
(a) Velocity profiles for various values of M. (b) Temperature profiles for various values of M.Figure 4
(a) Concentration profiles for various values of M. (b) Concentration profiles for various values of Le.The effect of the Prandtl number on the temperature and concentration is shown in Figures5(a) and 5(b), respectively. It is depicted that the temperature and concentration profiles and their boundary layer thickness reduce significantly as the Prandtl number increase. Because increasing the Prandtl number tends to reduce the thermal diffusivity of the fluid and causes weak penetration of heat inside the fluid. However, in the region near to the boundary surface, the heat transfer rate increases with an increase in Pr. This is due to the fact that the temperature gradient at the surface increase.Figure 5
(a) Temperature profiles for various values of Pr. (b) Concentration profiles for various values of Pr.The effect of viscous dissipation parameter Ec on the temperature and concentration profiles is presented in Figures6(a) and 6(b), respectively. The Eckert number expresses the conversion of kinetic energy into internal energy by work done against the viscous fluid stress. It is observed that the temperature increases significantly from the surface and attains a peak value around η=0.5 and then decreases in the rest of the region as given in Figure 6(a). This implies that the thermal boundary layer becomes thicker with large Eckert number. The concentration profile gradually reduces near the surface up to η=1 and then it increases with an increase in the viscous dissipation parameter Ec as highlighted in Figure 6(b).Figure 6
(a) Temperature profiles for various values of Ec. (b) Concentration profiles for various values of Ec.The influence of the Brownian motion parameter Nb on the nanofluid temperature and concentration profiles is presented in Figures7(a) and 7(b), respectively. Figure 7(a) reveals that the temperature profile increases with an increase in Nb, particularly in the region close to the surface. The physics behind this phenomenon is that the increased Nb increases the thickness of thermal boundary layer, which finally enhances the temperature. Figure 7(b) remarks that an increase in the values of Nb tends to decrease the concentration profile near the surface. Figures 8(a) and 8(b), respectively, reveal the usual temperature and concentration profiles for various values of thermophoresis parameter Nt. The thermophoresis force generated by the temperature gradient produces a fast flow and more fluid is heated away from the surface. Consequently, the higher the value of Nt increase temperature profile and its boundary layer thickness as given in Figure 8(a). Figure 8(b) reveals that concentration profile declines near the boundary surface until η=0.5 and afterward it increases with an increase in thermophoresis parameter Nt.Figure 7
(a) Temperature profiles for various values of Nb. (b) Concentration profiles for various values of Nb.Figure 8
(a) Temperature profiles for various values of Nt. (b) Concentration profiles for various values of Nt.


## 5. Conclusions

The problem of two-dimensional steady laminar MHD boundary layer wedge flow with heat and mass transfer of nanofluid past a porous media with viscous dissipation, Brownian motion, and thermophoresis effects has been studied. Using suitable similarity transformations, the governing equations are transformed to a system of nonlinear ordinary differential equations and solved numerically employing spectral quasilinearization method. From the above discussions the following conclusions are given:(1)
The thickness of velocity boundary layer reduces with an increase in pressure gradient, permeability, and magnetic parameters.(2)
Thermal boundary layer thicker with an increase in Eckert number, Brownian motion, and thermophoresis parameters.(3)
Greater values of Prandtl number, Lewis number, Brownian motion, and magnetic parameter reduce the nanofluid concentration profile.(4)
The skin-friction coefficient at the surface enhances with an increase in pressure gradient, permeability, and magnetic parameters.(5)
The local Nusselt number is a decreasing function, but a local Sherwood number is an increasing function of the pressure gradient parameter, magnetic parameter, permeability parameter, thermophoresis parameter, and Eckert number.



---

*Source: 4507852-2019-01-30.xml*