Next Article in Journal
Subjective and Objective Quality Assessments of Display Products
Next Article in Special Issue
MHD Double-Diffusive Carreau Fluid Flow through a Porous Medium with Variable Thermal Conductivity and Suction/Injection
Previous Article in Journal
A Novel Method to Determine Basic Probability Assignment Based on Adaboost and Its Application in Classification
Previous Article in Special Issue
Thermal Management and Modeling of Forced Convection and Entropy Generation in a Vented Cavity by Simultaneous Use of a Curved Porous Layer and Magnetic Field
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hybrid Nanofluids Flows Determined by a Permeable Power-Law Stretching/Shrinking Sheet Modulated by Orthogonal Surface Shear

Department of Mathematics, Faculty of Mathematics and Computer Science, Babeş-Bolyai University, 400084 Cluj-Napoca, Romania
*
Author to whom correspondence should be addressed.
Entropy 2021, 23(7), 813; https://doi.org/10.3390/e23070813
Submission received: 1 June 2021 / Revised: 23 June 2021 / Accepted: 24 June 2021 / Published: 25 June 2021
(This article belongs to the Special Issue Entropy Analysis in Nanofluids and Porous Media)

Abstract

:
The present paper studies the flow and heat transfer of the hybrid nanofluids flows induced by a permeable power-law stretching/shrinking surface modulated orthogonal surface shear. The governing partial differential equations were converted into non-linear ordinary differential equations by using proper similarity transformations. These equations were then solved applying a numerical technique, namely bvp4c solver in MATLAB. Results of the flow field, temperature distribution, reduced skin friction coefficient and reduced Nusselt number were deduced. It was found that increasing mass flux parameter slows down the velocity and, hence, decreases the temperature. Furthermore, on enlarging the stretching parameter, the velocity and temperature increases and decreases, respectively. In addition, that the radiation parameter can effectively control the thermal boundary layer. Finally, the temperature decreases when the values of the temperature parameter increases. We apply similarity transformation in order to transform the governing model into a system of ODEs (ordinary differential equations). Numerical solutions for particular values of involved parameters are in very good agreement with previous calculations. The most important and interesting result of this paper is that for both the cases of shrinking and stretching sheet flows exhibit dual solutions in some intervals of the shrinking and stretching parameter. In spite of numerous published papers on the flow and heat transfer over a permeable stretching/shrinking surface in nanofluids and hybrid nanofluids, none of the researchers studied the present problem. Therefore, we believe that the results of the present paper are new, and have many industrial applications.

1. Introduction

