# Analytical and Numerical Study on Magnetoconvection Stagnation-Point Flow in a Porous Medium with Chemical Reaction, Radiation, and Slip Effects

Authors: H. Niranjan; S. Sivasankaran; M. Bhuvaneswari
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

We investigate the effects of slip and radiation on magnetoconvection flow of a chemically reacting fluid near a stagnation-point towards a vertical plate embedded in a porous medium analytically and numerically. The governing partial differential equations are diminished into the coupled ordinary differential equations by similarity transformations. Then they are solved analytically by homotopy analysis method and solved numerically by shooting method with RK fourth-order method. In this study, the analytical and numerical results are compared for many combinations of parameters. The rates of heat and mass transfer are calculated. The velocity profile near the plate overshoots on increasing the slip parameter. The concentration and temperature are decreasing on increasing the slip parameter.


---

## Body

## 1. Introduction

The study of convective heat transfer through porous medium for an incompressible fluid on the heated surface has received major attention because of its diverse uses in the insulation of nuclear reactors, petroleum industry, geothermal problems, storage of nuclear waste, and several other areas. The study of slip condition on convective boundary layer flow beside a vertical plate embedded in a porous medium has accepted considerable practical and theoretical awareness. The research has been maintained out in this area to analyze the heat and mass transfer characteristics within the boundary layer flow. Hiemenz [1] is the pioneer to investigate the forced convection two-dimensional flow at stagnation-point by similarity transformation. Many researchers extended Crane’s [2] problem of steady flow of an incompressible fluid over a linearly stretched plate to analyze several conditions. Mohammadi [3] numerically analyzed on boundary layer natural convection flow along a vertical plate. Cheng et al. [4] developed an analytic solution by homotopy analysis method for an unsteady mixed convection flow near the stagnation-point towards a vertical surface in a porous medium. Layek et al. [5] studied the boundary layer flow near stagnation-point in a porous medium over a stretching sheet with heat generation and blowing/suction. Lee et al. [6] studied convective heat transfer over a persuaded plate in a porous medium using Lie group analysis. They perceived that the thickness of the thermal and momentum boundary layer increases on increasing the radiation parameter. The natural convection flow over a plate in a porous medium was examined by Ferdows et al. [7], Rosali et al. [8], and Hou [9]. Bhuvaneswari et al. [10] studied the internal heat generation effect on natural convective flow over an inclined surface in a porous medium. They observed that the profiles of velocity and temperature increase on increasing the heat generation parameter.The research on magnetohydrodynamic (MHD) flow and heat transfer has several applications in engineering fields such as power generator, design of MHD accelerators, cooling of nuclear reactors, and heat exchangers. The effect of radiation on heat and mass transfer is essential in some applications, such as hydrometallurgical industries, solar power technology, and astrophysical flows. Takhar et al. [11] and Aydin and Kaya [12] explored the mixed magnetoconvection of a viscous fluid in the surrounding area of a stagnation-point next to a heated vertical surface. Hayat et al. [13] considered the mixed convection stagnation-point flow towards a sheet in a porous medium with magnetic field and radiation effects using homotopy analysis method. Makinde [14] studied the effects of internal heat generation and radiation on MHD mixed convection stagnation-point flow in a perpendicular plate embedded in a porous medium.The study of chemical reaction on heat and mass transfer is very significant in chemical technology and food processing. Bhuvaneswari et al. [15] studied the effect of chemical reaction on a semi-infinite inclined surface using Lie group analysis. Bhattacharyya [16] and Afify and Elgazery [17] investigated the two-dimensional flow of a viscous fluid with chemical reaction over a stretching/shrinking sheet near a stagnation-point. Recently, boundary layer flow with slip condition has attracted many researchers. Abel and Mahesha [18] analyzed the characteristics of MHD flow of a non-Newtonian viscoelastic fluid over a stretching sheet in the presence of nonuniform radiation and heat source. Hari et al. [19] numerically analyzed the effect of chemical reaction on MHD mixed convective flow over a vertical plate embedded in a porous medium with thermal radiation. The velocity and concentration increase on increasing the values of the chemical reaction parameter in the case of generative case and they are decreasing in the case of destructive case. Makinde and Sibanda [20] investigated the effects of chemical reaction and internal heat generation on the boundary layer flow over a linearly stretching sheet. The local skin friction and mass transfer rate increase on increasing the Schmidt number. Bhattacharyya and Layek [21] analyzed steady, laminar stagnation-point flow in a porous shrinking sheet with radiation and suction. Wang [22] and Labropulu and Li [23] examined the effect of slip of a non-Newtonian fluid behavior at a stagnation-point flow over a plate. Harris et al. [24] and Bhattacharyya et al. [25] explored the steady mixed convection flow at a stagnation-point in a porous medium over a perpendicular surface with velocity slip condition. Merkin et al. [26] and Rohni et al. [27] considered temperature slip effect on an unsteady boundary layer flow in a porous medium. Aman et al. [28] and Seini and Makinde [29] studied the MHD flow of an incompressible fluid over a linearly shrinking/stretching sheet with slip effects. Singh and Chamkha [30] examined the heat transfer and fluid flow near the vertical linear shrinking sheet using second-order slip condition. They noticed that the skin friction falls on increasing the slip parameter. Roşca and Pop [31] explored steady forced flow and convection heat transfer over a vertical shrinking/stretching sheet using slip condition.Inspired by the above applications and surveys explained, the purpose of this present study is to explore the effects of chemical reaction and velocity slip on mixed convection stagnation-point flow over a vertical plate embedded in a porous medium in the presence of external magnetic field and thermal radiation.


