A Galerkin strategy for tri - hybridized mixture in ethylene glycol comprising variable di ﬀ usion and thermal conductivity using non - Fourier ’ s theory

: This research is conducted to investigate heat and mass transport past over a stretched surface having pores in a pseudo - plastic model. To study porosity e ﬀ ect, Darcy Forchheimer relation is used. Thermal and mass transport expressions are derived by engaging the double di ﬀ usion theories as extensively used by researchers pro posed by Cattaneo and Christov. Furthermore, the thermal performance is studied by mixing the tri - hybrid nano particles in a pseudo - plastic material. The phenomenon of boundary layer is used to derive the complex model. The correlation for tri - hybrid nanoparticles is used to con vert the model partial di ﬀ erential equations into ordinary di ﬀ erential equations ( ODE ) along with appropriate simi larity transformation. The trans ﬁ gured ODEs are coupled nonlinear in nature, and the exact solution is not possible. To approximate the solution numerically, ﬁ nite element scheme ( FES ) is used and code is developed in MAPLE 18.0 for the graphical results, grid independent survey, and tabular results. The obtained results are compared with the published ﬁ ndings that con ﬁ rm the accuracy and authenticity of the solution and engaged scheme. From the performed analysis, it is concluded that FES can be applied to complex engineering problems. Furthermore, it is monitored that nanoparticles are essential to boost the thermal performance and higher estimation of Schmidt number control the mass di ﬀ usion.


Introduction
Fluid flows past over a stretched and rotated surfaces have gained much attention of the researchers due to their wider applications. Diffusion theories play a vital role to control heat and mass transportation in fluid flows. Researchers have much attraction to study physical problems discussing thermal transport. To keep thermal stability and maintain high thermal performance, mixing of nanoparticles is essential, and it boosts the heat transfer in fluid flows. Several models of nanofluids have been proposed and extensively used. Buongiorno's model discusses the involvement of Brownian motion and thermophoresis, whereas the single phase model contains the involvement of volume fraction and also researchers have addressed numerous relations that present the involvement of particles shapes. Several remarkable attempts have been made to discuss the utilization of nanoparticles in fluid flows. For instance, Koriko et al. [1] discussed the growth of energy transfer inserting nanoparticles along with bioconvection flow in the presence of magnetic field and gyrotactic microorganisms towards a vertical surface. Ali et al. [2] studied the role of hybrid nanoparticles in the energy transfer under the action of viscous dissipation. [3] investigated the study of various nanoparticles in wavy wall in view of heat transfer phenomena. Pop et al. [4] discussed the role of magnetohydrodynamic in stagnation point flow along with hybrid nanoparticles considering melting aspects. Tlili et al. [5] estimated features of nanoparticles in an unsteady flow using various shapes of nanoparticles considering Oldroyd-B ferrofluid under the action of magnetic field. Qin [6] estimated the heat transfer features including the role of nanostructures in chamber. Nawaz et al. [7] discussed the comparative analysis among hybrid nanoparticles and nanoparticles into dust particles considering hyperbolic tangent liquid adopting the finite element approach. Nazir et al. [8] used the finite element approach to analyse comparative investigation between hybrid nanostructures and nanostructure in base liquid called ethylene glycol in Carreau Yasuda liquid. Nazir et al. [9] used non-Fourier's law in rheology of hyperbolic tangent liquid inserting vital impact of hybrid nanostructures towards a stretching surface. Nazir et al. [10] discussed the thermal character of hybrid nanostructures in hyperbolic tangent rheology in the base fluid (ethylene glycol) via finite element approach. Sheikholeslami and Farshad [11] concluded thermal aspects of the solar collector system inserting hybrid nanostructures. Swain et al. [12] used an exponentially surface to analyse a role of hybrid nanostructures in the presence of chemical reaction via slip conditions. Faraji et al. [13] developed a model related to phase change liquid in view of hybrid nanomaterials along with heat source in a rectangular enclosure. Abdelmalek et al. [14] investigated an inclination into heat energy using nanoparticles and hybrid nanostructures in view of base fluid in wedge. They have used finite element scheme (FES) to know comparison among hybrid nanostructures and nanostructures including chemical reaction into solute particles. Nawaz et al. [15] analysed rheology of power law liquid into solute particles and thermal energy inserting an impact of hybrid nanoparticles over a heated surface via finite element approach. Latocha et al. [16] estimated prediction of hydroxyapatite in nanoparticles towards 3D-printed reactors. Haider et al. [17] discussed the role of hybrid nanoparticles in Williamson liquid using heat flux model. They adopted heated surface to captures comparison performance of hybrid nanoparticles and hybrid nanoparticles in ethylene glycol. Nazir et al. [18] discussed features of thermal energy using nanoparticles towards a surface in Casson liquid. Dadheech et al. [19] analysed the comparison of heat energy performance among SiO MoS 2 2 − hybrid nanofluid and nanofluid MoS 2 ( ) in ethylene glycol under the presence of magnetic field. Shafiq et al. [20] discussed the heat energy features in Walters' B nanofluid over a Riga plate via statistical approach. Marzougui et al. [21] used the lid-driven cavity to capture the impacts of entropy generation and thermal transfer in the attendance of non-variable magnetic field. Pushpa et al. [22] captured flow and convective thermal transfer inserting nanoparticles in thin baffle. Dhif et al. [23] studied the performance of hybrid nanoparticles in a solar collector system. Abdelsalam et al. [24] discussed the study related to Rabinowitsch suspension in the presence of leveraging elasticity in a wavelike conduit. Abdelsalam et al. [25] conducted results of swimming sperms in the presence of electro-magnetically using a wavelike conduit. Raza et al. [26] adopted curved surface to achieve consequences in Williamson liquid under thermal radiation. Eldesoky et al. [27] integrated thermal features into conjunction along with slip conditions in view of peristaltic motion using a catheterized pipe. Bhatti and Abdelsalam [28] captured peristaltic motion adding hybrid nanoparticles along with magnetic field. Elkoumy et al. [29] performed the role of Maxwell liquid in terms of peristaltic motion inserting magnetic field and Hall force. Bhatti et al. [30] simulated the heat transfer based on intra-uterine motion in channel. Some important contributions covering the modelling aspects are reported in earlier studies [31][32][33][34][35][36].