During the last few years, hybrid nanofluids appeared as an extension of nanofluids and are believed to improve their thermophysical and rheological characteristics and also heat transfer attributes. Suresh et al. [1] supervised an analysis from an experimental point of view on the synthesis and characterization of Al 2 O 3 Cu / H 2 O nanoparticles for various concentrations: 0.1%; 0.33%; 0.75%; 1%; and 2%. Huminic and Huminic [2] have affirmed a review paper that the results concerning the thermo-physical properties and the heat transfer and flow characteristics of hybrid nanofluids are used in various heat exchangers and energy applications. Hybrid nanofluid is invented by combining a base fluid with a mixture or composite form of suspended dissimilar nanoparticles. This kind of nanofluid was found to have great efficiency in terms of its thermophysical properties. The reviews on hybrid nanofluids were briefly discussed in the references [3,4]. Mabood et al. [5] highlighted the use of Fe3O4-graphene/water hybrid nanoliquid in three-dimensional unsteady boundary layer flow with magnetic field and non-linear radiation. Devi and Devi [6] predicted the thermophysical properties of hybrid nanofluid by using the correlations which have been validated with the existing experimental data. Moreover, the correlations of hybrid nanofluids in Takabi and Salehi [7] were also practical and considered by few researchers. Later, Khashi’ie et al. [8] used both types of thermophysical properties (Devi and Devi [6] and Takabi and Salehi [7]). Experimental and numerical results presented in many papers indicate that the hybrid nano-sized particle fluids will enhance the heat transfer in heat exchangers significantly, but research is still needed regarding the study of distinct combinations of nanoparticles, their mixing ratio, the stability of the hybrid nanofluids, and the understanding the mechanisms which contribution to the heat transfer enhancement. Devi and Devi [6] engaged a numerical study of the Cu- Al 2 O 3 /water hybrid nanofluid in two-dimensional flow over a sheet which is stretched. It was found that by choosing various and suitable quantities of nanoparticles the best choice heat transfer rate of hybrid nanofluid will be obtained. Since the paper by Devi and Devi [6] was published, many researchers have studied the hybrid nanofluids past surfaces that are stretched or shrunk, such as Rostami et al. [9], Waini et al. [10,11], Khashi’ie et al. [12], Zainal et al. [13], to mention just a few authors. We can also mention the following published papers suggested by reviewers: Rahman et al. [14], Shafiq [15,16,17], Sanni et al. [18], Khan et al. [19], Khan [20] and Makinde [21], Mabood et al. [5]. The heat transfer applications of hybrid nanofluids are: industrial cooling applications; smart fluids (due to the lack of abundant sources of clean energy and widespread dissemination of battery operated devices, such as cell phones and laptop it is essential to use nanofluids as a smart fluid; biomedical applications; cancer therapeutics (this initiative involves the use of iron-based nanoparticles as delivery vehicles for drugs or radiation in cancer patients); liquid cooling of computer processors due to their high thermal conductivity. emulsions, spray pyrolysis; and thermal spraying, etc.
The boundary layer flows driven by stretching/shrinking sheets have been investigated by a lot of people due to their many applications. The heat transfer and flow characteristics generated by stretching/shrinking surfaces appear widely in engineering application processes (see Fisher [22]). Therefore, in the present paper, we intend to analyze the hybrid nanofluid flows induced by a permeable power-law stretching/shrinking surface modulated by orthogonal surface shear. The influences of various parameters over the main physical quantities of interest are provided in figures and are discussed in detail.

2. Mathematical Model

This paper extends the zero-pressure-gradient boundary-layer solutions for power-law stretching/shrinking sheet studied by Banks [23] to include an orthogonal power-law shearing motion with equal exponents to the important case of hybrid nanoparticle-sized fluids. Figure 1 presents the schematic of the problem, in which ( x ,   y ) axes are taken in the plane of the sheet and z axis is normal directions to the surface of the sheet, the flow being in the region z 0 . We suppose that the velocity of the stretching/shrinking sheet is u w ( x ) and the mass flux velocity is w w ( x , y , 0 ) , with w w ( x , y , 0 ) < 0 for suction case and w w ( x , y , 0 ) > 0 for injection case. We also assume that the constant surface temperature is T w , and T is the temperature of the ambient fluid.
Based on these hypotheses and using hybrid nanofluid model suggested by Takabi and Salehi [7], the governing conservation equations of continuity, streamwise momentum, the spanwise momentum and energy equations, are expressed in the Cartesian coordinates ( x ,   y ,   z ) as (see Weidman [24]; Devi and Devi [6]):
  u     x +   v     y +   w     z = 0
u     u   x + v     u   y + w     u   z = μ h n f ρ h n f   2   u   z 2
u     v   x + v     v   x + w     v   z = μ h n f ρ h n f   2 v   z 2
u     T   x + v     T   y + w   T   z = k h n f ( ρ   C p ) h n f   2 T   z 2
subject to the boundary conditions:
u ( x , y , z ) = u w ( x , y , z )   λ = a   x n λ ,       v ( x , y , z ) = v w ( x , y , z ) = b x n w = w w ( x , y , z ) , T = T w )   at   z = 0 u ( x , y , z ) 0 , v ( x , y , z ) 0 , w ( x , y , z ) 0 , T T as z }
Here u ( x , y , z ) ,   v ( x , y , z )   and   w ( x , y , z ) represent the components of the velocity along ( x , y ,   z ) axes, T is the temperature of the hybrid nanofluid, a is a positive constant having dimension t 1 , where a describes the stretching/shrinking strength of the sheet and λ is the stretching/shrinking parameter of the sheet which is constant, with λ > 0 for the stretching case, λ < 0 for the shrinking case and λ = 0 for the fixed sheet, respectively, and n is the power-law stretching/shrinking parameter not being required to be an integer number.
Furthermore, μ h n f represents the dynamic viscosity, ρ n h f represents the density, k h n f represents the thermal conductivity and ( ρ   C p ) h n f represents the heat capacity of the hybrid nanoparticles fluid, calculated accordingly to Takabi and Salehi [7].
μ h n f = μ f ( 1 ϕ A l 2 O 2 ϕ C u ) 2.5 ρ h n f = ϕ A l 2 O 3 ρ A l 2 O 3 + ϕ C u ρ C u + ( 1 ϕ h n f ) ρ f   k h n f k f = { ϕ A l 2 O 3 k A l 2 O 3 + ϕ C u k C u ϕ A l 2 O 3 + ϕ C u + 2 k f + 2 ( ϕ A l 2 O 3 k A l 2 O 3 + ϕ C u k C u ) 2 ( ϕ A l 2 O 3 + ϕ C u ) k f }   × { ϕ A l 2 O 3 k A l 2 O 3 + ϕ C u k C u ϕ A l 2 O 3 + ϕ C u + 2 k f ( ϕ A l 2 O 3 k A l 2 O 3 + ϕ C u k C u ) + ( ϕ A l 2 O 3 + ϕ C u ) k f } 1 ( ρ C p ) h n f = ϕ A l 2 O 3 ( ρ C p ) A l 2 O 3 + ϕ C u ( ρ C p ) C u + ( 1 ϕ h n f ) ( ρ C p ) f }
Here ϕ is the volume fraction of the nanoparticles ( ϕ = 0 corresponding to a normal fluid), ϕ A l 2 O 3 denotes Al2O3 particles and ϕ C u denotes Cu particles), ρ f and ρ s represents the densities of the base fluid and the hybrid nano-sized particles, respectively, kf and ks represents the thermal conductivities of the base fluid and the hybrid nano-sized particles, respectively, ( ρ C p ) f and ( ρ C p ) s represents the heat capacitance of the base fluid and the hybrid nano-sized particle, and C p is the heat capacity at constant pressure. The above relations are correct and realistic. They are determined starting from a physical hypothesis, and they agree with the mass and energy conservation. Table 1 presents in detail the physical properties of the base fluid (water) and nano-sized particles that are of great interest for us.
According to Weidman [24], the next variables which are dimensionless are introduced:
  u ( x , y , η ) = a   x n f ( η ) ,       v ( x , y , η ) = b   x n   g ( η ) , w =   a   ν f μ f { n   x n 1   f ( η ) + x n σ x μ f   [ η   f ( η ) f ( η ) ] } θ ( η ) = T   T T w     T ,           η = z   a / ν f   σ ( x ) }
