#### Abstract

The thermally developed flow of nanoparticles due to bidirectional moving space has been presented numerically. The additional impact of heat source and sink are also incorporated for enhancing the heat transfer rate. The copper and aluminum nanoparticles are selected for enhancing the thermal diffusion system. The thermal stability of nanoparticles is carefully addressed. The moving surface allows the shrinking and stretching phenomenon. The system of resulting equations with nonlinear nature is numerically swapped via shooting technique with help of MATLAB software. The dual numerical simulations for the thermally radiative problem with external heating phenomenon are being performed. It is noted from results that thermal stability of nanoparticles is more stable in the upper branch of shrinking/stretching surface. The presence of heat source provides extra heat and improves the thermal transportation phenomenon which is more progressive in upper branch. The enhancement in temperature profile due to copper nanoparticles is more stable as compared to aluminum nanoparticles. Moreover, the increasing change in the velocity is resulted upon increasing the rotation constant. The obtained results reflect applications in thermal sciences, heating systems, energy production, solar applications, nuclear reactions, biomedical applications, etc.

#### 1. Introduction

Recently, motivations of researchers are noticed towards the nanofluid which reports leading contributions in the thermal engineering and industrial processes. The growing interest in nanomaterials is justified due to highly attractive applications in solar energy, thermal transportation systems, engineering devices, cooling systems, heating objectives, nuclear reactions, medical sciences, emission processes, etc. With nanostructured (1–100 nm), the nanofluids report excellent thermal activities. The interaction between nanoparticles and some base fluids reports results enhanced thermal mechanism of base fluids. Choi [1] focused first time on the thermal aspect of nanofluids via experimental approach. Later on, diverse research is being performed by investigators to highlight the thermal mechanism and applications of nanomaterials. Turkyilmazoglu [2] developed a slip nanofluid model for the annuli flow by utilizing the dual nanofluid and single relations of the nanofluid model. Nayak et al. [3] exploited the viscous dissipation consequences for illustrating the nanofluids flow with Oldroyd-B base material via bidirectional space. Basir et al. [4] performed the novel aspect of slip flow of nanofluid with sensitivity insurance and bioconvection phenomenon. The Maxwell nanofluid thermal determinations subject to the interaction of multiple types of slip consequences have been identified by Ahmad et al. [5]. The Buongiorno nanofluid model implementation for inspecting the viscoelastic fluid thermal capacity was reported in the work of Nadeem et al. [6]. The analytical simulation for the nanofluid Riga surface problem was investigated by Vaidya et al. [7]. Wang et al. [8] focused on the migration of nanofluid material minichannels having parallel plates. Abdelaziz et al. [9] discussed the characteristics of nanofluid pattern via horizontal tubes. The continuation of Sabu et al. [10] determined the heat transfer improvement for the alumina nanoparticles with slip constraints due to heated rotatory disk. Ojjela et al. [11] addressed the viscoelastic behavior and thermal capability for the nanoparticles with non-Newtonian base fluid. The stability determination of nanofluids with experimental approach was determined by Dong et al. [12]. Sivasankaran and Bhuvaneswari [13]investigated the heat transfer mechanism with constraints of pressure loss and porous space. Acharya [14] numerically focused on the natural convective flow of nanofluid due to heated wall surface in square enclosure. Mondal et al. [15] addressed the optimized thermal mechanism due to the couple stress nanofluid with radiative phenomenon. Acharya [16] discussed the finite element numerical simulations for the spinning flow of nanoparticles with active and passive control approach. Shafiq et al. [17] addressed the bioconvection aspect of nanofluid with zero mass constraints. Haya et al. [18] addressed the Joule heating effects while inspecting the enhanced aspect of heat transfer in third grade nanofluid. Some more recent work on nanofluid flow is presented in [19–24].

The hybrid nanofluid is the modified form of nanofluids with more improved thermal activities. In fact, the hybrid nanofluid model is based on the interaction of more than one different types of nanoparticles along with base material. The preferences of the hybrid nanofluid model over simple nanomaterials are due to extra high activities and more thermal stability. Recently, different investigations are performed on hybrid nanofluid with different nanoparticles. Sowmya et al. [25] reported the inclusion of iron nanoparticles with existence of Lorentz force in cavity with two heated fins. Sundar et al. [26] claimed the hybrid nanofluid thermal properties with collector applications. Akhtar et al. [27] modeled a hybrid nanofluid problem in elliptic duct following the peristaltic mechanism. Eid and Mabood [28] noticed the ethylene glycol base fluid thermal inspection with hybrid nanofluid numerically. Eid and Mabood [28] presented the thermal research on the carbon nanotubes flow with heat generation features. Shahid et al. [29] implemented the Darcy porous law for studying the hybrid nanofluid prospective in parabolic surface. Madhukesh et al. [30] explored the thermal efficiencies of AA7072-AA7075 nanoparticles with Newtonian heating enrollment. Acharya [31] discussed the shape characteristics for the hybrid nanofluid in heated obstacles. Devi and Devi [32] visualized the copper and aluminum oxide nanoparticles thermal prospective accounted via stretching surface. Devi and Devi [33] observed the 3D thermal flow problem with Newtonian heating. Lund et al. [34] intended the rotating flow of hybrid nanofluid. Teh and Ashgar [35] determined the Joule heating for the rotating flow of hybrid nanofluid. Zaimi et al. [36] discussed the multiple solution of rotating viscoelastic fluid with hybrid nanofluid.