Saleem and Heidarshenas
Current model is developed using the rheology of pseudo-plastic material along with the presence of tri-hybrid nanoparticles. The theories regarding non-Fourier's and Darcy's Forchheimer within heat generation and chemical reaction are investigated. FES is considered to find numerical simulations. Investigations regarding current model is not achieved yet. • Section 1 presents literature survey; • Section 2 presents modelling; • Section 3 contains methodology procedure and meshfree investigations; and • Section 4 reports detailed discussion and description of obtained solution against numerous involved parameters, and Section 5 lists important findings. Developing approach of tri-hybrid nanoparticles is shown in Figure 1.

Flow analysis of tri-hybrid nanoparticles
The flow of tri-hybrid nanoparticles in pseudo-plastic liquid towards the heated surface is addressed. Solute particles and heat energy of particles are analysed considering the theory of Cattaneo-Christov model (CCM).
Processes of chemical reaction and heat generation are implemented in the presence of CCM. The concept of Forchheimer porous is modelled in the momentum equation. Layers associated with momentum and thermal are generated because of stretching of plate. Composite relation among nanofluid, nanoparticles, hybrid nanomaterials, and tri-hybrid nanoparticles is considered in equations (12)- (14). Physical configuration of current model is captured in Figure 2 ( Table 1). System of formulated partial differential equations [33][34][35] is as follows: Where heat generation is Q, temperature is T , thermal conductivity is k, fluid density is ρ, velocities are (U V ,̇), time relaxation is λ t , specific heat capacitance is denoted by C P , kinematic viscosity is ν, power law number is m, space coordinates are y x , , tri-hybrid nanoparticles are represented by Thnf , ambient temperature is T ∞ , permeability within porous medium is k ⁎ , inertia coefficient in term of porous medium is F D , C is concentration, ambient concentration is C ∞ , chemical reaction K a , and mass diffusion is D. No-slip theory is used to generate boundary conditions (BCs) and required conditions are as follows: Transformations of model are as follows: Variable thermal conductivity and variable mass diffusion in view of tri-hybrid nanoparticles [36] are defined as follows: Transformations deliver system of ODEs and nonlinear ODEs are as follows: Properties associated with thermal for correlations among tri-hybrid nanoparticles are as follows: The formulated skin friction coefficient is defined as follows: The temperature gradient is modelled as follows: The rate of mass diffusion is as follows: Reynolds number is Re .
FES is utilized to simulate numerical results. FES is very capable of simulating CFD problems. Detail steps are described below.
Step I: Equations (6) and (7) within BCs are called the strong form. It is noticed that collecting all terms of equations (6) and (7) on one side and integrating it over each elements of domain are residuals. Such procedure is known as weighted (residual method) for the development of weak forms. The residuals are as follows: Here, Q Q Q , , 1 2 3 , and Q 4 are named as weight functions. The shape functions are as follows: Step II: An approach associated with Galerkin finite element is imposed to obtain weak form in view of shape functions.
Step III: The assembly approach is utilized for the development of stiffness element, whereas assembly approach is performed via assembly procedure of FEA. Stiffness elements are as follows: Step IV: Picard linearization approach provides the transformed algebraic system (linear equations).
Step V: Finally, system of linear algebraic equations is numerically solved within the computational tolerance (10 5 − ). The stopping condition is listed below.
Step VI: Table 2 demonstrates study of mesh-free.
Step VII: A total of 300 elements are required to obtain convergence analysis. The programming of finite element method is generated using software related to MAPLE 18.

