1. Introduction {#sec1}
===============

The study of non-Newtonian fluids has attracted the attention of new researchers due to their huge range of practical applications in engineering and science. These fluids show a nonlinear relationship between the stress and the rate of strain which offer much complicated and nonlinear differential systems. Among the non-Newtonian fluids, Casson fluid attracted more attention of researchers due to its applications in the manufacturing of pharmaceutical products, coal in water, china clay, paints, synthetic lubricants, and biological fluids. Casson [@bib1] introduced the first Casson model for the behavior of fluid flow in pigment-oil suspensions. Casson fluid characteristics are different from Newtonian fluids due to the relationships between shear stress and the rate of shear strain. Human blood, Jelly, honey, concentrated fruit juices, and tomato sauce are some examples of Casson fluid. Recently, Khan et al. [@bib2] studied Casson fluid flow with homogeneous-heterogeneous reactions. The collective effects of fluid dynamics, thermodynamics and, electrodynamics give rise to the topic MFD. It is divided into two parts: MHD and MGD. MHD deals with electrically conducting fluids whereas MGD deals with ionized compressible gases. MHD has several applications, namely, fusion research, MHD accelerator, and power generator. Some recent studies concerning the effect of magnetic field on fluid flow can be mentioned through studies \[[@bib3], [@bib4], [@bib5]\]. Convection term is arising the mixing of one portion to another portion of the fluid due to gross moments of the bulk of the fluid. There are three types of convection mode which are divided in to free, force and mixed in fluid flow problems. Recently, the free convective MHD flow plays a significant role in the petrochemical industry, heat exchanger design, geophysics and, MHD power generation system. Kataria and Patel [@bib6] discussed unsteady free convective MHD flow of micropolar fluid. Selimefendigil and Öztop \[[@bib7], [@bib8]\] studied magnetic field effects on free convective nanofluid flow in a cavity. The concept of Mixed convection in fluid flow problems has attracted many researchers due to its importance in many engineering applications such as food processing, microelectronic devices, cooling of electronic devices and, nuclear reactors. Selimefendigil and. Öztop \[[@bib9], [@bib10], [@bib11], [@bib12], [@bib13], [@bib14], [@bib15]\] considered mixed convection MHD flow of nanofluid with different physical conditions. Hayat et al. [@bib16] performed three-dimensional mixed convection MHD flow with variable viscosity whereas, Waqas et al. [@bib17] studied mixed convection MHD flow due to the nonlinear stretched sheet. Currently, many researchers have scrutinized the influence of thermal radiation on the heat transport features of fluids across stretched surfaces. The heat transfer process has many applications in solar power technology, and photochemical reactors, etc. Shah et al. [@bib18] and Sheikholeslami et al. [@bib19] discussed unsteady MHD flow with thermal radiation. Hayat et al. [@bib20] performed solar radiation and Joule heating in the MHD flow of thixotropic nanofluid. Non-linear thermal radiation effects are also important in heat transfer fluid flow problems. The application of heat transfer problems with non-linear thermal radiation is summarized as gas cooled nuclear reactors, counting combustion, furnace design and, nuclear reactor protection. Kumar et al. \[[@bib21], [@bib22]\] considered MHD flow over a convective surface with non-linear thermal radiation effects. The porous medium is a hard containing random spaces (pores), joined or separate, distributed within it in either a fixed or random manner. If pores signify a certain portion of the bulk volume, a difficult complex network can be designed which is capable to carry fluids. The porous media play an important role in technology and engineering field. The significant of the MHD flow problems with a porous medium which is done by Kataria and Patel [@bib23]. Research work on MHD Casson fluid flow in a porous medium with radiation can be found in Das [@bib24]. The study of chemical reaction procedures is beneficial for improving a number of chemical technologies, such as food processing, polymer production, and manufacturing of ceramics or glassware. In these procedures, the molecular diffusion via chemical reaction within or at the boundary cannot be determined. Abbas et al. [@bib25] obtained a numerical solution for chemical reaction effects on stagnation point flow of Casson fluid over a stretching/shrinking sheet with thermal radiation. Analytical solution of MHD Casson fluid flow past an oscillating plate in the presence of chemical reaction and thermal radiation obtained by Kataria and Patel [@bib26]. Further, the phenomenon of heat generation/absorption plays a significant role in adjusting the heat transfer rate. The effects of heat generation and heat absorption on MHD flow with different physical situations mentioned in Refs. \[[@bib27], [@bib28]\]. Some other relevant studies with thermal radiation, chemical reaction and, heat generation effects on MHD flow with different fluids can be found in Refs. \[[@bib29], [@bib30]\]. The Soret and Dufour effect arises when simultaneous heat and mass transfer affecting each other. Due to the concentration gradient, heat transfer process produced which cause the Dufour effect whereas, the mass transfer produced by temperature gradient is called the Soret effect. Kataria and Patel [@bib31] considered Soret effects on free convective Casson fluid flow whereas, Sheikholeslami et al. [@bib32] defined analytical solution for the effect of thermal diffusion and heat-generation on MHD nanofluid flow in a porous medium. The cross diffusion effects are usually of a small order of magnitude. These combine effect shows an important role in the process of solar ponds, biological systems, and, the microstructure of oceans. The cross diffusion effects were used in free or mixed convection MHD flow with heat and mass transfer problems applications. Relevant research work is done in Refs. \[[@bib33], [@bib34], [@bib35], [@bib36], [@bib37], [@bib38]\]. Moreover, Kataria and Mittal [@bib39], Kumar et al. [@bib40], Kumaran and Sandeep [@bib41] studied the impact of Brownian motion and thermophoresis on MHD flow with different physical conditions.

The elementary concepts of the homotopy in topology to recommend analytic approximation method, namely the HAM. This technique has been used to solve different types of nonlinear differential equation successfully. The HAM is better compared to numerical methods because it does not include discretization of the variable, hence is free from rounding off errors. The homotopy analysis method was introducing first by Liao [@bib42]. The HAM method can be applied Liao \[[@bib43], [@bib44]\] for solving several kinds of a nonlinear system of differential equations problems. The HAM theory is applied effectively to a nonlinear problem arise from the real situation in lots of different fields by researchers, engineers and, scientists. Fallahzadeh and Shakibi [@bib45] to solve the convection-diffusion equation using the homotopy analysis method. Recently, Fagbade et al. [@bib46] obtained an analytical approximation of thermal radiation and heat generation effects on MHD flow of viscoelastic fluid using HAM methods.

The main theme of the present article is to develop the model of cross diffusion effects on MHD Casson fluid flow through a porous medium. A novel aspect of chemical reaction, heat generation and, nonlinear thermal radiation is present. The focal concern of this problem is to investigate the collective effect of cross-diffusion, heat generation, nonlinear radiation and chemical reaction on two-dimensional MHD flow of Casson fluid through a porous medium which is not yet studied. The governing dimensionless systems of ordinary differential equations with imposing initial and boundary conditions are solved using the Homotopy analysis method. The convergence of series solutions is useful for providing the developing linear functions of solutions which is an advantage of the HAM method. Application of these types of research work can be found in fire dynamics in insulations, solar collection systems and, recovery of petroleum products, etc.

2. Model {#sec2}
========

A graphical explanation of the physical problem of steady two dimensional laminar mixed convective MHD flow of Casson fluid past an infinite plate in a porous medium is depicted in [Fig. 1](#fig1){ref-type="fig"}. The uniform ambient temperature $T_{\infty}$ and ambient concentration $C_{\infty\ }$ is considered. The values of velocity, temperature and concentration distribution of the sheet to be considered as $u_{w}\left( x \right),\ T_{w}$ and $C_{w\ }$ respectively. The flow is kept in the region $\left( y\  > \ 0 \right)$ toward a surface according. Coordinate system chosen in such a way that, x-axis along the sheet and the y-axis perpendicular to it. In y-axis direction, the constant strength $B_{0}$ is applied.Fig. 1Physical sketch of the problem.Fig. 1

The constitutive equation for the Casson fluid can be written as$$\text{τ}_{\text{ij}} = \left\{ \begin{matrix}
{2\left( {\text{μB} + \frac{\text{Py}}{\sqrt{2\text{π}}}} \right)\text{e}_{\text{ij}}\quad\pi > \text{π}_{\text{c}}} \\
{2\left( {\text{μB} + \frac{\text{Py}}{\sqrt{2\text{π}_{\text{c}}}}} \right)\text{e}_{\text{ij}}\quad\pi < \text{π}_{\text{c}}} \\
\end{matrix} \right.$$Where, $e_{ij}$ Deformation rate at $\left( {\text{i},\ \text{j}} \right)^{\text{th}}$ component, $Py$ Yield stress, $\mu B$ Dynamic viscosity, $\pi$ Product of the component for deformation rate with itself and, $\pi_{c}$ is a critical value of this product based on the non-Newtonian model.

The viscous dissipation and induce magnetic are negligible effects. The governing continuity, momentum, energy and concentration equations under the Boussinesq′s approximation are given below,$$\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} = 0$$$$u\frac{\partial u}{\partial x} + v\frac{\partial u}{\partial y} = \nu\left( {1 + \frac{1}{\gamma}} \right)\frac{\partial^{2}u}{\partial y^{2}} - \sigma B_{0}^{2}u - \frac{\mu\varphi}{k_{1}}u + \rho g\beta_{T}\left( T - T_{\infty} \right) + \rho g\beta_{C}\left( {C - C_{\infty}} \right)$$$$u\frac{\partial T}{\partial x} + v\frac{\partial T}{\partial y} = \frac{k_{3}}{\rho c_{p}}\frac{\partial^{2}T}{\partial y^{2}} - \frac{1}{\rho c_{p}}\frac{\partial q_{r}}{\partial y} + \frac{DK_{T}}{c_{s}c_{p}}\frac{\partial^{2}C}{\partial y^{2}} + \frac{Q_{0}}{\rho c_{p}}\left( {T - T_{\infty}} \right)\ $$$$u\frac{\partial C}{\partial x} + v\frac{\partial C}{\partial y} = D\frac{\partial^{2}C}{\partial y^{2}} + \frac{DK_{T}}{T_{m}}\frac{\partial^{2}\theta}{\partial y^{2}} - K_{2}\left( {C - C_{\infty}} \right)$$$$\ \ u = u_{w} = ax,\ \ v = 0,\ \ T = T_{w},\ \ C = C_{w}\ \ at\ y = 0.$$$$\left. u\rightarrow 0,\ T\rightarrow T_{\infty},\ \ \ C\rightarrow C_{\infty};\ as\ \ y\rightarrow\infty \right.$$Here, $u$ velocity in the x-direction, $v$ velocity in the y-direction, $T$ Fluid temperature, $C\ $ Concentration, $C_{p}$ Specific heat at constant pressure, $q_{r}$ Radiative heat flux, $c_{s}$ Concentration susceptibility, $D$ Mass diffusion coefficient, $K_{T}\ $ thermal diffusion ratio, $\sigma$ Electrical conductivity, $\rho$ Fluid density, $\varphi$ Porosity of the porous medium, $\beta_{T}$ Thermal expansion coefficient, $\beta_{C}\ $ Concentration expansion coefficient, $g\ $ Acceleration due to gravity, $\nu\ $ Kinematic viscosity coefficient and, $k_{3}\ $ Thermal conductivity.

Using Rosseland [@bib48], radiation heat flux $q_{r}$ can be written as$$q_{r} = - \frac{16\sigma^{\ast}}{3k^{\ast}}\ T^{3}\frac{\ \partial T}{\partial y}$$where, $\sigma^{\ast}\ and\ k^{\ast}\ $ the Stefan--Boltzmann constant and the mean absorption coefficient respectively.

In view to Eqs. [(8)](#fd8){ref-type="disp-formula"}, [Eq. (4)](#fd4){ref-type="disp-formula"} reduced to$$u\frac{\partial T}{\partial x} + v\frac{\partial T}{\partial y} = \frac{k_{3}}{\rho c_{p}}\frac{\partial^{2}T}{\partial y^{2}} + \frac{16\sigma^{\ast}}{3k^{\ast}}\frac{1}{\rho c_{p}}\left\lbrack {T^{3}\frac{\ \partial^{2}T}{\partial y^{2}} + 3T^{2}\left( \frac{\ \partial T}{\partial y} \right)^{2}} \right\rbrack + \frac{DK_{T}}{c_{s}c_{p}}\frac{\partial^{2}C}{\partial y^{2}} + \frac{Q_{0}}{\rho c_{p}}\left( {T - T_{\infty}} \right)\ $$

Introducing non dimensional variables$$\psi = \left( \sqrt{\nu\ xu_{w}} \right)f\left( \eta \right),\ T = T_{\infty}\left( {1 + \left( {\theta_{w} - 1} \right)\theta} \right),\ \ \ C\left( \eta \right) = \frac{\text{C} - C_{\infty}}{C_{w} - C_{\infty}}\ ,\eta = \sqrt{\frac{u_{w}}{\nu}}\ y,\ \ \ u = \frac{\partial\psi}{\partial y},\ v = - \frac{\partial\psi}{\partial x}\ $$where $f$ dimensionless stream function, $u_{w}$ Stretching velocity and $\eta$ Dimensionless coordinator.

Therefore, the dimensionless form of governing equations is expressed as,$$\left( {1 + \frac{1}{\gamma}} \right)\ f^{\prime\prime\prime} + \frac{3}{4}f^{\prime\prime}f - \frac{1}{2}{f^{\prime}}^{2} + \lambda\left( {\theta + NC} \right) - \left( {M^{2} + \frac{1}{k}} \right)f^{\prime} = 0$$$$\theta^{\prime\prime} + Nr\left\lbrack {\left( {1 + \left( {\theta_{w} - 1} \right)\theta} \right)^{3}\theta^{\prime\prime} + 3\left( {\theta_{w} - 1} \right){\theta^{\prime}}^{2}\left( {1 + \left( {\theta_{w} - 1} \right)\theta} \right)^{2}} \right\rbrack + \frac{3Pr}{4}\theta^{\prime}f + Pr.\ H\theta + PrD_{f}C^{\prime\prime} = 0$$$$\frac{C^{\prime\prime}}{Sc} + \frac{3}{4}f\ C^{\prime} + Sr\theta^{\prime\prime} - KrC = 0$$

Subject to$$f = 0,\ f^{\prime} = 1,\ C = 1,\ \theta = 1\mspace{9mu}\ at\mspace{9mu}\ \eta = 0$$$$\left. \ f^{\prime} = 0,\ \theta = 0,\ C = 0\ \mspace{9mu} at\mspace{9mu}\ \eta\rightarrow\infty \right.$$where, $\lambda = \frac{Gr_{x}}{R{e_{x}}^{2}}$ Mixed convection parameter, $N = \frac{G{r_{x}}^{\ast}}{T_{x}}$ Buoyancy force parameter, $Gr_{x} = \frac{g\beta_{T}\left( {T_{w} - T_{\infty}} \right)}{\nu^{2}}$ Thermal grashof number, $G{r_{x}}^{\ast} = \frac{g\beta_{C}\left( {C_{w} - C_{\infty}} \right)}{,\nu^{2}}$ Mass grashof number, $Re_{x} = \frac{x\ u_{w}}{\nu}$ Local Reynolds number, $M^{2} = \ \frac{\sigma{B_{0}}^{2}}{\rho b}$ Magnetic parameter, $Nr = \frac{16\sigma^{\ast}}{3k^{\ast}}\frac{{T^{3}}_{\infty}}{k_{1}}$ Radiation parameter, $Pr = \frac{\rho c_{p}}{k_{3}}$ Prandtl number, $Sr = \frac{DK_{T}\left( {T_{w} - T_{\infty}} \right)}{v\left( {C_{w} - C_{\infty}} \right)}$ Soret number, $D_{f} = \frac{\left( {C_{w} - C_{\infty}} \right)\ DK_{T}}{vc_{s}c_{p}\left( {T_{w} - T_{\infty}} \right)}$ Dufour number, $Sc = \frac{v}{D}$ Schmidt number, $\frac{1}{k} = \frac{\mu\varphi}{k_{1}}$ permeability of porous medium, $Kr = \frac{K_{2}}{\nu^{2}}$ chemical reaction parameter, $H = \frac{Q_{0}v^{2}}{\rho c_{p}}$ Heat generation/absorption parameter.

Expressions of skin-friction, Nusselt number, and Sherwood number are defined using the relations.$$R{e_{x}}^{\frac{1}{2}}C_{f_{x}} = \left( {1 + \frac{1}{\gamma}} \right)f^{\prime\prime}\left( 0 \right)$$$$R{e_{x}}^{- \frac{1}{2}}Nu_{x} = - \left( {1 + Nr{\theta_{w}}^{3}} \right)\theta^{\prime}\left( 0 \right)\ $$$$\ R{e_{x}}^{- 1/2}Sh_{x} = - C^{\prime}\left( 0 \right)$$Where $Re_{x} = u_{w}\frac{x}{\nu}$ is local Reynolds number.

3. Methodology {#sec3}
==============

The HAM is a fundamental idea of topology, proposed by Liao [@bib42]. Homotopy analysis methods is an analytic approximation method with a guarantee of convergence, mainly for nonlinear differential equations. In the application point of view, the HAM has been applied to solve fluid mechanics problems. In this study, the governing equations [(10)](#fd10){ref-type="disp-formula"} -- (12) are coupled non-linear ordinary differential equations with initial and boundary conditions (13)--(14) whose exact solutions are not possible. Therefore, HAM suggested by Liao \[[@bib42], [@bib43], [@bib44]\] is suitable to solve these types of system. Here initial guess and auxiliary linear operators are required for HAM solution which is presented given below:$$f_{0}\left( \eta \right) = 1 - e^{- \eta};\ \theta_{0}\left( \eta \right) = e^{- \eta};\ g_{0}\left( \eta \right) = e^{- \eta}$$$$L_{f} = \frac{\partial^{3}f}{\partial\eta^{3}} - \frac{\partial f}{\partial\eta}\ ,\ \ L_{\theta} = \ \frac{\partial^{2}\theta}{\partial\eta^{2}} - \theta\ ,\ L_{C} = \frac{\partial^{2}C}{\partial\eta^{2}} - \ C\ $$

Satisfying$$L_{f}\left( {C_{1} + C_{2}e^{\eta} + C_{3}e^{- \eta}} \right) = 0,\ \ \ L_{\theta}\left( {C_{4}e^{\eta} + C_{5}e^{- \eta}} \right) = 0,\ \ \ \ L_{C}\left( {\ C_{6}e^{\eta} + C_{7}e^{- \eta}} \right) = 0.$$where $c_{1},\ \ c_{2},\ldots.,\ c_{8}$ are constants.

The deformation of zero^th^ order can be written as,$$\left( {1 - p} \right)L_{f}\left\lbrack {\mspace{1mu}\widehat{f}\left( {\eta;\ p} \right) - f_{0}\left( \eta \right)} \right\rbrack = p\hslash_{f}N_{f}\left\lbrack {\widehat{f}\left( {\eta;\ p} \right),\widehat{\theta}\left( {\eta;\ p} \right),\ \widehat{C}\left( {\eta;\ p} \right)} \right\rbrack,$$$$\left( {1 - p} \right)L_{\theta}\left\lbrack {\widehat{\theta}\left( {\eta;\ p} \right) - \theta_{0}\left( \eta \right)} \right\rbrack = p\hslash_{\theta}N_{\theta}\left\lbrack {\widehat{f}\left( {\eta;\ p} \right),\widehat{\theta}\left( {\eta;\ p} \right),\ \widehat{C}\left( {\eta;\ p} \right)} \right\rbrack,$$$$\left( {1 - p} \right)L_{C}\left\lbrack {\widehat{C}\left( {\eta;\ p} \right) - C_{0}\left( \eta \right)} \right\rbrack = p\hslash_{C}N_{C}\left\lbrack {\widehat{f}\left( {\eta;\ p} \right),\widehat{\theta}\left( {\eta;\ p} \right),\ \widehat{C}\left( {\eta;\ p} \right)} \right\rbrack,$$

Subject to the boundary conditions:$$\widehat{f}\left( {0;\ p} \right) = 0,\ {\widehat{f}}^{\prime}\left( {0;\ p} \right) = 1;$$$${\widehat{f}}^{\prime}\left( {\infty;\ p} \right) = 0\text{;}$$$$\widehat{\theta}\left( {0;\ p} \right) = 1,\ \widehat{\theta}\left( {\infty;\ p} \right) = 0,\ $$$$\widehat{C}\left( {0;\ p} \right) = 1,\ \widehat{C}\left( {\infty;\ p} \right) = 0;\ $$

The nonlinear operator is defined as$$N_{f}\left\lbrack {\widehat{f}\left( {\eta;\ p} \right),\widehat{\theta}\left( {\eta;\ p} \right),\widehat{C}\left( {\eta;\ p} \right)} \right\rbrack = \left( {1 + \frac{1}{\beta}} \right)\ \frac{\partial^{3}\widehat{f}}{\partial\eta^{3}} + \frac{3}{4}{\widehat{f}}^{\prime\prime}\widehat{f} - \frac{1}{2}{\widehat{f}}^{\prime 2} + \lambda\left( {\widehat{\theta} + N\widehat{C}} \right) - \left( {M^{2} + \frac{1}{k}} \right){\widehat{f}}^{\prime}$$$$N_{\theta}\left\lbrack {\widehat{f}\left( {\eta;\ p} \right),\widehat{\theta}\left( {\eta;\ p} \right),\widehat{C}\left( {\eta;\ p} \right)} \right\rbrack = \ \ \ \frac{\partial^{2}\widehat{\theta}}{\partial\eta^{2}} + Nr\left\lbrack \left( {1 + \left( {{\widehat{\theta}}_{w} - 1} \right)\widehat{\theta}} \right)^{3}{\widehat{\theta}}^{\prime\prime} + \left\lbrack {\left( {1 + \left( {\theta_{w} - 1} \right)\theta} \right)^{3}\theta^{\prime\prime} + 3\left( {{\widehat{\theta}}_{w} - 1} \right){\widehat{\theta}}^{\prime 2}\left( {1 + \left( {{\widehat{\theta}}_{w} - 1} \right)\widehat{\theta}} \right)^{2}} \right\rbrack + \frac{3Pr}{4}{\widehat{\theta}}^{\prime}f + Pr\ H\widehat{\theta} + PrD_{f}{\widehat{C}}^{\prime\prime} \right\rbrack$$$$N_{C}\left\lbrack {\widehat{f}\left( {\eta;\ p} \right),\widehat{\theta}\left( {\eta;\ p} \right),\widehat{C}\left( {\eta;\ p} \right)} \right\rbrack = \frac{1}{Sc}{\widehat{C}}^{\prime\prime} + \frac{3}{4}f\ {\widehat{C}}^{\prime} + Sr{\widehat{\theta}}^{\prime\prime} - Kr\widehat{C}$$where $\widehat{f}\left( {\eta;\ p} \right),\ \widehat{\ \theta}\left( {\eta;\ p} \right)\ $ and $\widehat{\ C}\left( {\eta;\ p} \right)$ are unknown functions. $\hslash_{f}$, $\hslash_{\theta}$ and $\hslash_{C}\ $ are auxiliary parameters and the nonlinear operators are $\ N_{f}$, $N_{\theta}$ and $N_{C}$.

Also $p \in \left( {0,\ \ 1} \right)$ is an implanting parameter.

For $p = 0$ and $p = 1$ we have$$\widehat{f}\left( {\eta;0} \right) = f_{0}\left( \eta \right),\ \widehat{f}\left( {\eta;\infty} \right) = f\left( \eta \right),$$$$\widehat{\theta}\left( {\eta;0} \right) = \ \theta_{0}\left( \eta \right),\ \widehat{\theta}\left( {\eta;\infty} \right) = \ \theta\left( \eta \right),\ $$$$\widehat{C}\left( {\eta;0} \right) = C_{0}\left( \eta \right),\ \widehat{C}\left( {\eta;\infty} \right) = C\left( \eta \right),$$

In other words, when the variation of $p$ is taken from 0 to 1 is considered then $\ \widehat{f}\left( {\eta;\ p} \right)$, $\ \widehat{\theta}\left( {\eta;\ p} \right)$ and $\widehat{C}\left( {\eta;\ p} \right)$ are change from $f_{0}\left( \eta \right),\ \theta_{0}\left( \eta \right)$ and $C_{0}\left( \eta \right)$ to $f\left( \eta \right),\ \theta\left( \eta \right)$ and $C\left( \eta \right)$. The series expansion using Taylor′s formula of these functions yields the following,$$\widehat{f}\left( {\eta;\ p} \right) = f_{0}\left( \eta \right) + \sum\limits_{m = 1}^{\infty}{f_{m}\left( \eta \right)p^{m}}\text{,}$$$$\widehat{\theta}\left( {\eta;\ p} \right) = \theta_{0}\left( \eta \right) + \sum\limits_{m = 1}^{\infty}{\theta_{m}\left( \eta \right)p^{m}}\text{,}$$$$\widehat{C}\left( {\eta;\ p} \right) = C_{0}\left( \eta \right) + \sum\limits_{m = 1}^{\infty}{C_{m}\left( \eta \right)p^{m}}\text{,}$$Where$$f_{m}\left( \eta \right) = \frac{1}{m!}\left\lbrack \frac{\partial^{m}f\left( {\eta;\ p} \right)}{\partial p^{m}} \right\rbrack_{p = 0}\ ,$$$$\theta_{m}\left( \eta \right) = \frac{1}{m!}\left\lbrack \frac{\partial^{m}\theta\left( {\eta;\ p} \right)}{\partial p^{m}} \right\rbrack_{p = 0},$$$$C_{m}\left( \eta \right) = \frac{1}{m!}\left\lbrack \frac{\partial^{m}C\left( {\eta;\ p} \right)}{\partial p^{m}} \right\rbrack_{p = 0}\ ,$$

It is seen that the above series convergence depends upon $\hslash_{f},\ \hslash_{\theta}$ and $\ \hslash_{C}$. Assuming the nonzero auxiliary parameters so that the convergence of Eqs. [(21)](#fd21){ref-type="disp-formula"}, [(22)](#fd22){ref-type="disp-formula"}, and [(23)](#fd23){ref-type="disp-formula"} at $p = 1$ is appear. Hence it can be obtained the following,$$f\left( \eta \right) = f_{0}\left( \eta \right) + \sum\limits_{m = 1}^{\infty}{f_{m}\left( \eta \right)}\text{,}$$$$\theta\left( \eta \right) = \theta_{0}\left( \eta \right) + \sum\limits_{m = 1}^{\infty}{\theta_{m}\left( \eta \right)}\text{,}$$$$C\left( \eta \right) = C_{0}\left( \eta \right) + \sum\limits_{m = 1}^{\infty}{C_{m}\left( \eta \right)}\text{,}$$

Differentiating Eqs. [(21)](#fd21){ref-type="disp-formula"}, [(22)](#fd22){ref-type="disp-formula"}, [(23)](#fd23){ref-type="disp-formula"}, [(24)](#fd24){ref-type="disp-formula"}, [(25)](#fd25){ref-type="disp-formula"}, [(26)](#fd26){ref-type="disp-formula"} and [(27)](#fd27){ref-type="disp-formula"} m times with respect to $p$, putting $p = 0$ , and dividing by $m!$ , it is obtained the $m^{th}$ order deformation $\left( m \geq 1 \right)$ which can written as,$$L_{f}\left\lbrack {f_{m}\left( \eta \right) - \text{χ}_{m}f_{m - 1}\left( \eta \right)} \right\rbrack = \hslash_{f}R_{f,\ m}\left( \eta \right)\text{.}$$$$L_{\theta}\left\lbrack {\theta_{m}\left( \eta \right) - \text{χ}_{m}\theta_{m - 1}\left( \eta \right)} \right\rbrack = \hslash_{\theta}R_{\theta,\ m}\left( \eta \right)$$$$L_{C}\left\lbrack {C_{m}\left( \eta \right) - \text{χ}_{m}C_{m - 1}\left( \eta \right)} \right\rbrack = \hslash_{C}R_{C,\ m}\left( \eta \right)$$

Subject to the initial and boundary conditions$$f_{m}\left( 0 \right) = f_{m}^{\prime}\left( \infty \right) = 0,$$$$f_{m}^{\prime}\left( 0 \right) = 1,$$$$\theta_{m}\left( 0 \right) = 1,\ \ \theta_{m}\left( \infty \right) = 0,$$$$C_{m}\left( 0 \right) = 1,\ \ C_{m}\left( \infty \right) = 0,\ $$with$$R_{f,\ m}\left( \eta \right) = \left( {1 + \frac{1}{\beta}} \right)f_{m - 1}^{\prime\prime\prime} + \frac{3}{4}\sum\limits_{j = 0}^{m - 1}f_{j}f_{m - 1 - j}^{\prime\prime} - \frac{1}{2}\sum\limits_{j = 0}^{m - 1}{f_{m - 1 - j}^{\prime 2} + \lambda\left( \theta_{m - 1} + NC_{m - 1} \right)} - \left( {M^{2} + \frac{1}{k}} \right)f_{m - 1}^{\prime}$$$$R_{\theta,\ m}\left( \eta \right) = \theta_{m - 1}^{\prime\prime} + Nr\left\lbrack {\theta_{m - 1}^{\prime\prime} + \left( {\theta_{w}^{3} - 3\theta_{w}^{2} + 3\theta_{w} - 1} \right)}{\left( {\sum\limits_{j = 0}^{m - 1}{\theta_{m - 1 - j}\ }\sum\limits_{l = 0}^{j}\theta_{j - l}\sum\limits_{p = 0}^{l}{\theta_{l - p}\theta_{p}^{\prime\prime}}} \right) + \left( {\theta_{w}^{2} - 2\theta_{w} + 1} \right)\left( {3\sum\limits_{j = 0}^{m - 1}{\theta_{j}\theta_{m - 1 - j}^{\prime\prime}}} \right)}{+ \left( {3\sum\limits_{j = 0}^{m - 1}\theta_{m - 1 - j}\sum\limits_{l = 0}^{j}{\theta_{j - l}\theta_{l}^{\prime\prime}}} \right)\left( {\theta_{w} - 1} \right) + \left( {3\left( {\theta_{w} - 1} \right)} \right)\left\{ {\sum\limits_{j = 0}^{m - 1}{\theta_{j}^{\prime}\theta_{m - 1 - j}^{\prime\prime}} + 2\left( {\theta_{w} - 1} \right)\sum\limits_{j = 0}^{m - 1}\theta_{m - 1 - j}^{\prime}\sum\limits_{l = 0}^{j}{\theta_{j - l}^{\prime}\theta_{l}} + \left( {\theta_{w} - 1} \right)^{2}\sum\limits_{j = 0}^{m - 1}\theta_{m - 1 - j}\sum\limits_{l = 0}^{j}\theta_{j - l}\sum\limits_{p = 0}^{l}{\theta_{l - p}^{\prime}\theta_{p}^{\prime}}} \right\}} \right\rbrack + \frac{3Pr}{4}\sum\limits_{j = 0}^{m - 1}{\theta_{m - 1 - j}^{\prime}\, f_{j}} + H.Nr\ \theta_{m - 1} + PrD_{f}C_{m - 1 - j}^{\prime\prime}$$$$C_{c,\ m}\left( \eta \right) = \frac{1}{Sc}C_{m - 1}^{\prime\prime} + \frac{3}{4}\sum\limits_{j = 0}^{m - 1}{f_{j}C_{m - 1 - j}} + Sr\theta_{m - 1 - j}^{\prime\prime} - KrC_{m - 1}$$$$\text{with}\ \chi_{m\ } = \left\{ \begin{matrix}
{0,} & {m \leq 1} \\
{1,} & {m \geq 1} \\
\end{matrix} \right.\text{,}$$

Solving the corresponding *m*^th^-order deformation equations,$$f_{m}\left( \eta \right) = f_{m}^{\ast}\left( \eta \right) + C_{1} + C_{2}e^{\eta} + C_{3}e^{- \eta}$$$$\theta_{m}\left( \eta \right) = \theta_{m}^{\ast}\left( \eta \right) + C_{4}e^{\eta} + C_{5}e^{- \eta}$$$$C_{m}\left( \eta \right) = C_{m}^{\ast}\left( \eta \right) + C_{6}e^{\eta} + C_{7}e^{- \eta}$$

Here the particular solutions of m^th^ -- order equations are $f_{m}^{\text{*}}\ ,\ \theta_{m}^{\text{*}}\ and\ C_{m}^{\text{*}}$. Using boundary conditions, the constants $C_{i}\left( \ i = 1,2,\ \ldots,8 \right)$ are to be determined.

4. Analysis {#sec4}
===========

It is well known that convergence of the HAM solutions depends on the auxiliary parameters $\hslash_{F},\hslash_{\theta}$ and $\ \hslash_{C}$. For this concept, the allied h-curves are plotted in Figs. [2](#fig2){ref-type="fig"}, [3](#fig3){ref-type="fig"}, and [4](#fig4){ref-type="fig"} respectively. These all figures clearly suggest admissible ranges for the auxiliary parameters.Fig. 2h-curve of $f^{\prime\prime}\left( \eta \right)$.Fig. 2Fig. 3h-curve of $\theta^{\prime}\left( \eta \right)$.Fig. 3Fig. 4h-curve of $C^{\prime}\left( \eta \right)$.Fig. 4

5. Results & discussion {#sec5}
=======================

In order to get a rich understanding of the behavior of the problem, a parametric study is completed and the found results are clarified with the support of graphical diagrams. The non-dimensional fluid velocity, fluid temperature, and concentration for several values of the different physical parameter are presented in Figs. [5](#fig5){ref-type="fig"}, [6](#fig6){ref-type="fig"}, [7](#fig7){ref-type="fig"}, [8](#fig8){ref-type="fig"}, [9](#fig9){ref-type="fig"}, [10](#fig10){ref-type="fig"}, [11](#fig11){ref-type="fig"}, [12](#fig12){ref-type="fig"}, [13](#fig13){ref-type="fig"}, [14](#fig14){ref-type="fig"}, [15](#fig15){ref-type="fig"}, [16](#fig16){ref-type="fig"}, [17](#fig17){ref-type="fig"}, [18](#fig18){ref-type="fig"}, [19](#fig19){ref-type="fig"}, [20](#fig20){ref-type="fig"}, [21](#fig21){ref-type="fig"}, [22](#fig22){ref-type="fig"}, and [23](#fig23){ref-type="fig"}.Fig. 5$f^{\prime}\left( \eta \right)$ for different values of $M$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 5Fig. 6$f^{\prime}\left( \eta \right)$ for different values of $\gamma$ at $k = 0.1,\ Nr = 0.6,\ M = 0.5,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 6Fig. 7$f^{\prime}\left( \eta \right)$ for different values of $k$ at $M = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 7Fig. 8$f^{\prime}\left( \eta \right)$ for different values of $Nr$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 8Fig. 9$\theta\left( \eta \right)$ for different values of $Nr$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ $.Fig. 9Fig. 10$f^{\prime}\left( \eta \right)$ for different values of $H$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ M = 0.5\ .$Fig. 10Fig. 11$\theta\left( \eta \right)$ for different values of $H$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ M = 0.5\ \ .$Fig. 11Fig. 12$f^{\prime}\left( \eta \right)$ for different values of $- H$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ M = 0.5\ .$Fig. 12Fig. 13$\theta\left( \eta \right)$ for different values of $- H$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ H = 0.1,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ M = 0.5\ \ .$Fig. 13Fig. 14$f^{\prime}\left( \eta \right)$ for different values of $D_{f}$ at $k = 0.1,\ M = 0.5,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ H = 0.1,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ M = 0.5.$Fig. 14Fig. 15$f^{\prime}\left( \eta \right)$ for different values of $Sr$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ M = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1.$Fig. 15Fig. 16$C\left( \eta \right)$ for different values of $Sr$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ M = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 16Fig. 17$f^{\prime}\left( \eta \right)$ for different values of $Kr$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ M = 0.5\ and\ H = 0.1\ .$Fig. 17Fig. 18$C\left( \eta \right)$ for different values of $Kr$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ M = 0.5\ and\ H = 0.1\ \ .$Fig. 18Fig. 19$f^{\prime}\left( \eta \right)$ for different values of $N$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ M = 0.5,\ Kr = 0.1\ and\ H = 0.1\ \ .$Fig. 19Fig. 20$f^{\prime}\left( \eta \right)$ for different values of $\lambda$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ M = 0.5,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ \ .$Fig. 20Fig. 21$\theta\left( \eta \right)$ for different values of $\ \theta_{w}$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ M = 0.5,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 21Fig. 22$\theta\left( \eta \right)$ for different values of $Pr$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ M = 0.5,\ \theta_{w} = 0.1,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 22Fig. 23$C\left( \eta \right)$ for different values of $Sc$ at $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ M = 0.5,\ D_{f} = 0.5,\ Sr = 0.5,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ .$Fig. 23

[Fig. 5](#fig5){ref-type="fig"} shows magnetic field parameter $M$ effects on velocity profiles. Physically, the Lorentz force becomes stronger with an increase in $M$. Due to this reason, Magnetic parameter tends to the reduced motion of the fluid which is displayed in [Fig. 5](#fig5){ref-type="fig"}. Therefore, in fluid flow problems, hydrodynamics situation is better in comparison to the hydromagnetic case. The influence of Casson fluid parameter $\gamma$ on velocity profiles are shown in [Fig. 6](#fig6){ref-type="fig"}. The yield stress decreases and the thickness of the momentum boundary layer increases, when $\gamma$ is increased. Due to this reason velocity distribution decrease with an increase in the values of $\gamma$. Physically, the fluid becomes more viscous with increasing Casson fluid parameter $\gamma$. It is further observed from this graph that when $\left. \gamma\rightarrow\infty \right.$, the fluid behaviours change from non-Newtonian to Newtonian. [Fig. 7](#fig7){ref-type="fig"} shows velocity distribution for several values of permeability parameter $\ k$. It is seen that velocity distribution increase with increase in $k$. This phenomena is clearly satisfied which is explains the physical situation that as $k$ increases. These result are clearly supported with Kataria and Patel [@bib26]. Figs. [8](#fig8){ref-type="fig"} and [9](#fig9){ref-type="fig"} shows the effects of thermal radiation parameter $Nr$ on velocity and temperature distribution. It is seen that, thermal radiation parameter $Nr$ tends to improve velocity and temperature distribution. The thermal buoyancy force increase and thickness of the thermal boundary layer decrease, when $Nr$ is increased. Physically, Due to increasing thermal radiation parameter $Nr$, heat is generated in fluid flow, which leads to improvement in velocity and temperature distribution throughout the fluid flow region. The effects of Heat generation and Heat absorption parameter $H$ on temperature and velocity distribution are demonstrated in Figs. [10](#fig10){ref-type="fig"}, [11](#fig11){ref-type="fig"}, [12](#fig12){ref-type="fig"}, and [13](#fig13){ref-type="fig"}. The positive values of $H$ indicate the heat is generated whereas negative values means heat is absorption. Above all figures illustrate that velocity and temperature distribution increase with positive values of H (heat source), whereas negative values of H (heat absorbed) have reverse effects on it. These outcomes are strongly supported because the hall of porous is increase due to the effects of heat source which rises the temperature in the flow field.

The influence of Dufour number $\left( Df \right)$ and Soret number $\left( Sr \right)$ on velocity, temperature and, concentration distribution is shown in Figs. [14](#fig14){ref-type="fig"}, [15](#fig15){ref-type="fig"}, and [16](#fig16){ref-type="fig"}. According to the definition, the values of $Df$ and $Sr$ to be considered in such a way that their product is constant. From the figure, it is observed that, motion of the fluid elevates with increase in $Sr$ (decrease in $Df$). Physically, the mass buoyancy force increase due to increase the values of $Sr$, which leads to improve heat transfer process. It is also seen form [Fig. 14](#fig14){ref-type="fig"} that, temperature distribution and its boundary layer thickness heighten with increase in $Df$. Figs. [17](#fig17){ref-type="fig"} and [18](#fig18){ref-type="fig"} shows effects of reaction parameter on $Kr$ on velocity and concentration profiles. The values of chemical reaction $Kr\  > \ 0$ leads to reduce the concentration distribution. This occurrence has a superior agreement with the physical realities because of the weakens the buoyancy effects due to concentration gradients. Finally, it is concluded that chemical reaction tends to reduced velocity and concentration distribution. Effects of buoyancy force parameter N and Mixed convection parameter $\lambda$ on velocity distribution are discussed in Figs. [19](#fig19){ref-type="fig"} and [20](#fig20){ref-type="fig"} respectively. It is observed that, both parameter tends to improve velocity distribution throughout the flow field. Physically, $\lambda > 0$ indicates that heat is converted from the surface. Therefore, velocity of the fluid increases. The buoyancy parameter N performs a ratio of the buoyancy with the inertia of the external forced on the heat and fluid flow. These results indicate that the velocity distribution increase with buoyancy parameter. The impact of temperature ratio parameter $\theta_{w}$ on temperature distribution is shown in [Fig. 21](#fig21){ref-type="fig"}. From the graph, it is seen that the heat transfer process improve with temperature ratio parameter $\theta_{w}$. The results of Figs. [5](#fig5){ref-type="fig"}, [6](#fig6){ref-type="fig"}, [16](#fig16){ref-type="fig"}, [19](#fig19){ref-type="fig"}, and [20](#fig20){ref-type="fig"} are in good agreement with the outcomes of work of Kumar et al. [@bib36]. Figs. [22](#fig22){ref-type="fig"} and [23](#fig23){ref-type="fig"} shows effects of Prandtl number $Pr$ and Schmidt number $Sc$ on temperature and concentration distribution. It is depict that the temperature and concentration distribution decreases with increase in $Pr$ and $Sc$ respectively. This result is justified because the increase value in $Pr$, thermal conductivity of the fluid decrease, due to this reason the thermal boundary layer thickness is decreases.

[Table 1](#tbl1){ref-type="table"} shows the HAM convergence analysis for governing dimensionless momentum, energy and concentration equations. It is seen that 10^th^ order of approximation is enough for momentum equation, 15^th^ order approximation is adequate for energy and concentration equations. The variation of the Skin friction, Nusselt number, Sherwood number are shown in Tables [2](#tbl2){ref-type="table"} and [3](#tbl3){ref-type="table"} for various values of the governing parameters. It is illustrated that the permeability of porous medium $k$, buoyancy force parameter N and mixed convection parameter $\lambda$ tends to improve Skin friction coefficient whereas, Magnetic field M and Casson fluid parameter $\gamma$ have reverse effects on it. The local Nusselt number increase with increase in thermal radiation parameter Nr, heat generation H and temperature ratio parameter $\theta_{w}$. Sherwood number decrease with increase in chemical reaction parameter Kr and Dufour number $Df$. Thermal-diffusion parameter Sr tends to improved skin friction coefficient, local Nusselt number and Sherwood number. To check the validation of numerical approach, local Nusselt number is computed and compared with previous published work for some special case which is presented in [Table 4](#tbl4){ref-type="table"}. It is seen that, the Nusselt number ($- \theta^{\text{'}}\left( 0 \right)\ $ at $Nr = 0$) is strong agreement with Kumar et al. [@bib36] and Chen [@bib47].Table 1Convergence of the homotopic solutions for different order of approximations at M = 0.5, $k = 0.5,\ Nr = 0.6,\ \gamma = 1,\ \lambda = 1,\ \theta_{w} = 0.1,\ Pr = 5,\ Sr = 0.5,\ D_{f} = 0.5,\ Sc = 0.66,\ N = 1,\ Kr = 0.1\ and\ H = 0.1\ $Table 1Order of approximation$- f^{\prime\prime}\left( 0 \right)$$- \theta^{\text{'}}\left( 0 \right)\ $$- C^{\prime}\left( 0 \right)$21.05161.60970.233441.05341.64380.203661.05381.63970.205681.05391.63380.2085101.05381.63110.2098151.05381.63090.2098201.05381.63110.2097251.05381.63110.2097301.05381.63110.2097Table 2Skin friction coefficient, Local Nusselt number and Sherwood number variation at $Nr = 0.1,\ Sr = 2,\ D_{f} = 1,\ Sc = 0.22,\ Kr = 2,\ H = 0.5,\ \theta_{w} = 1.2,\ Pr = 15\ \text{and}\ \gamma = 0.5$.Table 2$M$$k$$N$$\lambda$$\theta_{w}$$Pr$$\left( {1 + \frac{1}{\gamma}} \right)f^{\prime\prime}\left( 0 \right)$$- \left( {1 + Nr\theta_{w}^{3}} \right)\theta^{\prime}\left( 0 \right)$$- C^{\prime}\left( 0 \right)$0.50.510.11.215−2.8489991.9140610.2096400.6−2.9055041.9157040.2089050.7−2.9709271.9175950.2080640.6−2.6706981.9088130.2120050.7−2.5359311.9047840.2138430.6−2.8417231.9138310.2097480.7−2.8344491.9136020.2098560.2−2.7814391.9121590.2104690.3−2.7138091.9102580.2112971.3−2.8491062.0031530.2056471.4−2.8492272.1086310.20087316−2.8491041.9032500.21231017−2.8491361.8901280.216024Table 3Skin friction coefficient, Local Nusselt number and Sherwood number variation at $M = 0.5,\ k = 0.5,\ N = 0.5,\ \lambda = 0.1,\ \theta_{w} = 1.2\ \text{and}\ Pr = 15$.Table 3$Nr$$Sr$$Df$$Sc$$Kr$$H$$\gamma$$\left( {1 + \frac{1}{\gamma}} \right)f^{\prime\prime}\left( 0 \right)$$- \left( {1 + Nr\theta_{w}^{3}} \right)\theta^{\prime}\left( 0 \right)$$- C^{\prime}\left( 0 \right)$0.1210.2220.50.5−2.8489991.9140610.2096400.2−2.8494852.2386770.1978270.3−2.8499792.5760700.1856022.1−2.8480361.8626670.2020472.2−2.8470931.8133600.1955491.1−2.8497511.8868520.2201771.2−2.8504541.85984250.23132550.3−2.8512981.8182720.2121630.4−2.8532381.7026330.1764212.1−2.8500231.9476440.2193862.2−2.8510181.98089650.2288000.6−2.8500752.0154360.17530320.7−2.8509552.1111420.1424450.6−2.6837651.9190380.2074160.7−2.5594481.9230650.205632Table 4Comparison of $\left( - \theta^{\text{'}}\left( 0 \right) \right)$ with "Kumar et al. [@bib36]" and "Chen [@bib47]" at.$Nr = 0,\ \theta_{w} = 1,\ H = 0\ \text{and}\ D_{f} = 0\ .$Table 4$Pr$$- \theta^{\text{'}}\left( 0 \right)\ $ Ref. [@bib36]$- \theta^{\text{'}}\left( 0 \right)\ $ Ref. [@bib47]HAM solution1.0−0.58223−0.58199−0.582243.0−1.16522−1.16523−1.165205.0−1.56803-−1.5680010.0−2.30798−2.30796−2.30790

6. Conclusion {#sec6}
=============

The problem of two-dimensional mixed convection MHD Casson fluid flow in a porous medium is examined by considering cross-diffusion, heat generation, chemical reaction, and nonlinear thermal radiation effects. The numerical solution is obtained by homotopy analysis method. The important results of these research work are listed below:➢The velocity, temperature and concentration distributions attain maximum value in the neighborhood of the plate and then decreased properly to approach the free stream value.➢The fluid flow accelerated via mixed convection parameter $\lambda$, buoyancy force parameter $N$, the permeability of porous medium parameter $k$, heat generation parameter H \> 0 and thermo-diffusion parameter $Sr$.➢A large value of Casson fluid parameter $\gamma$, Magnetic field parameter M, heat absorption parameter H \< 0 and chemical reaction parameter $Kr$ have retard effects on velocity distribution.➢Heat and mass transfer process raise with Dufour number $Df$ and Soret number $Sr$ respectively.➢Chemical reaction parameter Kr tends to reduced concentration distribution.➢The thermal radiation $Nr$ and heat generation $H\  > \ 0$ tends to increase temperature distribution.➢In most of the Engineering and real-world problems, skin friction is not desirable, so it can be decreased by decreasing the values of magnetic field parameter M, Casson fluid parameter $\gamma$ and chemical reaction parameter $Kr$.➢The rate of heat and mass transfer process improved with $Sr$, $Nr$ and $Kr$.

Declarations {#sec7}
============

Author contribution statement {#sec7.1}
-----------------------------

Harshad R Patel: Conceived and designed the analysis; Analyzed and interpreted the data; Contributed analysis tools or data; Wrote the paper.

Funding statement {#sec7.2}
-----------------

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Competing interest statement {#sec7.3}
----------------------------

The authors declare no conflict of interest.

Additional information {#sec7.4}
----------------------

No additional information is available for this paper.
