Numerical Study on Generalized Heat and Mass in Casson Fluid with Hybrid Nanostructures

The rheological model for yield stress exhibiting fluid and the basic laws for fluid flow and transport of heat and mass are used for the formulation of problems associated with the enhancement of heat and mass due to dispersion of nanoparticles in Casson. The heat and mass transfer obey non-Fourier’s laws and the generalized Fick’s law, respectively. Model problems are incorporated by thermal relaxation times for heat and mass. Transfer of heat energy and relaxation time are inversely proportional, and the same is the case for mass transport and concentration relaxation time. A porous medium force is responsible for controlling the momentum thickness. The yield stress parameter and diffusion of momentum in Casson fluid are noticed to be inversely proportional with each other. The concentration gradient enhances the energy transfer, and temperature gradient causes an enhancement diffusion of solute in Casson fluid. FEM provides convergent solutions. The relaxation time phenomenon is responsible for the restoration of thermal and solutal changes. Due to that, the thermal and solutal equilibrium states can be restored. The phenomenon of yield stress is responsible for controlling the momentum boundary layer thickness. A porous medium exerts a retarding force on the flow, and therefore, a deceleration in flow is observed. The thermal efficiency of MoS2−SiO2− Casson fluid is greater than the thermal efficiency of SiO2− Casson fluid.


