Mixed Convection of Silica–Molybdenum Disulphide/Water Hybrid Nanoliquid over a Rough Sphere

: A steady combined convective motion over a rough sphere with hybrid nanoparticles is analyzed. We have considered silica (SiO 2 ) and molybdenum disulphide (MoS 2 ) nanoadditives which are added in H 2 O to form MoS 2 –SiO 2 /H 2 O hybrid nanoliquid. The partial differential equations describing the boundary layer ﬂow characteristics are reduced into non-dimensional form with appropriate non-similar reduction. It should be noted that the governing equations have been written using the conservation laws of mass, momentum and energy. These considered equations allow simulating the analyzed phenomenon using numerical techniques. Implicit ﬁnite difference approximation and technique of Quasilinearization are utilized to work out the dimensionless control equations. The inﬂuence of various physical characteristics included in this challenge, such as the velocity ﬁelds and temperature patterns, is investigated. The study of border gradients is performed, which deals with the skin friction and energy transport strength. The plots of computational outcomes are considered, which ascertain that velocity distribution reduces, whilst coefﬁcient of friction at the surface, energy transport strength and temperature distribution augment for enhancing values of hybrid nanoﬂuid. For enhancing magnitude of combined convection parameter, dimensionless velocity distribution, surface drag coefﬁcient and energy transport strength enhance, while temperature distribution diminishes. High impact of hybrid nanoﬂuid on energy transport strength and the surface friction compared to the host liquid and mono nanoﬂuid in presence/absence of surface roughness is shown. Velocity distribution enhances for rising values of velocity ratio parameter. Enhancing values of frequency parameter rise the friction at the surface and energy transport strength. It is also examined that the hybrid nanoﬂuid has a maximum temperature for the blade-shaped nanoparticles and has a low temperature for the spherical-shaped nanoparticles. for rising values of Ri , which increases the energy transport strength from wall to ﬂuid. Further, for enhancing values of volume fraction ( ϕ 2 ), the friction at the surface and the energy transport strength enhance. This is because addition of nanoparticles in the host liquid leads to intensifying its temperature and as a result one can ﬁnd an intensiﬁcation of the energy transport.


