Transportation of Hybrid MoS2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\textbf {MoS}}_2$$\end{document}–SiO2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\textbf {SiO}}_2$$\end{document}/EG Nanofluidic System Toward Radially Stretched Surface

A stochastic computing approach is implemented in the present work to solve the nonlinear nanofluidics system that occurs in the model of atomic physics. The process converts the partial differential nanofluidics system with suitable level of similarities transformation into nonlinear systems of differential equations. For the construction of datasets, finite difference scheme (Lobatto IIIA) is applied through different selection of collocation points for nonlinear nanofluidics system having accuracy of order four. Lobatto IIIA has a strong point to tackle extremely nonlinear systems of ordinary differential equations in smooth way. For different scenarios, datasets are well trained through computing scheme to investigate the heat transfer and thermal performance of nanofluidic transportation system of nanofluids and hybrid nanofluids toward stretching surfaces with variation of Biot number, Nusselt number and skin fraction. Furthermore, the reliability, accuracy and efficiency are endorsed through various statistical analysis and graphical illustrations of proposed computing scheme.


Introduction
Several base fluids like water, glycol and alcohol have rather low viscosity, thermal conductivity and may not be used independently in several applications, the suspension containing nanoparticles combined with the base fluid [1]. A special type of fluid known as nanofluid has been introduced, to overcome this weakness. Due to the extreme essential properties of heat transfer, such nanofluids are commonly used for macroscopic cooling [2]. In everyday life, implementations such as portable computer processors, lubricant processing, coolants, air conditioning, freezers and nanotechnology, the single-phase flow model of nanofluids is of considerable significance [3]. Narla et al. have presented the pumping of electroosmotic nanofluid past a curved channel [4]. A very simple and environment friendly conversion method of ultraviolet irradiance to heat energy is identified as a solar collector, used in the industrial sector [5]. In such collectors, different fluids like water and glycol are being used to improve their power and performance [7]. Sheikholeslami et al. have analyzed the hydrothermal analysis of nanofluids with wavy pipe absorber [6]. Recently, both numerically and experimentally, Li et al. [8] have presented the use of nanofluids in certain kinds of solar collectors. Oke [9] has discussed and numerically analyzed a special type of nanofluid known as gold-water nanofluid and presented the obtained results.
There is another classification of nanofluids obtained by dispersing two or more types of nanoparticles in conventional fluid, termed as hybrid nanofluids [10]. Hybrid nanofluids are essential because they are potential heat transfer fluids that possess more enhanced thermophysical properties and better thermal performance than base fluids and nanofluids having single type of nanoparticles [11]. Algehyne et al. [12] have presented the BVP of Maxwell hybrid nanofluid and its thermophysical performances. Khan et al. [13] have investigated a fractional model on hybrid nanofluids with porous medium. The use of nanoparticles for imaging and sensing purposes in bio-medical field has become very important for advance applications [12]. This growth in the field of nanoparticles is directly linked with designing new materials at nanoscale level along with innovations in imaging technologies to measure and manipulate nanomaterials [14]. The ternary MHD hybrid nanofluid flow past a rotating surface is treated numerically by Oke [15]. The production, usage and disposal of synthesized nanoparticles will cause expulsion of air, sand, soil and water systems. Harmful effects are to be expected; therefore, minimization and quantification of these negative effects on ecological health are of great importance [16]. Therefore, true information about chemical and physical properties of nanoparticles in consideration of realistic circumstances is necessary to predict their toxicity, fate and behavior in normal aquatic environment [17].
The phrase thin film is widely used to describe flow that is defined by the thickness of a layer of material ranging from nanometer scale to the several micrometers. In one dimension (one or two dimensions), the flow spectra are significantly less than the other [18]. It is a property that may be utilized to simplify Navier-Stokes equations [19]. The structural importance of flow problems with thin films (flow) is significant [20]. Thin-film flow problems in regular activities range from simple to sophisticated effects, such as transmitting a rain drop on an airplane window [21]. Siddiqui et al. [22] implemented the He's homotopy method for third-grade TFF fluid flow. It is true that thin-film flow is used in a variety of industrial processes [23]. The Casson hybrid nanofluid based on TFF is analyzed by Alhussain and Tassaddiq [24]. Researchers are interested in working on this topic because of its extensive applications in sectors such as coating material, industrial processes, chemical engineering, cooling and the lubrication of thermal changer fins into contact lens motion, dam breaking wave simulation and nuclear reactor fluid dynamics that are among the most significant industrial applications [25]. The TFF flow past a stretching sheet is numerically tackled by Kamis et al. [26]. The representation of TFF dynamics is similar to tear films in the eye and membranes in the biophysics field [27]. The most prevalent fluid properties are thin-film flow and flow characteristics. Furthermore, in this advanced and fast-paced era, thin film allows all those involved to monitor the output of paints in surface flow coverings, plan the layout of nuclear reactors and reduce the structure of equipment such as lubricating fluid [28]. TFF's biological uses are similar to red blood cells flowing via veins, eyes and lungs [18]. A neuroheuristic intelligent algorithm is utilized by Ahmad et al. [29] to evaluate TFF model.
Many engineering systems involving electro-magnetic accelerators and magnetohydrodynamic (MHD) power generators work on high temperatures; therefore, the understanding of MHD nanofluidic flow past has become very essential in production of appropriate equipment [31]. Thus, problems regarding nanofluids and magnetohydrodynamics (MHD) have become mechanically more important [32]. The modified Eyring Powell MHD fluid flow is presented and numerically solved by Oke [33]. The concept of using nanofluids in magnetohydrodynamic (MHD) flows over different flow geometries (specifically stretching sheets) has been presented due to the stability, enhanced heat transfer and effective thermal conductivity of these kinds of fluids [34]. Thermal conductivity of the nanofluids contributes in better performance, high energy proficiency and subordinate operating costs, discussed by Ijaz and Ayub in [35].
Lobatto IIIA based on finite difference method is a version of bvp4c method. The numerical code bvp4c in MATLAB, which is a finite difference technique, implements the 3-stage Lobatto IIIA formula that is a collocation technique having accuracy of order four. In the literature usually, the symbol III is associated with Lobatto methods because symbols I and II are held in reserve for two kinds of Radau methods [36]. A collocation method based on Lobatto IIIA technique is utilized by Oke et al. [37] for numerical solution of transformed ordinary differential equations. Lobatto IIIA method has a strong point to tackle extremely nonlinear systems of ordinary differential equations (ODEs) [38]. Implementation of Lobatto IIIA method to attain better approximation and rapid convergence is an inventive effort, used by Ahmad et al. [39].
The salient features of the proposed techniques are itemized as follows: • A novel investigation is presented to examine the heat transfer and thermal performance of nanofluidic transportation of nanofluids and hybrid nanofluids toward stretching surfaces by Lobatto IIIA techniques. • The governing physical flow systems are modeled mathematically in curvilinear coordinates and simplified by using appropriate transformations. • The simplified systems of nonlinear differential equations are numerically treated to investigate the influences of emerged physical parameters on temperature, pressure and velocity profiles. • A comparative dynamical analysis of nanofluids and hybrid nanofluids under the influence of emerging parameters is presented through tables and graphs which ascertain the worth of the scheme as an effective and viable computing solver. • The conceptual simplicity, smoothness and stability are the hallmarks of the proposed study.
The framework of rest of the manuscript is structured as follows: Mathematical development and modeling of the problems is provided in Section two, results obtained by computational technique are presented graphically as well as in tabulated form in Section three, while the Section four includes the conclusion and future research study.

