Influence of Cattaneo – Christov Heat Flux on MHD Jeffrey , Maxwell , and Oldroyd-B Nanofluids with Homogeneous-Heterogeneous Reaction

This research article deals with the determination of magnetohydrodynamic steady flow of three combile nanofluids (Jefferey, Maxwell, and Oldroyd-B) over a stretched surface. The surface is considered to be linear. The Cattaneo–Christov heat flux model was considered necessary to study the relaxation properties of the fluid flow. The influence of homogeneous-heterogeneous reactions (active for auto catalysts and reactants) has been taken in account. The modeled problem is solved analytically. The impressions of the magnetic field, Prandtl number, thermal relaxation time, Schmidt number, homogeneous–heterogeneous reactions strength are considered through graphs. The velocity field diminished with an increasing magnetic field. The temperature field diminished with an increasing Prandtl number and thermal relaxation time. The concentration field upsurged with the increasing Schmidt number which decreased with increasing homogeneous-heterogeneous reactions strength. Furthermore, the impact of these parameters on skin fraction, Nusselt number, and Sherwood number were also accessible through tables. A comparison between analytical and numerical methods has been presented both graphically and numerically.


Introduction
A fluid composed of nanoparticles is called nanofluid.Nanoparticles of materials such as metallic oxides, carbide ceramics, nitride metals, ceramics, semiconductors, single, double or multi walled carbon nanotubes, alloyed, nanoparticles, etc. have been used for the preparation of nanofluids.Nanofluids have many characteristics in heat transfer, including microelectronics, Symmetry 2019, 11, 439 2 of 18 local refrigerator, cooler, machining, and heat exchanger.The idea of nanofluid was introduced by Choi [1].The Fourier's [2] recommended law of heat conduction normally works for heat transmission features from the time it was presented in the literature.By including the relaxation time parameter Cattaneo [3] it has improved this law and this term overwhelms the paradox of heat conduction.Christov [4] has named this theory Cattaneo-Christov heat flux theory, by further modifying the Cattaneo theory by exchanging the time derivative with Oldroyd-B upper convicted derivative.Mustafa [5] scrutinized the model [4] for heat transmission in a rotating Maxwell nanofluid flow.Chen [6] probed the influence of heat transfer and viscous dissipation of nanofluid flow over a stretching sheet.Sheikholeslami et al. [7][8][9][10][11] deliberated the three-dimensional magnetohydrodynamics (MHD) nanofluid flow in parallel rotating plates.Sheikholeslami [11][12][13][14][15] analytically and numerically deliberated the applications of nanofluids with different properties, behavior, and influences.Dawar et al. [16] examined the flow Williamson nanofluid over a stretching surface.Shah et al. [17] examined the micropolar nanofluid flow in rotating parallel plates with Hall current impact.Maleki et al. [18] scrutinized the non-Newtonian nanofluids flow and heat transfer over a porous surface.Nasiri et al. [19] deliberated the smoothed particle by a hydrodynamics approach for numerical simulation of nanofluid flows.Rashidi et al. [20] used the nanofluids in a circular tube heat exchanger and examined the entropy generation.Safaei et al. [21] studied numerically and experimentally the nanofluids convective heat transfer in closed conduits.Mahian et al. [22,23] presented the advances in the simulation and modeling of the flows of nanofluids.
Due to its relaxation properties, Jeffrey, Maxwell, and Oldroyd-B nanofluids have significant applications in the area of fluid mechanics.Ahmad et al. [24] scrutinized the flow of Jeffrey nanofluid with Magnetohydrodynamic impact.Ahmad and Ishak [25] deliberated the flow of Jeffrey nanofluid with MHD and transverse magnetic field impacts in a porous medium.Hayat et al. [26] probed the Oldroyd-B nanofluid flow with heat transfer and thermal radiation impacts.Raju et al. [27] deliberated the flow of Jeffrey nanofluid with a homogenous-heterogeneous reaction and non-linear thermal radiation impacts.The articles that are related to Jeffrey nanofluid can be found in [28][29][30][31][32]. Hayat et al. [33] inspected the MHD Maxwell nanofluid flow using suction/injection.Raju et al. [34] presented the heat and mass transmission in three-dimensional non-Newtonian nanofluid and Ferrofluid.Sandeep and Sulochana [35] investigated the mixed convection micropolar nanofluid flow over a stretching sheet.Raju et al. [36] deliberated the impacts of an inclined magnetic field, thermal radiation and cross diffusion on the two-dimensional flow.Nadeem et al. [37] presented the heat and mass transfer in Jeffrey nanofluid.Makinde et al. [38] deliberated the unsteady fluid flow with convective boundary conditions.Sheikholeslami [39] discussed the hydro-thermal behavior of nanofluids flow because of its external heated plates.Shah et al. [40] presented the Darcy-Forchheimer flow of radiative carbon nanotubes in a rotating frame.Chai et al. [41] presented a review of the heat transfer and hydrodynamic characteristics of nano/microencapsulated phase.Shah et al. [42] examined the electro-magneto micropoler Casson Ferrofluid over a stretching/shrinking sheet.Dawar et al. [43] analyzed the MHD CNTs Casson nanofluid in rotating channels.Khan et al. [44] deliberated the Williamson nanofluid flow over a linear stretching surface.Imtiaz et al. [45] examined the unsteady MHD flow due to a curved stretchable surface with homogeneous-heterogeneous reactions.Hayat et al. [46] deliberated the flow of nanofluids with homogeneous-heterogeneous reaction impacts over a non-linear stretched sheet with variable thickness.The recent study about nanofluid with application can be seen [47][48][49][50].
The present work is based on an analysis of MHD flow of three combine nanofluids (Maxwell, Oldroyd-B, and Jeffrey) over a linear stretching surface.The present model composed of Cattaneo-Christov heat flux.The impact of homogeneous-heterogeneous reactions were taken in this model.A boundary layer methodology was used in the mathematical expansion.The impact of dimensionless parameters on the fluid flow have been presented through graphs and tables.