Introduction
The transport behavior of nanofluids can be described using, for example, Buongiorno's model [1] and Tiwari and Das approach [2]. In Buongiorno's approach, the nanoparticles' slip velocity is not zero considering seven possible slip phenomena including thermophoresis, Magnus effect, inertia, Brownian motion, diffusiophoresis, gravity and fluid drainage. Tiwari and Das model [2] can be considered as an effective model with formulae for heat conductivity and dynamic viscosity that are greater than in the case of regular liquids. Therefore, Tiwari and Das approach can be considered as a homogeneous model which is different from Buongiorno's model.
Nanofluid is a combination of host fluid and nano-sized solid particles [3]. Scientists are in search of better type of liquid in order to enhance the thermal efficiency of different engineering systems. Taking into account such investigation, the hybrid nanofluid as a combination of base liquid and several different types of nanoparticles, can have high thermal conductivity compared to mono nanofluid and hybrid nanofluids can be considered as the advanced forms of nanofluid. Hybrid nanofluids have heat transfer applications in numerous engineering fields such as, manufacturing, naval structures, transportation, acoustics, defense, micro fluidics, medical etc. [4][5][6]. Thus, Waini et al. [7] have analyzed the combined convection hybrid nanoliquid circulation about a vertical border and their results reveal that the energy transport strength diminishes in shrinking case. In the work of Dinarvand et al. [8], hybrid nanofluid flow past a wedge by considering titanium and copper nanoparticles has been analyzed. Rostami et al. [5] have scrutinized the combined convection of a hybrid silica-alumina nanofluid and their results ascertain that the hybrid nanofluid has high energy transport strength than the base fluid. Moreover, many other researchers have worked by considering the different combination of nanoparticles for hybrid nanofluids such as, Ag-TiO 2 , CuO-MgO, Cu-Al 2 O 3 and MWCNTs-Fe 2 O 4 -SWCNTs, etc. [9][10][11][12][13]. Since, many authors such as Maraj et al. [14], Ghadikolaei et al. [15], Khan et al. [16], etc., have studied the water containing Molybdenum disulphide and Silica nanoparticles, considering the either of one or both nanoparticles. Chemical attributes of molybdenum disulfide including capabilities of lubrication, heat conductivity and thermal capacitance are employed in various engineering systems [17]. In addition, Silicon oxide (silica) is widely applied in energy transport applications and it can be successfully use as nano-inclusions [18]. The work of these authors has motivated us to perform the present research and therefore, we have investigated the current problem.
Some kinds of nanoadditives have been employed by some scientists in theoretical and experimental researches. Among the various types of nanoadditives, one of them is known as molybdenum disulfide nanoadditives. Such nanoadditives are rarely employed in the papers particularly in theoretical researches. MoS 2 is characterized by large band gaps structure which is related to the structure of grapheme and as a result, MoS 2 is employed in different engineering systems [19][20][21][22][23][24][25][26]. Moreover, chemical attributes of MoS 2 can be applied in various mechanical systems [23]. MoS 2 particles have been coated with copper to improve their combination with the copper particles for the sintering [24]. Further, silicon is widely spread element on earth after oxygen [25]. Silicon has very useful attributes that illustrates many different applications of this element in various fields. Silicon chips have essentially changed different technical systems including electronics, space and aviation industry, and others. Silicon is used in photovoltaic panels that characterize the high efficiency and suitability. Silicon oxide (silica) is widely employed in energy transport applications and it can be successfully applied as a nanoadditive [18].
In many branches of science and technology, the notion of mixed convection flow can be applied. Correspondingly, combined convection flows play a crucial role in the operation of thermal processing, such as in heated or cooled enclosures and electrical power supply, ventilation and air conditioning systems in atmospheric and ocean circulation, electronic equipment and thermal distribution in buildings [27][28][29]. In particular, owing to its prevalence in engineering fields, the combined convection flow around the sphere has been stimulated by several researchers. The examination of heat transport in combined convective motion over a sphere helps in polymer and fiber industries. Different investigators have studied a combined convective boundary layer motion over a sphere. Combined convective motion around a solid sphere was studied by Salleh et al. [30], taking Newtonian warming into the consideration and their outcomes exhibit that the dimensionless velocity and temperature patters can be reduced with Pr. In addition, Rashad et al. [31] examined the combined convective motion past a sphere and authors have ascertained that the chemical effects are noticeable for energy transport strength. Further, many researchers also analyzed the combined convective water boundary layer flow over a sphere [32,33].
Many scientists have also examined the hybrid nanofluid circulation over different geometries. Recently, Waini et al. [34] have scrutinized combined convective motion with Cu-Al 2 O 3 /H 2 O nanoliquid past a surface and their results reveal that the surface drag coefficient lessens for opposing buoyancy flow, while it boosts for assisting buoyancy flow. Further, Muhammad et al. [35] have analyzed combined convection over a curved surface with MWCNT-Cu/H 2 O hybrid nanofluid and the obtained outcomes ascertain that the fluid velocity enhances for combined convection parameter. Bognar et al. [36] have studied the nanofluid flow over a uniformly moving plate considering different nanoparticles and their outcomes reveal that the values of velocity and temperature are greater for alumina nanoparticles compared to titania and magnetite particles. In addition, Bognar and Hriczo [37] have studied the ferrofluids behavior over a flat sheet under the variable magnetic field influence. They have revealed that a rise of the ferromagnetic parameter characterizes the velocity reduction and temperature augmentation. Many researchers have also worked on the combination of silica and molybdenum nanoparticles by considering various base fluids such as water, ethylene glycol, etc. [38,39].
The one feature that tests micro-irregularities presented in the surface texture, is the surface roughness. The roughness analysis was initiated by Nikuradse [40] where he explored the roughness effects on surface drag coefficient. Owen and Thompson [41] have examined the influence of roughness on energy transport strength, while Dawson and Trass [42] have investigated its influence on Sherwood number. Sedahmed and Shemilt [43] have examined the forced convection behavior owing to the border roughness. Patil et al. [44] have investigated the nanosuspension combined convection motion over a rough sphere and their outcome reveals that with suitable choice of roughness and combined convection parameter the point of separation can be delayed. In addition, Patil et al. [45] have analyzed nanoliquid motion about a moving plate and their outcome shows that the friction at the surface enhances.
Taking into account this brief literature review, it can be found that investigation of the combined convective flow of molybdenum silica/water hybrid nanofluid about a rough sphere has not been studied earlier. The novelty of the present analysis includes the following points: -combined convective flow over a rough sphere; -impact of molybdenum-silica/water hybrid nanofluid on flow structures and heat patterns; -effect of surface roughness on hydro-and thermodynamics.
The modeled governing equations and their boundary constraints are affected to the appropriate non-similar reduction. Further, the obtained relations are attempted for solution transformations by employing the method of Quasilinearization and the implicit finite difference approximation [44,45].