Introduction
The diverse rheological behavior of non-Newtonian fluid motivated researchers to propose diverse rheological models. Each model captures some specific rheological features. Casson rheological model is a non-Newtonian rheological model and captures yield stress rheological features. This means, to some extent, Casson fluid behaves like a solid and, after some specific value of applied stress, it starts behaving like a fluid. This specific value of applied stress is called yield stress. The simplest rheological model for the incompressible flow of Casson fluid is given by where β is the Casson fluid parameter and for β → ∞ , Equation (1) reduces to the rheological tensor of an incompressible Newtonian fluid. Numerical studies on various aspects related to Casson fluid have been published so far. For example, Hamid et al. [1] wrote a study on the existence of dual solutions for the problem related to heat transfer in Casson fluid over the stretchable surface. Nadeem et al. [2] discussed the 3D flow of Casson fluid induced by wall moving with the exponential surface. Mukhopadhyay et al. [3] theoretically studied the transfer of heat energy in Casson fluid over a surface moving with nonlinear velocity. Hayat et al. [4] examined the effects of temperature gradient and concentration gradients on mass and heat transfer in Casson fluid, respectively. Khan et al. [5] examined the simultaneous effects of homogeneous-heterogeneous chemical reactions on mass transfer in Casson fluid. The author has carefully reviewed the literature and has come to know that no study using hybrid nanostructures (MoS 2 and SiO 2 ) in Casson fluid with the novel model (given by Equations (2)- (5)) is published yet. The energy crisis is the biggest challenge to humankind. Energy storage and its efficient utilization is the need of the hour. Techniques for minimizing energy losses are highly required. With this need in mind, researchers working in the field of thermal and cooling systems have been striving for introducing novel techniques for efficient transport of heat energy and the efficient usage of stored energy. Recently, a technique of inclusion of nanosized solid structures in the fluid has been introduced to obtain the fluid mixtures of very good thermal conductivity. Such mixtures serve as better working fluids, and consequently, the thermal performance of the fluid increases. Because of numerous applications of fluids with nanoparticles, researchers have studied much of the dynamics of such fluids. For instance, Sheikholeslami et al. [6] investigated the effects of nanoparticles Fe 3 O 4 on heat transfer in water. Sheikholeslami et al. [7] conducted finite element numerical experiments for the analysis of thermal enhancement in fluid subjected to the variable magnetic field. Dogonchi et al. [8] also analyzed the impact of nanoparticles and thermal radiations on the transfer of heat in the mixed convection flow of fluid over a moving body. Dogonchi et al. [9] studied the role of Cu− nanoparticles on the enhancement of heat energy in the presence of homogeneous-heterogeneous chemical reactions in water subjected to the uniform heat source. Ayub et al. [10] used novel heat flux models for the modeling of heat transfer in viscoelastic fluid with nanoparticles. Rana et al. [11] used Koo-Kleinstreuer and Li models with novel heat flux theory to model heat transfer in Sutter by fluid with nanoparticles. Studies [6][7][8][9][10][11] and the references therein are limited to thermal enhancement via mono-nanoparticles. The recent advancement in the field reveals that the better method is the inclusion of nanoparticles of more than one kind. Therefore, several studies have been published. The reasoning and justification are given below.
The advancement in technologies and success in the synthesis of nanosized particles has brought a great revolution for introducing a class of special class of fluids called hybrid nanofluids. These are fluids that are homogeneous mixtures of nanoparticles of more than one kind. It is a proven fact that hybrid nanofluids are better coolants in comparison to nanofluids. Due to this factual reason, hybrid nanofluids have been frequently studied. The recent advancement in the field of thermal enhancement has revealed that the dispersion of nanoparticles of more than one kind is more effective than the dispersion of nanoparticles of a single kind. The fluid with more than one nanoparticle is a hybrid nanofluid. Several studies on this topic have been published. For example, Ahmed et al. [12] studied the impact of hybrid nanoparticles on heat transfer in a pulsating flow exposed to the magnetic field and thermal radiations. Iftikhar et al. [13] discussed the role of Cu and SiO 2 in enhancing heat transfer in water. Benkhedda et al. [14] studied the importance of hybrid nanoparticles in enhancing thermal transfer in fluid in a pipe. Aziz et al. [15] examined the role of hybrid nanoparticles on entropy generation in Erying-Powell fluid subjected linear thermal radiations and viscous dissipation. Subhani et al. [16] theoretically analyzed the impact of hybrid nanoparticles on heat-transfer enhancement in the micropolar fluid. Wiani et al. [17] discussed the role of hybrid nanoparticles on heat transfer in fluid moving on a curved surface. Kaneez et al. [18] numerically analyzed the influence of hybrid nanoparticles on thermal enhancement in Casson fluid with dust particles. It is important to mention that this study on an enhancement of heat in Casson fluid with hybrid nanoparticles. However, the present study is about the impact of hybrid nanoparticles on heat and mass transport in Casson fluid with hybrid nanoparticles. Further, the present work considers the flow and heat transfer in Casson fluid over a vertical surface. In this case, the buoyance force is significant and is considered here. This is another aspect due to which the present work differs from the work of Kaneez et. al. [18]. 3 of 18 The electrically conducting fluids moving in the presence of magnetic fluid have been found in many natural and man-made applications. The flows of such fluids are called MHD flows and have been investigated a lot in the present era. For instance, Nawaz et al. [19] examined the impact of magnetic fields on the flow of micropolar fluid composed of ionized particles. Nawaz et al. [20] studied the magnetic field on the mass and heat diffusion in the Newtonian fluid between two walls moving with nonuniform velocity in the radial direction. Asif et al. [21] investigated the thermal performance of the fluidexposed magnetic field. In light of these studies [19][20][21], the authors were inspired and considered how magnetic field applied to the flow because Casson fluid is an electrically conducting fluid and experiences Lorentz force arises due to the change in magnetic flux. Flows over vertical boundaries are greatly affected by a buoyancy force that arises due to density differences caused by temperature and concentration differences. Downward flows experience positive a buoyancy force, whereas upward flows experience a negative buoyancy force. Momentum equations are modified in terms of buoyancy force. Such flows are called mixed convective flows over moving boundaries. Such flows have been studied intensely in recent years. The suitable references are listed [22][23][24][25][26][27].
Brownian motion and thermophoretic effects affect flow characteristics and, therefore, are neglected. Several kinds of research have discussed these effects. For instance, Sheikholeslam et al. [28] studied the influence of thermophoresis and Brownian motion on heat and mass in a mixed convection flow of electrically conducting fluid containing copper nanoparticles and subjected to the magnetic field. Kandasamy et al. [29] also studied the influence of thermophoresis and Brownian motion on the simultaneous impact of heat and mass transfer subjected to the thermal stratification in fluid containing nanoparticles. Makinde et al. [30] examined the role of thermophoresis and Brownian motion on the biconvective flow of MHD fluid over a surface of paraboloid revolution. Lin et al. [31] discussed the influence of thermophoresis and Brownian motion on heat transfer in fluid in a groove rotating with constant angular velocity. Some further applications and related aspects are mentioned in recent published work. Non-newtonian casson fluid flow over a stretching, exponentially embedded surface was examined by Animasaun, I et al. [32] for the laminar free convective MHD boundary layer. The effect of all of the influencing factors on the velocity and temperature profiles of the flow at the boundary was studied with an approximate solution (HAM) approach. Alruwashid et al. [33] analyzed the effects of graphene concentration on the electrochemical properties of cobalt ferrite nanocomposite materials. Pal et al. [34] developed novel nanostructures using pentablock polymers. Nasr et al. [35] introduced the technique of surface enhancement for nanosheets. Zemtsova et al. [36] published a comprehensive study on the feature of the synthesis of the disposed TiC with a nickel nanostructures.
A comprehensive review of the literature has shown that no study on the impact of hybrid nanoparticles (MOS 2 − SiO 2 ) on thermal performance of Casson under Brownian and thermophoretic effects has been conducted so far. Secondly, the researches on thermal enhancement are conducted numerically using numerical techniques other than FEM. FEM is a powerful method and is less constrained. It converges rapidly, and the accuracy of results is quite sound. The stated problems are solved numerically by FEM. Further, no author has discussed the Brownian motion and thermophoretic mechanism in the presence of hybrid nanofluid so far. This means the impact of dispersion of nanoparticles on simultaneous heat and mass transfer has not been studied yet. This research is arranged in five sections. Section 1 is related to the background. Section 2 contains the formulation of problems. Section 3 discusses numerical methodology. Section 4 is about the discussion of results based on numerical simulations. Section 5 contains the key results of this study.