Mathematical Development
This section is held in reserve to the mathematical development and numerical simulation of axisymmetric flow of SiO 2 /EG-based nanofluid and hybrid MoS 2 -SiO 2 /EGbased nanofluid on a radially stretched surface under consideration of magnetic field effect along with the internal heat absorption/generation. Lobatto IIIA is employed for computational exploration of governing flow model. Here ethylene glycol is taken as base fluid, molybdenum sulfide (MoS 2 ) and silicon dioxide (SiO 2 ) as ultra-fine nanoparticles. Hybrid nanofluid is prepared by inserting the nanoparticles of MoS 2 and SiO 2 in ethylene glycol.

Mathematical Modeling
A porous infinite radially stretched disk is considered for electrically conducting and incompressible nanofluidic flow. This opted disk is placed at the plane wherez = 0. Nanocomposite MoS 2 -SiO 2 /EG fluid flow is taking place above a half spacez > 0. This hybrid nanofluidic mixture contains 1% of SiO 2 volumetric fraction (φ 1 ), while MoS 2 volumetric fraction (φ 2 ) varies from 1 up to 5%. Equation (1) determines the nanoparticles volumetric fraction φ hnf of the hybrid MoS 2 - Magnetic field B 0 having uniform intensity is applied in axial direction, and magnetic induction is unnoticed assuming the sufficiently weak magnetic field. The current analysis is done in cylindrical coordinates (r ,θ,z). Moreover, the flow is generated due to the reason of stretching of the surface along radial direction. This flow has uniform velocity U w (r ) = cr , where real number c is nonnegative. Surface of the porous disk is incorporated with heat transfer along with internal heat generation. Furthermore, lower surface is assumed to be heated through convection using a warm fluid having temperature T f . Temperature T f is more than that of constant ambient temperature T ∞ and gives a heat-transfer coefficient h f . The physical structure of governing flow model is demonstrated in Fig. 1.
The derivatives regardingθ are absent because of axial symmetry of the leading engineering problem. Using above considerations together with boundary-layer approximations, the principal equations can be coupled into the subsequent forms.
along with related boundary conditions where u and w are components of velocity alongr andz-axis, correspondingly. Furthermore, μ hnf , ρ hnf and σ hnf indicate hybrid nanofluidic dynamic viscosity, density and electrical conductivity, respectively. f w (r ) is injection/suction velocity over the radial disk's surface.

