Significance of nanoparticle’s radius, heat flux due to concentration gradient, and mass flux due to temperature gradient: The case of Water conveying copper nanoparticles

The performance of copper selenide and effectiveness of chemical catalytic reactors are dependent on an inclined magnetic field, the nature of the chemical reaction, introduction of space heat source, changes in both distributions of temperature and concentration of nanofluids. This report presents the significance of increasing radius of nanoparticles, energy flux due to the concentration gradient, and mass flux due to the temperature gradient in the dynamics of the fluid subject to inclined magnetic strength is presented. The non-dimensionalization and parameterization of the dimensional governing equation were obtained by introducing suitable similarity variables. Thereafter, the numerical solutions were obtained through shooting techniques together with 4th order Runge–Kutta Scheme and MATLAB in-built bvp4c package. It was concluded that at all the levels of energy flux due to concentration gradient, reduction in the viscosity of water-based nanofluid due to a higher radius of copper nanoparticles causes an enhancement of the velocity. The emergence of both energy flux and mass flux due to gradients in concentration and temperature affect the distribution of temperature and concentration at the free stream.

The roles of energy and mass fluxes due to temperature and concentration gradients play major roles in electrochemical processes, dynamics of gases, chemical catalytic reactors, production of copper selenide sequel to the pioneering work by Louis Dufour 1 ; see Ingle and Horne 2 , Hollinger and Lucke 3 , Korzhuev 4 , Kim et al. 5 , Mahdy 6 , Anjum Badruddin 7 . Conservation of energy is capable to unravel the nature of outward flux of energy. Meanwhile, Hort et al. 8 once remarked that concentration currents and heat are two driven forces of fluctuations in temperature as in the case of conservation of energy. The effect of energy flux due to concentration gradient on six different gases ( CO 2 + O 2 , N 2 + H 2 , H 2 + CO 2 , N 2 + O 2 , CO 2 + N 2 , H 2 + O 2 ) at various levels of pressure was examined experimentally by Rastogi and Madan 9 . It was remarked that kinetic theory is useful to determine the Dufour coefficient of each gas and temperature difference within the domain (T w − T ∞ ) is dependent on pressure that maintains the motion of each gas. Formation of energy flux due to temperature gradient is highly significant in gases than in liquids. In fact, the occurrence in a less viscous fluid is 10,000 larger than in high viscous fluid. This fact suggests that the magnitude of the Dufour number is larger for gases (Hort et al. 10 ). Dynamics of peristaltic flow through a channel of width 2a subject to energy flux due to concentration gradient and mass flux due to energy gradient was closely examined by Hayat et al. 11 . It was shown that temperature distribution is an increasing property of both Dufour and Soret effects. Although, the increase in the temperature distribution is more enhanced near the surface. The reverse is the observed effects of Dufour and Soret phenomena on the concentration as both properties decreases.
Further examination of energy and mass fluxes by Linz 12 led to a conclusion that a larger Dufour number was recommended for gaseous mixtures in which energy flux due to concentration gradient is significant. Following the suggestion by Partha et al. 15 , Lewis number is more appropriate than Schmidt number in characterizing heat and mass transfer due to the significance of energy flux due to concentration gradient and mass flux due to temperature gradient. This led to a robust analysis of the interconnectedness of Lewis number, Dufour number, and Soret number. The results show that when L e = 0.1 and S r = 0.6 , Nusselt number Nu x Ra −0.5 x decreases with Dufour effect at the rate of −0.076145399 . But, when Lewis number is more higher (i.e. L e = 1 ) and Soret number still remain S r = 0.6 , the same Nusselt number Nu x Ra −0.5 x decreases with Dufour effect at higher rate of −0.585851817 . Energy flux due to concentration gradient was pointed out by Garcia-Colin et al. 13 as the most important source of heat conduction within the region where magnetic fields do not exist or even weak. The further associated increase in effective thermal conductivity to growth in Dufour effect. The major reason why temperature increases due to higher magnetic strength were associated with the fact that heat conduction in the perpendicular direction diminishes whenever magnetic field is intensified. The conclusions mentioned above are not affected in the correction published as Ref. 14 .
Enhancement in the transfer of heat across fluid flow has made experts in thermal engineering embrace the efficiency of nanofluids. The earlier mentioned advancement is based on the nature of the base fluid and nanoparticles. Nanoparticle concentration and temperature effects on the ratio of mass to density and viscosity are some of the physical properties. However, thermal conductivity and specific heat capacity at different levels of concentration of nanoparticles, nanoparticles' size, and temperature are some of the thermal properties. The concentration of nanoparticles, pressure drop, friction factor, nanoparticles radius are also some of the outlined characteristics of nanofluids as pointed out by Mohamoud Jama 16 , Narayanan and Rakesh 17 , Lin and Yang 18 . Intrinsic magnetic related properties of magnetite vary as to its diameter changes (i.e. higher ratio of surface to volume). As the size of nanoparticles becoming smaller, the superparamagnetic nature of magnetic nano-sized particles even changes 19 . Ashraf et al. 20 remarked that changes in the radius of nanoparticles affect the characteristics of both interphase and nanoparticles. According to Yapici et al. 21 , the outcome of comparative analysis of ethylene glycol conveying SiO 2 (20-30 nm, 60-70 nm), TiO 2 (30 nm, 50 nm), ZnO (10-30 nm, 35-45 nm, 80-200 nm), CuO (40 nm, 80 nm), and MgO (20 nm, 40 nm) shows that the relative viscosity of these nanofluids is a decreasing property of particle sizes. Vishal et al. 22 noted that the viscosity of nanofluid can be greatly influenced by particle size, the nature of energy transfer between fluid's layer and the surface of the particles. There exist significant changes in the melting point of nanoparticles due to rises in the radius. This conclusion was based on the fact that the melting temperature was seen to be a decreasing property of higher particle size with a significant decreasing rate between spherical and nanoparticles; see Antoniammal and Arivuoli 23 .
The analysis of Namburu et al. 24 indicates that increasing the diameter of silicon dioxide nanoparticles in ethylene glycol and water corresponds to a decrease in the nanofluid's viscosity. It is worth remarking that the observed decrease in fluid's viscosity is more significant when the nanofluid is cold (negative temperature). At a larger temperature, the decrease in viscosity with particle size disappears. In the case of a nanofluid parallel to a moving stretchable sheet, platelet shape of molybdenum disulfide (MoS2) nanoparticles was found by Hamid et al. 25 to produce a unique heat transfer. In another study, Sheikholeslami et al. 26 illustrated the movement of multiple wall carbon nanotube and Iron(iii)oxide nanoparticles in a typical water based nanofluid through a porous medium when Lorentz force is predominant. The dynamics of water conveying alumina and copper nanoparticles over a split lid-driven square cavity was examined by Khan et al. 27 and it was shown that higher Nusselt number proportional to the heat transfer rate is attained at the point when both lids meet. In another related report by Khan et al. 28 , carbon nanotubes and heating of the wall are two factors capable to boost the local Nusselt number in the case of water-based carbon nanotubes over a right-angle trapezoidal cavity. Meanwhile, the local Nusselt numbers are increasing property of solid volume fraction of carbon nanoparticles; Hamid et al. 29 . The analysis of seven different hybrid nanofluids by Nehad et al. 30 shows that optimal Nusselt number is achievable when suction and stretching are significantly large but less dense nanoparticles like silicon dioxide and multiple wall CNT.
Sequel to the aforementioned reviews of related literatures, in the presence of Joule heating and spacedependent heating, it is noteworthy to examine the significance of increasing nanoparticle radius and inclined magnetic field on the dynamics of chemical reactive water conveying copper nanoparticles through a porous medium. The outcome of such study would be very helpful to experts dealing with the performance of copper selenide and effectiveness of chemical catalytic reactors. This study was designed to provide answers to the following related research questions:

Mathematical formulation
When energy flux due to concentration gradient and mass flux due to temperature gradient are significant, the dynamics of water conveying copper nanoparticles over a horizontal surface subject to magnetic field of strength B o inclined at an angle γ was formulated. Due to the stretching at the wall, stretching velocity is assumed to be u w = U o x where < 0 implies shrinking while > 0 implies stretching at the wall y = 0 . The chemical reaction that occurs in the transport phenomenon was model as first order where the rate of the chemical reaction is Kc * . Following Bachok et al. 39 , the governing equation suitable to investigate the aforementioned transport phenomenon is

Suitable boundary conditions are
The model proposed by Graham 31 and Gosukonda et al. 32 for the ratio of dynamic viscosity of the nanofluid to the dynamic viscosity of base fluid defined as where the radius of nanoparticle is d p and the inter-particle spacing is h was adopted. The effective nanofluid properties are given by where φ is the solid volume fraction, µ f is the dynamic viscosity of the base fluid, ρ f and ρ s are the densities, (ρc p ) f and (ρc p ) s are the heat capacitance, κ f and κ s are the thermal conductivities and σ f and σ s are the electrically conductivities of the base fluid and nanoparticles respectively. Prior to one of the conclusions by Buongiorno et al. 33 on the thermal conductivity of nanofluids, the model proposed by Maxwell 34 was adopted to incorporate the enhancement in the thermal conductivity of Cu-Water nanofluid. Next is to introduce the following variables where D f is the Dufour number, S r is the Soret number, Prandtl number P r , magnetic field parameter M, and porosity parameter P, heat source/sink parameter Q e , Lewis number L e , viscous dissipation term E c , and suction S are defined as   38 , and Bachok et al. 39 . The dimensionless local skin friction, heat transfer rate, and mass transfer rate are