## 2. Mathematical Model

Consider a steady, 2D laminar, mixed convection boundary layer flow of an incompressible viscous fluid near a stagnation-point along a vertical plate in the presence of magnetic field through a porous medium as exposed in Figure1. It is presumed that the chemical reaction is considered as homogeneous and of first order. It is supposed that the porous medium is heat generating or absorbing internally at a fixed rate. Constant strength is enforced on unvarying magnetic field along the y-axis. The induced magnetic field is negligible due to small magnetic Reynolds number. As the fluid hits the wall at the stagnation-point, the flow is divided into two equal and opposite forces. In the potential flow region, the velocity distribution is assumed as U∞=cx, where c is a positive constant. The porous medium has isotropic, homogeneous, and thermodynamic equilibrium with local fluid. Therefore, the governing equations are given by (1)ux+vy=0,(2)uux+vuy=υuyy+gβ∗C-C∞+gβT-T∞-u-U∞σeB02ρ+υK~+U∞dU∞dx,(3)uTx+vTy=αTyy-1ρCp∂qr∂y+QT-T∞ρCp,(4)uCx+vCy=DCyy-Γ0C-C∞.Subject to the boundary conditions (5)u=N1uy,v=0,C=Cw,T=Tw,aty=0,u⟶U∞=cx,C⟶C∞,T⟶T∞,asy⟶∞,the medium is taken to be gray, and emitting radiation and heat flux in the equation of energy due to thermal radiation are given by means of the Rosseland approximation:(6)qr=-4σ∗3K′∂T4∂y.Figure 1
Schematic diagram of the problem.It is presumed that the transformation in temperature within the flow is too small, so that the temperature termT4 might be stated as a linear combination of temperature (T) about free stream temperature (T∞) using Taylor series and ignoring higher terms; we may express it as T4≈4T∞3T-3T∞4; we get (7)qr=-16σ∗T∞33K′Ty.We declare the ensuing nondimensional variables to nondimensionalise the governing equations:(8)b=N1cν1/2,Cr=Γ0c,GrC=gβ∗Cω-C∞x3υ2,GrT=gβTω-T∞x3υ2,K=υcK~,M=σeB02cρ,Pr=υα,Rd=4σ∗T∞3kK′,Rex=U∞xυ,RiC=GrCRex2,RiT=GrTRex2,S=QυρCpc,Sc=υD,η=ycυ,θη=T-T∞Tω-T∞,ϕη=C-C∞Cω-C∞,ψx,y=υcxfη.Wherever the stream function (ψ) is defined in the normal notation as u=ψy and υ=-ψx, then the continuity equation (1) is satisfied. Replacing the expression in (7) together with (8) into (1)–(5), we attain the nonlinear ordinary differential equations:(9)f′′′+ff′′-f′2+RiTθ+RiCϕ-K+Mf′-1+1=0,(10)1+43Rdθ′′+Prfθ′+Sθ=0,(11)ϕ′′+Scfϕ′-ScCrϕ=0.By taking (K+M)=Mp = porous-magnetic parameter and in (10) also considering 1+4/3Rd=RD, (9) and (10) are modified as follows: (12)f′′′+ff′′-f′2+RiTθ+RiCϕ-Mpf′-1+1=0,(13)RDθ′′+Prfθ′+Sθ=0.The corresponding boundary conditions are(14)f=0,f′=bf′′0,θ=1,ϕ=1,atη=0,f′=1,θ=0,ϕ=0,asη⟶∞.The nondimensional numbersCfx, Nu, and Sh are defined as follows: (15)Cfx=2τωρU∞2,Nu=xqωkTω-T∞,Sh=xqmDCω-C∞,(16)τω=μ∂u∂yy=0,qω=-k∂T∂yy=0-4σ∗3K′∂T4∂yy=0,qm=-D∂C∂yy=0.Substituting (7), (8), and (15) into (16), we receive the terms for the skin friction and local Nusselt and Sherwood numbers as follows:(17)Rex1/2Cfx=f′′0,Rex-1/2Nu=-RDθ′0,Rex-1/2Sh=-ϕ′0.


