Inclusion of hybrid nanoparticles in hyperbolic tangent material to explore thermal transportation via finite element approach engaging Cattaneo-Christov heat flux

This report is prepared to examine the heat transport in stagnation point mixed convective hyperbolic tangent material flow past over a linear heated stretching sheet in the presence of magnetic dipole. Phenomenon of thermal transmission plays a vital role in several industrial manufacturing processes. Heat generation is along with thermal relaxation due to Cattaneo-Christov flux is engaged while modeling the energy equation. In order to improve the thermal performance, inclusion of hybrid nanoparticles is mixed in hyperbolic tangent liquid. The conservation laws are modeled in Cartesian coordinate system and simplified via boundary layer approximation. The modeled partial differential equations (PDEs) system are converted into ordinary differential equations (ODEs) system by engaging the scaling group transformation. The converted system of modeled equations has been tackled via finite element procedure (FEP). The efficiency of used scheme has been presented by establishing the grid independent survey. Moreover, accurateness of results is shown with the help of comparative study. It is worth mentioning that the inclusion of hybrid nanoparticles has significant higher impact on heat conduction as compared with nanoparticle. Moreover, hybrid nanoparticles are more efficient to conduct maximum production of heat energy as compared with the production of heat energy of nanoparticles. Hence, hybrid nanoparticles (MoS2/Ag) are observed more significant to conduct more heat energy rather than nanoparticle (Ag).


