In this work, we inspect and analyze a two dimensional, unsteady mixed convectional hybrid nanofluid hydromagnetic flow (Al2O3-Cu/H2O) over a convectional heated an extending/contracting surface with the influence of thermal radiation. Hybrid nanofluid (Al2O3-Cu/H2O) flows with magnetohydrodynamic and heat source or sink. Brownian motion and thermophoresis were incorporated using the Buongiorno model. Hybrid nanofluid with vol. fraction range limited to 1.5% and within the higher temperature range of 50°C to 70°C is considered for thermal conductivity and viscosity analysis. The proposed model is then converted into ODEs through similarity transformation with the help of homotopy analysis. The effect of embedded input factors on the temperature, velocity, and concentration profiles is visually demonstrated and explained. The magnetic field has inverse impact on velocity and temperature profiles. Velocity profile increases for both mixed convection and buoyancy ratio parameters. It has been noticed that the temperature profile increases with thermal radiation. For increasing values of Lewis number, the concentration of hybrid nanoparticles is considerably lowered. Moreover, we observed an increase in the concentration of hybrid nanoparticles through a destructive chemical reaction, whereas a generative chemical reaction has the reverse effect. It has been proved that skin friction is increasing function of and decreasing function of . On the other hand, Nusselt number increased with the increase of while Sherwood number is decreased, with the increase of .

1. Introduction

It is eminent that many studies have been carried out on two cases of nanoparticles floating in a base fluid known as “Hybrid Nanofluid” the forefront nanofluid. The key importance of the hybrid nanofluid is that by selecting the right mix of nanoparticles, favorable characteristics may be enhanced, and drawbacks can be mitigated owing to their interactive influence. It is reported that these hybrid nanofluids are new and have a number of application acoustics, defense, manufacturing, transportation, medical, microfluidics, and naval constructions. Nanofluid flow, in particular, is long familiar for its high heat conveyance when equated to regular fluid. The hybrid nanofluid is applied to raise it even further. Jamshed and Aziz [1] studied the Cattaneo-Christov heat flux impact and discovered that spherical shaped nanoparticles transmit heat at the fastest rate when compared to hexagon, platelet type nanoparticles. The authors in [2] examined squeezing flow in a hybrid base fluid with nanoparticles suspended in it. The hybrid nanofluid stagnation point flow past an extending sheet was investigated by the researchers [3]. The influence on flow between two riga-plates was investigated by Ahmed et al. [4]. The authors in [5] studied the impact of changing viscosity on a stretched sheet containing hybrid nanoparticles using the Runge–Kutta fourth order technique. The researchers in [6] interrogated the synthesis of hybrid nanofluid while in [7] examined the effects of a hybrid nanofluid of (Al2O3-Cu/H2O) on heat transmission.

There are numerous numerical methods in the literature to inspect the boundary layer flow and heat transport of a hybrid nanofluid. These concepts were further explored by the researchers in their research [9] to conceptualize the main idea. The problem was then expanded to a three-dimensional flow capable to the Newtonian heating consideration by the authors [10]. In both the above studies, they obtained key results about the heat transfer. The difficulty of a hybrid nanofluid consisting of a rotating flow was described by Hayat and Nadeem [11]. Zainal et al. [12] inspected the unsteady stagnation point flow of a hybrid nanofluid past a convectional heated an extending/contracting surface while accounting for velocity slip’s influence on heat transfer. Using the Buongiorno model, Daniel et al. [13] investigated an unstable blended convectional electrical magnetohydrodynamic (MHD) flow and heat conveyance generated by nanofluid across a permeable stretched surface. The friction coefficient and wear volume were studied by Xie et al. [14] to determine the tribological properties of hybrid nanoparticles. Devi and Anjali [10] investigated the three-dimensional flow of (Cu-Al2O3/water) copper alumina/water hybrid nanofluid using the RK-Fehlberg integration technique.