This analysis deals with the impact of heat source and thermal radiation effects in flow of hybrid nanofluid due to rotating space. The main motivation for performing this investigation is to highlight the enhanced thermal aspect of heat transfer phenomenon subject to the magnetized copper and aluminum nanoparticles. The moving surface attained the stretching and shrinking phenomenon. The thermal stability of hybrid nanoparticles is focused. The shooting numerical algorithm is followed for the simulations process. The dual simulations are performed in each branch of moving space.

#### 2. Problem Formulation

The stretching and shrinking phenomenon of 3D hybrid nanofluid is taken into account. The heat transfer analysis is considered in presence of external heating source and thermal radiation. The magnetic force impact is directed along the *z*-axis for 3D flow. The three-dimensional moving space attained the velocity The stretching and shrinking results are attained for and respectively. The surface temperature of nanoparticles is notified with while is for free stream temperature case as shown in Figure 1. The nanoparticles and surface present rotatory behavior with angular velocity along *z*-direction. Let , , and be representation of velocities along and directions, respectively. The assumptions of low magnetic Reynolds number reduce the role of induced magnetic force. The governing relations in view of such considerations are as follows [26–28]:

The related boundary conditions are as follows:

The relation for radiative flux is given by . Table 1 and Table 2 provide the thermophysical properties.

The following similarity variables (7) are used to obtain similarity solutions for equations (1)–(5):

The dimensionless set of equation issubject to BCs,where prime “” denotes the differentiation by , is the magnetic number, is Prandtl, is rotation parameter, shows the heat source/sink parameter, and denotes the thermal radiation constant and suction ()/injection () constant.

The dynamic of local Nusselt number and skin friction coefficient is

Substituting (7) into (12) yieldswhere is the local Reynolds number.

#### 3. Stability Analysis

Two branches of solution against the involved parameters are observed. In order to inspect the branch of solution with more physical relevance, the stability analysis has been performed for both branches. The governing equations (2)-(5) with unsteady form are

The unsteady similarity solution can be considered as follows [33–35]:where the is the same as in equation (7) and the dimensionless time is the , where the subscription represents the derivative of the variable subscripted. Remember that continuity (1) is also satisfied by the values and given above. Equation (18) is substituted into equations (14)-(17), and one can getsubject to BCs,

It is assumed that [28]

The solution associated to the formulated equations (8)–(10) with boundary conditions (11) is denoted with , and The decay and growth disturbance rate is defined with It is assumed that , and and derivatives of these expressions are smaller in contrast to the steady state situation. Such assumptions are followed due to the linear stability of flow pattern. Upon replacing equation (26) into set of expressions (19)–(21) with taking , , and are associated to , , and . On this end, the linear eigenvalue problem is formulated aswith

We can solve equations (24)–(26) along with BCs (27) by using the numerical values of , and which were obtained from the solutions of equations (8)–(10) with BCs (11). It is worth noting that the homogeneous linear equations (24)–(26) subject to constants (27) construct an eigenvalue system with eigenvalue . By solving equations (24)–(27), a set of eigenvalues is obtained.

#### 4. Results and Discussion

Following to the investigation of Devi and Devi [33, 34] and Lund et al. [35], here water is considered as base material with copper () and aluminum nanoparticles . Moreover, volume fraction of aluminum is kept constant against the variation of copper volume oxide (). The comparison validation for ensuring the validity task has been done in Table 3 with work of Zaimi et al. [37] with excellent comparison results.

Figures 2–4 are sketched to observe the physical aspect of on and . Results are visualized for suction case. The dual branches of solution are results for various values of suction case. The critical points for are while and are the critical points against and respectively. Two different categories of branches such as dual branch and no branch for suction parameter exist. Moreover, larger change in results increment in which preserve a separation. The dominant values of suction do not allow tension in moving particles, and vorticity is not smothered. In upper branch, with larger suction, the decrement in wall shear force and heat transfer rate is noted. In lower branch, the heat transfer rate is lower but wall shear force enhanced. In upper branch, no significant change in due to suction parameter is inspected while in lower branch, declined.