Validation of results
The comparison among present analysis and published results [33][34][35] are derived in Table 3. It is noticed that the present problem is reduced into published problems [33][34][35]

Analysis of velocity curves against physical parameters
Motion into nanoparticles, hybrid nanoparticles, and trihybrid nanoparticles are measured with respect to variation in power law number (m), Darcy's number (ε), and Forchheimer number (F r ). Figure 3a-c is prepared for the measurement of the velocity distribution versus m ε , , and F r . Figure 3a prescribes the motion into nanoparticles versus distribution in ε. It is noticed that the presentation of ε happened because of pores at the surface. The flow of nanoparticles, ternary hybrid nanoparticles, and hybrid nanoparticles is declined over the surface. Momentum layers are decreased when ε is enhanced. Physically, higher values of ε produce frictional force into fluid particles. Hence, the frictional force makes a reduction into motion among fluid layers. Hence, inverse proportional relation is existed among flow and ε. Impact of F r on velocity curves is measured by Figure 3b, inserting the role of ternary hybrid nanoparticles. Appearance of F r is model out due to impact of Darcy's Forchheimer law. A declination role into motion of particles is investigated when F r is enhanced. Layers associated with the momentum boundary are investigated decreasing function against impact of F r model related to Darcy-Forchheimer is based on Darcian motion into fluid particles. Mathematically, it appeared as a velocity squared in motion equations. A retardation motion is developed due to appearance of Forchheimer. Momentum layers have a decreasing function versus impact of F r . Hence, fluid is termed as thick when F r is increased. An effect of m on velocity curves is established in Figure 3c. It is estimated that the description of m is modelled out due to inserting impact of pseudo-plastic liquid. Moreover, the category of shear thinning, shear thickening, and Newtonian liquid is based on values m. For m 1, < fluid appeared as shear thickening among fluid particles. It is noticed that layers along with momentum layers are declined when m is increased.

Analysis of temperature curves against physical parameters
Thermal aspects including ternary hybrid nanoparticles are observed against variation in H ε , h 1 , and Ω a . These thermal aspects are measured with respect to change in H ε , h 1 , and Ω a by Figure 4a-c. The impact of variable thermal conductivity (ε 1 ) on temperature curves is examined by Figure 4a. It is noticed that thermal energy of particles is boosted when ε 1 is increased. The parameter related to ε 1 is produced because of variable thermal conductivity. Mathematically, ε 1 has a relationship related to direct proportional temperature. Hence, an increment into thermal energy is based on distribution in ε 1 . Moreover, layers of thermal at boundary are increased versus higher values of ε 1 . Physically, ε 1 has a directly proportional relation versus temperature difference. Hence, an increment in ε 1 results in a significant ability of thermal energy to conduct more heat energy into fluid particles is observed. The ability to conduct temperature into hybrid and nanoparticles is also increased when ε 1 is increased. Figure 4b demonstrates influence of H h on temperature curves. Dual role of heat energy is visualized on aspects of thermal energy while negative values of H h are due to the role of heat absorption and positive values of H h are due to the role of heat energy generation. Production into heat energy of particles is inclined into nanoparticles and tri-hybrid nanoparticles based on variation in H h . Mathematically, H h appeared in energy equation, which is a product of temperature difference. Basically, it occurred because of an external heat source. Thermal energy is boosted when an external heat source is placed at the wall. Boundary layers based on the thermal impact are increasing behaviour. This increasing effect of thermal energy is produced due to applying an external source of heat energy, which is placed at the wall of surface. Hence, the heat energy is established according to values of H h . The prediction of heat energy against distribution in Ω a is addressed by Figure 4c. Argumentation into heat energy is measured when Ω a is increased. The parameter is called Ω a is made because of non-Fourier's theory. Non-Fourier's theory is implemented in current analysis. Hence, heat energy is boosted when Ω a is increased. An increasing function is investigated among Ω a and thermal layers. Moreover, fluid particles take more time to restore more heat energy when Ω a is increased. The layers associated within thermal are increasing inserting higher values of Ω a . Figure 5 is observed as a vital impact on comparison of tri-hybrid nanoparticles, nanoparticles, and fluid and hybrid nanoparticles on the thermal performance. It is examined that a significant performance of ternary hybrid nanoparticles is addressed on the thermal performance. Maximum production of thermal energy is produced for tri-hybrid nanoparticles rather than production into thermal energy of hybrid nanoparticles, fluid and hybrid nanoparticles.

