Next Article in Journal
Transversal Jacobi Operators in Almost Contact Manifolds
Previous Article in Journal
The Existence and Uniqueness of the Solution of a Nonlinear Fredholm–Volterra Integral Equation with Modified Argument via Geraghty Contractions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Nanofluid Slip Flow over an Exponentially Stretching/Shrinking Permeable Sheet with Heat Generation

1
Department of Mathematics, Faculty of Science, Universiti Putra Malaysia, UPM, Serdang 43400, Selangor, Malaysia
2
Institute for Mathematical Research, Universiti Putra Malaysia, UPM, Serdang 43400, Selangor, Malaysia
3
Fakulti Teknologi Kejuruteraan Mekanikal dan Pembuatan, Universiti Teknikal Malaysia Melaka, Hang Tuah Jaya, Durian Tunggal, Melaka 76100, Malaysia
4
Department of Mathematics, Babeş-Bolyai University, R-400084 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
Submission received: 25 November 2020 / Revised: 14 December 2020 / Accepted: 19 December 2020 / Published: 24 December 2020
(This article belongs to the Section Engineering Mathematics)

Abstract

:
An investigation has been done on the hybrid nanofluid slip flow in the existence of heat generation over an exponentially stretching/shrinking permeable sheet. Hybridization of alumina and copper with water as the base fluid is considered. The mathematical model is simplified through the similarity transformation. A numerical solver named bvp4c in Matlab software is utilized to facilitate the problem-solving process and dual solutions are attained. The influences of several pertinent parameters on the main physical quantities of interest and the profiles are scrutinized and presented in the form of graphs. Through the stability analysis, only the first solution is considered as the physical solution. As such, the findings conclude that the upsurges of volume fraction on the copper nanoparticle could enhance the skin friction coefficient and the local Nusselt number.

Graphical Abstract

1. Introduction

Heat transfer becomes the most important process with the evolution of advanced technology. Applications and industries involving heat generation, such as manufacturing, thermal power generation, transportation, chemical processes and many others, necessitate an efficient heat transfer performance to produce the best outcome. Nanofluid applications have been thoroughly reviewed by Sidik et al. [1,2] and also by Huminic and Huminic [3]. Historically, various investigations have been conducted to improvise the heat transfer capability of the transferring medium which is perceived as a heat transfer fluid.
The idea of dispersing high-thermal conductivity solid particles of microsize into the conventional fluid was first pioneered by Maxwell [4] and then continued by Hamilton and Crosser [5] to intensify the fluid thermal conductivity. Nevertheless, some limitations and flaws still occurred, such as the coagulation in the fluid flow passage. Owing to that limitation, a rapid investigation is being made and nanofluid is established. Choi and Eastman [6] invented this novel heat transfer fluid by dispersing the nanometer-sized solid particles into the base fluid, and it is believed these could overcome the coagulation of the flow passage due to its special feature.
However, along with the development of technologies, another class of heat transfer fluid named hybrid nanofluid is established as the extension to nanofluid that is invented by the dispersion of dual or multiple kinds of nanosized solid particles with a good thermal conductivity into a base fluid. Looking through the reviews of the previous investigations by the researchers and scientists [7,8,9,10,11], hybrid nanofluid is claimed to have a preferable capability in terms of thermal conductivity and heat transfer due to its synergistic effects, contrasted to mono nanofluid and regular fluid. Because of that, this new heat transfer fluid captivated many researchers and encouraged them to conduct further investigation into this fluid to understand its rheological properties and the mechanisms behind it, as well as to cope with the exponential growth demand of technology and science.
It seems like Devi and Devi [12,13,14] are among the earliest ones who discovered the flow of hybrid nanofluid at the boundary layer, mathematically and numerically, and their study emphasizes alumina and copper as the hybrid nanoparticles. Aly and Pop [15] considered hybrid nanofluid in their research into the flow and heat transfer over a stretching/shrinking sheet. They concluded that the hybrid nanofluid causes more reduction in the temperature profile unlike the mono nanofluid when the magnetic field, suction and copper volume fraction parameters are amplified. Waini et al. [16] also reported that the presence of a huge volume fraction of copper intensifies the skin friction coefficient and the local Nusselt number. Many other recent interesting investigations towards hybrid nanofluid, focusing on alumina and copper as the hybrid nanoparticles, with various configurations of parameters, have also been revealed by Yashkun et al. [17], Lund et al. [18,19], Aladdin et al. [20], Roşca et al. [21], Khashi’ie et al. [22,23,24,25], Roy and Pop [26], Zainal et al. [27], Anuar et al. [28], Gangadhar et al. [29] and Wahid et al. [30].
Slip condition is one of the parameters that is extensively being considered among the researchers and is significantly important in the heat transfer process. This condition happens when the fluid flow velocity is different from the velocity at the boundary. In industrial application, slip is significantly used in the microelectronic cooling system, micro heat exchangers, internal cavities’ and artificial heart valves’ polishing, and drug delivery system [31,32,33]. The presence of slip could reduce the resistance of the flow in microchannels and nanochannels as well as enhancing the efficiency of microscale viscous pumps. The appropriate degree of slip should be considered for the effectiveness of the applications that are also concerned with the other related parameters as well, thus various investigations towards the slip condition should be examined specifically. Presumably, Andersson [34] was the first one who examined the slip effect of the flow in the boundary layer of viscous fluid. In this regard, the investigation on the slip effect towards the hybrid nanofluid flow is still an ongoing research area. More recently, Hayat et al. [35] considered the partial slip on the Ag-CuO/H2O hybrid nanofluid due to rotating flow. They found out that the skin friction coefficient is lessened as the velocity slip is augmented in the case of the stretching sheet. Nadeem and Abbas [36] scrutinized the micropolar hybrid nanofluid with slip effect. A comparative investigation was made by Hayat et al. [37] on the hybrid and mono nanofluid of multi-walled carbon nanotubes with the slip flow. Yan et al. [38] elucidated the multiple slip conditions on the magnetized hybrid nanofluid with the consideration of a few other parameters. It was found out that thermal slip diminishes the heat transfer rate. Later, the velocity slip and temperature slip were simultaneously modeled by Alharbi [39] in the study of hybrid nanofluid and the influence of velocity slip is consistent with the result obtained by Hayat et al. [35]. Such other studies on slip conditions towards hybrid nanofluid were also reported by Tlili et al. [40], Aly and Pop [41] and Eid and Nafe [42]. The reported results in these investigations could provide the initial insight for the engineers and researchers needed to predict and control the slip parameter towards the related applied area.
Besides this, heat generation/absorption also becomes an essential aspect when managing the performance of the heat transfer. Due to that, investigations into the effects of heat generation/absorption towards nanofluid have been carried out up until now. Abbas et al. [43] contemplated the slip flow of nanofluids with heat generation due to the curved stretching surface. Upreti et al. [44] explored the heat generation of nanofluids due to the flat porous plate. Jusoh et al. [45] elaborated on the nanofluid stagnation point flow with heat generation. Murugesan and Kumar [46] also elucidated the same effect but with the flow of nanofluid over the exponentially stretching sheet. With the evolution of nanofluid to hybrid nanofluid, this heat generation/absorption effect has been considered by Hayat and Nadeem [47] with regard to the model of hybrid nanofluid with the three-dimensional flow, which was then continued by Hayat et al. [48], who analyzed such effects on the rotating hybrid nanofluid. Eid and Nafe [42] also considered the presence of heat generation together with the slip effect in their research on hybrid nanofluid. Most of them indicated that the amplification of heat generation increases the temperature profile.
Therefore, in this present study, we intended to scrutinize the slip flow and heat transfer of a hybrid nanofluid of alumina and copper (hybrid nanoparticles) with water (base fluid) over an exponentially permeable stretching/shrinking sheet with the consideration of heat generation. Both velocity and thermal slips are being considered in the model. We believe the model of a hybrid nanofluid with such configurated effects and surface is still not being explored by the researchers, judging by the previous literature. The influences of various parameters over the main physical quantities are provided in the figure of the graph, and are discussed. An analysis of flow stability is also conducted, as the non-unique solutions are visible and the least eigenvalues that imply the stability of the flow are tabulated.

