Computational Assessment of Thermal and Solute Mechanisms in Carreau–Yasuda Hybrid Nanoparticles Involving Soret and Dufour Effects over Porous Surface

Engineers, scientists and mathematicians are greatly concerned about the thermal stability/instability of any physical system. Current contemplation discusses the role of the Soret and Dufour effects in hydro-magnetized Carreau–Yasuda liquid passed over a permeable stretched surface. Several important effects were considered while modelling the thermal transport, including Joule heating, viscous dissipation, and heat generation/absorption. Mass transportation is presented in the presence of a chemical reaction. Different nanoparticle types were mixed in the Carreau–Yasuda liquid in order to study thermal performance. Initially, governing laws were modelled in the form of PDEs. Suitable transformation was engaged for conversion into ODEs and then the resulting ODEs were handled via FEM (Finite Element Method). Grid independent analysis was performed to determine the effectiveness of the chosen methodology. Several important physical effects were explored by augmenting the values of the influential parameters. Heat and mass transfer rates were computed against different parameters and discussed in detail.


Introduction
The mechanism of transport phenomenon in different materials has received reasonable attention recently due to its wider applications in industry and different medical processes. Several important materials exist for the support of these mechanisms. Due to their different characteristics, these materials cannot be explained through one constitutive relation. Carreau-Yasuda is one such important material which has the following constitute relation.
For Y = 0 or n = 1, the Newtonian model is recovered. This model predicts the relation of shear stress with frequency. Several important contributions have been made by considering this material. For example, Zare et al. [1] discussed this model by experimentally considering the complex viscosity relationship. In their investigation, they found an excellent settlement of frequency data. They considered the involvement of carbon nanotubes in the mixture of Carreau-Yasuda material. Kayani et el. [2] reported on the behavior of wall properties on the peristaltic flow of the Carreau-Yasuda model in a sinusoidal channel by considering the Hall effect. Governing laws for the transport of species, heat and momentum were modeled under the low-Reynolds-number assumption along with the long wavelength approach. After implementing the scaling group transformation, the transformed problem was approximated numerically via the ND-Solve tool in the MATHEMATICA 15.0 symbolic package. The authors conducted a parametric analysis and their findings were shown in several graphs. They noticed a decline in the thermal field for the Biot number and recorded an enhancement in the mounting values of the Brinkman number. The unsteady rheology of the Carreau-Yasuda material in a circular tube was examined by Rana and Murthy [3]. In their investigation, they reported the wall absorption effect. They retrieved different flow behavior cases by considering different values of material parameters. Sochi [4] presented the modelling of Cross and Carreau liquid through a circular pipe. Analytical and numerical schemes were jointly implemented for the solution of flow equations and an excellent settlement was monitored. The rheology of the Carreau-Yasuda model in a cavity at high-Reynolds-number was examined by Shamekhi and Aliabadi [5] via the mesh-free algorithm. The phenomenon of blood flow via the Carreau-Yasuda model was reported by Jahangiri [6] through the FEM package. The involvement of nanoparticles enhanced the thermal performance and heat transportation rate. Several models of the nanomaterials are available and frequently used to study the thermal performance of different materials. Several researchers have paid attention to these materials due to their wider applications and usage. For instance, Gorla and Gireesha [7] developed the modeling of steady viscoelastic material with convective heat transport. The Buongiorno model is utilized to capture the characteristics of Brownian diffusion and thermophoresis. Modelling of heat transport is carried out by considering thermal radiation and heat generation. They solved boundary layer equations via a shooting procedure in the MATLAB symbolic package. The impact of several pertinent parameters were displayed through graphs, and tabular results were prepared to demonstrate the effectiveness and applicability of the shooting scheme for a large set of nonlinear data arising in the mechanical engineering problem. Muhammad et al. [8] modelled the squeezed flow of a viscous nanofluid with an updated mass and heat fluxes between parallel plates and handled the resulting expressions analytically via the OHAM in MATHEMATICA 15.0 computational tool. They noticed the enhancement in fluid velocity for the larger squeezing parameter. Rashid et al. [9] presented the exact solution of a water-based mixture of aligned nanoparticle materials under the radiation effect. They recorded the depreciation in heat transfer rate against the radiation and slip parameter. The double stratification phenomenon in the buoyancy-driven flow of a micro-polar viscous nanofluid was examined by Ramzan et al. [10]. They analyzed the depreciation in velocity for the buoyancy parameter and an enhancement was recorded against the micro-polar parameter. Hezma et al. [11] studied the behavior of SWCNTs in order to investigate the mechanical properties of polyvinyl chloride. Upadhay and Raju [12] examined the inclusion of dust particles in Eyring-Powell material over a stretched sheet. They studied the thermal and mass transport in dusty Eyring-Powell nanofluid by engaging the revised definitions of mass and heat fluxes. They found the numerical solutions for nonlinear modeled equations via the shooting scheme. Several important pieces of research on the transport phenomenon are reported in [13][14][15][16][17] and the references therein. Nazir et al. [18] studied comparison analysis among hybrid nanoparticles and nanoparticles in hyperbolic tangent liquid past a stretching surface. They adopted a finite element approach to conduct numerical results. Chu et al. [19] modeled correlations between nanoparticles and hybrid nanoparticles, considering activation energy and chemical reaction. They noticed thermal aspects past a parabolic surface using a finite element scheme. Cui et al. [20] simulated effects related to radius and roughness of inserting nanoparticles. Awais et al. [21] applied the KKL model in the transfer of energy using nanoparticles. Nazir et al. [22] discussed the numerical results of the Carreau-Yasuda liquid in heat/energy considering the hybrid nanoparticles and nanomaterials, numerically solved by the finite element approach.
In the above cited literature, no study deals with the combined behavior of the following: mass, heat transport in hydro-magnetized Carreau-Yasuda material using Joule heating, viscous dissipation, heat generation, chemical reaction and the Soret and Dufour influences in the Darcy-Forchheimer porous stretching sheet. This report fills the gap in this discussion and should be used as a foundation for researchers to work further on this model. The inclusion of nanoparticles in the Carreau-Yasuda material is attractive to researchers. Organization of this research is divided in the following way: the literature survey is reported in Section 1, modelling with important physical assumptions are covered in Section 2, Section 3 covers a detailed description of the finite element procedure with a grid independent survey, a detailed description of the solution and the influence of several emerging parameters are explained in Section 4 and important findings of the reported study are listed in Section 5. Figure 1  applied the KKL model in the transfer of energy using nanoparticles. Nazir et al. [22] discussed the numerical results of the Carreau-Yasuda liquid in heat/energy considering the hybrid nanoparticles and nanomaterials, numerically solved by the finite element approach.
In the above cited literature, no study deals with the combined behavior of the following: mass, heat transport in hydro-magnetized Carreau-Yasuda material using Joule heating, viscous dissipation, heat generation, chemical reaction and the Soret and Dufour influences in the Darcy-Forchheimer porous stretching sheet. This report fills the gap in this discussion and should be used as a foundation for researchers to work further on this model. The inclusion of nanoparticles in the Carreau-Yasuda material is attractive to researchers. Organization of this research is divided in the following way: the literature survey is reported in Section 1, modelling with important physical assumptions are covered in Section 2, Section 3 covers a detailed description of the finite element procedure with a grid independent survey, a detailed description of the solution and the influence of several emerging parameters are explained in Section 4 and important findings of the reported study are listed in Section 5. Figure 1 reveals the division of the base fluids, hybrid nanoparticles and nanoparticles. In this Figure, , , , and are known as nanoparticles whereas , ethylene glycol and oils are called base fluids. In this current analysis, ethylene glycol is considered as a base fluid. Mixtures of and are hybrid nanoparticles.