Thermal radiation’s effect on MHD blood flow and heat conveyance in a permeable capillary in extending motion was inspected. It has been examined that at high temperatures, heat transfer and thermal radiation are recognized to have a significant influence on numerous physiological processes, technological, and engineering industry equipment. Solar collector performance, plume dynamics, rocket propulsion, high-dose cancer therapy, fire propagation, material processing, and combustion systems are just a few examples. With the advancement of analytical and computational techniques, thermal convection flows with high radiative flux have received more attention recently. The effect of heat radiation on peristaltic transport of ionic nanoliquids in biomicrofluid channels had analyzed by Prakash et al. [15]. In context from melting heat transfer and thermal radiation, the researcher [16] proposed the dynamics of stagnation point flow of carbon nanotubes. To get series solution, the optimal homotopy technique was used. Muhammad et al. [17] make another notable attempt to address the key aspects of heat radiation and viscous dissipation impacts in a viscous nanofluid. Mixed convectional and slide impacts in viscous nanofluid across an extended sheet were described by Hsiao [18]. The solution was numerically defined using the convective boundary conditions [19]. The reader can further study about heat transfer through nanofluid flow in refs. [2027]. The relevant study has been seen in [2831].

The main aim of this research is to see how magnetic fields, thermal radiation, and heat generation/absorption affect unstable hybrid nanofluid hydromagnetic flow (Al2O3-Cu/H2O) across a convectional heated an extending/contracting surface. The thermophoresis and Brownian motion characteristics were included in the Buongiorno model. To minimize the autonomous variables in governance equations resulting from mathematical modelling, a suitable collection of dimensionless variables is employed. By implementing the homotopy approach, an analytical solution has been calculated. Comparison between HAM and ND solve has been shown in Table 1.

2. Mathematical Formulation

Assume the unsteady hybrid nanofluid hydromagnetic flow (Al2O3-Cu/H2O) past a convectional heated an extending/contracting surface in the bearing of the magnetic field, chemical reaction, heat generation/absorption, and thermal radiation. The flow problem plot is depicted in Figure 1, where , the extending/contracting is surface velocity, indicates a constant that agree to extending and contracting instances, and denotes the problem of unsteadiness. The velocity of the free stream is notified by , where is the strength of the stagnation flow. and are used to represent the ambient and reference temperature separately. Then, we consider that the bottom of the surface is heated by convectional from a hot fluid at a particular temperature and coefficient of heat transfer has noted by while the mass transfer coefficient is indicated by . The governing boundary layer equations may be recognized as [12, 13] based on all of the assumptions mentioned.

Applying concept of Rosseland approximation is as follows: where the parameters and are explained in nomenclature

Further simplification leads us to the following:

Here, we rewrite Equation (3) as follows: where denotes the factor of velocity in -axis, is the velocity factor in -axis, is the Al2O3-Cu/H2O dynamic viscosity, and and are the Brownian thermophoretic diffusion terms while is the density of Al2O3-Cu/H2O, is the Al2O3-Cu/H2O temperature, is the thermal/heat conductivity of Al2O3-Cu/H2O, and is the Al2O3-Cu/H2O heat capacity. The boundary conditions, as well as the velocity partial slip, are set to

where is the velocity slip term, in which stands for the starting merit of the velocity slip term. Table 2 lists the thermophysical characteristics of copper (Cu), as well as aluminium oxide (Al2O3) and water (H2O) nanoparticles. The thermophysical characteristics of hybrid nanofluid are shown in Table 3. The solid volume fraction of nanoparticles is described as , where denotes the density, is the density of the hybrid nanoparticle, is the constant pressure of heat capacity, represents the thermal conductivity of H2O, and is the thermal conductivity of the hybrid nanoparticle.

The follower similarity transformations are provided in order to represent the governed Equations (1), (2), (7), and (8) concerning the BCs (9) in a much easy way [12]. where is the stream function that can be specified as and and is the similarity variable. Thus, we attain

Equations (2), (7), and (8) are converted into the following set of nonlinear similarity ODEs by using the similarity variables (10) and (11) in light of the aforementioned relationships.

Now the starting and final constraints (9) are also changed into the following: where are unsteadiness parameter, magnetic parameter, thermal radiation parameter, Prandtl number, buoyancy ratio parameter, Biot number, local Reynolds number in -axis, ratio of velocity parameter, heat source/sink parameter, mixed convection parameter, Lewis number, thermophoretic parameter, and Brownian motion parameter.

2.1. Physical Quantities of Interest

For above model local Nusselt number , the skin friction coefficient and Sherwood number are clear as follows: where

In dimensional form, we have from above as

2.2. Idea of HAM for the Model

Here, we used HAM to solve Equations (12)–(14) using boundary conditions (15). The Mathematica software is utilized for this, and the proposed model can be solved using HAM in the following manner: where operators are specified as and ; moreover, we have

We also define the operators, such that

For Equations (12)–(14), the 0th-order system is shown as follows: with the BCs are

While the embedding constraint is , to regulate for the solution convergence , , and are used. At specific values , the following is obtained:

Expand the , , and through Taylor’s series for While BCs are While

3. Outcomes with Discussion

In the following discussion section, the outcomes of various parametric quantity on velocity, concentration profiles, and temperature are shown graphically and discussed. These parameters are unsteadiness parameter, buoyancy ratio parameter, magnetic parameter, thermal radiation parameter, Lewis number, mixed convection parameter, Prandtl number and thermophorotic parameter, and Brownian motion parameter. The geometry of the problem is shown in Figure 1. It can be seen that when the magnetic parameter assigns maximum values, the velocity distribution retains its declining tendency as shown in Figure 2. There is a periodic oscillation velocity with diminishing amplitude while the magnetic force produces a resistive nature force known as the Lorentz force which controls the flow of fluid particles. As a result, to demonstrate the value of reducing the velocity of moving particles, the effect of the unsteadiness parameter on the nanofluid velocity profile has seen in Figure 3. Because of the acceleration situation, brings down rate of fluid flow and a narrower momentum boundary barrier thickness; the behavior occurs. The velocity profile reduces for the higher acceleration.

Figure 4 illustrates the influence of the buoyancy ratio parameter on velocity drawings and also increasing levels of ; there is a noticeable increase in velocity. This increase in velocity drawings is more noticeable in stable condition than in the unsteadness situation. Figure 5 shows the effect of on . From Figure 5, it can be expected that has greater values of . Physically, this is because of the larger values of buoyancy force. The impact of heat generation/absorption on the temperature is plotted in Figure 6. The temperature of the fluid and the thickness of the thermal boundary level are both increased when a heat source is used. The heat sink , on the other hand, lowers the fluid’s temperature and thins the thermal boundary level thickness.

To denote , the lack of heat generation or absorption, the increase in thermal radiation raises the temperature but lowers the concentration profiles, as seen in Figure 7. This is because an increase in thermal radiation gives to a greater extent heat to the hybrid nanofluid, resulting in a rise in temperature and the thickness of the thermal boundary layer. Physically, radiative aspect creates the Brownian motion of minute ingredients faster than normal; thus, random migrated particles strike with one another, and the caused frictional energy transforms it to thermal energy. The effect from the Prandtl number along the temperature profile has seen in Figure 8. By increasing the Prandtl number, the fluid’s thermal conductivity decreases. This fact is due to inverse relation of with thermal diffusivity, and it is well-known fact that the fluid with higher values of has weaker thermal diffusion so that the temperature declines. Figures 9 and 10 demonstrate the consequence of the thermophoresis parametric quantity along the thermal profile. When the thermophoresis parametric quantity is increased which enhance the behavior of the system, thermophoresis effect is associated with movement of nanoparticles from a hot wall to a cold wall, and since it is generated by temperature gradients, it creates a fast flow away from the moving plate; thus, more fluid is heated away from the surface, and this leads to an increase in the temperature within the thermal boundary layer. Figure 10 demonstrates the consequence of the Brownian diffusion parameter on the thermal profile. Physically, the imperfect nature of the Brownian motion parametric quantity enables heating the physical setup. This heating causes nanoparticles to be transferred from the cooler stretched sheet area to the quiescent fluid zone. Because the particles migrate from a high up to a low concentration part, an increase in the Brownian motion parameter and leads decrease in the fluid concentration gradient, as seen in Figure 11.

The Brownian motion is the movement of fluid particles from higher to lower concentration so increases in the thermophoresis parameter produce a huge rise in the concentration dispersion of the fluid flow , which meet to zero at the boundary level, as shown in Figure 12. A minor change in the thermophoresis parameter causes the fluid particles to move quickly, releasing surplus heat energy and causing a huge rise in the concentration dispersion. As a result, with an increase in the value of , Figure 12 indicates a considerable rise in the concentration dispersion. Figure 13 also shows the influence of the Schmidt number on the concentration dispersion. It is noticed that for larger values of , the fluid concentration decreases. This fact is quite similar with that of Prandtl number impact on temperature. Actually, has inverse relation with mass diffusivity; therefore, larger values of correspond thinner concentration boundary layer.

4. Discussion on Tabulated Results

Table 4 shows that is increased when the values of are increased. The is decreased, when the values of are increased. Table 5 shows that is increased when the values of are increased. The is decreased, when the values of are increased. Table 6 shows that is decreased, when the values of are increased. Table 1 shows the very excellent agreement of the HAM and ND solve solution by computer-based package Mathematica 11.1.0.

5. Conclusions

In this research work, the uniform mixed convectional with the combined consequences of thermal radiation and heat generation/absorption along the magnetohydrodynamic (MHD) flow of hybrid nanofluid via a convectional heated the stretching/shrinking surface is studied. Brownian motion and thermophoresis were incorporated using the Buongiorno model. We came at the following key conclusions based on our findings with various factors in this study: (i)Increasing the combined convection and buoyancy ratio parameters improved nanofluid flow properties(ii)It is well established that an increase in the unsteadiness and magnetic parameters result in a decrease in the velocity profile(iii)Consequently, a rise in the thermophorotic parameter, radiation parameter, or Brownian motion parameter means an increase in the temperature profile(iv)It is noticed that the increase of the Prandtl number decreases the temperature profile(v)Increases in the values of the thermophoretic parameter and Brownian motion parameter have an inverse influence on the concentration profile(vi)The displays lessening tendency for rising values of (vii)Skin friction is augmented when the values of are improved. When the values of is enlarged, is reduced(viii)Nusselt number is improved when the values are increased. The is reduced, when the values of are augmented(ix)Sherwood number is declined, when the values of are enlarged


and :Components of velocity along - and -axes (m/s)
:Stretching/shrinking velocity (m/s)
:Ambient temperature (K)
:Heat transfer coefficient
:(kgm-1s-1) Viscosity of hybrid nanofluid
:Reynolds number
:Dimensionless temperature
:The density of the nanoparticle (kgm-3)
:Nanoparticle thermal conductivity (Wm-1K-1)
:Stream function
:Magnetic parameter
:Velocity slip factor
:Ratio of velocity parameter and heat
:Local Nusselt number
:Dimensionless velocity
:Kinematic viscosity (m2s-1)
:Wall mass flux
:Rosseland approximation
:Nanoparticle solid volume fraction
:Coefficient of mean absorption
:Plate concentration
:Concentration of ambient
:Mass transfer coefficient
:Brownian motion parameter
:Cartesian coordinate axis
:Strength of the stagnation flow
:Reference temperature (K)
:Hybrid nanofluid density (kgm-3)
:Hybrid nanofluid volumetric heat capacity (m2s-2K-1)
:Initial value of the velocity slip factor
:Constant pressure of heat capacity
:Base fluid density (kgm-3)
:Base fluid thermal conductivity (Wm-1K-1)
:Unsteadiness parameter
:Dimensionless concentration
:Prandtl number
:Source/sink parameter
:Skin friction coefficient
:Temperature of fluid (K)
:Dynamic viscosity (kgm-1s-1)
:Wall shear stress
:Transportation of heat
:Stefan-Boltzmann constant
:Porosity parameter
:Thermal radiation parameter
:Wall concentration
:Biot number
:Thermophoretic parameter.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.