Thus,
w w ( x , 0 ) = a   ν f μ f   ( n   x n 1 x n σ x μ f ) S
where the differentiation with respect to η is denoted by prime, S is the constant mass flux velocity parameter. Here S > 0 represents suction and S < 0 represents injection, respectively, and σ ( x ) = ( n + 1 ) / 2   x ( n 1 ) / 2 .
Substituting the similarity variables (7) into Equations (1)–(3), the following ordinary differential equations are obtained:
μ h n f / μ f ρ h n f / ρ f   f + f   f 2   n n + 1 f 2 = 0
μ h n f / μ f ρ h n f / ρ f   g + f   g 2   n n +   1 f g = 0
1 P r k h n f / k f ( ρ C p ) h n f / ( ρ C p ) f θ + 2   n n + 1   f θ = 0
with the corresponding boundary conditions:
f ( 0 ) = S ,       f ( 0 ) = λ ,         g ( 0 ) = 0 ,     g ( 0 ) = 1 ,     θ ( 0 ) = 1   f ( η ) 0 ,         g ( η ) 0 ,         θ ( η ) 0           as         η           }
where P r is the Prandtl number, which are defined as P r = ( ρ   C p ) f / k f .
From a physical point of view we are interested in the skin friction coefficients C f x and C f y , and the Nusselt local number N u , presented as follows:
C f x = μ h n f ρ f   (   u   z ) z = 0 ,         C f y = μ h n f ρ f   (   v   z ) z = 0 ,           N u x = k h n f k f   (   T   z ) z = 0
Using (7) and (13), we get
R e x 1 / 2 C f x = μ h n f μ f   f ( 0 ) ,         R e x 1 / 2 C f y = μ h n f μ f   g ( 0 ) ,   R e x 1 / 2   N u x =   k h n f k f   θ ( 0 ) }
where R e x = u w ( x )   x / ν f and R e y = u w ( y )   y / ν f are the local Reynolds numbers.
We can observe that for a normal fluid ( ϕ = 0 ) , Equation (9) is similar with Equation (6) from Mahapatra et al. [26], namely:
f + f   f 2   n n + 1 f 2 = 0
subject to:
f ( 0 ) = 0 ,       f ( 0 ) = 1 ,         f ( η ) 0         as       η
However, for a regular fluid ( ϕ = 0 ) and n = 1 , Equation (9) reduces to:
f + f   f f 2 = 0
subject to:
f ( 0 ) = S ,         f ( 0 ) = λ ,         f ( η ) 0           a s         η
The exact solution of the boundary value problem (17) and (18) is given by (Vajravelu and Cannon [27] or Cortell [28]), as:
f ( η ) = S + α   [ 1 exp ( β   η ) ] ,               ( β = S + α > 0 )
where:
β = 1 2   ( S   ±   S 2 + 4   λ )
Thus, we have:
f ( 0 ) = λ 2   ( S   ±   s 2 + 4   λ )
so that it gives, as expected, λ c = s 2 4 .