Development of the Flow Model
An enhancement in the thermal and solute performance of Carreau-Yasuda rheology, inserting the impact of nanoparticles and hybrid nanoparticles, is considered as shown in Figure 1. The flow runs towards the stretching surface under the action of a constant magnetic field. Heat takes place due to Joule heating and viscous dissipation. The Soret and Dufour influences are captured with heat generation and chemical reaction. Forchheimer's porous theory is imposed in the transport phenomenon. The geometrical flow diagram is considered in Figure 2 and the thermal properties of the nanoparticles are shown in Figure 3.

Development of the Flow Model
An enhancement in the thermal and solute performance of Carreau-Yasuda rheology, inserting the impact of nanoparticles and hybrid nanoparticles, is considered as shown in Figure 1. The flow runs towards the stretching surface under the action of a constant magnetic field. Heat takes place due to Joule heating and viscous dissipation. The Soret and Dufour influences are captured with heat generation and chemical reaction. Forchheimer's porous theory is imposed in the transport phenomenon. The geometrical flow diagram is considered in Figure 2 and the thermal properties of the nanoparticles are shown in Figure 3.
The non-linear PDEs are developed according to the physical happenings and the boundary layer approximations.
∂ u ∂x where ( u, v) is velocity componets, space coordinates are (x, y), kinamtic viscosity is ν, F s is inertia cofficient (porous medium), permeability (porou medium) is k a , power law index number is m, time constant is Λ, magnetic induction is B 0 , electrical conductivity is σ, temprature is T, C is the concentration, heat source is Q, T ∞ is ambient temprature, fluid density is ρ, C is concentarion, mass diffusion is D, T m is fluid mean temprature, k T is thermal diffusion, C s is concentration susceptibility and k 0 is chemical reaction number. The no-slip theory provides the required boundary conditions of the current model: Change in variables is constructed as: Transformations are used in Equations (1)-(5) and the system of non-linear PDEs are converted to following ODEs Here, A 1 , A 2 , A 3 , A 4 , A 5 and H 1 are involved parameters (representing the correlation of nanoparticles and hybrid nanostructures) in the above equations which are defined as Figure 3 demonstrates the thermal properties of density, electrical conductivity, thermal conductivity and specific heat capacitance for ethylene glycol, MoS 2 /SiO 2 and MoS 2 . Density of C 2 H 6 O 2 is 1113.5, density of MoS 2 is 2650, density of MoS 2 /SiO 2 is 5060, thermal conductivity of C 2 H 6 O 2 is 0.253, thermal conductivity of MoS 2 is 1.5, thermal conductivity of MoS 2 /SiO 2 is 34.5, electrical conductivity of MoS 2 is 0.0005, electrical conductivity of Here, the Weissenberg number is The surface force in attendance of the Carreau-Yasuda liquid at the wall of the melting surface is The temperature gradient because of the nano and hybrid nanoparticles is Concentration gradient at the surface of the melting sheet is The local Reynolds number is Re = ax 2 ν f . Figure 2. Geometry of transport phenomenon in Carreau-Yasuda fluid.

