ACombined Convection Carreau–YasudaNanofluidModel over a Convective Heated Surface near a Stagnation Point: A Numerical Study

Department of Mathematics, University of Gujrat, Gujrat 50700, Pakistan Department of Mathematics, Quaid-I-Azam University, Islamabad 44000, Pakistan Department of Mathematics, College of Sciences, King Khalid University, Abha 61413, Saudi Arabia Department of Mathematical and Computer Modeling, Al-Farabi Kazakh National University, Almaty, Kazakhstan Department of Mathematical and Computer Modeling, Kazakh British-Technical University, Almaty, Kazakhstan Department of Computer Science, University of Sahiwal, Sahiwal, Pakistan


Introduction
Today, enhancement in the heat transfer mechanism using cooling fluids has become very significant in industrial applications such as machining and electronics, energy production industries, petrochemical, aeronautics, and transportation industries. e use of irreversible strength resources and the creation of environmental provisions have supported this improvement. e objective of manufacturing heat transfer devices requires reducing cost and increasing heat transfer in each surface area unit to achieve high efficiency. In recent years, technical advances have made improvements in the rheological properties of cooling fluids, and the manufacture of a steady solid-liquid suspension called nanofluid has optimized the thermal performance of industrial instruments and heat exchangers [1]. erefore, several researchers have reported enhancement in heat transfer through nanofluids [2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21]. e combination of forced and free convection is known as mixed convection. Extensive cases of mixed-convection are characterized as internal mass forces, where the flow is determined collectively by a few external forcing systems (i.e., external power supplied to the fluid-streamlined body system).
is is a distribution of fluid in a gravity domain with differing densities. e acceleration of the Earth's temperature stratified mass of air and water region is the realistic presentation of mixed convection that is conventionally studied in geophysics. However, in many engineering gadgets, mixed convection is found within the system of much smaller scales. is will be encapsulated on account of certain cases alluded to in channel flows. Mixed convection is often perceived in normal and traditional situations, i.e., heating or cooling of channel walls, limited velocities of a fluid flow that are characteristics of a laminar flow. Turbulent channel flows studies with considerable gravity area impacts have evolved since the 1960s after its essential use in industrial applications in nature, e.g., electronic appliances, cooling through electric fans, solar panels exposed to wind currents, atomic reactor cooled during emergency shutoff, flow in the ocean and the atmosphere, heat transmission placed in a low-velocity environment, and so on. Cesarano [22] discussed several fractional generalizations of different types, starting from the heat equation, and suggested a method useful for analytical/numerical solutions. Cesarano [23] introduced a nonconventional approach of multi-dimensional Chebyshev polynomials. Mixed convection flow was also reported by many researchers [24][25][26][27][28][29][30].
Commonly, fluids are categorized into two main classifications: non-Newtonian and Newtonian fluids; the main difference between Newtonian and non-Newtonian fluid is the connection between the deformation rate tensor and the extra stress tensor [31]. In a 3D free-form extrusion printing process, the rheological property plays an important role. But no rheological model was available in the open literature that could accurately consider the impacts of both the concentration of nano/microparticles and the non-Newtonian viscosity in the paste [32]. Zuo and Liu [32] introduced a fractal rheological model for 3D print pastes by applying a fractional derivative and it can be also used for other non-Newtonian fluids. e experimental findings are in strong agreement with the theoretical forecast, suggesting that the model is accurate and realistic. Defining the behaviors of melts and many polymer solutions, Carreau-Yasuda model was established successfully. It was accurately applicable to many experimental viscosity patterns with a wide range of shear rates [33]. Khan et al. [34] investigated the dissimilar solutions of MHD Carreau-Yasuda fluid with slip conditions flow past a rotating disk. Andrade et al. [35] presented the turbulent flow of the non-Newtonian Carreau-Yasuda fluid. Salahuddin et al. [36] examined the influence of transverse magnetic field on the squeezed twodimensional flow of a free stream past a sensor surface of an electrically conducting Carreau-Yasuda fluid. Kumar et al. [37] discussed pulsatile blood flow through the human carotid artery by applying computational fluid dynamics.
Inspired by the above research work, the ambition of the current investigation is to study the impact of mixed convection on Carreau-Yasuda stagnation point nanofluid flow over a vertical elastic surface by applying a mathematical nanofluid model introduced by Buongiorno [38]. In our point of view, the problem is new and original. e outcomes are discussed through various parameters such as the Brownian motion parameter, mixed convection parameter, thermophoresis parameter, buoyancy ratio parameter, stretching parameter, Lewis, Weissenberg, Prandtl, and Biot numbers on the mass and heat transfer.