Method of solution, results and discussion
The method of superposition (Na 40 ) was used to obtain the system of first order IVP for Eqs. (10)- (14). The numerical solution of the corresponding IVP was obtained using MATLAB based bvp4c package suggested by Gokhan 41 , Kierzenka and Shampine 42 , and Ali Umit Keskin 43 . Reliability and validity of solution was established by comparing the limiting case of this present study with that of Hamad 44 at various values of Prandtl number. As shown in Table 1, good agreement is seen, hence further analysis is reliable. The simulation was carried out for fixed magnetic field strength M = 0.5 , inclination of the magnetic field γ = 30 o , inter-particle spacing h = 1 , Lewis number L e = 0.2 , Prandtl number P r = 6.2 , chemical reaction parameter K c = 0.5 , heat source parameter Q e = 0.5 , intensity of heat source n = 0.5 , viscous dissipation term-Eckert number E c = 0.3 , stretching related parameter = 0.5 , and suction S = 0.5 . The answer to the remaining research questions was obatined when .
Hamad 44 Wang 45 Gorla and Sidawi 46 51 , and Animasaun et al. 52 , the analysis presented as Tables 2, 3, 4 and 5 reveal that a larger radius of copper nanoparticles leads to negligible higher local skin friction coefficients. Firstly, when S r = 0.1 and D f = 0.1 , as d p increases, local skin friction coefficients increases at the rate of 0.046678714, higher heat transfer rate of 0.019499889, and higher mass transfer rate of 0.004874409. In order to capture the significance of energy flux   www.nature.com/scientificreports/ due to concentration gradient, the same analysis was carried out for D f = 10 ; see Table 3. Percentage increase in Nu x √ Re x with d p as D f changes from 0.1 to 10 was estimated 344.286288% . Comparative analysis of the data presented as Tables 4 and 5 and figures illustrated as Figs. 1 and 2 confirm that there exists a relationship between energy flux due to concentration gradient and mass flux due to temperature gradient.
When S r = 0.1 , the significance of increasing radius of copper nanoparticles d p and energy flux due to concentration gradient D f were examined. The outcome of the analysis presented as Figs. 3 and 4 show that the velocity of the transport phenomenon increases with a larger radius of nanoparticles. Meanwhile, an increment in the friction across the layers near the inclined surface is observed through an increase in the radius of nanoparticles as shown in Fig. 5. This can be associated with the fact that increasing the diameter of nanoparticles in water as pointed out by Namburu et al. 24 corresponds to a decrease in the nanofluid's viscosity. Not only that, the outcome of an examination of (i) ethylene glycol and alumina nanoparticles mixture (ii) water and CuO nanoparticles mixture by Pastoriza-Gallego et al. 48,49 shows that higher viscosity is bound to occur as particle size diminishes. It is worth deducing from Figs. 3, 4, and 5 that the velocity is the same at all the levels of heat energy flux due to concentration gradient D f . However, in the case of distribution of temperature illustrated as Fig. 6, decreases negligibly with d p when D f is small in magnitude. When D f = 5 and D f = 10 , the observed decrease in the temperature distribution due to the higher radius of nanoparticles is ascertained. Further exploration of temperature distribution for 0.1 ≤ D f ≤ 10 was achieved through kriging gridding method package in Surfer plot version 11.1.719 as shown in Figs. 7 and 8. When the transfer of species due to temperature gradient is small in magnitude optimal temperature is seen a few distances away from the inclined surface at larger values of D f . Now, when mass flux due to temperature gradient is sufficiently large, optimal temperature, although it is small in magnitude is seen at all the levels of energy flux due to the concentration gradient.

Conclusion
Attempt had been made to examine the significance of increasing radius of nanoparticles, energy flux due to concentration gradient, and mass flux due to temperature gradient in the dynamics of chemically reactive fluid subject to suction and inclined magnetic strength. Based on the analysis, it is worth concluding that 1. at all the levels of energy flux due to concentration gradient, reduction in the viscosity of water-based nanofluid due to a higher radius of copper nanoparticles causes an enhancement of the velocity. 2. significance decrease in distribution of temperature across the domain due to increasing radius of copper nanoparticles is achievable when energy flux due to concentration gradient is sufficiently large in magnitude. 3. when mass flux due to temperature gradient is highly negligible, optimal temperature is also observable when energy flux due to concentration gradient is sufficiently large. Reverse is the case when mass flux due to temperature gradient is sufficiently large as optimal temperature is ascertained at all levels of energy flux due to concentration gradient.    www.nature.com/scientificreports/ 4. reduction in the mass transfer rate Sh x √ Re x due to higher radius of copper nanoparticles is guaranteed when the transfer of species (mass) due to temperature gradient is sufficiently large but transfer of heat energy due to concentration gradient is highly negligible. 5. the emergence of both energy flux and mass flux due to gradients in concentration and temperature affect the distribution of temperature and concentration at the free stream.