Next Article in Journal
Unsupervised Clustering for Hyperspectral Images
Next Article in Special Issue
Effects of Stefan Blowing and Slip Conditions on Unsteady MHD Casson Nanofluid Flow Over an Unsteady Shrinking Sheet: Dual Solutions
Previous Article in Journal
A Hybrid Approach of Interval Neutrosophic Vague Sets and DEMATEL with New Linguistic Variable
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dual Solutions and Stability Analysis of a Hybrid Nanofluid over a Stretching/Shrinking Sheet Executing MHD Flow

1
School of Quantitative Sciences, Universiti Utara Malaysia, Sintok 06010, Kedah, Malaysia
2
KCAET Khairpur Mir’s Sindh Agriculture University, Tandojam Sindh 70060, Pakistan
3
Faculty of Mathematics and Statistics, Ton Duc Thang University, Ho Chi Minh City 72915, Vietnam
4
Center of Excellence for Research in Engineering Materials (CEREM), King Saud University, P.O. Box 800, Al-Riyadh 11421, Saudi Arabia
5
Electrochemistry and Corrosion Laboratory, Department of Physical Chemistry, National Research Centre, El-Behoth St. 33, Dokki, Cairo 12622, Egypt
*
Author to whom correspondence should be addressed.
Submission received: 16 January 2020 / Revised: 2 February 2020 / Accepted: 5 February 2020 / Published: 12 February 2020
(This article belongs to the Special Issue Turbulence and Multiphase Flows and Symmetry)

Abstract

:
In this paper, the unsteady magnetohydrodynamic (MHD) flow of hybrid nanofluid (HNF) composed of C u A l 2 O 3 /water in the presence of a thermal radiation effect over the stretching/shrinking sheet is investigated. Using similarity transformation, the governing partial differential equations (PDEs) are transformed into a system of ordinary differential equations (ODEs), which are then solved by using a shooting method. In order to validate the obtained numerical results, the comparison of the results with the published literature is made numerically as well as graphically and is found in good agreements. In addition, the effects of many emerging physical governing parameters on the profiles of velocity, temperature, skin friction coefficient, and heat transfer rate are demonstrated graphically and are elucidated theoretically. Based on the numerical results, dual solutions exist in a specific range of magnetic, suction, and unsteadiness parameters. It was also found that the values of f ( 0 ) rise in the first solution and reduce in the second solution when the solid volume fraction ϕ C u is increased. Finally, the temporal stability analysis of the solutions is conducted, and it is concluded that only the first solution is stable.

1. Introduction

Researchers are interested in studying the enhancement of heat transfer due to its significant applications in engineering and industries. The heat transfer of convectional fluids such as ethylene glycol, water, and oil can be utilized in various apparatus of engineering, for instance, devices of electrons and heat exchangers. However, these base liquids keep limited thermal conductivity or, in other words, they have low thermal conductivity. To overcome this shortcoming, engineers, mathematicians, and researchers of various fields attempt to improve thermal conductivity of previously mentioned liquids by including a solitary kind of nanosized particles to form a mixture called ‘nanofluid’, which was initially presented by Choi and Eastman [1]. It is easily shown from the previous studies that solid nanoparticles possess the ability increase thermal conductivity and the rate of heat transfer of convectional base fluids. Consequently, numerous analysts and thermal experts have carried out investigations in order to enhance the rate of heat transfer of nanofluid for various perspectives, numerically as well as experimentally. For instance, the examination of the single-phase models of nanofluids was established by Tiwari and Das [2]. This model has been constructed by considering a solid volume fraction of nanoparticles in the base fluid, and, later, the governing equations have been solved numerically by employing the finite volume method. As a result, this model has been widely considered by many scientists, engineers, and mathematicians, such as Benzema et al. [3], Dero et al. [4], Lund et al. [5,6], Dogonchi, et al. [7,8], Amini, et al. [9], Zaib et al. [10], Raza et al. [11], Rasool et al. [12,13], Dinarvand et al. [14], and Roşca, et al. [15] to investigate different types of flow. Furthermore, to keep the demand of the high heat transfer rate from industries and other sectors, researchers have introduced a new kind of the nanofluid by considering the two different types of the solid particles in the single convectional base fluid. This kind of nanofluid is known as a hybrid nanofluid [16]. It is worth mentioning that thermal conductivity of regular base fluid in hybrid nanofluid is higher than in simple nanofluid.
Due to this nice property, studies of different types of hybrid nanofluids have gained attention. Devi and Devi [17] examined the C u A l 2 O 3 / water hybrid nanofluid. The obtained numerical results were compared with the experimental results of Suresh et al. [18] and were found in excellent agreement. Toghraie et al. [19] further adopted this model for the Z n O T i O 2 / E G hybrid nanofluid. Meanwhile, Moghadassi et al. [20] studied the A l 2 O 3 Cu / water base hybrid nanofluid numerically and found that “for the hybrid nanofluids, the average Nusselt number increase was 4.73% and 13.46% in compared to A l 2 O 3 water and pure water, respectively”. Moreover, this model has been widely employed by many researchers, such as Hayat et al. [21], Saba et al. [22], Acharya et al. [23], Afridi et al. [24], Shafiq et al. [25,26], and Manh et al. [27]. Furthermore, Khashi’ie et al. [28] found dual solutions of magnetohydrodynamic (MHD) flow of a hybrid nanofluid in the presence of Joule heating and noticed that higher values of Eckert number do not affect boundary layer separation. Lund et al. [29] studied hybrid nanofluids by considering copper and alumina as solid particles with water as a base fluid. Dual solutions were noticed in the presence of high suction. Dual solutions of a hybrid nanofluid over a vertical thin needle were investigated by Waini et al. [30]. They performed a stability analysis on dual solutions and discovered that only the first solution is stable. Later, Waini et al. [31] considered an unsteady flow of hybrid nanofluid over a stretching/shrinking surface and found that dual solutions appeared in certain ranges of an unsteadiness parameter.
This paper is an extension of the study carried out by Waini et al. [31] by considering the unsteady magnetohydrodynamic (MHD) flow of the C u A l 2 O 3 / water hybrid nanofluid in the presence of thermal radiation over the stretching/shrinking parameter. A model of Tiwari and Das [2] is used to deal with governing equations by including two solid nanoparticles, namely copper ( C u ) and alumina ( A l 2 O 3 ), with water as a base fluid. To the best of authors’ knowledge and based on a survey of previous literature, this kind of study has not been investigated.