3. Discussion of the Results

Using the bvp4c solver from MATLAB programming language, Equations (9)–(11) and the conditions (12) at the boundaries are numerically solved for the value of the Prandtl number fixed to Pr = 6.2 corresponding to pure water at 25 °C. In particular, bvp4c is a finite-difference code that implements the three-stage Lobatto IIIa formula (see Shampine et al. [29]). This is a collocation formula that provides a continuous solution with fourth-order accuracy. Mesh selection and error control are based on the residual of the continuous solution. The effectiveness of this solver ultimately counts on our ability to provide the algorithm with an initial prediction for the numerical result. Due to the fact that our problem may have multiple (dual) solutions, the bvp4c solver has to be provided with an initial estimate of the solution for the problem (9)–(12). The size of our mesh and the tolerance of the error are 0.001 and 10 9 respectively, ensuring accurate enough numerical results. Based on this guess value, the velocity and temperature contours are determined satisfying the boundary conditions (12) at infinity. Because the bvp4c method will converge to the first solution even for poor guesses we can use a simple initial guess for the first solution. However, for the second solution of Equations (9)–(12), it is more complicated to determine a good guess which will be adequate. Also, this convergence issue is influenced by the value of the selected parameters. In this paper we studied both stretching case sheet ( λ > 0 ) and shrinking ( λ < 0 ) case sheet. The obtained numerical results for the reduced skin friction f ( 0 ) can be compared with those reported by Vajravelu and Cannon [27] or Cortell [28], and they are shown in Table 2. Since our results match very well we are confident that our solutions are correct and precise.
Graphs that display the influence of the governing parameters on the reduced skin friction coefficients f′(0) and g′(0), the reduced local Nusselt θ ( 0 ) , velocity f ( η ) , g ( η ) as well as the temperature θ ( η ) contours were obtained and displayed in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13. For numerical results we considered the non-dimensional parameter values as: Prandtl number P r , the power-law parameter n = 1 3 , the stretching/shrinking parameter λ c   λ < 1 in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 and λ = 0.95 in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13, and the constant mass flux velocity parameter S = 2 . Also, different nanofluid particles concentrations are considered in this paper. Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 show that for the studied problem (9)–(12) dual solutions (upper and lower branch solutions) exist for λ c λ (shrinking and stretching sheet) and there is no solution in the case that λ < λ c < 0 . Here λ c < 0 represents the critical value of λ < 0 starting at the ODE’s system (9) to (12) has no solutions. We should assert that in the situation that λ < λ c < 0 the full equations (Navier–Stokes and energy) must be resolved.
Finally, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 present the velocities f ( η ) and g ( η ) contours and the temperature θ ( η ) contours for λ = 0.95 (shrinking sheet situation). In these figures the upper branch solution is illustrated by solid lines, while the lower branch solution is depicted by the dotted lines, respectively. It is obvious looking at these figures that the conditions at infinity (12) are verified. This shows the correctness of our numerical results for the boundary-value problem (9)–(12). We deduce from Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13 that the dimensionless contours for velocity and temperature of the upper branch solution show a thinner boundary layer thickness in comparison to the lower branch one.