Modeling
We inspect the incompressible, 2D, steady flow on a stagnation point attaining the stretched surface in a Carreau--Yasuda nanofluid as shown in Figure 1. Sheet is stretched with velocity U w (x) � ax and ambient fluid velocity is U E (x) � bx while the origin is fixed at S. It is assumed that the surface is warmed up by convection from a hot fluid at the temperature T f and the sheet is warmed up by convection fluid at the temperature T ∞ which is by heat transport coefficient h f . e surface motion will source the evolution of the boundary layer. Consider the Cartesian coordinates system; fluid flow velocity will change through x and y axis in a manner that the x-axis is taken vertically and the y-axis is taken horizontally. e elemental mathematical form of Carreau-Yasuda fluid is [39] where Γ and d are Carreau-Yasuda fluid parameters, A 1 is first Rivlin Ericksen tensor, μ ∞ is infinite shear rate viscosity, μ 0 is zero shear rate viscosity, τ is extra stress tensor, and _ c is . Take into account that infinite shear rate viscosity μ ∞ � 0 and then equation (1) reduces to the form Under these considerations, governing equations can be written as follows:

Mathematical Problems in Engineering
Here, ] is the viscosity characteristics and u and v are the velocity components in the x and y directions, respectively, where (ρc p ) p /(ρc p ) f is the proportion of nanoparticles capacity of heat to that of the base fluid capacity of heat, ρ f is the density of the base fluid, α � K/(ρc p ) f is the thermal diffusibility of the base fluid, c p is the constant pressure specific heat, D B is the Brownian diffusion, D T is the thermophoresis diffusion, and the subscripts np, ∞, and f denote nanoparticles, a value very far away from the solid surface, and the base fluid, respectively. e subjected boundary conditions are defined as Appropriate similarity transformations are defined as Applying the above relationship, the continuity equation (3) is satisfied. Equations (4)-(6) take the following form: with boundary conditions Here, prime denotes derivation for η, and further dimensionless parameters are defined as Here, N r , λ, Bi, W e , and Le denote buoyancy ratio parameter, mixed convection parameter, Biot number, Lewis number, and Weissenberg number respectively. Pr represents the Prandtl number, r is the stretching characteristic, Nb is the Brownian motion characteristic, and Nt is the thermophoresis characteristic. e skin friction coefficient C f , local Nusselt number Nu x , and local Sherwood number Sh x are where q w presents surface heat flux, τ w denotes surface shear stress, and q m defines surface mass flux for Carreau-Yasuda fluid as follows: By using appropriate similarity transformations (8), the expressions for dimensionless Nusselt number, skin friction, and the Sherwood number become ,