Mathematical Modeling and Formulation
The incompressible electrically conducted three combined nanofluids (Jeffrey, Maxwell, and Oldroyd-B) were confined by a linear stretched surface.The fluid flow was taken in a two-dimensional steady state with stable surface temperature.x-axis was considered parallel to the surface, while y-axis was orthogonal to x-axis in the chosen coordinate system.The stretching velocity in x-axes direction was defined as U w (x) = ζx.The conclusion of homogeneous-heterogeneous reactions on the fluid flows of two chemical species I and J were taken in account.In the y-axis direction a uniform magnetic field B 0 was acting.The heat transmission procedure was applied through Cattaneo-Christov heat flux theory.
In case of cubic autocatalysis, the Homogeneous reaction is [45,46] While on the catalyst surface, the heterogeneous reaction has been defined by where k c , k s , I, J, i, j are the rate constants, chemical species, and concentrations of chemical species, respectively.
In the absence of viscous dissipation and thermal radiation, the boundary layer equations leading to the flow of viscoelastic fluids can be written as follows [47]: ρc p uT x + vT y = −∇.q,(5) Here u, v, µ, ρ f , υ f are velocity components in their respective directions, dynamic viscosity, density, and kinematic viscosity respectively.λ 1 , λ 2 , λ 3 are the relaxation time, a proportion of the relaxation to retardation times, respectively.T, σ f , B 0 indicated the temperature, electrical conductivity and the transverse magnetic field.
Symmetry 2019, 11, 439 4 of 18 The heat transfer equation proceeded as: where Φ E is given as: The accompanying boundary conditions were: where α = k ρc p indicated the thermal diffusivity, D I and D J indicated the diffusion coefficients, T w denoted the temperature at the surface, T ∞ for the surrounding fluid temperature and ζ for non-negative stretching rate constant with T −1 as the dimension.
Apparently the equation of continuity is satisfied and Equations ( 4)-( 13) become: with boundary conditions: In the above equations, M = Now Equations ( 16) and ( 17) yield: The subjected boundary conditions are: Skin friction coefficient through the dimensionless scale is: where Re x is called the local Reynolds number.
The dimensionless form of Nu x and Sh x were found as:
The initial assumptions were picked as below: The linear operators were taken as L F , L G and L φ : With the following properties: where r i (i = 1 − 7) were the constants: The resulting non-linear operators N F , N G and N φ were specified as: The zero th -order problem for Equations ( 14), ( 15) and (20) were: The related boundary conditions where: Symmetry 2019, 11, 439 where τ ∈ [0, 1] is the embedding parameter, F , G , φ that were used to control the solution convergence.When τ = 0 and τ = 1 we have: Expanding F(η; τ), G (η; τ) and φ(η; τ) by Taylor's series: where: The F , G and φ are taken in such a way that the series (35) converges at τ = 1, we have: The following are satisfied by the q th -order problem.
Which have the following boundary conditions: Here Symmetry 2019, 11, 439 7 of 18 where:

HAM Solution Convergences
In this segment we graphically discussed the superior effect of the concerned parameters.The convergence of Equation ( 36) was subjected entirely through the auxiliary constraints F , G , φ .This is a collection in a way that it controls and converges the series solutions.The optional division of , was plotted through -curves F (0), G (0), φ (0) for the 2 nd ordered approximated solution of HAM.The operational region of is −2.2 < F < 0.2, −2.1 < G < −0.1, −2.4 < φ < 0.1.The convergence of HAM through the -curve on velocity profile F (0), temperature profile G (0) and concentration profile φ (0) is presented in Figure 1.