Numerical Approach and Convergence Analysis
The finite element method is an effective method in the view of accuracy and the convergence of a problem compared with other numerical approaches. There are many advantages to FEM but some are discussed here: FEM has the ability to handle various complex geometries; This numerical method is thought to be most significant in solving physical problems with wide ranges; FEM requires less investment in the view of time and resources; A main advantage of FEM is its handling of various types of boundary conditions and It has a good ability with regards to the discretization (of derivative) problems into small elements; The Working scheme of finite element method has been shown with the help of Figure 4.
The numerical approach called finite element scheme is used to simulate the numerical results of highly non-linear PDEs and numerous applications of FEM are found in CFD (computational fluid dynamics) problems. The FEM approach is explained in the following steps: Step I: The division of a problem domain into a finite number of elements and residuals. The weak form is captured from the strong form due to residuals. The approximation result is simulated using shape functions and the approximation simulations of the variables are: Here F ξ is H and the shape function is defined as: Step II: In this step, the matrices are stiffness, vector and boundary (integral vector). The global stiffness (matrix) is obtained whereas the Picard (linearization approach) is utilized to obtain a linear system of equations that are defined as: Here F l and H l are variables (nodal values).
Step III: The algebraic equations (non-linear) resulting from the assembly process are: Step IV: The computational domain is considered as [0, 8] while mesh-free analysis is computed along with 270 elements. The problem is converged at mid of each of the  Table 1 reveals the study of convergence analysis. The solution to the problem is converged after simulations of 210 to 270 elements. It is observed that 270 elements are ensured for the convergence of the problem. Outcomes are provided for velocity, concentration and temperature at mid of each of the 270 elements. All numerical simulations related to tables and graphs are captured for the 270 elements.
Comparative analysis: The numerical result of the current problem is verified with published results [23] by the disappearing effects of We = = F r = H s = Ec = M = D f = ϕ 1 and ϕ 2 = 0. Numerical values of the Nusselt number are computed against the distribution in Prandtl number. Good agreements among the results of the present problem and published work [23] are presented in Table 2.