Numerical Simulation of Model
Numerical exploration of temperature and velocity profiles for the hybrid MoS 2 -SiO 2 /EG nanofluidic flow over a porous infinite radially stretched surface is done by employing Lobatto IIIA technique according to process illuminated  in Fig. 2. Moreover, thermo-physical properties of ethylene glycol, SiO 2 and MoS 2 are listed in Table 1. Similarity transformations as well as dimensionless variable ζ and f (ζ ) are accustomed to more assist and simplify the current investigation.
After simplification, system of PDEs Eqs. (2)(3)(4)(5) for proposed physical problem is reduced into equivalent transformed system of nonlinear coupled ODEs Eqs. (7-9) in which Local Nusselt number Nur and skin friction coefficient C f are concerned physical quantities, which are defined as where q w and τ w are In sight of Eq. (6), expression, which is labeled in Eqs. (14) and (15), gives the dimensionless local Nusselt Nur and skin friction C f as where Nur and C f are reduced Nusselt number and skin friction, respectively. The system of nonlinear coupled ODEs Eqs. (7-9) is then converted into system of first-order ODEs for Lobatto IIIA. In this investigation, mesh points are fixed at 900, while relative error tolerance varies as 10 −3 , 10 −6 , 10 −9 and 10 −12 . Approximated velocity and temperature profiles of the hybrid nanofluidic system are determined under consideration of dimensionless quantities, e.g., Hartman number M, heat generation parameter λ, Prandtl number Pr and Biot number Bi-based scenarios. At the end stability, accuracy and convergence analysis are performed.