4. Conclusions

Our main conclusions are the following:
  • For the potential stretching and shrinking flows there exist two solutions (first/upper branch and second/lower branch), when some of the governing parameters lie in specific intervals.
  • We can control the separation of the boundary layer by using lower magnitude of λ < 0 .
  • The increment of the copper or alumina nano-sized particles solid volume concentration has a large influence on the critical value λ c < 0 of λ < 0 .
  • The increase of concentration ϕ 2 ( C u ) raises the velocity of the fluid along the x direction for the shrinking case, while the velocity of the fluid along the y direction shows an opponent result. The increment of concentration ϕ 1 ( A l 2 O 3 ) amplify the velocity of the fluid in the y direction for the shrinking case, whereas the velocity of the fluid in the x direction shows a contrary result. The above conclusions are true for the upper branch solution. For the lower branch solution the results show contrary behaviors.
  • The increase in n comes along with a contraction of the temperature contours for the shrinking sheet situation for both upper and lower branch solutions.
  • The increase in concentration ϕ 1 ( A l 2 O 3 ) comes along with a small dilatation of the temperature contours for the shrinking case for the first solution. A contrary behavior is seen for the second solution.

Author Contributions

Conceptualization, N.C.R. and I.P.; methodology, N.C.R. and I.P.; software, N.C.R. and I.P.; validation, N.C.R. and I.P.; formal analysis, N.C.R. and I.P.; investigation, N.C.R. and I.P.; resources, N.C.R. and I.P.; data curation, N.C.R. and I.P.; writing—original draft preparation, N.C.R. and I.P.; writing—review and editing, N.C.R. and I.P.; visualization, N.C.R. and I.P.; supervision, N.C.R. and I.P.; project administration, N.C.R. and I.P.; funding acquisition, N.C.R. and I.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors wish to express their thanks to the very competent reviewers for their excellent comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Roman letters
a positive constant
b share strength constant
C f x , C f y skin friction coefficient along the x- and y- directions, respectively
C p specific heat at constant pressure ( Jkg 1 K 1 )
f ( η ) , g ( η ) dimensionless stream function in the x- and y- directions, respectively ( ms 1 )
k thermal conductivity of the fluid ( Wm 1 K 1 )
n is the power-law stretching/shrinking parameter
N u x local Nusselt number
P r Prandtl number
Re x , R e y local Reynolds number in the x- and y- directions, respectively
S suction parameter
t time ( s )
T fluid temperature ( K )
T w wall temperature ( K )
T ambient temperature ( K )
u , v , w velocities in the x-, y- and z- directions ( ms 1 )
u w , v w velocities of the stretching/shrinking sheet in the x- and y- directions, respectively ( ms 1 )
x , y , z Cartesian coordinates ( m )
Greek symbols
θ dimensionless temperature
λ stretching/shrinking parameter
η similarity variable
μ dynamic viscosity of the fluid ( kgm 1 s 1 )
ν kinematic viscosity of the fluid ( m 2 s 1 )
ρ fluid density ( kgm 3 )
ρ C p heat capacity +( JK 1 m 3 )
ϕ 1 alumina volume fractions
ϕ 2 copper volume fractions
Subscripts
f base fluid
n f nanofluid
h n f hybrid nanofluid
s 1 solid component for alumina (Al2O3)
s 2 solid component for copper (Cu)
Superscript
differentiation with respect to η