Results and Discussion
In this segment the impact of emerging parameters on velocity function () G  and concentration function ()  within the defined domain have been discussed.
The impact of  on () 2. The Lorentz force theory deliberated that the magnetic field grows at a reversed force to the fluids flow.This force reduced the momentum boundary layer while it improved the thickness of the boundary layer.Therefore, with the escalating magnetic field  the velocity profile () F   declined.From here we concluded that Jeffrey nanofluid was greatly subjected by the magnetic field compared to the other two.In Figures 3 and 4 the impact of Pr and  on () G  were presented respectively.In Figure 3 we perceived that ()

G 
diminished with the rise in Pr .Physically the thickness of the boundary layer increased with the reduction in thermal diffusion.In addition, it can also be seen from the figure that Pr is more effective on Jeffrey and Maxwell nanofluids compared to the Oldroyd-B nanofluid.In Figure 4 the effect of thermal relaxation parameter  on () G  has been described.From here we saw that () G  reduced with the escalation in  .This was attributable to the fact that as we escalate  , the material particles need more time for heat transmission to its nearest particles.In addition, it can be stated that this material shows a non-conducting behavior with higher values of  which results in a reduction in ()  are schemed in Figures 5-7 respectively.In Figure 5 the effect of Sc on ()  has been described.Schmidt number is the ratio of momentum diffusivity to Sc .In Figure 6 the impact of  on ()  has been described.From here we concluded that

Results and Discussion
In this segment the impact of emerging parameters on velocity function F (η), temperature function G(η) and concentration function φ(η) within the defined domain have been discussed.The impact of M on F (η) is deliberated in Figure 2. The Lorentz force theory deliberated that the magnetic field grows at a reversed force to the fluids flow.This force reduced the momentum boundary layer while it improved the thickness of the boundary layer.Therefore, with the escalating magnetic field M the velocity profile F (η) declined.From here we concluded that Jeffrey nanofluid was greatly subjected by the magnetic field compared to the other two.In Figures 3 and 4 the impact of Pr and Ω on G(η) were presented respectively.In Figure 3 we perceived that G(η) diminished with the rise in Pr.Physically the thickness of the boundary layer increased with the reduction in thermal diffusion.In addition, it can also be seen from the figure that Pr is more effective on Jeffrey and Maxwell nanofluids compared to the Oldroyd-B nanofluid.In Figure 4 the effect of thermal relaxation parameter Ω on G(η) has been described.From here we saw that G(η) reduced with the escalation in Ω.This was attributable to the fact that as we escalate Ω, the material particles need more time for heat transmission to its nearest particles.In addition, it can be stated that this material shows a non-conducting behavior with higher values of Ω which results in a reduction in G(η).The impact of Sc, K and K s on φ(η) are schemed in Figures 5-7 respectively.In Figure 5 the effect of Sc on φ(η) has been described.Schmidt number is the ratio of momentum diffusivity to mass diffusivity.Physically, the Schmidt number is related to hydrodynamic layer's thickness and boundary layer.Sc intensifies the momentum of the boundary layer flow which results in an increase in concentration profile.It is clear from the figure that φ(η) upsurges with the rise in Sc.In Figure 6 the impact of K on φ(η) has been described.From here we concluded that larger K results in a reduction in φ(η).This may have been caused by the fact that the reaction rates controlled the diffusion coefficients.To a certain extent similar results are displayed in Figure 7.In Figure 7 the impact of K s on φ(η) has been described.From this figure we have concluded that the growing values of K s showed a drop in behavior in φ(η).This results from an agreement with the general physical behavior of the homogeneous reaction K and the heterogeneous reaction K s .In Figures 8 and 9 the impact of M on C f x and Nu x for Jeffrey, Maxwell and Oldroyd-B nanofluids have been described.It is clear from the figures that the growing values of M were decreasing for both C f x and Nu x .The magnetic field was applied perpendicular to the flow of fluids and had an inverse variation with the skin friction of the fluid flow.This is the reason why the increasing magnetic field reduced the skin friction of the fluids flow.Similarly, the behavior of the heat transfer rate was due to the growing magnetic force on the fluids flow phenomena, with the fluid particles requiring more time to transfer the heat to the nearest particle.This was because the heat transfer rate reduced with the escalating magnetic field.The impact of Pr and Ω on Nu x for the nanofluids flow has been described in Figures 10 and 11.From here we have concluded that the escalation in Pr increased the heat transfer rate while the increased Ω reduced the heat transfer rate for the nanofluids flow.Figure 12 shows the Total Residual error for the three types of nanofluid flow.

