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.

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 Figure 1. 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 -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 , where is a positive constant. The porous medium has isotropic, homogeneous, and thermodynamic equilibrium with local fluid. Therefore, the governing equations are given by Subject to the boundary conditions 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:

It is presumed that the transformation in temperature within the flow is too small, so that the temperature term might be stated as a linear combination of temperature () about free stream temperature () using Taylor series and ignoring higher terms; we may express it as ; we get

We declare the ensuing nondimensional variables to nondimensionalise the governing equations:

Wherever the stream function () is defined in the normal notation as and , 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:By taking = porous-magnetic parameter and in (10) also considering , (9) and (10) are modified as follows: The corresponding boundary conditions are

The nondimensional numbers , , and are defined as follows: Substituting (7), (8), and (15) into (16), we receive the terms for the skin friction and local Nusselt and Sherwood numbers as follows:

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 asThe auxiliary linear operators , , and aresatisfying the following properties:where ( = 1–7) denote the arbitrary constants. The zeroth-order deformation problems are The -order deformation problem is of the formwhere The general solutions of (22)–(24) arewhere , , and are the special solutions. The symbolic calculations are obtained by MATLAB. Figure 2 indicates that the respective admissible values of , , and are and . We choose the values of auxiliary parameter from this range; we will get the more accurate results.

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 to by taking consecutive steps . The assumed initial values of , , and are compared with the measured values of , , and at . If the dissimilarity is present, assume another set of primary values for , , and and the procedure is repeated. This process is repeated up to the arrangement between considered value and the agreed condition at 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.

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 parameters , , , , and are fixed throughout the study. Tables 24 depict the values of , , and from the analytical and numerical methods. Table 2 depicts that skin friction increases on increasing the internal heat generation () and radiation (Rd) parameters. It is inferred that the skin friction decreases on increasing the chemical reaction (Cr) and slip () parameters. Tables 3 and 4 show that the heat transfer rate decreases and the mass transfer rate increases on increasing , Rd, and Cr parameters. It is noticed that the mass and heat transfer rates increase on increasing the slip parameter ().

Figure 3(a) illustrates profiles of velocity along the boundary layer for , , and 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 (), internal heat generation (), and slip parameter (). 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.

It can be observed from Figures 6(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.

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.

Nomenclature

Slip parameter
:Strength of the magnetic field
:Species concentration
Cr:Chemical reaction parameter
:Diffusion coefficient
:Dimensionless stream function
:Gravitational acceleration
:Solutal Grashof number
:Thermal Grashof number
:Thermal conductivity
:Permeability parameter
:Porous medium permeability
:Mean absorption coefficient
:Magnetic field parameter
:Navier slip coefficient
Nu:Nusselt number
Pr:Prandtl number
:Heat generation/absorption
Rd:Thermal radiation parameter
:Reynolds number
:Thermal Richardson number
:Solutal Richardson number
:Heat generation parameter
Sc:Schmidt number
Sh:Sherwood number
:Temperature
:Velocity components
:Cartesian coordinates.
Greek Symbols
:Coefficient of thermal expansion
:Coefficient of solutal expansion
:Chemical reaction rate
:Similarity variable
:Dimensionless temperature
:Dynamic viscosity
:Kinematic viscosity
:Density
:Stefan-Boltzmann constant
:Electrical conductivity
:Dimensionless concentration
:Stream function.
Subscripts
:At wall
:At free stream.

Competing Interests

The authors declare that they have no competing interests.