Results and Discussion
In proceeding section, numerically computed impacts for evolving physical constraints on nanofluidic motion and temperature of SiO 2 /EG nanofluid and hybridized MoS 2 -SiO 2 /EG nanofluid are explained. Hartman number M, suction parameter S, nanoparticles volumetric fractions φ 1 and φ 2 , heat absorption/generation parameter λ, Biot number Bi and shape factors m are several influential constraints of interest in this exploration. Graphical and numerical outcomes are designed for both categories of nanofluids to provide comparative study. Dotted curves illustrate outcomes achieved for SiO 2 /EG nanofluid, whereas solid curves explain outcomes for hybridized MoS 2 -SiO 2 /EG nanofluid. Six scenarios each with three cases of the evolving constraints for both nanofluids are shown in Tables 2, 3, 4, 5, 6, 7 and 8. In these tables different values are given to sundry variables in different cases. We have 6 scenarios and 3 cases for each scenario, which demonstrates variation of each parameter. Numerical calculations for different tolerance rates like (10 −3 , 10 −6 , 10 −9 and 10 −12 ) show the adequate accuracy and better convergence obtained by computing scheme. Mesh points used for solutions of each case of 6 scenarios are revealed in Table 2. Table 3 indicates stopping criteria depending on maximum residual obtained for different cases of each scenario, during problem evaluation. Complexity analysis of numerical method Lobatto IIIA to obtain required values of residual error is performed in terms of number of evaluated ODEs and BCs. Number of evaluated ODEs to develop the desired residual is given in Table 4 for different scenarios and cases. The evaluated boundary conditions for the specified accuracy given in Table 3 during different cases of each scenario are presented in Table 5.  Tables 7  and 8.
The impacts of suction and magnetic constraints on fluidic motion are displayed in Fig. 3. From Fig. 3a it is determined that suction parameter decelerates flow past of both nanofluids and works as an opposing force. This decreasing trend is caused by the suction taking place at porous stretched surface that interrupts the boundary layer fluid flow and leads to descend fluidic velocity. Similarly rising values of magnetic parameter reduce rate of fluid's flow as can be perceived in Fig. 3b. Existence of magnetic field in directed flows is significant. Moreover, perpendicularly applied magnetic field to the nanofluidic flow direction has dominant effects in opposing manner. Therefore, parameter M acts as a drag force and retards the movement of nanofluids. Deceleration of flow due to suction and magnetic constraints is observed leading for SiO 2 /EG nanofluid. Furthermore, in Table-6 skin friction coefficient which is also termed as velocity gradient rises using positively increasing values of both suction and magnetic constraints.
Biot number (sheet convection) gives influence on nanofluidic velocity and f (ζ ) which is presented in Fig. 4. Nanofluidic flow under up surging values of Biot number possesses ascending behavior, while graphs for f (ζ ) exhibit descending behavior as demonstrated in Fig. 4a, b, respectively. Moreover, leading rise is appeared for hybridized MoS 2 -SiO 2 /EG nanofluid, while leading drop is found for SiO 2 /EG nanofluid. Figure 5 reveals the effects for volumetric concentrations of SiO 2 and MoS 2 -SiO 2 in ethylene glycol on nanofluidic velocity and f (ζ ). It is noticed that a marginal rise in volumes of nanoparticles accelerates the flow and raises the graphs for f (ζ ) as given in Fig. 5a, b, respectively. Since usage of nanoparticles in conventional fluid decreases its viscosity which helps the fluid flow to accelerate. Faster increase in velocity and f (ζ ) is observed for hybridized nanofluid. Furthermore, greater values of φ 1 and φ 2 increase the skin friction coefficient which is more leading for MoS 2 -SiO 2 /EG as given in Table 6. Figures 6, 7 and 8 are plotted to explain the influence of some evolving constraints such as volumetric fractions, Biot number and heat absorption/generation parameter with each against five shape factors on profiles of temperature for both fluids. Sphere, hexahedron, tetrahedron, column and lamina are five considered shapes of nanoparticles. Figure 6 is devoted to inspect the thermal performance of nano and hybrid nanofluid by submerging nanoparticles of both categories up to 5% of volume. Figure 6a, b demonstrates that increasing volume of nanoparticles raises the temperature of both fluids. For hybrid nanofluid 1% volume of        SiO 2 nanoparticles is reserved fix. Moreover, both nanofluids holding lamina shape of nanoparticles offer maximum temperature for variation of φ 1 and φ 2 . On comparing Fig. 6a, b, it is also visualized that the temperature of hybrid nanofluid leads the temperature of nanofluid. Thermal behavior of nano and hybrid nanofluid subject to rising sheet convection constraint is pictured in Fig. 7. Biot number (sheet convection) is defined as a quotient of heat transmission resistances within and above the conductive surface that measures the temperature difference by means of applied thermal gradient over the conductive surface. Therefore, upsurge in this quotient raises the temperature curves for both SiO 2 /EG and MoS 2 -SiO 2 /EG nanofluids. Moreover, high temperature is measured by submerging nanoparticles with lamina shape into the both fluid categories. Tables 7 and 8 validate that Nusselt numbers for both fluids are increased using greater values of Bi while decreased using greater values of shape factors.
Impressions of heat absorption/generation constraint on temperature curves are provided in Fig. 8. Since λ < 0 is heat absorption case, so descending λ values decrease the temperature graphs as well. This behavior is perceived for both SiO 2 /EG and MoS 2 -SiO 2 /EG nanofluids. Lowest temperature is achieved for spherical nanoparticles. In Tables 7  and 8, it is witnessed that Nusselt numbers for both fluids are increased using smaller values of λ while decreased using greater values of shape factors.

Conclusion
In this work, a novel approach is presented to numerically examine the nanofluidic and hybrid nanofluidic sys-tems MoS 2 -SiO 2 /EG, through designed computing scheme. In the proposed model, axisymmetric flow of SiO 2 /EG nanofluid as well as MoS 2 -SiO 2 /EG hybrid nanofluid over a radially stretched disk under influential impacts of magnetic field, suction, internal heat absorption/generation, sheet convection and volumetric fractions is investigated in detail. The foremost judgments of presented work are: • It is determined that MoS 2 -SiO 2 /EG hybrid nanofluid holds better thermal abilities because its temperature upsurges more rapidly with rising Biot number and volume fractions as matched to SiO 2 /EG nanofluid. • High magnetic field and suction of porous disk decelerate flow past both nanofluids inside the boundary; on the other hand, volumetric fractions accelerate it. • Lowest nondimensional temperature is witnessed for nanoparticles of sphere shape, whereas peak temperature is perceived for lamina shape. • In case of heat absorption constraint, negative impacts on thermal performance of both nanofluids are observed. • Skin friction coefficient of both nanofluids gain rises, when subject to the influence of suction, magnetic and volumetric fraction constraints and leads in hybrid nanofluidic case. • Constraints of heat absorption and sheet convection raise the Nusselt number of both nanofluidic categories, but shape factors are a cause to decrease its values.
In future, one may exploit artificial intelligence-based numerical computing paradigm [40][41][42][43] for dynamical analysis of hybrid nanoparticles (MoS 2 -SiO 2 , TiO 2 -CuO) on the fluid toward stretched surfaces with the influential impacts of evolving constraints such as slip condition, nanoparticles  shape factors, heat generation, thermal radiation and applied transverse magnetic field.
Funding Open access funding provided by Università degli Studi di Palermo within the CRUI-CARE Agreement.

Conflict of Interest
The authors declare that there is no conflict of interest, in any aspect.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.