Mathematical Analysis
Herein, we investigate the viscous, laminar, incompressible combined convection motion around a sphere. Here R is the sphere size. Let T w and T ∞ be the temperature at the border and away from the surface. The geometry and the configuration of the flow system are schematically presented in Figure 1. Here x is the curvilinear coordinate over the sphere border, while coordinate y is perpendicular it. The section r(x) explains the contour of the considered object. Further, Muhammad et al. [35] have analyzed combined convection over a curved surface with MWCNT-Cu/H2O hybrid nanofluid and the obtained outcomes ascertain that the fluid velocity enhances for combined convection parameter. Bognar et al. [36] have studied the nanofluid flow over a uniformly moving plate considering different nanoparticles and their outcomes reveal that the values of velocity and temperature are greater for alumina nanoparticles compared to titania and magnetite particles. In addition, Bognar and Hriczo [37] have studied the ferrofluids behavior over a flat sheet under the variable magnetic field influence. They have revealed that a rise of the ferromagnetic parameter characterizes the velocity reduction and temperature augmentation. Many researchers have also worked on the combination of silica and molybdenum nanoparticles by considering various base fluids such as water, ethylene glycol, etc. [38,39]. The one feature that tests micro-irregularities presented in the surface texture, is the surface roughness. The roughness analysis was initiated by Nikuradse [40] where he explored the roughness effects on surface drag coefficient. Owen and Thompson [41] have examined the influence of roughness on energy transport strength, while Dawson and Trass [42] have investigated its influence on Sherwood number. Sedahmed and Shemilt [43] have examined the forced convection behavior owing to the border roughness. Patil et al. [44] have investigated the nanosuspension combined convection motion over a rough sphere and their outcome reveals that with suitable choice of roughness and combined convection parameter the point of separation can be delayed. In addition, Patil et al. [45] have analyzed nanoliquid motion about a moving plate and their outcome shows that the friction at the surface enhances.
Taking into account this brief literature review, it can be found that investigation of the combined convective flow of molybdenum silica/water hybrid nanofluid about a rough sphere has not been studied earlier. The novelty of the present analysis includes the following points: -combined convective flow over a rough sphere; -impact of molybdenum-silica/water hybrid nanofluid on flow structures and heat patterns; -effect of surface roughness on hydro-and thermodynamics.
The modeled governing equations and their boundary constraints are affected to the appropriate non-similar reduction. Further, the obtained relations are attempted for solution transformations by employing the method of Quasilinearization and the implicit finite difference approximation [44,45].