Tables Discussion
In this section we have demonstrated the effect of emerging dimensionless parameters on the presented model of nanofluids.Table 1 displayed

Tables Discussion
In this section we have demonstrated the effect of emerging dimensionless parameters on the presented model of nanofluids.Table 1 displayed the conclusions associated with emerging parameters of skin fraction coefficients.This shows the impression of magenatic field parameter M on skin fraction coefficients.The magenatic field parameter shows a reduction in the skin fraction coefficient.Table 2 demonstrated the conclusion of incipient parameters on local Nusselt numbers.The heat transfer rate decreases with the rise in thermal relaxation parameter Ω while escalates with the increase in Prandtl number Pr.Table 2 shows that the thermal relaxation parameter has more effect on Jefferey nanofluids in comparison to Maxwell and Oldroyd-B nanofluids.Table 3 demonstrated the conclusion of an emerging parameter on the Sherwood number.The Sherwood number reduces with its rise, which upsurges with the escalation of the strength of homogeneous reaction K and the strength of heterogeneous reaction K s .

Comparison of Analytical Solutions and Numerical Solutions
An analytical solution means an exact solution.To study the behavior of systems, an analytical solution can be used with varying properties.Regrettably there are many practical systems that lead to an analytical solution, and analytical solutions are often of limited use.This is why we use a numerical approach to generate answers that are closer to practical results.These solutions which cannot be used as complete mathematical expressions are numerical solutions.In the natural worldthere are almost no problems that are exactly solvable, which makes the problem more difficult than all the exactly solvable problems.There are three or four of them in nature that have already been solved, unfortunately even numerical methods cannot always give an exact solution.Numerical techniques can handle any completed physical geometries which are often impossible to solve analytically.In this article both analytical and numerical approaches are tested to solve the modeled problem.A comparison of HAM and ND-Solve technique for F (η), G(η) and φ(η) are deliberated in Figures 13-15 and Tables 4-6, respectively.
solution can be used with varying properties.Regrettably there are many practical systems that lead 247 to an analytical solution, and analytical solutions are often of limited use.This is why we use a 248 numerical approach to generate answers that are closer to practical results.These solutions which 249 cannot be used as complete mathematical expressions are numerical solutions.In the natural 250 worldthere are almost no problems that are exactly solvable, which makes the problem more difficult 251 than all the exactly solvable problems.There are three or four of them in nature that have already

Conclusions
In this article the MHD flow of three combined nanofluids (Jefferey, Maxwell, and Oldroyd-B) over a linear stretched surface have been scrutinized.The problem was solved analytically by HAM.
The convergence of HAM has been presented through graphical presentations.The concluding remarks are as follows: The upsurges in magnetic field diminishes the velocity field.The upsurges in Prandtl number and thermal relaxation parameters diminish the temperature field.The upsurges in Schmidt number upsurges the concentration field.The larger homogeneous reaction and heterogeneous reaction strengths falloff from the concentration field.
Author Contributions: A.S. and Z.S. modeled the problem and wrote the manuscript.S.I., A.D. and P.K. thoroughly checked the mathematical modeling and English corrections.W.K. and A.S. solved the problem using Mathematica software.Z.S., S.I. and P.K. contributed to the results and discussions.All authors finalized the manuscript after its internal evaluation.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature
were the Debora numbers with respect to relaxation and retardation time, Pr = υ f α represented the Prandtl number, Ω = ζλ 1 indicated the thermal relaxation time, Sc = υ f D I is the Schmidt number, K = the diffusion coefficient, When D I = D J then δ = 1 and as a result:
mass diffusivity.Physically, the Schmidt number is related to hydrodynamic layer's thickness and boundary layer.The escalating Sc intensifies the momentum of the boundary layer flow which results in an increase in concentration profile.It is clear from the figure that ()  upsurges with the rise in

Figure 9 .
Figure 9. Impact of  on

Figure 8 .Figure 8 .
Figure 8. Impact of M on C f x .

Figure 9 .
Figure 9. Impact of  on x Nu .
the conclusions associated with emerging parameters of skin fraction coefficients.This shows the impression of magenatic field parameter M on skin fraction coefficients.The magenatic field parameter shows a reduction in the skin fraction coefficient.Table 2 demonstrated the conclusion of incipient parameters on local Nusselt numbers.

Table 1 .
Distinction in −C f x for different M.

Table 2 .
Distinction in Nu x for different Ω and Pr.

Table 3 .
Distinction in Sh x for different Sc, K and K s .

Table 4 .
Symmetry of HAM versus numerical solutions for F (η), when Sc = Pr