2. Mathematical Formulation

A steady boundary layer of hybrid nanofluid slip flow induced by an exponentially permeable stretching/shrinking sheet with heat generation is scrutinized. According to the physical model in Figure 1, u and v are the elements of velocity in x and y directions, respectively. The surface velocity is given as u w ( x ) = c e x / L and the wall mass transfer velocity is specified as v w ( x ) = v 0 e x / 2 L ( v 0 < 0 is the mass suction and v 0 > 0 is the mass injection), λ > 0 is the stretching constant, λ < 0 is the shrinking constant, λ = 0 refers to a motionless sheet, T is temperature, T w = T + T 0 e x / 2 L is the sheet of varying temperature with constant T 0 and q = q 0 e x / L is the heat generation rate constant.
Therefore, concerning the mentioned boundary layer assumptions, the equations of continuity, momentum and energy are formulated as (Murugesan and Kumar [46], Ghosh and Mukhopadhyay [49], Waini et al. [16], Devi and Devi [14], Eid and Nafe [42])
u x + v y = 0
u u x + v v y = μ h n f ρ h n f 2 u y 2
u T x + v T y = k h n f ( ρ C p ) h n f 2 T y 2 + q ( ρ C p ) h n f ( T T )
due to the boundary conditions (Ghosh and Mukhopadhyay [49], Mukhopadhyay and Andersson [50])
u = u w ( x ) λ + A 1 μ h n f ρ h n f u y ,   v = v w ,   T = T w ( x ) + B 1 T y   at   y = 0 u 0 ,   T T   as   y
in which A 1 and B 1 are the slip factors for velocity and thermal, respectively. The correlation of the physical properties (Tayebi and Chamkha [51], Takabi and Salehi [52], Babu et al. [11], Ghalambaz et al. [53]) and the hybrid nanofluid’s thermo-physical properties (Oztop and Abu-Nada [54]) is depicted in Table 1 and Table 2, respectively. Here, the hnf, f, s1 and s2 subscripts refer to the hybrid nanofluid, base fluid and solid nanoparticle for the alumina and solid nanoparticle for copper, accordingly, and ϕ s 1 and ϕ s 2 symbolize the alumina and copper volume fraction parameters, separately.
As to reduce the aforementioned governing equations, the similarity variables are employed (Waini et al. [16], Eid and Nafe [42])
ψ = e x / 2 L 2 ν f L c f ( η ) ,   u = ψ y ,   v = ψ x ,   θ ( η ) = T T T w T ,   η = y e x / 2 L c 2 ν f L
which consequently gratify Equation (1) and transform Equations (2) and (3) into the subsequent equations
( μ h n f / μ f ρ h n f / ρ f ) f + f f 2 f 2 = 0
1 Pr ( k h n f k f ) θ + ( ρ C p ) h n f ( ρ C p ) f ( f θ f θ ) + β θ = 0
due to the transformed boundary conditions
f ( 0 ) = S ,   f ( 0 ) = λ + A f ( 0 ) ,   θ ( 0 ) = 1 + B θ ( 0 ) f ( η ) 0 ,   θ ( η ) 0   as   η
such that Pr = ( μ C p ) f / k f is the Prandtl number, β = 2 q 0 L / c ( ρ C p ) f is the heat generation parameter, A = A 1 ( μ h n f / ρ h n f ) e x / 2 L c / 2 ν f L is the dimensionless velocity slip parameter and B = B 1 e x / 2 L c / 2 ν f L is the dimensionless thermal slip parameter, and S = ν 0 / ν f c / 2 L is the wall mass transfer parameter in which S > 0   ( ν 0 < 0 ) equates to mass suction and S < 0   ( ν 0 > 0 ) equates to mass injection.
The main physical quantities of interest (local skin friction coefficient C f and the local Nusselt number N u x ) are mathematically computed as (Waini et al. [16], Eid and Nafe [42])
C f = μ h n f ρ f 1 u 2 w ( u y ) y = 0 ,   N u x = k h n f k f 2 L ( T w T ) ( T y ) y = 0
or by introducing the transform, we have
Re x 1 / 2 C f = ( μ h n f μ f ) f ( 0 ) ,   Re x 1 / 2 N u x = ( k h n f k f ) θ ( 0 )
where Re x = 2 L u w / ν f is the local Reynolds number.

3. Flow Stability

