Computational analysis of water based Cu - Al2O3/H2O flow over a vertical wedge

Theoretical and numerical investigation of the fluctuating mixed convection of hybrid nanofluid flow over a vertical Riga wedge is considered in this analysis. Two kinds of solid nanoparticles with base fluid at vertical Riga wedge is studied. Thermal and velocity slip impacts on vertical Riga wedge are investigated in the current study. We discussed both the unsteady and steady cases. The water has low thermal conductivity. We added the nanoparticle Cu and A l 2 O 3 which increases the thermal conductivity of the base fluid. This phenomena increase the heat transfer rate at the surface Riga plate. Partial differential equations are reduced into an ordinary differential equation by means of dimensionless similarity variables. The resulting ordinary differential equations are further solved through numerical and perturbation methods. Thickness of momentum boundary layer is reduced because of the solid nanoparticle rises in all cases of ξ = 0 , ( 2 i ξ ) 0 , ( 2 i ξ ) 1 and ( 2 i ξ ) n . Our results are more agreeing with the decay results of Bachok et al. and Yacob et al. when rest of the physical parameters dimensions.


Introduction
Heat transfer and boundary layer flow at vertical wedge have been considered much attention in recent days because of various physical applications of metallurgical and engineering, namely: plastic and metal extrusion, hot rolling, manufacturing of the fiberglass, crystal growing, wire drawing, manufacturing of paper and continuous casting. In the early time, Crane 1 pioneered boundary layer flow caused by stretching the surface. Carragher and Crane 2 highlighted that the temperature difference between the ambient fluid and surface was proportional to the power of the distance from the fixed point. The conduct of Newtonian nanofluids could be valuable in assessing the chance of heat transfer improvement in different procedures of these industries. Wu 3 studied the boundary layer flow of water entry of twin wedges numerically. Micropolar fluid flow at convective stretching sheet discussed by Haq et al. 4 Ramesh 5 explored the results of slip effects with convective condition of peristaltic flow in a porous medium. Ramesh 6 investigated the effects of inclined MHD in the inclined asymmetric channel. Gireesha et al. 7 investigated the hall impacts with MHD nanofluid and also worked on irregular heat consumption/generation. Gireesha et al. 8 discussed the impacts of hall and thermal radiation impacts at porous medium in the presence of suspension of fluid particles. Mahanthesh 9 explored the influence of dusty liquid with hall effects at the stretching sheet. Gireesha et al. 10 analyzed the hall effects on nanomaterial fluid at linear stretching sheet under the hall effects. Gireesha et al. 11 examined the dusty fluid flow with hall impacts under the nonlinear radiation at the stretching plate. Shalini and Mahanthesh 12 explained the nanomaterial dusty fluid with Rayleigh-Benard convection without/with Coriolis force. Recently, few investigators have discussed the behavior of fluid at various stretching surfaces under various assumptions at sheet see Ramesh et al., 13 Mahanthesh et al. 14 and Wakif et al. 15 Within the past few a protracted time, researchers have confronted a substantial issue of thermal efficiency in numerous engineering and industrial applications. Advance technology such as microfluidics, optical, chemical synthesis, high speed microelectronics, transportation, microsystems including electrical components and mechanical bears high thermal loads. In such days, cooling specific systems is a key problem. The rate of heat transfer is typically enhanced by adding the region accessible for heat transfer. A further approach is to use a thermally effective fluid, which is typically distributed inside traditional heat transfer fluids, such as water, propylene glycol, ethylene glycol, etc. Heat transfer is most important in these days due to an energy crisis in the whole world. Choi and Eastman 16 was the first person who worked on nanofluid. He used the nanoparticle with base fluid like water, oil, and ethylene glycol. Before this idea, several authors analyzed the heat transfer rate enhancement. But Choi and Eastman 16 achieved heat transfer rate more than two times as compared to conventional flow because low thermal conductivity. Lee et al. 17 used the oxide nanoparticle. The properties of the nanofluid discussed by Das et al. 18 They additionally brought a physically and chemically method for synthesizing nanofluids and also defined the methods for making a stable suspension of nanoparticles with the base fluids. They also analyzed the stagnation point flow at the convective stretching sheet. Mixed convection of nanomaterial fluid flow at porous wedge explored by Ellahi et al. 19 Most interesting research has done by few authors on nanofluid in Ellahi et al., 19 Usman et al., 20 Sheikholeslami et al., 21 Ramesh and colleague, 22,23 Zeeshan et al., 24 Liu et al. 25 Hybrid nanofluid more important role in the heat transfer enhancement because of their lots of applications namely, biomedical, aselectronic cooling, defense, generator cooling, automobile radiators, nuclear system cooling, lubrication, coolant inmachining, welding, thermal storage, solar heating, cooling and heating in buildings, drug reduction and biomedical etc. At the mechanical level, there are positive highlights like chemical soundness and tall thermal efficiency, which empowers them to perform efficiently as compared to nanofluids. Furthermore, hybrid words were introduced by Suresh et al. 26 He claimed that two different kinds of particles are used with fluid. He proved that the heat transfer rate enhancement was two times more than nanofluid. The experimental results of hybrid nanofluid explored by Suresh et al. 26 and Labib et al. 27 The turbulent flow of hybrid nanofluid analyzed at a low volume concentration by Sundar et al. 28 Stagnation point flow of hybrid nanomaterial fluid flow explored by Nadeem et al. 29 Slip effects of hybrid nanomaterial fluid flow at cylinder analyzed by Nadeem and Abbas. 30 Several scientists are being researched the heat transfer of hybrid nanofluid for different physical assumptions reported in Nadeem and colleagues, 31,32 Mahanthesh and colleagues, 14,33 Nehad et al. 34 Mixed convection of hybrid nanofluid flow dusty fluid over a vertical Riga wedge plate is considered. Mathematical model has been constructed on the assumptions of fluid flow. The partial differential equations have been converted into systems of ordinary differential equations using the suitable similarity transformations. Nonlinear ordinary differential equations are solved through perturbation scheme and numerically through BVP4C method. The solutions are deliberated for numerous physical parameters.