Introduction
Non-Newtonian fluids have diverse applications and usage in several engineering disciplines. There constitutive relation is different than Newtonian materials. For non-Newtonian materials, the stress strain relationship is complicated. These materials occur frequently in different phenomenon and it has applications in polymers extrusion, medical, oil-pipeline friction reduction and several others. An important non-Newtonian model is hyperbolic tangent [1][2][3][4][5], whose constitutive expression is expressed as for m 1 ¼ 0; G 1 _ a < 1; we get Naseer et al. [1] examined the exponentially stretched flow of mixed convective hyperbolic tangent model in horizontal cylinder with thermal transportation. They engaged shooting approach to solve transformed nonlinear boundary layer equations (BLEs). They have discussed the contribution of several influential variables on velocity and temperature. They observed the depreciation in temperature field for higher values of Reynolds number and Prandtl number. Khan et al. [2] used BVPh2.0 package to handle the solution for chemically reactive hyperbolic tangent model undergoing radiation and heat generation phenomenon. They have found that velocity is a diminishing function of Weissenberg number and temperature profile grows against radiation parameter. Nawaz et al. [3] explored the contribution of hybrid nanoparticles in hyperbolic tangent model in a stretched cylinder using finite element approach. They modeled the thermal phenomenon under radiation effect and heat generation. They observed that heat transfer rate of hybrid nanoparticles is higher than nanofluid and it is recommended for better thermal performance. Mathematical analysis for inclined MHD hyperbolic tangent model was studied by Ali et al. [4] via BVP4C package. They recorded the increase in Nusselt number for Prandtl number. Patil and Raju. [5] used the convective conditions in hyperbolic tangent model past over a nonlinear exponential porous stretching sheet. They modeled the energy transport in the presence of viscous dissipation and radiation effect. They used shooting approach to compute solution for the arising boundary layer equations. They concluded that higher values of Eckert number reduce the skin friction and enhances the Sherwood number.
Nanofluids are a hot topic of research due to their wider use in different industrial mechanisms. It is used in several engineering systems to control/enhance the thermal performance. The inclusion of nanoparticles in base fluids enhances the thermal conductivity which enhance in thermal process. Several researchers have a look on this hot topic due to their wider applications and usage. For instance, Aman et al. [6] studied the thermal performance of unsteady convective Casson model by mixing SWCNTs and MWCNTs in a vertical channel. They used the thermal conductivity model which was proposed by Maxwell. They have considered radiation effect in thermal profile and mixed convection along with MHD effect is taken in momentum transport. They presented the perturbation solution for the arising problem. Several graphs and tables are presented by them to notice the influence of emerging parameters on obtained solution. They reported the increase in thermal field by enhancing the values of radiation parameter and volume fraction because it upsurges the conduction. Moreover, they recorded the enhancement in velocity against Richardson parameter. Mixed convective stratified unsteady electrically conducting viscous nanofluid was investigated by Daniel et al. [7]. They modeled a strong problem over a melting stretching sheet with several important physical effects. They modeled a strong problem by taking several important effects like mixed convection, viscous dissipation, chemical reaction, radiation and heat generation which are used in the modeled laws. They presented the numerical solution via finite difference technique and established a comparative analysis to judge the scheme efficiency and results validity. They have shown an increase in momentum and thermal profile against electric parameter. Besthapu et al. [8] worked on thermally stratified dissipated stretched nanofluid past over a heating stretching porous sheet numerically. They plotted the solution against different physical parameters and noticed the decline in fluid velocity against magnetic parameter and increase in concentration and thermal fields. Also, they reported the rise in concentration and temperature profiles against mounting values of thermophoresis parameter. Ramzan et al. [9] discussed the optimal solution foe Maxwell nanofluid model with thermal and mass transportation having soret and dufour effects past over a porous heated sheet. They presented the error analysis graphically against different approximations order. They noticed the diminishing impact of porous parameter on velocity. Also, an enhancement in heat and mass transfer coefficients is monitored by them against suction parameter and viscoelastic parameter. Modelling of water based nanoparticles in rotating disc having Darcy Forchhiemer medium with slip effects was analyzed by Hayat et al. [10] analytically. They have mixed different type of nanoparticles in water to investigate thermal performance. They used the relation proposed by Xue for thermos-physical properties. They monitored the enhancement in heat transfer coefficient against Forchhiemer number and retardation in fluid velocity for higher volume fraction. Khanafer and Vafai [11] deliberated the experimental analysis for thermophysical features of nanofluids. They have critically studied different effective viscosity models and thermal expansion coefficient in their reported study. Kempannagari et al. [12] discussed characteristics of heat transfer in micropolar liquid including the effects of thermal conductivity (variable), thermal radiation and Joule heating past a heated curved sirafce. Kumar et al. [13] discussed the aspects of thermal energy in hybrid ferrofluid considering the impacts of thermal radiation and heat generation. Kumar et al. [14] addressed the phenomenal contrast in view of heat transfer and flow behavior over wedge and heated cone. They considered the theory of non-Fourier's under the action of magnetic field via shooting approach. They also investigated that flow phenomena and heat transfer for the case of cone are higher than flow phenomena and heat transfer for the case of wedge. Kumar et al. [15] studied influences of thermal energy along with thermal radiation in the attendance of magnetic field past a heated curved surface via shooting scheme. Kumar et al. [16] captured the behavior of micropolar fluid including the BCs (2 nd order and first order slips) whereas they considered variable heat (sink) and Lorentz force. Kumar et al. [17] investigated the flow situation in the attendance of magnetic field. They used slandering heated surface to know the behavior considering effects along with viscous dissipation. Anantha Kumar et al. [18] observed rheology of micropolar fluid along with stagnation point and magnetic field. They used second-order (velocity slip) of considered flow model past a heated vertical surface. They also considered the thermal conductivity (variable), Joule heating and heat source sink and electrical energy. Ashwinkumar [19] visualized the influences of thermal energy and solute particles along with action of magnetic field inserting nanoparticles (alloy/silver-water) past an elongated heated surface. Ashwinkumar and Sulochana et al. [20] simulated consequences buoyancy force and thermal radiation under the action of magnetic field inserting Al 50 Cu 50 -nanoparticles in the rheology of Casson fluid. Ashwinkumar et al. [21] securitized features of ferrous nanoparticles under the action of magnetic field over melting vertical plate (semi-infinite). Sulochana et al. [22] discussed the aspects of heat transfer phenomena in Maxwell liquid including the nanoparticles over a melting elongated surface. They considered various kinds of influences (heat generation, thermophoresis, Brownian moment and magnetic field) simulated the results via shooting scheme. Mabood et al. [23] discussed an enhancement in thermal energy inserting the hybrid nanoparticles (Fe3O4/graphene-H 2 0) in the attendance of magnetic field. Mabood et al. [24] simulated the maximum production of heat energy due the role hybrid nanoparticles under the action of thermal radiation and magnetic field past a slandering melting surface. They used shooting scheme to simulate the simulations of considered model. Samrat et al. [25] investigated the rheology of Casson liquid inserting nanoparticles in the presence of thermal rediation past a stretching surface. Sulochana et al. [26] discussed the phenomena of thermal energy in the rheology of ferrofluid involvement of nanoparticles towards the horizontal needle via magnetic field. They studied the role of ferrous nanoparticles in solute particles and thermal energy including magnetic field and thermal radiation subjected to chemical reaction over an elongated plate. Sulochana et al. [27] studied the impact of 2D flow, thermal energy and mass diffusion inserting nanoparticles under action of thermal radiation and chemical reaction over a heated elongated plate via shooting scheme. Some important contributions are reported in [28][29][30][31][32][33][34][35][36].
Above cited studies ensure that no analysis has been done for hyperbolic tangent model with thermal relaxation time and radiation phenomenon. This contribution fills this gap. This contribution is organized as: section 1 consists of literature survey, modeling of considered flow situation with physical quantities are listed in section 2, solution methodology with advantages is reported in section 3, graphical and tabular results are analyzed in section 4 and section 5 contains the important findings. The sketched scheme of hybrid nanoparticles and nanoparticles are illustrated by Fig 1.