2. Problem Formulation

Let us consider an unsteady MHD incompressible flow of the C u A l 2 O 3 / water nanofluid flow on a stretching/shrinking sheet in the presence of thermal radiation effect. Figure 1 shows the physical model and coordinate system of the concerned problem. The velocity of wall mass transfer is v w ( x ) = ( ϑ f c ( 1 ε t ) ) 1 2 f ( η ) where c is a positive constant. It is also assumed that the flow is subjected to a transverse magnetic field of strength B = B 0 ( 1 ε t ) 1 / 2 where B 0 is the constant applied magnetic field. By including the considered assumptions, the equation of mass, momentum, and energy in the form of a boundary layer are expressed in Tiwari and Das’ model [2], as below:
u x   +   v y = 0
u t   +   u u x   +   v u y = μ h n f ρ h n f 2 u y 2 σ * B 2 u ρ h n f
T t   +   u T x   +   v T y = [ k h n f ( ρ c p ) h n f   +   16 σ 1 T 3 3 k * ( ρ c p ) h n f ] 2 T y 2
subject to the following boundary conditions:
{ t < 0 ,   u = v = 0 , T = T t 0 , v = v w ,   u = λ u w , T = T w   a t   y = 0   u 0 ,   T T   a s   y
where u w ( x , t ) =   = c x ( 1 ε t ) is the velocity of the surface. In the current study, thermophoresis properties of Waini et al. [31] are adopted. These properties of hybrid nanofluid are presented in Table 1 and Table 2.
In order to reduce the governing equations into a system of ordinary differential equations (ODEs), the following similarity transformation variables are employed:
η = ( c ϑ f ( 1 ε t ) ) 1 2 y ,   u = c x ( 1 ε t ) f ( η ) ,   v = ( ϑ f c ( 1 ε t ) ) 1 2 f ( η ) , θ ( η )   =   T T T w T
Substituting Equation (5) in Equations (1)–(3), Equation (1) is definitely fulfilled, and Equations (2)–(3) take the accompanying dimensionless form of ODEs:
f   +   ξ 1 { f f f 2 A ( 0.5 η f     +   f ) } ξ 2 M f = 0
1 ξ 3 P r   [ ( k h n f / k f )   +   4 R d 3 ] θ   +   f θ 0.5 A η θ = 0
subject to the following boundary conditions:
{ f ( 0 ) = S , f ( 0 ) = λ ,   θ ( 0 ) = 1 f ( η ) 0 , θ ( η ) 0   as   η
The reduced quantities are expressed as A = ε c , M = σ * B 0 2 c   ρ f ,   P r = ϑ f α f ,   R d = 4 σ 1 T 3 k * k f .
{ ξ 1 = ( 1 ϕ A l 2 O 3 ) 2.5 ( 1 ϕ C u ) 2.5 { ( 1 ϕ C u ) [ 1 ϕ A l 2 O 3   +   ϕ A l 2 O 3 ( ρ A l 2 O 3 ρ f ) ]   +   ϕ C u ( ρ C u ρ f ) } ξ 2 = ( 1 ϕ C u ) 2.5 ( 1 ϕ A l 2 O 3 ) 2.5 ξ 3 = { ( 1 ϕ C u ) [ 1 ϕ A l 2 O 3   +   ϕ A l 2 O 3 ( ρ c p ) A l 2 O 3 ( ρ c p ) f ]   +   ϕ C u ( ρ c p ) C u ( ρ c p ) f }
The physical quantities of attention are the skin friction coefficient C f and the local Nusselt number N u x , which are expressed as follows:
C f = μ h n f ρ f u w 2 ( u y ) | y = 0 ,     N u x = x k h n f k f ( T w T ) ( T y ) | y = 0
Applying Equation (5) in Equation (10), leads to:
R e C f = 1 ( 1 ϕ A l 2 O 3 ) 2.5 ( 1 ϕ C u ) 2.5 f ( 0 ) ;   1 R e N u x = k h n f k f [ 1   +   4 R d 3 ] θ ( 0 )
where R e = x u w ϑ f is the local Reynold number.

3. Stability Analysis

In order to perform a temporal analysis of the solutions’ stability, introducing a new dimensionless time-dependent similarity transformation variable is required, as recommended by Merkin [32], Dero et al. [33,34], and Lund et al. [35,36,37]. Letting τ = c t ( 1 ε t ) yields the following new similarity transformation variables:
{ u = c x ( 1 ε t ) f ( η ,   τ ) η ,   v = ( ϑ f c ( 1 ε t ) ) 1 2 f ( η , τ ) , θ ( η , τ )   =   T T T w T η = ( c ϑ f ( 1 ε t ) ) 1 2 y ;   τ = c t ( 1 γ t )
Substituting Equation (12) in Equations (2) and (3) leads to:
3 f η 3   +   ξ 1   { f 2 f η 2 ( f η ) 2 A ( 0.5 η 2 f η 2   +   f η ) ( 1   +   A τ ) 2 f τ η } ξ 2 M f η = 0
1 ξ 3 P r   [ ( k h n f / k f )   +   4 R d 3 ] 2 θ η 2   +   f θ η 0.5 A η θ η ( 1   +   A τ ) θ τ = 0
subject to the following boundary conditions:
{ f ( 0 , τ ) = S , f ( 0 ,   τ ) η = λ , θ ( 0 , τ ) = 1 f η ( ,   τ ) η =   θ ( η , τ ) = 0   a s   η
To check the stability of the steady flow solutions, where f ( η ) = f 0 ( η ) , θ ( η ) = θ 0 ( η ) , and ( η ) = 0 ( η ) , the following boundary value problem (8)-(11) must be satisfied:
{ f ( η , τ ) = f 0 ( η )   +   e γ τ F ( η , τ ) θ ( η , τ ) = θ 0 ( η )   +   e γ τ G ( η , τ )
where F ( η ) and G ( η ) are small relatives of f 0 ( η ) and θ 0 ( η ) , respectively, and γ is an unknown eigenvalue parameter, which need to be determined. Substituting Equation (16) in Equations (13)–(14) by considering τ = 0 , gives the following linearized eigenvalue problem:
F 0   +   ξ 1 { f 0 F 0 2 f 0 F 0   +   F 0 f 0 A ( 0.5 η F 0   +   F 0 )   +   F 0 } ξ 2 M F 0 = 0
  1 ξ 3 P r   [ ( k h n f / k f )   +   4 R d 3 ] G 0   +   f 0 G 0   +   F 0 θ 0 0.5 η A G 0   +   G 0 = 0
subject to the following boundary conditions:
{ F 0 ( 0 ) = 0 ,   F 0 ( 0 ) = 0 ,   G 0 ( 0 ) = 0 F 0 ( η ) 0 ,   G 0 ( η ) 0   a s   η
According to Haris et al. [38] and Khashi’ie et al. [39], one boundary condition has to be relaxed in order to obtain the value(s) of the smallest eigenvalue. It should be noted that this relaxation in boundary condition does not affect the results [40]. In this particular problem, F 0 ( η ) 0   as   η is relaxed into F 0 ( 0 ) = 1 .

4. Result and Discussion

In this section, Equations (6) and (7) with boundary conditions (8) have been solved numerically by employing the shooting method in Maple software. The shooting method has been employed by many researchers [41,42,43,44]. As suggested by Devi and Deci [17], the solid volume fraction ϕ A l 2 O 3 = 0.1 is kept as a constant in many graphs. Furthermore, the solid volume fraction ϕ C u is added in order to make the C u A l 2 O 3 / water hybrid nanofluid. As recommended by Dero et al. [4] and Iqbal et al. [45], the Prandtl number P r = 6.2 is kept as constant and employed to produce the results of the current study. Table 3 shows the comparison of the numerical values of f ( 0 ) and θ ( 0 ) with the published results of waini et al. [31] for Cu–water nanofluid. It is worth mentioning that Waini et al. [31] used the 3-stage Labatto IIIa formula in a BVP4C solver to solve the resultant similarity equations. It is found that the current results show outstanding agreement with the results of Waini et al. [31]. Therefore, it can be concluded that the accuracy of the present method is fine, and it can be used in this study. Figure 2 is demonstrated for the comparison purpose in order to check the accuracy of the used method. It is observed that the critical values of the current study have a good agreement with the critical values of published paper (allude Figure 6 in Waini et al. [31]). After comparison, this method can be used confidently in this problem.
Figure 3 and Figure 4 show the variation of f ( 0 ) and θ ( 0 ) against A for ϕ C u = 0.001 ,   0.01 , and 0.1 . By keeping A fixed, an expansion of ϕ C u creates a reduction in the heat transfer rate, which implies that the heat transfer coefficient is decreased on the sheet, while an opposite nature is observed for the f ( 0 ) in the second solution. On the other hand, the coefficient of skin friction decreases in the first solution. Moreover, it is noticed that multiple solutions exist for A   > A c , whereas no solution exists for A < A c . It is worth mentioning here that A c is the critical value where multiple solutions exist. Variations of f ( 0 ) and θ ( 0 ) along M for numerous values of ϕ C u were drawn in Figure 5 and Figure 6, respectively. It is observed from them that when the solid volume fraction ϕ C u increases the values of the critical point of M are also increased. Furthermore, M c 1 = 0.4492 ,   M c 2 = 0.3912 are the respective critical values of M for ϕ C u = 0.001 and 0.01 . It is examined that the critical values of M get smaller for higher values of the solid volume fraction ϕ C u , which is the physical reason of the extension of the thickness of boundary layer separation. Moreover, it can be seen that the values of f ( 0 ) rises in the first solution and reduces in the second solution when the solid volume fraction ϕ C u is increased. In the meantime, the rate of heat transfer decreases for the unstable solution and increases for the stable solution with the increasing of ϕ C u .
In detail, Figure 7 and Figure 8 show the variation of f ( 0 ) and θ ( 0 ) with A for the numerous values of the mass suction parameter S . It is observed that the domain of the solutions increases along the critical value of A , moving to one side when the mass suction parameter S is decreased. In light of our calculations, A c 1 = 10.1708 ,   A c 2 = 7.2207 , and A c 3 = 5.2042 are the respective critical values of A for S = 1.8 ,   1.75 and 1.7 . At the same time, it is examined that the smaller critical value is for S = 1.8 as compared to S = 1.75 and S = 1.7 . Henceforth, boundary layer separation is delayed for the higher values of the mass suction S . In addition, the coefficient of skin friction increases in the first solution with the developing of S , whereas it declines in the second solution. The expanding of S prompts to rise the effect of the heat transfer rate θ ( 0 ) in both solutions. Figure 9 and Figure 10 demonstrate the variation of f ( 0 ) and θ ( 0 ) with S for the numerous values of ϕ C u . It can be found that the increments in ϕ C u lead to smaller critical values of S c . According to our computation, S c 1 = 1.7886 ,   S c 2 = 1.78234 , and S c 3 = 1.7496 are the respective critical values of S for ϕ C u = 0.001 ,   0.01 , and 0.1 . Furthermore, the heat transfer rate rises in both solutions for the higher values of the mass suction S . This is due to fact that mass suction S helps the molecules of the hybrid nanofluid to transfer the heat effectively. On the other hand, f ( 0 ) increases (decreases) in the first (second) solution when ϕ C u is enhanced.
Figure 11 and Figure 12 illustrate the profiles of velocity f ( η ) and the temperature θ ( η ) for S. It was demonstrated that the rise in S decreases the velocity of the fluid in the first solution but it increases the velocity of the fluid in the second solution. On the other hand, it is perceived that the temperature of fluid decreases with the expansion of the mass suction S . The effect of the thermal radiation parameter on θ ( η ) is shown in Figure 13. It can be noticed that the thickness of the thermal boundary layer, as well as temperature of the fluid, increases with the increasing values of R d .
Finally, the linearized eigenvalue Equations (17) and (18) have been solved using a 3-stage Labatto IIIa formula in a BVP4C solver in MATLAB 2017b. The results of the stability analysis of the solutions are demonstrated graphically in Figure 14. It can be easily observed from the figure that the initial decay of disturbance is denoted by the positive value of γ and a stable flow, while the initial growth of disturbance is indicated by the negative values of γ and an unstable flow. Additionally, it is noticed that γ tends to zero at the critical values of A from both solutions.
Table 1. Thermophysical properties of the hybrid nanofluid [46].
Table 1. Thermophysical properties of the hybrid nanofluid [46].
PropertiesHybrid Nanofluid
Dynamic viscosity μ h n f = μ f ( 1 ϕ A l 2 O 3 ) 2.5 ( 1 ϕ C u ) 2.5
Density ρ h n f = ( 1 ϕ C u ) [ ( 1 ϕ A l 2 O 3 ) ρ f   +   ϕ A l 2 O 3 ρ A l 2 O 3 ]   +   ϕ C u ρ C u
Thermal conductivity k h n f = k C u   +   2 k n f 2 ϕ C u ( k n f k C u ) k C u   +   2 k n f   +   ϕ C u ( k n f k C u ) × ( k n f )
where k n f = k A l 2 O 3   +   2 k f 2 ϕ A l 2 O 3 ( k f k A l 2 O 3 ) k A l 2 O 3   +   2 k f   +   ϕ A l 2 O 3 ( k f k A l 2 O 3 ) × ( k f )
Heat capacity ( ρ c p ) h n f = ( 1 ϕ C u ) [ ( 1 ϕ A l 2 O 3 ) ( ρ c p ) f   +   ϕ A l 2 O 3 ( ρ c p ) A l 2 O 3 ]   +   ϕ C u ( ρ c p ) C u
Table 2. The thermo physical properties of the base fluid (water) and the nanoparticles [47,48].
Table 2. The thermo physical properties of the base fluid (water) and the nanoparticles [47,48].
Fluids ρ c p ( J / kg   K ) k (W/m K)
Alumina ( A l 2 O 3 )397076540
Copper (Cu)8933385400
Water ( H 2 O )997.141790.613
Table 3. Values of f ( 0 ) and θ ( 0 )   for the Cu–water nanofluid ( ϕ C u = 0.2 ) with various values of A   when M = R d = 0 , S = 2.1 , P r = 6.2 , and λ = 1 .
Table 3. Values of f ( 0 ) and θ ( 0 )   for the Cu–water nanofluid ( ϕ C u = 0.2 ) with various values of A   when M = R d = 0 , S = 2.1 , P r = 6.2 , and λ = 1 .
A ϕ A l 2 O 3 f ( 0 ) θ ( 0 )
Waini[31]PresentResultsWaini[31]PresentResults
1st Soln2nd Soln1st Soln2nd Soln1st Soln2nd Soln1st Soln2nd Soln
−102.194247−1.4912812.19424658−1.491281197.0736806.8845487.07367986.88454793
−10.1----1.60888878−0.69818494----5.22381505.0737595
−301.521197−4.1447461.52119229−4.144745727.4971517.2961767.49715087.2961762
−30.1----0.83697819−2.53135605----5.63167945.4925121
−500.844435−6.4315070.84443506−6.431507387.8584467.6578017.85844577.6578012
−50.01----0.77597819−6.20435575----7.63756957.4426329
−90−0.517287−10.58983−0.51728551−10.58983038.4733168.2776768.47331638.27767616
−90.01----−0.60021181−10.2417356----8.247096398.05819003

5. Conclusions

The problem of the unsteady MHD flow of the C u A l 2 O 3 / water hybrid nanofluid in the presence of the thermal radiation effect over the stretching/shrinking sheet is examined. Two methods are adopted to carry out the numerical computations of the current problem, a shooting method is used to get the dual solutions of ODEs (6–8) in MAPLE software 2018, and the computation of the stability analysis is carried out by employing a 3-stage Labatto IIIa formula in a BVP4C solver in MATLAB 2017b. The impacts of ϕ C u unsteadiness parameter S and suction parameter S on the f ( 0 ) ,   θ ( 0 ) ,   f ( η ) , and θ ( η ) were exhibited graphically and examined. The results show that the development (reduction) of f ( 0 ) and the decline (reduction) of θ ( 0 ) past a shrinking sheet is found with the expansion of ϕ C u in the first (second) solution. The temperature of fluid in both solutions increase (decrease) for the higher values of R d and S . The fluid velocity decreases when the effect of suction parameter S enhanced in the first solution, whereas it increases in the second solution. Furthermore, it is found that the skin friction coefficient increases with an increase in solid volume friction of ϕ C u and suction parameter S in the first solution, while it decreases in the second solution. On the other hand, an improvement in the heat transfer rate is noticed when suction parameter S is increased for both solutions, whereas the totally opposite trend is noticed for the case of solid volume friction of ϕ C u . For the certain ranges of the suction parameter S , unsteadiness parameter A , and magnetic parameter M , dual solutions are noticed. Furthermore, the results of the temporal stability analysis confirm that the first solution is stable, whereas the second solution is unstable. In future, this work can be continued to assess the impact of viscous dissipation, joule heating, boundary slip conditions, and convective boundary conditions. Other than these impacts, it can be extended for three-dimensional boundary layer flows by conserving various possible effects.

Abbreviation

Nomenclature
T   0 a constant ϕ A l 2 O 3 nanoparticle volume fraction of the iron oxide.
T ambient temperatureMmagnetic parameter
differentiation with respect to η K 1 Porous parameter
R d Thermal radiationcconstant
ρ h n f effective density of hybrid nanofluid P r Prandtl number
ρ n f effective density of nanofluid C f skin friction coefficient
μ h n f effective dynamic viscosity of hybrid nanofluid γ 1 smallest eigen value
μ n f effective dynamic viscosity of nanofluid τ Stability transformed variable
σ * electrical conductivity v w suction/injection velocity
f fluid fractionTTemperature
MHartmann/magnetic number k h n f thermal conductivity of the hybrid nanofluid
( ρ c p ) h n f heat capacitance of the hybrid nanofluid k n f thermal conductivity of the nanofluid
( ρ c p ) n f heat capacitance of the nanofluidttime
hnfHybrid nanofluid η transformed variable
N u local Nusselt number A Unsteadiness parameter
R e local Reynolds number T   w variable temperature at the sheet
BMagnetic fieldu, vvelocity components
n f nanofluid fraction λ shrinking/stretching parameter
k * mean absorption coefficient σ 1 Stefan–Boltzmann constant
ϕ C u nanoparticle volume fraction of the copper S S < 0 for suction parameter and S > 0 for blowing parameter

Author Contributions

L.A.L. derived the equations and generated the results and wrote the paper. Z.O. formulated the model and proof read the manuscript. I.K. checked the whole manuscript. I.K. and E.-S.M.S. derived the equations of stability analysis and wrote the introduction section. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by Researchers Supporting Project number (RSP-2019/33), King Saud University, Riyadh, Saudi Arabia. This research is also supported by the Universiti Utara Malaysia (UUM).

Acknowledgments

Researchers Supporting Project number (RSP-2019/33), King Saud University, Riyadh, Saudi Arabia. The first author is thankful to the School of Quantitative Sciences for providing a good environment to conduct this research in the postgraduate lab, special thanks to Z.O. and I.K.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; No. ANL/MSD/CP-84938; CONF-951135-29; Argonne National Lab: Lemont, IL, USA, 1995. [Google Scholar]
  2. Tiwari, R.K.; Das, M.K. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int. J. Heat Mass Transf. 2007, 50, 2002–2018. [Google Scholar] [CrossRef]
  3. Benzema, M.; Benkahla, Y.K.; Labsi, N.; Ouyahia, S.E.; El Ganaoui, M. Second law analysis of MHD mixed convection heat transfer in a vented irregular cavity filled with Ag-MgO/water hybrid nanofluid. J. Therm. Anal. Calorim. 2019, 137, 1113–1132. [Google Scholar] [CrossRef]
  4. Dero, S.; Rohni, A.M.; Saaban, A. The dual solutions and stability analysis of nanofluid flow using tiwari-das modelover a permeable exponentially shrinking surface with partial slip conditions. J. Eng. Appl. Sci. 2019, 14, 4569–4582. [Google Scholar] [CrossRef] [Green Version]
  5. Lund, L.A.; Omar, Z.; Khan, I.; Dero, S. Multiple solutions of Cu-C6H9NaO7 and Ag-C6H9NaO7 nanofluids flow over nonlinear shrinking surface. J. Cent. South Univ. 2019, 26, 1283–1293. [Google Scholar] [CrossRef]
  6. Lund, L.A.; Omar, Z.; Khan, U.; Khan, I.; Baleanu, D.; Nisar, K.S. Stability analysis and dual solutions of micropolar nanofluid over the inclined stretching/shrinking surface with convective boundary condition. Symmetry 2020, 12, 74. [Google Scholar] [CrossRef] [Green Version]
  7. Dogonchi, A.S.; Chamkha, A.J.; Hashemi-Tilehnoee, M.; Seyyedi, S.M.; Ganji, D.D. Effects of homogeneous-heterogeneous reactions and thermal radiation on magneto-hydrodynamic Cu-water nanofluid flow over an expanding flat plate with non-uniform heat source. J. Cent. South Univ. 2019, 26, 1161–1171. [Google Scholar] [CrossRef]
  8. Dogonchi, A.S.; Ismael, M.A.; Chamkha, A.J.; Ganji, D.D. Numerical analysis of natural convection of Cu-water nanofluid filling triangular cavity with semicircular bottom wall. J. Therm. Anal. Calorim. 2019, 135, 3485–3497. [Google Scholar] [CrossRef]
  9. Amini, Y.; Akhavan, S.; Izadpanah, E. Vortex-induced vibration of a cylinder in pulsating nanofluid flow. J. Therm. Anal. Calorim. 2019. [Google Scholar] [CrossRef]
  10. Zaib, A.; Khan, M.; Shafie, S. Boundary-layer flow of a cu-water nanofluid over a permeable shrinking cylinder with homogenous-hetrogenous reactions: Dual solutions. Therm. Sci. 2019, 23, 295–306. [Google Scholar] [CrossRef] [Green Version]
  11. Raza, J.; Rohni, A.M.; Omar, Z.; Awais, M. Heat and mass transfer analysis of MHD nanofluid flow in a rotating channel with slip effects. J. Mol. Liq. 2016, 219, 703–708. [Google Scholar] [CrossRef]
  12. Rasool, G.; Zhang, T.; Shafiq, A. Marangoni effect in second grade forced convective flow of water based nanofluid. J. Adv. Nanotechnol. 2019, 1, 50. [Google Scholar] [CrossRef]
  13. Rasool, G.; Shafiq, A.; Khalique, C.M.; Zhang, T. Magneto-hydrodynamic Darcy-Forchheimer nanofluid flow over nonlinear stretching sheet. Phys. Scr. 2019, 94, 10. [Google Scholar] [CrossRef]
  14. Dinarvand, S. Nodal/saddle stagnation-point boundary layer flow of CuO-Ag/water hybrid nanofluid: A novel hybridity model. Microsyst. Technol. 2019, 25, 2609–2623. [Google Scholar] [CrossRef]
  15. Roşca, N.C.; Roşca, A.V.; Pop, I. Unsteady separated stagnation-point flow and heat transfer past a stretching/shrinking sheet in a copper-water nanofluid. Int. J. Numer. Methods Heat Fluid Flow 2019. [Google Scholar] [CrossRef]
  16. Ahmed, N.; Saba, F.; Khan, U.; Mohyud-Din, S.T.; Sherif, E.S.M.; Khan, I. Nonlinear thermal radiation and chemical reaction effects on a Cu-CuO/NaAlg hybrid nanofluid flow past a stretching curved surface. Processes 2019, 7, 962. [Google Scholar] [CrossRef] [Green Version]
  17. Devi, S.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]
  18. Suresh, S.; Venkitaraj, K.P.; Selvakumar, P.; Chandrasekar, M. Synthesis of Al2O3-Cu/water hybrid nanofluids using two step method and its thermo physical properties. Colloids Surf. A Physicochem. Eng. Asp. 2011, 388, 41–48. [Google Scholar] [CrossRef]
  19. Toghraie, D.; Chaharsoghi, V.A.; Afrand, M. Measurement of thermal conductivity of ZnO-TiO2/EG hybrid nanofluid. J. Therm. Anal. Calorim. 2016, 125, 527–535. [Google Scholar] [CrossRef]
  20. Moghadassi, A.; Ghomi, E.; Parvizian, F. A numerical study of water based Al2O3 and Al2O3-Cu hybrid nanofluid effect on forced convective heat transfer. Int. J. Therm. Sci. 2015, 92, 50–57. [Google Scholar] [CrossRef]
  21. 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]
  22. Saba, F.; Ahmed, N.; Khan, U.; Waheed, A.; Rafiq, M.; Mohyud-Din, S. Thermophysical analysis of water based (Cu-Al2O3) hybrid nanofluid in an asymmetric channel with dilating/squeezing walls considering different shapes of nanoparticles. Appl. Sci. 2018, 8, 1549. [Google Scholar] [CrossRef] [Green Version]
  23. Acharya, N.; Bag, R.; Kundu, P.K. Influence of hall current on radiative nanofluid flow over a spinning disk: A hybrid approach. Phys. E Low-Dimens. Syst. Nanostruct. 2019, 111, 103–112. [Google Scholar] [CrossRef]
  24. Afridi, M.I.; Alkanhal, T.A.; Qasim, M.; Tlili, I. Entropy Generation in Cu-Al2O3-H2O Hybrid Nanofluid Flow over a Curved Surface with Thermal Dissipation. Entropy 2019, 21, 941. [Google Scholar] [CrossRef] [Green Version]
  25. Shafiq, A.; Khan, I.; Rasool, G.; Sherif, E.S.M.; Sheikh, A.H. Influence of single-and multi-wall carbon nanotubes on magnetohydrodynamic stagnation point nanofluid flow over variable thicker surface with concave and convex effects. Mathematics 2020, 8, 104. [Google Scholar] [CrossRef] [Green Version]
  26. Shafiq, A.; Zari, I.; Rasool, G.; Tlili, I.; Khan, T.S. On the MHD casson axisymmetric marangoni forced convective flow of nanofluids. Mathematics 2019, 7, 1087. [Google Scholar] [CrossRef] [Green Version]
  27. Manh, T.D.; Nam, N.D.; Abdulrahman, G.K.; Moradi, R.; Babazadeh, H. Impact of MHD on hybrid nanomaterial free convective flow within a permeable region. J. Therm. Anal. Calorim. 2019. [Google Scholar] [CrossRef]
  28. Khashi’ie, N.S.; Arifin, N.M.; Nazar, R.; Hafidzuddin, E.H.; Wahi, N.; Pop, I. Magnetohydrodynamics (MHD) axisymmetric flow and heat transfer of a hybrid nanofluid past a radially permeable stretching/shrinking sheet with joule heating. Chin. J. Phys. 2019. [Google Scholar] [CrossRef]
  29. Lund, L.A.; Omar, Z.; Khan, I.; Seikh, A.H.; Sherif, E.S.M.; Nisar, K.S. Stability analysis and multiple solution of Cu-Al2O3/H2O nanofluid contains hybrid nanomaterials over a shrinking surface in the presence of viscous dissipation. J. Mater. Res. Technol. 2019, 9, 421–432. [Google Scholar] [CrossRef]
  30. Waini, I.; Ishak, A.; Pop, I. Hybrid nanofluid flow and heat transfer past a vertical thin needle with prescribed surface heat flux. Int. J. Numer. Methods Heat Fluid Flow 2019. [Google Scholar] [CrossRef]
  31. Waini, I.; Ishak, A.; Pop, I. Unsteady flow and heat transfer past a stretching/shrinking sheet in a hybrid nanofluid. Int. J. Heat Mass Transf. 2019, 136, 288–297. [Google Scholar] [CrossRef]
  32. Merkin, J.H. On dual solutions occurring in mixed convection in a porous medium. J. Eng. Math. 1986, 20, 171–179. [Google Scholar] [CrossRef]
  33. Dero, S.; Rohni, A.M.; Saaban, A.; Khan, I. Dual solutions and stability analysis of micropolar nanofluid flow with slip effect on stretching/shrinking surfaces. Energies 2019, 12, 4529. [Google Scholar] [CrossRef] [Green Version]
  34. Dero, S.; Uddin, M.J.; Rohni, A.M. Stefan blowing and slip effects on unsteady nanofluid transport past a shrinking sheet: Multiple solutions. Heat Transf. Asian Res. 2019, 6, 2047–2066. [Google Scholar] [CrossRef]
  35. Lund, L.A.; Omar, Z.; Khan, I. Mathematical analysis of magnetohydrodynamic (MHD) flow of micropolar nanofluid under buoyancy effects past a vertical shrinking surface: Dual solutions. Heliyon 2019, 5, e02432. [Google Scholar] [CrossRef] [Green Version]
  36. Lund, L.A.; Omar, Z.; Dero, S.; Khan, I. Linear stability analysis of MHD flow of micropolar fluid with thermal radiation and convective boundary condition: Exact solution. Heat Transf. Asian Res. 2019. [Google Scholar] [CrossRef]
  37. Lund, L.A.; Omar, Z.; Khan, I.; Kadry, S.; Rho, S.; Mari, I.A.; Nisar, K.S. Effect of viscous dissipation in heat transfer of MHD flow of micropolar fluid partial slip conditions: Dual solutions and stability analysis. Energies 2019, 12, 4617. [Google Scholar] [CrossRef] [Green Version]
  38. Lund, L.A.; Omar, Z.; Khan, I.; Raza, J.; Sherif, E.S.M.; Seikh, A.H. Magnetohydrodynamic (MHD) flow of micropolar fluid with effects of viscous dissipation and joule heating over an exponential shrinking sheet: Triple solutions and stability analysis. Symmetry 2020, 12, 142. [Google Scholar] [CrossRef] [Green Version]
  39. Khashi’ie, N.S.; Arifin, N.M.; Rashidi, M.M.; Hafidzuddin, E.H.; Wahi, N. Magnetohydrodynamics (MHD) stagnation point flow past a shrinking/stretching surface with double stratification effect in a porous medium. J. Therm. Anal. Calorim. 2019. [Google Scholar] [CrossRef]
  40. Waini, I.; Ishak, A.; Pop, I. Transpiration effects on hybrid nanofluid flow and heat transfer over a stretching/shrinking sheet with uniform shear flow. Alex. Eng. J. 2019. [Google Scholar] [CrossRef]
  41. Asma, M.; Othman, W.A.M.; Muhammad, T. Numerical study for darcy-forchheimer flow of nanofluid due to a rotating disk with binary chemical reaction and arrhenius activation energy. Mathematics 2019, 7, 921. [Google Scholar] [CrossRef] [Green Version]
  42. Alarifi, I.M.; Abokhalil, A.G.; Osman, M.; Lund, L.A.; Ayed, M.B.; Belmabrouk, H.; Tlili, I. MHD flow and heat transfer over vertical stretching sheet with heat sink or source effect. Symmetry 2019, 11, 297. [Google Scholar] [CrossRef] [Green Version]
  43. Basir, M.F.M.; Kumar, R.; Ismail, A.I.M.; Sarojamma, G.; Narayana, P.S.; Raza, J.; Mahmood, A. Exploration of thermal-diffusion and diffusion-thermal effects on the motion of temperature-dependent viscous fluid conveying microorganism. Arab. J. Sci. Eng. 2019, 44, 8023–8033. [Google Scholar] [CrossRef]
  44. Wang, Q. Optimal strokes of low Reynolds number linked-sphere swimmers. Appl. Sci. 2019, 9, 4023. [Google Scholar] [CrossRef] [Green Version]
  45. Iqbal, Z.; Akbar, N.S.; Azhar, E.; Maraj, E.N. Performance of hybrid nanofluid (Cu-CuO/water) on MHD rotating transport in oscillating vertical channel inspired by Hall current and thermal radiation. Alex. Eng. J. 2018, 57, 1943–1954. [Google Scholar] [CrossRef]
  46. Khashi’ie, N.S.; Arifin, N.M.; Hafidzuddin, E.H.; Wahi, N. Thermally stratified flow of Cu-Al2O3/water hybrid nanofluid past a permeable stretching/shrinking circular cylinder. J. Adv. Res. Fluid Mech. Therm. Sci. 2019, 63, 154–163. [Google Scholar]
  47. Rahman, A.N.H.; Bachok, N.; Rosali, H. Numerical solutions of MHD stagnation-point flow over an exponentially stretching/shrinking sheet in a nanofluid. In Journal of Physics: Conference Series; IOP Publishing: Bristol, UK, 2019; Volume 1366. [Google Scholar]
  48. Zaib, A.; Khan, U.; Shah, Z.; Kumam, P.; Thounthong, P. Optimization of entropy generation in flow of micropolar mixed convective magnetite (Fe3O4) ferroparticle over a vertical plate. Alex. Eng. J. 2019, 58, 1461–1470. [Google Scholar] [CrossRef]
Figure 1. Physical models and coordinate systems.
Figure 1. Physical models and coordinate systems.
Symmetry 12 00276 g001
Figure 2. Comparison with the 6th figure of Waini et al. [31].
Figure 2. Comparison with the 6th figure of Waini et al. [31].
Symmetry 12 00276 g002
Figure 3. Variation of f ( 0 ) with A for various values of ϕ C u .
Figure 3. Variation of f ( 0 ) with A for various values of ϕ C u .
Symmetry 12 00276 g003
Figure 4. Variation of θ ( 0 ) with A for various values of ϕ C u .
Figure 4. Variation of θ ( 0 ) with A for various values of ϕ C u .
Symmetry 12 00276 g004
Figure 5. Variation of f ( 0 ) with M for various values of ϕ C u .
Figure 5. Variation of f ( 0 ) with M for various values of ϕ C u .
Symmetry 12 00276 g005
Figure 6. Variation of θ ( 0 ) with M for various values of ϕ C u .
Figure 6. Variation of θ ( 0 ) with M for various values of ϕ C u .
Symmetry 12 00276 g006
Figure 7. Variation of f ( 0 ) with A for various values of S .
Figure 7. Variation of f ( 0 ) with A for various values of S .
Symmetry 12 00276 g007
Figure 8. Variation of θ ( 0 ) with A for various values of S .
Figure 8. Variation of θ ( 0 ) with A for various values of S .
Symmetry 12 00276 g008
Figure 9. Variation of f ( 0 ) with S for various values of ϕ C u .
Figure 9. Variation of f ( 0 ) with S for various values of ϕ C u .
Symmetry 12 00276 g009
Figure 10. Variation of θ ( 0 ) with S for various values of ϕ C u .
Figure 10. Variation of θ ( 0 ) with S for various values of ϕ C u .
Symmetry 12 00276 g010
Figure 11. Variation of f ( η ) with η for various values of S .
Figure 11. Variation of f ( η ) with η for various values of S .
Symmetry 12 00276 g011
Figure 12. Variation of θ ( η ) with η for various values of S .
Figure 12. Variation of θ ( η ) with η for various values of S .
Symmetry 12 00276 g012
Figure 13. Variation of θ ( η ) with η for various values of R d .
Figure 13. Variation of θ ( η ) with η for various values of R d .
Symmetry 12 00276 g013
Figure 14. Smallest eigenvalues γ for various values of A .
Figure 14. Smallest eigenvalues γ for various values of A .
Symmetry 12 00276 g014

Share and Cite

MDPI and ACS Style

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. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12020276

AMA Style

Lund LA, Omar Z, Khan I, Sherif E-SM. Dual Solutions and Stability Analysis of a Hybrid Nanofluid over a Stretching/Shrinking Sheet Executing MHD Flow. Symmetry. 2020; 12(2):276. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12020276

Chicago/Turabian Style

Lund, Liaquat Ali, Zurni Omar, Ilyas Khan, and El-Sayed M. Sherif. 2020. "Dual Solutions and Stability Analysis of a Hybrid Nanofluid over a Stretching/Shrinking Sheet Executing MHD Flow" Symmetry 12, no. 2: 276. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12020276

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