Mathematical formulations
Consider the mixed convection fluctuating flow of hybrid nanofluid over a vertical Riga wedge in the presence of solid nanoparticle. The coordinate system and flow configuration reveals in Figure 1. The thermal and momentum boundary layer thickness is represented by d M and d T . It is considered that both the free stream and surface temperature present small amplitude oscillations in time about steady, non-zero free stream velocity and mean temperature. The electro-magnetohydrodynamic applied on the vertical wedge see Wakif et al. 15 Under the assumption and usual boundary layer approximation, the following mathematical system has been considered as below Where, T = T À T ' , the boundary conditions are as following The coefficient of the volumetric expression for temperature expresses as b, V e X , t ð Þ be the fre stream velocity, V and W be the local velocity components of the fluid, r hnf , m hnf , a hnf and rC p hnf is the density, viscosity, thermal conductivity and heat capacity of hybrid nanofluid. Under these conditions, we combine the energy equations (3) and (6), momentum equations (2) and (5) for the two phases which we expressed as Corresponding boundary conditions are Where The mean surface temperature and mean velocity represent as T 0 X ð Þ and V 0 X ð Þ, e is the amplitude of oscillation which is e ( 1 while the frequency of oscillation is v:

Solution procedure
On the behalf of the boundary conditions given in equation (9), we deliberate solutions of the equations (6)(7)(8) in the following form Where V 1 , V 2 are the fluctuating parts and W 1 , W 2 are the velocity components. V , W and T 1 , T 2 are the real and fluctuating part of the temperature function T . The predictable solution functions are the real parts of the functions known in equation (10). The thermophysical properties are presented in Tables 1 and 2. After applying the equation (10) in equations (6-9) and getting the terms and conditions to 0 e ð Þ, we obtained succeeding equations With the boundary conditions are And Relevant boundary conditions are Now we consider functions V 0 X ð Þ and T 0 X ð Þ like as V 0 X ð Þ= V c X 1=2 and T 0 X ð Þ= T c while V c and T c are the constants. We introduced the similarity following transformations are Where C 1 and C 2 are the stream function which satisfy the equations 11Þ and 15 ð Þ ð correspondingly. The similarity and local frequency variables are h and j. By applying the similarity transformation equation (23) on equations (20)(21)(22), we get the dimensionless form The steady part of the flow presented as follows Respected boundary conditions Perturbation solutions for small j To gain the impacts of mixed convection near the leading edge, the results based on the series 'finite numbers of terms are only valid for the small frequency range'. So Putting these values in equations (25)(26), than equating the coefficients of 2ij Respected boundary conditions are Equating the coefficients of 2ij ð Þ 1 Relevant boundary conditions are Equating the coefficients of 2ij ð Þ n 1 + Pr ð Þn hnf G n 000 + 3 = 4 FG n 00 + 3 = 4 + n = 2 n o F 00 G n Relevant boundary conditions are