References

  1. Suresh, S.; Venkitaraj, K.P.; Selvakumar, P. Synthesis, characterization of Al2O3-Cu nanocomposite powder and water-based nanofluids. Adv. Mater. Res. 2011, 328–330, 1560–1567. [Google Scholar] [CrossRef]
  2. 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]
  3. Mahian, O.; Kolsi, L.; Amani, M.; Estellé, P.; Ahmadi, G.; Kleinstreuer, C.; Marshall, J.S.; Siavashi, M.; Taylor, R.A.; Niazmand, H.; et al. Recent advances in modeling and simulation of nanofluid flows-Part I: Fundamentals and theory. Physic Rep. 2019, 790, 1–48. [Google Scholar] [CrossRef]
  4. Mahian, O.; Kolsi, L.; Amani, M.; Estellé, P.; Ahmadi, G.; Kleinstreuer, C.; Marshall, J.S.; Taylor, R.A.; Abu-Nada, E.; Rashidi, S.; et al. Recent advances in modeling and simulation of nanofluid flows-Part II: Applications. Physic Rep. 2019, 791, 1–59. [Google Scholar] [CrossRef]
  5. Mabood, F.; Khan, W.; Ismail, A. MHD boundary layer flow and heat transfer of nanofluids over a nonlinear stretching sheet: A numerical study. J. Magn. Magn. Mater. 2015, 374, 569–576. [Google Scholar] [CrossRef]
  6. Devi, S.U.; Devi, S.P.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]
  7. Takabi, B.; Salehi, S. Augmentation of the heat transfer performance of a sinusoidal corrugated enclosure by employing hybrid nanofluid. Adv. Mech. Engng. 2014, 6, 147059. [Google Scholar] [CrossRef]
  8. 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]
  9. Rostami, M.N.; Dinarvand, S.; Pop, I. Dual solutions for mixed convective stagnation-point flow of an aqueous silica–alumina hybrid nanofluid. Chin. J. Phys. 2018, 56, 2465–2478. [Google Scholar] [CrossRef]
  10. Waini, I.; Ishak, A.; Pop, I. MHD flow and heat transfer of a hybrid nanofluid past a permeable stretching/shrinking wedge. Appl. Math. Mech.-Engl. Ed. 2020, 41, 507–520. [Google Scholar] [CrossRef]
  11. Waini, I.; Ishak, A.; Pop, I. Hybrid nanofluid flow towards a stagnation point on an exponentially stretching/shrinking vertical sheet with buoyancy effects. Int. J. Numer. Methods Heat Fluid Flow 2021, 31, 216–235. [Google Scholar] [CrossRef]
  12. Khashi’Ie, N.S.; Arifin, N.M.; Pop, I.; Nazar, R.; Hafidzuddin, E.H.; Wahi, N. Three-dimensional hybrid nanofluid flow and heat transfer past a permeable stretching/shrinking sheet with velocity slip and convective condition. Chin. J. Phys. 2020, 66, 157–171. [Google Scholar] [CrossRef]
  13. Zainal, N.A.; Nazar, R.; Naganthran, K.; Pop, I. Unsteady three-dimensional MHD non-axisymmetric Homann stagnation point flow of a hybrid nanofluid with stability analysis. Mathematics 2020, 8, 784. [Google Scholar] [CrossRef]
  14. Rahman, M.M.; Saghir, Z.; Pop, I. Free convective heat transfer efficiency in Al2O3–Cu/water hybrid nanofluid inside a recto-trapezoidal enclosure. Int. J. Numer. Methods Heat Fluid Flow 2021. [Google Scholar] [CrossRef]
  15. Shafiq, A.; Sindhu, T.N.; Al-Mdallal, Q.M. A sensitivity study on carbon nanotubes significance in Darcy–Forchheimer flow towards a rotating disk by response surface methodology. Sci. Rep. 2021, 11, 1–26. [Google Scholar] [CrossRef] [PubMed]
  16. Shafiq, A.; Mebarek-Oudina, F.; Sindhu, T.N.; Abidi, A. A study of dual stratification on stagnation point Walters’ B nanofluid flow via radiative Riga plate: A statistical approach. Eur. Phys. J. Plus 2021, 136, 1–24. [Google Scholar] [CrossRef]
  17. Anum, S.; Tabassum, N.S.; Chaudry, M.K. Numerical investigation and sensitivity analysis on bioconvective tangent hyperbolic nanofluid flow towards stretching surface by response surface methodology. Alex. Eng. J. 2020, 59, 4533–4548. [Google Scholar]
  18. Sanni, K.M.; Asghar, S.; Rashid, S.; Chu, Y.-M. Nonlinear radiative treatment of hydromagnetic non-Newtonian fluid flow induced Kehinde nonlinear convective boundary driven curved sheet with dissipations and chemical reaction effects. Front. Phys. 2021, 9, 283. [Google Scholar] [CrossRef]
  19. Khan, W.A.; Khan, Z.H.; Rahi, M. Fluid flow and heat transfer of carbon nanotubes along a flat plate with Navier slip boundary. Appl. Nanosci. 2013, 4, 633–641. [Google Scholar] [CrossRef] [Green Version]
  20. Khan, W.; Culham, J.R.; Yovanovich, M.M. Fluid flow and heat transfer from a cylinder between parallel planes. J. Thermophys. Heat Transf. 2004, 18, 395–403. [Google Scholar] [CrossRef] [Green Version]
  21. Makinde, O.D.; Khan, Z.H.; Ahmad, R.; Khan, W.A. Numerical study of unsteady hydromagnetic radiating fluid flow past a slippery stretching sheet embedded in a porous medium. Phys. Fluids 2018, 30, 083601. [Google Scholar] [CrossRef]
  22. Fisher, E.G. Extrusion of Plastics; Wiley: New York, NY, USA, 1976. [Google Scholar]
  23. Banks, W.H.H. Similarity solutions of the boundary-layer euations for a stretching wall. J. Mécnique Théorique Appl. 1983, 2, 375–392. [Google Scholar]
  24. Weidman, P.D. Flows induced by power-law stretching surface motion modulated by transverse or orthogonal surface shear. Comptes Rendus Mec. 2017, 345, 169–176. [Google Scholar] [CrossRef]
  25. 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]
  26. Mahapatra, T.R.; Nandy, S.K.; Gupta, A.S. Oblique stagnation-point flow and heat towards a shrinking sheet with thermal radiation. Meccanica 2012, 47, 1623–1632. [Google Scholar] [CrossRef]
  27. Vajravelu, K.; Cannon, J.R. Flow over a nonlinearly stretching sheet. Appl. Math. Comput. 2006, 181, 609–618. [Google Scholar] [CrossRef]
  28. Cortell, R. A nonlinearly stretching sheet. Appl. Math. Comput. 2007, 184, 864–873. [Google Scholar]
  29. Shampine, L.F.; Gladwell, I.; Thompson, S. Solving ODEs with MATLAB; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