Results and Discussion
Mechanisms of velocity, thermal energy and diffusion of mass influenced by chemical reaction are addressed over a stretched melting surface. Correlations between silicon dioxide and Molybdenum dioxide in EG (ethylene glycol) are used in the presence of the Carreau-Yasuda liquid. Various kinds of influences (Soret, Dufour, viscous dissipation, Joule heating and magnetic field) are also addressed. As such, the complex transport phenomenon is simulated with the help of a numerical approach (FEM). The graphical computational investigations are captured in graphs and tables. The detailed outcomes are discussed below: Graphical investigations of velocity against distribution in various parameters: The change in Weissenberg, power law index, Forchheimer numbers and Carreau-Yasuda variables are addressed in the motion of fluid particles considered in Figures 5-8. Figure 5 is plotted to measure the role of We in the motion of hybrid nanoparticles. It is estimated that the motion of hybrid and fluid -nanoparticles is slowed down by applying higher We Values. The Weissenberg number is constructed in the current model due to the consideration of the rheology of the Carreau-Yasuda fluid while We is defined as a ration of elastic and viscous forces. An increase in We brings the declination in motion of fluid particles in the presence of nanoparticles and hybrid nanoparticles. Hence, a reduction is noticed versus the change in We. Moreover, the thickness of the momentum boundary layers decline when We is increased. The flow for Newtonian fluid is the dominated flow for a case of non-Newtonian fluid. The relationship between the velocity and power law index number is shown in Figure 6. The decreasing phenomenon of motion in fluid particles is captured and m is created due to tensor of the Carreau-Yasuda liquid. The numerical values of m are decided by the category of fluids (shear thinning or shear thickening). The fluid becomes thick in the case of large m values. Hence, the power law index number is not a significant parameter in the case of an enhancement in flow involvement of nanoparticles and hybrid nanoparticles. Parameter related to the power law number has a significant impact on adjusting the momentum boundary layer thickness. The role of F r is noticeable in the flow of nanoparticles and hybrid nanoparticles (see Figure 7). It is demonstrated that the parameter related to F r occurs in the momentum equation because of the Forchheimer porous. This kind of parameter behaves like a non-linear-type function in the flow of nanoparticles. In this case, the retardation force is created in fluid motion and brings resistance of the fluid particles into motion. Moreover, the thickness related to the boundary layer is reduced when large values of F r are applied. Further, the motion created by the Forchheimer porous is less than the motion created in the particles, excepting the involvement of Forchheimer porous media. The parameter associated with d is called the fluid variable and the change in d versus the velocity is captured in Figure 8.          Figures 9-12 reveal the characterization of thermal energy phenomenon versus the variation of H s , F r , Ec and D f . The thermal energy performance is measured with respect to the variation in heat generation number while this phenomena is shown in Figure 9. The production of heat energy is at its maximum when using higher values of H s . The external heat source at the sheet surface results in the maximum production into heat energy. to Thickness related to thermal layers is also enhanced due to the large values of heat generation number. Thus, the heat generation number is a significant parameter for the maximum production of heat energy. Moreover, an inclination in thermal energy is created due to the direct relation to thermal energy. It is demonstrated that the positive values for H s are present due to the phenomena of heat generation. As such, the impact of heat generation is visualized in our analysis. Figure 10 plots the enhancement in heat energy against the change in Eckert number. Physically, an enhancement in heat energy due to viscous dissipation is simulated. A direct relationship between viscous dissipation and the kinetic energy phenomenon was found. The temperature of the fluid particles is enhanced due to the Eckert number. The work-done rate is enhanced by the particle heat energy when viscous dissipation occurs. The relationship between F r and temperature profile is visualized in Figure 11. This figure captures the better performance in heat energy, including the appearance of F r . In fact, F r is generated due to the Forchheimer porous while F r generates more heat energy in the presence of nanoparticles and hybrid nanoparticles. Moreover, the concept of Dufour's number is characterized in the dimensionless heat equation due to the first law of thermodynamics while terms related to thermodynamics refer to the concept of Dufour (heat energy) because of the concentration gradient. The concentration gradient is enhanced using large values of D f . The fluid particles absorb more heat energy when D f is increased. Hence, the temperature profile is increased when there are higher values of D f (see Figure 12). Further, thickness of the thermal boundary layers are controlled by the impacts of F r and D f . Figures 13-15 have been plotted to visualize the transport of diffusion against the change in Sc, K c and Sr. The measurement of mass diffusion is captured in Figure 13, considering the influence of Sc. The diffusion of mass decreases when Sc is enhanced. We can observe that this reduction into mass diffusion happens due to the definition of Sc. Physically, Sc is rationed among mass and momentum diffusivities. According to the concept of Sc, the concentration curves are reduced when Sc is inclined. A similar situation occurs in terms of the boundary layer thickness in relation to concentration. Figure 14 visualizes the effect of the chemical reaction number on the transport of mass diffusion We found that the parameter related to K c revealed the coefficient of thermal energy along with the chemical reaction. The positive values of K c correspond to the destructive chemical reaction. In this case, a reduction is captured in the diffusion of the mass species. In the current flow model, the case related to a destructive chemical reaction is used. Thickness associated with the concentration layers is reduced when K c is increased. The decreasing graph is plotted between the Soret number and diffusion of mass as shown in Figure 15. The concept of Sr (fractioned between the difference in temperature and concentration) appears due to the temperature gradient in the concentration equation. Using Soret's theory, the solute diffusion is enhanced due to thermal energy.