Modeling development
The modeling of 2D hyperbolic tangent liquid [1][2][3][4][5] inserting the thermal properties of nano and hybrid nanostructures is addressed past a melting surface. The non-Fourier's law [14,36] in the presence of heat generation is studied in energy equation. The movement in fluid particles is responsible due to velocity (bx) of surface. The role of magnetic induction is captured and magnetic field is also taken at away from the surface considered by H 0 x. The base fluid is considered as a EG (ethylene glycol) along with hyperbolic tangent liquid while nanoparticles (Ag) and hybrid nanoparticles are called MoS 2 /Ag. The geometrical phenomenon is demonstrated by Fig 2. The considerations of current model are addressed below: ➢ Two dimensional flow in the presence of hyperbolic tangent liquid is considered; ➢ Cattaneo-Christove heat flux is addressed; WhereŨ ;Ṽ (velocity components), x,y (space coordinates), m (power law index number), H (induced magnetic number), H a , H b (components of magnetic induction),m (magnetic permeability), μ e (magnetic diffusivity), T (fluid temperature), ρ (fluid density), T 1 (ambient temperature), Q (heat generation), k (thermal conductivity), λ A (relaxation number), G (gravitation force), μ f (fluid dynamic viscosity), Γ 1 (time constant number) and β (thermal expansion) are mentioned above.
No-slip theory provides required boundary conditions (BCs) arẽ The variables are used to obtain ODEs and these variables are: The system of ODEs along with desired BCs are addressed [35] as The correlations of thermo-physical properties in nano and hybrid nanoparticles are mentioned below and their values are listed in Table 1.

Solution scheme and convergence analysis
The non-linear ODEs are numerically simulated using numerical scheme called finite element approach [33,34]. Such numerical approach has ability to compute numerical solution of current complex problem. The main steps related this scheme is discussed here: ❖ The residuals (weighted integral) are made whereas the residuals are made using the integrating approach in set of non-linear ODEs (ordinary differential equations); ❖ In next step, computational domain is achieved through problem domain considering [0, 1). It is demonstrated that numerical computations reveal that asymptotic BCs are agreed at η 1 = 7. Hence, [0, 7] is considered as computation domain; ❖ The stiffness matrix is computed using GFEA (Galerikin finite element approximation) in weak form of residuals (weighted integrals). The stiffness elements along with 300 elements in integral residuals are made. Finally, assembly scheme is captured; ❖ A Picard linearization scheme is utilized to make linearization in algebraic equations. The tolerance is captured as 10 −8 for solving equations; ❖ The convergence approach is addressed considering 300 elements. Hence, current problem is converged at mid of each 300 elements. This convergence study is simulated in Table 2. Each outcome is simulated at the mid of each element. The validation of numerical results related to gradient temperature is simulated by Table 3.

Results and discussion
The code related to finite element approach called Galerikin finite element scheme is developed in Maple 18. This numerical approach has ability to simulate the complex flow problems. The current problem is addressed physical happening in view of PDEs including characterizations magnetic induction and non-Fourier's theory. In this section, prime graphical outcomes of motion and thermal energy in fluid particles in the presence of heat enhancement role of nano and hybrid nanoparticles versus variation parameters are simulated. The prime discussion is addressed below:

Graphical outcomes of flow situation
The measurement of flow situation is addressed versus the variation in We (Weissenberg number), H s (heat generation number), λ 1 (mixed convection number) and m (power law index number) inserting the nanoparticles and hybrid-nanostructures is captured. It is investigated    Meanwhile, bouncy force has a direct relation versus the pressure gradient which generates the maximum magnitude in MBL (momentum boundary layer thickness). The values of λ 1 >0 and λ 1 <0 are associated along with the role of cooling and heating processes in fluid particles. In graphical view, cooling process is considered for decrement in flow and heating process is assumed for increasing flow. Moreover, nanoparticles are not enough significant to obtain maximum flow as compared maximum flow for hybrid nanoparticles. Hence, the change in λ 1 becomes favorable physical parameter against temperature gradient and bouncy force. This happening derives the fluid significant over heated surface. Here, dual role of λ 1 is notices against flow. The thickness related (BL) is increased using the higher values of λ 1 . The variation in H s against the primary velocity is conducted by Fig 6. The flow is enhanced using enlargement in H s . The internal heat in fluid particles are increased due to H s and internal heat has direct relation versus kinetic energy. Hence, direct relation is simulated versus H s . Therefore, MBL become significant. The large values of H s produces more heat energy into fluid particles. Therefore, heat energy of hybrid nanoparticles and nanoparticles is increased. It is noticed that property related to generative of heat energy is not suitable when hybrid nanoparticles are considered as a coolant.