## 3. Method of Solution

### 3.1. Analytical Solution by HAM

The analytical solutions of (11)–(13) along with the boundary conditions (14) are obtained by homotopy analysis method (HAM). The initial approximations for homotopy analysis solutions are chosen as(18)f0η=η+e-η1+b-11+b;θ0η=e-η;ϕ0η=e-η.The auxiliary linear operators Lf, Lθ, and Lϕ are(19)Lf=f′′′-f′;Lθ=θ′′-θ,Lϕ=ϕ′′-ϕ,satisfying the following properties:(20)LfC1+C2eη+C3e-η=0;LθC4eη+C5e-η=0;LϕC6eη+C7e-η=0,where Ci(i = 1–7) denote the arbitrary constants. The zeroth-order deformation problems are (21)1-pLff-η,p-f0η=phfNff-η,p,θ-η,p,ϕ-η,p,1-pLθθ-η,p-θ0η=phθNθf-η,p,θ-η,p,ϕ-η,p,1-pLϕϕ-η,p-ϕη=phϕNϕf-η,p,θ-η,p,ϕ-η,p,f-0,p=0,f-′0,p=bf-′′0,p,f-′∞,p=1,θ-0,p=1,θ-∞,p=0,ϕ-0,p=1,ϕ-∞,p=0.The mth-order deformation problem is of the form(22)Lffmη-χmfm-1η=hfRmfη(23)Lθθmη-χmθm-1η=hθRmθη(24)Lϕϕmη-χmϕm-1η=hϕRmϕη(25)fm0=0,fm′0=bfm′′0,θm0=0,ϕm0=0,fm′∞=0,θm∞=0,ϕm∞=0,where (26)Rmfη=fm-1′′′η+∑k=0m-1fm-1-kηfk′′η-fm-1-k′ηfk′η+RiTθ+RiCϕ-Mpfm-1′η+Mp+11-χm,Rmθη=1+43Rdθm-1′′η+Pr∑k=0m-1fm-1-kηθk′η+Sθm-1η,Rmϕη=ϕm-1′′η+Pr∑k=0m-1fm-1-kηϕk′η-ScCrϕ,χm=0,m≤11,m>1.The general solutions of (22)–(24) are(27)fmη=f∗mη+C1+C2η+C3e-η,θmη=θ∗mη+C4eη+C5e-η,ϕmη=ϕ∗mη+C6eη+C7e-η,where f∗mη, θ∗mη, and ϕ∗mη are the special solutions. The symbolic calculations are obtained by MATLAB. Figure 2 indicates that the respective admissible values of hf, hθ, and hϕ are -0.5≤hf≤-0.2 and -0.6≤hθ,hϕ≤-0.2. We choose the values of auxiliary parameter (h=-0.35) from this range; we will get the more accurate results.Figure 2
h curves of f′′0, θ′(0), and ϕ′0.