As it is attainable to compute dual solutions, the flow stability analysis is executed. The problem is regarded as unsteady in this analysis where the continuity equation (Equation (1)) is held, and the time derivative is introduced to Equations (2) and (3)
u u x + v u y + u t = μ h n f ρ h n f 2 u y 2
u T x + v T y + T t = k h n f ( ρ C p ) h n f 2 T y 2 + q ( ρ C p ) h n f ( T T )
where the variable t here denotes time.
Referencing the work reported by Merkin [55], Weidman et al. [56], Waini et al. [16] and Yan et al. [38], the dimensionless time-dependent variable transformation for Equations (11) and (12) is
η = y e x / 2 L c 2 ν f L ,   ψ = f ( η , τ )   e x / 2 L 2 ν f L c ,   θ ( η , τ ) = T T T w T , τ = c 2 L t e x / L
As such, with the utilization of the transformation, the subsequent equations are procured
( μ h n f / μ f ρ h n f / ρ f ) 3 f η 3 + 2 f η 2 f 2 ( f η ) 2 2 f η τ = 0
1 Pr ( k h n f k f ) 2 θ η 2 + ( ρ C p ) h n f ( ρ C p ) f ( θ η f θ f η ( θ τ ) ) + β θ = 0
restricted to the conditions
f ( 0 , τ ) = S ,   f η ( 0 , τ ) = λ + A 2 f η 2 ( 0 , τ ) ,   θ ( 0 ) = 1 + B θ η ( 0 , τ ) f η ( η , τ ) 0 ,   θ ( η , τ ) 0   as   η
The succeeding perturbation equation is introduced for the stability analysis of the similarity solutions, f ( η ) = f 0 ( η ) and θ ( η ) = θ 0 ( η ) ,
f ( η , τ ) = f 0 ( η ) + e γ τ F ( η ) θ ( η , τ ) = θ 0 ( η ) + e γ τ G ( η )
where γ is an unidentified eigenvalue, while F ( η ) and G ( η ) are small relative to f 0 ( η ) and θ 0 ( η ) , correspondingly (Weidman et al. [56], Waini et al. [16], Yan et al. [38]).
Substituting Equation (17) into Equations (14)–(16) and setting the value of τ to zero, the initial decay or growth of disturbance in Equation (17) can be identified and the subsequent linearized equations are attained
( μ h n f / μ f ρ h n f / ρ f ) F + f 0 F + F f 0 4 f 0 F + γ F = 0
1 Pr ( k h n f k f ) G + ( ρ C p ) h n f ( ρ C p ) f ( θ 0 F + G f 0 θ 0 F G f 0 + γ G ) + β G = 0
subject to the linearized conditions
F ( 0 ) = 0 ,   F ( 0 ) = A F ( 0 ) ,   G ( 0 ) = B G ( 0 ) F ( ) 0 ,   G ( ) 0
Since most of the linearized conditions (Equation (20)) are equal to zero, it is essential for F ( ) 0 to relax and be replaced with F ( 0 ) = 1 to produce a realizable range of eigenvalues. Thus, to attain the unknown eigenvalues γ which act as the indicator to ascertain the stability of the flow, the bvp4c solver is utilized. The flow is contemplated to be stable when γ > 0 , and vice versa.

4. Result and Discussions