Results and Discussion
Nonlinear differential equations (9)-(11) with boundary conditions (12) and (13) are solved by using the MATLAB program bvp4c method. ese ODEs (9)-(11) with boundary conditions (12) and (13) are also solvable by Homotopy perturbation and variational iteration method [40,41]. For physical significance, the graphical and numerical outcomes of the solution will be taken into account. Figure 2 displays the effect of λ on the velocity patterns f ′ (η) for distinct values of λ. It may be noticed that field of velocity increases when (λ > 0) and decreases (λ < 0). Figure 3 demonstrates the influence of Nt on f ′ (η). It has been observed that the velocity curveincreases while incrementing the values of Nt. e contrasting behaviour is shown for opposing flow in Figure 3. Figure 4 shows the effect of N r on f ′ (η), inclined value of N r causing the boundary layer thickness upward and downward for reverse behavior. It is revealed that the value of n increases, and the flow of f ′ (η) expands. e observation is the same for opposing flow in Figure 5. Figure 6 shows that the profile of velocity increases when the values of Biot number enhancing. In general, Prandtl number has an impact on f ′ (η) reversed flow for the larger value of Pr in Figure 7. Figure 8 show that velocity distribution of the stretching parameter increases for expanding the value of r. It may be noted that W e is enhanced, and the flow of velocity moves upward in Figure 9.
Deviation of temperature flow is discussed in Figures 10-14. Figure 10 represents the behavior of Biot number, the field of temperature inclined for a larger value of Bi. Figure 11 depicts that the larger value of the Prandtl number profile of temperature goes down rapidly. Figures 12  and 13 represent Le and Nb, deeper presentation shown for assisting flow and upper for opposing flow on the temperature profile. Figure 14 shows the influence of the thermophoresis parameter on the thermal field; temperature flow moves downward for a larger value of Nt. Nanoparticle concentration field dispute is shown in Figures 15-19. Figure 15 shows that boundary layer thickness declines by increasing the values of Le,and on the other hand, boundary layer thickness inclines for opposing flow. Figure 16 describes the impact of Biot number on the concentration field. Concentration field declines by increasing the values ofBi. Figure 17 demonstrates the impact of Pr on φ(η) field; field rises with the increment of Pr. Figure 18 shows that concentration curve goes up for the Mathematical Problems in Engineering inclined value of Nb and in the reverse flow curve reduces for enlarge value of Nb. It has been noted that in Figure 19 concentration distribution declines to expand the values of Nt. e opposite behavior is shown for opposing flow. e results of the Nusselt number, skin friction, and Sherwood number are described numerically. Table 1 investigates the heat transfer parameter expanding if the values of N r , Nb, Le, Pr, λ, Bi, Nt inclined. Table 2 presents the impact of N r , Nb, Le, Pr, λ, Bi, and Nt on the Sherwood number. e values of Nb increase, if the mass transfer number rapidly inclines. e numerical amount of the mass transport coefficient diminishes if the values of these parameters N r , Nt, Le, Pr, λ, Bi rise. Table 3 represents the effect of a distinct amount of the stretching parameter and other parameters on the skin friction coefficient numerically. Skin friction enlarges when the values of N r , W e , Le, and λ increase and declines if the amount of Nb expands. In this paper, we cannot figure out the effect of the particle's size and distribution on the flow properties; the two-scale fractal calculus has to be adopted for this purpose [32,42]. Furthermore, streamlines for distinct values of stretching parameter are shown in Figures 20-22, which show that fluid flows in the same direction.

Concluding Remarks
In this manuscript, the problem of combined convection Carreau-Yasuda stagnation point nanofluid flow enclosed by a stretchable surface is discussed. e problem is solved numerically by applying the bvp4c technique through MATLAB. e main findings of the highlighted results are as follows: (i) Velocity profile accelerates by increasing the amount of λ, Nt, N r , n, W e , r, Bi and decelerates, when decreasing the amount of Prandtl number. (ii) It is observed that the concentration profile declined when Bi, Le, and Nt decreased. It is also revealed that, with the enlargement in Pr and Nb, the concentration pattern expands.
(iii) e temperature profile increases with an increase in the value of Biot number (Bi) and falls when the values of Pr, Le, Nb, and Nt increase.

Data Availability
All data are included within this manuscript.

Conflicts of Interest
e authors declare that they have no conflicts of interest.