### 3.2. Numerical Method

The coupled equations (11)–(13) along with the boundary conditions (14) are transformed into the linear differential equations by defining dependent variables of first order. The set of coupled linear differential equations of first-order system are integrated using Runge-Kutta fourth-order method from η=0 to η=ηmax by taking consecutive steps Δη. The assumed initial values of f′′(0), θ′(0), and ϕ′(0) are compared with the measured values of f′′(0), θ′(0), and ϕ′(0) at η=ηmax. If the dissimilarity is present, assume another set of primary values for f′′(0), θ′(0), and ϕ′(0) and the procedure is repeated. This process is repeated up to the arrangement between considered value and the agreed condition at η=ηmax within the degree of exactness. The analytical and numerical results are compared with the results available in the literature and they are shown in Table 1. This provides confidence in our analytical and numerical results to be reported subsequently.Table 1
Comparison of analytical and numerical values off′′, -θ′, and -ϕ′ for different values of S when RiT=1, RiC=0.5, Pr = 1, Sc = 0.5, M=0, K=0, Cr = 0, and Rd = 0 with Makinde [14].
S
f
′
′
0
-
θ
′
0
-
ϕ
′
0
HAM
RK4
Makinde [14]
HAM
RK4
Makinde [14]
HAM
RK4
Makinde [14]
−1
1.9164
1.8501
1.8444
1.1215
1.3897
1.3908
0.4744
0.4618
0.4631
0
2.0158
1.9989
1.9995
0.6333
0.6401
0.6932
0.4814
0.4761
0.4789
1
2.0705
2.1299
2.1342
−0.1055
−0.0727
−0.0730
0.5062
0.4902
0.4917


## 3.1. Analytical Solution by HAM

The analytical solutions of (11)–(13) along with the boundary conditions (14) are obtained by homotopy analysis method (HAM). The initial approximations for homotopy analysis solutions are chosen as(18)f0η=η+e-η1+b-11+b;θ0η=e-η;ϕ0η=e-η.The auxiliary linear operators Lf, Lθ, and Lϕ are(19)Lf=f′′′-f′;Lθ=θ′′-θ,Lϕ=ϕ′′-ϕ,satisfying the following properties:(20)LfC1+C2eη+C3e-η=0;LθC4eη+C5e-η=0;LϕC6eη+C7e-η=0,where Ci(i = 1–7) denote the arbitrary constants. The zeroth-order deformation problems are (21)1-pLff-η,p-f0η=phfNff-η,p,θ-η,p,ϕ-η,p,1-pLθθ-η,p-θ0η=phθNθf-η,p,θ-η,p,ϕ-η,p,1-pLϕϕ-η,p-ϕη=phϕNϕf-η,p,θ-η,p,ϕ-η,p,f-0,p=0,f-′0,p=bf-′′0,p,f-′∞,p=1,θ-0,p=1,θ-∞,p=0,ϕ-0,p=1,ϕ-∞,p=0.The mth-order deformation problem is of the form(22)Lffmη-χmfm-1η=hfRmfη(23)Lθθmη-χmθm-1η=hθRmθη(24)Lϕϕmη-χmϕm-1η=hϕRmϕη(25)fm0=0,fm′0=bfm′′0,θm0=0,ϕm0=0,fm′∞=0,θm∞=0,ϕm∞=0,where (26)Rmfη=fm-1′′′η+∑k=0m-1fm-1-kηfk′′η-fm-1-k′ηfk′η+RiTθ+RiCϕ-Mpfm-1′η+Mp+11-χm,Rmθη=1+43Rdθm-1′′η+Pr∑k=0m-1fm-1-kηθk′η+Sθm-1η,Rmϕη=ϕm-1′′η+Pr∑k=0m-1fm-1-kηϕk′η-ScCrϕ,χm=0,m≤11,m>1.The general solutions of (22)–(24) are(27)fmη=f∗mη+C1+C2η+C3e-η,θmη=θ∗mη+C4eη+C5e-η,ϕmη=ϕ∗mη+C6eη+C7e-η,where f∗mη, θ∗mη, and ϕ∗mη are the special solutions. The symbolic calculations are obtained by MATLAB. Figure 2 indicates that the respective admissible values of hf, hθ, and hϕ are -0.5≤hf≤-0.2 and -0.6≤hθ,hϕ≤-0.2. We choose the values of auxiliary parameter (h=-0.35) from this range; we will get the more accurate results.Figure 2
h curves of f′′0, θ′(0), and ϕ′0.