Graphical outcomes of thermal energy
The measurement of temperature is essential part in this research article. In fact, comparative simulations are captured inserting nano and hybrid nano-structures in base fluid versus H s , O a and λ 1 . In Figs 8-10, it is essential to mention that hybrid nanoparticles are more efficient to obtain maximum enhancement in thermal conductivity rather than the role of nanoparticles. Therefore, hybrid nanoparticles are efficient strong to achieve maximum production of heat energy and cooling process. Fig 8 conducts the distribution of thermal energy against λ 1 . The reduction in the production of heat energy is achieved versus the higher λ 1 . In physical point of view, temperature gradient enhances due to large values of λ 1 and more acceleration is produced in velocity. Physically, flow of fluid particles accelerates and temperature gradient increases due to higher values of λ 1 . These dual influences enhance the flow and thermal energy transfer. Further, These dual influences decrees internal energy of particles. Hence, these happenings create decline in temperature profile. The solid curves are higher than curved related to das dot lines. The graphical study with respect to H s is addressed by Fig 9. The fluid particles are heated due internal energy source. So, more heat is generated in fluid particles using H s . The large values of H s produces more heat energy into fluid particles. Therefore, heat energy of hybrid nanoparticles and nanoparticles is increased. It is noticed that external heat source is considered at surface of sheet generate more heat into fluid particles. Hence, fluid particles are heated due to enhancement of heat generation number. Therefore, increasing function is conducted of heat energy. Moreover, maximum amount of heat energy is captured inserting hybrid nanoparticles rather than inserting the nanoparticles.

Outcomes related to gradient temperature and divergent flow
The divergent flow and gradient temperature are noticed versus We, λ 1 , H s and O a simulating in Table 4. The reduction in gradient called skin friction coefficient versus variation in We and λ 1 but increasing values of divergent flow is simulated against We. Physically, flow of fluid particles accelerates and temperature gradient increases due to higher values of λ 1 . These dual influences enhance the flow and thermal energy transfer. Further, These dual influences decrees internal energy of particles. In comparative view, hybrid nanoparticles make the more  enhancements in divergent flow as compared nanoparticles. Hence, nanoparticles are more essential rather than study of nanoparticles in view of divergent velocity. The gradient temperature is enhanced versus heat generation number due to gravitational force. But inverse trend is simulated in temperature gradient due to using the large values of O a . Meanwhile, it is observed that study of hybrid nanoparticles make the more impact to attain maximum temperature gradient called Nusselt number rather than the study of nanoparticles.

Conclusion and prime findings
The features of magnetic induction in tangent hyperbolic liquid suspending with nanoparticles and hybrid nano-structures is addressed. The role of non-Fourier's theory along with heat generation is considered. This complex model is simulated with strong technique called finite element approach (FEM). The prime outcomes related to current model is captured below: • The hybrid nanoparticles are more efficient to conduct maximum production of heat energy as compared the production of heat energy of nanoparticles. Hence, hybrid nanoparticles (MoS 2 /Ag) are observed more significant to conduct more heat energy rather than nanoparticles (Ag); • The flow is significantly enhanced for the case of hybrid nanomaterials as compared for the case nanomaterials. The flow of hybrid nanoparticles is higher than flow of nanoparticles considering various physical parameters; • The convergence of current model is analyzed up to 300 elements. The implementation of FEM is considered very useful for discretization of derivatives and this approach handles various types BCs (boundary conditions); • The power law and Weissenberg numbers indicate the reduction in growth of velocity but flow is enhanced versus the large values of heat generation and mixed convection number. The thickness of boundary layers is declined considering higher values of Weissenberg number while thickness of boundary layers in inclined versus the change in power law number; • The maximum amount of heat energy is simulated using large value heat generation number but less production of heat energy is attained against the variation of time relaxation and mixed convection numbers. Further, thickness related to thermal layers becomes higher versus the values of heat generation number. The variation of time relaxation and mixed convection numbers make the reduction in thermal boundary layer thickness.
• The temperature gradient is enhanced versus enlargement in time relaxation and mixed convection numbers while decimation in temperature gradient is captured using higher values of heat generation and Weissenberg numbers. Moreover, hybrid nanoparticles are also considered useful to maximized production of temperature gradient rather than the production of temperature gradient for the case of nanoparticles; • The skin friction is reduced versus enhancement in heat generation mixed convection numbers but reverse role is addressed versus the higher values of Weissenberg and mixed convection numbers. The surface force becomes higher for hybrid nanoparticles than surface force for the case of nanoparticles.