Chemically reactive flow of upper-convected Maxwell fluid with Cattaneo–Christov heat flux model

This attempt concentrates on impact of chemically reactive flow of upper-convected Maxwell liquid. Nonlinear slip condition for Maxwell fluid is employed. The processes of heat and mass transfer through theory of Cattaneo–Christov flux are studied. Ordinary differential systems have been considered. Convergent solutions are constructed for the governing equations. Incoming nonlinear modeled problems have been computed for the velocity, temperature and concentration. The impact of emerging variables, namely Deborah number (β), Schmidt number (Sc), Thermal relaxation parameter (γ), Prandtl number (Pr) and chemical reaction parameter (δ) on quantities of interest is graphically investigated. Both temperature and concentration fields decay when thermal relaxation and chemical reaction parameters are increased.


Introduction
There is a wide range of chemical reactions in nature which have widespread practical applications. These reactions are involved in various processes, especially in fog formation and dispersion, food processing, hydrometallurgical industry, air and water pollutions, atmospheric flows, fibers insulation and crops damage due to freezing, etc. In these processes the molecular diffusion of species on the boundary or inside the chemical reaction is very intricate. Some of the reactions have the capacity to proceed gradually or do not react at the moment without catalyst. In this direction Merkin [ [12]. Impact of nonlinear thermal radiation and induced magnetic field in flow of viscoelastic material with heterogeneous and homogeneous reactions is reported by Animasaun et al. [13]. Raju et al. [14] examined induced magnetic field effects in flow of Casson liquid with heterogeneous and homogeneous reactions. Qayyum et al. [15] examined heterogeneous-homogeneous reactions flow of silver-water and copper-water nanoparticles in the presence of nonlinear thermal radiation.
It is greatly acknowledged that in circumstances comprising extremely small times, maximal temperatures or thermal gradients nearby absolute zero, heat diffusion concept provided by Fourier becomes imprecise and non-Fourier consequence becomes decisive in characterizing the diffusion mechanism and anticipating temperature distribution. Practical prospects where deviance from the Fourier's model turns noteworthy may be encountered, for example, in microelectronic materials including IC chips, heating of laser pulse with high heat flux or exceptionally short duration for hardening of semiconductors, impulse drying and laser surgery in biomedical engineering. Numerous investigations have been reported in order to present new formulation for heat conduction. For instance Straughan [16] disclosed thermal convection phenomenon utilizing non-Fourier heat conduction concept. Non-Fourier heat conduction effectiveness in stretching flow of Maxwell liquid is presented by Han et al. [17]. Waqas et al. [18] extended the idea presented in [17] considering temperature-dependent conductivity. Analysis of radiation and non-Fourier flux in differentially heated two-dimensional square cavity is developed by Sasmal and Mishra [19]. Chemical reaction and stagnation point characteristics in stratified variable thermal conductivity stretched flow of Eyring-Powell material in presence of non-Fourier heat conduction is explored by Hayat et al. [20]. Upper-convected Maxwell liquid flow with non-Fourier heat flux is studied by Saleem et al. [21]. Makinde et al. [22] presented magnetohydrodynamic flow over various geometries with Cattaneo-Christov heat flux.
The aforestated investigations witness that slip effects are not explored correctly. No doubt wall slip appears in complex liquids comprising suspensions, emulsions, foams and polymer analysis. The liquids which communicate boundary slip characteristics possess considerable demands in internal cavities and artificial heart valves [23]. MHD slip flow near a stagnation point in the presence of variable thicked surface is investigated by Khan et al. [24] and Babu and Sandeep [25]. Thermal radiation and mixed convection effects in stagnation point slipped flow of viscous liquid is explored by Rashad et al. [26]. Aly and Sayed [27] presented a comparative analysis considering four types of nanoparticles in MHD radiative stretched flow of viscous liquid. Simultaneous influences of mixed convection and magnetohydrodynamics in axisymmetric stretched flow of viscous material through slip effects and convective conditions are analyzed by Ganesh et al. [28]. Hayat et al. [29] established numerical solutions for radiative viscous nanomaterial considering melting effects.
It has been found from the existing information that mostly the flow of non-Newtonian fluid in the presence of slip condition valid for viscous fluid is considered. This is not correct. No doubt the stress in non-Newtonian material even for incompressible case is different than the viscous fluid. Thus, in reality the slip conditions for viscous and non-Newtonian fluids are distinct. Motivated in such fact our prime intention here was to report the slip effects in Maxwell material induced by stretchable surface. This model is capable of predicting relaxation time characteristics [30][31][32][33][34]. Polymer having low molecular weight is the best example for Maxwell material. Note that slip condition in viscous fluid is linear whereas it becomes nonlinear in Maxwell fluid situation. Also the generalized concept of heat and mass fluxes are imposed. Such concept has been used in view of Cattaneo-Christov theory. Another important feature is related to the consideration of stretching surface with chemical reaction. Besides this homotopy concept [35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50] is implemented for arising nonlinear problems. Velocity, temperature and concentration are described for meaningful discussion considering important variables.

Formulation
Let us consider flow of an incompressible Maxwell fluid over a stretching sheet with slip effect. Heat and mass fluxes have been characterized using Cattaneo-Christov theory. Further, we assume T w and C w as the temperature and concentration of stretched surface at y ¼ 0 while T 1 and C 1 denote the ambient temperature and concentration. The present flow is governed by the following basic expressions: dC dt ¼ Àr:j; ð4Þ Now (4) gives in which C represents the concentration of species, D * the mass diffusivity and j the mass flux. The Cauchy stress tensor for Maxwell fluid model is where the extra stress tensor S satisfies the following relation: The velocity, temperature and concentration fields are V ¼ ½uðx; yÞ; vðx; yÞ; 0; T ¼ Tðx; yÞ; C ¼ Cðx; yÞ: The boundary layer equations now are u ! 0 when y ! 1: To study the characteristics of heat and mass transfer we have For k 2 = 0 the above equation reduces to Fourier law of heat conduction. For incompressible fluid r:v ¼ 0 and then Eq. (14) gives The energy equation is Equations (15) and (16) with boundary conditions The concentration equation is Employing the transformations the mass conservation law (11) is identically fulfilled and Eqs. (12,13,17,19,20) are reduced to the forms: The dimensionless variables are Homotopic solutions and convergence analysis. We consider f(g), h(g) and u(g) via set of base functions in the forms where a k; n ; b k; n and c k; n are coefficients. Thus all the approximations of f g ð Þ; h g ð Þ and u g ð Þ must obey the above expressions. The initial approximations and operators are ð32Þ with Þare the arbitrary constants. The general solutions f m ; h m and u m ð Þin terms of where Àb Àb The convergence of series solutions in HAM procedure is quite necessary. Such convergence analysis heavily depends upon the auxiliary variables h. Thus hÀ curves for f ; h and u have been displayed in Fig. 1. Figure clearly   ð Þ will produce a retarding force between two adjacent layers in the flow. Due to this fact there is a reduction in the velocity and associated layer thickness. The results of viscous fluid are obtained when b ¼ 0: Figure 3 is sketched for the effect of ðcÞ on temperature hðgÞ: It is scrutinized that temperature hðgÞ shows recessive behavior for higher estimation of thermal relaxation parameter c ð Þ: It is found that larger values of c decreases both temperature field and thermal layer thickness. Physically for higher estimation of thermal relaxation parameter c ð Þ a material requires more time to transfer heat from more energetic particles to low energetic particles, i.e., it demonstrates the features of a non-conducting material. Therefore, temperature hðgÞ decays. Further, it is noted that for c ¼ 0 ð Þ the heat transfers without any delay through the whole material. Thus temperature hðgÞ dominants for Fourier law (i.e., for c ¼ 0) in comparison to Cattaneo-Christov heat flux model. Behavior of Prandtl number Pr ð Þ on temperature hðgÞ is examined in Fig. 4. Higher estimation of Pr decays the temperature field and thermal layer thickness. It is due the fact that Pr is the combination of thermal diffusivity to momentum diffusivity. Therefore, small values of Pr; ðPr ( 1Þ mean the thermal diffusivity dominates. For higher estimation of Pr; ðPr ) 1Þ the momentum diffusivity dominates. Brownian diffusion coefficient decreases due to which concentration boundary layer is reduced. Influence of chemical reaction parameter d ð Þ on concentration uðgÞ is displayed in Fig. 5. As anticipated, a reduction in concentration uðgÞ is observed when the chemical reaction parameter d [ 0 ð Þ is increased. Figure 6 illustrates the behavior of Schmidt number ðScÞ on concentration uðgÞ: Ratio of viscous to molecular diffusion rate is known as Schmidt number. Larger viscous diffusion rate is observed when Schmidt number is increased and consequently fluid concentration enhances. Table 1 is constructed to show the order of convergence. f 00 ð0Þ; h 0 ð0Þ and u 0 ð0Þ converge at 15th and 16th order of approximations, respectively. To examine the correctness of flow problem, we have compared our results with the published results of Makinde et al. [51] in Table 2. The results are found in an excellent agreement.

Closing remarks
In this paper we employ the upper-convected Maxwell model and non-Fourier heat flux model to investigate heat and mass transfer above a stretching plate with velocity slip. The numerical results suggest the following: • Velocity decays for higher estimation of b: • Temperature field decreases for larger Pr and c:  • Both h and u decay by increasing c and d: • Concentration of fluid decreases for higher estimation of Sc: Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.