## 3.2. Numerical Method

The coupled equations (11)–(13) along with the boundary conditions (14) are transformed into the linear differential equations by defining dependent variables of first order. The set of coupled linear differential equations of first-order system are integrated using Runge-Kutta fourth-order method from η=0 to η=ηmax by taking consecutive steps Δη. The assumed initial values of f′′(0), θ′(0), and ϕ′(0) are compared with the measured values of f′′(0), θ′(0), and ϕ′(0) at η=ηmax. If the dissimilarity is present, assume another set of primary values for f′′(0), θ′(0), and ϕ′(0) and the procedure is repeated. This process is repeated up to the arrangement between considered value and the agreed condition at η=ηmax within the degree of exactness. The analytical and numerical results are compared with the results available in the literature and they are shown in Table 1. This provides confidence in our analytical and numerical results to be reported subsequently.Table 1
Comparison of analytical and numerical values off′′, -θ′, and -ϕ′ for different values of S when RiT=1, RiC=0.5, Pr = 1, Sc = 0.5, M=0, K=0, Cr = 0, and Rd = 0 with Makinde [14].
S
f
′
′
0
-
θ
′
0
-
ϕ
′
0
HAM
RK4
Makinde [14]
HAM
RK4
Makinde [14]
HAM
RK4
Makinde [14]
−1
1.9164
1.8501
1.8444
1.1215
1.3897
1.3908
0.4744
0.4618
0.4631
0
2.0158
1.9989
1.9995
0.6333
0.6401
0.6932
0.4814
0.4761
0.4789
1
2.0705
2.1299
2.1342
−0.1055
−0.0727
−0.0730
0.5062
0.4902
0.4917


## 4. Results and Discussion