Figures 5–7 show the variation of and magnetic parameter against copper solid volume fraction . As grows, the critical point of shifts to the left, assisting in the separation of the boundary layer. , and referred to critical points of for , and . As a result, the dual-branch representation is possible for where ; however, no branch exists for . The variation in and get dominant change against in the upper branch for a fixed value of , while the reversal phenomenon can be seen in the bottom branch. Furthermore, as the magnetic field is increased, the skin friction coefficient decreases by maintaining constant values of in the upper branch while slower trend is predicted in the lower branch.

The influences of positive constant on , and is reported in Figures 8 and 9The obtained results are prepared in view of stretching/shrinking parameter . A lower observation offor stretching constant is noted. In upper zone, upon increasing the change in is declining while increases. The thermal boundary layer separation is greater for linear stretching as compared to nonlinear case. Moreover, the heat transfer rate for stretching case is more progressive as compared to shrinking case. From Figure 10, a declining change in is observed for greater

Figures 11–13 show the effect of copper solid volume fraction on velocity profiles and in the - and -directions, as well as a temperature profile . Different profiles in both branches fulfill the far-field BCs for various values of and hence validate the obtained numerical results. In Figure 11, the field of velocity decreases as parameter increases in the upper surface while contradictory results are noted in the lower branch. This is because the volume fraction slows the velocity of the fluid, which reduces the thickness of the momentum layer and, as a result, the velocities in the -direction. On the other hand, when rises, velocity rate increases in the -direction for both branches (refer to Figure 12). Furthermore, as expected, heat transfer rate gets arisen for copper solid volume percentage rises in both branches.

The graphical outcomes for rotation constant parameter on velocity along -directions () have been worked out in Figure 14. The solution has a symmetrical nature. When the value of is positive, velocity increases in both branches. When the value of is negative, the velocity of the hybrid nanofluid decreases. It is also discovered that when , two branches are possible. Figure 15 shows the effect of heat sink/source parameter on temperature profiles. It has been observed that temperature of hybrid nanofluid increases in both solutions when the parameter increases. This enhancement is due to external heat source. This aids in the formation of the thermal boundary layer thickness.

It is well justified fact that upon existing of different solution branches, the importance of stability analysis cannot be denied identifying the more realistically branch. Usually, the upper branch is referred to the physical branch as it meets with the boundary assumptions. However, the stability of lower branch is also importance and can be observed via stability procedure. After carefully noticing, it is claimed that upper branch is more stable in view of The positive sign of ensures that stability of branch and reflects the decay and growth results. The negative values of show that the branch is unstable. Both positive and negative values of upper and lower branches are presented via Table 4. Figure 16 shows the lower numerical values of eigenvalue for upper and lower surface branches.

#### 5. Conclusions

In a rotating hybrid nanofluid, the laminar steady boundary flow on a stretching/shrinking sheet was investigated. The governing PDEs were numerically solved in *bvp4c* using the well-known collocation approach known as the 3-stage-Lobatto-3a method. Numerical results for various values of physical application parameters have been graphically presented and thoroughly analyzed. The values of and are compared to published publications in order to validate the obtained numerical results, which show a positive agreement. In the precise range of the magnetic, suction, and shrinking/stretching parameters, dual branches have been discovered. The following are a few key points from the current study:(i)For suction case, the heat transfer rate is larger as compared to injection phenomenon in the upper branch.(ii)In upper branch, no significance change in y-directional velocity is observed while declining change in velocity is noticed for lower branch.(iii)The presence of external heat source enhanced the heat transfer rate in both branches. However, the increment in heat transfer is relatively larger in the upper branch.(iv)No dual branches exist outside the magnetic parameter and stretching/shrinking constant.(v)Upon increasing the copper nanoparticle volume fraction, the enhancement in temperature rate is measured.(vi)The fluctuation in thermal boundary layer is more progressive for linear stretching instead of nonlinear case.

#### Nomenclature

: | Alumina |

: | Ambient temperature |

: | Copper |

: | Dimensionless velocity |

M: | Hartmann number |

: | Injunction/suction parameter |

: | Local Nusselt number |

: | Local Reynolds number |

B: | Magnetic field |

: | Mass flux of velocity |

Pr: | Prandtl number |

: | Skin friction coefficient |

T :: | Temperature |

: | Sink/source parameter |

: | Variable temperature at the sheet |

u, v, : | Velocity components |

: | Velocity of surface |

: | Water |

: | Angular velocity |

: | Density |

: | Dimensionless temperature |

: | Dynamic viscosity |

: | Electrical conductivity |

: | Heat capacity |

: | Smallest eigenvalue |

: | Stability transformed variable |

: | Stretching/shrinking parameter |

: | Thermal conductivity, |

: | Transformed variable |

: | Unknown eigenvalue. |

#### Data Availability

All data are given in the paper.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.