Formulation of Problem
The boundary layer approximations give the following governing models: ∂u ∂x The conditions implemented to solve the above system of PDEs are where [u, v, 0] is the velocity; g is the gravitational acceleration; ρ is the density; µ is the kinematic viscosity and c p specific heat constant; k is the thermal conductivity; the subscript hn f stands for hybrid nanofluid; D is the mass diffusion coefficient; τ 0 is the thermal relaxation time; τ 1 is the concentration relaxation time; and γ 0 is the parameter associated with thermal slip. The flow diagram of the problem is represented in Figure 1. The following model for nanoparticles is incorporated Equations (2)-(7) in dimensionless form can be changed via the following change of variable: Hence, one obtains: The associated boundary conditions are dimensionalized, and hence, one obtains: where is the wall shear stress and = is the Reynolds number. Thermophysical properties are tabulated in Table 1.

Numerical Procedure
The different numerical schemes can be applied for the solution of Equations (8)-(12 However, FEM is powerful and provides the most accurate numerical solutions. The ma steps are mentioned below. Step 1: Residuals are multiplied by weight functions and are integrated over th Derivatives concerning η are denoted by prime. M is a magnetic parameter; (Gr) t and (Gr) c are the Grashof numbers for temperature and concentration, respectively; β * is heat generation parameter; K is the porous medium parameter; Pr is the Prandtl-number; Sc is the Schmidt number; Ec is the Eckert number; γ is the thermal relaxation time parameter; γ 1 is the concentration relaxation time parameter; and B i and E i are Biot numbers. These are expressed as The skin friction and Nusselt and Sherwood numbers are obtained as where τ w is the wall shear stress and Re x = ax 2 ν f is the Reynolds number.
Thermophysical properties are tabulated in Table 1. Table 1. Values of properties of pure fluid and nanostructure.

Numerical Procedure
The different numerical schemes can be applied for the solution of Equations (8)- (12). However, FEM is powerful and provides the most accurate numerical solutions. The main steps are mentioned below.
Step 1: Residuals are multiplied by weight functions and are integrated over the typical element.
Linear weight and shape functions are given by The dependent unknowns are approximated over the element [η e , η e+1 ] by the finite element approximations where f j , h j , θ j and φ j are to be computed. S j is the shape function.
Step 2: The weighted integrals are evaluated only for second-order linear terms to obtain weak forms of residuals.
Step 3: The weak forms are approximated via Galerkin approximations.
where f and h are nodal values that are computed.
Step 4: Stiffness elements are derived, and the stiffness matrix is constructed under the assembly procedure.
The nonlinear system is linearized by Picard linearization proposed as where {π} r−1 is the nodal value at (r − 1)th the iteration and {π} r is the nodal value at the rth iteration.
Step 6: A grid-independent study is performed.
The numerical solutions are also checked to be grid independent. This mesh-free analysis is recorded in Table 2. The recorded values in Table 2 show that the numerical results become grid independent when the computational domain [0, 12] is broken into two elemental line segments. Table 2. Local skin friction coefficient, Nusselt number and Sherwood number for hybrid nanofluid (MOS 2 − SiO 2 ) and nanofluid (SiO 2 ) when M = 0.9, Pr = 0.9, Sc = 0.8, (Gr) c = 0.1, (Gr) t = 0.2,

Results and Discussion
The yield stress rheological model, non-Fourier law of heat conduction, and generalized Fick's law are simultaneously used in the formulation of problems dealing with the simultaneous transfer of heat energy and mass in Casson fluid with hybrid nanoparticles over a vertical surface. The buoyancy forces due to the density differences under Boussinesq approximations are considered of a significant order of magnitude, and Lorentz force is measured from Ohm's law when the magnetic Reynolds number is small. The numerical experiments with three samples of parametric values are performed. The simulations are visualized, and the outcomes are displayed in the form of graphs and numerical data. The mesh-free analysis is performed, and the results are validated. A detailed discussion about visualized outcomes is provided below.
Fluid particles motion and magnetic field: The Lorentz force is incorporated in the momentum equation through Ohm's law and the term involved Hartmann number in the dimensionless equation (Equation (8)). This dimensionless number (M) is capable of determining the impact of the change of intensity of magnetic field on the motion of fluid particles. Since the Lorentz force is negative (in this case), an increase in M through values

Fluid particles motion and resistance by the porous medium:
The parameter K is called the porous medium parameter, and its variation through increasing values ( K = 0.05, 0.1, 0.15, 0.2) depicts the decrease in voids between the particles of the porous medium. Therefore, the fluid particles of MoS 2 − SiO 2 − Casson f luid and SiO 2 − Casson f luid face more a resistive force by the medium, and as a consequence, the motion of the particles of the fluid slows down. This fact is also visualized in numerical simulations (see Figure 3). It is also observed that MoS 2 − SiO 2 − Casson f luid experiences more resistive force than SiO 2 − Casson f luid. Numerical experiments have also predicted that the presence of a porous medium in the fluid regime may play a vital role in controlling the boundary layer thickness.
Fluid particle motion and buoyancy force: The parameters (Gr) c and (Gr) t are the coefficients of terms associated with buoyancy forces that arise due to the density differences caused by the compositional and temperature differences under Boussinesq approximation. For (Gr) c > 0 and (Gr) t > 0, the buoyancy force is posture and assists the flow (see Figures 4 and 5), whereas when (Gr) c < 0 and (Gr) t < 0, the Buoyancy force is negative and opposes the flow. This opposing behavior can be noticed during numerical experiments. However, the observations are not included here (due to the restrictions on the length of the article). It is also found that the buoyancy force in MoS 2 − SiO 2 − Casson f luid is stronger than the buoyancy forces in SiO 2 − Casson f luid.
Casson fluid face more a resistive force by the medium, and as a consequence, the motion of the particles of the fluid slows down. This fact is also visualized in numerical simulations (see Figure 3). It is also observed that − − Casson fluid experiences more resistive force than − Casson fluid. Numerical experiments have also predicted that the presence of a porous medium in the fluid regime may play a vital role in controlling the boundary layer thickness.

Temperature of fluid particles and intensity of magnetic field:
The MoS 2 − SiO 2 − Casson and SiO 2 − Casson fluids are electrically conducting fluids. When electric current passes through MoS 2 − SiO 2 − Casson and SiO 2 − Casson fluids, some of the electrical energy converts into heat, and this phenomenon is called Ohmic dissipation. This phenomenon is also considered here in this study. The term in dimensionless energy equation involving M is due to Joule heating (Ohmic heating) and determines the impact of Ohmic dissipation on the temperature of the fluid. This dissipated heat raises the kinetic energy of the fluid particles, and therefore, their temperature increases. This effect can be noticed in Figure 6. This figure also demonstrates that the Ohmic dissipation phenomenon in MoS 2 − SiO 2 − Casson f luid is stronger than that in SiO 2 − Casson. Therefore, the rise in temperature of fluid particles of MoS 2 − SiO 2 − Casson is greater than the rise of temperature of fluid particles of SiO 2 − Casson f luid. Therefore, the influence of Ohmic dissipation on boundary layer thickness associated with the flow of MoS 2 − SiO 2 − Casson f luid is greater than the influence of Ohmic dissipation on boundary layer thickness associated with the flow of SiO 2 − Casson f luid.

Temperature of fluid particles and viscous dissipation:
The friction among the particles of fluid causes the dissipation of heat which is equal to the rate at which the work is performed by friction force. This effect is called viscous dissipation. It is noticed during simulations that the inclusion of hybrid nanostructures ( and ) increases the friction among the particles of the fluid. It can also be observed that the rise in friction among the fluid particles due to the inclusion of hybrid nanoparticles ( − ) is greater than the rise in friction among the particles of fluid due to the inclusion of in fluid (see Figure 7).

Temperature of fluid particles and viscous dissipation:
The friction among the particles of fluid causes the dissipation of heat which is equal to the rate at which the work is performed by friction force. This effect is called viscous dissipation. It is noticed during simulations that the inclusion of hybrid nanostructures (MoS 2 and SiO 2 ) increases the friction among the particles of the fluid. It can also be observed that the rise in friction among the fluid particles due to the inclusion of hybrid nanoparticles (MoS 2 − SiO 2 ) is greater than the rise in friction among the particles of fluid due to the inclusion of SiO 2 in fluid (see Figure 7).

Temperature of fluid particles and Biot number:
The Biot number B i appears in the dimensionless form of convective boundary conditions. The role of B i on the temperature of fluid particles of MoS 2 − SiO 2 − Casson f luid is examined. The comparative analysis between the role of B i on temperatures of SiO 2 − Casson f luid given by Figure 8. In both types of fluids (mono-nanofluid and hybrid nanofluid), temperature fluid particles have an increasing behavior for Biot number B i .

Temperature of fluid particles and Prandtl number:
The Pr is the dimensionless parameter, and an increase in the Pr relates to a decrease in thermal conductivity corresponding to the lower ability of fluid to conduct heat. Consequently, the temperature decreases. This decrease in temperature of MoS 2 − SiO 2 − Casson f luid and SiO 2 − Casson f luid can be noticed (see Figure 9). The decreasing impact of Pr on the temperature of SiO 2 − Casson f luid is greater than the impact on the temperature of MoS 2 − SiO 2 − Casson f luid. performed by friction force. This effect is called viscous dissipation. It is noticed during simulations that the inclusion of hybrid nanostructures ( and ) increases the friction among the particles of the fluid. It can also be observed that the rise in friction among the fluid particles due to the inclusion of hybrid nanoparticles ( − ) is greater than the rise in friction among the particles of fluid due to the inclusion of in fluid (see Figure 7).

Temperature of fluid particles and Prandtl number:
The is the dimensionless parameter, and an increase in the Pr relates to a decrease in thermal conductivity corresponding to the lower ability of fluid to conduct heat. Consequently, the temperature decreases. This decrease in temperature of − − and − can be noticed (see Figure 9). The decreasing impact of on the temperature of − is greater than the impact on the temperature of − − .

Temperature of fluid particles and Prandtl number:
The is the dimensionless parameter, and an increase in the Pr relates to a decrease in thermal conductivity corresponding to the lower ability of fluid to conduct heat. Consequently, the temperature decreases. This decrease in temperature of − − and − can be noticed (see Figure 9). The decreasing impact of on the temperature of − is greater than the impact on the temperature of − − .

Temperature of fluid particles and thermal memory effects:
The parameter is called the thermal relaxation parameter and has an analogy with the Deborah number proposed for viscoelastic fluid. This parameter determines the restoring characteristics of thermal equilibrium. It is important to mention here for = 0, the non-Fourier's law re-

Temperature of fluid particles and thermal memory effects:
The parameter γ is called the thermal relaxation parameter and has an analogy with the Deborah number proposed for viscoelastic fluid. This parameter determines the restoring characteristics of thermal equilibrium. It is important to mention here for γ = 0, the non-Fourier's law reduces to the classical Fourier's law, and an increase of γ replicates the increase in ability of fluid to restore thermal equilibrium. Therefore, the temperature of fluid decreases (see Figure 10).

Concentration of fluid particles and Schmid number (Sc):
The Sc is the number that determines the impact of diffusion coeffect on the diffusion solute in the fluids. Here, in this work, the behavior of Sc on the diffusion of solute in MoS 2 − SiO 2 − Casson f luid and SiO 2 − Casson f luid is examined, and the outcomes are presented by Figure 11.
Concentration of fluid particles and solutal relaxation time: The effects of thermal relaxation time on the temperature of fluid particles are similar to the impact of solutal relaxation time on the concentration field (compare Figures 10 and 12).
Behaviors of wall shear rate of heat transfer and mass flux against selected parameters: The behaviors of selected parameters Bi, Pr, (Gr) t , γ and γ 1 on wall shear rate, heat-transfer rate and mass-transfer rate are examined. The results so obtained are recorded in the tabular numerical data (see Table 4) for both SiO 2 − Casson f luid (nanofluid) and MoS 2 − SiO 2 − Casson f luid. This table clearly depicts that heat transfer is enhanced due to simultaneous dissipation of nanostructures MoS 2 and SiO 2 . Biot number Bi compels the wall shear stress to decrease in both types of fluids (nano and hybrid nanofluids). It is noticed that wall shear stress for the case of hybrid nanofluid (MoS 2 − SiO 2 − Casson f luid) is greater than that for the case of nanofluid (SiO 2 − Casson f luid). An increasing rate of Bi on wall heat-transfer rate is observed. However, it has a decreasing trend for the case of mass transfer rate for both SiO 2 − Casson and MoS 2 − SiO 2 − Casson f luids. Pr has shown an increasing behavior on wall shear stress and wall heat-transfer rate. However, the opposite trend is noticed for mass flux. For positive values of Grashof number (Gr) t , the wall shear stress decreases; however, opposite trend is noticed for (Gr) t < 0. It is also found that for the assisting flow ((Gr) t > 0), the rate of heat transfer and rate of mass transfer increase. The reverse behavior is noted for (Gr) t < 0. The trend of shear stress, mass flux and the rate of heat transfer for various values of γ and γ 1 can be observed in Table 4.

Concentration of fluid particles and Schmid number ( ):
The is the number that determines the impact of diffusion coeffect on the diffusion solute in the fluids. Here, in this work, the behavior of on the diffusion of solute in − − and − is examined, and the outcomes are presented by Figure 11.

Behaviors of wall shear rate of heat transfer and mass flux against selected parameters:
The behaviors of selected parameters , , ( ) , and on wall shear rate, heattransfer rate and mass-transfer rate are examined. The results so obtained are recorded in the tabular numerical data (see Table 4 is greater than that for the case of nanofluid ( − Casson fluid). An increasing rate of on wall heat-transfer rate is observed. However, it has a decreasing trend for the case of mass transfer rate for both − Casson and − − Casson fluids . has shown an increasing behavior on wall shear stress and wall heat-transfer rate. However, the opposite trend is noticed for mass flux. For positive values of Grashof number ( ) , the wall shear stress decreases; however, opposite trend is noticed for ( ) < 0. It is also found that for the assisting flow (( ) > 0), the rate of heat transfer and rate of mass transfer increase. The reverse behavior is noted for ( ) < 0. The trend of shear stress, mass flux and the rate of heat transfer for various values of and can be observed in Table 4.

Conclusions
The formulated models are solved numerically with the implementation of FEM. The Cattaneo-Christov generalized laws for heat and mass diffusion are helpful for analyzing thermal and solutal relaxation phenomena. The numerical samples were obtained for various values of related parameters to visualize the dynamics of flow quantities, including the rate of heat and mass transport flux. The main observations related to this study are 1.
FEM provides convergent solutions.

2.
The relaxation time phenomenon is responsible for the restoration of thermal and solutal changes, due to which, the thermal and solutal equilibrium states can be restored.

3.
The phenomenon of yield stress is responsible for controlling momentum boundary layer thickness.

4.
Porous medium exerts retarding force on the flow, and therefore, the deceleration inflow is observed. 5.
The thermal efficiency of MoS 2 − SiO 2 − Casson fluid is greater than the thermal efficiency of SiO 2 − Casson fluid.

6.
The change in magnetic flux by the motion of the fluid particles results in a Lorentz force that opposes the flow, and the boundary layer thickness becomes shorter. It can be noticed also noticed that the Lorentz force associated with the motion of SiO 2nanofluid has a smaller magnitude relative to the Lorentz force induced by the motion of hybrid nanofluid (MoS 2 -SiO 2 -nanofluid). Thus, the viscous region for momentum boundary layer flow nanofluid is wider than that for the flow of hybrid nanofluid. 7.
The vortex viscosity is responsible for decreasing the macromotion of fluid particles. Thus, velocity of Newtonian fluid is greater than the velocity of micropolar fluid. Further, the impact of vortex motion of solid structures in the fluid on velocity of nanofluid is negligible in comparison with the case of hybrid nanofluid. 8.
Joule heating associated with hybrid nanofluid flow has greater value in comparison with Joule heating in nanofluid. Thus, the temperature of hybrid nanofluid increases more rapidly than the temperature of nanofluid. Moreover, the Ohmic heating process in hybrid nanofluid is stronger than that in nanofluid. 9.
Viscous dissipation results an increase in temperature of the fluids (both mono-and hybrid nanofluids). 10. The thermal relaxation time is responsible for decreasing the temperature of fluid. This is because of the ability of the fluid to restore thermal changes in order to attain thermal equilibrium.
Future aspects: In this study, MoS 2 and SiO 2 are used to enhance the thermal conductivity of the Casson fluid, and the optimized heat transfer is noticed. Based on the simulation, it is claimed that thermal performance is significantly increased. This observation is related to Casson fluid with the nanoparticles. However, it is very important to mention that this analysis should be performed for other non-Newtonian fluids. A similar analysis is also required for other nanoparticles. The comprehensive comparison between present study and the studies for other nanoparticles in other non-Newtonian fluids will give a comprehensive picture. Authors have the intension to do so in the future.