Figure 1. Physical scheme the proposed problem.
Figure 1. Physical scheme the proposed problem.
Entropy 23 00813 g001
Figure 2. Change in f ( 0 ) with λ for some values of n for the hybrid nanoparticles fluid having ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025 for S = 2 .
Figure 2. Change in f ( 0 ) with λ for some values of n for the hybrid nanoparticles fluid having ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025 for S = 2 .
Entropy 23 00813 g002
Figure 3. Change in g ( 0 ) with λ for some values of n for the hybrid nanoparticles fluid having ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025 for S = 2 .
Figure 3. Change in g ( 0 ) with λ for some values of n for the hybrid nanoparticles fluid having ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025 for S = 2 .
Entropy 23 00813 g003
Figure 4. Change in θ ( 0 ) with λ for some values of n for the hybrid nanoparticles fluid having ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025 for S = 2 .
Figure 4. Change in θ ( 0 ) with λ for some values of n for the hybrid nanoparticles fluid having ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025 for S = 2 .
Entropy 23 00813 g004
Figure 5. Change in f ( 0 ) with parameter λ for the case of hybrid nanoparticles fluid having distinct concentrations of Al2O3 and Cu nanoparticles when n = 1 , S = 2 .
Figure 5. Change in f ( 0 ) with parameter λ for the case of hybrid nanoparticles fluid having distinct concentrations of Al2O3 and Cu nanoparticles when n = 1 , S = 2 .
Entropy 23 00813 g005
Figure 6. Change in g ( 0 ) with parameter λ for the case of hybrid nanoparticles fluid having distinct concentrations of Al2O3 and Cu nanoparticles when n = 1 , S = 2 .
Figure 6. Change in g ( 0 ) with parameter λ for the case of hybrid nanoparticles fluid having distinct concentrations of Al2O3 and Cu nanoparticles when n = 1 , S = 2 .
Entropy 23 00813 g006
Figure 7. Change in θ ( 0 ) with parameter λ for the case of hybrid nanoparticles fluid having distinct concentrations of Al2O3 and Cu nanoparticles when n = 1 , S = 2 .
Figure 7. Change in θ ( 0 ) with parameter λ for the case of hybrid nanoparticles fluid having distinct concentrations of Al2O3 and Cu nanoparticles when n = 1 , S = 2 .
Entropy 23 00813 g007
Figure 8. Change on the dimensionless velocity f ( η ) for some n when λ = 0.95 , S = 2 for concentrations ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025.
Figure 8. Change on the dimensionless velocity f ( η ) for some n when λ = 0.95 , S = 2 for concentrations ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025.
Entropy 23 00813 g008
Figure 9. Change on the dimensionless velocity g ( η ) for some n when λ = 0.95 , S = 2 for concentrations ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025.
Figure 9. Change on the dimensionless velocity g ( η ) for some n when λ = 0.95 , S = 2 for concentrations ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025.
Entropy 23 00813 g009
Figure 10. Change on the dimensionless temperature θ ( η ) for some n when λ = 0.95 , S = 2 for concentrations ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025.
Figure 10. Change on the dimensionless temperature θ ( η ) for some n when λ = 0.95 , S = 2 for concentrations ϕ1 (Al2O3) = 0.025 and ϕ2 (Cu) = 0.025.
Entropy 23 00813 g010
Figure 11. Change in the dimensionless velocity f ( η ) for λ = 0.95 , S = 2 and n = 1 in the case of distinct concentrations of nanoparticles.
Figure 11. Change in the dimensionless velocity f ( η ) for λ = 0.95 , S = 2 and n = 1 in the case of distinct concentrations of nanoparticles.
Entropy 23 00813 g011
Figure 12. Change in the dimensionless velocity g ( η ) for λ = 0.95 , S = 2 and n = 1 in the case of distinct concentrations of nanoparticles.
Figure 12. Change in the dimensionless velocity g ( η ) for λ = 0.95 , S = 2 and n = 1 in the case of distinct concentrations of nanoparticles.
Entropy 23 00813 g012
Figure 13. Change in the dimensionless temperature θ ( η ) for λ = 0.95 , S = 2 and n = 1 in the case of distinct concentrations of nanoparticles.
Figure 13. Change in the dimensionless temperature θ ( η ) for λ = 0.95 , S = 2 and n = 1 in the case of distinct concentrations of nanoparticles.
Entropy 23 00813 g013
Table 1. Thermo-physical characteristics (Oztop and Abu-Nada [25]).
Table 1. Thermo-physical characteristics (Oztop and Abu-Nada [25]).
Physical CharacteristicsWater A l 2 O 3 C u
ϱ   ( kg / m 3 ) 997.039708933
C p ( J / kgK ) 4180765385
k   ( W / mK ) 0.607140400
Table 2. Comparisons of f ( 0 ) for some n, in the case ϕ 1 = ϕ 2 = 0 , S = 0 and λ = 1 .
Table 2. Comparisons of f ( 0 ) for some n, in the case ϕ 1 = ϕ 2 = 0 , S = 0 and λ = 1 .
n f ( 0 )
Vajravelu and Cannon [27]Cortell [28]Present Study
0-0.6275470.627554
0.2-0.7667580.766837
0.5-0.8894770.889543
0.75-0.9537860.953956
111.01.0
1.5-1.0615871.061600
3-1.1485881.148593
51.1945-1.194487
7-1.2168471.216850
101.23481.2348751.234874
20-1.2574181.257423
100-1.2767681.276773
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Roşca, N.C.; Pop, I. Hybrid Nanofluids Flows Determined by a Permeable Power-Law Stretching/Shrinking Sheet Modulated by Orthogonal Surface Shear. Entropy 2021, 23, 813. https://doi.org/10.3390/e23070813

AMA Style

Roşca NC, Pop I. Hybrid Nanofluids Flows Determined by a Permeable Power-Law Stretching/Shrinking Sheet Modulated by Orthogonal Surface Shear. Entropy. 2021; 23(7):813. https://doi.org/10.3390/e23070813

Chicago/Turabian Style

Roşca, Natalia C., and Ioan Pop. 2021. "Hybrid Nanofluids Flows Determined by a Permeable Power-Law Stretching/Shrinking Sheet Modulated by Orthogonal Surface Shear" Entropy 23, no. 7: 813. https://doi.org/10.3390/e23070813

APA Style

Roşca, N. C., & Pop, I. (2021). Hybrid Nanofluids Flows Determined by a Permeable Power-Law Stretching/Shrinking Sheet Modulated by Orthogonal Surface Shear. Entropy, 23(7), 813. https://doi.org/10.3390/e23070813

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