Analysis of concentration curves against physical parameters
Diffusion of solute particles is observed versus impacts of K , Ω , Sc c c , and ε 2 considering Figure 6a-d. Figure 6a illustrates that the variation into thermal energy is verified against chemical reaction number (K c ). In this figure,  dual trends of K c is measured by Figure 6a. Two kinds of trend are called destructive and generative chemical reactions for K 0 c > and K 0 c < , respectively. It is estimated that the chemical reaction is known as non-reactive reaction for K 0. c = Diffusion into species is decreased versus large values of K c . Further, the generative chemical reaction is higher than the reaction related to destructive and non-reactive reactions. Diffusion into species slows down because chemical reaction among fluid particles occurred. Hence, layers associated with concentration are also reduced when chemical reaction becomes fast among nanoparticles. An influence of Schmidt number (Sc) on concentration curves is addressed by Figure 6b. Diffusion of species slows down with respect to change in Sc. Physically, Sc has inverse proportional versus mass diffusivity. Therefore, an inclination into mass diffusion is addressed when mass diffusivity is increased. Hence, a declination role is predicted versus impact of Sc. Moreover, concentration layers are inclined with respect to change in Sc. It is the ratio of viscous diffusion and mass diffusion rates. Therefore, mass diffusion rate is decreased versus higher impact of Sc. Concentration layers are declined against higher value of Sc. Appearance of Ω c is noticed on concentration curves plotted by Figure 6c. Three kinds of nanoparticles are inserted into diffusion of species, whereas higher values of Ω c are observed as a significant performance of species rate. Significant performance of mass species is observed versus the role of Ω c . Figure 6d addresses characterizations of solute particles versus the impact of Ω c Existence of Ω c is formulated due to the impact of Cattaneo--Christov theory. A significant achievement into solute particles is observed when Ω c is enhanced including tri-hybrid nanoparticles.

Analysis of Sherwood number, temperature gradient, and drag force coefficient
Temperature gradient, Sherwood number, and drag force coefficient are simulated by incorporating influences of Forchheimer number, heat generation number, chemical reaction number, and Sc including the study of ternary hybrid nanoparticles. Impacts related to these parameters on temperature gradient, Sherwood number, and drag force coefficient are simulated as shown in Table 4. Role of F r boosts surface force but the concentration and temperature gradients are declined when F r is increased. H h diminishes the rate of thermal energy, mass diffusion, and surface force. It is estimated that the dual trends (heat absorption and heat generation) are addressed. In both cases for H h are observed significantly for establish maximum production in rate of concentration and temperature gradients. A declination is predicted in view of surface force versus an impact of chemical reaction (generative and destructive) but in terms of Sherwood and Nusselt, numbers are inclined. For the case of Sc, rates of mass diffusion and thermal energy are significantly increased.

Conclusion
Solute and thermal characterizations in pseudo-plastic liquid are investigated along with the ternary hybrid nanoparticles. An expanding surface is considered to measure the flow, thermal energy, and solute particles into particles against the distribution of various physical parameters. Darcy's Forchheimer law is incorporated along with non-Fourier's law in the presence of heat generation and heat absorption. A term related to chemical reaction is also considered into solute particles. The finite element approach has used to characterize graphical outcomes. Main findings are addressed below.
• Flow slows down versus higher values of Darcy's number, power law number, and Forchheimer number; • Argumentation is investigated into heat energy of particles for the case of ternary hybrid nanoparticles rather than for case of fluid, nanoparticles, and hybrid nanoparticles. • Maximum achievement of heat energy for tri-hybrid nanoparticles as compared than nanofluid and hybrid nanoparticles. • Diffusion of species slows down when Sc is enhanced.
• Applications of ternary hybrid nanoparticles based on present problem are utilized in medicines, cancer cells, Funding information: The authors state no funding involved.
Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

Conflict of interest:
The authors state no conflict of interest.
Data availability statement: All data generated or analysed during this study are included in this published article.