key: cord-1024364-pe2oxvql authors: Ahmad, Iftikhar; Cheema, Tahir Nawaz; Raja, Muhammad Asif Zahoor; Awan, Saeed Ehsan; Alias, Norma Binti; Iqbal, Sana; Shoaib, Muhammad title: A novel application of Lobatto IIIA solver for numerical treatment of mixed convection nanofluidic model date: 2021-02-24 journal: Sci Rep DOI: 10.1038/s41598-021-83990-8 sha: 3572852212be5a615e292557aab318ce03293c1d doc_id: 1024364 cord_uid: pe2oxvql The objective of the current investigation is to examine the influence of variable viscosity and transverse magnetic field on mixed convection fluid model through stretching sheet based on copper and silver nanoparticles by exploiting the strength of numerical computing via Lobatto IIIA solver. The nonlinear partial differential equations are changed into ordinary differential equations by means of similarity transformations procedure. A renewed finite difference based Lobatto IIIA method is incorporated to solve the fluidic system numerically. Vogel's model is considered to observe the influence of variable viscosity and applied oblique magnetic field with mixed convection along with temperature dependent viscosity. Graphical and numerical illustrations are presented to visualize the behavior of different sundry parameters of interest on velocity and temperature. Outcomes reflect that volumetric fraction of nanoparticles causes to increase the thermal conductivity of the fluid and the temperature enhances due to blade type copper nanoparticles. The convergence analysis on the accuracy to solve the problem is investigated viably though the residual errors with different tolerances to prove the worth of the solver. The temperature of the fluid accelerates due the blade type nanoparticles of copper and skin friction coefficient is reduced due to enhancement of Grashof Number. The suspension of nanoparticles into the base fluid to increase the thermal conductivity and rate of heat transfer is termed as nanofluid. The term mixed convection is used for the transfer of heat. In the last few decades, it is observed that nanofluids are vastly used to increase the thermal conductivity of the traditional fluids i.e. water, machine oil, fuel oil etc. The uses of nanofluids are in nanotechnology, automatic cooling, temperature exchanger and targeted medication transfer etc. Innovative work by Choi and Eastman conclude that the thermal conductivity of the base fluid can be improved by adding the small number of nanoparticles 1 . Metals, oxides, nitrates and carbides are mostly used nanoparticles. Only 1-5% amount of nanoparticles is added to the base fluid to improve the thermal conductivity of the base fluid. After these preliminary considerations, many investigations have been performed for the nanofluidic problems along with transfer of heat. Advanced heat conduction is important for nanofluids due to large surface area of nanoparticles which permits for additional heat transfer 2 . The important flow phenomenon arises in the field of boundary layer flow problems along with its extending surface. Aerodynamics, plastic sheet extrusion and artificial extrusion are some of the applications for different fluids. Initially, Sakiadis 3 analyzed the boundary layer singularity on a solid plate. Crane 4 observed the flow over a parallel stretching sheet. Later, many researchers worked on the flow for extended as well as moving surfaces under different situations. The flow on a stagnation point past a stretching surface was examined by Banks 5 . Khan and Pop 6 investigated the boundary layer for nanofluid under the impact of Brownian wave theory and concluded that the both phenomena have significant impact on heat flux and the skin friction at the surface. Das et al. 7 examined the thermal energy influence on the nanofluid. The researchers concentrated on the nanofluid transference by taking base fluid that has invariant fluidic features. The impact of viscosity on the transfer of heat and nanofluid flow was observed by mixing copper and silver nanoparticles that are antimicrobial agents and considered by field emission scanning electron microscope, transmission electron microscope and scanning electron microscope into the water. The experimental results are useful for reproduction filaments 8 . Elbashbeshy 9 investigated the results of transfer of heat analysis past a stretching sheet. In the recent articles, too much focus was given on the transport of nanofluid and in these articles the base fluid was considered to have invariant rheological characteristics 10 . The influence of variation in viscosity and micro variation on oblige transport of copper-water nanofluid is also examined. Lobatto IIIA method is considered for examining the stability of properties for boundary value problems (BVP). This method is named after Rehuel Lobatto and it is used for the numerical integration of differential equations. Lobatto IIIA method is applied for nonlinear couple system of differential equations arises in the fields of mechanics and electrical circuits 11 . Generally, industrialized liquids come across having a variable viscosity, which can also be observed as compression, shear or dependent on temperature. Zehra et al. 12 investigated the flow of non-Newtonian fluid over an inclined channel. The influence of viscosity on drift and warmth transfer with Ag-water and Cu-water nanofluids over a shifting surface is examined by Vajravelu 13 . In recent times, Tabassum et al. 14 observed the impact temperature is enhanced by increasing the variable viscosity for the nanofluidic system under. Copper steel and silver nanoparticles are frequently used in electronic equipment and thermal conduction 15 . Main functions of copper are in the electrical cable, tiling and pipes used in manufacturing. Numerical work for the flow of Darcy Forchhieimer was studied for the analysis of the sisko nanoparticle by compelling nonlinear thermal radiation with Lobatto IIIA method 16 . The study of Lobatto IIIA was investigated on the porous medium for the flow of two-dimensional flow of magneto hydro dynamic fluid taking the conditions of Navier's slip and activation energy 17 . Nanofluid flow pass on a curved surface using Joule heating was considered to study the influence of copper and silver nanoparticles. The nonlinearity in thermal radiation on the curved elongating sheet was examined 18 . The effect of silver and copper nanomaterials on an inclined stretching sheet was investigated to study the effect of viscous fluid with mixed convection on the sheet. Resultant ordinary differential equations are numerically solved by a software Mathematica with Euler's Explicit Method (EEM) 19 . The effort for the enlargement of heat transfer was analyzed for enlightening the performance of electrical devices by manipulating the porous media for enhancing the heat transfer rate by using different boundary conditions and structures 20 . Ghadikolaei et al. 21 investigated the mixed convection flow of hybrid nanoparticles of silver with the 50-50 percent quantity of ethylene-glycol water on the stretching surface in the presence of variable viscosity and different shapes of nanoparticles by a well-known numerical method Runge-Kutta Fehlberg fifth order (RKF-5) and also analyzed the flow of GO-MoS2 hybrid nanoparticles in H2O-(CH2O) hybrid base fluid under the effect of H 2 bond with (RKF-5) numerical technique 22 . Zangooee 23 examined the hydrothermal flow of magneto hydro dynamic flow of titanium dioxide with glycol nanofluid flow on the two radiative stretchable revolving disks using AGM. Similarly, many fluidic problems are investigated with numerical and analytical solver such as magnetic field effect on nanofluid flow using AGM and ADM [24] [25] [26] , heat transfer rate in nuclear waste 27 , fins arrangement in cubic enclosure 28 , heat transfer simulation in a channel with rectangular cylinder 29, 30 and 3D optimization of baffle arrangement 31 . Besides these, the research community has exploited different numerical schemes in diversified applications [32] [33] [34] [35] [36] [37] [38] [39] . The current article is an attempt to investigate the impact of viscosity in the transfer of heat and thermal conductivity by implementation of numerical scheme based on Lobatto IIIA method. The tabular and diagrammatical To analyze the effect of nanofluidic system under the impact of heat transfer along a stretching sheet, the rectangular coordinate structure is considered, where x-axis and y-axis are perpendicular and parallel to the stretching sheet in Cartesian coordinates system. The mixed convection with copper and silver in the water is constricted above the surface and stretching sheet is taken along y-axis. V w (x) = ax, where a is the constant, V w (x) is the velocity of a linearly stretching surface and the solid-state two-dimensional stagnation flow is assumed on the surface which is V e (x) = cy, where c is constant. Tabassum 14 worked on the influence of variation in the resistance. Figure 1 represent the physical geometry of the proposed problem. Moreover, we supposed the thickness of nanoparticle in the Vogel's model. In the direction of y-axis, a uniform magnetic field B o is applied in the normal direction of the stretching surface. The uniform magnetic field is induced due to magnetic Reynolds number which is taken very small that can be neglected as compared to magnetic field. Maraj 24 express the system of governing equations in the following form: where, T z and T ∞ are temperatures at the surface and away from the surface. The Boundary conditions are where u and v are components of velocity along x and y direction, respectively σ nf , k nf , ρ nf , µ nf and (cp) nf are electrical conductivity, thermal conductivity representing the rate at with heat permits through the fluid, density, dynamic viscosity and specific heat capacity of nanofluid. The properties of base fluid and nanoparticles is represented in Table 1 . Table 2 represents the thermophysical properties of σ nf , k nf , µ nf ,ρ nf . The variable viscosity measures the resistance occurs due to deformation of the fluid represented in the form of Vogel's model as The physical quantities alike rate of heat flux Nu x and shear strain Cf x and at the surface can be expressed as Density Thermal conductivity Heat capacity www.nature.com/scientificreports/ where, τ z represent the shear stress of wall and q z is heat flux at wall. The dimensionless form of above expression by invoking similarity transformation is To find the solution of coupled nonlinear Eqs. (7) (8) with the boundary conditions in Eq. (12) . Lobatto IIIA is incorporated with MATLAB routine 'bvp4c' . The strategy for the solution of problem is presented in Fig. 2 . The first step in the diagram provides the fluid flow system dynamics, the transformation system for PDEs to ODEs is given in the second step, the third step is consisting of two parts i.e., non-linear higher orders ODEs of the fluidic system and their conversion into first order system of ODEs. In the fourth step of the fluid flow diagram, velocity and temperature profiles are attained based on parameters of interest of the flow system, while the last step in the block diagram represents the accuracy, convergence, and stability analysis in terms of residual errors. The nonlinear ordinary differential equations are solved with Lobatto IIIA method by using MATLAB software. The numerical solution obtained for the distribution of temperature and velocity are observed and showed through diagrams. Figures 3, 4 , 5, 6, 7, 8, 9, 10, 11 and 12 shows the consequence of extending ratio parameter, www.nature.com/scientificreports/ Figure 3a ,b it is examined that when A < 1 the velocity profile is concave down for variation in ratio parameter A for both silver and copper particles with water based nanofluids and velocity profile is concave up for the increasing value of elongating ratio parameter. It is observed that the decreasing value of A represent the condition when stagnation velocity is less than elongating ratio and the increasing value of A represent the state when the stagnation velocity becomes greater than the ratio parameter A. Figure 4a ,b shows that the viscosity www.nature.com/scientificreports/ parameter decreases for A < 1 and goes on increasing when A > 1. It is observed that the velocity decreases with increasing value of viscosity parameter. In Fig. 5a ,b behavior of volume fraction for Cu and Ag is shown. It shows that the flow accelerates for increasing value of volume fraction when the value of ratio parameter is decreasing, and the flow behaves oppositely when the value of elongating ratio parameter is increasing. Figure 6a ,b shows the result of Grashof number and it is observed that the velocity increase for increasing the values of Grashof number. The velocity increases because Grashof number is proportion of resistance forces to the viscid forces. Figure 7a ,b represent the impact of magnetic field and it shows that the magnetic field constantly turns as a resistive energy so the velocity decreases for increasing values of magnetic field for both silver and copper based nanofluids. Figures 8, 9 , 10, 11 and 12 shows the effect of temperature for different values of parameters such as nanoparticle shape m, Grashof number Gr, magnetic field M, viscosity parameter and viscosity parameter ϕ for www.nature.com/scientificreports/ both the nanoparticles i.e. copper and silver. In Fig. 8a ,b it is observed that the presence of nanoparticle enhance the base fluid capacity to bearing heat and therefore, temperature droplets so that is why due to high volumetric fraction ϕ the temperature of the fluid is declined. Figure 9a ,b shows that for increasing the value of particle shape factor m temperature decreases and maximum temperature is illustrious for blade formed nanoparticle whereas minimum temperature is observed for bricked formed particles. Figure 10a ,b shows the influence of magnetic parameter M. It is observed that in the presence of external magnetic field temperature rises. Impact of fluid viscosity on temperature is represented in Fig. 11a ,b. It is observed that temperature of fluidic syste upsurges due to dependence of viscosity on temperature. Figure 12a In the present article, we discussed the results of mixed convection on MHD nanofluid and temperature dependent viscosity of the fluid by incorporating the strength of numerical computational approach. In this flow the stagnation point was towards the extending sheet. Copper-silver nanoparticles with base fluid were used for complete analysis. Vogel's model was taken to determine the impact of viscosity which is depending on temperature. Problem was expressed in the x-y coordinate system. The results of present analysis give the following key findings: • The fluid flow is decreased due to temperature dependent viscosity. • The fluid flow of assorted convection is affected by accelerating it. • The rise in volumetric fraction of nanoparticles causes to increase the thermal conductivity of the traditional fluid. • The impact of temperature was maximum for blade type nanoparticles of copper. • Skin friction coefficient decreases with the increasing value of Grashof Number Gr and increased with the increasing value of viscosity factors and volume fraction while these parameters reduced the Nusselt Number. In future, one may implement the Lobatto IIIA scheme for numerical treatment of many potential application arising in the fields of bioinformatics [42] [43] [44] , astro/plasma/atomic physics [45] [46] [47] , nonlinear circuit models [48] [49] [50] , fluid mechanics 51-56 , financial mathematics 57,58 and COVID-19 virus models 59,60 . www.nature.com/scientificreports/ Enhancing thermal conductivity of fluids with nanoparticles (No. ANL/MSD/CP-84938; CONF-951135-29) Embedded system for high performance digital temperature sensor on distributed parallel computer system Boundary layer behaviour on continuous solid surfaces-1 Flow past a stretching plate Similarity solutions of the boundary-layer equations for a stretching wall Boundary-layer flow of a nanofluid past a stretching sheet Nanofluid flow over an unsteady stretching surface in presence of thermal radiation Extraordinary role of hydrogen possessions and viscosity variation in electrically conducting copper and silver nanoparticles inspired by mixed convection Heat transfer over a stretching surface with variable surface heat flux A comprehensive shape factor analysis using transportation of MoS2-SiO2/H2O inside an isothermal semi vertical inverted cone with porous boundary On the numerical solution of stiff IVPs by Lobatto IIIA Runge-Kutta methods Numerical solutions of Williamson fluid with pressure dependent viscosity The effect of variable viscosity on the flow and heat transfer of a viscous Ag-water and Cu-water nanofluids Impact of viscosity variation and micro rotation on oblique transport of Cu-water fluid HPCL: integrated software for parameter characterization and mechanical properties of rubber nanocomposite Numerical treatment for Darcy-Forchheimer flow of Sisko nanomaterial with nonlinear thermal radiation by lobatto IIIA technique Numerical treatment for fluidic system of activation energy with non-linear mixed convective and radiative flow of magneto nanomaterials with Navier's velocity slip Comparative study of silver and copper water nanofluids with mixed convection and nonlinear thermal radiation Modern developments about statistical declaration and probable error for skin friction and Nusselt number with copper and silver nanoparticles. Chin Nanofluid flow and heat transfer in porous media: a review of the latest developments Natural convection MHD flow due to MoS2-Ag nanoparticles suspended in C2H6O2H2O hybrid base fluid with thermal radiation 3D mixed convection MHD flow of GO-MoS2 hybrid nanoparticles in H2O-(CH2OH) 2 hybrid base fluid under the effect of H2 bond Hydrothermal analysis of MHD nanofluid (TiO2-GO) flow between two radiative stretchable rotating disks using AGM. Case Stud Magnetic field effect on nanofluid flow between two circular cylinders using AGM Nanofluid heat transfer between two pipes considering Brownian motion using AGM. Alexandria Eng Application of He's variational iteration method to nonlinear Jaulent-Miodek equations and comparing it with ADM Evaluation of the heat transfer rate increases in retention pools nuclear waste Numerical study of fins arrangement and nanofluids effects on three−dimensional natural convection in the cubical enclosure Hybrid conduction, convection and radiation heat transfer simulation in a channel with rectangular cylinder Numerical analysis of secondary droplets characteristics due to drop impacting on 3D cylinders considering dynamic contact angle 3D optimization of baffle arrangement in a multi-phase nanofluid natural convection based on numerical simulation Numerical treatment for hydro-magnetic unsteady channel flow of nanofluid with heat transfer Numerical treatments to analyze the nonlinear radiative heat transfer in MHD nanofluid flow with solar energy Numerical computing paradigm for investigation of micropolar nanofluid flow between parallel plates system with impact of electrical MHD and Hall current Fuzzy conformable fractional differential equations: Novel extended approach and new numerical solutions Application of reproducing kernel algorithm for solving Dirichlet time−fractional diffusion-Gordon types equations in porous media Application of residual power series method for the solution of time−fractional Schrödinger equations in one− dimensional space Numerical solutions of systems of first-order, two-point BVPs based on the reproducing kernel algorithm Numerical algorithm for the solutions of fractional order systems of Dirichlet function types with comparative analysis Nanotechnology theory used for simulation of emerging big data systems on high performance computing: a conceptual framework High performance nanotechnology software (HPNS) for parameter characterization of nanowire fabrication and nanochip system Novel applications of intelligent computing paradigms for the analysis of nonlinear reactive transport model of the fluid in soft tissues and microvessels A stochastic computational intelligent solver for numerical treatment of mosquito dispersal model in a heterogeneous environment Integrated neuro-evolution-based computing solver for dynamics of nonlinear corneal shape model numerically A new heuristic computational solver for nonlinear singular Thomas-Fermi system using evolutionary optimized cubic splines Design of neuro-swarming-based heuristics to solve the third-order nonlinear multi-singular Emden-Fowler equation Integrated intelligent computing with neuro-swarming solver for multi-singular fourth-order nonlinear Emden-Fowler equation Two-dimensional brain microtubule structures behave as memristive devices Design of nature−inspired heuristic paradigm for systems in nonlinear electrical circuits Design and application of nature inspired computing approach for nonlinear stiff oscillatory problems Heat and mass transfer phenomenon for the dynamics of Casson fluid through porous medium over shrinking wall subject to Lorentz force and heat source/sink Effects of Gyro-Tactic organisms in bio-convective nano-material with heat immersion, stratification, and viscous dissipation. Arab Integrated intelligent computing paradigm for the dynamics of micropolar fluid flow with heat transfer in a permeable walled channel Dynamical analysis for nanofluid slip rheology with thermal radiation, heat generation/absorption and convective wall properties Hydromagnetic mixed convective flow over a wall with variable thickness and Cattaneo-Christov heat flux model: OHAM analysis Nanoparticles and nonlinear thermal radiation properties in the rheology of polymeric material Wavelets optimization method for evaluation of fractional partial differential equations: an application to financial modelling Fractional neuro-sequential ARFIMA-LSTM for financial market forecasting A stochastic intelligent computing with neuro-evolution heuristics for nonlinear SITR system of novel COVID-19 dynamics Intelligent computing with Levenberg-Marquardt artificial neural networks for nonlinear system of COVID-19 epidemic model for future generation disease control The authors declare no competing interests. Correspondence and requests for materials should be addressed to M.A.Z.R. Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.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://creat iveco mmons .org/licen ses/by/4.0/.