Mathematical Analysis
Herein, we investigate the viscous, laminar, incompressible combined convection motion around a sphere. Here R is the sphere size. Let Tw and T∞ be the temperature at the border and away from the surface. The geometry and the configuration of the flow system are schematically presented in Figure 1. Here x is the curvilinear coordinate over the sphere border, while coordinate y is perpendicular it. The section r(x) explains the contour of the considered object.   The surface roughness is modeled in terms of large frequency and low amplitude sinusoidal waves. The control equations representing the hybrid nanosuspension flow variations are [2,7,[42][43][44][45]: Prescribed border constraints are Here, u indicates the nanosuspension velocity component over x coordinate, v indicates the nanosuspension velocity component over y axis and hybrid nanofluid temperature is indicated T. Further, ρ hnf , µ hnf , (ρc p ) hnf , (ρβ) hnf and k hnf denote the density, viscosity, thermal capacitance, heat expansion coefficient and heat conductivity of the hybrid nanosuspension. Following Khan et al. [36], Devi and Devi [46], Rostami et al. [5] and Waini et al. [7], the chemophysical attributes for mono and hybrid nanosuspensions are provided in Table 1. Table 1. Chemo-physical attributes of hybrid nanosuspension [2,5,36,46].

Properties
Water Meanwhile, Table 3 shows the different shape impacts on governing parameters. Further, the solid components of SiO 2 and MoS 2 are denoted by the subscripts s 1 and s 2 , respectively. The subscripts hnf, f and nf and denotes hybrid nanoliquid, base liquid and nanoliquid respectively. Following non-similarity parameters as in Patil et al. [44,47], Rajakumar et al. [32] and Roy et al. [33] are employed. The non-similarity transformations and resulting flow variables are as follows: In view of Equation (5), we have and Now, utilizing the transformations (5), Equations (1)-(3) are reduced as below and Equation (1) is satisfied identically.
Symmetry 2021, 13, 236 The employed border conditions are given below: The non-dimensional characteristics arising in this research are: The velocity profile in the mainstream current is given below: Therefore, ξ and γ(ξ) can be expressed in x as below: The above equations, that are expressed in terms of ξ, are rewritten in x, using the following relation: In view of Equations (11) and (12), Equations (6) and (7) become as follows The relevant border constraints are All parameters and non-dimensional characteristics that arise in the analysis are described in Nomenclature. The gradients at the sphere are: surface drag parameter and energy transport strength

Method of Analysis
Equations (13) and (14) are linearized employing the Quasilinearization method in the form: Here the parameters at the (i + 1)th iteration are found by using the known ith iterative coefficients. The non-dimensional border conditions become where boundary layer edge is denoted by η ∞ . Parameters in Equations (18) and (19) are Numerical solutions to the Equations (18) and (19) with the border constraints (20) were worked out through implicit finite difference approximation, where the backward and central differences are employed. The obtained algebraic set is solved employing the Varga's algorithm [48]. The step size ∆x = 0.01 is chosen along x direction and ∆η = 0.01 is chosen along η direction. When the present and previous iterations difference is low than 0.0001, the iteration process is stopped.

Results and Discussion
Molybdenum disulphide and silica/water combination of nanoparticles is taken for the present investigation. In the present model, 0.1 volume fraction of silica (ϕ 1 ) is included to the host liquid (H 2 O) and accordingly molybdenum disulphide (ϕ 2 ) is introduced with numerous concentrations to obtain MoS 2 -SiO 2 /water hybrid nanosuspension. By considering the different combination of nanoparticles, the similar approach has been investigated by many researchers [46,49,50]. The surface roughness is modeled in terms of large frequency and low amplitude sinusoidal waves. The frequency characteristic (n) and low variable (α) are employed for 5 ≤ n ≤ 100 and 0 ≤ α ≤ 0.3. The smooth border is denoted using α = 0, whilst the rough border is denoted using α = 0 [44,51]. As curvilinear x and y coordinates are perpendicular one another, the angle between these coordinates is taken as π/2. The magnitude ε = 0 characterizes the no-slip condition at the border, i.e., smooth border. The empirical shape factor s for nanoparticles is calculated using the formula s = ψ/3, where ψ is the sphericity. As per the Hamilton and Crosser approach [52], sphericity is the proportion of sphere area to the border area of physical particles with the same volumes. In the numerical computation, we have considered Pr = 7 that demonstrates the water as the host liquid, and it is constant throughout the investigation [30,53]. Figure 2 exhibits the effect of combined convection characteristic (Ri) and velocity ratio characteristic (ε) on velocity pattern {F(x, η)}. Rising values of combined parameter boost the nanosuspension velocity. Enhancing magnitudes of combined convection parameter specify higher significant impact of forced convection. Thus, the velocity of the fluid augments with the forced convection. In addition, in existence of rough surface (α = 0.1), velocity slip is noticed close to the surface. That is, abrupt jump arises close to sphere with roughness. Moreover, the characteristic (ε) is determined as the ratio between the slip velocity and the ambient circulation velocity. Thus, as the velocity ratio characteristic (ε) augments, the velocity of the liquid is enhanced. the wall temperature and ambient flow temperature increases, for increasing values of Ri, which reduces the fluid temperature. Further, for escalating magnitude of volume fraction (φ2), the temperature of the liquid augments. The liquid heat conductivity enhances for enhancing magnitude of volume fraction (φ2) and rises the capacity of the hybrid nanofluid to conduct more heat, which enhances the temperature of the liquid.       locity of the liquid diminishes. It is clear that the velocity of the fluid is larger in absence of nanoparticles compared its presence, for each corresponding value of ε. For increasing values of φ2 the hybrid nanofluid becomes denser and concentrated. Moreover, velocity of the fluid enhances with enhancing values of velocity ratio characteristic. This is caused because the larger value of velocity ratio characteristic indicates more slip at the border that raises the liquid velocity.    Figure 5 presents the impact of small parameter (α) and volume fraction (ϕ 2 ) on temperature profile {G(x, η)}. Temperature of the liquid diminishes with enhance in this small parameter. As we seen that the fluid temperature diminishes more in rough surface (α = 0) case compared to smooth surface (α = 0). Such phenomenon can be explained by the packets of fluid that are captured in the troughs formed owing to the surface roughness. It experiences slip at the wall of the sphere, when the liquid flows about the rough sphere and this enhances the velocity of the liquid. The growth of the liquid velocity yields more fluid movement across the wall of the sphere and this diminishes the liquid's temperature. In addition, the liquid's temperature enhances with rise in the volume fraction (ϕ 2 ). It is very clear that hybrid nanofluid has greater thermal conductivity than that of base fluid. Physically, it means that the nanoparticles with high thermal conductivity allow intensifying the heat transference and an inclusion of such nanoparticles enhances energy transport with a rise of the temperature. more fluid movement across the wall of the sphere and this diminishes the liquid's temperature. In addition, the liquid's temperature enhances with rise in the volume fraction (φ2). It is very clear that hybrid nanofluid has greater thermal conductivity than that of base fluid. Physically, it means that the nanoparticles with high thermal conductivity allow intensifying the heat transference and an inclusion of such nanoparticles enhances energy transport with a rise of the temperature. In addition, the drop between the border temperature and free stream temperature enhances for rising values of Ri, which increases the energy transport strength from wall to fluid. Further, for enhancing values of volume fraction (φ2), the friction at the surface and the energy transport strength enhance. This is because addition of nanoparticles in the host liquid leads to intensifying its temperature and as a result one can find an intensifi- for rising values of Ri, which increases the energy transport strength from wall to fluid. Further, for enhancing values of volume fraction (ϕ 2 ), the friction at the surface and the energy transport strength enhance. This is because addition of nanoparticles in the host liquid leads to intensifying its temperature and as a result one can find an intensification of the energy transport.
values of combined convection characteristic. The positive magnitudes of Ri illustrate essential impact of free convection over the forced one. Thus, the buoyancy force ins the liquid intensifies the liquid motion and as a result it augments frictions at the surf In addition, the drop between the border temperature and free stream temperature hances for rising values of Ri, which increases the energy transport strength from wa fluid. Further, for enhancing values of volume fraction (φ2), the friction at the surface the energy transport strength enhance. This is because addition of nanoparticles in host liquid leads to intensifying its temperature and as a result one can find an inten cation of the energy transport.      Figure 8 portrays the influence of small parameter (α) and volume fraction (ϕ 2 ) on energy transport strength Re −1/2 Nu . The energy transport strength boosts with enhancing nanoparticle volume fraction (ϕ 2 ). The physical reason is that these nanoparticles which are characterized by large heat conductivity display the augmentation of the energy transport coefficient because this parameter is directly proportional to heat conductivity. Further, heat transfer rate enhances with enhance in the small parameter. The sinusoidal waves are observed in presence of surface roughness. This is because the sphere has more area owing to the presence of roughness and thereby intensifies the heat transfer in a larger quantity from the sphere into the ambient fluid. transport coefficient because this parameter is directly proportional to hea Further, heat transfer rate enhances with enhance in the small parameter. T waves are observed in presence of surface roughness. This is because the sp area owing to the presence of roughness and thereby intensifies the hea larger quantity from the sphere into the ambient fluid.  for the surface drag coefficient Re 1/2 C f and energy transport strength Re −1/2 Nu . Energy transport strength and coefficient of friction at the border are enhanced for hybrid nanosuspension compared to mono nanofluid and regular liquid. In addition, the friction and rate of heat transfer high for rough surface (α = 0) case compared to smooth surface (α = 0) and sinusoidal waves are seen in presence of rough surface (α = 0). Hence, we can conclude that the hybrid nanofluid has more thermal conductivity than that of nanofluid and base fluid due to the rising magnitudes of nanoparticles concentration fundamentally designate the quicker heat transmission from surface of the cylinder to ambient fluid, this boosts the thermal conductivity and consequently temperature of hybrid nanofluid, which causes the intensive enhance in the magnitudes of temperature profile and rate of heat transfer coefficient, significantly.
Symmetry 2021, 13, x FOR PEER REVIEW conclude that the hybrid nanofluid has more thermal conductivity than that of na and base fluid due to the rising magnitudes of nanoparticles concentration fundam designate the quicker heat transmission from surface of the cylinder to ambient fl boosts the thermal conductivity and consequently temperature of hybrid na which causes the intensive enhance in the magnitudes of temperature profile and heat transfer coefficient, significantly.      Their impact can be observed in the chemo-physical attributes of the hybrid nano sion. Among the considered shapes, blade-shaped particles give the high heat mance compared to the platelets, cylindrical and bricks, and sphere-shaped partic the low heat performance. Further, enhancing values of small parameter, the temp profile lessens, while Nusselt number augments. Physically, addition of nanopar the regular liquid leads to the temperature growth, so an intensification of the number is quite obvious. For instance, at 0.5 x = the heat transfer rate enhance 16%, from spherical to blade shaped nanoparticles type in presence of rough surf   Taking into account an influence of the considered shapes of nanoparticles on velocity profiles (see Figure 15) it is possible to highlight that the blade-shaped particles give the highest growth of the velocity in comparison with other studied shapes, and sphericalshaped particles illustrate the formation of the smallest velocity. Further, a rise of small parameter reflects a growth of the velocity profile. Figure 16 explains the impact of non-similarity parameter and combined convection parameter on velocity distribution Velocity profile reduces while, temperature profile boosts in absence of non-similarity variable in opposing buoyancy flow case. Further, opposite outcomes are observed in presence of non-similarity parameter in assisting buoyancy motion case. Reason is that x directs as a favorable pressure drop. Intensifying magnitudes of combined convection coefficient augments the velocity distribution while, reduces the temperature distribution. It is clear from figure that the velocity of liquid is more pronounced for larger Ri in assisting motion variant. Physically, more buoyancy force is generated for greater values of Ri, which gave the maximum energy to the moving fluid, and such heat creates the confrontation in the circulation. The opposing model was scrutinized for the velocity in opposing flow version. In addition, it ascertains that the temperature and thermal boundary layer are diminished owing to Ri in the considered circulations. Taking into account an influence of the considered shapes of nanoparticles on velocity profiles (see Figure 15) it is possible to highlight that the blade-shaped particles give the highest growth of the velocity in comparison with other studied shapes, and sphericalshaped particles illustrate the formation of the smallest velocity. Further, a rise of small parameter reflects a growth of the velocity profile.   Figure 16 explains the impact of non-similarity parameter and combined convection parameter on velocity distribution {F(x, η)} and temperature distribution {G(x, η)}. Velocity profile reduces while, temperature profile boosts in absence of non-similarity variable in opposing buoyancy flow case. Further, opposite outcomes are observed in presence of non-similarity parameter in assisting buoyancy motion case. Reason is that x directs as a favorable pressure drop. Intensifying magnitudes of combined convection coefficient augments the velocity distribution while, reduces the temperature distribution. It is clear from figure that the velocity of liquid is more pronounced for larger Ri in assisting motion variant. Physically, more buoyancy force is generated for greater values of Ri, which gave the maximum energy to the moving fluid, and such heat creates the confrontation in the circulation. The opposing model was scrutinized for the velocity in opposing flow version. In addition, it ascertains that the temperature and thermal boundary layer are diminished owing to Ri in the considered circulations.  To verify the accuracy of the computational approach, the obtained outcomes found for several values of Ri, by setting all characteristics values to zero except Pr = 0.7 are compared with the outcomes of Rajakumar et al. [32] in Figure 17 and it is obtained that the outcomes are well agreed. To verify the accuracy of the computational approach, the obtained outcomes found for several values of Ri, by setting all characteristics values to zero except Pr = 0.7 are compared with the outcomes of Rajakumar et al. [32] in Figure 17 and it is obtained that the outcomes are well agreed.  Table 4 shows the influence of different values of combined convection parameter for mono nanofluid and hybrid nanofluid for skin-friction coefficient and energy transport strength. A rise of Ri reflects an increase in heat transport strength and friction at the surface for both the mono nanofluid and hybrid nanofluid. In addition, it is clear that, the hybrid nanofluid values are comparatively larger than that of nanofluid, for skin-friction coefficient as well as energy transport strength in smooth surface case. In addition, it is seen that the same trend is observed in rough surface case. Table 4. Comparison of hybrid nanofluid (φ1 = 0.1 and φ2 = 0.1) and mono nanofluid (φ1 = 0.1 and φ2 = 0.0) for different magnitudes of combined convection coefficient (Ri) at smooth and rough surfaces, for skin-friction coefficient and energy transport strength when 0.5 x = , n = 10, ε = 0.1.

Combined
Mono Nanofluid Hybrid Nanofluid  Table 4 shows the influence of different values of combined convection parameter for mono nanofluid and hybrid nanofluid for skin-friction coefficient and energy transport strength. A rise of Ri reflects an increase in heat transport strength and friction at the surface for both the mono nanofluid and hybrid nanofluid. In addition, it is clear that, the hybrid nanofluid values are comparatively larger than that of nanofluid, for skin-friction coefficient as well as energy transport strength in smooth surface case. In addition, it is seen that the same trend is observed in rough surface case. Table 4. Comparison of hybrid nanofluid (ϕ 1 = 0.1 and ϕ 2 = 0.1) and mono nanofluid (ϕ 1 = 0.1 and ϕ 2 = 0.0) for different magnitudes of combined convection coefficient (Ri) at smooth and rough surfaces, for skin-friction coefficient and energy transport strength when x = 0.5, n = 10, ε = 0.1.

Combined Convection Parameter (Ri)
Mono Nanofluid (ϕ 1 = 0.1 and ϕ 2 = 0.0) Hybrid Nanofluid (ϕ 1 = 0.1 and ϕ 2 = 0.1)  Table 5 demonstrates the comparison with previous published data illustrating a very good agreement. It presents the values of heat transfer rate against x and Ri. The negative magnitudes of mixed convection coefficient (Ri < 0) indicate the opposing buoyant flow, while Ri > 0 is for assisting buoyant flow. From Table 5, it is clear that Nusselt number is a reducing relation of x.

Conclusions
The current investigation explores the combined convective flow of silica-molybdenum disulphide/water nanosuspension over a rough sphere. The influences of varying small parameter, mixed convection, concentration of nanoadditives and frequency characteristic on the non-dimensional velocity and temperature distributions, surface drag coefficient and energy transport strength are investigated and explained with the help of different graphs. From the detailed analysis of the obtained results, the following vital observations are made: the hybrid nanofluid augments the temperature, as well as Nusselt number at the wall; -velocity distribution is reduced, while friction at the surface enhances in the case of hybrid nanofluid; -slip exists near the surface of the sphere owing to roughness, which yields steep jump in the liquid velocity close to the border; -for enhancing magnitudes of small parameter, the friction at the surface and the Nusselt number are enhanced; -skin-friction parameter and the Nusselt number are boosted for hybrid nanoliquid on comparison with mono nanofluid and the base fluid; -skin-friction parameter and the Nusselt number have the highest values for bladeshaped nanoadditives and lowest value is found for spherical-formed nanoparticles. Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available in this paper.

Conflicts of Interest:
The authors declare no conflict of interest.