The solutions for Equations (6) and (7) conditioned to Equation (8) are numerically enumerated by employing the bvp4c solver in Matlab software. The suitable preliminary guesses, preferable boundary layer thickness and various values of parameters need to be selected and well-adjusted in the coding function of the solver to compute the most precise results. Stretching ( λ > 0 ) and shrinking ( λ < 0 ) surfaces are considered in this study. The Prandtl number and nanoparticle volume fraction for alumina are specified to be fixed throughout this study as Pr = 6.2 and ϕ s 1 = 0.01 , respectively, while the other parameters, such as velocity A and thermal B slips, nanoparticle volume fraction for copper ϕ s 2 , heat generation β and suction S , are set to be varied in order to study their impacts with regard to the boundary layer flow and heat transfer.
For verification purposes, the results presented in this study are contrasted with the past results revealed by Ghosh and Mukhopadhyay [49], Hafidzuddin et al. [57] and Waini et al. [16]. The numerical values of f ( 0 ) and θ ( 0 ) are computed and compared in Table 3 for the case of viscous fluid ( ϕ s 1 = ϕ s 2 = 0 ) and shrinking surface ( λ = 1 ) with S = 3 and Pr = 0.7 . We observe that the present and previous results are found to be in a reasonable correlation. Thus, we guarantee that the technique, as well as the results provided in this present study, are valid and acceptable.
The non-unique or the dual solutions of Equations (6) and (7) are perceived to exist in this study for certain values of parameters, as can be observed in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9 for the plots of skin friction coefficient Re x 1 / 2 C f and the local Nusselt number Re x 1 / 2 N u x , and also in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17 for the plots of velocity and temperature profiles. The plots of the skin friction coefficient Re x 1 / 2 C f and local Nusselt number Re x 1 / 2 N u x against the stretching/shrinking parameter λ , with varied values of velocity slip parameter A   ( = 0.2 ,   0.4 ,   0.6 ) when S = 2.4 ,   ϕ s 2 = 0.01 ,   B = 0.2 and β = 0.02 , are exhibited in Figure 2 and Figure 3. It is discovered from the figures that the presence of dual solutions is visible for A = 0.2 ,   0.4 ,   0.6 when λ c 1 = 1.4419 ,   λ c 2 = 1.7322 ,   λ c 3 = 2.0362 respectively. This also indicates that the bifurcation point of the boundary layer only appears at the shrinking surface ( λ < 0 ) region and the surges of velocity slip parameter A delay the bifurcation of the boundary layer.
In Figure 2, the skin friction coefficient Re x 1 / 2 C f is noted to be in a decreasing trend as the velocity slip parameter A escalates for the first solution in a certain range ( 1.4419 λ < 0 ) in the shrinking surface ( λ < 0 ) region, but the opposite trend occurs when approaching the stretching surface ( λ > 0 ) region and onwards. However, for the local Nusselt number Re x 1 / 2 N u x , it is noticed from Figure 3 that the enhancement in the value of the velocity slip parameter A promotes the increment in the local Nusselt number Re x 1 / 2 N u x for the first solution, but when approaching the stretching surface ( λ > 0 ) region, the opposite trend occurs. Meanwhile, for the second solution of the skin friction coefficient Re x 1 / 2 C f and the local Nusselt number Re x 1 / 2 N u x , a similar trend is also noticed to occur like the first solution, but a change of trend happens within a certain range of shrinking surface ( λ < 0 ) region and the stretching surface region ( λ > 0 ) onwards, unlike in the first solution. The skin friction coefficient Re x 1 / 2 C f is also perceived to increase as the value of the stretching/shrinking parameter λ diminishes, while the obverse occurs for the local Nusselt number Re x 1 / 2 N u x for both the first and second solutions.
The effects of the thermal slip and heat generation parameters only impact the local Nusselt number Re x 1 / 2 N u x . The plot of the local Nusselt number Re x 1 / 2 N u x against the stretching/shrinking parameter λ for varied values of thermal slip B   ( = 0.2 ,   0.22 ,   0.23 ) and heat generation β   ( = 0.02 ,   0.2 ,   0.3 ) when S = 2.4 ,   ϕ s 2 = 0.01 and A = 0.2 is displayed in Figure 4 and Figure 5. Dual solutions are observed to exist, starting at almost the same critical values for the varied values of both the thermal slip B and heat generation β parameters, which are λ c = 1.4419 . The existence of the thermal slip B and heat generation β parameters causes the local Nusselt number Re x 1 / 2 N u x to decrease, but it is seen to increase as the stretching/shrinking parameter rises λ for both the first and second solutions.
Moreover, for the influence of copper volume fraction ϕ s 2   ( = 0.001 ,   0.005 ,   0.01 ) over the skin friction coefficient Re x 1 / 2 C f , and the local Nusselt number Re x 1 / 2 N u x against the stretching/shrinking parameter λ when S = 2.4 ,   A = B = 0.2 ,   β = 0.02 , is shown in Figure 6 and Figure 7, and the same against the suction parameter S when λ = 1 ,   A = B = 0.2 ,   β = 0.02 is shown in Figure 8 and Figure 9. It is noted that in Figure 6 and Figure 7, the existence of dual solutions is visible for ϕ s 2   ( = 0.001 ,   0.005 ,   0.01 ) when λ c 1 = 1.3686 ,   λ c 2 = 1.4016 ,   λ c 3 = 1.4419 , which indicates that the bifurcation point of the boundary layer eventuates at the shrinking surface ( λ < 0 ) region and the presence of the copper volume fraction delays the bifurcation of the boundary layer.
The upsurge in the copper volume fraction parameter ϕ s 2 leads to the increase in the skin friction coefficient Re x 1 / 2 C f at the shrinking surface ( λ < 0 ) region for the first solution, but leads to a reduction for the second solution. We also notice a change in the trend of the skin friction coefficient Re x 1 / 2 C f for the first solution towards the copper volume fraction parameter ϕ s 2 when λ > 0 . Meanwhile, the plot of the local Nusselt number Re x 1 / 2 N u x is observed to rise when there is an upsurge in the copper volume fraction parameter ϕ s 2 for the first and second solutions. The plots of the skin friction coefficient Re x 1 / 2 C f and local Nusselt number Re x 1 / 2 N u x are noticed to increase and decrease, respectively, when the stretching/shrinking parameter λ decreases for both solutions.
In Figure 8 and Figure 9, dual solutions are perceived to exist for ϕ s 2   ( = 0.001 ,   0.005 ,   0.01 ) when S c 1 = 2.0775 ,   S c 2 = 2.0551 ,   S c 3 = 2.0290 , which indicates that only certain appropriate values of suction are required to achieve dual solutions for the shrinking surface. In these figures, the increases in copper volume fraction ϕ s 2 and suction S parameters augment the skin friction coefficient Re x 1 / 2 C f for the first solution but reduce the skin friction coefficient Re x 1 / 2 C f for the second solution. Concurrently, the local Nusselt number Re x 1 / 2 N u x is noted to augment when the copper volume fraction ϕ s 2 and suction S parameters increase for both of the solutions. From all of the figures (Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8 and Figure 9), we noted that the plots of the skin friction coefficient and the local Nusselt number against both the stretching/shrinking and the suction parameters are all moving towards the left side of the domain, either up or down on the left side, to achieve the critical point.
The velocity f ( η ) and temperature θ ( η ) profiles are also provided as in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17. The increment of velocity slip A   ( = 0.2 ,   0.4 ,   0.6 ) and suction S   ( = 2.2 ,   2.3 ,   2.4 ) parameters causes the velocity profile f ( η ) to increase for the first solution, but it reduces for the second solution and for both solutions of the temperature profile θ ( η ) as depicted in Figure 10, Figure 11, Figure 12, Figure 13, separately. For the first realizable solution, physically, the presence of velocity slip at the boundary of the shrinking surface reduces the flow resistance which causes more flow to slip through the surface and increase the flow velocity, thus a high degree of velocity slip is required to speed up the flow velocity and vice versa. For the copper volume fraction parameter ϕ s 2   ( = 0.001 ,   0.005 ,   0.01 ) , the upsurge in value amplifies both the velocity f ( η ) and temperature θ ( η ) profiles for both of the solutions, as exemplified in Figure 14 and Figure 15, but in the velocity profile f ( η ) for the second solution, dissimilar behavior is noticed at a certain range. In Figure 16 and Figure 17, the influences of the thermal slip parameter B   ( = 0.2 ,   0.4 ,   0.6 ) and heat generation β   ( = 0.02 ,   0.04 ,   0.06 ) parameter over the temperature profile θ ( η ) are portrayed. Since these parameters do not affect the velocity profile f ( η ) , it is only necessary to provide the temperature profile θ ( η ) in this present study. The temperature profile θ ( η ) is noticed to increase when the thermal slip parameter B decreases and when the heat generation parameter β increases. The appropriate amounts of thermal slip and heat generation parameters need to be controlled to obtain the required heat transfer performance, as both parameters produce a different impact that could offset each other.
As dual solutions exist, it is essential to execute the stability analysis of the flow. The stability analysis is conducted as has been discussed in the previous section. The computations of the lowest eigenvalues γ for A = 0.2 and A = 0.6 for different λ when S = 2.4 ,   ϕ s 2 = 0.01 ,   B = 0.2 and β = 0.02 are tabulated in Table 4. It is realized that the first solution gives the positive value of eigenvalues γ 1 , meanwhile the second solution gives the negative value of eigenvalues γ 2 . As according to Merkin [55] and Weidman et al. [56], the positive eigenvalues γ 1 imply that the flow is real and stable due to the initial decay of disturbance as the time passes, and the negative eigenvalues γ 2 imply that the flow is unreal and unstable due to the initial growth of disturbance. It is also noticeable from the table that as the value of λ λ c , the value of γ 0 . So, as supported by the statement above, it can be affirmed that the flow is stable for the first solution whilst it is unstable for the second solution.

5. Conclusions

The slip flow of a hybrid nanofluid over a permeable exponentially stretching/shrinking sheet with the existence of heat generation has been numerically scrutinized by using the bvp4c solver. The present findings can be summarized as follows:
  • Dual solutions are achievable at the shrinking surface region as well as with the appropriate amount of suction, and only the first solution is stable.
  • The increment of the velocity slip parameter depreciates the skin friction coefficient but enhances the local Nusselt number in the shrinking surface region (first solution).
  • The skin friction coefficient is intensified (first solution) when the values of the copper volume fraction and suction parameters increase in the shrinking surface region.
  • The local Nusselt number is augmented (first and second solutions) when the copper volume fraction and suction parameters increase, but is diminished (first and second solutions) as the thermal slip and heat generation parameters increase.
  • The rise in the stretching/shrinking parameter reduces the skin friction coefficient but augments the local Nusselt number.
  • The augmentation of the velocity slip, suction and copper volume fraction parameters augments the velocity profile (first solution).
  • The temperature profile can be enhanced (first and second solutions) with the enhancement of the copper volume fraction and heat generation parameters, but can be diminished (first and second solutions) with the increment of the velocity and thermal slips and suction parameters.