Graphical investigations of mass diffusion against distribution in various parameters:
Mechanisms of gradient temperature, surface force and mass diffusion rate versus the distribution of various parameters: The computational analysis of surface force (skin friction coefficient), gradient temperature (Nusselt number) and rate of mass diffusion versus the variation of D f , Sc, H s , F r and We is simulated in Table 3. Table 3 reveals that the drag force (skin friction coefficient) declines when We and F r are increased. However, we see a reduction in the skin friction coefficient when the heat generation number is increased. When D f is enhanced, the constant variation is simulated in surface force at the wall. The temperature gradient decreases using the large values of D f , H s , F r and We. The parameter related to Sc plays a vital impact in maximizing the rate of thermal energy. As for the concentration gradient, it shows the same behavior as the temperature gradient.

Prime Consequences of the Problem
The transport features in the rheology of Carreau-Yasuda liquid and involvement of nanostructures and hybrid nanoparticles over a heated surface have been visualized. The Dufour and Soret effects under the action of a magnetic field have been addressed. Forchheimer porous media was also considered. The simulations of the current model were computed by finite element approach. The prime findings are captured below:

•
Convergence of the problem is ensured at 270 elements; • The motion of nanoparticles and hybrid nanoparticles in ethylene glycol is boosted versus the enhancement in fluid variable, power law index number, Weissenberg number and Forchheimer porous number;

Prime Consequences of the Problem
The transport features in the rheology of Carreau-Yasuda liquid and involvement of nanostructures and hybrid nanoparticles over a heated surface have been visualized. The Dufour and Soret effects under the action of a magnetic field have been addressed. Forchheimer porous media was also considered. The simulations of the current model were computed by finite element approach. The prime findings are captured below:

•
Convergence of the problem is ensured at 270 elements; • The motion of nanoparticles and hybrid nanoparticles in ethylene glycol is boosted versus the enhancement in fluid variable, power law index number, Weissenberg number and Forchheimer porous number; • Significant production of heat energy versus higher values of heat generation, Eckert, Dufour and Forchheimer porous numbers; • The transportation of solute particles declines versus the large values of Schmidt and chemical reaction numbers, but solute particles accelerate against higher values of Soret number; • Surface force is increased via large values of Weissenberg and Forchheimer porous numbers but surface force is decreased versus the large values of heat generation number and The role of the Schmidt number is significant in the development of temperature and concentration gradient.