Analytical solution by using HAM and numerical solution by using shooting method were obtained for different values of radiation, chemical reaction, and velocity slip parameters. The values of the parametersRiT=2.0, RiC=2.0, Mp=2.0, Pr=7.2, and Sc=1 are fixed throughout the study. Tables 2–4 depict the values of f′′(0), -θ′(0), and -ϕ′(0) from the analytical and numerical methods. Table 2 depicts that skin friction increases on increasing the internal heat generation (S) and radiation (Rd) parameters. It is inferred that the skin friction decreases on increasing the chemical reaction (Cr) and slip (b) parameters. Tables 3 and 4 show that the heat transfer rate decreases and the mass transfer rate increases on increasing S, Rd, and Cr parameters. It is noticed that the mass and heat transfer rates increase on increasing the slip parameter (b).Table 2
Analytical and numerical solutions for various values off′′(0) when K=1, M=1, Pr = 7.2, and Sc = 1.0.
S
Rd
Cr
b
f
′
′
0
Error (%)
HAM
RK4
−1.5
1.0
0.5
0.4
1.741935
1.742816
0.05055
−1.0
1.745303
1.746196
0.05114
0.0
1.748974
1.749733
0.04338
1.0
1.752804
1.753440
0.03627
1.5
1.756711
1.757330
0.03522
0.5
0.0
0.5
0.4
1.696237
1.697773
0.09047
0.3
1.718371
1.719537
0.06780
0.5
1.731982
1.731025
0.05528
0.8
1.744971
1.745330
0.02057
1.5
1.770112
1.770255
0.08078
2.0
1.782315
1.783635
0.07400
0.5
1.0
−2.0
0.4
2.061214
2.060439
0.03761
−1.5
1.925079
1.924475
0.03139
−0.5
1.805891
1.805358
0.02952
0
1.774912
1.775401
0.02754
0.5
1.752906
1.753440
0.03045
1.5
1.723012
1.722443
0.03303
2.0
1.711073
1.710801
0.01589
0.5
1.0
0.5
0.0
3.158507
3.157450
0.03347
0.2
2.273258
2.272768
0.02156
0.6
1.422014
1.421517
0.03496
0.8
1.194172
1.193269
0.07567
1.0
1.028569
1.027359
0.11777Table 3
Analytical and numerical solutions for various values of-θ′(0) when  K=1, M=1,  Pr = 7.2, and Sc = 1.0.
S
Rd
Cr
b
-
θ
′
0
Error (%)
HAM
RK4
−1.5
1.0
0.5
0.4
1.503127
1.502298
0.05518
−1.0
1.439975
1.439975
0.00000
0.0
1.375810
1.375810
0.00000
1.0
1.309663
1.309663
0.00000
1.5
1.241382
1.241382
0.00000
0.5
0.0
0.5
0.4
1.912091
1.911775
0.01653
0.3
1.646183
1.645690
0.02995
0.5
1.523071
1.522510
0.03684
0.8
1.383101
1.382877
0.01619
1.5
1.169021
1.169727
0.06035
2.0
1.068943
1.068377
0.05297
0.5
1.0
−2.0
0.4
1.445169
1.444713
0.03156
−1.5
1.387088
1.386636
0.03259
−0.5
1.334001
1.333199
0.06015
0
1.319987
1.319596
0.02963
0.5
1.309897
1.309663
0.01786
1.5
1.296015
1.295783
0.01790
2.0
1.291136
1.290635
0.03881
0.5
1.0
0.5
0.0
1.018166
1.017375
0.07775
0.2
1.209972
1.209753
0.01810
0.6
1.370195
1.369418
0.05673
0.8
1.409185
1.408859
0.02313
1.0
1.437043
1.436746
0.02067Table 4
Analytical and numerical solutions for various values of-ϕ′(0) when K=1, M=1, Pr = 7.2, and Sc = 1.0.
S
Rd
Cr
b
-
ϕ
′
0
Error (%)
HAM
RK4
−1.5
1.0
0.5
0.4
1.044089
1.043630
0.04398
−1.0
1.045321
1.044155
0.11166
0.0
1.045172
1.044706
0.04460
1.0
1.046004
1.045285
0.06878
1.5
1.046148
1.045894
0.02428
0.5
0.0
0.5
0.4
1.036874
1.036006
0.08378
0.3
1.039876
1.039445
0.04146
0.5
1.042083
1.041357
0.06971
0.8
1.044115
1.043834
0.02692
1.5
1.049270
1.048402
0.08279
2.0
1.051298
1.050986
0.02968
0.5
1.0
−2.0
0.4
−1.559173
−1.564133
0.31710
−1.5
−0.489101
−0.490619
0.30940
−0.5
0.513916
0.514218
0.05872
0
0.808952
0.808147
0.09961
0.5
1.046013
1.045285
0.06964
1.5
1.425416
1.424446
0.06809
2.0
1.585091
1.584309
0.04935
0.5
1
0.5
0.0
0.949062
0.948406
0.06916
0.2
1.012067
1.011676
0.03864
0.6
1.066135
1.065652
0.04532
0.8
1.079006
1.079208
0.01871
1.0
1.089222
1.088848
0.03434Figure3(a) illustrates profiles of velocity along the boundary layer for Rd=0.3, S=1.5, and Cr=0.3 and different values of slip parameter. The velocity increases as the slip parameter is increased. The boundary layer thickness increases when the slip parameter is increased. Figures 3(b) and 3(c) show the influence of slip parameter on temperature as well as concentration, respectively. These two profiles decrease on increasing the slip parameter. This is the reason why the viscosity decreases in the medium or the free stream velocity increases, and it results in the thermal and solutal boundary layer thickness decrease. The velocity and temperature profiles for various values of radiation parameter (Rd) are illustrated in Figures 4(a) and 4(b). As the radiation increases, the velocity as well as temperature increases significantly. The thickness of the thermal boundary layer increases on increasing the radiation parameter because the intensity of the electromagnetic radiation in the medium decreases and the Stefan-Boltzmann constant of the medium increases. Figures 5(a) and 5(b) depict the chemical reaction effect on both the velocity and concentration profiles along the surface for constant values of radiation (Rd=1.0), internal heat generation (S=1.0), and slip parameter (b=0.4). The velocity and concentration increase on increasing the chemical reaction parameter. The concentration profile overshoots in the generative chemical reaction case, that is, the negative values of the chemical reaction parameter.Figure 3
Velocity (a), temperature (b), and concentration (c) profiles for different values of slip parameter withS=1.5, Rd=0.3, and Cr=0.3.
(a)
(b)
(c)Figure 4
Velocity (a) and temperature (b) profiles for different values of radiation parameter withS=1.0, Cr=0.3, and b=0.4.
(a)
(b)Figure 5
Velocity (a) and concentration (b) profiles for different values of chemical reaction parameter withS=1.0, Rd=1.0, and b=0.4.
(a)
(b)It can be observed from Figures6(a)–6(c) that the local skin friction decreases on increasing the chemical reaction parameter and slip parameter. As the radiation increases, local skin friction increases due to temperature gradient causes in the fluid. In Figures 7(a)–7(c), it is interesting to note that Nusselt number increases on increasing the slip and radiation parameters. However, as the chemical reaction parameter increases, the local Nusselt number decreases. It is interesting to note that the local heat transfer rate decreases on increasing the internal heat generation parameter. Figures 8(a)–8(c) show that the local Sherwood number increases on increasing the chemical reaction, radiation, and slip parameters. Further scrutinising these figures, it is found that the local mass transfer rate increases with slip parameter. However, the effect of internal heat generation on local mass transfer rate is negligible.Figure 6
Local skin friction for different (a) chemical reaction parameter withS=0.5 and Rd=1.0, (b) radiation parameter with S=0.5 and Cr=0.5, and (c) slip parameter with Cr=0.5 and Rd=1.0.
(a)
(b)
(c)Figure 7
Local Nusselt numbers for different (a) slip parameter withCr=0.5 and Rd=1.0, (b) radiation parameter with S=0.5 and Cr=0.5, and (c) chemical reaction parameter with S=0.5 and Rd=1.0.
(a)
(b)
(c)Figure 8
Local Sherwood numbers for different (a) chemical reaction parameter withS=0.5 and Rd=1.0, (b) radiation parameter with S=0.5 and Cr=0.5, and (c) slip parameter with Cr=0.5 and Rd=1.0.
(a)
(b)
(c)


## 5. Conclusions

The effects of chemical reaction and slip parameters on magnetoconvection stagnation-point flow in a porous medium over a vertical plate with thermal radiation are investigated. The governing equations are solved analytically using HAM and numerically by shooting technique along with Runge-Kutta fourth-order integration. The results show a good agreement between the analytical and numerical results of the considered problem. The velocity increases on increasing the slip and radiation parameters and it increases on decreasing the chemical reaction parameter. The temperature decreases on increasing the slip parameter and it increases on increasing the radiation parameter. The concentration increases on decreasing the chemical reaction parameter and it decreases on increasing the slip parameter. The skin friction decreases on increasing the slip and chemical reaction parameters. The local mass and heat transfer rates increase on increasing the slip parameter.



---

*Source: 4017076-2016-04-26.xml*