Author Contributions

Research design, N.S.W., N.M.A., N.S.K. and I.P.; formulation and methodology, N.S.W., I.P.; result analysis, N.S.W.; validation, N.M.A. and I.P.; article preparation, N.S.W.; review and editing, N.S.W., N.M.A., N.S.K. and I.P. All authors have read and agreed to the published version of the manuscript.

Funding

Funded by Universiti Putra Malaysia.

Acknowledgments

The authors appreciatively acknowledge the support received from Ministry of Education (Malaysia) through the Fundamental Research Grant Scheme (FRGS)-5540309, Universiti Putra Malaysia and Universiti Teknikal Malaysia Melaka.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

u , v elements of velocity in x and y directions, respectively
u w surface velocity
v w wall mass transfer velocity
L characteristic length
λ stretching/shrinking parameter
T fluid temperature
T w varying temperature of the surface
T reference temperature
q heat generation rate constant
A 1 velocity slip factor
B 1 thermal slip factor
ρ f fluid density
ρ h n f fluid density for hybrid nanofluid
( ρ C p ) f heat capacity for fluid
( ρ C p ) h n f heat capacity for hybrid nanofluid
μ f dynamics viscosity for fluid
μ h n f dynamics viscosity for hybrid nanofluid
k f thermal conductivity for fluid
k h n f thermal conductivity for hybrid nanofluid
ϕ s 1 , ϕ s 2 nanoparticles volume fraction for alumina and copper, respectively
f dimensionless velocity
θ dimensionless temperature
Pr Prandtl number
β heat generation parameter
A dimensionless velocity slip parameter
B dimensionless thermal slip parameter
S wall mass transfer parameter
C f local skin friction coefficient
N u x local Nusselt number
Re x local Reynolds number
t time (s)
τ dimensionless time variable
γ eigenvalue
γ 1 , γ 2 least eigenvalue for first and second solutions, respectively