Result and discussions
The influence of the physical parameters on the heat transfer and skin fraction is presented in Table 3. The skin fraction improved because of improving the solid nanoparticle. Heat transfer rate declines as well as solid nanoparticle enhances which reveals in Table 3. The heat transfer rate 2ij ð Þ n (Unsteady flow) achieve much higher than that of the j = 0 (Steady flow) in case of the Re 1=2 X C f and Re 1=2 X N u X . Skin friction reduces when the slip parameter rises. The thickness of momentum is reduced by cause of the enhancing velocity slip parameter. While the heat transfer improved with improving the velocity slip parameter. Because thickness of thermal boundary rises by cause of velocity slip enhances. Incase of j = 0, variation of the modified Hartman number on the skin friction and heat transfer revealed in Table 3. Skin friction reduces when the modified Hartman number rises. The thickness of momentum is reduced by cause of the enhancing modified Hartman number but opposite results are achieved for the 2ij ð Þ n . While the heat transfer reduced with improving the modified Hartman number because thickness of thermal boundary declines by cause of velocity slip enhances but opposite results are achieved for the 2ij ð Þ n . In case of j = 0, variation of the Richardson number on the skin friction and heat transfer revealed in Table 3. Skin friction improves when the Richardson number rises. The thickness of momentum is improves by cause of the enhancing Richardson number but opposite results are achieved for the 2ij ð Þ n . In case of j = 0, the heat transfer improved with improving the Richardson number because thickness of thermal boundary improves by cause of velocity slip enhances but opposite results are achieved for the 2ij ð Þ n . In case of j = 0, variation of the dimensionless number on the skin friction and heat transfer revealed in Table 3. Skin friction improves when the dimensionless number rises. The thickness of momentum improves by cause of the enhancing dimensionless number, but opposite results are achieved for the 2ij ð Þ n . In case of j = 0, the heat transfer improved with improving the dimensionless number because thickness of thermal boundary improves by cause of velocity slip enhances but opposite results are achieved for the 2ij ð Þ n . Impacts of the physical parameters on the modified Hartman number, thermal slip, velocity slip and solid nanoparticle on temperature gradient highlighted in Figures 2-5 for the various cases j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . Figure 2 reveals the variation of thermal thickness and thermal slip. Results are achieved that thermal thickness reduced by cause of improving the thermal slip parameter for all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . But in case of the thermal thickness reduced 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n respecticly, as the values of n enhances at the surface of Riga vertical wedge. Table 4 presents the comparative results of decay literature. Figure 3 reveals the variation of thermal thickness and solid nanoparticle concentration F 2 . Results are achieved that thermal thickness rised by cause of improving the solid nanoparticle for all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n because thermal thickness increases due to enhancing the thermal conductivity of the fluid increases. When the solid nanoparticles added in the base fluid which enhances the thermal conductivity. But in case of the thermal thickness improved 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n respecticly, as the values of n enhances at the surface of Riga vertical wedge. Figure 4 reveals the variation of thermal thickness and modified Hartman number. Results are achieved that thermal thickness rised by cause of improving the modified Hartman number for all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . But in case of the thermal thickness improved 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n respecticly, as the values of n enhances at the surface of Riga vertical wedge. This phenomena exist due to Lorentz forces F m = ÀF m , 0, 0 ð Þ which explained physically.
The electromagnetic forces F m = 1 worked as opposite along to the motion of fluid as drag-like forces which successful extents are rotting exponentially along the y axis. Figure 5 reveals the variation of thermal thickness and velocity slip. Results are achieved that thermal thickness reduced by cause of improving the velocity slip parameter for all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . But in case of the thermal thickness reduced 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n respecticly, as the values of n enhances at the surface of Riga vertical wedge. Impacts of the physical parameters on    the modified Hartman number, thermal slip, velocity slip and solid nanoparticle on the velocity distribution highlighted in Figures 6-9 for the various cases j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . Variation of the thermal slip and velocity profile presented in Figure 6. Thickness of momentum reduced by cause of the thermal slip rises in all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . Momentum thickness decreases by cause of enhancing the values of n. Variation of the velocity slip and velocity profile presented in Figure 7. Thickness of momentum improved by cause of the velocity slip rises in all cases of j = 0., 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . Momentum thickness decreases by cause of enhancing the values of n. Variation of the solid nanoparticle and velocity profile presented in Figure 8. Thickness of momentum reduced by cause of the solid nanoparticle rises in all cases of j = 0., 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . because thermal thickness increases due to enhancing the thermal conductivity of the fluid increases. Momentum thickness decreases by cause of enhancing the values of n. Variation of the modified Hartman number and velocity profile presented in Figure 9. Thickness of momentum improved by cause of the modified Hartman number rises in all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . Momentum thickness decreases by cause of enhancing the values of n.

Final remarks
We considered the mixed convection fluctuating flow of hybrid nanofluid over a vertical Riga wedge. Major    effects of the physical parameters on the temperature, velocity, heat transfer and skin friction on the Riga vertical wedge are highlighted as below: The thickness of momentum is reduced by cause of the enhancing velocity slip parameter while the heat transfer improved with improving the velocity slip parameter. In case of j = 0, the heat transfer improved with improving the Richardson number because thickness of thermal boundary improves by cause of velocity slip enhances but opposite results are achieved for the 2ij ð Þ n . Momentum thickness decreases by cause of enhancing the values of n ( 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n ). Thickness of momentum reduced by cause of the solid nanoparticle rises in all cases of j = 0, 2ij ð Þ 0 , 2ij ð Þ 1 and 2ij ð Þ n . Our results are more agreeing with the decay results of Bachok et al. 35 and Yacob et al. 36 when and rest of the physical parameters dimensions.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.