References

  1. Sidik, N.A.C.; Adamu, I.M.; Jamil, M.M.; Kefayati, G.H.R.; Mamat, R.; Najafi, G. Recent Progress on Hybrid Nanofluids in Heat Transfer Applications: A Comprehensive Review. Int. Commun. Heat Mass Transf. 2016, 78, 68–79. [Google Scholar] [CrossRef]
  2. Sidik, N.A.C.; Jamil, M.M.; Japar, W.M.A.A.; Adamu, I.M. A Review on Preparation Methods, Stability and Applications of Hybrid Nanofluids. Renew. Sustain. Energy Rev. 2017, 80, 1112–1122. [Google Scholar] [CrossRef]
  3. Huminic, G.; Huminic, A. Hybrid Nanofluids for Heat Transfer Applications—A State-of-the-Art Review. Int. J. Heat Mass Transf. 2018, 125, 82–103. [Google Scholar] [CrossRef]
  4. Maxwell, J.C. A Treatise on Electricity and Magnetism; Clarendon Press: Oxford, UK, 1873. [Google Scholar]
  5. Hamilton, R.L.; Crosser, O.K. Thermal Conductivity of Heterogeneous Two-Component Systems. Ind. Eng. Chem. Fund. 1962, 1, 187–191. [Google Scholar] [CrossRef]
  6. Choi, S.U.S.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles. ASME Fluids Eng. Div. 1995, 231, 99–106. [Google Scholar]
  7. Sarkar, J.; Ghosh, P.; Adil, A. A Review on Hybrid Nanofluids: Recent Research, Development and Applications. Renew. Sustain. Energy Rev. 2015, 43, 164–177. [Google Scholar] [CrossRef]
  8. Sundar, L.S.; Sharma, K.V.; Singh, M.K.; Sousa, A.C.M. Hybrid Nanofluids Preparation, Thermal Properties, Heat Transfer and Friction Factor—A Review. Renew. Sustain. Energy Rev. 2017, 68, 185–198. [Google Scholar] [CrossRef]
  9. Nabil, M.F.; Azmi, W.H.; Hamid, K.A.; Zawawi, N.N.M.; Priyandoko, G.; Mamat, R. Thermo-Physical Properties of Hybrid Nanofluids and Hybrid Nanolubricants: A Comprehensive Review on Performance. Int. Commun. Heat Mass Transf. 2017, 83, 30–39. [Google Scholar] [CrossRef]
  10. Sajid, M.U.; Ali, H.M. Thermal Conductivity of Hybrid Nanofluids: A Critical Review. Int. J. Heat Mass Transf. 2018, 126, 211–234. [Google Scholar] [CrossRef]
  11. Babu, J.R.; Kumar, K.K.; Rao, S.S. State-of-Art Review on Hybrid Nanofluids. Renew. Sustain. Energy Rev. 2017, 77, 551–565. [Google Scholar] [CrossRef]
  12. Devi, S.P.A.; Devi, S.S.U. Numerical Investigation of Hydromagnetic Hybrid Cu–Al2O3/Water Nanofluid Flow over a Permeable Stretching Sheet with Suction. Int. J. Nonlinear Sci. Numer. Simul. 2016, 17, 249–257. [Google Scholar] [CrossRef]
  13. Devi, S.S.U.; Devi, S.P.A. Numerical Investigation of Three-Dimensional Hybrid Cu–Al2O3/Water Nanofluid Flow over a Stretching Sheet with Effecting Lorentz Force Subject to Newtonian Heating. Can. J. Phys. 2016, 94, 490–496. [Google Scholar] [CrossRef]
  14. Devi, S.U.; Devi, S.A. Heat Transfer Enhancement of Cu-Al2O3/Water Hybrid Nanofluid Flow over a Stretching Sheet. J. Niger. Math. Soc. 2017, 36, 419–433. [Google Scholar]
  15. Aly, E.H.; Pop, I. MHD Flow and Heat Transfer over a Permeable Stretching/Shrinking Sheet in a Hybrid Nanofluid with a Convective Boundary Condition. HFF 2019, 29, 3012–3038. [Google Scholar] [CrossRef]
  16. Waini, I.; Ishak, A.; Pop, I. Hybrid Nanofluid Flow Induced by an Exponentially Shrinking Sheet. Chin. J. Phys. 2019. [Google Scholar] [CrossRef]
  17. Yashkun, U.; Zaimi, K.; Abu Bakar, N.A.; Ishak, A.; Pop, I. MHD Hybrid Nanofluid Flow over a Permeable Stretching/Shrinking Sheet with Thermal Radiation Effect. HFF 2020. [Google Scholar] [CrossRef]
  18. Lund, L.A.; Omar, Z.; Khan, I.; Seikh, A.H.; Sherif, E.-S.M.; Nisar, K.S. Stability Analysis and Multiple Solution of Cu–Al2O3H2O Nanofluid Contains Hybrid Nanomaterials over a Shrinking Surface in the Presence of Viscous Dissipation. J. Mater. Res. Technol. 2020, 9, 421–432. [Google Scholar] [CrossRef]
  19. Lund, L.A.; Omar, Z.; Khan, I.; Sherif, E.-S.M. Dual Solutions and Stability Analysis of a Hybrid Nanofluid over a Stretching/Shrinking Sheet Executing MHD Flow. Symmetry 2020, 12, 276. [Google Scholar] [CrossRef] [Green Version]
  20. Aladdin, N.A.L.; Bachok, N.; Pop, I. Cu-Al2O3/Water Hybrid Nanofluid Flow over a Permeable Moving Surface in Presence of Hydromagnetic and Suction Effects. Alex. Eng. J. 2020, 59, 657–666. [Google Scholar] [CrossRef]
  21. Roşca, N.C.; Roşca, A.V.; Jafarimoghaddam, A.; Pop, I. Cross Flow and Heat Transfer Past a Permeable Stretching/Shrinking Sheet in a Hybrid Nanofluid. HFF 2020. [Google Scholar] [CrossRef]
  22. Khashi’ie, N.S.; Md Arifin, N.; Pop, I.; Nazar, R. Melting Heat Transfer in Hybrid Nanofluid Flow along a Moving Surface. J. Therm. Anal. Calorim. 2020. [Google Scholar] [CrossRef]
  23. Khashi’ie, N.S.; Arifin, N.M.; Pop, I.; Nazar, R.; Hafidzuddin, E.H.; Wahi, N. Flow and Heat Transfer Past a Permeable Power-Law Deformable Plate with Orthogonal Shear in a Hybrid Nanofluid. Alex. Eng. J. 2020, 59, 1869–1879. [Google Scholar] [CrossRef]
  24. Khashi’ie, N.S.; Arifin, N.M.; Pop, I.; Wahid, N.S. Flow and Heat Transfer of Hybrid Nanofluid over a Permeable Shrinking Cylinder with Joule Heating: A Comparative Analysis. Alex. Eng. J. 2020, 59, 1787–1798. [Google Scholar] [CrossRef]
  25. Khashi’ie, N.S.; Md Arifin, N.; Pop, I. Mixed Convective Stagnation Point Flow towards a Vertical Riga Plate in Hybrid Cu-Al2O3/Water Nanofluid. Mathematics 2020, 8, 912. [Google Scholar] [CrossRef]
  26. Roy, N.C.; Pop, I. Flow and Heat Transfer of a Second-Grade Hybrid Nanofluid over a Permeable Stretching/Shrinking Sheet. Eur. Phys. J. Plus 2020, 135, 768. [Google Scholar] [CrossRef]
  27. Zainal, N.A.; Nazar, R.; Naganthran, K.; Pop, I. MHD Flow and Heat Transfer of Hybrid Nanofluid over a Permeable Moving Surface in the Presence of Thermal Radiation. HFF 2020. [Google Scholar] [CrossRef]
  28. Anuar, N.S.; Bachok, N.; Pop, I. Radiative Hybrid Nanofluid Flow Past a Rotating Permeable Stretching/Shrinking Sheet. HFF 2020. [Google Scholar] [CrossRef]
  29. Gangadhar, K.; Bhargavi, D.N.; Kannan, T.; Venkata Subba Rao, M.; Chamkha, A.J. Transverse MHD Flow of Al2O3–Cu/H2O Hybrid Nanofluid with Active Radiation: A Novel Hybrid Model. Math. Methods Appl. Sci. 2020, 1–19. [Google Scholar] [CrossRef]
  30. Wahid, N.S.; Md Arifin, N.; Turkyilmazoglu, M.; Hafidzuddin, M.E.H.; Abd Rahmin, N.A. MHD Hybrid Cu-Al2O3/Water Nanofluid Flow with Thermal Radiation and Partial Slip Past a Permeable Stretching Surface: Analytical Solution. JNanoR 2020, 64, 75–91. [Google Scholar] [CrossRef]
  31. Bataineh, K.M.; Al-Nimr, M.A. 2D Navier–Stokes Simulations of Microscale Viscous Pump With Slip Flow. J. Fluids Eng. 2009, 131, 051105. [Google Scholar] [CrossRef]
  32. Sharatchandra, M.C.; Sen, M.; Gad-el-Hak, M. Thermal Aspects of a Novel Viscous Pump. J. Heat Transf. 1998, 120, 99–107. [Google Scholar] [CrossRef]
  33. Nandal, J.; Kumari, S.; Rathee, R. The Effect of Slip Velocity on Unsteady Peristalsis MHD Blood Flow through a Constricted Artery Experiencing Body Acceleration. Int. J. Appl. Mech. Eng. 2019, 24, 645–659. [Google Scholar] [CrossRef] [Green Version]
  34. Andersson, H.I. Slip Flow Past a Stretching Surface. Acta Mechanica 2002, 158, 121–125. [Google Scholar] [CrossRef]
  35. Hayat, T.; Nadeem, S.; Khan, A.U. Rotating Flow of Ag-CuO/H2O Hybrid Nanofluid with Radiation and Partial Slip Boundary Effects. Eur. Phys. J. E 2018, 41, 75. [Google Scholar] [CrossRef]
  36. Nadeem, S.; Abbas, N. On Both MHD and Slip Effect in Micropolar Hybrid Nanofluid Past a Circular Cylinder under Stagnation Point Region. Can. J. Phys. 2019, 97, 392–399. [Google Scholar] [CrossRef]
  37. Muhammad, K.; Hayat, T.; Alsaedi, A.; Ahmad, B.; Momani, S. Mixed Convective Slip Flow of Hybrid Nanofluid (MWCNTs + Cu + Water), Nanofluid (MWCNTs + Water) and Base Fluid (Water): A Comparative Investigation. J. Therm. Anal. Calorim. 2020. [Google Scholar] [CrossRef]
  38. Yan, L.; Dero, S.; Khan, I.; Mari, I.A.; Baleanu, D.; Nisar, K.S.; Sherif, E.-S.M.; Abdo, H.S. Dual Solutions and Stability Analysis of Magnetized Hybrid Nanofluid with Joule Heating and Multiple Slip Conditions. Processes 2020, 8, 332. [Google Scholar] [CrossRef] [Green Version]
  39. Alharbi, S.O. Influence of Wall Slip and Jump in Wall Temperature on Transport of Heat Energy in Hybrid Nanofluid. J. Therm. Anal. Calorim. 2020, 1–8. [Google Scholar] [CrossRef]
  40. Tlili, I.; Khan, W.A.; Khan, I. Multiple Slips Effects on MHD SA-Al2O3 and SA-Cu Non-Newtonian Nanofluids Flow over a Stretching Cylinder in Porous Medium with Radiation and Chemical Reaction. Results Phys. 2018, 8, 213–222. [Google Scholar] [CrossRef]
  41. Aly, E.H.; Pop, I. MHD Flow and Heat Transfer near Stagnation Point over a Stretching/Shrinking Surface with Partial Slip and Viscous Dissipation: Hybrid Nanofluid versus Nanofluid. Powder Technol. 2020, 367, 192–205. [Google Scholar] [CrossRef]
  42. Eid, M.R.; Nafe, M.A. Thermal Conductivity Variation and Heat Generation Effects on Magneto-Hybrid Nanofluid Flow in a Porous Medium with Slip Condition. Waves Random Complex Media 2020, 1–25. [Google Scholar] [CrossRef]
  43. Abbas, Z.; Naveed, M.; Sajid, M. Hydromagnetic Slip Flow of Nanofluid over a Curved Stretching Surface with Heat Generation and Thermal Radiation. J. Mol. Liquids 2016, 215, 756–762. [Google Scholar] [CrossRef]
  44. Upreti, H.; Pandey, A.K.; Kumar, M. MHD Flow of Ag-Water Nanofluid over a Flat Porous Plate with Viscous-Ohmic Dissipation, Suction/Injection and Heat Generation/Absorption. Alex. Eng. J. 2018, 57, 1839–1847. [Google Scholar] [CrossRef]
  45. Jusoh, R.; Nazar, R.; Pop, I. Impact of Heat Generation/Absorption on the Unsteady Magnetohydrodynamic Stagnation Point Flow and Heat Transfer of Nanofluids. HFF 2019, 30, 557–574. [Google Scholar] [CrossRef]
  46. Murugesan, T.; Kumar, M.D. Effects of Thermal Radiation and Heat Generation on Hydromagnetic Flow of Nanofluid over an Exponentially Stretching Sheet in a Porous Medium with Viscous Dissipation. World Sci. News 2019, 128, 130–147. [Google Scholar]
  47. Hayat, T.; Nadeem, S. Heat Transfer Enhancement with Ag–CuO/Water Hybrid Nanofluid. Results Phys. 2017, 7, 2317–2324. [Google Scholar] [CrossRef]
  48. Hayat, T.; Nadeem, S.; Khan, A.U. Numerical Analysis of Ag–CuO/Water Rotating Hybrid Nanofluid with Heat Generation and Absorption. Can. J. Phys. 2019, 97, 644–650. [Google Scholar] [CrossRef]
  49. Ghosh, S.; Mukhopadhyay, S. Stability Analysis for Model-Based Study of Nanofluid Flow over an Exponentially Shrinking Permeable Sheet in Presence of Slip. Neural Comput. Appl. 2020, 32, 7201–7211. [Google Scholar] [CrossRef]
  50. Mukhopadhyay, S.; Andersson, H.I. Effects of Slip and Heat Transfer Analysis of Flow over an Unsteady Stretching Surface. Heat Mass Transf. 2009, 45, 1447–1452. [Google Scholar] [CrossRef]
  51. Tayebi, T.; Chamkha, A.J. Entropy Generation Analysis during MHD Natural Convection Flow of Hybrid Nanofluid in a Square Cavity Containing a Corrugated Conducting Block. HFF 2019, 30, 1115–1136. [Google Scholar] [CrossRef]
  52. Takabi, B.; Salehi, S. Augmentation of the Heat Transfer Performance of a Sinusoidal Corrugated Enclosure by Employing Hybrid Nanofluid. Adv. Mech. Eng. 2015, 6, 147059. [Google Scholar] [CrossRef]
  53. Ghalambaz, M.; Roşca, N.C.; Roşca, A.V.; Pop, I. Mixed Convection and Stability Analysis of Stagnation-Point Boundary Layer Flow and Heat Transfer of Hybrid Nanofluids over a Vertical Plate. Int. J. Numer. Methods Heat Fluid Flow 2019, 30, 3737–3754. [Google Scholar] [CrossRef]
  54. Oztop, H.F.; Abu-Nada, E. Numerical Study of Natural Convection in Partially Heated Rectangular Enclosures Filled with Nanofluids. Int. J. Heat Fluid Flow 2008, 29, 1326–1336. [Google Scholar] [CrossRef]
  55. Merkin, J.H. On Dual Solutions Occurring in Mixed Convection in a Porous Medium. J. Eng. Math. 1986, 20, 171–179. [Google Scholar] [CrossRef]
  56. Weidman, P.D.; Kubitschek, D.G.; Davis, A.M.J. The Effect of Transpiration on Self-Similar Boundary Layer Flow over Moving Surfaces. Int. J. Eng. Sci. 2006, 44, 730–737. [Google Scholar] [CrossRef]
  57. Hafidzuddin, M.E.H.; Nazar, R.; Arifin, N.M.; Pop, I. Boundary Layer Flow and Heat Transfer over a Permeable Exponentially Stretching/Shrinking Sheet with Generalized Slip Velocity. JAFM 2016, 9, 2025–2036. [Google Scholar] [CrossRef]
Figure 1. Physical model.
Figure 1. Physical model.
Mathematics 09 00030 g001
Figure 2. Plot of Re x 1 / 2 C f with λ for varied A .
Figure 2. Plot of Re x 1 / 2 C f with λ for varied A .
Mathematics 09 00030 g002
Figure 3. Plot of Re x 1 / 2 N u x with λ for varied A .
Figure 3. Plot of Re x 1 / 2 N u x with λ for varied A .
Mathematics 09 00030 g003
Figure 4. Plot of Re x 1 / 2 N u x with λ for varied B .
Figure 4. Plot of Re x 1 / 2 N u x with λ for varied B .
Mathematics 09 00030 g004
Figure 5. Plot of Re x 1 / 2 N u x with λ for varied β .
Figure 5. Plot of Re x 1 / 2 N u x with λ for varied β .
Mathematics 09 00030 g005
Figure 6. Plot of Re x 1 / 2 C f with λ for varied ϕ s 2 .
Figure 6. Plot of Re x 1 / 2 C f with λ for varied ϕ s 2 .
Mathematics 09 00030 g006
Figure 7. Plot of Re x 1 / 2 N u x with λ for varied ϕ s 2 .
Figure 7. Plot of Re x 1 / 2 N u x with λ for varied ϕ s 2 .
Mathematics 09 00030 g007
Figure 8. Plot of Re x 1 / 2 C f with S for varied ϕ s 2 .
Figure 8. Plot of Re x 1 / 2 C f with S for varied ϕ s 2 .
Mathematics 09 00030 g008
Figure 9. Plot of Re x 1 / 2 N u x with S for varied ϕ s 2 .
Figure 9. Plot of Re x 1 / 2 N u x with S for varied ϕ s 2 .
Mathematics 09 00030 g009
Figure 10. Velocity profile f ( η ) for varied A .
Figure 10. Velocity profile f ( η ) for varied A .
Mathematics 09 00030 g010
Figure 11. Temperature profile θ ( η ) for varied A .
Figure 11. Temperature profile θ ( η ) for varied A .
Mathematics 09 00030 g011
Figure 12. Velocity profile f ( η ) for varied S .
Figure 12. Velocity profile f ( η ) for varied S .
Mathematics 09 00030 g012
Figure 13. Temperature profile θ ( η ) for varied S .
Figure 13. Temperature profile θ ( η ) for varied S .
Mathematics 09 00030 g013
Figure 14. Velocity profile f ( η ) for varied ϕ s 2 .
Figure 14. Velocity profile f ( η ) for varied ϕ s 2 .
Mathematics 09 00030 g014
Figure 15. Temperature profile θ ( η ) for varied ϕ s 2 .
Figure 15. Temperature profile θ ( η ) for varied ϕ s 2 .
Mathematics 09 00030 g015
Figure 16. Temperature profile θ ( η ) for varied B .
Figure 16. Temperature profile θ ( η ) for varied B .
Mathematics 09 00030 g016
Figure 17. Temperature profile θ ( η ) for varied β .
Figure 17. Temperature profile θ ( η ) for varied β .
Mathematics 09 00030 g017
Table 1. Hybrid nanofluid physical properties.
Table 1. Hybrid nanofluid physical properties.
PropertiesHybrid Nanofluid Correlations
Density ρ h n f = ρ s 1 ϕ s 1 + ρ s 2 ϕ s 2 + ρ f ( 1 ϕ h n f )
where ϕ h n f = ϕ s 1 + ϕ s 2
Heat Capacity ( ρ C p ) h n f = ( ρ C p ) s 1 ϕ s 1 + ( ρ C p ) s 2 ϕ s 2 + ( ρ C p ) f ( 1 ϕ h n f )
Dynamic Viscosity μ h n f μ f = 1 ( 1 ϕ h n f ) 2.5
Thermal Conductivity k h n f k f = [ 2 k f + ( ϕ s 1 k s 1 + ϕ s 2 k s 2 ϕ h n f ) + 2 ( ϕ s 1 k s 1 + ϕ s 2 k s 2 ) 2 ϕ h n f k f 2 k f ( ϕ s 1 k s 1 + ϕ s 2 k s 2 ) + ( ϕ s 1 k s 1 + ϕ s 2 k s 2 ϕ h n f ) + ϕ h n f k f ]
Table 2. Thermo-physical properties.
Table 2. Thermo-physical properties.
Physical PropertiesWaterAl2O3Cu
ρ ( kg / m 3 ) 997.139708933
C p ( J / kgK ) 4179765385
k ( W / mK ) 0.61340400
Table 3. Tabulation of f ( 0 ) and θ ( 0 ) for viscous fluid ( ϕ s 1 = ϕ s 2 = 0 ) when Pr = 0.7 ,   S = 3 and A = B = β = 0 for the shrinking surface ( λ = 1 ) case.
Table 3. Tabulation of f ( 0 ) and θ ( 0 ) for viscous fluid ( ϕ s 1 = ϕ s 2 = 0 ) when Pr = 0.7 ,   S = 3 and A = B = β = 0 for the shrinking surface ( λ = 1 ) case.
f ( 0 ) θ ( 0 )
First SolutionSecond SolutionFirst SolutionSecond Solution
Ghosh and Mukhopadhyay [49]2.39082−0.972231.77120.84832
Hafidzuddin et al. [57]2.3908−0.97221.77120.8483
Waini et al. [16]2.390814−0.9722471.7712370.848316
Present Results2.390813634−0.9721288681.7712373070.847748272
Table 4. The least eigenvalues for preferred values of A and λ as S = 2.4 ,   ϕ s 2 = 0.01 ,   B = 0.2 and β = 0.02 .
Table 4. The least eigenvalues for preferred values of A and λ as S = 2.4 ,   ϕ s 2 = 0.01 ,   B = 0.2 and β = 0.02 .
A λ Least Eigenvalues
First   Solution   ( γ 1 ) Sec ond   Solution   ( γ 2 )
0.2−1.440.08790−0.08778
−1.4410.06126−0.06120
−1.44190.01326−0.01326
0.6−2.030.13447−0.13450
−2.0360.02518−0.02518
−2.03620.00711−0.00713
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wahid, N.S.; Arifin, N.M.; Khashi’ie, N.S.; Pop, I. Hybrid Nanofluid Slip Flow over an Exponentially Stretching/Shrinking Permeable Sheet with Heat Generation. Mathematics 2021, 9, 30. https://0-doi-org.brum.beds.ac.uk/10.3390/math9010030

AMA Style

Wahid NS, Arifin NM, Khashi’ie NS, Pop I. Hybrid Nanofluid Slip Flow over an Exponentially Stretching/Shrinking Permeable Sheet with Heat Generation. Mathematics. 2021; 9(1):30. https://0-doi-org.brum.beds.ac.uk/10.3390/math9010030

Chicago/Turabian Style

Wahid, Nur Syahirah, Norihan Md Arifin, Najiyah Safwa Khashi’ie, and Ioan Pop. 2021. "Hybrid Nanofluid Slip Flow over an Exponentially Stretching/Shrinking Permeable Sheet with Heat Generation" Mathematics 9, no. 1: 30. https://0-doi-org.brum.beds.ac.uk/10.3390/math9010030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop